第二章 线性方程组的数值解法ppt课件.ppt

《第二章 线性方程组的数值解法ppt课件.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第二章线性方程组的数值解法§2.1消元法§2.2直接分解法§2.3向量和矩阵的范数§2.4雅可比迭代§2.5高斯-赛德尔迭代§2.6松弛迭代直接法:经过有限次运算后可求得方程组精确解的方法(不计舍入误差!)（Gauss消去法及其变形、矩阵的三角分解法）迭代法：从解的某个近似值出发，通过构造一个无穷序列去逼近精确解的方法。（一般有限步内得不到精确解）直接法比较适用于中小型方程组。对高阶方程组，既使系数矩阵是稀疏的，但在运算中很难保持稀疏性，因而有存储量大，程序复杂等不足。迭代法则能保持矩阵的稀疏性，具有计算简单，编制程序容

2、易的优点，并在许多情况下收敛较快。故能有效地解一些高阶方程组。线性方程组的数值解法AX=b（2.1）1．三角形方程组的解法---回代法（对角型方程组略）（2.2）（2.3）§2.1消元法2．顺序高斯消去法基本思想：通过消元将上述方程组化为三角形方程组进行求解。（1）消元过程其中第一步：若用乘第一行加到第i行中，得到第二步：若用…….……第k步：若用乘第k行加到第i行中，得到其中第n-1步：……（2）回代过程若则顺序高斯消去算法高斯消去法计算复杂度顺序Gauss消去法可执行的前提定理1给定线性方程组，如果n阶方阵的所有顺序

3、主子式都不为零，即则按顺序Gauss消去法所形成的各主元素均不为零，从而Gauss消去法可顺利执行。注：当线性方程组的系数矩阵为对称正定或严格对角占优阵时，按Gauss消去法计算是稳定的。顺序Gauss消去法计算过程中的akk(k)称为主元素，在第k步消元时要用它作除数,则可能会出现以下几种情况1、若出现akk(k)＝0，消元过程就不能进行下去。2、akk(k)≠0，消去过程能够进行，但若

4、akk(k)

5、过小，也会造成舍入误差积累很大导致计算解的精度下降。例2-1在四位十进制的限制下，试用顺序Gauss消去法求解如下方程

6、组此方程组具有四位有效数字的精确解为x1＝17.46，x2＝－45.76，x3＝5.546解用顺序Gauss消去法求解，消元过程如下经回代求解得x3＝5.546，x2＝100.0，x1＝－104.0和此方程组的精确解相比x3＝5.546，x2＝－45.76，x1＝17.46有较大的误差。对于此例，由于顺序Gauss消去法中的主元素绝对值非常小，使消元乘数绝对值非常大，计算过程中出现大数吃掉小数现象，产生较大的舍入误差，最终导致计算解x1＝－104.0和x2＝100.0已完全失真。为避免这种现象发生，可以对原方程组作等价变

7、换，再利用顺序Gauss消去法求解。写出原方程组的增广矩阵：针对第一列找出绝对值最大的元素，进行等价变换：求得方程的解为：x3＝5.546，x2＝－45.76，x1＝17.46精确解为：x3＝5.546，x2＝－45.76，x1＝17.46由此可见，第二种Gauss消去法的精度明显高于顺序Gauss消去法，我们称它为列主元Gauss消去法。列主元Gauss消去法与顺序Gauss消去法的不同之处在于：后者是按自然顺序取主元素进行消元前者在每步消元之前先选取主元素然后再进行消元3、列主元Gauss消去法计算步骤：1、输入矩阵

8、阶数n，增广矩阵A(n,n+1);2、对于(1)按列选主元：选取l使(2)如果，交换A(n,n+1)的第k行与第l行元素(3)消元计算：3、回代计算从§2.1中讨论可知，顺序Gauss消去法的消元过程是将增广矩阵[A，b]＝[A（1），b（1）]逐步化为矩阵[A（n），b（n）]。现在说明，在消元过程中，系数矩阵A＝A（1）是如何经矩阵运算化为上三角矩阵A（n）。即用矩阵运算的观点来看，消元的每一步计算等价于用一个单位下三角矩阵左乘前一步化得到的矩阵。§2.2直接分解法若，令，i＝2，3，…，n，得到下三角矩阵施行第一步

9、消元，我们得到若，令，i＝2，3，…，n，则有施行第二步消元，我们得到如此下去，施行第n－1步消元，得到(n)由此可见，在顺序Gauss消去法的过程中，系数矩阵A＝A（1）经过一系列单位下三角矩阵的左乘运算化为上三角矩阵A（n），即这时由得令容易验证则从顺序Gauss消去法的矩阵运算表示式可知，系数矩阵A可分解为一个单位下三角矩阵L和一个上三角矩阵U的乘积，即其中第一个方程组的系数矩阵为下三角矩阵，第二个方程组的系数矩阵为上三角矩阵，两个方程组都非常容易求解，具体求解结果如下：我们将A=LU称为矩阵A的三角分解,这时线性

10、方程组为：令则有对于由解得对于由求得可以看出对于方程组：只要对系数矩阵作了三角分解：由这个简单的计算过程可知，系数矩阵的三角分解很关键，如何进行三角分解更容易？下面介绍几种方法。通过如下两组公式很容易求解：注：（1）L为单位下三角阵而U为一般上三角阵的分解称为Doolittle分解（2）L为一般下三角阵而U为单位上三