版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、第五章 控制系统仿真5.2 微分方程求解方法以一个自由振动系统实例为例进行讨论。如下图1所示弹簧-阻尼系统,参数如下:M=5 kg, b=1 N.s/m, k=2 N/m, F=1N x b M F k 图1 弹簧-阻尼系统假设初始条件为:时,将m拉向右方,忽略小车的摩擦阻力, 求系统的响应。)用常微分方程的数值求解函数求解包括ode45、ode23、ode113、ode15s、ode23s等。 wffc1.m myfun1.m一、常微分方程的数值求解函数ode45求解解:系统方程为 这是一个单变量二阶常微分方程。将上式写成一个一阶方程组的形式,这是函数ode45调用规定的格式。令: (位移)
2、 (速度)上式可表示成: 下面就可以进行程序的编制。%写出函数文件myfun1.mfunction xdot=myfun1(t,x)xdot=x(2);1-10*x(2)-20*x(1);% 主程序wffc1.m t=0 30;x0=0;0;tt,xx=ode45(myfun1,t,x0);plot(tt,yy(:,1),:b,tt,yy(:,2),-r)legend(位移,速度)title(微分方程的解 x(t) 二、方法2:%用传递函数编程求解ksys1.mnum=1;den=5 1 2;%printsys(num,den)%t=0:0.1:10;sys=tf(num,den);figur
3、e(1)step(sys)figure(2)impulse(sys)figure(3)t=0:0.1:10;ramp=t;lsim(sys,ramp,t);figure(4)tt=size(t);noise=rand(tt,1);lsim(sys,noise,t)figure(5)yy=0.1*t.2;lsim(num,den,yy,t)w=logspace(-1,1,100);m p=bode(num,den,w);figure(6)subplot(211);semilogx(w,20*log10(m);grid onsubplot(212);semilogx(w,p)grid ongm,p
4、m,wpc,wgc=margin(sys)figure(7)margin(sys)figure(8)nyquist(sys)figure(9)nichols(sys)方法3:% 主程序wffc1.m t=0 30;x0=0;0;tt,yy=ode45(myfun1,t,x0);figure(1)plot(tt,yy(:,1),:b,tt,yy(:,2),-r)hold onplot(tt,0.2-0.2*yy(:,2)-0.4*yy(:,1),-.k)legend(位移,速度,加速度)title(微分方程的解) figure(2)plot(yy(:,1),yy(:,2)title(平面相轨迹)%写出函数文件my
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 人教版九年级物理 16.1电 压(学习、上课课件)
- 护理服务概念
- 带形基础(砖、毛石、砼基础并附图)
- 面向生命的思考作文
- 安徽省亳州市谯城区2024-2025学年三年级上学期第一次月考数学试卷
- 黄南藏族自治州同仁县2024-2025学年六上数学期末质量检测试题含解析
- 黄石市黄石港区2024-2025学年六上数学期末预测试题含解析
- 会理县2024年数学六年级第一学期期末调研模拟试题含解析
- 湖北省十堰市竹溪县第二高级中学2024-2025学年高一上学期9月月考物理试题
- DB34-T 4880.1-2024 食安安徽品牌评价 第1部分:通则
- 城市管理监控指挥中心管理制度
- 胃肠动力学常用检查知识
- 医院信息管理系统软件项目计划书
- 电子体重秤毕业论文
- igcse-物理大纲
- 机床设备搬迁施工方案,北京京城起重吊装搬运公司完整
- 三生事业六大价值
- PE管拉管施工方案
- 屋面防水工程竣工验收报告
- 名著导读星星离我们有多远
- 墙面钢架焊接施工技术交底
评论
0/150
提交评论