线性方程组直接方法(迭代法).ppt

第六章 解线性方程组的直接方法;;;§2 Gauss消去法;;;;;;;;;;二 矩阵的三角分解;;;;;;三 Guass消去法的计算量;;用Gauss消去法解AX=b时，设A非奇异，但可能出现 这时必须进行行交换的Gauss消去法，但在实际计算中即使 但其绝对值很小时，用 作除数，会导致中间结果矩阵 的元素数量级严重增长和舍入误差的扩散，使最后的结果不可靠。;; 这个例子还告诉我们，对同一个数值问题，用不同的计算方法，得到的精度大不一样，一个计算方法，如果用此方法的计算过程中舍入误差得到控制，对计算结果影响较小，称此方法为数值稳定的；否则，如果用此计算方法的计算过程中舍入误差增长迅速，计算结果受舍入误差影响较大，称此方法为数值不稳定。因此，我们解数值问题时，应选择和使用数值稳定的计算方法，否则，如果使用数值不稳定的计算方法去解数值计算问题，就可能导致计算失败。;一 完全主元素消去法;;（1）选主元素：选取 使;（4）回代求解。 经过上面的过程，即从第1步到第n-1步完成选主元，交 换两行，交 换两列，消元计算，原方程组约化为 ：;二 列主元素消去法(1); 完全主元素消去法是解低阶稠密矩阵方程组的有效方法，但完全主元素方法在选主元时要花费一定的时间。现介绍一种在实际计算中常用的部分选主元，（即列主元）消去法。列主元消去法，即是每次选主元时，仅依次按列选取绝对值最大的元素作为主元素，且仅交换两行，再进行消元计算。; 第k步计算如下： 对于k=1，2……..，n-1做到（4） （1）按列选主元，即确定 使;解：精确解为（舍入值）：; 本例是具有舍入的4位浮点数进行运算，所得的计算解还是比较准确的。 下面是完全主元素消去法框图（图6-1）;;消元计算 ;①;②—; 与原方程组准确解 比较，可知，本题用3位有效数字计算的列元素法是相当准确的。 大量实践表明：列主元法为解线性方程组的精确方法。;;列主元素消去法步骤：;6, 消元计算： aij - mik aij→ aij (i,j=k+1,…,n) bi - mik bk→ bi (i=k+1,…,n);⑴;则由排列阵性质（左乘矩阵是对矩阵进行行变换。）; Gauss消去法总是消去对角线下方的元素。现考虑一种修正，即消去对角线下方和上方的元素。这即为Gauss-Jordan（G-J）消去法。; 在第k步计算时，考虑对上述矩阵的第k行上、下都进行消元计算（k=1,2…n）;上述过程全部执行完后有：;定理;例;;故;4 Gauss消去法的变形;;;;（B) 选主元的三角分解法;;;二 平方根法;;;;;;;;例;在一些实际问题中常有解三对角线性方程组Ax=f的问题，即：; 对于具有条件（2）的方程组（1），我们介绍下面的追赶法求解。追赶法具有计算量少，方法简单，算法稳定的特点。; 现假设定理对n-1阶的满足条件（2）的三对角矩阵成立，求证对满足条件（2）的n阶三对角矩阵定理亦成立。由条件 ，则利用消元法的第1步有：;证明：由于A是满足条件（2）的n阶三对角阵。因此A的任一个顺序主子式 亦满足条件（2）的n阶三对角矩阵。由上一个定理;根据这一结论以及三角分解定理知，这种矩阵A可进行三角分解： A=LU。在这里特别的有：;;上面已经验证(5)对i=1成立。现设(5)对i-1成立。求证(5)对i成立。 由假设;由此可得求解三对角线性方程组的追赶法:; 我们将计算 及 的过程称为追的过程。计算方程组的解 的过程称为赶的过程。追赶法求解 Ax=f 仅需5n-4次乘除运算。工作量较小。;（2）计算 ：;需要对 （n维 列向量空间）中向量及 中的矩阵引进某种度; ||x+y|| ||y||