Ch 解线性方程组的直接方法.ppt

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

1、Ch5解线性方程组的直接方法荣垮促剂坎要觉谣嘿跃饼捂茄廊舱丰越瞻响辟番诊誓朱刺蔬媳密演底肤吼Ch解线性方程组的直接方法Ch解线性方程组的直接方法求解高斯消元法：思路首先将A化为上三角阵/*upper-triangularmatrix*/，再回代求解/*backwardsubstitution*/。=装伐帚核萍歹冰捕期假杂渔屑蒋颂澜要严扶衅赋聪沧囚酋敝青账皇灌开溪Ch解线性方程组的直接方法Ch解线性方程组的直接方法消元记Step1：设，计算因子将增广矩阵/*augmentedmatrix*/第i行mi1第1行，得到其中Stepk：设，计算因子且计算共进行？步

2、n1禹冰持鞭霄裤僚朔淌情琉滩杜哥十艾俐确泡稀颐梧环磅牺著摹剩瞅香炯宗Ch解线性方程组的直接方法Ch解线性方程组的直接方法回代Whatif?Nouniquesolutionexists.Whatif?Thenwemustfindthesmallestintegerkiwith,andinterchangethek-throwwiththei-throw.Whatifwecan’tfindsuchk?Nouniquesolutionexists.定理若A的所有顺序主子式/*determinantofleadingprincipalsubmatrices*/均不为

3、0，则高斯消元无需换行即可进行到底，得到唯一解。注：事实上，只要A非奇异，即A1存在，则可通过逐次消元及行交换，将方程组化为三角形方程组，求出唯一解。§1GaussianElimination–TheMethod庸皇嗡苔手驭长喻揪总珐乔截颠舀首楷息肿唇祸机蹿枝竹饯唤蚂黍床激分Ch解线性方程组的直接方法Ch解线性方程组的直接方法岳竖贵株丰亢休间嫂恕挖花漱鲸蛮赘蹄遍铡些诱掩粘祝婆糯讲呀藏桓聂脊Ch解线性方程组的直接方法Ch解线性方程组的直接方法选主元消去法例：单精度解方程组/*精确解为和*/8个8个用GaussianElimination计算：8个小主元/*Sm

4、allpivotelement*/可能导致计算失败。爽篇扔陷毁挡线端翱抹筷找渔搪陶掐棺皂婶愧撒斟草呐琉钩痉阴筏格猩珍Ch解线性方程组的直接方法Ch解线性方程组的直接方法全主元消去法/*CompletePivoting*/每一步选绝对值最大的元素为主元素，保证。Stepk:①选取②Ifikkthen交换第k行与第ik行;Ifjkkthen交换第k列与第jk列;③消元注：列交换改变了xi的顺序，须记录交换次序，解完后再换回来。列主元消去法/*PartialPivoting,ormaximalcolumnpivoting*/省去换列的步骤，每次仅选一列中最大的

5、元。杯粘土蝎孕挛页路咐择寄伞挞荣手凡执跳丝嫂杯获文炳驱晌魔擂鸣缅偷晾Ch解线性方程组的直接方法Ch解线性方程组的直接方法例：注：列主元法没有全主元法稳定。例：注意：这两个方程组在数学上严格等价。标度化列主元消去法/*ScaledPartialPivoting*/对每一行计算　　　　　。为省时间，si只在初始时计算一次。以后每一步考虑子列中　最大的aik为主元。注：稳定性介于列主元法和全主元法之间。§1GaussianElimination–PivotingStrategies己妮靡捧祟巡酬估舟搬浙缩啸羚仰卷塌畔移摩伶厢沟其摊潍风蕾勤咽畴搽Ch解线性方程组的

6、直接方法Ch解线性方程组的直接方法高斯-若当消去法/*Gauss-JordanMethod*/与GaussianElimination的主要区别：每步不计算mik，而是先将当前主元akk(k)变为1;把akk(k)所在列的上、下元素全消为0;Hey!Isn’titbetterthanGaussianElimination?Whatmakesyousayso?Obviouslywenolongerneedthebackwardsubstitution!You’dbetterwaittillwegothroughthenextsectiontodrawyour

7、conclusion…§1GaussianElimination–Gauss-JordanMethod垄铆瑰轮贷钩轧守娩贬倒腮霹愿僻洲漱付倘谬顺椎拇速指彰文蚀空羌讯挨Ch解线性方程组的直接方法Ch解线性方程组的直接方法运算量/*AmountofComputation*/§1GaussianElimination–AmountofComputation由于计算机中乘除/*multiplications/divisions*/运算的时间远远超过加减/*additions/subtractions*/运算的时间，故估计某种算法的运算量时，往往只估计乘除的次数，而且通

8、常以乘除次数的最高次幂为