6.4.1 追赶法

三角法在对角占优阵中的应用
从三角分解法的角度来看,追赶法实际上是一种特殊的LU分解。对于三对角矩阵,我们可以将其分解为一个特殊的下三角矩阵L和上三角矩阵U的乘积,这种分解具有独特的稀疏结构。

三对角方程组的特点

三对角方程组的系数矩阵A为三对角矩阵,形式如下:
方程组可以写为:

将三对角矩阵A分解为LU的形式:
通过比较系数可得:
可以看出:
  • L矩阵是一个单位下三角矩阵,只在次对角线上有非零元素
  • U矩阵是一个上三角矩阵,只在主对角线和上对角线上有非零元素
  • 这种特殊的分解形式大大减少了计算量和存储空间
追赶法中的α实际上就是U的主对角线元素,这种LU分解的形式帮助我们理解追赶法的本质。

追赶法的计算步骤

追赶法分为两个步骤:

1. 追的过程(前代换)

计算系数:
计算常数项:
在实际计算中,将A写成增广矩阵 的形式,最后一列就是上面算的“常数项”

2. 赶的过程(回代)

实际上,"赶"的过程就是求解上三角方程组Ux=β。注意, 就是原来的A中的 前面已经通过迭代(追)解出。
当我们拿到β向量后,由于U是上三角矩阵,我们可以从最后一个方程开始,依次向上求解每个未知数。
这种回代过程之所以称为"赶",是因为我们是从下往上逐个消除未知数,就像把未知数一个个"赶"出方程一样。
从后向前依次求解:

算法特点

  • 计算量小:运算次数为O(n)
  • 存储空间少:只需要存储三条对角线上的元素
  • 适用于大型三对角方程组的求解

算法流程图

收敛条件

追赶法的收敛条件要求系数矩阵满足:
  • 系数矩阵为严格对角占优矩阵
  • 所有主对角元素不为零(也就是严格占优,不能取等号)
 

6.4.2 改进平方根法

三角法在对称正定阵中的应用

1. 基本概念

改进平方根法是一种用于求解对称正定线性方程组的高效算法。它是平方根法的改进版本,通过利用矩阵的对称性质来减少计算量。

2. 适用条件

  • 系数矩阵A必须是对称矩阵
  • 所有顺序主子式均大于零(正定性质)
  • 方程组形式为 Ax = b

3. 矩阵分解

对于满足条件的矩阵A,可以将其分解为:
其中:
  • L是单位下三角矩阵(主对角线上元素均为1)
  • U是上三角矩阵
分解的具体形式为:
 
notion image
notion image
例3.6 用改进的Cholesky分解法解方程组:
解题步骤:
第1步:矩阵构造
第2步:计算
得到等价的三角形方程组(LUx=y)为:
回代求解,得:
Loading...
Z_cosy
Z_cosy
浙江大学电气工程学院本科生
最新发布
Dairy
2025-4-21
MATH
2025-3-23
English
2025-3-21
Courses
2025-3-21
公告
🎉Welcome to Z-cosy🎉
-- 食用指南 ---
目前只有课程笔记以及电控学习笔记
陆续会整理更多内容!