第三章 线性方程组的解法(上) 1.ppt

* 式中 * 第k步：（k=1，2，…，n-1）继续上述消去过程，设第1步～第k-1步计算已经完成，得到与原方程组等价的方程组 记为A(k)x=b(k)；现进行第k步消元计算，设 ，计算乘数 用（mik）乘上式的第k个方程加到第i个方程，消去第i个方程的未知数xk，得到与原方程组等价的方程组 3.2.2 高斯消去法的算法构造（续） 第i行 第k行 （i=k+1，…，n） * 简记为A(k＋1)x=b(k＋1)；其中A(k＋1)，b(k＋1)元素计算公式为： 3.2.2 高斯消去法的算法构造（续） * 最后，重复上述约化过程，即k=1，2，…，n-1且设 （k=1，2，…，n-1）共完成n-1步消元计算，得到与原方程组（3.1）等价的三角形方程组 3.2.2 高斯消去法的算法构造（续） 消元过程 * 3.2.2 高斯消去法的算法构造（续） 第1步 在方程（3.4）的最后一个方程中解出xn，得 （2）回代过程 * 3.2.2 高斯消去法的算法构造（续） 第3步 依次继续下去，一般可得xk的计算公式 第2步 将xn的值代入式（3.4）的倒数第二个方程，解出xn-1，得 当k=1时，就完成了回代过程，得到所求的解。 (k=n-1, n-2, …, 1) * 将（3.1)约化为(3.4)的过程称为消元过程， （3.4）的求解过程称为回代过程 由消元过程和回代过程求解线性方程组的方法称为高斯消去法。 3.2.2 高斯消去法的算法构造（续） * ① 消元计算 k=1，2，…，n-1 3.2.2 高斯消去法的算法构造（续） 定理3.1 设线性方程组Ax=b，其中A∈Rn×n(Rn×n表示n阶方阵的集合）。 （1）如果 则可通过高斯消去法将Ax=b化为等价的上三角方程组，且有计算公式： * ② 回代计算 (2)如果A为非奇异矩阵时，可能有某 ，则在第k列存在有元素 ， 于是可能通过交换（A，b）的第k行和第ik行元素，将 调到（k，k）位置，然后再进行消元计算。 于是，在A为非奇异矩阵时，只要引进行交换，则高斯消去法可将 化为上三角方程组，且通过回代即可求得方程组解。 3.2.2 高斯消去法的算法构造（续） * 高斯消去法要求 第k步消元的主元素 3.2.2 高斯消去法的算法构造（续） 判断主元素 的一个充要条件： 定理3.2 对矩阵A=(aij)n×n消元时，主元素 的一个充要条件是矩阵的各阶顺序主子式 * 3.2.2 高斯消去法的算法构造（续） 即 * 3.2.3 高斯消去法算法分析 消元过程的计算量，高斯消去法消去过程分n-1步 第k步的计算工作量为： 1）计算乘数：需要作（n-k）次除法运算； 2）消元：需作(n-k)2次乘法运算和(n-k)2次加减法运算； 3）计算b(k)：需作（n-k）次乘法运算和(n-k)次加减法运算； k=1，2，…，n-1 （1）高斯消去法的计算量 （n-k）行,（n-k）列 （n-k）行 （n-k）列 * 乘除法运算 加减法运算 于是完成n-1步运算全部消元计算共需要作 求和公式 3.2.3 高斯消去法算法分析(续） * （2）回代计算： 共需要作 n(n+1)/2 乘除法运算，n,…,1 n(n-1)/2 加减法运算 ，n-1,…,1 用高斯消去法解 （其中A∈Rn×n）的计算量为共需作 乘除法运算 加减法运算 3.2.3 高斯消去法算法分析(续） * 在计算机上用高斯消去法解低阶稠密矩阵线性方程组时要注意几点： （1）要用一个二维数组A(n,n)存放系数矩阵A的元素，用一维数组b(n)存放常数项b分量。 （2）需要输入的数据： A ，b。 （3）约化的中间结果A(k)元素冲掉A元素，b(k)冲掉b。 例如，计算 （a）mik←-A(i,k)/ A(k,k)，（i=k+1，…，n）； （b）A(i,j)←A(i,j)+mik*A(k,j)，(i=k+1,…,n；j=k+1,…,n）; （c） b(i)←b(i) +mik* b(k)，（i=k+1，…，n）。 （4）如果不存在ik，使 ，输出方程没有唯一解的信息。 3.2.3 高斯消去法算法分析(