天津理工大学数值计算第三章 戴敏.ppt

计算机科学与工程系 计算机科学与工程系 数值计算方法 第三章 线性方程组的数值解法 许多科学技术问题要归结为解含有多个未知量x1, x2, …, xn的线性方程组 方程组的矩阵形式 AX = b 线性方程组的克莱姆解法 如果矩阵A非奇异，即A的行列式det A≠0，根据克莱姆法则，方程组有唯一的解 其中△表示det A，△i表示A中第i列换成b后所得的行列式 例：已知二元线性方程组 方程组的解为 对于较高阶的情况，用这种方法求解是不现实的，一个n阶行列式有n!项，每一项又是n个数的乘积 例如求解一个10阶线性方程组，必须求出11个10阶行列式的值，计算量为 Md0＝10！＝3628800 Md＝11×Md0＝11×10！对于n阶线性方程组，其计算工作量为 (n+1)! 项 线性方程组的解法 直接法 将方程组的形式，根据某些计算规律，解析成另外一种形式（例如将系数矩阵变为单位矩阵），直接求得线性方程组的解 直接法要求存储的数据很多，因此必须应用存储容量大的计算机 在计算过程中，会产生积累误差，因此计算精度是有限的 直接法的计算时间是固定的、有限的 当线性方程组的维数不大，计算机的容量又较大时，适宜采用直接解法 间接解法 假定一组线性方程组的解，将这组解带入方程组（迭代形式的线性方程组），然后得到一组新解，再把这组新解代入方程组，然后又得到一组更新的解，反复这样代入求解下去，直到这组解满足一定精度的时候，才结束计算 间接法在计算过程中，只需将系数矩阵中的非零元素存入，因此可使计算机的存储容量小一些 采用反复迭代求解，逐步逼近最佳解，每次迭代都可看成是重新开始计算，因此在计算机中无积累误差 迭代法计算时间较长 所以对于维数较大的线性方程组，在电子计算机容量有限的情况下，选用间接解法是适宜的 3.1 高斯消去法 高斯消去法是一种直接法 在这种方法中，只包含有限次的四则运算，计算有限步，能直接得到方程组的精确解 在实际计算过程中由于舍入误差的存在和影响，只能求得近似解 高斯消去法的计算步骤就是“消元－回代”的两个过程 3.1.1 顺序高斯消去法 消元过程 已知线性方程组AX＝b，将A与b合并为一个n×(n+1)的增广矩阵，然后对这个增广矩阵进行初等行变换，将它变为一个上三角矩阵 例：对三阶线性方程组的消元过程 3.1.1 顺序高斯消去法 把方程（I）乘（ ）后加到方程（II）上去，把方程（I）乘（ ）后加到方程（III）上去，即可消去方程（II）、（III）中的x1，得同解方程组 将方程（II）乘（ ）后加于方程（III），得同解方程组 3.1.1 顺序高斯消去法 由回代过程得 x3 = 2 x2 = 8 x1 = -13 一般形式的线性方程组的消元过程 如果a11 ? 0，从其它n –1个方程中消x1，使它变成如下形式 3.1.1 顺序高斯消去法 3.1.1 顺序高斯消去法 如果 ，则以 为主元素，又可以把方程组化为 3.1.1 顺序高斯消去法 继续消元，重复同样的手段，第k步所要加工的方程组是 3.1.1 顺序高斯消去法 设 ，从其它（n - k）个方程中消xk，消元公式为 3.1.1 顺序高斯消去法 按照上述步骤进行n－1次后，将原方程组加工成下列形式 称为各次消元的主元素； 称为各次消元的乘数，主元素所在的行称为主行 3.1.1 顺序高斯消去法 例：四阶方程组的消元过程 3.1.1 顺序高斯消去法 回代过程 下三角方程组 如果aii? 0，i = 1, 2,…, n，则解为 3.1.1 顺序高斯消去法 上三角方程组 若aii ? 0, i = 1, 2,…, n，则解为 3.1.1 顺序高斯消去法 高斯消去法的回代公式 高斯消去法分为消元过程与回代过程，消元过程将所给方程组加工成上三角形方程组，再经回代过程求解 3.1.1 顺序高斯消去法 高斯消去法的流程图 3.1.1 顺序高斯消去法 算法的特点 按消元规则进行运算后，对角线以下元素为0，故运算中对于对角线以下元素不作计算，以减小计算量 对角线以下元素对回代求解没有影响，故将乘数放在该处，即 以节省存储单元 3.1.1 顺序