上机实习二【DOC精选】.docVIP

上机实习二 例1已知 方程组的准确解为x* = [1，1，1]T，试用迭代法产生向量序列逐步逼近准确解。 解：（1）雅可比（Jacobi）迭代 A=[9 -1 -1;-1 10 -1;-1 -1 15]; b=[7;8;13]; x=[0;0;0];y=x; for k=1:4 for i=1:3 s=0;t=x(i);x(i)=0; for j=1:3 s=s+A(i,j)*x(j); end x(i)=t; y(i)=(b(i)-s)/A(i,i); end er=norm(x-y,1);x=y end （2）塞德尔（Seidel）迭代 A=[9 -1 -1;-1 10 -1;-1 -1 15]; b=[7;8;13];x=[0;0;0]; er=1;k=0; while er0.01 er=0;k=k+1; for i=1:3 s=0;t=x(i);x(i)=0; for j=1:3 s=s+A(i,j)*x(j); end x(i)=(b(i)-s)/A(i,i); er=max(abs(x(i)-t),er); end x end 例2 平面温度场的计算问题。 10℃、15