matlab迭代算法求解四元线性方程组,数值分析实验之线性方程组的迭代求解(MATLAB实现)...

详细实验指导见上上一篇，此处只写内容啦

实验内容： 求解如下4元线性方程组的近似解。

Gauss-Seidel迭代

A=[10,-1,2,0;-1,11,-1,3;2,-1,10,-1;0,3,-1,8];

b=[6;25;-11;15];

x=zeros(4,1);

temp=zeros(4,1);

x0=[1;2;-1;1];

for k=1:100

for l=1:4

temp(l,1)=abs(x(l,1)-x0(l,1));

end

if max(max(temp))>0.00001

for i=1:4

x(i,1)=1/A(i,i)*b(i,1);

for j=1:i-1

if i-1~=0

x(i,1)=x(i,1)-1/A(i,i)*A(i,j)*x(j,1);

end

end

for j=i+1:4

if j<=4

x(i,1)=x(i,1)-1/A(i,i)*A(i,j)*x(j,1);

end

end

end

else

break

end

end

fprintf(‘迭代次数为%g\n‘,k);

运行结果：