计算方法 第七章 线性代数方程组的直接解法.ppt

7.2.4 高斯主元素消去法 交换原则：通过方程或变量次序的交换，使在对角线位置上获得绝对值尽可能大的系数作为akk(k)，称这样的akk(k) 为主元素，并称使用主元素的消元法为主元素法 根据主元素选取范围分为：列主元素法、行主元素法、全主元素法 主元素法的意义 全主元素法不是按列选主元素，而是在全体待选系数中选取，则得全主元素法。 例7.3 用全主元素法解下列线组 计算m21=-19/40=0.475，m31=4/40=0.1 (5)- m21(4), (6)- m31(4) 消去x2 得 保留有主元素的方程 7.2.4.1 列主元素法 列主元素法就是在待消元的所在列中选取主元，经方程的行交换，置主元素于对角线位置后进行消元的方法。 例7.4 用列主元素法解下列线性方程组 (5)- m21(4), (6)- m31(4)得 保留有主元素的方程 7.4 平方根法 在应用数学中，线性方程组大多数的系数矩阵为对称正定这一性质，因此利用对称正定矩阵的三角分解式求解对称正定方程组的一种有效方法，且分解过程无需选主元，有良好的数值稳定性。矩阵的这一特性使它的三角分解也有更简单的形式，从而导出一些特殊的解法，如平方根法与改进的平方根法。 A对称：AT=A A正定：A的各阶顺序主子式均大于零。即 由Doolittle分解，A有唯一分解 A=LDLT分解，既适合于解对称正定方程组，也适合求解A为对称，而各阶顺序主子式不为零的方程组 而对A=LLT只适合于对称正定方程组 LDLT分解法解对称方程组所含的乘除法运算约为n3/4次， LLT分解法解对称正定方程组约需乘除法n3/6次。 LDLT法虽然增加了计算量，但避免了开方运算且扩大了使用范围，优点是明显的。 三对角方程组求解的追赶法 三对角方程组求解的追赶法 三对角方程组求解的追赶法 追赶法 在线性代数方程组的数值解法中，经常需要分析解向量的误差，需要比较误差向量的“大小”或“长度”。那么怎样定义向量的长度呢？我们在初等教学里知道，定义向量的长度，实际上就是对每一个向量按一定的法则规定一个非负实数与之对应，这一思想推广到维线性空间里，就是向量的范数或模。 一、内积的定义及性质 二、向量的长度及性质 定义7.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* ||，其相对误差可表示成 定义7.5（矩阵的范数）如果矩阵 的某个非负的实值函数 ，满足 相容性 （1）矩阵范数与矩阵范数的相容:‖AB‖≤‖A‖‖B‖ 条件数的定义 条件数的性质 条件数的计算 条件数的计算（续1） 条件数的计算（续2） 条件数的计算 “病态”方程的经验判断 “病态”方程的处理 对角占优矩阵 定理 7.13 (A的扰动对解的影响) 设A非奇异， 若Ax=b≠0，且 则 证 ：见p154 此定理说明：若 充分小且 时，解X的相对误差限约为矩阵A的相对误差的||A-1|| ||A||倍。同样有 ||A-1|| ||A||越大，解的相对误差限也就越大，近似程度就越差。 我们还可证明更为一般的结论： 当方程组的系数矩阵A非奇异和常数项b为非零向量时，且同时有扰动δA，δb，满足 ，若x和x+δx分别是方程组 Ax=b 及 的解 则 证 ：见p155 从以上几个定理可以看出，从某种意义上讲，||A-1|| ||A||刻划了方程组的解对原始数据的灵敏度，或者说刻划了方程组的“病态”程度，为此先引入矩阵条件数的概念。 定义7.9 （矩阵条件数）设A为非奇异矩阵，称