【matlab】 幂法 求解最大特征值

一. 算法：

1.输入矩阵A，初始向量x误差限ep,最大迭代次数N

2.置 k = 1, m1 = 0;

3.求Xr-> norm(x)   abs(Xr)=max[Xi] 1<=i<=n

4.计算    y = x/norm(u)

5. 若m1-m 小于误差限 输出结果，停止 否则 to6

6.若k<N k++ norm(x) = m1

二.程序：

A = [-6.9,,;
     -,10.1,;
     -,,-0.1];
N=;
ep=1e-;
n=length(A);
y=ones(n,);
k=;
m1=;
while k<=N
   x=A*y;
   m=max(abs(x));
   y=x/m
   if abs(m-m1)<ep
        break;
   end
   m1=m;
   k=k+;
end

final_answer_u = m;
final_answer_x = x;