matlab lu分解求线性方程组_计算方法（二）直接三角分解法解线性方程组

封面是WH2里春希在编辑部的上司麻理前辈，有一说一，这条线的第一次H有点恶趣味，不是很喜欢。

---

一：概述

矩阵分解我学过的挺多种，比如极分解，谱分解，满秩分解，正交三角分解还有这里的直接三角分解大部分我都没有具体运用的经验。但是这里的三角分解的应用就很直白了，就是把矩阵分解为规律的三角矩阵后，我们就能用上次上篇文章里那种解出一个值然后不断回代的方式得到方程组的解。

二：具体步骤

计算方法课上老师只讲了一个很机械的方法，先用一个例子说明吧（带下标的通式打起来费劲看着也费劲）。

图一

我们可以通过如下步骤使之变成两个三角矩阵的乘积。

图二

矩阵L

矩阵U

已知

使用两次回代公式，

第一次求出y，第二次求出x。

三：原理分析

图二中的计算步骤很容易理解，线性代数课程中我们学习过求标准正交基的施密特方法，于是很容易联想到上述计算步骤的原理：

（1）首先提取出矩阵列向量：

（2）然后从第一行开始把主元素上方的元素消去化为0，并把主元素化为1：

（3）再从上式中反解出

的关系，很容易得出矩阵L即为

为列空间组成的正交矩阵，矩阵U为从

的过渡矩阵。

（实际上正交三角分解即为把上述步骤替换为施密特方法，而上文这种正交化方法叫啥，评论区有大神可以告诉我一下）

（查到了，叫Dolittle分解，如果矩阵L对角元不进行单位化，称为Courant分解。）

四：算法实现（MATLAB）

function

五：总结

刚学的时候对这种计算方法只是机械的记忆，今天要写这篇文章的时候才对计算的原理思考了一下，其实类比施密特正交化，我们很容易就能想到。

（插个旗子：计算方法课程笔记结束后，开一篇文章，专门研究一下各种矩阵分解的应用，应该很有趣）

上一篇

水法天天：计算方法（一）：GAUSS消去法与主元素法解线性方程组​zhuanlan.zhihu.com

下一篇

水法天天：计算方法（三）平方根法及其改进解线性方程组​zhuanlan.zhihu.com