矩阵三角分解法
3.2矩阵三角分解,平方根
2014-12-29
北京信息科技大学
4
2 3 1 0 2 3 1 0 2 3 例如A= 4 1 2 1 0 5 2 3 0 5 / 3
这里有A的两种不同的三角分解，类似可举出很多，一般， 若A=LU是一个三角分解，任取与A同阶的非奇异对角矩阵D， 则
xn yn / unn xi ( yi i n 1,
k i 1
u
n
ik
xk ) / uii
, 2,1
2014-12-29
北京信息科技大学
23
例 用杜利特尔分解法求解方程组
2 x1 2 x2 3 x3 3 4 x1 7 x2 7 x3 1 2 x 4 x 5 x 7 1 2 3
， 则矩 , n 1 ）
U 为上三角阵。 其中 L 为单位下三角阵，
其中
A LU
1 l21 1 L l31 l32 l n1 ln 2
1 ln ,n 1
u11 u12 u22 ，U 1
u1n u2 n unn
矩阵三角分解法
矩阵直接三角分解法是高斯消去法的变形方法。 高斯消去法 有多种变形，有的是高斯消去法的改进，有的是用于某种特殊系 数矩阵的化简。 分解原理 1．概述 高斯消去法解线性方程组先消元，然后再回代。当用矩阵描 述时，是对系数矩阵分解为一个上三角阵和一个下三角阵的乘 积，即 LU 分解。因此，高斯消去法与矩阵的 LU 分解是一致的。
2014-12-29 北京信息科技大学 15
紧凑格式
2 2 3 A 4 7 7 5 2 4 (2) 2 (4) 2 (-2) -1 (2) 2 (7) 3 (4) 2 (3) 3 (7) 1 (5) 6
矩阵分解——精选推荐
矩阵分解矩阵分解矩阵分解是将矩阵拆解为数个矩阵的乘积，可分为三⾓分解、满秩分解、QR分解、Jordan分解和SVD（奇异值）分解等，常见的有三种.矩阵的三⾓分解、正交三⾓分解、满秩分解将矩阵分解为形式⽐较简单或性质⽐较熟悉的⼀些矩阵的乘积,这些分解式能够明显地反映出原矩阵的许多数值特征,如矩阵的秩、⾏列式、特征值及奇异值等. 另⼀⽅⾯, 构造分解式的⽅法和过程也能够为某些数值计算⽅法的建⽴提供了理论依据. 接下来就讨论⼀下矩阵的三⾓分解.1 矩阵的三⾓分解1.1 矩阵的三⾓分解基本概念与定理定义1.1[]5设m n∈和上三⾓矩L C?A C?∈,如果存在下三⾓矩阵m n阵n m∈, 使得A=LU, 则称A可作三⾓分解或LU分解.U C?定义1.2设A为对称正定矩阵, D为⾏列式不为零的任意对⾓矩阵,则T=成⽴：A A=, U为⼀个单位上三⾓矩阵, 且有A LDU1) 如果L是单位下三⾓矩阵, D是对⾓矩阵, U是单位上三⾓矩阵, 则称分解D=为LD U分解.A L U2) 如果L=LD是下三⾓矩阵, ⽽U是单位上三⾓矩阵, 则称三⾓分解A LUCrout分解;= 为克劳特()3) 如果U DU是单位下三⾓矩阵, U 为上三⾓矩阵, 则称三⾓=分解A LUDoolittle分解;= 为杜利特()U --=== , 称为不带平⽅根的乔累斯基()Cholesky 分解;5) 如果12L D L = , 12D U U= , 则1122A LD U LD D U LU=== , 由于T UL = , 则T A LL= , 称为带平⽅根的乔累斯基()Cholesky 分解. 定理 1.1 n阶⾮奇异矩阵A可作三⾓分解的充要条件是k 0A ≠()1,2,,1k n =- ,这⾥A k为A 的k 阶顺序主⼦阵, 以下同.证明必要性. 设⾮奇异矩阵A 有三⾓分解A L U=, 将其写成分块形式k12k122122212222A L 0U =A A 0U kA U L L这⾥A k ,k L 和k U 分别为A, L和U 的k 阶顺序主⼦阵. ⾸先由0⽽L 0k ≠,U 0k ≠; 因此A =L U0kkk ≠()1,2,,1k n =-.充分性. 对阶数n 作数学归纳法. 当n=1时, 1A =（11a ）=（1）（11a ）,结论成⽴. 设对n k =结论成⽴, 即k =k k A L U , 其中k L 和k U 分别是下三⾓矩阵和上三⾓矩阵. 若k 0A ≠,则由kA =L k k U 易知L k 和k U 可逆. 现证当1n k =+时结论也成⽴, 事实上-1k k k k1TT 1T 1-1k+1,1k 1,1k k k A c 0c A =10c kkk T kk k k k k L U L r a r U a r U L +--+++??= ? ?-.由归纳法原理知A 可作三⾓分解.定理 1.1 给出了⾮奇异矩阵可作三⾓分解的充要条件, 由于不满⾜定理1.1的条件, 所以它不能作三⾓分解. 但110000110011211011202A ?????????? ?===.上例表明对于奇异矩阵,它还能作三⾓分解未必要满⾜定理1.1的条件.⾸先指出,⼀个⽅阵的三⾓分解不是唯⼀的, 从上⾯定义来看,杜利特分解与克劳特分解就是两种不同的三⾓分解,其实,⽅阵的三⾓分解有⽆穷多, 这是因为如果D 是⾏列式不为零的任意对⾓矩阵, 有1()()A LU C D D U LU-== ,其中,LU 也分别是下、上三⾓矩阵, 从⽽A LU = 也使A 的⼀个三⾓分解. 因D 的任意性, 所以三⾓分解不唯⼀. 这就是A 的分解式不唯⼀性问题, 需规范化三⾓分解.定理 1.2 （LD U 基本定理）设A 为n 阶⽅阵,则A 可以唯⼀地分解为A =LD U(1.1)的充分必要条件是A 的前1n -个顺序主⼦式k 0A ≠()1,2,,1k n =- .其中L,U分别是单位下、上三⾓矩阵, D是对⾓矩阵D=diag ()12,,,n d d d ,1k k k A d A -=()1,2,,kn = , 01A =.证明充分性. 若k 0A ≠()1,2,,1k n =- , 则由定理1.1, 即实现⼀个杜利特分解A LU= , 其中L 为单位下三⾓矩阵, U 为上三⾓矩阵,记1112122==()()()()()()1111112122222n n n nn a a a a a a ??=()n A , 因为()u 0i ii ii a ≡≠()1,2,,1i n =- .下⾯分两种情况讨论：1) 若A ⾮奇异,由式(1)有n ?=()()() 121122n nn a a a =A ≠, 所以()n nn nna u =≠,这时令()()()()121122diag n nn D a a a = , 则() ()()1121122111,,,n nn D diag a a a -??= ?.LD D U LDU -=== (1.2)是A 的⼀个LD U 分解.2)若A 奇异,则()u 0i iiii a ≡=,此时令()()()12111221,1(,,,,0)n n n D diag a a a ---= ,()()()()121n-111221,1,,,n n n D diag a a a ---= , α=()1n1u,,,Tn u n - ,则10n T UU α-??≡ =1111110=DU 0001n n n n T T U D U D α------,因此不论哪种情况, 只要k0A ≠()1,2,,1k n =- , 总存在⼀个LD U分解式(1.1),1a kk k kk k A d A -==()1,2,,1kn =- ,01A =.均⾮奇异.若还存在另⼀个LD U 分解111A L D U =, 这⾥1L ,1D , 1U 也⾮奇异,于是有111L D U L D U =(1.3)上式两端左乘以11L -以及右乘以1U -和1D -, 得111111L L D U U D---=, (1.4)但式(1.4)左端是单位下三⾓矩阵, 右端是单位上三⾓矩阵, 所以都应该是单位阵, 因此1LL I-=,1111D U UDI--=,即1L L =,111--=. 由后⼀个等式类似地可得11U UI-=,11D D I-=,即有1U U=,1D D=.2) 若A 奇异, 则式(1.3)可写成分块形式1111100001000110001T T T T T L D U L D U ααββ= ? ? ? ? ? ???????????, 其中1L, 1L 是1n -阶单位下三⾓阵; U , 1U 是1n -阶上三⾓阵; D,1D 是1n -阶对⾓阵; α, 1α,β, 1β是1n -维列向量. 由此得出111111=D U D DUD ααββαββα???? ? ???, 其中1L, 1D , 1U 和L ,D, U均⾮奇异, 类似于前⾯的推理, 可得1L =L ,1D =D , 1U =U ,1=αα,T T1=ββ.必要性. 假定A 有⼀个唯⼀的LD U 分解, 写成分块的形式便是1111A 00=0101n n n n T T nn n x D L U ya d αβ----,(1.5)其中1n L -,1D n -, 1n U -, 1n A -分别是L,A的1n -阶顺序主⼦矩阵;x , y, α,β为1n -维列向量. 由式(1.5)有下⾯的矩阵⽅程:1111n n n n A L D U ----=, (1.6)11TTn n yD U β--=,(1.7)11n n x L D α--=, (1.8)1Tnn n na D d βα-=+. (1.9)否则, 若10n A -=, 则由式(1.6)有111110n n n n n A L D U D -----===.于是有1110n n n L D D ---==, 即11n n L D --奇异. 那么对于⾮其次线性⽅程组(1.8)有⽆穷多⾮零解, 不妨设有α', 使11n n L D x α--'=, ⽽α'=α.同理, 因11n n D U --奇异, ()1111TTT n n n n L D U D ----=也奇异,故有ββ'≠, 使11TTn n U D yβ--=, 或11TTn n D U yn nn n d a D βα-'''=-, 则有1111000101n n n n T T nn nA x D L U y a d αβ----'= ? ? ? ?'',这与A 的LD U 分解的唯⼀性⽭盾, 因此10n A -≠.考察1n -阶顺序主⼦矩阵1n A -由式(1.6)写成分块形式, 同样有2222n n n n A L D U ----=. 由于10n D -≠, 所以20n D -≠, 可得222220n n n n n A L D U D -----==≠, 从⽽20n A -≠. 依此类推可得0k A ≠()1,2,,1k n =- .综上所述, 定理证明完毕.推论 1[]3 设A 是n 阶⽅阵, 则A 可惟⼀进⾏杜利特分解的充分必要条件是A 的前1n -个顺序主⼦式11110k k k kka a A a a =≠,1,2,,1k n =- , 其中L 为单位上三⾓矩阵, 即有11121212223132121111n nnn n n n n u u u l u u l l A u l l l -=并且若A 为⾮奇异矩阵, 则充要条件可换为: A的各阶顺序主⼦式全不为零, 即:0k A ≠,1,2,,k n = .推论 2[]3 n 阶⽅阵A 可惟⼀地进⾏克劳特分解111212122212111n nn n nnl u u ll u A LUl l l==的充要条件为11110k k k kka a A a a =≠, 1,2,,1k n =- .若A 为奇异矩阵, 则0nn l =, 若A 为⾮奇异矩阵, 则充要条件也可换为0k A ≠, 1,2,,k n = .定理 1.3[]3 设A 为对称正定矩阵, 则A 可惟⼀地分解为T A LDL =, 其中L 为下三⾓矩阵, D 为对⾓矩阵, 且对⾓元素是L 对⾓线元素的倒数. 即2212n n nnl l l L l l l ?? ?=, 1122111nn l l D l ?? ? ? ? ?=. 其中11/j ijij ik jk kkk l a l l l -==-∑,1,2,,ni = , 1,2,,j i = .。
矩阵论中不同形式的三角分解的优缺点
矩阵论是数学领域中的一个重要分支，它研究的是矩阵和线性方程组的理论和方法。
在矩阵论中，三角分解是一种常见的矩阵分解方法，它可以将一个复杂的矩阵分解为一个或多个简单的三角形矩阵的乘积。
不同形式的三角分解有着各自的优缺点，本文将从几种不同的角度来讨论这些优缺点。
一、LU分解LU分解是将一个矩阵分解为一个下三角矩阵L和一个上三角矩阵U的乘积，即A=LU。
LU分解的优点是计算简单，因为它只需要进行一次分解即可得到L和U两个矩阵，后续的线性方程求解可以直接使用LU 分解后的矩阵进行计算。
然而，LU分解的缺点是当原始矩阵A的某些主对角线元素接近于零时，LU分解可能会失效，需要采取一些特殊的技巧来解决这个问题。
二、Cholesky分解Cholesky分解是将一个对称正定矩阵分解为一个下三角矩阵L和其转置矩阵的乘积，即A=LL^T。
Cholesky分解的优点是计算量较小，而且分解出的L矩阵的元素都是实数，因此在存储和计算上都有一定的优势。
然而，Cholesky分解的缺点是它只适用于对称正定矩阵，对于非对称矩阵或不正定矩阵是无法进行Cholesky分解的。
三、QR分解QR分解是将一个矩阵分解为一个正交矩阵Q和一个上三角矩阵R的乘积，即A=QR。
QR分解的优点是适用范围广，对于任意矩阵都可以进行QR分解，并且分解出的Q和R矩阵都具有一些良好的性质，比如Q是正交矩阵，R是上三角矩阵。
然而，QR分解的缺点是计算量较大，尤其是对于大型矩阵来说，QR分解的计算时间会比较长。
不同形式的三角分解都有各自的优缺点，选择合适的分解方法需要根据具体的问题来决定。
在实际应用中，可以根据矩阵的特点和计算需求来选择最合适的三角分解方法，以达到最优的计算效果。
研究和探索更加高效的矩阵分解方法也是矩阵论研究的重要方向之一。
四、SVD分解SVD分解是将一个矩阵分解为三个矩阵的乘积，即A=UΣV^T，其中U和V分别是正交矩阵，Σ是对角矩阵。
SVD分解的优点是适用于所有的矩阵，无论是否为方阵，而且SVD分解是唯一的，即对于每一个矩阵都存在唯一的SVD分解。
矩阵的三角分解法
矩阵的三角分解法矩阵的三角分解法是一种用于将一个矩阵分解为上三角矩阵和下三角矩阵的方法。
这种分解方法可以帮助我们更好地理解和解决矩阵相关的问题。
下面我将按要求逐段解释这个问题。
1. 什么是三角分解法三角分解法是一种将矩阵分解为上三角矩阵和下三角矩阵的方法。
在三角分解中，我们将原始矩阵分解为两个三角矩阵，一个是上三角矩阵，另一个是下三角矩阵。
上三角矩阵的主对角线以下的元素全为零，而下三角矩阵的主对角线以上的元素全为零。
这种分解法在解线性方程组、计算矩阵的行列式和求逆等问题中非常有用。
2. 如何进行三角分解三角分解的具体过程是通过一系列的行变换将原始矩阵转化为上三角矩阵或下三角矩阵。
这些行变换包括行交换、行缩放和行替换等操作。
首先，我们选择一个主元素，通常是第一行第一列的元素。
如果主元素为零，则需要进行行交换，将一个非零元素移动到主元素的位置。
然后，我们使用行缩放操作，将主元素所在列的其他元素变为零。
具体操作是将主元素所在行的每个元素除以主元素的值，然后将结果乘以其他行的主元素所在列的元素，并将其减去相应的行。
重复以上步骤，直到得到上三角矩阵或下三角矩阵。
最后，我们可以将得到的上三角矩阵和下三角矩阵合并为一个新的上三角矩阵或下三角矩阵。
3. 三角分解的应用领域有哪些三角分解法在数值计算和线性代数中有广泛的应用。
它可以用于求解线性方程组、计算矩阵的行列式和求逆等问题。
在求解线性方程组时，我们可以将系数矩阵分解为上三角矩阵和下三角矩阵，然后使用回代法或前代法来求解方程组。
这样可以简化计算过程，提高求解的精度和效率。
在计算矩阵的行列式时，我们可以通过三角分解将矩阵转化为上三角矩阵或下三角矩阵，然后将主对角线上的元素相乘即可得到行列式的值。
这种方法比直接计算行列式的方法更简单、高效。
在求解矩阵的逆时，我们可以将矩阵分解为上三角矩阵和下三角矩阵，然后通过对分解得到的上三角矩阵和下三角矩阵进行反向的行变换，得到原始矩阵的逆矩阵。
数值计算_矩阵的三角分解算法
数值计算_矩阵的三角分解算法矩阵的三角分解是一种将矩阵分解为下三角矩阵和上三角矩阵的方法。
三角分解在数学和计算机科学中都有广泛的应用，特别是在线性代数、数值计算和优化问题中。
在本篇文章中，我们将介绍几种常见的矩阵三角分解算法。
一、LU分解LU分解是矩阵三角分解中最常见的一种方法。
它将一个矩阵分解为一个下三角矩阵L和一个上三角矩阵U，使得原始矩阵A可以表示为A=LU。
其中，L矩阵的主对角线元素全为1，而U矩阵的主对角线元素是A矩阵的主对角线元素。
实际上，LU分解可以看作是高斯消元法的矩阵形式。
在进行LU分解时，我们可以通过对原始矩阵A进行一系列的行变换来得到上三角矩阵U。
同时，我们可以记录每一次行变换的乘积以及主元元素的倒数，从而得到下三角矩阵L。
因此，LU分解可以通过高斯消元法来直接实现。
二、Cholesky分解Cholesky分解是一种仅适用于对称正定矩阵的三角分解方法。
它将一个对称正定矩阵分解为一个下三角矩阵L和其转置矩阵的乘积，即A=LL^T。
Cholesky分解非常有效率，尤其适用于解线性方程组和进行矩阵的逆运算。
由于分解结果是一个下三角矩阵，因此Cholesky分解可以减少计算量并提高计算速度。
三、QR分解QR分解是一种将矩阵分解为一个正交矩阵Q和一个上三角矩阵R的方法，即A=QR。
其中，Q矩阵是正交矩阵，其列向量是正交的，而R矩阵是上三角矩阵。
QR分解可以看作是对矩阵A进行一系列的正交变换，使其变为上三角形式。
其中，每一次正交变换可以通过Givens旋转来实现，即通过矩阵的乘积来实现矩阵的旋转。
QR分解在多元线性回归分析、奇异值分解和特征值分解等领域有广泛的应用。
四、LUP分解LUP分解是LU分解的一个变种，并增加了行交换的步骤。
LUP分解将一个矩阵分解为一个下三角矩阵L、一个上三角矩阵U和一个置换矩阵P，使得PA=LU。
其中，L和U的构造方式与LU分解相同，而置换矩阵P是一个与单位矩阵相似的矩阵，用于记录行交换的信息。
5.3 矩阵的三角分解法
8
解： (1)分解A LU，令 2 5 6 1 4 13 19 l 21 6 3 6 l31 0 1 l32 0 u11 0 1 u12 u22 u13 u23 u33
24
由A L( DLT ) 1 l 21 l31 ... l n1 1 l32 ... ln 2 1 ... ... lnn 1 1 d1 ... d1l21 d2 ... d1l31 d 2 l32 d3 ... ... ... ... d1 l n1 d 2 ln2 d 3 ln 3 dn
25
由 i j时aij = l ik d k l jk l ij d j , 知
k =1
j -1
L, D元素计算公式
lij =
aij lik d k l jk
k =1
j -1
dj
j -1
( j 1, 2, ,i 1)
2 d i =aii l ik d k ( i 1, 2, , n) k =1
y1 b1 i -1 y y b l i ij i j j 1
i 2, 3, , n
( i n 1, , 1)
7
或 用 Doolittle 分解法
例:用矩阵的直接三角分解法解方程组
5 6 x1 10 2 4 13 19 x 19 2 6 3 6 30 x3
27
d1 a11
改进平方根法解方程组
1. 分解计算A=LDLT ,
d1 a11 对于i 2, 3, ..., n j 1 c a cik l jk ij ij k 1 cij ( j 1, 2, ..., i 1) lij dj i 1 d i aii cik l ik k 1
矩阵的三角分解法
矩阵的三角分解法矩阵的三角分解法，是一种重要的矩阵分解方法之一。
它可以将一个矩阵分解为上三角矩阵和下三角矩阵的乘积，这样的分解在许多数值计算和线性代数问题中都有广泛的应用。
首先，我们来介绍一下矩阵的基本概念。
矩阵是由数个数按照矩阵的规则排列形成的矩形阵列。
矩阵可以表示线性变换，用于解线性方程组，并在数据科学、图像处理等领域中有广泛应用。
矩阵由m行和n列组成，形式为m×n.。
在进行矩阵运算时，我们经常需要对矩阵进行分解，以便更好的进行计算和处理。
接下来，让我们具体了解一下矩阵的三角分解法。
矩阵的三角分解法，是一种将一个矩阵分解为上三角矩阵和下三角矩阵的方法。
具体来说，对于一个n×n的矩阵A，可以找到一个上三角矩阵U和一个下三角矩阵L，使得A=LU。
其中，U的主对角线元素全为1，而L的副对角线元素全为0。
这样的分解称为LU分解。
矩阵的三角分解法具有以下的重要性和指导意义。
首先，矩阵的三角分解法可以简化矩阵计算。
在进行矩阵乘法、求逆、解线性方程组等计算时，三角分解法能够将原始矩阵转化为上三角矩阵和下三角矩阵的乘积形式，从而简化了计算过程。
这极大地提高了计算的效率，并降低了计算的复杂性。
其次，矩阵的三角分解法对于解线性方程组非常有用。
对于一个线性方程组Ax=b，我们可以将矩阵A进行三角分解为LU形式，然后通过前代和回代的方法，可以快速地求解出x的值。
这种方法称为三角分解解法，它在求解大型线性方程组时具有较高的效率和稳定性。
此外，矩阵的三角分解法在数学理论研究中也发挥着重要作用。
矩阵的三角分解法能够将一个矩阵A分解为上三角矩阵和下三角矩阵的乘积形式，这种分解方式是独一无二的。
同时，矩阵的三角分解法还能够帮助我们理解和证明线性代数的一些基本定理，如矩阵的秩、可逆性等。
它为我们研究和深入理解矩阵及其性质提供了重要的工具和思路。
综上所述，矩阵的三角分解法是一种重要的矩阵分解方法，它可以将一个矩阵分解为上三角矩阵和下三角矩阵的乘积形式。
矩阵三角分解
矩阵三角分解1. 什么是矩阵三角分解矩阵三角分解是将矩阵分解成一个上三角矩阵和一个下三角矩阵的过程。
这种分解形式可以用于简化矩阵的乘法计算，并且对于求解线性方程组等问题也有很大的帮助。
2. 矩阵三角分解的应用矩阵三角分解在科学计算和工程应用中具有广泛的应用。
其中包括求解线性方程组、矩阵求逆、矩阵特征值计算等。
在这些应用中，矩阵三角分解能够有效地简化计算，并且提高计算速度和精度。
3. 矩阵三角分解的实现方法矩阵三角分解的实现可以采用多种方法，例如高斯消元法、Crout分解法、Cholesky分解法等。
其中，高斯消元法是最常用的方法之一，其基本思想是通过初等行变换将矩阵转化为上三角矩阵。
而Crout分解法和Cholesky分解法则是把矩阵分解成一个下三角矩阵和一个上三角矩阵。
4. 矩阵三角分解的优势相对于直接计算矩阵乘法的方法，矩阵三角分解具有一定的优势。
首先，矩阵三角分解能够避免矩阵乘法中的数值溢出问题，提高计算的稳定性和精度。
其次，矩阵三角分解的计算量通常比直接计算矩阵乘法要少，能够提高计算速度和效率。
5. 矩阵三角分解在深度学习中的应用矩阵三角分解在深度学习中也有广泛应用。
例如，在神经网络中，卷积层的计算可以通过矩阵相乘来实现。
因此，将卷积层中的卷积核分解成一个上三角矩阵和一个下三角矩阵，就能够简化卷积层的计算，并且减少计算量和存储空间。
这种方法被称为低秩近似技术，能够提高卷积层的计算速度和效率。
6. 总结矩阵三角分解是将矩阵分解成一个上三角矩阵和一个下三角矩阵的过程，其应用广泛且计算效率高。
在科学计算和工程应用中具有广泛的应用，同时在深度学习中也有着重要的作用。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
L3 ( I 3,i3 L2 I 3,i3 )( I 3,i3 I 2,i2 L1 I 2,i2 I 3,i3 )
~ ~ ~ ( I 3,i3 I 2,i2 I1,i1 ) A L3 L2 L1 PA,
（3.Байду номын сангаас）
1
其中
~ L1 I 3,i3 I 2,i2 L1 I 2,i2 I 3,i3 , ~ L2 I 3,i3 L2 I 3,i3 , ~ L3 L3 ,
~ Lk (k 1,2,3) 为单位下三角阵，其元素的绝对值不超过1.
P I 3,i3 I 2,i2 I1,i1 .
记
~ ~ ~ L L3 L2 L1 ,
1
由(3.2)得到
PA LU .
其中 P为排列矩阵， L 为单位下三角阵， U为上三角阵.
2
这说明对(2.1)应用列主元素消去法相当于对 ( A b)先 进行一系列行交换后对 PAx Pb再应用高斯消去法. 而在实际计算中只能在计算过程中做行的交换. 定理8 (列主元素的三角分解定理) 如果 A为非奇异矩阵， 则存在排列矩阵 P使
AX I n 其中 I n 为单位矩阵.
7
将 X按列分块
X x1 x2 xn , I e1 e2 en ,
于是求解 AX I n 等价于求解 n个方程组
Ax j e j , ( j 1,2, , n).
可用高斯-若当方法求解 AX I n .
8
例4 用高斯-若当方法求
6 5 3
0 0 1
0 1 0
1 0 0
1 第1次消元 0 0
1/ 3 2/3 1/ 3 c3
9
1 第2次消元 0 0 1 第3次消元 0 0
0 1 0
1/ 2 3/ 2 1/ 2
0 0 1
5/ 2 3/ 2 1/ 2 c2 3 3 1
2 1 0
0 1 0
0 0 1
1 3 2
2 1 ( I n A1 ). 0
c1
且
m1 (m11 , m21 , m31 )T c3 , m2 (m12 , m22 , m32 )T c2 ,
) 方程组 A( k ) x b ( k，其中
( A( k )
1 (k ) b )
0 1
0 0 1
a1k a2 k ak 1,k ak ,k an ,k
a1n a2 n ak 1,n ak ,n an ,n
3. 计算乘数
mik aik / akk mkk 1/ akk . (k 1,2,,n且i k ),
( mik可保存在存放 aik 的单元中).
5
4.
消元计算
aij aij mik akj k 1,2,, n且i k , j k 1,, n
m3 (m13 , m23 , m33 )T c1.
10
为了节省内存单元，可不必存放单位矩阵， c3存放在
A 的第1列， c2 存放在 A 的第二列位置， c1 存放在 A的
第3列. 经消元计算，最后再调整一下列就可在 A 的位置得到
A1 . 注意第 k 步消元时，由 A 的第 k 列
ak (a1k , , akk , , ank )T
利用(3.1)得到
Ln1I n1,in1 L2 I 2,i2 L1I1,i1 A A( n ) U .
若记 则有
~ P Ln1 I n1,in1 L2 I 2,i2 L1 I1,i1 .
~ P A U ,
~ Pb b ( n ) ,
~. 考虑 n 4 时的 P
U A( 4) L3 I 3,i3 L2 I 2,i2 L1 I1,i1 A
bi bi mik bk
k 1,2,,n且i k .
5.
计算主行
akj akj mkk ( j k , k 1,, n),
bk bk mkk .
上述过程结束后有
( A b) ( A( k 1)
1 b ( k 1) )
1 1
ˆ b 1 ˆ b 2 . ˆ b n
1 A 2 3 2 4 5 3 5 6
的逆矩阵. 解
1 C 2 3 2 4 5
3 5 6
1 0 0
0 1 0 5/ 3 2/3 1/ 3
0 3 0 2 r 1 r3 1 1 2 1 1 0 0 1 0 1 0
5 4 2
b1 b2 bk 1 . bk bn
4
在第 k步计算时 (k 1,2,,n), 对上述矩阵第 k 行的上、下 都进行消元.
1. 按列选主元素，即确定 ik 使
aik ,k max aik .
k in
2. 换行(当 ik k 时)交换 ( A b)第 k行与第 ik 行元素.
PA LU ,
其中 L 为单位下三角阵，U 为上三角阵. 编程时， L 元素存放在数组 A的下三角部分， U元素 存放在 A上三角部分，由记录主行的整型数组 Ip(n) 可知
P 的情况.
3
5.3.2
高斯-若当消去法
高斯消去法中，若同时消去对角线下方和上方的元素， 这种方法称为高斯-若当(Gauss-Jordan)消去法. 设用高斯-若当消去法已完成 k 1步， Ax b 化为等价
6
用高斯-若当方法将 A约化为单位矩阵，计算解就在 常数项位置得到，用不着回代求解，计算量大约需要 n 3 / 2 次乘除法，比高斯消去法大，但用高斯-若当方法求矩阵 的逆矩阵还是比较合适的.
定理9(高斯-若当法求逆矩阵) 设 A为非奇异矩阵，
方程组 AX I n的增广矩阵为 C ( A I n ) . 如果对 C应 用高斯-若当方法化为 ( I n T ) ， 则 A1 T . 事实上，求 A的逆矩阵 A1 ，即求 n阶矩阵 X，使