




版權說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權,請進行舉報或認領
文檔簡介
1、數(shù)學建模與數(shù)學實驗數(shù)學建模與數(shù)學實驗 微微 分分 方方 程程微分方程的解析解微分方程的解析解 求微分方程(組)的解析解命令:dsolve(方程方程1, 方程方程2,方程方程n, 初始條件初始條件, 自變量自變量) 記號: 在表達微分方程時,用字母 D 表示求微分,D2、D3 等表示求高階微分.任何 D 后所跟的字母為因變量,自變量可以指定或由系統(tǒng)規(guī)則選定為確省.例如,微分方程 022dxyd應表達為:D2y=0.例例 1 求 21 udtdu 的通解.解解 輸入命令:dsolve(Du=1+u2,t)To Matlab(ff1) 結(jié) 果:u = tan(t-c)例例 2 求微分方程的特解. 1
2、5)0( , 0)0(029422yyydxdydxyd 解解 輸入命令: y=dsolve(D2y+4*Dy+29*y=0,y(0)=0,Dy(0)=15,x)結(jié) 果 為 : y =3e-2xsin(5x)To Matlab(ff2) 例例 3 求微分方程組的通解. zyxdtdzzyxdtdyzyxdtdx244354332解解 輸入命令 : x,y,z=dsolve(Dx=2*x-3*y+3*z,Dy=4*x-5*y+3*z,Dz=4*x-4*y+2*z, t); x=simple(x) % 將x化簡 y=simple(y) z=simple(z)結(jié) 果 為:x = (c1-c2+c3+
3、c2e -3t-c3e-3t)e2t y = -c1e-4t+c2e-4t+c2e-3t-c3e-3t+c1-c2+c3)e2t z = (-c1e-4t+c2e-4t+c1-c2+c3)e2t To Matlab(ff3)返 回微分方程的數(shù)值解微分方程的數(shù)值解(一)常微分方程數(shù)值解的定義(一)常微分方程數(shù)值解的定義 在生產(chǎn)和科研中所處理的微分方程往往很復雜且大多得不出一般解。而在實際上對初值問題,一般是要求得到解在若干個點上滿足規(guī)定精確度的近似值,或者得到一個滿足精確度要求的便于計算的表達式。因此,研究常微分方程的數(shù)值解法是十分必要的因此,研究常微分方程的數(shù)值解法是十分必要的。的相應近似值求
4、出準確值,值處,即對的若干離散的開始其數(shù)值解是指由初始點,:對常微分方程nnnyyxyxyxxxxxy,y )(,),(),y(x x )y(xy)f(x,y 2121210000返 回(二)建立數(shù)值解法的一些途徑(二)建立數(shù)值解法的一些途徑001i)y(xy)f(x,y , 1, 2 , 1 , 0 , xynihxi解微分方程:可用以下離散化方法求設1、用差商代替導數(shù)、用差商代替導數(shù) 若步長h較小,則有hxyhxyxy)()()( 故有公式:1-n,0,1,2,i )(),(001xyyyxhfyyiiii此即歐拉法歐拉法。2、使用數(shù)值積分、使用數(shù)值積分對方程y=f(x,y), 兩邊由xi
5、到xi+1積分,并利用梯形公式,有:)(,()(,(2)(,()()(11111iiiiiixxiixyxfxyxfxxdttytfxyxyii實際應用時,與歐拉公式結(jié)合使用:, 2 , 1 , 0 ),(),(2),()(11)1(1)0(1kyxfyxfhyyyxhfyykiiiiikiiiii的計算。然后繼續(xù)下一步,取時,當滿足,對于已給的精確度)( y y 2i111i)(1)1(1kikikiyyy此即改進的歐拉法改進的歐拉法。故有公式:)(),(),(200111xyyyxfyxfhyyiiiiii3、使用泰勒公式、使用泰勒公式 以此方法為基礎,有龍格龍格-庫塔法庫塔法、線性多步法
6、線性多步法等方法。4、數(shù)值公式的精度、數(shù)值公式的精度 當一個數(shù)值公式的截斷誤差可表示為O(hk+1)時(k為正整數(shù),h為步長),稱它是一個k階公式階公式。k越大,則數(shù)值公式的精度越高。歐拉法是一階公式,改進的歐拉法是二階公式。龍格-庫塔法有二階公式和四階公式。線性多步法有四階阿達姆斯外插公式和內(nèi)插公式。返 回(三)用(三)用Matlab軟件求常微分方程的數(shù)值解軟件求常微分方程的數(shù)值解t,x=solver(f,ts,x0,options)ode45 ode23 ode113ode15sode23s由待解方程寫成的m-文件名ts=t0,tf,t0、tf為自變量的初值和終值函數(shù)的初值ode23:組合
7、的2/3階龍格-庫塔-芬爾格算法ode45:運用組合的4/5階龍格-庫塔-芬爾格算法自變量值函數(shù)值用于設定誤差限(缺省時設定相對誤差10-3, 絕對誤差10-6),命令為:options=odeset(reltol,rt,abstol,at), rt,at:分別為設定的相對誤差和絕對誤差. 1、在解n個未知函數(shù)的方程組時,x0和x均為n維向量,m-文件中的待解方程組應以x的分量形式寫成. 2、使用Matlab軟件求數(shù)值解時,高階微分方程必須等價地變換成一階微分方程組.注意注意:例例 4 0)0( ; 2)0(0)1 (1000222xxxdtdxxdtxd解解: 令 y1=x,y2=y1則微分
8、方程變?yōu)橐浑A微分方程組:0)0(, 2)0()1 (1000211221221yyyyyyyy1、建立m-文件vdp1000.m如下: function dy=vdp1000(t,y) dy=zeros(2,1); dy(1)=y(2); dy(2)=1000*(1-y(1)2)*y(2)-y(1); 2、取t0=0,tf=3000,輸入命令: T,Y=ode15s(vdp1000,0 3000,2 0); plot(T,Y(:,1),-)3、結(jié)果如圖050010001500200025003000-2.5-2-1.5-1-0.500.511.52To Matlab(ff4) 例例 5 解微分方程組. 1)0(, 1)0(, 0)0(51. 0321213312321yyyyyyyyyyyy解解 1、建立m-文件rigid.m如下: function dy=rigid(t,y) dy=zeros(3,1); dy(1)=y(2)*y(3); dy(2)=-y(1)*y(3); dy(3)=-0.51*y(1)*y(2);2、取t0=0,tf=12,輸入命令: T,Y=ode45(rigid,0 12,0 1 1);
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
- 4. 未經(jīng)權益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責。
- 6. 下載文件中如有侵權或不適當內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 假期怎么寫學習計劃
- 湖北省武漢市漢南區(qū)職教中心2025年高二物理第二學期期末學業(yè)質(zhì)量監(jiān)測試題含解析
- 企業(yè)年度述職報告?zhèn)€人工作總結(jié)
- 從學校畢業(yè)后參加工作的心得體會
- 路易體癡呆的健康宣教
- 貿(mào)易強國的發(fā)展歷程與趨勢
- 先天性喉隔護理查房
- 扁桃體白喉的健康宣教
- 軀干創(chuàng)傷性切斷的查房
- 腦干梗死的健康宣教
- 工會審計實務課件
- 預防艾滋病、梅毒和乙肝母嬰傳播相關報表、上報流程和要求
- 《鐵路技術管理規(guī)程》(普速鐵路部分)-14年新版
- 食用油儲存期品質(zhì)變化的太赫茲光譜無損識別
- 胎盤早剝預案演練腳本
- 五山文學全集第一卷
- 聚磷腈功能高分子材料的合成及應用
- 中國鐵路總公司《鐵路技術管理規(guī)程》(高速鐵路部分)2014年7月
- 鈣加維生素Dppt課件(PPT 14頁)
- TRD深基坑止水帷幕施工方案(22頁)
- 企業(yè)組織架構圖模板
評論
0/150
提交評論