实际问题中解线性方程组的经典解法.docVIP

在实际问题中解线性方程组的经典解法 直接法与三角形方程组的求解分析 线性方程组求解问题在许多科学计算问题中都会遇到，如应力分析、电学网络、自由振动问题等。在计算机数值方法的课程中，2线性方程组求解在样条插值、数据拟合的最小二乘法以及常微分方程边值问题中都要用到.产生的线性方程组的类型有很多，如按系数矩阵含零元素多少分类，有稠密和稀疏（零元素占80%以上）线性方程组之别；如按阶数的高低分类，有高阶（阶数在1000阶以上）和低阶之别；如按系数矩阵的形状和性质分类，又有对称正定、三对角线对角占优等之别，因为在电子计算机上求解，必须要考虑算法的计算复杂性以及算法的数值稳定性问题。所以针对不同类型的线性方程组有不同的解法，但是，基本的方法可归结为两大类，即为直接法和迭代法。 本章介绍的经典解法，都把原方程组化为一个或者两个三角形方程组来求解，主要包括Gauss4消去法和它的变形----直接三角分解法 设有线性方程组 （1，1） 其中 根据线性代数知识，当时，方程组（1.1）的解存在且唯一，对增广矩阵施行行初等变换，化为上三角形矩阵，同时化为，这时与增广矩阵相应的线性方程组为上三角形方程组 （1.2） 其中 设 则（1.2）的解为 （1.3） 它便是原方程组（1.1）的解，实现上述求解过程的方法称为Gauss消去法 如果方程组（1.1）的系数矩阵可分解为两个形式简单的三角形矩阵和的乘积，即 (1.4) 若为下三角形矩阵，则为上三角形矩阵，反之亦然 从而求解的问题转化为解三角形方程组 (1.5) 和 （1.6） 其中 则为下三角形方程组，它的第个方程为 （1.7） 假定按的顺序解得 （1,8） 上三角形方程组的第个方程为 （1.9） 假定按的顺序求解得 （1, 10） 直接法的求解过程，在计算过程中无舍入误差的前提下，都可以经有限步算数运算而得到精确解。由于计算机的字长有限，初始数据取浮点数以及运算过程都不断地产生舍入误差，这些误差的传播和积累，会影响计算精度，所以，如何避免舍入误差的增长是设计算法时必须考虑的问题。 直接法通常需要存储系数矩阵的全部元素，当方程组的阶数很高时，需要相当大的内存空间，因此，在算法设计上应当注意节省内存，比如，对称矩阵可以只存其下三角形部分于一维数组中。 综上所述，如果矩阵非奇异，总可以通过带有行交换或不带行交换的消元过程，将化为非奇异上三角形矩阵因此，回代求解过程（1,3）也可以进行到底，但是，在实际应用中，常常难于事先判断系数矩阵的奇异性，因此，需要进一步考虑当为奇异矩阵时计算过程可能发生的情况。一是消元过程的某一步找不到非零的，于是计算中断；二是虽然消元过程能进行到底，但，使回代求解过程无法进行下去，因此，在计算设计中必须考虑到上述两种可能发生的情形，此时应在算法设计中给出计算中断的信息。 Gauss列主元素消去法 在Gauss消元过程中位于矩阵的主对角线位置上的元素称为主元素，因为在计算解的分量时，都做除数，但应当避免用小的数做除数，即避免用较在数量级上相对小的主元做除数，以防止舍入误差的扩大，降低解得精度，下面的例子说明“小主元”对解得精度的影响。 例1用Gauss消去法解线性方程组 用8位十进制尾数的浮点数计算 解 类似的算出我们看到，由于小主元做除数，使得行乘数的数量级变得很大，这样在计算时，即与在计算机上做代数和时，由于阶码升为使尾数右移变成了机器零，这里揭示了数量级相对于分子数做除数，会使Gauss消去法在计算机有限字长的环境下造成较大舍入误差的原因。 经第一步消元得 经第二步消元得 这一结果表明原方程组有无穷多个解，在计算机上，将会给出“无唯一解”的信息，停止计算，但实际上，由于，所以原方程组有唯一解，是计算过程的舍