| 查看: 837 | 回复: 11 | |||
| 当前主题已经存档。 | |||
[交流]
【讨论】【Matlab】|稀疏矩阵的处理|“十金散尽还复来”
|
|||
最近用Matlab处理比较大矩阵的本征值问题，迫于机器条件限制（内存和速度不够），想跟大家讨论一下如下几个稀疏矩阵的存储和计算问题      计算的时候，前几个矩阵相乘得到最后一个矩阵，然后求解最后一个矩阵的本征值问题。使用 sparse 讲矩阵转化为稀疏矩阵后，运算速度快多了。问题是， 1°稀疏矩阵的本征值好像只能用 eigs 这个命令来求，想要求最小的几个本征值，用了 eig(A, k, sigma)这个命令，sigma 取0，结果在边界出现了畸变（肉眼可见）； 2°运行时 command window 老提示 "A few Ritz values of..."，刷屏很厉害； 3°如果对稀疏矩阵的某些零元素进行赋值，是不是要先将其还原为满矩阵呢？ 怎么解决以上问题，恳请各位达人指教，谢谢！ [ Last edited by Yaniel on 2009-5-21 at 15:05 ] |
回复此楼» 猜你喜欢
中国科学院东莞材料科学与技术研究所-2026年博士招生-吴昊研究员-磁学与自旋电子学
已经有0人回复
-
《电磁学》教材推荐
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有209人回复
-
【急招】合肥工大核聚变材料计算方向2026级工程博士生
已经有4人回复
-
大豆异黄酮分离
已经有0人回复
-
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标
已经有10人回复
-
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名
已经有1人回复
-
中国科学院物理研究所谌志国研究员团队招收2027年博士研究生
已经有3人回复
-
2026年中德博士后交流项目 - 新型量子和磁性材料：材料制备表征和中子散射研究
已经有12人回复
2楼2009-05-21 15:06:26
zeoliters
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 贵宾: 0.902
- 金币: 2088.4
- 散金: 16
- 帖子: 445
- 在线: 60.6小时
- 虫号: 772181
- 注册: 2009-05-16
- 专业: 计算化学，量子化学
3楼2009-05-21 17:12:09
4楼2009-05-21 17:57:40
redcrag
木虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 1778.2
- 红花: 1
- 帖子: 1163
- 在线: 45小时
- 虫号: 254694
- 注册: 2006-05-27
- 性别: GG
- 专业: 光学
★ ★ ★ ★ ★ ★
Yaniel(金币+2,VIP+0):谢谢 ， 我先去试试看 5-22 09:09
Yaniel(金币+2,VIP+0):options.disp=0 解决了刷屏问题，谢谢！ 5-27 14:00
zhaohscas1976(金币+2,VIP+0):谢谢 ， 6-10 11:52
Yaniel(金币+2,VIP+0):谢谢 ， 我先去试试看 5-22 09:09
Yaniel(金币+2,VIP+0):options.disp=0 解决了刷屏问题，谢谢！ 5-27 14:00
zhaohscas1976(金币+2,VIP+0):谢谢 ， 6-10 11:52
|
我见别人光波导模式的差分程序中eigs的用法如下： options.tol = 1e-8; options.disp = 0; % suppress output [V,D]=eigs(P,speye(size(P)),nmodes,shift,options); P是稀疏矩阵，第二个是单位矩阵，第三个是特征值数，第四个是偏移值（在这个值附近求nmodes个特征值）。 因此针对你的问题： 1 好像eigs是求最大特征值，如果求最小的，可能需要对矩阵进行适当变换，将最小特征值变为最大的。 2 我贴的那个第二行可能是不让eigs输出一些信息的 3 那个赋值应该不用转换矩阵，否则你的内存会不够的 |
5楼2009-05-22 08:54:21
6楼2009-05-22 09:13:01
7楼2009-05-22 15:59:21
8楼2009-05-23 13:29:05
9楼2009-05-24 14:35:00
10楼2009-05-26 09:02:08
10