《数值分析》-- 高斯消去法与矩阵三角分解法(LU分解)

Posted 2022-12-08 胜天半月子

文章目录

• 引言
• 一、预备知识
• • 1.1 矩阵基本运算和和行列式
  • 1.2 矩阵的特征值和谱半径
  • 1.3 特殊矩阵
• 二、Gauss消去法
• • 2.1 不消主元的Gauss消去法
  • 2.2 高斯主元素消去法
  • • 1. Gauss列主元消去法
    • 2. Gauss全主元消去法
• 三、矩阵的三角分解法
• • 3.1 用直接三角分解法求解方程组
  • 3.2 Doolittle分解法(杜利特尔)
• 题目

---

引言

1. 研究数值解法的必要性
   
   
2. 线性代数方程组的常用解法

• 直接法
  只包含有限次四则运算。在计算过程中不发生舍入误差的假定下，计算结果就是原方程组的精确解。（适用于小规模的n阶稠密线性方程组）
• 迭代法
  从解的某个近似值出发，通过构造一个无穷序列去逼近精确解的方法(一般有限步内得不到精确解)。实现这一极限过程每一步的结果是把前一步所得的结果施行相同的演算步骤得到的。（适用于大规模的n阶稀疏线性方程组）
  Remark：由于运算过程中舍入误差的存在，实际上
  直接方法得到的解也是方程组的近似解。

---

一、预备知识

1.1 矩阵基本运算和和行列式

1.2 矩阵的特征值和谱半径

设A是n阶方阵，如果数λ和n维非零列向量x使关系式Ax=λx成立，那么这样的数λ称为矩阵A特征值，非零向量x称为A的对应于特征值λ的特征向量。

• 特征值性质
  
  $tr(A)$为A 的迹。
  

• 注意
  

• 谱半径
  
  习题

1.3 特殊矩阵

• 对称正定矩阵
  

---

二、Gauss消去法

• 方程组的矩阵表示
  
  高斯消元的基本思想：将矩阵A的下三角部分消为零，即化为上三角形

2.1 不消主元的Gauss消去法

首先将方程组Ax=b化为上三角方程组，此过程称为消去过程，再求解上三角方程组，此过程称为回代过程.

例.


总结：Gauss消去法解Ax=b分两步:消元、回代.

• 不消主元的Gauss消去法
  将方程组Ax=b的系数矩阵与右端项合并为：
  
  第一步：
  
  第二步：
  
  第k步：
  
  第n-1步：
  
  
  回代：
  
• 定理
  
  

2.2 高斯主元素消去法

在计算过程中的舍入误差增大能得到控制，该方法就是稳定的。
小主元是不稳定的根源，这就需要采用“选主元素”技术，即选取绝对值最大的元素作为主元。

1. Gauss列主元消去法

$a_{kk}^{(k)}$ 称为第k步的主元

• Gauss列主元消去法:
  第一步：在第一列中选取最大值的元素，并将该行$a_{i,1}$与第一行$a_{1,1}$进行交换。(若该列中第一行$a_{1,1}$最大，无需交换)
  
  
  第k步：设第 步将增广矩阵(A,b)已化为
  
  
  结论：当$k$从$1$变化到$n-1$时,$(A,b)$成为上三角形式，再用回代公式求$x$.
  例. 用Gauss列主元解方程组
  

2. Gauss全主元消去法

概念不好理解，直接上例题+解析。
例.

通常用列主元消去法。

---

三、矩阵的三角分解法

3.1 用直接三角分解法求解方程组

1. 高斯消元与LU分解的等价性
   
   
   
   
   
   

• L与U ⭐⭐
  

2. 矩阵三角分解的基本定理
   
   

3.2 Doolittle分解法(杜利特尔)

• 一般计算公式⭐
  
  
  习题

---

题目

• 例题

1. 
   谱半径，就是特征值绝对值（复数取模）中的最大值，先求特征值。再取模，分别得到√5，√5，因此谱半径是√5。
2. 
   

• 例题

1. 
2. 
   
3.