NC-2-1.pptVIP

Numerical Calculation College of Science，Chongqing University 第二章 线性方程组的解法 掌握线性方程组的直接解法； 掌握线性方程组的迭代解法； 理解向量与矩阵的范数、条件数的概念。 第二章 线性方程组的解法 第二章 线性方程组的解法 §2.1 消元法 2.1.1 高斯消元法 2.1.1 高斯消元法 2.1.2 列主元Gauss消元法 2.1.2 列主元Gauss消元法 2.1.2 列主元Gauss消元法 2.1.3 Gauss-Jordan消元法 2.1.3 Gauss-Jordan消元法 2.1.3 Gauss-Jordan消元法 2.1.3 Gauss-Jordan消元法 §2.2 矩阵三角分解 2.2.1 Doolittle分解 2.2.1 Doolittle分解 2.2.1 Doolittle分解 2.2.1 Doolittle分解 2.2.2 Crout分解 2.2.2 Crout分解 2.2.3 平方根法 2.2.3 平方根法 2.2.3 平方根法 2.2.4 追赶法 2.2.4 追赶法 2.2.4 追赶法 2.2.4 追赶法 2.2.4 追赶法 2.2.4 追赶法 2.2.4 追赶法 就一般方程而言 （4） 以 ， 表示经过k步消元后相应的系数和右端项 令 （5） 先将方程（1）中 的系数化为1，并从其余方程中 消去 ，使之变为： 第一步 （6） 为此，先用 遍除方程 ，即进行运算 （7） 然后， 即令 （8） 第二步 再将方程（6）中 的系数化为1，并从其余方程中 消去 。 设经过k-1 步以后，方程组为： （9） 第k 步 再将方程（9）k中 的系数化为1，并从其余方程中 消去 。 （10） 为此，要实施运算 （11） 和 （12） 这样，经过 n 步以后，所给方程组最终被加工成： （13） 即为所给方程组的解 约当消去法的总计算量为： 约当消去法在加工过程中，用老系数计算新系数，而进 一步的计算中只用到新系数，故老系数可以不再保留。 故可动态地计算（11）、（12） 高斯-若当消元法解方程组并不比高斯消元法优越,但用于矩阵求逆是适宜的,实际上它是初等变换方法求逆的一种规范化算法: 例3 已知 ，求 所以 Gauss-Jordan法的Matlab程序见课本P13页。 设 ，若存在下三角矩阵L和一个上三角阵U使得A=LU，则称LU为矩阵A的三角分解。为保证三角分解的唯一，应该使分解式规范化。若限定L是单位下三角矩阵（对角线为1的下三角矩阵），则称该分解Doolittle分解，若限定U为单位上三角矩阵（对角线为1的上三角矩阵），则称该分解为Crout分解。 定理2.2 若矩阵A的各阶顺序主子式非零，则A可进行Doolittle分解和Crout分解，且这种分解是唯一的。 当A进行LU分解后， 就容易解了．这时 可改写为LUx=b ，令Ux=y，有Ly=b ．即 等价于 。即 分解为解两个三角形方程组，三角形方程组是极易求解的． 矩阵A的Doolittle分解可以用不带行交换Gauss消元法完成，其中L的下三角元素就是由消元运算的乘数因子组成，U是消元运算结束后的上三角矩阵。但也可以用矩阵乘法原理推出计算公式，当然结果是完全一致的。设： 则由矩阵的乘法原理有计算公式： 1、计算U 的第一行： ； 2、计算L 的第 一列， ； 3、对 进行以下的计算： （2.14） （2.15） 4、解 （2.16） 5、解 （2.17） 例2.3 用LU分解方法求解方程组