| 查看: 950 | 回复: 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 ] |
回复此楼» 猜你喜欢
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后
已经有673人回复
-
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有130人回复
-
物理学I论文润色/翻译怎么收费?
已经有74人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有6人回复
2楼2009-05-21 15:06:26
zeoliters
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 贵宾: 0.902
- 金币: 2088.4
- 散金: 16
- 帖子: 445
- 在线: 60.9小时
- 虫号: 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
