第5讲线性方程组的数值解法.pptVIP

计算方法 Computational Method Lecture 2 线性方程组的数值解法 计算方法 Computational Method Lecture 2 线性方程组的数值解法 计算方法 Computational Method 2.6 平方根法(square root method) 正定方程组——系数矩阵为对称正定矩阵的方程组 Lecture 2 线性方程组的数值解法 计算方法 Computational Method Lecture 2 线性方程组的数值解法 计算方法 Computational Method Lecture 2 线性方程组的数值解法 计算方法 Computational Method Lecture 2 线性方程组的数值解法 计算方法 Computational Method Lecture 2 线性方程组的数值解法 计算方法 Computational Method 自学：改进平方根法 Lecture 2 线性方程组的数值解法 2. 线性方程组的数值解法（linear equation） 计算方法 Computational Method 克莱姆（Cramer）法则 一个n阶行列式有n！项，每一项又是n个数的乘积！显然，让计算机用克莱姆法则求解线性方程组是不现实的！ Lecture 2 线性方程组的数值解法 2.1 高斯消去法（Gaussian elimination） 计算方法 Computational Method Lecture 2 线性方程组的数值解法 计算方法 Computational Method 归纳：首先逐步消去变元的系数，把原方程组化为系数矩阵为上三角阵的等价方程组，然后解系数矩阵为上三角阵的方程组。 消元过程——将原方程组化为系数矩阵为上三角阵的等价方程组 回代过程——求解系数矩阵为上三角阵的等价方程组 Lecture 2 线性方程组的数值解法 计算方法 Computational Method 消元过程 Lecture 2 线性方程组的数值解法 计算方法 Computational Method 回代过程 Lecture 2 线性方程组的数值解法 计算方法 Computational Method 高斯消去法的计算工作量 消去过程乘法次数： 消去过程除法次数： 消去过程乘法次数： 回代过程乘法次数： 回代过程加法次数： 回代过程除法次数： 乘除法总次数： 加法总次数： Lecture 2 线性方程组的数值解法 计算方法 Computational Method Lecture 2 线性方程组的数值解法 计算方法 Computational Method Lecture 2 线性方程组的数值解法 计算方法 Computational Method 第一组解显然是错误的，第二组解是准确解的具有2位有效数字的近似值。 导致第一组解错误的原因是消元过程中主元素很小 2.2 选主元素消去法（elimination with pivoting） 在系数矩阵中按列选取元素绝对值的最大值作为主元素，然后交换所在行与主元素所在行的位置，再按顺序消去法进行消元。 2.2.1 选列主元素消去法 Lecture 2 线性方程组的数值解法 计算方法 Computational Method Lecture 2 线性方程组的数值解法 计算方法 Computational Method 在系数矩阵中首先选取元素绝对值最大者作为主元素，交换第一行和此元素所在行，再交换第一列及此元素所在列，使主元素移到第一行第一列；在进行k－1次消元后准备进行第k次消元时，先从最后n-k+1行和最后n-k+1列所形成的矩阵中选出元素绝对值的最大者作为主元素，交换第k行和此元素所在行，再交换第k列和此元素所在列，然后再按在顺序消去法进行消元。 2.2.2 选全主元素消去法 Lecture 2 线性方程组的数值解法 计算方法 Computational Method Lecture 2 线性方程组的数值解法 计算方法 Computational Method 一般来说，列主元素消去法已有足够的精确度，但全主元素消去法往往会有更高的精确度，但全主元素消去花费在寻找绝对值为最大的元素以及交换行列上的时间较多，程序比较复杂。因此，实际应用时更偏向用列主元素消去法，理由为： ★只要系数矩阵的行列式不等于0，列主元素消去法总可以进行下去 ★lik的绝对值小于或等于1，有利于消元过程中对舍入误差的控制 ★列主元素消去法的选元过程只交换行而不交换列 讨 论 Lecture