§23直接三角分解法 Gaussi消去法消元过程的矩阵描述 n A,b(1)= n1 n2 nn n 第一次消元相当于用初等矩阵 左乘(4(①),b(1) ),其中 i=2.3
一、Gauss消去法消元过程的矩阵描述 ( , ) (1) (1) A b = (1) (1) (1) 2 (1) 1 (1) 2 (1) 2 (1) 2 2 (1) 2 1 (1) 1 (1) 1 (1) 1 2 (1) 1 1 n n n n n n n a a a b a a a b a a a b 第一次消元相当于用初等矩阵 − − = 1 1 1 1 21 1 n l l L ( , ) (1) (1) 左乘 A b ，其中 i n a a l i i 2,3, , (1) 11 (1) 1 1 = = §2-3 直接三角分解法
第k次消元相当于用矩阵 k+1,k nk (k) 左乘(6),b),其中k=i=k+1,3…,n lk 因此 1 L2·L1(4),b()=(40,b0 从而Ln1…L2·L1·A=A
第k次消元相当于用矩阵 − − = + 1 1 1 1 , 1, n k k k k l l L ( , ) (k ) (k ) 左乘 A b ，其中 i k n a a l k k k i k i k 1,3, , ( ) 1 ( ) = = + ( , ) (1) (1) 1 L A b L2 Ln−1 ( , ) (n) (n) 因此 = A b 从而 Ln−1 L1 A (n) L2 = A
故A=L1L2…Ln L=1,L2…L, 1,1 1,21n-1,2 n,1 n.2 n.n-1 为单位下三角矩阵 12 U=A(n)-= 27为上三角矩阵
1 ( ) 1 1 2 1 1 n L L L n A − − − − 故 A = = LU = = − − − − −− − − 1 1 1 1 1 . . . ,1 ,2 ,3 , 1 1,1 1,2 1,2 3 1 3 2 2 1 11 12 11 n n n n n n n n n l l l l l l l l l l L L L L 即 为单位下三角矩阵 = ( ) (2) 2 (2) 2 2 ( 1 ) 1 ( 1 ) 1 2 ( 1 ) 1 1 n n nnn a a a a a a ( n ) U = A 为上三角矩阵
、矩阵的三角分解 定义1设A为n阶方阵,若存在下三角方阵L和上三角 方阵U,使得A=LU,则称方阵A有三角分解或LU分解。特 别,若L为单位下三角方阵,则称它为杜利特尔 Doolittle) 分解;若U为单位上三角方阵,则称它为克劳特( Crout)分 解 定理1若n阶方阵4=(an)m的顺序主子式D≠0, k=1,2,…,n则存在n阶单位下三角方阵L和上三角方阵U, 使得A=LU,即A有杜利特尔分解. 且detA=detU=∏a
二、矩阵的三角分解 定义1 设A为n阶方阵，若存在下三角方阵L和上三角 方阵U，使得A=LU，则称方阵A有三角分解或LU分解。特 别，若L为单位下三角方阵，则称它为杜利特尔(Doolittle) 分解；若U为单位上三角方阵，则称它为克劳特(Crout)分 解. 定理1 则存在n阶单位下三角方阵L和上三角方阵U， 使得A=LU，即A有杜利特尔分解. = ( ) 0, 若n阶方阵A ai j nn 的顺序主子式Dk 且 det A = detU = = n i i aii 1 ( ) k = 1,2, ,n
11 因此A=a 1 1 根据矩阵的乘法原理,A的第一行元素a1,为 j=1 A的第r行元素主对角线以右元素an(j=r…,n)为 n ∑ 丿∫
根据矩阵的乘法原理, A的第一行元素a1 j 为 a1 j = u1 j j = 1,2, ,n A的第r行元素主对角线以右元素ar j ( j = r, ,n)为 = = r k rj rkukj a l 1 r = 1,2, ,n j = r, ,n = 1 1 1 1 1 n nr r l l l nn r r r n r n u u u u u u 11 1 1 = n nr nn r r r r n r n a a a a a a a a a A 1 1 11 1 1 因此