第五节 线性方程组的直接方法(modified by WQ)aaaa.ppt

§5.3 高斯主元素消去法 交换原则：通过方程或变量次序的交换，使在对角线位置上获得绝对值尽可能大的系数作为akk(k)，称这样的akk(k) 为主元素，并称使用主元素的消元法为主元素法 根据主元素选取范围分为：列主元素法、行主元素法、全主元素法 主元素法的意义 全主元素法不是按列选主元素，而是在全体待选系数中选取，则得全主元素法。 例5.3 用全主元素法解下列线组 计算m21=-19/40=0.475，m31=4/40=0.1 (5)- m21(4), (6)- m31(4) 消去x2 得 保留有主元素的方程 5.3.2 列主元素法 列主元素法就是在待消元的所在列中选取主元，经方程的行交换，置主元素于对角线位置后进行消元的方法。 例5.4 用列主元素法解下列线性方程组 (5)- m21(4), (6)- m31(4)得 保留有主元素的方程 为了研究线性方程组近似解的误差估计 和迭代法的收敛性, 有必要对向量及矩阵的 “大小”引进某种度量----范数的概念。向量范 数是用来度量向量长度的,它可以看成是二、 三维解析几何中向量长度概念的推广。用Rn 表示n维实向量空间。 定义5.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* ||，其相对误差可表示成 定义5.5（矩阵的范数）如果矩阵 的某个非负的实值函数 ，满足 定义5.7（矩阵的谱半径）设 的特征 值为 , 称 为A的 谱半径。 例 5.12 计算方阵 的三种常用范数 例5.12 计算方阵 的三种范数 解 先计算 所以 ,从而 定理5.8.1 设A为n阶方阵, 则对任意矩阵范数 都有 证: 设 为A的特征值，x是 对应于的特征向 量,则 x=Ax。两端取范数并依据其性质 得 由于x≠0，故 ，所以 5.8 误差分析 5.8.1 方程组的性态 在建立方程组时，其系数往往含有误差（如观测误差或计算误差），就是说，所要求解的运算是有扰动的方程组，因此需要研究扰动对解的影响。 例5.13 考察方程组 和 上述两个方程组尽管只是右端项有微小扰动, 但解大不相同, 第1个方程组的解是 第2个方程组的解是 。这类方程组称为病态的。 定义5.8 A或b的微小变化(又称扰动或摄动)引起方程组Ax=b解的巨大变化，则称方程组为病态方程组，矩阵A称为病态矩阵。否则方程组是良态方程组，矩阵A也是良态矩阵 为了定量地刻画方程组“病态”的程度，要对方程组Ax=b进行讨论，考察A（或b）微小误差对解的影响。为此先引入矩阵条件数的概念。 定义5.9 （矩阵条件数）设A为非奇异矩阵，称 为矩阵A条件数。 再确定U的第k行元素与L的第k列元素,对于k=2,3, …,n计算：① 计算U的第k行元素 （j=k,k+1,…,n） ② 计算L的第k列元素 （i=k,k+1,…,n） 利用上述计算公式便可逐步求出U与L的各元素 求解 Ly=b , 即计算: 求解 Ux=y , 即计算： 显然, 当 时, 解Ax=b直接三角分解法计算才能完成。设A为非奇异矩阵, 当 时计算将中断或者当 绝对值很小时，按分解公式计算可能引起舍入误差的积累，因此可采用与列主元消去法类似的方法，对矩阵进行行交换，则再实现矩阵的三角分解。 用直接三角分解法解Ax=b大约需要 次乘除法。 例6.8 用三角分解法解方程组