数值线性方程组直接解法介绍.pptxVIP

第三 章;第三章目录; ; 求解Ax = b，曾经学过高斯（Gauss）消元法，克莱姆（Cramer）法则，矩阵变换法等，但已远 远满足不了实际运算的需要，主要体现两个方面：一是运算的快速和准确，其次是方程组的个数增大时的计算问题。如何建立能在计算机上可以实现的有效而实用的解法，具有极其重要的意义，我们也曾指出过，Cramer法则在理论上是绝对正确的，但当n较大时，在实际计算中却不能用。;线性方程组的数值解法;§1 Gauss消元法 ;例1（续）; 为能更清楚地得到算法，下面以4阶线性方程组为例总结求解步骤，并且很容易地可推广至一般的n阶线性方程组。; 可以检查，分别以?li1乘第一个方程加到第i个方程上可以完成第一次消元，得同解方程组：;Gauss消元法的基本步骤3（4阶）;第三步：消元（4阶方程组需进行3次消元） 将上述 A (3) X = b(3) 中最后一个方程中???x3消为零： ;Gauss消元法的消元过程1、2（n阶）;第k 步：设第k?1步消元后得原方程组的同解方程组为： ;照此消元下去，完成n?1次 消元后，可将原方程组化成 同解的上三角形方程组如下： ;Gauss消元法的回代过程（n阶）;Gauss消元法的计算量 ;Gauss法与Cramer法则的计算量比较;§2 主元素法 ;例2（续1）; 若在解方程组前，先交换方程的次序，如将（2-8）交换一行与二行改写成如下所示：;例2两种解法的误差分析;2.2 列主元素法;列主元素法;2.3 全主元素法 ;例6;2.4 解三对角方程组的追赶法 ; 首先进行顺序消元， 且每步将主元系数 化为1，将方程组化为：;追赶法举例;§3 矩阵分解法 ;第二次消元相当于用初等矩阵： ;经过n?1步消元后得到： ;杜利特尔（Doolittle）分解——LU分解 ;3.2 矩阵的三角分解 ;对A进行LU分解;一般情况下，可由：;对A进行LU分解的具体步骤;矩阵A的LU分解举例;3.3 直接三角分解法 ;直接三角分解法（续）;解线性方程组举例;三角分解法的几点说明;三角分解法的几点说明（续）;解线性方程组举例;§4 平方根法与改进的平方根法 ;4.1 平方根法 ;定理 2.2（续）;Choleskg分解1 ;Choleskg分解2;平方根法举例;4.2 改进的平方根法 ;改进的平方根法说明;改进的平方根法举例