数值分析实验之线性方程组的迭代求解(MATLAB实现)
Posted ynly
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了数值分析实验之线性方程组的迭代求解(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 ‘,k);
运行结果:
以上是关于数值分析实验之线性方程组的迭代求解(MATLAB实现)的主要内容,如果未能解决你的问题,请参考以下文章