数值分析解线性代数方程组直接解法.pptVIP

第二节 高斯消元法及其计算机实现;线性代数方程组的一般形式;数值分析解线性代数方程组直接解法;MATLAB实现: x=A\b; 数值求解方法有以下三条途径（三种框架） 直接法：利用Gauss消元或矩阵分解，通过有限次运算 可求出精确解。; ; A b; 三角形方程组包括上三角形方程组和下三角形方程组，是最简单的线性方程组之一。上三角方程组的一般形式是： ;数值分析解线性代数方程组直接解法;;;数值分析解线性代数方程组直接解法;for j = n : – 1 : 2 b ( j ) = b ( j ) / A ( j , j ); b (1: j - 1 ) = b (1: j - 1 ) - b ( j ) *A (1: j - 1 , j ) ; end b ( 1 ) = b ( 1 ) / A ( 1 ,1 );;数值分析解线性代数方程组直接解法;数值分析解线性代数方程组直接解法; 高斯消元法是一个古老的直接法,由它改进得到的选主元法,是目前计算机上常用于求低阶稠密矩阵方程组的有效方法,其特点就是通过消元将一般线性方程组的求解问题转化为三角方程组的求解问题。;解：;数值分析解线性代数方程组直接解法;将方程组Ax=b的系数矩阵与右端项合并为; ;数值分析解线性代数方程组直接解法;数值分析解线性代数方程组直接解法;进行到第k步消元时;数值分析解线性代数方程组直接解法; 用回代过程求解上三角方程组，即可得解向量 ( x1*,x2*, …,xn* )T.;求解的全过程包括两个步骤：消元和回代 1 . 顺序消元;function X=gauss(A,b) %Input—A is an n×n nonsingullar matrix % ---b is an n×1 matrix %Output—X is the solution to the system AX=b;function X=gauss(A,b) %Input—A is an n×n nonsingullar matrix % ---b is an n×1 matrix %Output—X is the solution to the system AX=b; 高斯消元法的计算量分析 高斯消元法的乘除总运算分析为 消元次数k 消元乘法次数 消元除法次数 回代乘除法次数 1 n(n-1) n-1 2 (n-1)(n-2) n-2 . k (n-k+1)(n-k) n-k . n-1 2*1 1 n(n+1)/2 高斯消元法的计算量为 ; 消元法是解线性方程组的基本方法，具有计算简 单的优点，但有时由于主元过小，使得计算???果严重 失真,实际中常采用选主元高斯消元法。;例4:讨论下面方程组的解法; 选主元基本思想 k步消元前,应该在第k列元素 (i=k,…,n)中找出第一个出现的绝对值最大者,例如 ， 再把第ik个方程与第k个方程组进行交换,使 成为主元.我们称这个过程为选主元.由于只在第k列元素中选主元,通常也称为按列选主元. ;再交换第k,ik两个方程和第k,jk列,使 成为主元. 称这个过程为完全选主元. ;;算法 列主元高斯消元法解线性方程组 Ax = b;数值分析解线性代数方程组直接解法;