| 查看: 3760 | 回复: 11 | |||
[求助]
如何快速求解大型矩阵全部特征值与特征向量？ 已有3人参与
|
|||
|
矩阵是带状，还是旋转对称的方阵，规模大于10k*10k。 发自小木虫Android客户端 |
回复此楼» 猜你喜欢
26考研一志愿中国石油大学(华东)305分求调剂
已经有4人回复
-
286求调剂
已经有5人回复
-
材料与化工（0856）304求 B区 调剂
已经有3人回复
-
268求调剂
已经有8人回复
-
265求调剂
已经有14人回复
-
350 本科985求调剂，求老登收留
已经有3人回复
-
【考研调剂】化学专业 281分，一志愿四川大学，诚心求调剂
已经有9人回复
-
初试 317
已经有3人回复
-
279分求调剂 一志愿211
已经有15人回复
-
279求调剂
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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
