【圖文】第六章 常微分方程的數(shù)值解_第1頁
【圖文】第六章 常微分方程的數(shù)值解_第2頁
【圖文】第六章 常微分方程的數(shù)值解_第3頁
【圖文】第六章 常微分方程的數(shù)值解_第4頁
全文預(yù)覽已結(jié)束

下載本文檔

版權(quán)說明：本文檔由用戶提供并上傳，收益歸屬內(nèi)容提供方，若內(nèi)容存在侵權(quán)，請進(jìn)行舉報或認(rèn)領(lǐng)

文檔簡介

1、 例2 y c(1 y 2 y + y = 0 H=0,50 t = 0, y = 2, y = 1 （1）轉(zhuǎn)化方程。 （2）建立文件，設(shè)文件名為dzdt02.m。 function dz=f(t,z c=1; dz(1=z(2; 2 3 dz(2=c*(1-z(12*z(2-z(1; 1 dz=dz(1;dz(2; （3）調(diào)用求解。 H=0,50;z0=2;1; t,z=ode45(dzdt02,H,z0; plot(t,z(:,1,m-,t,z(:,2,g-. 0 -1 -2 -3 0 5 10 15 20 25 30 35 40 45 50 例3 x = y, y = 0.2 y sin

2、 t （1）轉(zhuǎn)化方程。 H=0,100 t = 0, x = 0, y = 2 （2）建立文件，設(shè)文件名為dzdt03.m。 function dz=f(t,z a=-0.2; 1.5 2 dz(1=z(2; 1 dz(2=a*z(2-sin(t; dz=dz(1;dz(2; （3）調(diào)用求解。 H=0,100;z0=0;2; 0.5 0 -0.5 t,z=ode45(dzdt03,H,z0; -1 plot(z(:,1,z(:,2,m- 0 1 2 3 4 5 6 x + 2ty = 0 例4 y 2tx = 0 H=-30，30 t = 0, x0 , y0 = 0,0, x0 , y0 =

3、 1,0 （1）轉(zhuǎn)化方程。 （2）建立文件，設(shè)文件名為dzdt04.m。 function dz=f(t,z dz(1=z(2; dz(2=-2*t*z(4; dz(3=z(4; dz(4=2*t*z(2; -0.5 0.5 0 dz=dz(1;dz(2;dz(3;dz(4; （3）調(diào)用求解。 H=-30,30;z0=0;1;0;0; t,z=ode45(dzdt04,H,z0; plot(z(:,1,z(:,3,m- -2 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 -1 -1.5 例5 x = ax + yz, y = b( y z , z = cy z xy （1）轉(zhuǎn)化方程。 （2）建立文件，設(shè)文件名為dzdt05.m。 function dz=f(t,z a=8/3;b=10;c=28; dz(1=-a*z(1+z(2*z(3; dz(2=-b*(z(2-z(3; 30 20 10 0 dz(3=c*z(2-z(3-z(1*z(2; dz=dz(1;dz(2;dz(3; （3）調(diào)用求解。 -10 -20 -30 20 10 60 0 -10 -20 -20 0 40 20 H=0,80;z0=0;0;0.2

溫馨提示

  • 1. 本站所有資源如無特殊說明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙，網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽，若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間，僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理，對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯，并不能對任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容，請與我們聯(lián)系，我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論