用MATLAB函数编写并求解微分方程.doc

《用MATLAB函数编写并求解微分方程.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、实验八：用MATLAB函数编写并求解微分方程一、实验原理为了对连续系统进行方针，首先需要建立其数学模型，然后利用计算机求这些数学模型，从而得出数学模型的数值解。由于连续系统是通过微分方程老建模的，因此对此的计算机方针本质上解释微分方程的数值求解问题。二、实验步骤1.建立数学模型，利用电路的拓扑元件的属性，列出网孔方程或节点方程，并简化为标准形式的计算机可求解的一组微分方程组的过程。2.选择适合的计算机求解方法求解仿真模型。3.编写MATLAB仿真程序式建立Simulink模块方块图，调试并运行程序。4.的初始

2、直接，及仿真结果，对仿真结果进行分析，以确定结果的可靠性和有效性。三、实验内容对下图的二阶系统，利用MATLAB编程仿真求解：1.零状态响应2.零输入相应3.全响应四、实验数据根据电路得微分方程:可得该微分方程的MATLAB函数(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.求解程序：零状态响应%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.零输入响应修改funcforex123.m中的in=0,将系统状态改为[0,1](修改ode45)3.全响应保持funcforex123.m不变(in=(t..0)*1),将系统状态改为[0,1]