matlab在常微分方程的应用,MATLAB在常微分方程上应用.ppt
您所在位置:網(wǎng)站首頁(yè) > 海量文檔
 > 計(jì)算機(jī) > matlab
MATLAB在常微分方程上應(yīng)用.ppt11頁(yè)
本文檔一共被下載:次,您可全文免費(fèi)在線閱讀后下載本文檔。
下載提示
1.本站不保證該用戶上傳的文檔完整性,不預(yù)覽、不比對(duì)內(nèi)容而直接下載產(chǎn)生的反悔問(wèn)題本站不予受理。
2.該文檔所得收入(下載+內(nèi)容+預(yù)覽三)歸上傳者、原創(chuàng)者。
3.登錄后可充值,立即自動(dòng)返金幣,充值渠道很便利
* * MATLAB在常微分方程上的應(yīng)用 2005/8/30 楊惠如 老師:王天楷教授 數(shù)值分析法 Taylor series Runge-Kutta Method Taylor’s series Runge-Kutta Method 一階 二階 A first-order Runge-Kutta method: uses a first-order Taylor series expansion y後=y前+(x後-x前)y’前+ y’’前 一階 二階 y後=y前+(x後-x前)y’前 MATLAB在常微分方程上的應(yīng)用 函數(shù)解法的特點(diǎn) M檔案 ODE範(fàn)例 函數(shù)解法的特點(diǎn) 採(cǎi)用TR-BDF2方法,和ode23s一樣,對(duì)於誤差容許範(fàn)圍較寬的情況,比ode15s效果好,可解剛性問(wèn)題 ode23tb 可變階的數(shù)值微分公式法(NDFs),屬於多步解法,如果用在ode45效果很差或失敗,可考慮用此方式。可解剛性問(wèn)題。 ode15s 採(cǎi)用改進(jìn)的2階Rosenbrock法,採(cǎi)單步法,比ode15s更適於誤差範(fàn)圍容許較寬的情況。 ode23s 採(cǎi)自由內(nèi)插法實(shí)現(xiàn)的梯形規(guī)則,適用剛性並要求無(wú)數(shù)值衰減的問(wèn)題 ode23t 採(cǎi)可變階的Adams-Bashforth-Moulton PECE法,屬於多步解法,比ode45更適合於誤差容許範(fàn)圍要求較嚴(yán)格時(shí) ode113 2階/3階Rung-Kutta(2,3)法。對(duì)於誤差容許範(fàn)圍較寬且稍帶剛性問(wèn)題比ode45效果好 ode23 4階/5階Rung-Kutta(4,5)法。採(cǎi)單步法,速度快,作為解問(wèn)題的初試法 ode45 資料來(lái)源:蔡宏睿 M檔案 以MATLAB程式碼所撰寫(xiě)的檔案通常以 “ m ” 為副檔名,所以這些檔案又稱為M檔案 M檔案又可分為兩類(lèi)—”函式” 和”主程式” 函式: 也是M檔案的一種,可以自行定義一個(gè)函式,定義好之後經(jīng)過(guò)存檔,不必執(zhí)行,當(dāng)對(duì)此函式輸入變數(shù)時(shí),程式會(huì)自動(dòng)讀取此已定義的函式 主程式: 是一個(gè)副檔名為m的檔案,其內(nèi)容包含MATLAB的各種指令,我們只要輸入主檔名,即可逐一執(zhí)行在此檔案內(nèi)的所有指令 資料來(lái)源:蔡宏睿 function dy=ff(x,y) dy=3*x.^2 ode45(‘ff’,[2 7],0.5) 函式 執(zhí)行 y’=3x2 檔名 X的範(fàn)圍 y的初值 ODE範(fàn)例 ODE範(fàn)例 函式 function yprime=cc(x,y) yprime=[y(2);-sin(y(1))] 執(zhí)行 ode45(‘cc’,[0 10],[1;1.5]) y(1)=1 y(2)=3
發(fā)表評(píng)論
請(qǐng)自覺(jué)遵守互聯(lián)網(wǎng)相關(guān)的政策法規(guī),嚴(yán)禁發(fā)布色情、暴力、反動(dòng)的言論。
用戶名:
驗(yàn)證碼:
匿名?
發(fā)表評(píng)論
總結(jié)
以上是生活随笔為你收集整理的matlab在常微分方程的应用,MATLAB在常微分方程上应用.ppt的全部?jī)?nèi)容,希望文章能夠幫你解決所遇到的問(wèn)題。
- 上一篇: 信息学奥赛一本通 2040:【例5.7】
- 下一篇: matlab1分钟数据转5分钟,把股票的