function [l,v,s]=rpmethod(A,x0,eps)
if nargin==2
eps = 1.0e-6;
end
v = x0; %v为主特征向量
M = 5000; %迭代步数限制
m = 0;
l = 0;
for(k=1:M)
y = A*v;
m = (y'*v)/(v'*v); %m为瑞利商
if(rank(m)>1)
m=(y*v')/(v*v');
end
v = y/m;
if(abs(m - l)<eps)
l = m; %到所需精度,退出,l为主特征值
s = k;
return;
else
if(k==M)
disp('迭代步数太多,收敛速度太慢!');
l = m;
s = M;
else
l = m;
end
end
end
MATLAB实现矩阵特征值计算,案例丰富【数学建模、科学计算算法】.zip
版权申诉
145 浏览量
2023-04-21
18:08:14
上传
评论
收藏 7KB ZIP 举报
![avatar](https://profile-avatar.csdnimg.cn/6b5319dc792644f0aea19d1ec195f6a2_m0_38106923.jpg!1)
![avatar-vip](https://csdnimg.cn/release/downloadcmsfe/public/img/user-vip.1c89f3c5.png)
不脱发的程序猿
- 粉丝: 24w+
- 资源: 5757
最新资源
- 666662222277777
- 图标生成工具 Mac 版本(支持.ico和.icns)IconBuilder v1.1
- 求最大公约数(简单).cpp
- elasticsearch数据库下载、配置、使用案例
- springboot的概要介绍与分析
- C语言的概要介绍与分析
- 第一个较大的Android项目,基于Android平台的图书管理系统(Android studio).zip
- Cisco Packet Tracer 6.2 for Windows Instructor Version
- 使⽤pyIAST计算⽓体吸附选择性
- tmp_b056727e59b8123365486983f32baa9732607ec3c6137b12.pdf
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
![feedback](https://img-home.csdnimg.cn/images/20220527035711.png)
![feedback](https://img-home.csdnimg.cn/images/20220527035711.png)
![feedback-tip](https://img-home.csdnimg.cn/images/20220527035111.png)