高斯消去法-Read.PPT

* * } 第3章 线性代数方程组的解法 矩阵形式: AX=b,其中 当A非奇异时，det(A)≠0，方程组有唯一解。 本章解决的问题：线性代数方程组的基本解法。 n个未知量n个方程的线性代数方程组 } (3.1) } 两类解法： 直接解法: 经过有限步运算就能求得精确解的方法，但实际计算中由于舍入误差的影响，这类方法也只能求得近似解； 迭代解法:构造适当的向量序列，用某种极限过程去逐步逼近精确解。 例如：高斯消去法、三角分解法等 例如：雅可比迭代法、高斯赛德尔迭代法等 } 高斯消去法的基本思想： 先消元，即按一定的规律逐步消去未知量，将方程组化为等价的上三角形方程组；然后进行回代，即由上三角形方程组逐个求出 } 3.1 顺序高斯消去法 (消元过程按方程和未知量的自然顺序进行.) 3.1.1 顺序高斯消去法举例 3.1.2 一般情况的计算过程 } 3.1.1 顺序高斯消去法举例 例 用顺序高斯消去法求解方程组 解 消元过程： [ A，b ]= 对增广矩阵[ A，b ]做消元计算， } [ A，b ]= 2 4 0 1 1 0 2 2 0.5 1.5 0 1 3 -0.5 5.5 0 1 2 1 2 2 4 0 1 1 0 0 0 1 0.75 1.25 2 4 0 1 1 0 2 2 0.5 1.5 0 0 2 -0.75 4.75 0 0 0 1.125 -1.125 消元过程结束， 得上三角形方程组 0 2 2 0.5 1.5 0 2 - 0.75 4.75 } 回代过程： 求解上三角形方程组 所以方程组的解为 } 基本的顺序高斯消去法,其消元过程遵守如下规则: 依从左到右、自上而下的次序将主对角元下方的元素化为零。 不作行交换，也不用非零数乘某行。 第k列消元（将主对角元下方化零）时，将下面各行分别减去第k行的适当倍数，不作其他变换。 } 练习： 用顺序高斯消去法求解方程组 解 2 3 0 1 1 1 1 5 5 6 2 11 2 3 0 1 0 -1/2 1 9/2 0 -3/2 2 17/2 2 3 0 1 0 -1/2 1 9/2 0 0 - 1 - 5 得上三角形方程组，再回代得 x1= –1 , x2= 1, x3= 5 } 3.1.2一般情况的计算过程 为统一公式，把常数项 bk 换名为 ak n+1，(k=1,2,…,n) 先看4阶的情形 第1 列消元（设a11≠0），作 } 化为 其中 令 ，则计算公式合写为 } 第2列消元（设a22（1）≠0）， 化为 其中 即 } 第3 列消元（设a33（2）≠0）， 化为上三角形方程组 其中 ，计算公式写为 } 对于n 阶的一般情形 } 一般，对于n 阶的方程组，需要对前n - 1列进行消元。 设已进行了第1列至第k–1列消元， 将第i 行 减去第k 行的l i k 倍， (i=k+1,…,n) 第k列消元， } } 第n-1列消元后 得到与原方程组同解的上三角形方程组 } 若 ，则 一般由第k个方程得 消元过程结束。进行回代过程，由上三角形方程组依次求出xn ,xn-1 ,…,x1 。 点击可显示方程组→ 顺序高斯消去法的计算步骤： (1)消元过程。对 k=1,2,… ,n-1 计算 (2)回代过程。对 k=n,n-1,… ,1 计算 ①若akk=0,则算法失败, 结束计算；否则做下步。 ②对i=k+1,k+2,…,n计算 称为消元过程的主元素。其中 } } 问题1.消元过程能按顺序进行到底的充要条件是什么？ 问题2.在原方程组的系数矩阵中如何反映出这个条件呢? 答： A的k阶顺序主子矩阵Ak的行列式值是: 考察矩阵A的k阶顺序主子矩阵Ak ， 即A的左上角k行k 列的子矩阵，例如： 1 2 3 A= 4 5 6 则A1=(1), A2= 1 2 , A3=A 7 8 9 4