| 查看: 685 | 回复: 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 ] |
回复此楼» 猜你喜欢
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有184人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
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