数值分析教教案21.docVIP

5.2.2 Gauss－Seidel（高斯-赛德尔迭代法） 1. Gauss－Seidel迭代原理 Jacobi迭代公式（5-9）用方程组表示为 （5-13） 因此，在Jacobi迭代法的计算过程中，要同时保留两个近似解向量和。如果把迭代公式改成以下形式 （5-14） 即每算出新近似解的一个分量，再算下一个分量时，用新分量代替老分量进行计算。这样，在整个计算过程中，只需要个单元存贮近似解分量。 选取初始向量，用迭代公式（5-14）产生近似解序列，这种方法叫Gauss－Seidel迭代法，式（5-14）为Gauss－Seidel迭代法的计算公式。 公式（5-14）用矩阵表示为 （5-15） 其中 移项可得 因为，故存在，上式可改写成 （5-16） 如果用矩阵来表法，记 则 于是 （5-17） 将式（5-17）代入式（5-16） 这是Gauss－Seidel迭代公式的矩阵表示，式中矩阵为迭代矩阵。 2. Gauss－Seidel算法 （1）输入,维数,误差限，最大容许迭代次数。 （2）置。 （3）计算： （4）若，输出，停机；否则转（5）。 （5）若，置，转（3）；否则输出失败信息，停机。 用Gauss－Seidel迭代法解例5-4。仍取，按式（5-14）计算得 如此继续下去，计算结果见表5-2。 0 0.0000 0.0000 0.0000 1 7.20