线性方程组的直接方法迭代法公开课一等奖优质课大赛微课获奖课件.pptVIP

第六章解线性方程组的直接方法1引言第1页第1页

第2页第2页

第3页第3页

§2Gauss消去法一Gauss消去法第4页第4页

第5页第5页

第6页第6页

第7页第7页

第8页第8页

第9页第9页

第10页第10页

第11页第11页

第12页第12页

第13页第13页

二矩阵的三角分解第14页第14页

第15页第15页

第16页第16页

第17页第17页

第18页第18页

第19页第19页

三Guass消去法的计算量第20页第20页

第21页第21页

用Gauss消去法解AX=b时，设A非奇异，但也许出现这时必须进行行互换Gauss消去法，但在实际计算中即使但其绝对值很小时，用作除数，会造成中间结果矩阵元素数量级严重增长和舍入误差扩散，使最后结果不可靠。例：设有方程组办法一：用Gauss消去法求解。（用含有舍入3位浮点数进行运算）解：方程组解3Gauss主元素消去法第22页第22页

办法二：用含有行互换Gauss消去法（避免小主元）这个解对于含有舍入3位浮点数进行计算,是一个较好结果。回代得解=1.00=0.00，与准确解比较，这结果很差。这个例子告诉我们,在采用高斯消去法解方程组时,小主元也许造成计算失败,故在消去法中应避免采用绝对值很小主元素。办法1计算失败原因,是用了一个绝对值很小数作除数,乘数很大,引起约化中间结果数量误差很严重增长,再舍入就使得计算结果不可靠了。第23页第23页

这个例子还告诉我们,对同一个数值问题,用不同计算方法,得到精度大不同,一个计算方法,假如用此方法计算过程中舍入误差得到控制,对计算结果影响较小,称此方法为数值稳定；不然,假如用此计算方法计算过程中舍入误差增加快速,计算结果受舍入误差影响较大,称此方法为数值不稳定。因此,我们解数值问题时,应选择和使用数值稳定计算方法,不然,假如使用数值不稳定计算方法去解数值计算问题,就可能造成计算失败。对普通矩阵方程组,需要引进主元技巧,即在高斯消去法每一步应当选取系数矩阵或消元后低阶矩阵中绝对值最大元素作为主元素,保持乘数,以便减少计算过程中舍入误差对计算解影响。第24页第24页

一完全主元素消去法第25页第25页

第26页第26页

（1）选主元素：选取使于是第k步计算：对于k=1，2，……..n-1做到（3）第k步选主元区域第27页第27页

（4）回代求解。通过上面过程，即从第1步到第n-1步完毕选主元，交换两行，交换两列，消元计算，原方程组约化为:第28页第28页

二列主元素消去法(1)第29页第29页

完全主元素消去法是解低阶稠密矩阵方程组有效办法,但完全主元素办法在选主元时要花费一定期间。现简介一个在实际计算中惯用部分选主元,（即列主元）消去法。列主元消去法,即是每次选主元时,仅依次按列选取绝对值最大元素作为主元素,且仅互换两行,再进行消元计算。第k步选主元区域设列主元消去法已经完毕第1步到第k-1步按列选主元，互换两行，消元计算得到与原方程组等价方程组：，其中第30页第30页

第k步计算下列：对于k=1，2……..，n-1做到（4）（1）按列选主元，即拟定使（3）假如，则互换第行第k行元素（2）假如=0，则A为奇异矩阵，停止计算。（4）消元计算：**第31页第31页

解：准确解为（舍入值）：例：用列主元消去法去解方程组计算解在常数项内得到（5）回代计算：第32页第32页

本例是含有舍入4位浮点数进行运算，所得计算解还是比较准确。下面是完全主元素消去法框图（图6-1）第33页第33页

图（6-1）stop输入n,A,b,cI(Z)I(I=1,2,…….n)k=1,2…….n是输出detA=0否互换行是否转下页第34页第34页

消元计算互换列且接上页否Stop第35页第35页

①②③第36页第36页

②—①③—①②*③*第二步，选4.12为列主元，不需换行，③*—②*③**由①②*③**回代求