矩阵代数（五）- 矩阵因式分解

小结

$\boldsymbol{LU}$ 分解
$\boldsymbol{LU}$ 分解算法

矩阵 $\boldsymbol{A}$ 的因式分解是把 $\boldsymbol{A}$ 表示为两个或更多个矩阵的乘积。

$\boldsymbol{LU}$ 分解

设 $\boldsymbol{A}$ 是 $m \times n$ 矩阵，它可以行换简为阶梯形而不必行对换（此后，我们将处理一般情形），则 $\boldsymbol{A}$ 可写成 $\boldsymbol{A}=\boldsymbol{LU}$ ， $\boldsymbol{L}$ 是 $m \times m$ 下三角矩阵，主对角线全是1， $\boldsymbol{U}$ 是 $\boldsymbol{A}$ 的一个 $m \times n$ 阶梯形矩阵。这样一个分解称为 $\boldsymbol{LU}$ 分解，矩阵 $\boldsymbol{L}$ 是可逆的，称为单位下三角矩阵。

当 $\boldsymbol{A}=\boldsymbol{LU}$ 时，方程 $\boldsymbol{Ax}=\boldsymbol{x}$ 可写成 $\boldsymbol{L}(\boldsymbol{Ux})=\boldsymbol{b}$ 。把 $\boldsymbol{Ux}$ 写成 $\boldsymbol{y}$ ，可以由解下面一对方程来求解 $\boldsymbol{x}$ ： $\begin{cases} \boldsymbol{Ly}=\boldsymbol{b} \\ \boldsymbol{Ux} = \boldsymbol{y} \end{cases}$

可以证明 $\boldsymbol{A}=\begin{bmatrix}3 & -7 & -2 & 2 \\ -3 & 5 & 1 & 0 \\ 6 & -4 & 0 & -5 \\ -9 & 5 & -5 & 12 \end{bmatrix}=\begin{bmatrix}1 & 0 & 0 & 0 \\ -1 & 1 & 0 & 0 \\ 2 & -5 & 1 & 0 \\ -3 & 8 & 3 & 1 \end{bmatrix}\begin{bmatrix}3 & -7 & -2 & 2 \\ 0 & -2 & -1 & 2 \\ 0 & 0 & -1 & 1 \\ 0 & 0 & 0 & -1 \end{bmatrix}=\boldsymbol{LU}$
应用 $\boldsymbol{A}$ 的 $\boldsymbol{LU}$ 分解来解 $\boldsymbol{Ax}=\boldsymbol{b}$ ，其中 $\boldsymbol{b}=\begin{bmatrix}-9 \\ 5 \\ 7 \\ 11\end{bmatrix}$ 。
解：解 $\boldsymbol{Ly}=\boldsymbol{b}$ 。
$\begin{bmatrix} \boldsymbol{L} & \boldsymbol{b} \end{bmatrix}=\begin{bmatrix}1 & 0 & 0 & 0 & -9 \\ -1 & 1 & 0 & 0 & 5 \\ 2 & -5 & 1 & 0 & 7 \\ -3 & 8 & 3 & 1 & 11 \end{bmatrix}$ ～ $\begin{bmatrix} 1 & 0 & 0 & 0 & -9 \\ 0 & 1 & 0 & 0 & -4 \\ 0 & 0 & 1 & 0 & 5 \\ 0 & 0 & 0 & 1 & 1 \end{bmatrix}=\begin{bmatrix} \boldsymbol{I} & \boldsymbol{y}\end{bmatrix}$
对 $\boldsymbol{Ux}=\boldsymbol{y}$ 进行行化简的向后步骤。
$\begin{bmatrix} \boldsymbol{U} & \boldsymbol{y} \end{bmatrix}=\begin{bmatrix}3 & -7 & -2 & 2 & -9 \\ 0 & -2 & -1 & 2 & 4 \\ 0 & 0 & -1 & 1 & 5 \\ 0 & 0 & 0 & -1 & 1 \end{bmatrix}$ ～ $\begin{bmatrix} 1 & 0 & 0 & 0 & 3 \\ 0 & 1 & 0 & 0 & 4 \\ 0 & 0 & 1 & 0 & -6 \\ 0 & 0 & 0 & 1 & -1 \end{bmatrix}$
故 $\boldsymbol{x}=\begin{bmatrix}3 \\ 4 \\ -6 \\ -1\end{bmatrix}$ 。

$\boldsymbol{LU}$ 分解的计算依赖于如何求 $\boldsymbol{L}$ 和 $\boldsymbol{U}$ 。

$\boldsymbol{LU}$ 分解算法

设 $\boldsymbol{A}$ 可以化为阶梯形 $\boldsymbol{U}$ ，化简过程中仅用行倍加变换，即把一行的倍数加于它下面的另一行。这样，存在单位下三角初等矩阵 $\boldsymbol{E_1},\cdots,\boldsymbol{E_p}$ 使 $\boldsymbol{E_p}\cdots\boldsymbol{E_1}\boldsymbol{A}=\boldsymbol{U}$ 。于是 $\boldsymbol{A}=(\boldsymbol{E_p}\cdots\boldsymbol{E_1})^{-1}\boldsymbol{U}=\boldsymbol{LU}$ ，其中 $\boldsymbol{L}=(\boldsymbol{E_p}\cdots\boldsymbol{E_1})^{-1}$ 。可以证明 $\boldsymbol{L}$ 是单位下三角矩阵。
注意将 $\boldsymbol{A}$ 化为阶梯形 $\boldsymbol{U}$ 过程中的行变换，它把 $\boldsymbol{A}$ 化为 $\boldsymbol{U}$ 。这写行变换也把 $\boldsymbol{L}$ 化为 $\boldsymbol{I}$ ，这是因为 $\boldsymbol{E_p}\cdots\boldsymbol{E_1}\boldsymbol{L}=\boldsymbol{E_p}\cdots\boldsymbol{E_1}(\boldsymbol{E_p}\cdots\boldsymbol{E_1})^{-1}=\boldsymbol{I}$

$\boldsymbol{LU}$ 分解的算法：

如果可能的话，用一系列的行倍加变换把 $\boldsymbol{A}$ 化为阶梯形 $\boldsymbol{U}$ 。
填充 $\boldsymbol{L}$ 的元素使相同的行变换把 $\boldsymbol{L}$ 变为 $\boldsymbol{I}$ 。

求下列矩阵的 $\boldsymbol{LU}$ 分解： $\boldsymbol{A}=\begin{bmatrix} 2 & 4 & -1 & 5 & -2 \\ -4 & -5 & 3 & -8 & 1 \\ 2 & -5 & -4 & 1 & 8 \\ -6 & 0 & 7 & -3 & 1 \end{bmatrix}$
解：因 $\boldsymbol{A}$ 有4行，故 $\boldsymbol{L}$ 应为 $4 \times 4$ 矩阵。 $\boldsymbol{L}$ 的第一列应该是 $\boldsymbol{A}$ 的第一列除以它的第一行主元素： $\boldsymbol{L}=\begin{bmatrix} 1 & 0 & 0 & 0 \\ -2 & 1 & 0 & 0 \\ 1 && 1 & 0 \\ -3 & & & 1\end{bmatrix}$
比较 $\boldsymbol{A}$ 和 $\boldsymbol{L}$ 的第一列。把 $\boldsymbol{A}$ 的第一列的后三个元素变成零的行变换同时也将 $\boldsymbol{L}$ 的第一列的后三个元素变成0。
$\boldsymbol{A}=\begin{bmatrix} \color{red}{2} & 4 & -1 & 5 & -2 \\ \color{red}{-4} & -5 & 3 & -8 & 1 \\ \color{red}{2} & -5 & -4 & 1 & 8 \\ \color{red}{-6} & 0 & 7 & -3 & 1 \end{bmatrix}$ ～ $\begin{bmatrix} 2 & 4 & -1 & 5 & -2 \\ 0 & \color{red}{3} & 1 & 2 & -3 \\ 0 & \color{red}{-9} & -3 & -4 & 10 \\ 0 & \color{red}{12} & 4 & 12 & -5 \end{bmatrix}$ ～ $\begin{bmatrix} 2 & 4 & -1 & 5 & -2 \\ 0 & 3 & 1 & 2 & -3 \\ 0 & 0 & 0 & \color{red}{2} & 1 \\ 0 & 0 & 0 & \color{red}{4} & 7 \end{bmatrix}$ ～ $\begin{bmatrix} 2 & 4 & -1 & 5 & -2 \\ 0 & 3 & 1 & 2 & -3 \\ 0 & 0 & 0 & 2 & 1 \\ 0 & 0 & 0 & 0 & \color{red}{5} \end{bmatrix}$
上式中标出的元素确定来将 $\boldsymbol{A}$ 化为 $\boldsymbol{U}$ 的行化简。在每个主元列，把标出的元素除以主元后将结果放入 $\boldsymbol{L}$ ： $\boldsymbol{L}=\begin{bmatrix} 1 & 0 & 0 & 0 \\ -2 & 1 & 0 & 0 \\ 1 & -3 & 1 & 0 \\ -3 & 4 & 2 & 1\end{bmatrix}$ 。
容易证明，所求出的 $\boldsymbol{L}$ 和 $\boldsymbol{U}$ 满足 $\boldsymbol{LU}=\boldsymbol{A}$ 。

人面猴
序言：七十年代末，一起剥皮案震惊了整个滨河市，随后出现的几起案子，更是在滨河造成了极大的恐慌，老刑警刘岩，带你破解...
沈念sama阅读 225,226评论 6赞 524
死咒
序言：滨河连续发生了三起死亡事件，死亡现场离奇诡异，居然都是意外死亡，警方通过查阅死者的电脑和手机，发现死者居然都...
沈念sama阅读 96,509评论 3赞 405
救了他两次的神仙让他今天三更去死
文/潘晓璐我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
开封第一讲书人阅读 172,523评论 0赞 370
道士缉凶录：失踪的卖姜人
文/不坏的土叔我叫张陵，是天一观的道长。经常有香客问我，道长，这世上最难降的妖魔是什么？我笑而不...
开封第一讲书人阅读 61,181评论 1赞 302
港岛之恋（遗憾婚礼）
正文为了忘掉前任，我火速办了婚礼，结果婚礼上，老公的妹妹穿的比我还像新娘。我一直安慰自己，他们只是感情好，可当我...
茶点故事阅读 70,189评论 6赞 401
恶毒庶女顶嫁案：这布局不是一般人想出来的
文/花漫我一把揭开白布。她就那样静静地躺着，像睡着了一般。火红的嫁衣衬着肌肤如雪。梳的纹丝不乱的头发上，一...
开封第一讲书人阅读 53,642评论 1赞 316
城市分裂传说
那天，我揣着相机与录音，去河边找鬼。笑死，一个胖子当着我的面吹牛，可吹牛的内容都是我干的。我是一名探鬼主播，决...
沈念sama阅读 41,993评论 3赞 431
双鸳鸯连环套：你想象不到人心有多黑
文/苍兰香墨我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
开封第一讲书人阅读 40,977评论 0赞 280
万荣杀人案实录
序言：老挝万荣一对情侣失踪，失踪者是张志新（化名）和其女友刘颖，没想到半个月后，有当地人在树林里发现了一具尸体，经...
沈念sama阅读 47,527评论 1赞 326
护林员之死
正文独居荒郊野岭守林人离奇死亡，尸身上长有42处带血的脓包…… 初始之章·张勋以下内容为张勋视角年9月15日...
茶点故事阅读 39,547评论 3赞 347
白月光启示录
正文我和宋清朗相恋三年，在试婚纱的时候发现自己被绿了。大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
茶点故事阅读 41,661评论 1赞 355
活死人
序言：一个原本活蹦乱跳的男人离奇死亡，死状恐怖，灵堂内的尸体忽然破棺而出，到底是诈尸还是另有隐情，我是刑警宁泽，带...
沈念sama阅读 37,250评论 5赞 351
日本核电站爆炸内幕
正文年R本政府宣布，位于F岛的核电站，受9级特大地震影响，放射性物质发生泄漏。R本人自食恶果不足惜，却给世界环境...
茶点故事阅读 42,991评论 3赞 340
男人毒药：我在死后第九天来索命
文/蒙蒙一、第九天我趴在偏房一处隐蔽的房顶上张望。院中可真热闹，春花似锦、人声如沸。这庄子的主人今日做“春日...
开封第一讲书人阅读 33,422评论 0赞 25
一桩弑父案，背后竟有这般阴谋
文/苍兰香墨我抬头看了看天上的太阳。三九已至，却和暖如春，着一层夹袄步出监牢的瞬间，已是汗流浃背。一阵脚步声响...
开封第一讲书人阅读 34,571评论 1赞 277
情欲美人皮
我被黑心中介骗来泰国打工，没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留，地道东北人。一个月前我还...
沈念sama阅读 50,241评论 3赞 382
代替公主和亲
正文我出身青楼，却偏偏与公主长得像，于是被迫代替她去往敌国和亲。传闻我的和亲对象是个残疾皇子，可洞房花烛夜当晚...
茶点故事阅读 46,737评论 2赞 366

矩阵代数（五）- 矩阵因式分解

小结

分解

分解算法

推荐阅读更多精彩内容

$\boldsymbol{LU}$ 分解

$\boldsymbol{LU}$ 分解算法