第四章线性代数方程组的迭代解法.ppt

;4.2 迭代法的基本思想 迭代法的基本思想是将线性方程组转化为便于迭代的等价方程组，对任选一组初始值 ，按某种计算规则，不断地 对所得到的值进行修正，最终获得满足精度要求的方程组的近似解。 ;设 非奇异， ，则线性方程组 有惟一解 ，经过变换构造出一个等价同解方程组 将上式改写成迭代式; 如果 存在极限 则称迭代法是收敛的，否则就是发散的。 收敛时，在迭代公式 中当 时， , 则 , 故 是方程组 的解。 对于给定的方程组可以构造各种迭代公式， 但并非全部收敛 。;例4.1 用迭代法求解线性方程组 ;4.3 雅可比（Jacobi）迭代法 4.3.1雅可比迭代法算法构造 ;取初始向量 进行迭代, 可以逐步得出一个近似解的序列： （k=1, 2, …） 直到求得的近似解能达到预先要求的精度， 则迭代过程终止，以最后得到的近似解作为线 性方程组的解。 当迭代到第10次有 计算结果表明，此迭代过程收敛于方程组的精 确解x*= (3, 2, 1)T。 ;考察一般的方程组，将n元线性方程组 ;4.3.2 雅可比迭代法的矩阵表示 设方程组 的系数矩阵A非奇异，且主对 角元素 ，则可将A分裂成 ;;;雅可比迭代矩阵表示法，主要是用来讨论其收敛性，实际计算中，要用雅可比迭代法公式的分量形式。即 ;4.3.3 雅可比迭代法的算法实现 ;4.4 高斯-塞德尔（Gauss-Seidel）迭代法 4.4.1 高斯-塞德尔迭代法的基本思想 在Jacobi迭代法中，每次迭代只用到前一次的迭代值，若每次迭代充分利用当前最新的迭代值，即在求 时用新分量 代替 旧分量 , 就得到高斯-赛德尔迭代法。 其迭代法格式为： ;例4.3 用Gauss?Seidel 迭代格式解方程组 ;4.4.2 Gauss—Seidel 迭代法的矩阵表示 将A分裂成A =L+D+U，则 等价于 ( L+D+U )x = b 于是,则高斯—塞德尔迭代过程 ;4.5 迭代法的收敛性 我们知道, 对于给定的方程组可以构造成简单迭代公式、???可比迭代公式、高斯-塞德尔迭代公式和超松弛迭代公式，但并非一定收敛。现在分析它们的收敛性。 对于方程组 经过等价变换构造出的等价方程组 ; 由此定理可知，不论是雅可比迭代法、还是高斯—塞德尔迭代法，它们收敛的充要条件是其迭代矩阵的谱半径 。 ;定理4.2 (迭代法收敛的充分条件) 若迭代矩阵B的一种范数 ,则迭代公式 收敛。 由定理知，当 时，其值越小，迭代收敛越快，在程序设计中通常用相邻两次迭代 （ε为给定的精度要求）作为控制迭代结束的条件。 ; 三种常用的向量范数：;三种常用的矩阵范数：;例4.5 已知线性方程组 ;⑵ 将系数矩阵分解 ;定理4.3 设n阶方阵 为对角占优阵, 则 A非奇异。 证: 因A为对角占优阵, 其主对角元素的绝对值大 于同行其它元素绝对值之和, 且主对角元素 全不为0, 故对角阵 为非奇异。 作矩阵 ;;定理4.4 对角占优线性方程组 的雅可比 迭代公式和高斯-赛德尔迭代公式均收敛。 证: 雅可比迭代公式的迭代矩阵为 ;;例4.6 设求解线性方程组 的雅可比迭代 ;例4.7 设 ,证明, 求解方程组 ;例4.7 设 ,证明, 求解方程组 ;;;;;;;;;;;;;;; 迭代法不存在误差累积问题。使用迭代法的 关键问题是其收敛性与收敛速度，收敛性与迭代 初