根据matlab--n阶非奇异方阵lu分解实现.pdfVIP

基于 Matlab 的 n 阶方阵的 LDU 分解实现 1. 引言 矩阵的 LDU 分解是 “矩阵理论与方法”课程中非常重要的一部分 LDU 分 解在实际工程应用中也非常广泛 LDU 分解可以将一个矩阵分解为一个下三角 矩阵和一个对角矩阵和一个上三角矩阵的乘积 LDU 分解主要应用在数值分析 中，用来解线性方程、求反矩阵或计算行列式 将系数矩阵 A 转变成等价两个矩阵 L 和 U 和对角矩阵的乘积 ，其中 L 和 U 分别是下三角和上三角矩阵，D 为对角矩阵。当 A 的所有顺序主子式都不为 0 时，矩阵 A 可以分解为 A=LDU 。即： Matlab 是很好的处理矩阵的工具。它的功能非常强大，包括创建矩阵，对 矩阵求逆，转置等操作非常简单，使其成为图像处理，信号分析等领域常用的工 具 Matlab 官方已经包括了对非奇异矩阵的 LU 分解函数[L,U]=lu(A)，为了加深 对矩阵分解的理解，本文不采用 Matlab 官方的 LU 分解函数对矩阵 A 进行 LDU 分解，而是根据理论推导和编程实现 LDU 分解 2. 程序设计 2.1. 输入合法检验 LU 分解需要被分解矩阵 A 满足如下条件： 1) 矩阵 A 为方阵 2) A 的顺序主子式 故 LU 分解需先检验 A 为 n 阶方阵，然后检验 A 的 n-1 个顺序主子式k 全 不为 0，才可进行 LU 分解。而检验主子式k 可以在 n-1 次循环 LU 分解中 进行，故先检验矩阵是否为方阵 代码如下 2.2. n-1 次循环 LDU 分解 LDU 分解本质上是高斯消元法。实质上是将 A 通过初等行变换变成一个 上三角矩阵，其变换矩阵就是一个单位下三角矩阵。从下至上地对矩阵 A 做初等行变换，将对角线左下方的元素变成零，然后再证明这些行变换的 效果等同于左乘一系列单位下三角矩阵，这一系列单位下三角矩阵的乘积 的逆就是 L 矩阵，它也是一个单位下三角矩阵 LDU 分解主要分为两步： i i i1 1 根据高斯消元法对A   消元，消元矩阵为L 1 ；2 计算L 1 A   A   以产 i1 i1 i 生下一步迭代的A   2.2.1. 根据 构造L ,L 1(j  i  1) j j i i 高斯消元A   ，使A   第 i+1 列从第 i+2 行至 n 行都为 0 。构造消元矩阵 L ,L 1 。首先判断 是否为 0，为 0 则无法继续分解，退出；否则继续 j j k 代码如下 i i1 L 1 A   A   2.2.2. 计算 i1 i i1