版權說明:本文檔由用戶提供并上傳,收益歸屬內容提供方,若內容存在侵權,請進行舉報或認領
文檔簡介
1、第六章 常微分方程的數(shù)值解常見的近似數(shù)值求解方法有歐拉折線法、阿當姆斯法、龍格-庫塔法與吉爾法。不受方程類型的限制,可以求任何形狀常微分方程的特解,但求出的解只能是數(shù)值的解函數(shù)。一、解法步驟: 將高階方程轉化為一階方程組。 (兩種情況) 建立相應的函數(shù)文件。 調用求解函數(shù)。 t,z=odeij(dzdtk,H,z0,tol) 剛性:設有一階常系數(shù)線性微分方程組y=Ay+f,如果它的Jacobian矩陣的特征值相差十分懸殊。odeij問題類型精度適用對象ode45非剛性中等多數(shù)情況下可優(yōu)先選用,但不能用來解剛性問題ode23非剛性較低可用來解中等剛性問題,或誤差允許范圍較寬的問題ode113非剛
2、性低到高不能解剛性問題,當誤差容限要求嚴格時效果較ode45好ode15s剛性低到高可用于解剛性問題,當采用ode45失敗或效果很差時,可考慮使用ode23s剛性低可用于解剛性問題,當誤差容限較寬時效果比ode15s好ode23t適度剛性低可用于解剛性問題,但要求無數(shù)值衰減ode23tb剛性低可用于解剛性問題,當誤差容限較寬時效果比ode15s好二、例題例1 在區(qū)間H=0.1,30上滿足初值條件x=0.1時, 的特解(1轉化原方程為一階方程組。(2建立描述上面微分方程組的函數(shù)文件,并設文件名為dzdt01.m。function dz=f(t,z)dz(1)=z(2);dz(2)=z(3);dz
3、(3)=(1/t)*z(3)-(2/t2)*z(2)+(2/t3)*z(1)+10*t3*sin(t);dz=dz(1);dz(2);dz(3)xxyyxyxyxsin1022623 0, 0, 0 yyy(3調用求解。H=0.1,30;z0=0;0;0;t,z=ode45(dzdt01,H,z0);plot(t,z(:,1),m-)plot(t,z(:,1),m-)051015202530-2.5-2-1.5-1-0.500.511.5x 105例2 H=0,50(1轉化方程。(2建立文件,設文件名為dzdt02.m。function dz=f(t,z)c=1;dz(1)=z(2);dz(2
4、)=c*(1-z(1)2)*z(2)-z(1);dz=dz(1);dz(2);(3調用求解。H=0,50;z0=2;1;t,z=ode45(dzdt02,H,z0);plot(t,z(:,1),m-,t,z(:,2),g-.)05101520253035404550-3-2-101230)1 (2 yyycy1, 2, 0yyt例3 H=0,100(1轉化方程。(2建立文件,設文件名為dzdt03.m。function dz=f(t,z)a=-0.2;dz(1)=z(2);dz(2)=a*z(2)-sin(t);dz=dz(1);dz(2);(3調用求解。H=0,100;z0=0;2;t,z=
5、ode45(dzdt03,H,z0);plot(z(:,1),z(:,2),m-)0123456-1-0.500.511.52tyyyxsin2 . 0,2, 0, 0yxt例4 H=-30,30(1轉化方程。(2建立文件,設文件名為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);dz=dz(1);dz(2);dz(3);dz(4);(3調用求解。H=-30,30;z0=0;1;0;0;t,z=ode45(dzdt04,H,z0);plot(z(:,1),z(:,3),m-)-0.4-0.200.20.40.60.811.21.41.6-2-1.5-1-0.500.5 0202x tyy tx0 , 1 ,0 , 0, , 00000yxyxt例5 (1轉化方程。(2建立文件,設文件名為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);dz(3)=c*z(2)-z(3)-z(1)*z(2);dz=dz(1);dz(2);dz(3);(3調用求解。H=0,80;z0=0;0;0.222/1015;t,z=ode23(dzdt05
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網頁內容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
- 4. 未經權益所有人同意不得將文件中的內容挪作商業(yè)或盈利用途。
- 5. 人人文庫網僅提供信息存儲空間,僅對用戶上傳內容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內容本身不做任何修改或編輯,并不能對任何下載內容負責。
- 6. 下載文件中如有侵權或不適當內容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 房產繼承和贈與合同(2篇)
- 2024年煙塵、粉塵自動采樣器及測定儀項目建議書
- 2024年抗滴蟲病藥合作協(xié)議書
- 2024年專用系統(tǒng)集成電路項目建議書
- 2024年嘧菌酯合作協(xié)議書
- 2024年新型電子封裝材料項目合作計劃書
- 中小學-陳述句與反問句互改-課件
- 食品行業(yè)運輸外包合同樣本
- 2024版服裝訂購合同范本
- 軟裝設計合同模板
- 中國心血管病風險評估和指南
- 心臟外科手術簡介
- 先進核電技術研發(fā)項目實施方案
- 結直腸癌中西醫(yī)結合治療護理干預的效果評估
- 雙方合作合同模板
- (完整版)二年級乘加乘減口算100題
- 公路工程質量檢驗評定標準 第二冊機電工程 JTG 2182-2020
- 小學各年級做好時間管理安排充實的生活主題班會
- 兒童商場游樂設置方案
- 如何做好研究生導師參考幻燈片
- 小學各年級“紅領巾獎章”獲章標準
評論
0/150
提交評論