数值分析 幂法求矩阵A按模最大的特征值和相应的特征向量

2019-04-13 12:25发布

%function [x,y]=jingmi(A,u0,ep)%x特征值,y特征向量,u0初始向量,ep允许误差。
%幂法求矩阵A按模最大的特征值和相应的特征向量 %参考数值分析第四版 颜庆津 P49
%2015.11.7  密密编写  (*^__^*) 嘻嘻…… function [x,y]=jingmi(A,u0,ep)%x特征值,y特征向量 %xx=1:10;i=1;     m=sqrt(u0'*u0);     y=u0/m;     u1=A*y;     x0=y'*u1; %xx(i)=x0; %i=i+1;     err=1;     u0=u1; while(err>ep)     m=sqrt(u0'*u0);     y=u0/m;     u1=A*y;     x=y'*u1;     err=abs(x-x0)/abs(x);     u0=u1;     x0=x; %xx(i)=x0; %i=i+1; end
end