




下載本文檔
版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
1、9.2對(duì)初值問題試用歐拉方法取步長(zhǎng)h=0.1和h=0.2計(jì)算其近似解,并與精確解比較。解: 歐拉方法:h=0.1;a=0;b=1;n=(b-a)/h;1y=zeros(1,n+1);xx=0;yy=0;y(1)=yy+h*(1.0/(1+xx.2)-2*(yy.2);xx=xx+h;yy=y(1);for i=1:ny(i+1)=y(i)+h*(1.0/(1+xx.2)-2*(yy.2);xx=xx+h;yy=y(i+1);end; yy = 0.1000 0.1970 0.2854 0.3609 0.4210 0.4656 0.4957 0.5137 0.5219 0.5227 0.5180
2、當(dāng)h=0.2時(shí)y = 0.2000 0.3763 0.4921 0.5423 0.5466 0.5271精確計(jì)算的結(jié)果:xx=0;h=0.1;a=0;b=1;n=(b-a)/h;y=zeros(1,n);for i=1:nxx=xx+h;y(i)=xx/(1+xx.2);end; yy =0.0990 0.1923 0.2752 0.3448 0.4000 0.4412 0.4698 0.4878 0.4972 0.5000誤差的發(fā)展情況: h=0.1;jsy= 0.1970 ;0.2854 ; 0.3609; 0.4210 ; 0.4656 ; 0.4957 ;0.5137 ; 0.5219
3、; 0.5227;0.5180;jqy=0.0990;0.1923;0.2752 ; 0.3448 ; 0.4000 ; 0.4412 ; 0.4698; 0.4878;0.4972 ; 0.5000; x=zeros(10,1); xx=0; x(1)=xx+h;for i=2:10x(i)=x(i-1)+h; end;yc=zeros(10,1);for j=1:10yc(j)=jsy(j)-jqy(j);end; plot(x,yc,b)誤差yc= 0.0980 0.0931 0.0857 0.0762 0.0656 0.0545 0.0439 0.0341 0.0255 0.0180誤
4、差的走勢(shì)圖為:9.3利用歐拉預(yù)測(cè)-校正法和經(jīng)典龍格庫(kù)塔方法,取步長(zhǎng)為h=0.1,求解方程并與精確解進(jìn)行比較。解: 歐拉預(yù)測(cè)-矯正方法h=0.1;a=0;b=1;n=(b-a)/h;yy=1;x=zeros(1,n+1);y=zeros(1,n);for i=2:11x(i)=x(i-1)+h;end;yb=yy+h*(x(1)+yy);y(1)=yy+(h/2)*(x(1)+yy+x(2)+yb);for j=2:nyb=y(j-1)+h*(x(j)+y(j-1);y(j)=y(j-1)+(h/2)*(x(j)+y(j-1)+x(j+1)+yb);end; yy =1.1100 1.2421
5、1.3985 1.5818 1.7949 2.0409 2.3231 2.6456 3.0124 3.4282四階經(jīng)典龍格庫(kù)塔方法:h=0.1;y0=1;a=0;b=1;n=(b-a)/h;y=zeros(1,n);x=zeros(1:10);k1=h*(0+1); k2=h*(0+h/2+1+k1/2); k3=h*(0+h/2+1+k2/2); k4=h*(0+h+1+k3);y(1)=y0+(k1+2*k2+2*k3+k4)/6; x(1)=0.1;for i=2:10x(i)=x(i-1)+h;end;for j=2:10k1=h*(x(j-1)+y(j-1); k2=h*(x(j-1
6、)+h/2+y(j-1)+k1/2);k3=h*(x(j-1)+h/2+y(j-1)+k2/2); k4=h*(x(j-1)+h+y(j-1)+k3);y(j)=y(j-1)+(k1+2*k2+2*k3+k4)/6;end; yy =1.1103 1.2428 1.3997 1.5836 1.7974 2.0442 2.3275 2.6511 3.0192 3.4366精確計(jì)算: x=0;h=0.1;a=0;b=1;n=(b-a)/h;y=zeros(1,n);for i=1:nx=x+h;y(i)=-x-1+2*exp(x);end; yy =1.1103 1.2428 1.3997 1.5
7、836 1.7974 2.0442 2.3275 2.6511 3.0192 3.4366誤差發(fā)展:h=0.1;jsy= 1.1100 1.2421 1.3985 1.5818 1.7949 2.0409 2.3231 2.6456 3.0124 3.4282;jqy=1.1103 1.2428 1.3997 1.5836 1.7974 2.0442 2.3275 2.6511 3.0192 3.4366; x=zeros(10,1); xx=0; x(1)=xx+h;for i=2:10x(i)=x(i-1)+h; end;yc=zeros(10,1);for j=1:10yc(j)=abs
8、(jsy(j)-jqy(j);end; ycplot(x,yc,b)ans =0.0003 0.0007 0.0012 0.0018 0.0025 0.0033 0.0044 0.0055 0.0068 0.0084誤差的走勢(shì)圖:9.5用Adams預(yù)測(cè)校正法求解并與精確解進(jìn)行比較。解:當(dāng)取步長(zhǎng)為0.01時(shí):h=0.01;y0=1;a=0;b=1;n=(b-a)/h;y=zeros(1,n);x=zeros(1,n);k1=h*(-1); k2=h*(-(1+k1/2).2); k3=h*(-(1+k2/2).2); k4=h*(-(1+k3).2);y(1)=y0+(k1+2*k2+2*k3+
9、k4)/6; x(1)=0.01;for i=2:nx(i)=x(i-1)+h;end;for j=2:nk1=h*(-y(j-1)*y(i-1); k2=h*(-(y(j-1)+k1/2).2);k3=h*(-(y(j-1)+k2/2).2); k4=h*(-(y(j-1)+k3).2);y(j)=y(j-1)+(k1+2*k2+2*k3+k4)/6;end;yb=zeros(1,100);yb(4)=y(3)+(h/24)*(55*(-y(3)*y(3)-59*(-y(3-1).2)+37*(-y(3-2).2)-9*(-(y0).2);y(4)=y(3)+(h/24)*(9*(-yb(4
10、).2)+19*(-y(3).2)-5*(-y(3-1).2)+(-y(3-2).2);for k=4:99yb(k+1)=y(k)+(h/24)*(55*(-y(k)*y(k)-59*(-y(k-1).2)+37*(-y(k-2).2)-9*(-y(k-3).2);y(k+1)=y(k)+(h/24)*(9*(-yb(k+1).2)+19*(-y(k).2)-5*(-y(k-1).2)+(-y(k-2).2);end;y(:,1:10) %100個(gè)校正值的前十個(gè)輸出;y(:,91:100) %100個(gè)校正值的最后是個(gè)輸出;ans = 0.9901 0.9820 0.9740 0.9646 0
11、.9554 0.9464 0.9375 0.9288 0.9202 0.9118ans =0.5245 0.5217 0.5190 0.5163 0.5137 0.5111 0.5085 0.5059 0.5034 0.5008精確解:h=0.01;y0=1;a=0;b=1;n=(b-a)/h;y=zeros(1,n);x=zeros(1,n);x(1)=0.01;for i=2:nx(i)=x(i-1)+h;end; h=0.01;y0=1;a=0;b=1;n=(b-a)/h;y=zeros(1,n);x=zeros(1,n);x(1)=0.01;for i=2:nx(i)=x(i-1)+h;end;for j=1:ny(j)=1/(x(j)+1);end;y(:,1:10) %100個(gè)校正值的前十個(gè)輸出;y(:,91:1
溫馨提示
- 1. 本站所有資源如無(wú)特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 員工三級(jí)安全教育培訓(xùn)
- 保險(xiǎn)綜合內(nèi)勤培訓(xùn)
- 鄉(xiāng)土樹種繁育基地行業(yè)深度調(diào)研及發(fā)展項(xiàng)目商業(yè)計(jì)劃書
- 風(fēng)箏沖浪教學(xué)與體驗(yàn)企業(yè)制定與實(shí)施新質(zhì)生產(chǎn)力項(xiàng)目商業(yè)計(jì)劃書
- 新型冠狀病毒感染的肺炎流行醫(yī)療救治實(shí)施方案
- 籃球興趣班課程教學(xué)計(jì)劃
- 初中學(xué)校教研教材本土化計(jì)劃
- 2024-2025日語(yǔ)基礎(chǔ)班教學(xué)計(jì)劃
- 部編四年級(jí)上冊(cè)《道德與法治》題庫(kù)建設(shè)計(jì)劃
- 管理咨詢師個(gè)人成長(zhǎng)計(jì)劃
- GB/T 39866-2021建筑門窗附框技術(shù)要求
- GB/T 3836.2-2021爆炸性環(huán)境第2部分:由隔爆外殼“d”保護(hù)的設(shè)備
- 關(guān)節(jié)脫位患者的護(hù)理-關(guān)節(jié)脫位患者的護(hù)理(外科護(hù)理ppt)
- 產(chǎn)品合格證模板-合格證模板樣本
- 領(lǐng)導(dǎo)干部的決策力與執(zhí)行力
- 水泵試運(yùn)行調(diào)試記錄
- 半導(dǎo)體中載流子的統(tǒng)計(jì)分布和計(jì)算
- 組織部處級(jí)干部培訓(xùn)審計(jì)財(cái)經(jīng)紀(jì)律課件
- 史上最全最權(quán)威婦產(chǎn)科icd編碼培訓(xùn)【版】課件
- 心血管診治與搶救標(biāo)準(zhǔn)操作規(guī)程(SOP)
- 中南大學(xué)自主招生綜合素質(zhì)測(cè)試面試試題答題技巧
評(píng)論
0/150
提交評(píng)論