使用Matlab实现：Jacobi、Gauss-Seidel迭代

例题

方程组$\{\begin{array}{l}5x_1+2x_2+x_3=-12\\ -x_1+4x_2+2x_3=20\\ 2x_1-3x_2+10x_3=3\end{array}$ 求解，当$max|x_i^{(k+1)}-x_i^{(k)}|\le 10^{-5}$ 时候迭代终止。

以下解答过程，上标表示迭代次数，下标表示序号。

Jacobi迭代

定义变量：
$$\begin{gathered} D=diag(a_{11},a_{22},...,a_{nn}), \\ L=\left[\begin{array}{cccccc}0& & & & & \\ -a_{21}& 0& & & & \\ ...& & & & & \\ -a_{i1}& ...& -a_{i,i-1}& 0& & \\ ...& & & & & \\ -a_{n1}& ...& -a_{n,i-1}& ...& -a_{n,n-1}& 0\end{array}\right], \\ U=\left[\begin{array}{cccccc}0& -a_{12}& ...& -a_{1,i}& ...& -a_{1,n}\\ ...& & & & & \\ & & 0& -a_{i-1,i}& ...& -a_{i-1,n}\\ ...& & & & & \\ & & & & 0& -a_{n-1,n}\\ ...& & & & & \\ & & & & & 0\end{array}\right] \end{gathered}$$

其矩阵迭代形式为：
$$\begin{gathered} x^{(k+1)}=B_J\cdot x^{(k)}+f_J \\ {B}_J=D^{-1}\cdot (L+U),f_J=D^{-1}\cdot b \end{gathered}$$。

写出分量形式：

$\{\begin{array}{l}{x}_1^{(k+1)}=\frac{1}{5}(-12-2x_2^{(k)}-x_3^{(k)})\\ x_2^{(k+1)}=\frac{1}{4}(20+x_1^{(k)}-2x_3^{(k)})\\ x_3^{(k+1)}=\frac{1}{10}(3-2x_1^{(k)}+3x_2^{(k)})\end{array}$

写成矩阵形式：

$\left[\begin{array}{c}{x}_1^{(k+1)}\\ x_2^{(k+1)}\\ x_3^{(k+1)}\end{array}\right]=\left[\begin{array}{ccc}0& -0.4& -0.2\\ 0.25& 0& -0.5\\ -0.2& 0.3& 0\end{array}\right]\cdot \left[\begin{array}{c}{x}_1^{(k)}\\ x_2^{(k)}\\ x_3^{(k)}\end{array}\right]+\left[\begin{array}{c}-2.4\\ 5\\ 0.3\end{array}\right]$

取初始向量：$x^0=(0,0,0{)}^T$，依次按照上式进行迭代。使用Matlab进行编程求解。

a=[0,-0.4,-0.2;0.25,0,-0.5;-0.2,0.3,0];
b = [-2.4;5;0.3];
x = [0;0;0];
xx = a * x + b;
i = 0;
while norm(x - xx, inf) >= 1e-5
	x = xx;
	xx = a * x + b;
	i = i +1;
end

以上代码，最终 $x=x^i,xx=x^{(i+1)}$ ，最终迭代次数位 $i+1$ 次，如果你需要看到更长的小数位置，可以使用以下Matlab代码，表示使用15位浮点或定点数。

format long g

运行结果为：

即精确解为 $x=(-4,3,2{)}^T$ 。

Gauss-Seidel迭代

其矩阵迭代形式为：
$$\begin{gathered} x^{(k+1)}=B_G\cdot x^{(k)}+f_G \\ {B}_G=(D-L{)}^{-1}\cdot U,f_G=(D-L{)}^{-1}\cdot b \end{gathered}$$

使用Matlab编程求解：

d = [5,0,0;0,4,0;0,0,10];
l = [0,0,0;1,0,0;-2,3,0];
u = [0,-2,-1;0,0,-2;0,0,0];
b = [-12;20;3];
t = inv(d - l);
bg = t * u;
fg = t * b;
x = [0;0;0];
xx = [-2.4;4.4;2.1];
i = 0;
while norm(x - xx, inf) >= 1e-5
	x = xx;
	xx = bg * x + fg;
	i = i +1;
end

运行结果为：

同样求得精确解为 $x=(-4,3,2{)}^T$ 。