高等工程数学 —— 第三章（1）三角分解和满秩分解

最新推荐文章于 2024-03-30 20:05:04 发布

梦里一声何处鸿

最新推荐文章于 2024-03-30 20:05:04 发布

阅读量1.5k

点赞数 2

分类专栏：高等工程数学文章标签：算法机器学习人工智能

本文链接：https://blog.csdn.net/Deam_swan_goose/article/details/127258173

版权

13 篇文章 44 订阅

订阅专栏

在这里插入图片描述感觉不用特别抠概念，看例题会方法就行。

简单讲就是将其进行如下分解：
在这里插入图片描述
计算公式如下：

哈哈哈哈，看不懂是吧。那就别看了，看例题解法就行了。

简单说就是这样：

在这里插入图片描述
例：

这里A的第一行不变就作为U的第一行。L的对角元全写1，第一列的其他元素 $l_{k1} = \frac{a_{k1}} {u_{11} }\quad$ .
心里去想那个分层的框，然后先算行。算u就是用原来的a的值去减去它每一层左边和上边的框外元素的乘积。算 $l$ 就是在此基础上在除一个 $u_{kk}$ 。例如：对于 $l_{32}$ 而言就是原来的 $a_{32}$ 减去框外左边的 $l_{31}$ 与框外上边的 $u_{12}$ 的乘积，再处于对应列的 $u_{kk}$ ，这里为 $u_{22}$ 。

对于上述Doolittle分解存在以下问题：
在这里插入图片描述

在真实算法计算时可能会用一些误差，如果 $u_{kk}$ 的值特别小那就会将误差放大影响结果。当 $u_{kk} = 0$ 时那就直接不能运行了。所以这里我们选列主元其实就是选每一列的最大元素行变换到对角元作为 $u_{kk}$ 。注意看列元素最大值时 $l$ 的值不用除以 $u_{kk}$ ，因为这里 $u_{kk}$ 并没有确定。
我们可以用一个置换矩阵 $P$ 来记录我们行变换的结果，就是先让 $P$ 作为单位阵，然后跟着前面做同样的行变换就好了。