31第一节 高斯消元法.pptVIP

数学学院 信息与计算科学系 直接方法 　　若在计算过程中没有舍入误差，经过有限次的算术运算，可求得方程组 Ax=b （|A|?0）的精确解的方法，称为直接方法。它主要用于解系数矩阵是低阶稠密矩阵的线性方程组。 　　这里的基本思想是：把一般线性方程组化为等价的系数矩阵为三角矩阵的线性方程组来求解。（为什么？） 第三章 解线性方程组的数值方法 解为 例如： 上三角矩阵所对应的线性方程组 下三角矩阵所对应的线性方程组 计算量（乘除法的主要部分）为 n2/2. 解为 　　因为计算量少，故需要化为这种形状. 记 Ax=b 为 A(1)x=b (1) ，即 第一节 消元法 一、 Gauss顺序消元—按自然顺序进行的消元法 1、基本思想：将其化为与之等价的三角形方程组来求解，具体步骤如下： 若a11(1)?0，记 li1= ai1(1) /a11(1), i=2 , 3 , …, n , 第i个方程加上第一个方程乘以- li1得 若用矩阵表示，相当于用一个初等矩阵 L1 左乘 A(1) 和 b(1) , 即 L1A(1) =A(2) , L1b(1)=b(2)， 其中 若 akk(k)?0, 记 lik= aik(k) /akk(k), i=k+1, k+2,…,n 第i个方程加上第k个方程乘以 -lik，用矩阵表示， b(k) , 分别得到 LkA(k) =A(k+1) , Lk b(k)=b(k+1)， 其中 相当于用一个初等 矩阵Lk 左乘A(k)和 方程组的的解为： 前面称为消元的过程，后面称为回代的过程。 最后得 即 A(n)x=b (n) L 为单位下三角矩阵 令 以上n-1步消元过程就是 U为上三角矩阵 即有 A=LU 其中 L 为单位下三角矩阵 , U为上三角矩阵。从而得求矩阵行列式的计算公式 其中lik= aik(k) /akk(k), k=1,2,…,n ,i=k+1 , k+2 , … , n,计算量为 n3/3。 2. 高斯消元法可用矩阵表示为： 例1 用高斯消元法求解方程组 解 用第一个方程消去后两个方程中的 得 再用第1个方程消去第2个方程中的 得 最后，经过回代求得原方程组的解为 例2 解方程组 解：消元 回代得解为 引理 A的主元素 (k=1,2,…,n) 的充要条件是矩阵A的各阶顺序主子式不为零，即 定理1 如果在消元过程中A的主元素 (k=1,2,…,n) ,则可通过高斯消元法求出Ax=b 的解. 二、高斯消元法的条件 定理2 Ax=b 可用高 斯消元法求解的充分必要条件是：系数矩阵 A 的各阶顺序主子式均不为零. 数学学院 信息与计算科学系