数值分析实验之线性方程组的迭代求解(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实现)的主要内容,如果未能解决你的问题,请参考以下文章

数值分析实验之线性方程组的迭代求解(java实现)

数值分析实验之非线性方程求根(MATLAB实现)

数值分析实验之矩阵的LU分解及在解线性方程组中的应用(MATLAB 代码)

数值分析实验之非线性方程求根(Python 现)

怎么用matlab解方程啊?

Matlab追赶法和迭代法解线性方程组