一. 算法:
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;