现代数值计算方法(matlab版)第3章(2)

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

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

2、{JIBackClose1nÙ)5§

3、{福建师范大学数计学院§3.2ÌGauss{9§S¢y2/16(k)^SGauss{OL§´Ø,Ñya=0,Okk(k)ÒÃ{?1e.=¦é¤kk=1,2,···,n,a6=0,ØUykkOL§´ê½.~3.4k5§

4、JJ0.0001x1+1.0x2=1.0,II1.0x1+1.0x2=2.0.JÙ°()I100009998Backx1=≈1.00010,x2=≈0.9999099999999Closey3b½^ê4?i2:ê5¦).)L§µâ42:ê$5K

5、1.0−10000.0=(0.00001−福建师范大学550.1)10=(0.0000−0.1)10=−10000.0(),Ón,2.0−10000.0=数计学院−10000.0,3/160.00011.01.0r2−104r10.00011.01.0−−−−−→1.01.02.001.0−10000.02.0−10000.0.00011.01.0−−−→.0−10000.0−10000.0JJ£L§µIIJ0.0001x1+1.0x2=1.0,x2=1.0,I=⇒−10000.0x2=−10000.0.x1=0.0.BackClose§

6、

7、,uy(Jîý.©Û(JýÏuy,du1Ì0.0001ýéLu,§3L§¥©1r¥mL§êâ10000,¦福建师范大学数计学院¥m(J/¯0K©êâ,lE¤êؽ.4/16é±þ¯K,ÄÀ^ýéêÌ.L§µ0.00011.01.0r1↔r21.01.02.0−−−→1.01.02.00.00011.01.0r2−0.0001r11.01.02.0JJ−−−−−−→II01.0−0.00011.0−0.0002J1.01.02.0I−−−→.Back01.01.0Closeùp,L§1.0−0.0001=(0.1−0.0000

8、1)101(),Ón1.0−0.0002=1.0.£L§µ福建师范大学数计学院1.0x1+1.0x2=2.0,x2=1.0,5/16=⇒1.0x2=1.0.x1=1.0.§

9、,uy(JÄÜn.þã~f`²ÀÌ5.e¡ãÌGauss{Äg.PA(1)=A,3GaussL§11Ú,11(1)JJ¥ýéa,=r11II(1)(1)a=max

10、a

11、r11i1J1≤i≤nIÌ.er1>1,1r11Ú111.BackClose/,3GaussL§1kÚ,(k)(k)a=max

12、a

13、(3.6)rkkikk≤i≤n福建师范大学数计学院Ì.er

14、k>k,1rk1Ú1k1.6/16ÌGauss{{Ú½Xeµ{3.2(ÌGauss{)Ú1ÑXêÝA,màb,k:=1¶Ú2ék=1,···,n−1?1Xeöµ(1)ÀÌ,(½rk,¦JJ(k)(k)a=max

15、a

16、,rkkikIIk≤i≤nea(k)=0,KÊO,ÄK,?1eÚ¶Jrkk(2)(k)(k)Ierk>k,(A,b)1k,rkü1¶Back(3)µéi,j=k+1,···,n,OClose(k)(k)(k+1)mik=aik/akk,aik=0,(k+1)(k)(k)(k+1)(k)(k)aij=aij−mikakj,bi=bi−mikbk.福建师范大

17、学Ú3£数计学院(n)(n)xn=bn/ann,7/16ék=n−1,···,1,OPn(k)(k)(k)xk=bk−akjxjakk.j=k+1â{3.2,ÑMATLAB§SXeµ•ÌGauss{MATLAB§SJJII%magauss2.mJfunctionx=magauss2(A,b,flag)I%^å:ÌGauss{)5§

18、Ax=bBackClose%ª:x=magauss(A,b,flag),AXêÝ,bmà,%eflag=0,KØw«¥mL§,ÄKw«¥mL§,%%@0,x)þ福建师范大学数计学院ifnargin<3,flag=0;end8/16n=length

19、(b);fork=1:(n-1)%ÀÌ[ap,p]=max(abs(A(k:n,k)));p=p+k-1;JJifp>kIIt=A(k,:);A(k,:)=A(p,:);A(p,:)=t;Jt=b(k);b(k)=b(p);b(p)=t;IendBackClose%m=A(k+1:n,k)/A(k,k);A(k+1:n,k+1:n)=A(k+1:n,k+1:n)-m*A(k,k+1:n);福建师范大学数计