| 查看: 4028 | 回复: 11 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
如何快速求解大型矩阵全部特征值与特征向量？ 已有3人参与
|
||
|
矩阵是带状，还是旋转对称的方阵，规模大于10k*10k。 发自小木虫Android客户端 |
回复此楼» 猜你喜欢
今天系统多次维护，明天很可能放榜！
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有74人回复
-
岩土工程学报什么时候才能终审完呐
已经有4人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- ANSYS14.0
已经有15人回复
virtualzx
木虫 (著名写手)
- 应助: 263 (大学生)
- 金币: 7161.3
- 红花: 54
- 帖子: 1605
- 在线: 317.6小时
- 虫号: 2069080
- 注册: 2012-10-18
- 性别: GG
- 专业: 理论和计算化学
【答案】应助回帖
感谢参与，应助指数 +1
|
dsyevd是求解密集对称矩阵特征值的函数。matlab的eig函数只是一个空壳，它会根据你的矩阵类型选择一个lapack函数调用。但是matlab只给出了密集和三对角矩阵的接口。你用的是带状矩阵，lapack里有专门的函数，但是你无法从matlab里调用。这样就没法利用带状矩阵的特性，速度稍慢，但也应该完全有能力处理。 发自小木虫IOS客户端 |
yueshuai97710楼2016-04-15 22:14:26
cooooldog
铁杆木虫 (著名写手)
ส็็็
- 应助: 237 (大学生)
- 金币: 6101.9
- 散金: 1114
- 红花: 39
- 帖子: 1380
- 在线: 553.8小时
- 虫号: 506699
- 注册: 2008-02-18
- 专业: 模式识别
【答案】应助回帖
感谢参与，应助指数 +1
|
这类问题对硬件要求高，对算法要求也高。 首先保证硬件足够好，通常要多核心、64位、高配置的内存的电脑 其次，软件不建议使用matlab,除非你的应该足够强，比如可以用超算计算机，否则还是用C++/Fortran之类的库； 如果用C++, 有很多现成的库可以尝试。 10K的规模，还是可以尝试用Eigen C++ template library的（不过不一定行） 专门用来解决稀疏特征值问题的开源的库，及其例子： The current version of SLEPc is 3.6, released in June, 2015 SLEPc, the Scalable Library for Eigenvalue Problem Computations http://slepc.upv.es/ Eigenvalue Problem Solver - EPS: Examples http://slepc.upv.es/documentatio ... ages/EPS/index.html Eigenvalue Problem Solver - EPS: Examples http://slepc.upv.es/documentatio ... utorials/index.html |
2楼2016-04-14 18:01:30
3楼2016-04-14 22:40:46
virtualzx
木虫 (著名写手)
- 应助: 263 (大学生)
- 金币: 7161.3
- 红花: 54
- 帖子: 1605
- 在线: 317.6小时
- 虫号: 2069080
- 注册: 2012-10-18
- 性别: GG
- 专业: 理论和计算化学
|
10k带状矩阵不大，正常电脑的话用lapack带的DSBTRD子程序转化为三对角矩阵，然后用DSTERF解特征值就是最快的选项了。这么小的矩阵不需要使用大矩阵的方法。大矩阵方法都不适合解全部特征向量 发自小木虫IOS客户端 |
4楼2016-04-15 09:15:19
