| 查看: 380 | 回复: 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; |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有56人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有9人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有3人回复
-
关于Filecode分析方法
已经有13人回复
-
时间戳又变了8-15
已经有9人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有8人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有9人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有6人回复
4楼2009-12-24 23:44:47
3楼2009-12-24 23:43:29
