Matrix_parallel_trans_pow.m

function [u_normal,lambta1]=Matrix_parallel_trans_pow(A,p,u)
% 幂法计算矩阵最大特征值及相应特征向量+原点平移

      [m,n]=size(A);I=eye(n);
      B=A-p*I;

      i=1;dellab=1;miu(1)=1;
while abs(dellab)>1e-6
      u_normal(i,1:n)=u;
      i=i+1;
      v=B*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+p;