%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