数值分析第五章线性方程组-数值分析课件.ppt

3.2.4 高斯主元素消去法 基本思想：每次消元之前在系数矩阵中按一定的范围选取绝对值最大的元素作为主元素，以便减少舍入误差的影响。 交换原则：通过方程或变量次序的交换，使在对角线位置上获得绝对值尽可能大的系数作为akk(k)，称这样的akk(k) 为主元素，并称使用主元素的消元法为主元素法 根据主元素选取范围分为：列主元素法、行主元素法、全主元素法 列主元素法：在待消元的所在列中选择主元，经方程的变换，置主元素于对角线位置后进行消元的方法。 全主元素：在全体待选系数中选取主元，则得全主元素法。 主元素法的意义 例3.3 用全主元素法解下列线组 计算m21=-19/40=0.475，m31=4/40=0.1 (5)- m21(4), (6)- m31(4) 消去x2 得 保留有主元素的方程 3.2.4.1 列主元素法 列主元素法就是在待消元的所在列中选取主元，经方程的行交换，置主元素于对角线位置后进行消元的方法。 即： 例3.4 用列主元素法解下列线性方程组 (5)- m21(4), (6)- m31(4)得 保留有主元素的方程 为了研究线性方程组近似解的误差估计 和迭代法的收敛性, 有必要对向量及矩阵的 “大小”引进某种度量----范数的概念。向量范 数是用来度量向量长度的,它可以看成是二、 三维解析几何中向量长度概念的推广。用Rn 表示n维实向量空间。 定义3.2 对任一向量X?Rn, 按照一定规则确定一个实 数与它对应, 该实数记为||X||, 若||X||满足下面三个 性质: (1) ||X||?0；||X||=0当且仅当X=0； (2) 对任意实数?， || ? X||=| ? | ||X||； 对任意向量Y?Rn，||X+Y|| ? ||X||+||Y|| 则称该实数||X||为向量X的范数 在Rn中，常用的几种范数有： 当不需要指明使用哪一种向量范数时，就用记号||.||泛指任何一种向量范数。 有了向量的范数就可以用它来衡量向量的大小和表示向量的误差。 设x*为Ax=b的精确解，x为其近似解，则其绝对误差可表示成||x- x* ||，其相对误差可表示成 定义3.5（矩阵的范数）如果矩阵 的某个非负的实值函数 ，满足 本文观看结束！！！ 祝各位身体健康！万事如意！！ 3.4 平方根法 工程实际计算中,线性方程组的系数矩阵常常具有对称正定性，其各阶顺序主子式及全部特征值均大于0。矩阵的这一特性使它的三角分解也有更简单的形式，从而导出一些特殊的解法，如平方根法与改进的平方根法。 定理3.6 设A是对称正定矩阵，则存在惟一的对角元素均为正数的下三角阵L，使A=LLT 证：因A是正定矩阵, A的顺序主子式 ?i0, i=1,2,…,n 因此存在惟一的分解 A=LU L是单位下三角阵, U是上三角阵, 将U再分解 其中D为对角阵, U0为单位上三角阵，于是 A = L U = L D U0 又 A = AT = U0TD LT 由分解惟一性, 即得 U0T=L A=L D LT 记 又因为det(Ak)＞0,(k=1,2,…,n), 故 于是对角阵D还可分解 其中 为下三角阵,令L=L1，定理得证。 将A=LLT展开，写成 按矩阵乘法展开，可逐行求出分解矩阵L的元素，计算公式是对于i=1,2,…,n j=i+1, i+2,…,n 这一方法称为平方根法,又称乔累斯基(Cholesky)分解,它所需要的乘除次数约 为数量级,比LU分解节省近一般的工作量。 例3.9 平方根法求解方程组 解: 因方程组系数矩阵对称正定,设A= ,即： 由Ly=b解得 由 解得 由此例可以看出，平方根法解正定方程组的缺点是需要进行开方运算。为避免开方运算，我们改用单位三角阵作为分解阵，即把对称正定矩阵A分解成 的形式，其中 为对角阵，而 是单位下三角阵,这里分解公式为 据此可逐行计算 运用这种矩阵分解方法,方程组Ax=b即 可归结为求解两个上三角方程组 和 其计算公式分别为 和 求解方程组的上述算法称为改进的平方根法。这种方法总的计算量约为 ，即仅为高斯消去法计算量的一半。 3.5 追赶法 在数值计算中,有一种系数矩阵是三对角方程组 简记为Ax=f,A满足条件 （1） （2） （3） 用归纳法可以证明，满足