数值分析03-线性方程组直接解法

《数值分析03-线性方程组直接解法》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第三章线性方程组直接解法3-1阜师院数科院第三章线性方程组直接解法第三章目录§1.Gauus消元法§2.主元素法2.1引入主元素法的必要性2.2列主元素法2.3全主元素法2.4解三对角方程组的追赶法§3.矩阵分解法3.1Gauss消去法的矩阵形式3.2矩阵的三角分解3.3直接三角分解法§4.平方根法与改进的平方根法§5.矩阵求逆§6.方程组的性态和条件数2阜师院数科院第三章线性方程组直接解法线性方程组的概念在科学研究和工程技术中所提出的计算问题中，线性方程组的求解问题是基本的，常见的，很多问题如插值函数，最小二乘数据拟合

2、，构造求解微分方程的差分格式等，都包含了解线性方程组问题，因此，线性方程组的解法在数值计算中占有较重要的地位。设n阶线性方程组：其矩阵形式为：Ax=b(3-2)其中：3阜师院数科院第三章线性方程组直接解法求解Ax=b，曾经学过高斯（Gauss）消元法，克莱姆（Cramer）法则，矩阵变换法等，但已远远满足不了实际运算的需要，主要体现两个方面：一是运算的快速和准确，其次是方程组的个数增大时的计算问题。如何建立能在计算机上可以实现的有效而实用的解法，具有极其重要的意义，我们也曾指出过，Cramer法则在理论上是绝对正确的，但

3、当n较大时，在实际计算中却不能用。线性方程组的概念（续）如果线性方程组Ax=b的系数行列式不为零，即det(A)0，则该方程组有唯一解。4阜师院数科院第三章线性方程组直接解法线性方程组的数值解法解线性方程组的数值方法大致分为两类：请注意：由于在计算中某些数据实际上只能用有限位小数，即不可避免地存在着舍入误差的影响，因而即使是确解法，也只能求到近似解。直接法在求解中小型线性方程组（≤100个），特别是系数矩阵为稠密型时，是常用的、非常好的方法。直接法：指假设计算过程中不产生含入误差，经过有限步四则运算可求得方程组准确解的

4、方法。2.迭代法：从给定的方程组的一个近似值出发，构造某种算法逐步将其准确化，一般不能在有限步内得到准确解。这一章介绍计算机上常用的直接法，它们都是以Gauss消元法为基本方法，即先将线性方程组化为等价的三角形方程组，然后求解。5阜师院数科院第三章线性方程组直接解法§1Gauss消元法Gauss消元法是最基本的一种方法，下例说明其基本思想:例1解线性方程组：解：消去x1，进行第一次消元：首先找乘数，以-12乘第一个方程加到第二个方程，以18乘第一个方程加到第三个方程上可得同解方程组：6阜师院数科院第三章线性方程组直接解法

5、例1（续）上述Gauss消元法的基本思想是：先逐次消去变量，将方程组化成同解的上三角形方程组，此过程称为消元过程。然后按方程相反顺序求解上三角形方程组，得到原方程组的解，此过程称为回代过程。再消一次元得：二次消元后将方程化为倒三角形式，然后进行回代容易解出：x3=3,x2=2,x1=1。我们的目的，是要总结归纳出一般情况下的n阶线性方程组的消元公式和回代求解公式，从而得到求解n阶线性方程组的能顺利在计算机上实现的行之有效的算法。7阜师院数科院第三章线性方程组直接解法Gauss消元法的基本步骤1（4阶）为能更清楚地得到算法

6、，下面以4阶线性方程组为例总结求解步骤，并且很容易地可推广至一般的n阶线性方程组。8阜师院数科院第三章线性方程组直接解法可以检查，分别以li1乘第一个方程加到第i个方程上可以完成第一次消元，得同解方程组：变化以后的方程组系数及右边的常数项可总结出如下的计算公式：完成第一次消元之后的方程组记为：A(2)x=b(2)Gauss消元法的基本步骤2（4阶）9阜师院数科院第三章线性方程组直接解法Gauss消元法的基本步骤3（4阶）以方程组中第i个方程减去第二个方程乘li2(i=3,4)，完成第二次消元。上标为3的系数和右端项可

7、由下面公式计算：10阜师院数科院第三章线性方程组直接解法Gauss消元法的基本步骤4（4阶）第三步：消元（4阶方程组需进行3次消元）将上述A(3)X=b(3)中最后一个方程中的x3消为零：然后可回代求解：由于A(4)为上三角形，所以可按变量的逆序逐步回代求原方程组的解：上述消元、回代求解过程很容易推广到一般的n阶线性方程组。经过上述消元步骤，得到同解的上三角形方程组：A(4)x=b(4)11阜师院数科院第三章线性方程组直接解法Gauss消元法的消元过程1、2（n阶）一般地，设n阶方程组：消元过程为：将上方程组中第i个方程

8、减去第2个方程乘以li2(i=3,…,n)，完成第二步消元。12阜师院数科院第三章线性方程组直接解法Gauss消元法的消元过程3（n阶）第k步：设第k1步消元后得原方程组的同解方程组为：第k步消元后同解方程组中上标为k+1的元素的计算公式见下屏13阜师院数科院第三章线性方程组直接解法照此消元下去，完成n1次消元后