itermulteight3.m

来自「数值分析方法,主要里面涵盖了计算方法方面matlab程序！」· M 代码 · 共 28 行

clear;
clc;
%function iterMult(A,u,KM,eps)     %u为初始列向量，KM是最大迭代次数，eps是容许误差。
A=xlsread('data8_3.xls');
eps=1e-5;
KM=50;
N=size(A,1);
u=ones(N,1);
k=0;r1=0;
while k<KM
    k=k+1;
    v=A*u;
    [r2,p]=max(abs(v));
    r2=v(p);
    u=v/r2;
    if(abs(r1-r2)<eps)
        break;
    end
    r1=r2;
end
if k>=KM
    fprintf('The Method is disconvergent\n');
else
    fprintf('r=%f\n',r2);
    for i=1:N
        fprintf('x[%d]=%f\n',i,u(i));
    end
end