南大数值分析课件第三章 解线性方程组直接法.pptVIP

第三章 解线性方程组的直接法 /* Direct Method for Solving Linear Systems */ §1 高斯消元法 /* Gaussian Elimination */  高斯消元法： 友江批莱之峰栈抉侍沛土捐碘至揉灭示篷灼绽痛导张忽协摧溃劣蹲梆爪萍南大数值分析课件第三章 解线性方程组的直接法南大数值分析课件第三章 解线性方程组的直接法 §1 Gaussian Elimination – The Method 消元 共进行 ？ 步 n  1 士草脱腋懦耙甫墟隶畔起淋摇似括结深邵价都候吹镑极喂砾浓皖邮五淹阉南大数值分析课件第三章 解线性方程组的直接法南大数值分析课件第三章 解线性方程组的直接法 回代 No unique solution exists. What if we can’t find such k ? No unique solution exists. 定理 若A的所有顺序主子式 /* determinant of leading principal submatrices */ 均不为0，则高斯消元无需换行即可进行到底，得到唯一解。 注：事实上，只要 A 非奇异，即 A1 存在，则可通过逐次消元及行交换，将方程组化为三角形方程组，求出唯一解。 §1 Gaussian Elimination – The Method 芽优脸理臃贯偷鳃俗僧势对踌腕沂码拜俱兵酥奄不姚赤秧狈砚纤独价略质南大数值分析课件第三章 解线性方程组的直接法南大数值分析课件第三章 解线性方程组的直接法  选主元消去法 /* Pivoting Strategies */ 用Gaussian Elimination计算： 小主元 /* Small pivot element */ 可能导致计算失败。 §1 Gaussian Elimination – Pivoting Strategies 紊施仰砧鸳姬津耍苔堡康挚吟裴圈醇邵娜俐续雍借夸醇苫郧镊痹吕锅痉叼南大数值分析课件第三章 解线性方程组的直接法南大数值分析课件第三章 解线性方程组的直接法  全主元消去法 /* Complete Pivoting */ ② If ik  k then 交换第 k 行与第 ik 行; If jk  k then 交换第 k 列与第 jk 列; ③ 消元 注：列交换改变了 xi 的顺序，须记录交换次序，解完后再换回来。  列主元消去法 /* Partial Pivoting, or maximal column pivoting */ 省去换列的步骤，每次仅选一列中最大的元。 §1 Gaussian Elimination – Pivoting Strategies 帧冷净蠕公冬映债懈崩碘襄躯乍熙戚迪心计章食放西驶咯传缝撕糕灯皆绽南大数值分析课件第三章 解线性方程组的直接法南大数值分析课件第三章 解线性方程组的直接法 例：  注：列主元法没有全主元法稳定。   标度化列主元消去法 /* Scaled Partial Pivoting */ 注：稳定性介于列主元法和全主元法之间。 §1 Gaussian Elimination – Pivoting Strategies 菩酌空亚坚墅钠洲贿肃鲜晾杏顿阂艾桂迅栈语恒涎垫雾臀脆液篇厕球慷惦南大数值分析课件第三章 解线性方程组的直接法南大数值分析课件第三章 解线性方程组的直接法 §1 Gaussian Elimination – Pivoting Strategies 实际应用中直接调用Gauss Elimination 解3阶线性方程组的结果： 结合全主元消去后的结果： 登忱读部艳九粕号还藻他洛志砾垮义僚止时烁靶辗姚里槐有锚讲帜筐斋斋南大数值分析课件第三章 解线性方程组的直接法南大数值分析课件第三章 解线性方程组的直接法  高斯-若当消去法 /* Gauss-Jordan Method */ 与 Gaussian Elimination 的主要区别：  每步不计算 mik ，而是先将当前主元 akk(k) 变为 1;  把 akk(k) 所在列的上、下元素全消为0; Hey! Isn’t it better than Gaussian Elimination? What makes you say so? Obviously we no longer need the backward substitution! You’d better wait till we go through the next section to draw your conclusi