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

第二节 高斯消元法及其计算机实现;线性代数方程组的一般形式;;MATLAB实现: x=A\b; 数值求解方法有以下三条途径（三种框架） ; ; A b;三角形方程组包括上三角形方程组和下三角形方程组，是最简单的线性方程组之一。上三角方程组的一般形式是： ;;;;;; 高斯消元法是一个古老的直接法,由它改进得到的选主元法,是目前计算机上常用于求低阶稠密矩阵方程组的有效方法,其特点就是通过消元将一般线性方程组的求解问题转化为三角方程组的求解问题。;解：;;将方程组Ax=b的系数矩阵与右端项合并为; ;;;进行到第k步消元时;; 用回代过程求解上三角方程组，即可得解向量 ( x1*,x2*, …,xn* )T.;求解的全过程包括两个步骤：消元和回代;;;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;;;例5 用列主元消去法解???程组; 由此回代,得 x=(1.9272,-0.69841,0.90038)T与精确解 x=(1.9273,-0.69850,0.90042)T相比较是比较准确的. ;二版习题 P113----6（1）