文章目录
引言一、预备知识1.1 矩阵基本运算和和行列式1.2 矩阵的特征值和谱半径1.3 特殊矩阵二、Gauss消去法2.1 不消主元的Gauss消去法2.2 高斯主元素消去法1. Gauss列主元消去法2. Gauss全主元消去法三、矩阵的三角分解法3.1 用直接三角分解法求解方程组3.2 Doolittle分解法(杜利特尔)题目引言
研究数值解法的必要性线性代数方程组的常用解法 直接法
只包含有限次四则运算。在计算过程中不发生舍入误差的假定下,计算结果就是原方程组的精确解。(适用于小规模的n阶稠密线性方程组)迭代法
从解的某个近似值出发,通过构造一个无穷序列去逼近精确解的方法(一般有限步内得不到精确解)。实现这一极限过程每一步的结果是把前一步所得的结果施行相同的演算步骤得到的。(适用于大规模的n阶稀疏线性方程组)
Remark:由于运算过程中舍入误差的存在,实际上
直接方法得到的解也是方程组的近似解。
一、预备知识
1.1 矩阵基本运算和和行列式
1.2 矩阵的特征值和谱半径
设A是n阶方阵,如果数λ和n维非零列向量x使关系式Ax=λx成立,那么这样的数λ称为矩阵A特征值,非零向量x称为A的对应于特征值λ的特征向量。
特征值性质
tr(A)tr(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列主元消去法
akk(k)a_{kk}^{(k)}akk(k) 称为第k步的主元
Gauss列主元消去法:
第一步:在第一列中选取最大值的元素,并将该行ai,1a_{i,1}ai,1与第一行a1,1a_{1,1}a1,1进行交换。(若该列中第一行a1,1a_{1,1}a1,1最大,无需交换)
第k步:设第 步将增广矩阵(A,b)已化为
结论:当kkk从111变化到n−1n-1n−1时,(A,b)(A,b)(A,b)成为上三角形式,再用回代公式求xxx.
例. 用Gauss列主元解方程组
2. Gauss全主元消去法
概念不好理解,直接上例题+解析。
例.
通常用列主元消去法。
三、矩阵的三角分解法
3.1 用直接三角分解法求解方程组
高斯消元与LU分解的等价性L与U⭐⭐
矩阵三角分解的基本定理
3.2 Doolittle分解法(杜利特尔)
一般计算公式⭐
习题
题目
例题谱半径,就是特征值绝对值(复数取模)中的最大值,先求特征值。再取模,分别得到√5,√5,因此谱半径是√5。
例题