第四章 解线性方程组的直接方法.ppt

4.2.4 高斯主元素消去法 交换原则：通过方程或变量次序的交换，使在对角线位置上获得绝对值尽可能大的系数作为akk(k)，称这样的akk(k) 为主元素，并称使用主元素的消元法为主元素法 根据主元素选取范围分为：全主元素法、列主元素法、行主元素法 例4.4 用全主元素法解下列线组 计算l21=-19/40=0.475，l31=4/40=0.1 (5)- l21(4), (6)- l31(4) 消去x2 得 保留有主元素的方程 2)列主元素法 列主元素法就是在待消元的所在列中选取主元，经方程的行交换，置主元素于对角线位置后进行消元的方法。 例4.5 用列主元素法解下列线性方程组 (5)-l21(4),(6)-l31(4)得 保留有主元素的方程 为了研究线性方程组近似解的误差估计 和迭代法的收敛性, 有必要对向量及矩阵的 “大小”引进某种度量----范数的概念。向量范 数是用来度量向量长度的,它可以看成是二、 三维解析几何中向量长度概念的推广。用Rn 表示n维实向量空间。 定义4.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* ||，其相对误差可表示成 定义4.5（矩阵的范数）如果矩阵 的某个非负的实值函数 ，满足 记笔记 或 例4.10 证明对任意同维向量x, y 有 证： 即 例4.11 设x=(1, 0, -1, 2)T, 计算 解: =1+0+|-1|+2=4 二维向量空间的几何意义： 定理4.7 对于任意向量x ,有 证: ∵ ∴ 即 当 p→∞， ∴ 定义4.4( 向量序列的极限 ) 设 为 中的 一向量序列， , 记 。如果 (i =1,2,…, n)，则称 收敛于向量 ，记为 定理4.8 （向量范数的等价性）设 为 上任意两种向量范数, 则存在常数C1,, C20, 使得对任意 恒有 （证:P133） 则称 是 上的一个矩阵范数(或模) 4.3 矩阵三角分解法 矩阵三角分解法是高斯消去法解线性方程组的一种变形解法 4.3.1 矩阵的LU分解原理 应用高斯消去法解n阶线性方程组Ax=b, 经过n步消元之后, 得出一个等价的上三角型方程组A(n)x=b(n), 对上三角形方程组用逐步回代就可以求出解来。 上述过程可通过矩阵分解来实现。即将非奇异阵A分解成一个下三角阵L和一个上三角阵U的乘积 。 将非奇异阵A分解成一个下三角阵L和一个上三角阵U的乘积 A=LU 称为对矩阵A的三角分解，又称LU分解 也称为杜利特尔（Doolittle）分解。其中 detLk=1, Lk invertible, k=1,2,…, n-1 k=1,2,…, n-1 上三角阵 4.3.2 用三角分解法解方程组 单位下三角阵 单位下三角阵 上三角阵 LU分解定理 A的各阶顺序主子式 单位下三角阵 上三角阵 存在性不必证,见前面的分析,仅证唯一性(采用反正法) L,U,L1U1都可逆 单位下三角 上三角阵 证明: L A的第一行 U的第一行 A的第一列 L的第一列 ? U Lrk=0, k=r+1,…, n; Lrr=1 设已知U的前r-1行，L的前r-1列，求U的第r行，L的第r列 ? A的第r行 U的第r行 Ukr=0, k=r+