直接三角解法.ppt

第四章方程组的直接解法 4.2 直接三角分解法 4.2.3 平方根法 4.2.1 一般矩阵的直接三角分解法 4.2.2 三对角方程组的追赶法 第四章方程组的直接解法 4.2 直接三角分解法 4.2.1 一般矩阵的直接三角分解法 本节讨论矩阵 A 的三角分解法的直接计算以及直接利用 A 的三角分 解式来求解方程组。 1. 不选主元的三角分解法 ), ( ), ( ), ( ij ij ij u U l L a A ? ? ? 设 A=LU ，记 其中 L 为单位下三角阵， U 为上三角阵。我们可直接给出 L 和 U 的元素的计算公式。 由 A 的第 1 行和第 1 列可计算出 U 的第 1 行和 L 的第 1 列，即 （ 4.2.1 ） （ 4.2.2 ） 如果 U 的第 1 至 k-1 列和 L 的第 1 至 k-1 列已经算出，则由 , , , 1 , , 1 n k k j u l a k r rj kr kj ? ? ? ? ? ? 1 1 1 1 11 , 1,2, , , , 2,3, , . j j k k u a j n a l k n u ? ? ? ? 第四章方程组的直接解法 可得 U 的第 k 行元素 同理，由 u kj =a kj - ,j =k,k+1, · · · ,n 。 （ 4.2.3 ） ? ? ? 1 1 k r rj kr u l a kj = ,i=k+1,k+2,· · · ,n, ? ? k r rk ir u l 1 可得 L 的第 k 列元素 交替使用（ 4.2.3 ）和 （ 4.2.4 ），就能逐次计算出 U （按行）和 L （按列）的全部元素，而且可以把它们存放在矩阵 A 对应的位置上 （ L 的对角线元素不必存放）。这就完成了 A 的 LU 分解。 ? ? k r rk ir u l 1 l ik =(a ik - )/ u kk ,i =k+1,k+2, · · · ,n 。 由（ 4.2.1 ） - （ 4.2.4 ）求得 L 和 U 后，解方程组 Ax=b 接化接为求 解 LUx=b ，若记 Ux=y ，则有 Ly=b 。于是可分两部解方程组 LUx=b ， 只要琢次向前代入的方法即可求得 y 。第二步求解 Ux=y ，只要琢次 第四章方程组的直接解法 用向后回代的方法即可求得 x 。设 x= （ x1 ， x2 ， · · · xn) T, y= （ y1, y2, · · · yn) T,b= （ b1 ， b2 ， · · · bn) T, 则有计算公式 ? ? ? ? ? ? ? ? ? ? ? ? 1 1 1 1 ,..., 2 , 1 , y i r r ir i i n i y l b y b ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? n i r ii r ir i i nn n n n i u x u y x u y 1 n 1 ,..., 2 , 1 , / / x （ 4.2.5 ） （ 4.2.6 ） 以上解方程组的计算与顺序 Gauss 消去法相当。如果有一系列方 程组，其系数距阵都是相同的，右端向量 b 不同，则只须进行一次 LU 分解计算。上述解方程的方法称为 LU 分解法 ，也称 Doolittle 方法 。 例 4.5 用 LU 分解法求解 第四章方程组的直接解法 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 5 5 1 6 3 1 0 1 1 4 1 1 0 1 4 2 1 1 2 6 4 3 2 1 x x x x 解 由（ 4.2.1 ） - （ 4.2.4 ）计算可得 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?