现代数值计算方法(matlab版)第3章(4)课件

《现代数值计算方法(matlab版)第3章(4)课件》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、福建师范大学yêO{数计学院1/18JJII1nÙ)5§

2、{JIBackClose1nÙ)5§

3、{福建师范大学数计学院§3.5)顽§

4、Cholesky©){2/18顽§

5、3ó§O¥kX2 A^.§

6、XêÝA顽,3¢ÛÉenÝL¦TA=LL,(3.18)½Lé,ù«©)´.JJIIl11Jl21l22IL=.......,..Backln1ln2···lnnClosedA=LLT'AÚLLTéA,¦LlXeµijda=l2,a=ll,1111i1i111√福建师范大学l11=

7、a11,数计学院li1=ai1/l11,i=2,···,n.3/18bL1k−1®²¦,e¡¦L1klik,i=k,···,n.5¿Xk−1aik=lirlkr+liklkk,r=1kP−11/2JJl=a−l2,kkkkkrIIr=1kP−1(3.19)Jlik=aik−lirlkrlkk,i=k+1,···,n.Ir=1Backþã©){¡Cholesky©){.duOéIClosengm²$,Cholesky©){q¡²{.ØJy,Cholesky©){¦ØOoþn3/6+O(n2),ÝLU©)Oþ.,Xd

8、,ÙOnm$´~Ø

9、福建师范大学数计学院.4/18;m$,rÝA©)¤TA=LDL(3.20)Ù¥LüenÝ,D´éÝ,éþØ".

10、^(3.20)ü>éA{,ÑOLÚDOúª.bLÚD11k−1®²¦,±e¦§1kJJTIIlik,i=k+1,···,nÚdk,ùplkk=1.'A=LDLü>1Jk,Xk−1Iaik=lirdrlkr+likdk,i=k,···,n.Backr=1CloseddkP−1d=a−l2d,kkkkrrr=1(3.21)福建师范大学kP−1数计学院lik=aik−lirdr

11、lkrdk,i=k+1,···,n.r=15/18U(3.21)O,LDLT©);m$,dOúªØJuy,¦Ø$oþO,q¡En3/3+O(n2).~¦{$þ,Ú9ÏCþtik=likdk,¿òOúª(3.21)nXeµék=1,2,···,n,OJJkP−1tik=aik−tijlkj,i=k+1,···,n,IIj=1Jlik=tik/dk,i=k+1,···,n,(3.22)IkP−1dk=akk−tkjlkj.Backj=1CloseN´wÑ,U?LDLT©)¦Ø$þn3/6+O(n2),ØIm$.T{;Cþtik,;þAO

12、.e¡·ïá^Cholesky©){¦)顽§

13、{Ú福建师范大学数计学院½.Ly=b,6/18TAx=b⇒LDLx=b⇒Dz=y,(3.23)TLx=z.{3.4(Cholesky©){)Ú1Ñ顽ÝAÚmàþb¶Ú2Cholesky©)µJJIId1=t11=a11,li1=ai1/d1,i=2,···,n,Jék=2,···,nOµIkP−1dk=akk−tkjlkj,Backj=1ClosekP−1tik=aik−tijlkj,lik=tik/dk,i=k+1,···,n.j=1Ú3^c{)en§

14、Ly=bµ福建师范大学

15、数计学院y1=b1,kP−17/18ék=2,···,nOyk=bk−lkjyj;j=1Ú4)é/§

16、Dz=yµék=1,···,n,Oµzk=yk/dk¶Ú5^£{)þn§

17、LTx=zµxn=zn,JJPnék=n−1,···,1Oµxk=zk−ljkxj.IIj=k+1Jâ{3.4,?MATLAB§SXeµIBack•Cholesky©){MATLAB§SClose%machol.mfunction[x,l,d]=machol(A,b)%^å:^Cholesky©){)§

18、Ax=b福建师范大学数计学院%LDL’©)8/18n=length(b);d=

19、zeros(1,n);l=eye(n,n);d(1)=A(1,1);l(2:n,1)=A(2:n,1)/d(1);d(2)=A(2,2)-l(2,1)*l(2,1)*d(1);fori=3:nforj=2:(i-1)s=0;JJIIfork=1:(j-1)s=s+d(k)*l(i,k)*l(j,k);endJl(i,j)=(A(i,j)-s)/d(j);IendBackCloses=0;forj=1:(i-1)s=s+d(j)*l(i,j)*l(i,j);endd(i)=A(i,i)-s;福建师范大学数计学院end9/18%¦)en§

20、Ly=b(c{)y

21、=zeros(n,1);y(1)=b(1);fori