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

%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