第二章线性代数方程组.pptVIP

第二章;确定小行星轨道;a1xj2 + 2a2xjyj + a3 yj2 +2a4 xj + 2a5 yj + 1 = 0; 在科学计算中，经常需要求解含有n个未知量 的n个方程构成的线性方程组;克莱姆法则需{[(n+1)(n-1)] n!+n}次乘除法运算 (1)? 输入系数矩阵A和右端向量b;;线性方程组数值解法的分类：;线性方程组数值解法的分类：;2.1 高斯消元法;1．三角形方程组的解法---回代法;2.顺序Gauss（高斯）消元法是一种规则化的加减消元法。 ;例1 用Gauss消元法解方程组;这样，对于方程组;顺序Gauss消元法的消元过程可表述如下： ;其中 ;其中 ;增广矩阵[A（n），b（n）]对应如下上三角形方程组 ;对于等价方程组;首先写出增广矩阵;然后进行消元，采用公式;在编程计算时，最后的增广矩阵存放的元素是：;算法. ;顺序Gauss消元法可执行的前提;例2 用Gauss消元法求解方程组：;例3 用Gauss消元法求解线性方程组;3、列主元Gauss消元法;此方程组具有四位有效数字的精确解为;经回代求解得 x3＝5.546，x2＝100.0，x1＝－104.0;写出原方程组的增广矩阵：;求得方程的解为：x3＝5.546，x2＝－45.76，x1＝17.46;下面将列主元Gauss消元法的计算步骤叙述如下：; 3. 再进行第二步消元得到增广矩阵[A（3），b（3）]。按此方法继续进行下去，经过n－1步选主元和消元运算，得到增广矩阵[A（n），b（n）]，它对应的方程组;选列主元过程：;例5 用列主元法解;x3 = 6.2/6.2 = 1; 高斯消元法是一种顺序消元法。消元过程按方程和未知数的顺序依次作消元计算，缺乏一定的灵活性。计算中有除数为零时，算法便无法实现；即使所有的除数都不为零，可以计算出方程组的解，也不能保证计算结果有很高的精度。;m21 = 2/0.00001 = 2 × 105;m21 = 0.00001/2 = 0.5 × 10-5;2.2 矩阵的直接分解法;? 高斯消元法的矩阵形式 ;;A 的 LU 分解 ( LU factorization );第一个方程组的系数矩阵为下三角矩阵，第二个方程组的系数矩阵为上三角矩阵，两个方程组都非常容易求解。;定理2：（矩阵的三角分解）设A为n ? n实矩阵，如果 解AX = b用高斯消去法能够完成（限制不进行行的交 换，即 ），则矩阵A可分解 为单位下三角矩阵L与上三角知阵U的乘积。 A = LU 且这种分???是唯一的。 ;注： （1） L 为单位下三角阵而 U 为一般上三角阵的分解称为多利特尔（Doolittle） 分解； （2）L 为一般下三角阵而 U 为单位上三角阵的分解称为库朗（Courant ）分解。 ;? Doolittle分解法 ： ;根据 A=LU 有等式成立：;n;可以解得：;下面，我们对具体矩阵进行Doolittle 三角分解。;例7 利用Doolittle三角分解法分解矩阵;1;如果我们要求解方程组;例8：利用Doolittle三角分解方法解线性方程组;得到;例9： 用矩阵分解方法解 AX = b;不选主元的LU分解;? 平方根法 ： ;设A为对称正定矩阵，那么A的所有顺序主子式均大于零，根据定理2.2，存在惟一三角分解 A＝LU，即; 记 Ak（1≤k≤n）为A的 k 阶顺序主子阵，则det（Ak）为A的k阶顺序主子式。由上式，利用矩阵分块运算规则，容易验证 ;即： A=LDM ，其中 DM=U， M=D-1U。;当A＝AT为对称矩阵时，根据 A＝LDM, 得到; 如果对称正定矩阵A具有如下分解 A＝GGT，其中G为下三角矩阵，则称其为对称正定矩阵的 Cholesky（乔列斯基）分解。 ;Ly＝b LTx＝ y;比较对应元素：;当 i = j 时;于是，根据计算公式; 关于方程组 Ax=b , 如果对系数矩阵进行了平方根分解 A＝LLT，则将方程组化为： Ly＝b , LTx＝y; 于是，关于系数矩阵是对称正定矩阵的线性方程组Ax=b的求解，分两步进行：;例10 用平方根法求解对称正定方程组 ;得 y1＝2，y2＝3.5 ，y3＝1 ;? 解三对角方程组的追赶法;考虑如下形式的三对角方程组问题;该三对角方程组的增广矩阵 为:;第一步：第一行元素除以b1，取 ?1 = c1 / b1，y1 = f1 / b1，可得