清华大学高等数值分析第三次作业第四题.docxVIP

矩阵阶数N=10

b-Axm

x

矩阵阶数N=50

b-A

x

矩阵阶数N=100

b-A

x

矩阵阶数N=500

b-A

x

计算的源代码如下：

N=10

O=3*ones(N,1);

P=-2*ones(N-1,1);

Q=-1*ones(N-1,1);

A=diag(O)+diag(P,1)+diag(Q,-1);

b=zeros(N,1);

b(1)=1;

b(N)=2;

e=zeros(N,1);

e(1)=1;

v=e;

x0=zeros(N,1);

r=b-A*x0;

B=zeros(N);

B(1,1)=1/3;

xm=x0+B*r;

K=linspace(1,N,N);

errorb=zeros(N,1);

errorx=zeros(N,1);

xx=ones(N,1);

errorb(1)=norm(b-A*xm);

errorx(1)=norm(xm-xx);

e(1)=0;

fori=2:N

e(i)=1;

v=[ve];

B(1:i,1:i)=invzhuigan(A(1:i,1:i));

xm=x0+B*r;

rk=b-A*xm;

errorb(i)=norm(rk);

rkx=xm-xx

errorx(i)=norm(rkx);

e(i)=0;

end

计算过程中，利用追赶法进