清华第五版数值分析第5章课件.pptVIP

证明：高斯消去法能进行下去的充要条件是 而消元过程不会改变顺序主子式 的值。所以 解的唯一性又要求 高斯列主元消去法的适用范围 定理：系数矩阵非奇异，这高斯列主元消去法可行。 高斯消元法每一消元步骤相当于对系数矩阵左乘一个行初等矩阵，具体为 一 矩阵的三角分解 若矩阵A有分解：A=LU，其中L为下三角阵，U为上三角阵，则称该分解为A的LU分解。 若矩阵A有分解A=LU，则解线性方程组 就等价于求解 当系数矩阵A相同，有很多不同的右端向量时，用此方法很有效。如求 即求方程 的解，用此法很合适。 下面看什么情况下，可对系数矩阵做三角分解。 三角分解的存在唯一性定理 定理 引理：若 则 存在唯一性证明 证明：存在性：高斯消去法能进行下去，则三角分解存在。 唯一性：设 则 下三角阵=上三角阵，所以 所以 =LU 求解 时，A 和 的误差对解 有何影响？ ? 设 A 精确， 有误差　 ，得到的解为 ，即 绝对误差放大因子 又 相对误差放大因子 §5 线性方程组的性态和解的误差分析 §6 Error Analysis for . ? 设 精确，A有误差 　 ，得到的解为 ，即 Wait a minute … Who said that ( I + A?1? A ) is invertible? (只要? A充分小，使得 是关键 的误差放大因子，称为 A的条件数，记为cond (A) , 越 则 A 越病态， 难得准确解。 大 例：Hilbert 阵 cond (H2)? = 27 cond (H3)? ? 748 cond (H6)? = 2.9 ? 106 cond (Hn)?? ? as n ? ? 注：现在用Matlab数学软件可以很方便求矩阵的状态数! 定义2: 设线性方程组的系数矩阵是非奇异的,如果cond(A) 越大,就称这个方程组越病态.反之,cond(A)越小,就称这个方程组越良态. 一般判断矩阵是否病态，并不计算A?1，而由经验得出。 ? 行列式很大或很小（如某些行、列近似相关）； ? 元素间相差大数量级，且无规则； ? 主元消去过程中出现小主元； ? 特征值相差大数量级。 ? 近似解的误差估计及改善： 设　　 的近似解为　 ，则一般有 cond (A) 误差上限 ? 改善方法(1) ： Step 1:　　　　　近似解 Step 2: Step 3: Step 4: 若　 可被精确解出，则有 就是精确解了。 经验表明：若 A 不是非常病态（例如： ），则如此迭代可达到机器精度；但若 A 病态，则此算法也不能改进。 ? 改善方法(2) ? 对方程组进行预处理，即适当选择非奇异对角阵D，C,使求解 Ax=b 的问题转化为求解等价方程组 DAC[C-1x]=Db,且使DAC 的条件数得到改善。 ?用双精度进行计算，以便改善和减轻病态矩阵的影响。 二 矩阵的直接三角分解 直接计算 A 的 LU 分解(例) （续） 1) 2) 3) 一般计算公式 计算量与 Gauss 消去法同. LU 分解求解线性方程组 a11 a12 a13 … a1n u11 u12 u13 … u1n a21 a22 a23 … a2n l21 u22 u23 … u2n a31 a32 a33 … a3n l31 l32 u33 … u3n an1 an2 an3 … ann ln1 ln2 ln3 … unn … … … ... … … … … ... ... (1) (3) (5) (2n-1) (2) (4) (6) (2n) 紧凑格式 ①旧元素减去左边行与 顶上列向量的点积 ②计算行不用除法 ③计算列要除主对角元 注： 例 求矩阵的Doolittle分解 ? ? ? ? 例:用矩阵的直接三角分解法解方程组