用MATLAB函数编写并求解微分方程_第1页
用MATLAB函数编写并求解微分方程_第2页
用MATLAB函数编写并求解微分方程_第3页
用MATLAB函数编写并求解微分方程_第4页
用MATLAB函数编写并求解微分方程_第5页
已阅读5页，还剩1页未读， 继续免费阅读

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

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

2、靠性和有效性。三、实验内容 对下图的二阶系统,利用 MA TLAB 编程仿真求解:1. 零状态响应 2. 零输入相应 3. 全响 应四、实验数据根据电路得微分方程 :可得该微分方程的 MA TLAB 函数 (funcforex123.h Function xdot+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;function in=f(t in=(t.0 *1; -0.200.20.40.60.811.21.41.61.8time secf +-u C (t

3、 =+-=12211111x C x fLx Lx LR x -0.2-0.15-0.1-0.0500.050.10.150.20.250.3time sec1. 求解程序:零状态响应 %filename ex123.mL=1;c=0.1;for R=1.5 3 5t,x=ode 45( funcforex123 ,-1,10,0;0,R,L,C; Figure(1;plot(t,x(:,1;hold onl xlibel( time sec ; Text(0.9,0.17, leftarrow i_L(t ;grid;Figure(2;plot(t,x(:,2;hold on; xlabel( time sec ; Text(0.5,0.3, leftarrow u_C(t ;grid; -0.50.511.52time sec -20246810time sec iL (t2. 零输入响应修改 funcforex123.m 中的 in=0,将系统状态改为 0,1(修改 ode45 time sec -0.50.51time sec3. 全响应保持 funcfore

人人文库> 全部分类> 教育资料 > 辅导培训

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论