版權(quán)說(shuō)明:本文檔由用戶(hù)提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
1第9章
MATLAB微分方程求解
29.1常微分方程初值問(wèn)題常微分方程的初值問(wèn)題是科學(xué)計(jì)算中的常見(jiàn)問(wèn)題,MATLAB提供了求解常微分方程(ODE)初值問(wèn)題的一系列ode函數(shù)。ode函數(shù)主要采用Runge-Kutta法求解常微分方程。常用ode函數(shù)函數(shù)名含義函數(shù)名含義ode23普通2-3階法解ODEode23s變階法解剛性odeode45普通4-5階法解ODEode23t解適度剛性odeode113普通變階法解ODEode23tb低階法解剛性ode39.1.1ODE的機(jī)理將高階常微分方程的初值問(wèn)題,轉(zhuǎn)變?yōu)?階常微分方程組的初值問(wèn)題。9.1.2普通2-3階法解ODE例9-1編制常微分方程的函數(shù)文件dfun.mfunctiondy=dfun(x,y)dy=y-2*x/y編制如下計(jì)算程序clear;clc;X12=linspace(0,5,50);[X,Y]=ode23('dfun',X12,2);plot(X,Y)運(yùn)行結(jié)果如圖9-1所示圖9-149.1.3普通4-5階法解ODE例9-2編制微分方程組的函數(shù)文件dfun.mfunctiondy=dfun(x,y)dy(1,1)=-2*y(1)+y(2)+2*sin(x);dy(2,1)=10*y(1)-9*y(2)+9*(cos(x)-sin(x));編制如下程序clear;clc;X12=linspace(0,9,30);ode45('dfun',X12,[2,3]);運(yùn)行結(jié)果如圖9-2所示圖9-259.2常微分方程邊值問(wèn)題9.2.1邊值問(wèn)題簡(jiǎn)介一般來(lái)說(shuō),微分方程邊值問(wèn)題可能有解、也可能無(wú)解,可能有唯一解、也可能有無(wú)數(shù)解。在假定有唯一解的前提下,邊值問(wèn)題有三種基本解法:迭加法:假如微分方程和邊界條件是線(xiàn)性的,邊值問(wèn)題可以轉(zhuǎn)化為初值問(wèn)
題。試射法:將問(wèn)題轉(zhuǎn)化為對(duì)漏缺初值的搜索,一旦漏缺初值確定,就可以按
初值問(wèn)題求解。松弛法:先猜測(cè)滿(mǎn)足邊界條件的區(qū)間網(wǎng)點(diǎn)上的解值,然后利用微分方程進(jìn)
行迭代改善。69.2.2bvp4c求解思路為求解常微分方程的邊值問(wèn)題,MATLAB提供了bvp4c函數(shù)文件。該文件是依據(jù)有限元法中的配置法編寫(xiě)的。其解題的基本思路如下:
把待解問(wèn)題轉(zhuǎn)化為標(biāo)準(zhǔn)邊值問(wèn)題為期望解指定初始猜測(cè)利用原微分方程,構(gòu)造殘差函數(shù)用原微分方程和邊界條件,借助迭代使殘差不斷減小為了bvp4c使用的方便,MATLAB提供了能產(chǎn)生較好猜測(cè)解的bvpinit指令,減輕了用戶(hù)進(jìn)行初始猜測(cè)的負(fù)擔(dān)。此外解算指令bvp4c也能處理包含未知參數(shù)的微分方程的邊值問(wèn)題。79.2.3bvp4c配套指令例9-3編制微分方程組的函數(shù)文件dfun.mfunctiondy=dfun(x,y,c)dy(1,1)=y(2);dy(2,1)=-c*abs(y(1));編制殘差方程組的函數(shù)文件bfun.mfunctionres=bfun(ya,yb,c)res(1,1)=ya(1);res(2,1)=yb(1)+2編制如下計(jì)算程序sol=bvp4c('dfun','bfun',sinit,[],c);plot(sol.x,sol.y(1,:),'r--')holdonc=1.5;sol=bvp4c('dfun','bfun',sinit,[],c);plot(sol.x,sol.y(1,:),'b-')legend('c=1','c=1.5')運(yùn)行結(jié)果如圖9-3所示clear;clc;x=linspace(0,4,5);v=[1;0];sinit=bvpinit(x,v);c=1;圖9-389.3常微分方程的解析解9.3.1微分方程的通解例9-4編制如下程序clear;clc;fun='Dx=-a*t'x=dsolve(fun,'t')運(yùn)行結(jié)果為fun=Dx=-a*t
x=
-1/2*a*t^2+C1例9-5編制如下程序clear;clc;fun1='Df=f+g';fun2='Dg=g-f';[f,g]=dsolve(fun1,fun2,'t')運(yùn)行結(jié)果為f=
-exp(t)*(C1*cos(t)-C2*sin(t))
g=
exp(t)*(C1*sin(t)+C2*cos(t))99.3.2微分方程的特解例9-6編制如下程序clear;clc;fun='D2y=4'x=dsolve(fun,'y(0)=1','y(1)=3','x')運(yùn)行結(jié)果為fun=D2y=4
x=
2*x^2+1例9-7編制如下程序clear;clc;fun1='Df=f+g';fun2='Dg=g-f';v1='f(0)=2';v2='g(0)=2';[f,g]=dsolve(fun1,fun2,v1,v2,'t')運(yùn)行結(jié)果為f=
exp(t)*(2*cos(t)+2*sin(t))
g=
exp(t)*(-2*sin(t)+2*cos(t))10例9-8編制如下程序clear;clc;fun='D3y=-y';v1='y(0)=3';v2='Dy(0)=0';v3='D2y(0)=0';y=dsolve(fun,v1,v2,v3,'t')運(yùn)行結(jié)果為y=
exp(-t)+2*exp(1/2*t)*cos(1/2*3^(1/2)*t)119.4偏微分方程求解9.4.1偏微分方程簡(jiǎn)介很多科學(xué)計(jì)算問(wèn)題都可以歸結(jié)為求解含定解條件(初始條件和邊界條件)的偏微分方程問(wèn)題。MATLAB提供了專(zhuān)門(mén)求解偏微分方程問(wèn)題工具箱(PDEToolbox)。MATLAB可以求解的線(xiàn)性偏微分方程類(lèi)型包括:
橢圓型方程拋物線(xiàn)方程雙曲型方程特征值方程129.4.2偏微分方程的求解過(guò)程利用MATLAB求解偏微分方程的過(guò)程如下:1.完成幾何造型2.輸入定解條件3.選擇方程類(lèi)型4.劃分網(wǎng)格5.數(shù)值求解6.結(jié)果可視化13例9-91415169.4.3PDE圖形用戶(hù)界面簡(jiǎn)介例9-10179.5綜合實(shí)例9.5.1有阻尼自由振動(dòng)例9-11編制標(biāo)準(zhǔn)一階常微分方程組的函數(shù)文件dfun.mfunctiondz=dfun(t,z)a=-2;b=-100;dz(1,1)=z(2);dz(2,1)=a*z(2)+b*z(1);編制如下程序clear;clc;t12=linspace(0,2,100);[t,z]=ode45('dfun',t12,[1,0]);plot(t,z(:,1))gridon運(yùn)行后振動(dòng)規(guī)律如圖9-19所示圖9-19189.5.2點(diǎn)的運(yùn)動(dòng)軌跡例9-12編制標(biāo)準(zhǔn)微分方程組的函數(shù)文件dfun.mfunctiondu=dfun(t,u)du(1,1)=u(3);du(2,1)=u(4);du(3,1)=0.2*u(1);du(4,1)=0.2*u(2);編制如下程序clear;clc;t12=linspace(0,2,100);[t,u]=ode45('dfun',t12,[2,0,0,1])plot(u(:,1),u(:,2))gridon運(yùn)行后,得到軌跡如圖9-20所示圖9-20199.5.3梁的撓曲線(xiàn)編制標(biāo)準(zhǔn)微分方程組的函數(shù)文件dfun.mfunctiondu=dfun(x,u)du(1,1)=u(2);du(2,1)=-0.02*x+0.01*x^3;根據(jù)定解條件,編制殘差方程組的函數(shù)文件bfun.mfunctionres=bfun(ya,yb)res(1,1)=ya(1);r
溫馨提示
- 1. 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶(hù)所有。
- 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ì)用戶(hù)上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶(hù)上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶(hù)因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 2024物業(yè)公司承擔(dān)住宅小區(qū)垃圾清運(yùn)的合同
- 2025年度留置車(chē)輛處置借款合同4篇
- 2025年grc構(gòu)件生產(chǎn)線(xiàn)投資建設(shè)與運(yùn)營(yíng)合同3篇
- 年度PAPTFE競(jìng)爭(zhēng)策略分析報(bào)告
- 年度童書(shū)產(chǎn)業(yè)分析報(bào)告
- 2024-2025學(xué)年新教材高中語(yǔ)文基礎(chǔ)過(guò)關(guān)訓(xùn)練15諫逐客書(shū)含解析部編版必修下冊(cè)
- 二零二五版白糖倉(cāng)儲(chǔ)物流服務(wù)合同范本2篇
- 2025年理療項(xiàng)目合作協(xié)議范本:特色理療項(xiàng)目合作框架協(xié)議3篇
- 2025年度中小企業(yè)間資金周轉(zhuǎn)互助合同范本
- 二零二五年度商業(yè)地產(chǎn)租賃合同中情勢(shì)變更處理辦法及責(zé)任劃分4篇
- 骨科手術(shù)后患者營(yíng)養(yǎng)情況及營(yíng)養(yǎng)不良的原因分析,骨傷科論文
- GB/T 24474.1-2020乘運(yùn)質(zhì)量測(cè)量第1部分:電梯
- GB/T 12684-2006工業(yè)硼化物分析方法
- 定崗定編定員實(shí)施方案(一)
- 高血壓患者用藥的注意事項(xiàng)講義課件
- 特種作業(yè)安全監(jiān)護(hù)人員培訓(xùn)課件
- (完整)第15章-合成生物學(xué)ppt
- 太平洋戰(zhàn)爭(zhēng)課件
- 封條模板A4打印版
- T∕CGCC 7-2017 焙烤食品用糖漿
- 貨代操作流程及規(guī)范
評(píng)論
0/150
提交評(píng)論