2.5 行阶梯形变换等价于矩阵乘法—LU分解.pdf

2.5 行阶梯形变换等价于矩阵乘法——LU 分解 在行阶梯形变换中，如果不变换成最简形 U ，只变换成行阶梯形矩阵U ，也就是把式 0 (2.4.1) 中的所有数乘变换矩阵Dk 都去掉，只保留消元变换E 和行交换矩阵 P ，则式(2.4.1) 可写成： U=(行阶梯变换中所有P 及E 的连乘积)*A 由于初等变换矩阵P 和E 的逆矩阵存在，其乘积也是可逆的，故可令行阶梯形变换中 所有P 和E 的连乘积矩阵的逆矩阵为L ，即 (行阶梯变换中所有P 及E 的连乘积)-1=L (2.5.1) 在式(2.5.1)两端同乘以L ，可写成 LU=A (2.5.2) 这种把矩阵A 通过初等矩阵左乘分解为一个对角元素全为 1 的下三角矩阵和一个上三 角矩阵乘积的变换称为 LU 变换。MATLAB 提供了矩阵的三角分解函数 lu.m，其调用格 式为： [L,U]=lu(A) 它返回的结果是一个有关行交换的下三角矩阵L 和一个上三角矩阵 U。这个变换程序实质 上是高斯消元法的另一种形式。 例2.12 用例1.4 中矩阵A ，b 的数据，用矩阵乘法求其行阶梯形变换的解。 3 2 2 4    解 该题的增广矩阵为：C A b ，按该例题解题过程可知三次消 [ , ]  3 3 1 5   2 2 1 4    元所需的消元元素e 应为e(2,1)=-3/3，e(3,1) =-2/3 及e(3,2)= -2/3，相应的初等消元矩阵为：  1 0 0  1 0 0 1 0 0       E1  1 1 0 , E2  0 1 0 , E3  0 1 0        0 0 1 2 / 3 0 1  0 2 / 3 1        用矩阵乘法可求出  1 0 03 2 2 4   3 2 2 4     