线性方程组的顺序高斯消去法.doc

2. 线性方程组的顺序高斯消去法. 设 A, x, b (7.3) 考察线性方程组 Ax=b . (7.4) (1) 顺序高斯消去法的消元过程: 对方程组(7.4)的增广矩阵 作初等行变换如下: 第一步, 记, , , . 设 , 对计算 , 第二步, 设, 对计算 , , 如此继续, 第步结束, 得矩阵 (n) 然后进行回代, 即可解出: , , . (2) 上述消元求解的整个过程, 需要做多少次乘除法运算? (3) 顺序消元过程的乘除法将会产生大量的舍入误差, 往往使结果很不可靠. 例6. 方程组 的准确解是, . 请你按顺序高斯消去法编制程序计算, 取五位有效数字, 得出什么样的解? 与准确解比较一下. (4) 舍入误差的产生, 主要来自于除法, 当时, 计算过程无法进行, 当接近于0时, 产生很大的舍入误差. 为使过程继续, 且减少舍入误差, 就导致下述改进了的列主元消去法. 3. 列主元高斯消去法. (1) 计算过程与顺序高斯消去法的不同之处仅在于增加一个选行号与交换行的过程. 即对, 执行 a) 选行号, 使, 称为第k个主元. b) 交换第k行与行, 即交换与及与, . 对例6按列主元高斯消去法编制程序计算, 看看得什么结果? (3) 按列主元消去法的计算量与顺序消去法的计算量, 有什么变化? (4) 列主元消去法的矩阵解释. 我们知道, 交换矩阵A的第i行与k行, 等价于A左乘一个初等换法阵Iik, 其中Iik,是交换单位矩阵I的第i行与k行的结果. 同理, 把A的第i行的c倍加到第k行上等价于A左乘一个初等消法阵Mik (c), 其中Mik (c)是把I的第i行的c倍加到第k行上的结果. 于是求解线性方程组Ax=b的列主元消去法的前半部分, 就是A依次左乘一系列换法阵与单位下三角阵(因消去过程用到的消法阵都是单位下三角阵, 即主对角元为1的下三角阵, 它们的乘积也是单位下三角阵), 结果化为上三角阵U, 即 其中Mk均为单位下三角阵. 事实上, 换行与消元的过程可以交换, 我们可以先对A作n-1次行交换, 即依次左乘矩阵, 然后再逐个消元. 于是, 记, 可得 , 其中M与L均为单位下三角阵, U为上三角阵.一般地，我们有 定义.若矩阵G = LU, 其中L为单位下三角阵, U为上三角阵, 则称G有LU 分解. 利用矩阵的LU分解, 方程组(7.4)化为 PAx=Pb LUx =Pb . 由此即可解得x . (5) 矩阵的LDU分解. 为使非奇异矩阵的LU分解标准化, 且有对称性, 我们将U写成如下形式 即U=DU1, 其中D是非奇异对角阵, U1是单位上三角阵. 于是, 若矩阵A=LDU,其中L与U分别为单位下三角阵与单位上三角阵, D为非奇异对角阵, 则称A有LDU分解. 定理7.1. 设A=(aij)是n阶方阵, 令 , 表示A的各阶主子矩阵. 则A有唯一LDU分解, 即A= LDU, 当且仅当Ak均为非奇异矩阵 证. 设A=LDU, 其中L, D, U如前述, 则L, D, U均为非奇异阵. 因Ak=LkDkUk,, 其中Lk, Dk, Uk分别为L, D, U的k阶主子矩阵, 均非奇异, 故Ak非奇异, 对. 反之, 设Ak均非奇异, . 我们对归纳证明Ak有LDU分解. 首先A1显然有LDU分解. 今设Ak有LDU分解: Ak=LkDkUk,, 令 其中x, y和a均为已知. 由 , , 即可解出唯一的u, v与d: , , . 故Ak+1有LDU分解. 这就归纳证明了A=An有LDU分解. 最后, 我们证明LDU分解的唯一性. 设A有LDU分解 A=L1D1U1= L2D2U2 . 因A非奇异, 故Li, Di, Ui均非奇异, i=1, 2. 得 , 左端为单位下三角阵, 右端则是上三角阵, 可见应是n阶单位阵I, 即,便推出L2=L1, 同理可证U2=U1, 并且D2=D1. (6) 如何利用矩阵A的LDU分解去求解方程组Ax=b . 4. 线性方程组迭代解法的一般形式 (1) 迭代法的一般形式. 用迭代法求线性方程组Ax=b的解的基本思想与§4介绍的迭代法求方程根是类似的. 如果Ax=b等价于方程组