清华第五版数值分析第6章课件.pptVIP

第五章 线性代数方程组的数值解法 ;一. 迭代法的基本思想;向量和矩阵收敛性定义;;五 迭代法的收敛性;定理：设有线性方程组 及迭代格式 ， 若有 的某种算子范数 ，则 （1）迭代法收敛，即对任意初值 ，有 且 （2） ;; 迭代法研究的主要问题; ; 一般地，设方程组 AX=b , 通过分离变量的过程建立Jacobi迭代公式，即 ;;将系数矩阵分裂为：;;Jacobi迭代矩阵;例1：设方程组为 ;Jacobi迭代算法描述; 2.2 k=k+1 2.3 if then 2.3.1 2.3.2 goto 2 else 输出 结束。 else 2.4 输出 迭代次数太大。 3 结束 ;;;例1的Jacobi 迭代公式为;G-S迭代算法描叙; 2.1.4 if then 2.2 k=k+1 2.3 if ,输出结果，结束。 else 2.4 输出迭代次数太大。 3 结束 ;Remark: Gauss-Seidel迭代法的计算过程比Jacobi迭代法更简单。计算过程中只需用一个一维数组存放迭代向量。 Gauss-Seidel迭代不一定比Jacobi迭代收敛快。;思想;?超松弛(SOR)迭代法的分量形式：;;例1的SOR迭代公式为;迭代格式的构造; 若从(3) 中精确解出 ,则作校正 即 就是方程组（1）的准确解。 注意: 从残差方程组(3)中精确解出 的工作量跟解方程组(1)是一样的。 ;此路不通， 怎么办?!; 我们转而考虑其次：用A 的某个简单的易求逆的近似矩阵 代替（3）中的A ，得到  (注意：因为 ，故 ） 从（4）中精确解出 （设这容易办到）， 作校正 将 代入方程组（1）计算残差 再解近似残差方程组 再作近似校正 如此下去得到一般形式的校正近似的迭代框架：;1）任取初始值 ， 于是可设 为已知； 2）计算残差 3）用 （设易求逆）代替A 精确解残差方程组 4)校正近似 若达到某种预先设定的终止条件则终止，否则转2）。 ;; 推论;3) 设A是有正对角元的n阶对称矩阵，则Jacobi迭代收?? A和2D-A同为正定矩阵。 证：记 则 即 ， 从而有相同的谱半径。; 对方程组 通过调整方程的次序，建立收敛的Jacobi迭代格式和Gauss-Seidel迭代格式. 解 将第二个方程调到第一行、第三个方 程调到第二行、第一个方程调到第三行后 有同解方程组; 其系数矩阵为严格对角占优矩阵，故Jacobi和 Gauss-Seidel迭代法都一定收敛. Jacobi迭代格式为 ; Gauss-Seidel迭代格式为 ;SOR 迭代法收敛性的结论：