| 查看: 3977 | 回复: 11 | ||
[求助]
如何快速求解大型矩阵全部特征值与特征向量？ 已有3人参与
|
||
|
矩阵是带状，还是旋转对称的方阵，规模大于10k*10k。 发自小木虫Android客户端 |
回复此楼» 猜你喜欢
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
-
5个%2F，非固定段22个字母
已经有13人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有5人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有6人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有7人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有5人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- ANSYS14.0
已经有15人回复
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
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
virtualzx
木虫 (著名写手)
- 应助: 263 (大学生)
- 金币: 7161.3
- 红花: 54
- 帖子: 1605
- 在线: 317.6小时
- 虫号: 2069080
- 注册: 2012-10-18
- 性别: GG
- 专业: 理论和计算化学
7楼2016-04-15 11:18:34
8楼2016-04-15 11:46:03
9楼2016-04-15 20:18:14
virtualzx
木虫 (著名写手)
- 应助: 263 (大学生)
- 金币: 7161.3
- 红花: 54
- 帖子: 1605
- 在线: 317.6小时
- 虫号: 2069080
- 注册: 2012-10-18
- 性别: GG
- 专业: 理论和计算化学
【答案】应助回帖
感谢参与，应助指数 +1
|
dsyevd是求解密集对称矩阵特征值的函数。matlab的eig函数只是一个空壳，它会根据你的矩阵类型选择一个lapack函数调用。但是matlab只给出了密集和三对角矩阵的接口。你用的是带状矩阵，lapack里有专门的函数，但是你无法从matlab里调用。这样就没法利用带状矩阵的特性，速度稍慢，但也应该完全有能力处理。 发自小木虫IOS客户端 |
» 本帖已获得的红花（最新10朵）
yueshuai97710楼2016-04-15 22:14:26
3楼2016-04-14 22:40:46
5楼2016-04-15 10:07:48
6楼2016-04-15 10:08:44
10