| 查看: 303 | 回复: 3 | |||
| 当前主题已经存档。 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 wuhaoshrek 的 3 个金币 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】大家帮我看看这个程序啊？改进的乘幂法
|
|||
|
我编的改进的乘幂法，算一般的矩阵特征值和特征向量是没有问题的，但是求三对角矩阵，如 2 -1 0 0 0 0 -1 2 -1 0 0 0 0 -1 2 -1 0 0 0 0 -1 2 -1 0 0 0 0 -1 2 -1 0 0 0 0 -1 2 时，就总是会出错啊，不知道怎么回事啊？谢谢大家了！ 源程序： function [mt,my]=maxtr2(A,ep) n=length(A); v0=ones(n,1); [p0,I0]=max(abs(v0)); p0=p0*sign(v0(I0)); v=A*v0; [p,I]=max(abs(v)); p=p*sign(v(I)); k=1; while abs(p-p0)>ep; v0=v; p0=p; u=v/p; v=A*u; [p,I]=max(abs(v)); p=p*sign(v(I)); k=k+1; end mt=p; my=u; |
回复此楼» 猜你喜欢
假如你的研究生提出不合理要求
已经有9人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
实验室接单子
已经有4人回复
-
全日制（定向）博士
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
3楼2009-12-24 23:43:29
4楼2009-12-24 23:44:47