A=[1,0.8,0.8;0.8,1,0.8;0.8,0.8,1];
b=[3;2;1];
经过计算,得到题目 2 的迭代阵 B 的谱半径 R=0.7155<1,所以能够得到收
敛解。经过 12 次迭代后,得到该线性方程组的解为:
x=[5.7692;0.7692;-4.2308]
题目 3 在 Command Window 中的程序:
A=[1,3;-7,1];
b=[4;6];
x0=[0;0];
经过计算,得到题目 3 的迭代阵 B 的谱半径 R=21>1,所以该线性方程组不
能得到用 Gauss-Seidel 迭代方法求得收敛解。
通过使用上面三个例题求解的过程发现,在 MATLAB中编写的 M 文件能够
很好地求解线性方程组,Gauss 消元法能够得到上面三个例题的全部收敛解,而
Gauss-Seidel 虽然无法求解到题目 3 的收敛解,但是整个 M 文件程序的使用性能
还是非常的不错,能够真实反映出线性方程组系数矩阵内在的一些性质。