矩阵的LU分解(自编MATLAB)实验报告.docVIP

. . 1矩阵的LU分解 1.1 LU 分解原理 定理：设AC?n?n，如果 A 的顺序主子式 A11≠0, a11a 则存在唯一的主对角线上元素全为 1 的下三角矩阵L与唯一的上三角矩阵 U，使得 A=LU. 证明：对矩阵A的阶数使用数学归纳法. 显然，当 n=1 时，A11=1 ?A11 就是唯一的分解式。现假定对 n-1 阶矩阵，定理的结论成立。对 A A=A 其中α1,α2∈Cn-1.由于 n-1 阶矩阵 An-1的 k 阶顺序主子式就是 A 的 k 阶主子式（k=1,2，…，n-2），故它们都不为零. A 其中Ln-1的主对角线上的元素都1. An-1=a11a 所以Ln-1及U 先假设已有 A=LU，其中 L=Ln-10γT β，γ LU = Ln-1U 则β，γ Ln-1β= 注意到Ln-1及 β=L 这就证明了 A 的 LU 分解的存在性和唯一性. 1.2 LU分解算法 当 n 阶矩阵满足定理的条件时，可以用初等变换的方法求出 L 和 U. 因为当 A=LU 时，由于 L 可逆，故必存在可逆矩阵 P 使得 PL 即 PA=PLU=U.也就是说，可以先对 A 施行行的初等变换得出上三角矩阵U，而矩阵 P可以通过对单位矩阵I进行相同的行初等变换得出，即 P(A,I) =(PA,PI) =(U,P) 于是A=P-1U,为保持P 1.3 LU分解用