Matrix_pow.m

function [u_normal,lambta1]=Matrix_pow(A,u)
%幂法计算矩阵最大特征值及相应特征向量
      [m,n]=size(A);
      i=1;dellab=1;miu(1)=1;
while abs(dellab)>1e-6
      u_normal(i,:)=[i-1,u'];
      i=i+1;
      v=A*u;
      [vmax,j]=max(abs(v));miu(i)=vmax*sign(v(j));
      u=v/miu(i);
      dellab=miu(i)-miu(i-1);
end
   lambta1=miu;