线性方程组求解__直接法详解.pptVIP

解线性方程组的直接法 引言 消元法：Gauss消元法 消元法：列主元素消元法 LU分解：求解三对角方程组的追赶法 误差分析 2.1 引言 小行星轨道问题： 天文学家要确定一小行星的轨道,在轨道平面建立以太阳为原点的直角坐标系.在坐标轴上取天文测量单位(一天文单位为地球到太阳的平均距离：9300万哩),对小行星作5次观察, 测得轨道上5个点的坐标数据如下： x 5.7640 6.2860 6.7590 7.1680 7.4800 y 0.6480 1.2020 1.8230 2.5260 3.3600 椭圆的一般方程: a1x2 + a2xy + a3y2 + a4x + a5y + 1 = 0 将数据逐个代入,可得五个方程的方程组，求解该线性方程组即可得行星轨道方程。 2.1 引言 对一般线性方程组: A x = b, 其中 (1)? 输入系数矩阵A和右端向量b； (2)计算系数矩阵A的行列式值D，如果D=0，则输出错误信息，结束，否则进行第(3)步； (3) 对k=1,2,···,n，用b替换A的第k列数据，并计算替换后矩阵的行列式值Dk； (4) 计算并输出x1 = D1 / D，x2 = D2 / D，····,xn=Dn/D， 结束。 但克莱姆法则只适用于低阶方程组，高阶方程组工作量太大，故一般用数值方法求解。数值方法分两类： 1. 直接法 2. 迭代法 2.2 Gauss消元法 基本思想：逐步消去未知元，将方程组化为与其等价的上三角方程组求解。 分两步： 第一步: 消元过程，将方程组消元化为等价的上三角形方程组; 第二步: 回代过程， 解上三角形方程组,得原方程组的解。 消元过程(化一般方程组为上三角方程组) 以四阶为例： 其系数增广矩阵为： 消元过程(化一般方程组为上三角方程组) 第一轮消元： 计算3个数: [m21 m31 m41]T = [a21 a31 a41]T / a11 用-m21乘矩阵第一行后加到矩阵第二行; 用-m31乘矩阵第一行后加到矩阵第三行; 用-m41乘矩阵第一行后加到矩阵第四行; 其系数增广矩阵变为： 消元过程(化一般方程组为上三角方程组) 第二轮消元： 计算2个数：[m32 m42]T = [a32(1) a42(1)]T / a22(1) 用-m32乘矩阵第二行后加到矩阵第三行; 用-m42乘矩阵第二行后加到矩阵第四行; 其系数增广矩阵变为： 消元过程(化一般方程组为上三角方程组) 第三轮消元： 计算: m43=a43(2)/a33(2) 用-m43乘矩阵第三行后加到矩阵第四行; 其系数增广矩阵变为： 消元过程(化一般方程组为上三角方程组) 其对应的上三角方程组为 消元过程(化一般方程组为上三角方程组) 若对于一般的线性方程组Ax=b，其消元过程的计算公式为： （k=1,2,…,n-1） 回代过程(解上三角方程组) 上三角方程组的一般形式为： 其中a11…ann≠0 回代过程(解上三角方程组) 工作量计算： 消去过程： “÷”：第k步，n-k次，共 (n-1)+(n-2)+……+1=n(n-1)/2 “×”：第k步，(n-k)(n-k+1)次，共 (n-1)n+(n-2)(n-1)+……+1×2= (n3-n)/3 总工作量： s1=n(n-1)/2+ (n3-n)/3 回代过程： “÷”：n “×”：1+2+……+(n-1)=n(n-1)/2 总工作量：s2=n+ n(n-1)/2= n(n+1)/2 故Gauss消元法的总工作量为： s=s1+s2= n2+(n3-n)/3 克莱姆法则求解的工作量为： “×”：（n+1个n阶行列式的值）(n+1)(n-1)n! “÷”：n 故总工作量为： [(n+1)(n-1)] n!+n 当n=6时, Gauss消元法工作量为106 ;而克莱姆法则求解工作量为25206。 定理： 约化的主元素ak+1,k+1(k) ≠ 0 (k=0,1,···,n-1)的充分必要条件是 矩阵A的各阶顺序主子式不为零。即 推论: 如果A的顺序主子式Dk ≠ 0 (k=1,···,n-1)，则Gauss消元法中的约化主元可以表示为 例 用高斯消元法求解方程组 矩阵的三角分解： 定理：设A为n阶矩阵，若A的顺序主子式Di ≠0（i=1,2,…n-1）,则A可分解为一个单位下三角矩阵L和一个上三角矩阵U的乘积，且这种分解是唯一的。 由Gauss消元过程可推得 L＝ U即为Gauss消元后所得的上三角方程的系数矩阵。 例 对矩阵A＝