| 查看: 414 | 回复: 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; |
回复此楼» 猜你喜欢
论文的挂名怎么评判贡献？基金申请中，排名第一的挂名文章
已经有3人回复
-
越来越担心硕士毕业后没有竞争力了
已经有13人回复
-
申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请
已经有3人回复
-
闲聊
已经有6人回复
-
申博求助
已经有4人回复
孤鸿飘逸
金虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 5861.4
- 散金: 504
- 红花: 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