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

科学计算与MATLAB 主讲：唐建国 中南大学材料科学与工程学院 2010.12 第九讲　解线性方程组的直接解法 内容提要 引言 Gauss消元法 列主元素消元法 误差分析 MATLAB的线性方程组求解函数1 小结 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 将数据逐个代入,可得五个方程的方程组，求解该线性方程组即可得行星轨道方程。 对一般线性方程组: 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.1 基本思想：逐步消去未知元，将方程组化为与其等价的上三角方程组求解。 2.2 分两步： 第一步: 消元过程，将方程组消元化为等价的上三角形方程组; 第二步: 回代过程， 解上三角形方程组,得原方程组的解。 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 回代过程 解上三角方程组 上三角方程组的一般形式为： 其中a11…ann≠0 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消元法中的约化主元可以表示为 例 用高斯消元法求解方程组 矩阵的三角分解： 例 对矩阵A＝ 作LU分解。 解 由Gauss消去法可得， m21 0， m31 2， m32 －1。 故 A＝ ＝LU 3、Gauss列主元素消元法 （用四位浮点数计算，精确解舍入到4位有效数字为：x1* -0.4904,x2* -0.05104,x3* 0.3675） 其中， m21 -1.000/0.001