幂法求特征值

%幂法
%用于计算一个可对角化的n*n矩阵的绝对值最大的特征值
%返回值m为绝对值最大特征值，而x为该特征值对应的特征向量
function [m,x]=PM(A)
    n=max(size(A));
    u=ones(n,1);
    y=ones(n,1);%初始化，初始值是多少不重要
    beta1=0;

    eta=norm(u,2);
    y=u./eta;
    u=A*y;
    beta2=y'*u;

    while abs((beta2-beta1)/beta1)>1e-12
        beta1=beta2;
        eta=norm(u,2);
        y=u./eta;
        u=A*y;
        beta2=y'*u;
    end
    m=beta2;
    x=y;
end