版權(quán)說(shuō)明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
1、龍格-庫(kù)塔方法是一種經(jīng)典方法,具有很高的精度,它間接的利用了泰勒級(jí)數(shù)展開(kāi),避免了高階偏導(dǎo)數(shù)的計(jì)算。此處以最為經(jīng)典的四級(jí)四階龍格-庫(kù)塔方法為例,計(jì)算格式如下 1龍格-庫(kù)塔法解一階ODE對(duì)于形如 的一階ODE初值問(wèn)題,可以直接套用公式,如今可以借助計(jì)算機(jī)方便的進(jìn)行計(jì)算,下面給出一個(gè)實(shí)例 取步長(zhǎng)h=0.1,此處由數(shù)學(xué)知識(shí)可得該方程的精確解為 。在這里利用MATLAB編程,計(jì)算數(shù)值解并與精確解相比,代碼如下:(1)寫出微分方程,便于調(diào)用和修改function val = odefun( x,y ) val = y-2*x/y; end(2)編寫runge-kutta方法的函數(shù)代碼function y
2、= runge_kutta( h,x0,y0 )k1 = odefun(x0,y0);k2 = odefun(x0+h/2,y0+h/2*k1);k3 = odefun(x0+h/2,y0+h/2*k2);k4 = odefun(x0+h,y0+h*k3); y = y0+h*(k1+2*k2+2*k3+k4)/6; end(3)編寫主函數(shù)解微分方程,并觀察數(shù)值解與精確解的差異clear allh = 0.1;x0 = 0;y0 = 1;x = 0.1:h:1;y(1) = runge_kutta(h,x0,y0);for k=1:length(x) x(k) = x0+k*h; y(k+1)
3、 = runge_kutta(h,x(k),y(k);endz = sqrt(1+2*x);plot(x,y,*);hold onplot(x,z,'r');結(jié)果如下圖,數(shù)值解與解析解高度一致2龍格-庫(kù)塔法解高階ODE對(duì)于高階ODE來(lái)說(shuō),通用的方法是將高階方程通過(guò)引入新的變量降階為一階方程組,此處仍以一個(gè)實(shí)例進(jìn)行說(shuō)明。 這是一個(gè)二階ODE,描述的是一個(gè)物體的有阻尼振動(dòng)情況。初始條件為 ,將方程降階,引入一個(gè)向量型變量Y故有 記則至此,二階方程降階為一階方程組。值得注意的是此時(shí)再用龍格-庫(kù)塔法進(jìn)行求解時(shí),代入的將是一個(gè)Y向量。同樣利用MATLAB進(jìn)行計(jì)算,步長(zhǎng)h=0.05,時(shí)間周
4、期為0,20.(1) 編寫ODE函數(shù)function Y = odefun1( ,Y0 )% 此處Y0為一個(gè)列向量,因?yàn)闀r(shí)間t未顯含在一階方程組中% 所以ode函數(shù)的第一個(gè)參數(shù)為空,要根據(jù)具體情況而定。Y = Y0(2); (2000-200*Y0(2)-750*Y0(1)/500;end(2) 編寫runge-kutta函數(shù)function Y = rkfa( h,t0,Y0 )k1 = odefun1(t0,Y0);k2 = odefun1(t0+h/2,Y0+h/2*k1);k3 = odefun1(t0+h/2,Y0+h/2*k2);k4 = odefun1(t0+h,Y0+h*k3);Y = Y0+h*(k1+2*k2+2*k3+k4)/6;end(3) 編寫主函數(shù)clear allh = 0.05;t = 0.05:h:20;t0 = 0;Y0 = 0; 0;%初值Y = cell(1,length(t);Y1 = rkfa( h,t0,Y0 );z = zeros(2,length(t);for k=1:length(t) Yk+1 = rkfa( h,t0,Yk); z(1,k) = Yk(1);z(2,k) = Yk(2
溫馨提示
- 1. 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝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ù)覽,若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
- 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ì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 2023年攀枝花市仁和區(qū)事業(yè)單位引才筆試真題
- 2025屆遼寧省大連高三適應(yīng)性調(diào)研考試物理試題含解析
- 廣東省2024年教師資格證 單項(xiàng)選擇題+材料分析題+寫作題真題沖刺卷上半年B卷
- 北京市2023年教師資格證 單項(xiàng)選擇題+材料分析題+寫作題真題沖刺卷9月份A卷
- 小學(xué)班主任工作經(jīng)驗(yàn)交流(10篇)
- 小貓咪咪作文4篇
- 采礦課程設(shè)計(jì)
- 軟件開(kāi)發(fā)課程設(shè)計(jì)方向
- 2024-2025學(xué)年云南省富寧縣數(shù)學(xué)九上開(kāi)學(xué)檢測(cè)試題【含答案】
- 2024年稀土鋁合金、銅合金材料項(xiàng)目發(fā)展計(jì)劃
- 牛津樹單詞卡(可放大彩打)
- 硬筆字常用筆畫練習(xí)綠色米字格
- (完整版)水泥土攪拌樁規(guī)范
- 兒童急救小常識(shí)(課堂PPT)
- 口腔診所器材清單
- 碳纖維粘貼加固施工方案匯總
- 坪塘小學(xué)2020-2021學(xué)年度第二學(xué)期教學(xué)工作總結(jié)
- 東林論文要求
- 中醫(yī)診所規(guī)章制度61728
- 肩難產(chǎn)的處理ppt課件
- 2021年保育員培訓(xùn)記錄
評(píng)論
0/150
提交評(píng)論