计算方法第5.6章.ppt

定理11 设 , 如果A为严格对角占 优阵，则解 的Jacobi方法，G-S迭 代法都收敛，且G-S迭代法收敛比Jacobi 方法为快。 证明 （1）首先证明解 的Jacobi方法收 敛。由假设有 (8.15) 解 的Jacobi方法迭代矩阵为: 故由定理9知解 的Jacobi方法收敛。 （2）由设 有唯一解记为 ，即 且由（8.15）有 其中 或 (8.16) G—S迭代公式：: (8.17) 引进误差向量： 且记 由（8.17）减去（8.16）即得G—S迭代近似解的误差递推公式： （8.18） 记 由设有 即 由（8.18）得 （8.19） （8.20） 记 在（8.18）中取 , 则有 即 其中 且由于 ，所以 。 利用递推关系（8.20），则有 故解 的G—S迭代法收敛。 所以， 成立。（8.21）说明解方程组 （其中设A为严格对角占优阵）G—S方法收敛比Jacobi方法为快。 （8.21） （3） 说明 事实上，对每一个 有 定理12 （1）设 其中A为对称正定阵； （2） ； 则解 的SOR方法收敛。 例13 设有方程组 试考察用Jacobi及G—S迭代法解此方程组的收敛性。 解 由于 为严格对角占优阵，于是由定理11可知解 的Jacobi迭代、G—S迭代均收敛。 §9 解非线性方程组的迭代法 在科学技术领域里常常提出求解非线性方程组问题，例如，用非线性函数拟合实验数据问题，非线性网络问题，用差分法求解非线性微分方程问题等。 设有非线性方程组 （9.1） 于是方程组（8.3）等价于方程组 （8.5） (8.6) 其中M为可选择的一个非奇异矩阵，应选择M使 容易求解。 对应于方程(8.5) 可构造一个迭代过程： 8.1 雅可比迭代法 选取M=D，于是 ，方程（8.3）转化为等价方程组 (8.7) J 称为 J a c o b i 迭代法的迭代矩阵。 于是得到雅可比迭代公式： Jacobi迭代公式的分量形式： 记 为第 k 次近似，则（8.7）式可写出为 或 （8.8） Jacobi 迭代法公式简单,由公式(8.7)或 (8.8)可知,每次迭代只需计算一次矩阵与向量的乘法，例9的迭代法就是解 的 Jacobi 迭代法。电算时Jacobi方法需要两组工作单元来保存 及 且可用 来控制迭代终止。由迭代法计算公式可知，迭代法一个重要特点是计算过程中原来矩阵A数据始终不变。 8.2 高斯-塞德尔迭代法 在（8.5）式中选取 M = D – L （下三角矩阵）。于是， N = M – A = U 方程（8.3）转化为等价方程组 于是得到高斯-塞德尔（G – S ）迭代公式： （8