PM.m

% 幂法
% 注释的地方需要修改

function mifa
A = [3,1,2;1,3,4;2,4,6];           % 所求矩阵
[m,u]=pow(A,1e-4,500) ;            % 第二个参数为精度，第三个为最大迭代次数
vpa(m,7)                           % 最大特征值有效数字位数
vpa(u,7)                           % 对应的特征向量有效数字位数

function [m,u]=pow(A,ep,Nmax)
n=length(A);u=ones(n,1);k=0;m1=0;
while k<=Nmax
    v=A*u;
    [vmax,i]=max(abs(v));
    m=v(i);u=v/m;
    if abs(m-m1)<ep
        break;
    end
    m1=m;k=k+1;
end