4 第五章 解线性方程组的直接方法.pptVIP

这个解对于具有舍入的3位浮点数进行计算，是一个很好的结果。这个例子告诉我们，在采用高斯消去法解方程组时，小主元可能导致计算失败，故在消去法中应避免采用绝对值很小的主元素。方法1计算失败的原因，是用了一个绝对值很小的数作除数，乘数很大，引起约化中间结果数量误差很严重增长，再舍入就使得计算结果不可靠了。这个例子还告诉我们，对同一个数值问题，用不同的计算方法，得到的精度大不一样，一个计算方法，如果用此方法的计算过程中舍入误差得到控制，对计算结果影响较小，称此方法为数值稳定的；否则，如果用此计算方法的计算过程中舍入误差增长迅速，计算结果受舍入误差影响较大，称此方法为数值不稳定。因此，我们解数值问题时，应选择和使用数值稳定的计算方法，否则，如果使用数值不稳定的计算方法去解数值计算问题，就可能导致计算失败。对一般矩阵方程组，需要引进主元的技巧，即在高斯消去法的每一步应该选取系数矩阵或消元后的低阶矩阵中的绝对值最大的元素作为主元素，保持乘数，以便减少计算过程中的舍入误差对计算解的影响。（4）回代求解。经过上面的过程，即从第1步到第n-1步完成选主元，交换两行，交换两列，消元计算，原方程组约化为：完全主元素消去法是解低阶稠密矩阵方程组的有效方法，但完全主元素方法在选主元时要花费一定的时间。现介绍一种在实际计算中常用的部分选主元，（即列主元）消去法。列主元消去法，即是每次选主元时，仅依次按列选取绝对值最大的元素作为主元素，且仅交换两行，再进行消元计算。完全主元素消去法在选主元素时要花费多机器时间。下面介绍另一种常用的方法即列主元素消去法，它仅考虑依次按列选主元素，然后换行使之变到主元素位置上，再进行消元计算。设用列主元素消去法解Ax=b已完成k-1步计算。即有：列主元素消去法步骤：设Ax=b。对于具有行交换的列主元素消去法，消元结果冲掉A，乘数mij冲掉aij，计算解X冲掉常数项b,行列式存在det.1,1→det,k=1（对k=1—n-1做2—7步。）3，若aik,k=0,则0→det.计算终止。5，计算乘数mik,mik=aik/akk→aik.(i=k+1—n)(|mik|≤1)6,消元计算：aij-mikaij→aij(i,j=k+1,…,n)bi-mikbk→bi(i=k+1,…,n)下面用矩阵运算来描述列主元素法。⑴其中Lk的元素满足|mik|≤1（k=1,2,……n-1）,由⑴得：则由排列阵性质（左乘矩阵是对矩阵进行行变换。）这表明：对Ax=b应用列主元素法相当于对（A|b）先进行一系列行交换后对PAX=Pb再应用Gauss消去法。在实际计算中我们只能在计算过程中做行的变换。有结论：定理：（列主元素三角分解定理）若A为非奇异性矩阵，则存在排列矩阵P使PA=LU。其中L为单位下三角阵，U为上三角阵。L存放在A的下三角部分，U存放在A的上三角部分。由整数型数组IP（n）记录可知P的情况。（例子见前）Gauss消去法总是消去对角线下方的元素。现考虑一种修正，即消去对角线下方和上方的元素。这即为Gauss-Jordan（G-J）消去法。在第k步计算时，考虑对上述矩阵的第k行上、下都进行消元计算（k=1,2…n）2，换行（当ik≠k）。交换[A,b]第k行与第ik行元素。3，计算乘数mik=-aik/akk(i=1—n,i≠k),mkk=1/akk(mik可存放在aik的单元中)（|mik|≤1）上述过程全部执行完后有：在一些实际问题中常有解三对角线性方程组Ax=f的问题，即：对于具有条件（2）的方程组（1），我们介绍下面的追赶法求解。追赶法具有计算量少，方法简单，算法稳定的特点。定理：设有三对角线性方程组Ax=f，且A满足条件（2），则A为非奇异矩阵。证明：用归纳法证明。显然n=2时，有：根据这一结论以及三角分解定理知，这种矩阵A可进行三角分解：A=LU。在这里特别的有：由此可得求解三对角线性方程组的追赶法:因此追赶法中不会出现中间结果迅速增长和全入误差产生积累现象。这就是矩阵的Frobenius范数,明显它满足上面定义。在计算上，（1），（2）比较容易，而（3）比较困难，但（3）在理论分析上十分有用。下面考察方程组的数据A(或b)的微小误差对解的影响。即估计x-y。其中y为（A+δA）y=b的解。可见方程组常数项分量只有微小变化（1%），而方程组的解却有较大的变化。即方程组的解对常数项b很灵敏。这样的方程组称为病态方程组。定义:如果矩阵A或常数项b的微小变化引起方程组Ax=b的解的巨大变化。则称此方程组为“病态”方程组，矩阵A称为“病态”矩阵。否则称此方程组为“良态”方程组，矩阵A称为“良态”矩阵。矩阵的“病态”性质是矩阵本身的特征。下面我们希望找出刻划矩阵“病态”性质的量。先考察常数项b的微小误差对解的影响。由⑷δx=-(A+δA)-