9-解线性方程组直接解法.pptVIP

科学计算与MATLAB;第九讲　解线性方程组的直接解法;内容提要;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， 结束。 克莱姆法则只适用于低阶方程组，高阶方程组工作量太大，故一般用数值方法求解。数值方法分两类： 直接解法：在计算过程中，如果所有运算都是精确的，在理论上，经过有限次运算就可以得到精确解，适用于变量较少的方程组。 迭代解法：近似解法，运算次数因要求的计算精度而变化。;迭代法也称辗转法，是一种不断用变量的旧值递推新值的过程。;2、Gauss消元法;2.3 Gauss消元的目的;2.4.1 消元过程(化一般方程组为上三角方程组);第一轮消元： 计算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）;2.4.2 回代过程(解上三角方程组);回代过程的计算公式：;2.5 工作量计算： 消去过程： “÷”：第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消元法中的约化主元可以表示为;例 用高斯消元法求解方程组; clear all; A=[2 3 4;3 5 2;4 3 30]; b=[6 5 32]; [L,U]=lu(A); x=U\(L\b) x = -13.0000 8.0000 2.0000;矩阵的三角分解：;将A 分解为单位下三角矩阵 L 和上三角矩阵 U的乘积的算法称为矩阵A的三角分解算法。;例 对矩阵A＝ 作LU分解。 解 由Gauss消去法可得， m21=0， m31=2， m32=－1。 故 A＝ ＝LU;如果已经有 A = L U 则 AX = b = L U X = b，;3、Gauss列主元素消元法;（用四位浮点数计算，精确解舍入到4位有效数字为：x1*=-0.4904,x2*=-0.05104,x3*=0.3675）;其中， m21=-1.000/0.001=-1000 m31=-2.000/0.001=-2000 m32=4001/2004=1.997 解为x1=-0.4000,x2=-0.09980,x3=0.4000 (x1*=-0.4904,x2*=-0.05104,x3*=0.3675) 显然，此解并不准确。;方法二：交换行，避免绝对值小的主元作除数。;其中， m21=0.5000