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
版权申诉
144 浏览量
2023-04-21
18:08:14
上传
评论
收藏 7KB ZIP 举报
不脱发的程序猿
- 粉丝: 24w+
- 资源: 5763
最新资源
- 高分项目,PID-小车类-PID算法控制小车直线行驶(制作步骤+程序+PID库)
- Solution1.rarSolution1.rarSolution1.rar
- 高分项目,PID-增量式PID和位置式PID算法实现和PID库(51单片机+STM32+arduino完成)
- 1414114147141
- 高分项目,PID-倒立摆-旋转倒立摆源码+PID库下载
- culane figure
- electron+vue3 开发引入软件平台窗口,可通过配置地址的打开配置的服务地址网站或系统
- 简鹿水印助手 1.0.2.0 离线包.zip
- centerpoint ganet hbm model
- CMMI课件,对软件能力成熟度感兴趣的、从事质量管理的可以关注下
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈