資源描述：
《用MATLAB函數(shù)編寫(xiě)并求解微分方程.doc》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在教育資源-天天文庫(kù)。
1、實(shí)驗(yàn)八：用MATLAB函數(shù)編寫(xiě)并求解微分方程一、實(shí)驗(yàn)原理為了對(duì)連續(xù)系統(tǒng)進(jìn)行方針，首先需要建立其數(shù)學(xué)模型，然后利用計(jì)算機(jī)求這些數(shù)學(xué)模型，從而得出數(shù)學(xué)模型的數(shù)值解。由于連續(xù)系統(tǒng)是通過(guò)微分方程老建模的，因此對(duì)此的計(jì)算機(jī)方針本質(zhì)上解釋微分方程的數(shù)值求解問(wèn)題。二、實(shí)驗(yàn)步驟1.建立數(shù)學(xué)模型，利用電路的拓?fù)湓膶傩?，列出網(wǎng)孔方程或節(jié)點(diǎn)方程，并簡(jiǎn)化為標(biāo)準(zhǔn)形式的計(jì)算機(jī)可求解的一組微分方程組的過(guò)程。2.選擇適合的計(jì)算機(jī)求解方法求解仿真模型。3.編寫(xiě)MATLAB仿真程序式建立Simulink模塊方塊圖，調(diào)試并運(yùn)行程序。4.的初始
2、直接，及仿真結(jié)果，對(duì)仿真結(jié)果進(jìn)行分析，以確定結(jié)果的可靠性和有效性。三、實(shí)驗(yàn)內(nèi)容對(duì)下圖的二階系統(tǒng)，利用MATLAB編程仿真求解：1.零狀態(tài)響應(yīng)2.零輸入相應(yīng)3.全響應(yīng)四、實(shí)驗(yàn)數(shù)據(jù)根據(jù)電路得微分方程:可得該微分方程的MATLAB函數(shù)(funcforex123.h)Functionxdot+funcforex123(t,x,flag,R,L,C)xdot+zeros(2,1);xdot(1)=-R/L*x(1)-1/L*x(2)+1/L*f(t);xdot(2)=1/C*x(1);functionin=f(t)in
3、=（t.0）*1;1.求解程序：零狀態(tài)響應(yīng)%filenameex123.mL=1;c=0.1;forR=[1.535][t,x]=ode45(‘funcforex123’,[-1,10],[0;0],[],R,L,C);Figure(1);plot(t,x(:,1));holdonlxlibel(‘timesec’);Text(0.9,0.17,’leftarrowi_L(t)’);grid;Figure(2);plot(t,x(:,2));holdon;xlabel(‘timesec’);Text(0.5
4、,0.3,’leftarrowu_C(t)’);grid;2.零輸入響應(yīng)修改funcforex123.m中的in=0,將系統(tǒng)狀態(tài)改為[0,1](修改ode45)3.全響應(yīng)保持funcforex123.m不變(in=(t..0)*1),將系統(tǒng)狀態(tài)改為[0,1]
