| 查看: 369 | 回复: 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邪修
已经有6人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有14人回复
-
2026国自然放榜时间
已经有4人回复
-
关于代码变化问题，想知道的进来
已经有22人回复
-
这样的filecode谁见过
已经有21人回复
-
filecode与中标关系的预测
已经有5人回复
-
fileCode有新解读？
已经有18人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
孤鸿飘逸
金虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 5861.4
- 散金: 499
- 红花: 17
- 帖子: 617
- 在线: 432.5小时
- 虫号: 864807
- 注册: 2009-10-07
- 性别: GG
- 专业: 机器人学及机器人技术
2楼2009-12-24 18:57:45
3楼2009-12-24 23:43:29
4楼2009-12-24 23:44:47
20