第二章线性代数方程组的直接解法解析.ppt

§2 线性代数方程组的直接解法 §2.1 引论 §2.2 Gaoss消去法 2.2.1 顺序消去与回代过程 2.2.2 顺序消去能实现的条件 2.2.3 矩阵的三角分解 2.2.4 列主元素消去法 §2.2.1 顺序消去与回代过程 在上面的消元过程中,应特别注意到元素(称之为主 元素).因为如果主元素为0,消元过程将无法进行下去.如 果主元 的绝对值较小,作为除数也不合适.所以这种 方法是不恰当的. 另外,我们总是通过乘数 来消元,从而称 为乘 数. §2.2.2 顺序消去实现的条件 §2.2.3 矩阵的三角分解 引理3.1 用m阶单位下三角方阵L左乘m×n阶矩阵A所得矩阵B,则A与B的p阶主子式相等. §2.2.4 列主元素消去法 解决问题的方法是:在系数矩阵中按列选取绝对值 最大的元素作为主元素进行消元.以后每一步都作一次 这样的调整.直至将系数矩阵化为上三角矩阵. 称这种消元法为: 列主元消元法. 另外还有: 全主元消元法. §2.3 直接三角分解方法 2.3.1 Doolittle分解方法 2.3.2 三对角方程组的追赶法 2.3.3 对称正定矩阵的Cholesky分解、 平方根法 §2.4 矩阵的条件数与病态方程组 2.4.1 扰动方程组、病态现象 2.4.2 矩阵的条件数与扰动方程组的 误差分析 2.4.3 病态方程组的解法 §2.4.1 扰动方程组、病态现象 定义7:设A是非奇异矩阵,若Cond(A)远大于1,则称方程组AX=b是病态的,否则称方程组是良态的. 对于先前的例子 §2.4.3 病态方程组的解法 一个方程组是否病态，由A的性质决定，与用什么数值 解法无关。对于病态的方程组。数值求解要小心进行，否 则可能得不到所要求的精确度。 解病态方程组，可以采用高精度运算，例如双倍或更多 倍字长的运算，使由于误差放大而损失若干有效数位之 后，还能保留一些有效位。 谢谢！ 数学与统计学院 主页 * 数 值 分 析 武 汉 大 学 数学与统计学院基础数学系 刘丁酉 liudingyou487@163.com 主页 数学与统计学院 主页 2.1 引论 2.2 Gauss消去法 2.3 直接三角分解方法 2.4 矩阵的条件数与病态方程组 数学与统计学院 数学与统计学院 主 AX = b 数学与统计学院 主页 数学与统计学院 主页 数学与统计学院 主页 消元公式 我们通过初等变换得到一个与原方程组等价的方程组 回代公式 数学与统计学院 主页 数学与统计学院 主页 定理2.2.1 设 ，若对k=1,2,……，n，A的顺序主子式 ，则每步消去过程的对角元 定理2.2.1的逆定理也成立，即若 则 数学与统计学院 主页 定义3.1 m×n阶矩阵A的前p行,前p列所组成的p(1≤p ≤min(m,n))阶方阵称为A的p阶顺序主子矩阵,A的p阶顺序主子矩阵的行列式被称为p阶顺序主子式 . 简称p阶主子式. 数学与统计学院 数学与统计学院 用矩阵的分块理论可设: 其中 分别为L,A,B的p阶顺序主子矩阵.依条件有: 数学与统计学院 因此有 即有 引理得证 数学与统计学院 定理3.1 若方阵A的各阶主子式都不为零,则存在唯一的单位下三角矩阵L和上三角矩阵U,使得LU=A 对A进行初等变换 数学与统计学院 若记 则 所以,当n=2时,定理成立. 设当n=k时,定理成立. 则,当n=k+1时,有: 数学与统计学院 于是 数学与统计学院 由引理3.1知 的各阶主子式不为零.(注)再由归纳法假设知:存在单位下三角矩阵和上三角矩阵 ,使 . 又 所以 数学与统计学院 设A还有另外的三角分解,即有 于是有 进一步有 由于此式左边上三角矩阵,右边为下三角矩阵,所 以必有 进而有 数学与统计学院 此定理是矩阵分解的重要定理,称作矩阵的LU分 解,也叫杜利特尔分解. 它是用直接法解线性方程组 的理论基础.此外,可利用LU分解来计算A的行列式的 值.即: 数学与统计学院 （4） （5） 进而,由 LY=b, UX=Y, 有: 数学与统计学院 主页 数学与统计学院 主页 ? 杜立特分解法 /* Doolittle Factorization —— LU 分解的紧凑格式 反复计算, 很浪费哦 …… 通过比较法直接导出L 和 U 的计算公式。 思路 §2.3.1 Doolittle分解方法 数学与统计学院