第一章 解线性方程组的直接法.pptVIP

解线性方程组的直接法 例题 例题 例题 * 理学院 University of Shanghai for Science and Technology College of Science 上海理工大学 实际中，存在大量的解线性方程组的问题。很多数值方法到最后也会涉及到线性方程组的求解问题： 如求样条插值的M和m关系式，解曲线拟合的法方程， 求矩阵特征值的反幂法等问题。 对线性方程组： 或者： 我们有Gram法则：当且仅当 时，有唯一的解，而且解为： 但在实际工程运算中Gram法则不能用于计算方程组的解， 如 n＝100，1033次/秒的计算机要算10120年，这是可想象，也 是不现实的。 一般的，在工程领域解线性方程组的方法可以分为2类： ①直接法：准确，可靠，理论上得到的解是精确的; ②迭代法：速度快，但有误差. 本章讲解直接法 我们知道，虽然方程组系数矩阵为一般比较“完整”的矩阵时， 解方程组非常困难，但如果方程的系数矩阵为下面有3种时， 我们可以直接求出其解： ① n次运算 ② （n＋1）n/2次运算 一、消去法 对角矩阵 下三角矩阵 ③ （n＋1）n/2次运算 问题是：我们如何将系数矩阵为一般形式的方程组，化为具有以上特殊系数矩阵的方程组？ 下三角矩阵 我们知道，对方程组作如下的变换，解不变 ①交换两个方程的次序； ②一个方程的两边同时乘以一个非0的数； ③一个方程的两边同时乘以一个数，加到另一个方程。 因此，对相应的增广矩阵(A，b) 作初等行变换方程组解不变，即 ①交换矩阵的两行； ②某一行乘以一个非0的数； ③某一行乘以一个数，加到另一行； 消元法就是对增广矩阵作上述的行变换，使之变为我们已知的 3种类型之一，而后求解的方法。 思路 首先将A化为上三角阵，再回代求解 。 = 1.高斯(Gauss)消元法 第一过程为消去过程——将系数矩阵化为上三角矩阵： 运算量： (n-1)?(1+n) 步骤分为两个过程 第一步：第i行＋第1行×(-ai1)/a11, i=2,3,…,n. 运算量： (n-2)?(1+n -1) =(n-2)n 第二步：第i行＋第2行×(-ai2(2))/a22(2), i=3,4,…,n. 第k步： 类似的做下去，我们有： 运算量： (n－k) ?(1＋n－k＋1)=(n－k)(n－k＋2). n－1步以后，我们可以得到变换后的矩阵为： 第二过程为回带过程——对系数矩阵为上三角矩阵方程组，逐步求解： ………… 因此，总的运算量为： 加上 解上述上三角阵的运算量(n+1)n/2，总共为： 用Gauss消去法解方程组: 1消元：对增广矩阵进行初等行变换 2回带：解以上矩阵为增广矩阵的方程组： 注意到，计算过程中 处在分母的位置，因此整个计算过程要保证它不为0 所以，Gauss消元法的可行条件为： 就是要求 A 的所有顺序主子式均不为 0，即 因此，有些有解的问题，不能用Gauss消元求解 另外，如果某个 很小的话，会引入大的误差 例：单精度解方程组 /* 精确解为 和 */ 8个 8个 用Gaussian 消元法计算： 8个 小主元可能导致计算失败。 在 Gauss 消元第 k 步之前，做如下的事情： 若 交换k行和j行 行的交换，不改变方程组的解，同时又有效地克服了Gauss消元地缺陷 例： ? 2、列主元消元法 用Gauss消去法解方程组 1消元：对增广矩阵进行初等行变换 2回带：解以上矩阵为增广矩阵的方程组： 选主元 消元 选主元 消元 将在Gauss消元第k步，变为 将该行上三角的部分也变为0 最后变为一个对角阵。 它地运算次数比Gauss消元多。使用于计算多个系数一样地方程组，如 X，B均为矩阵 3、Jordan消元法 用Jordan消去法解方程组 1消元：对增广矩阵进行初等行变换 2回带：解以上矩阵为增广矩阵的方程组： Gauss消元法的第k步： 从矩阵理论来看，相当于左乘矩阵 二、三角分解法 分析Gauss消元法的过程： 因此，整个Gauss消元法相当于左乘了一个单位下三角阵 所以有 L为单位下三角阵，U为上三角阵 因此 我们可以通过2次反代过程求解方程组 最终使系数矩阵化为上三角矩阵即A=L-1U L为下三角，U为单位上三角 比较第1行： 比较第1列： 1、Doolittle分解——直接三角分解 比较第2行： 比较第2列： 比较第k行： 比较第k列： K-1次 K-1＋1次 分解过程完毕，加上两次反代过程 总运算量为： 存储在矩阵的原来位