- 1、本文档共5页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
微分方程几种求解方法
第五章 控制系统仿真
§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调用规定的格式。
令: (位移)
(速度)
上式可表示成:
下面就可以进行程序的编制。
%写出函数文件myfun1.m
function xdot=myfun1(t,x)
xdot=[x(2);0.2-0.2*x(2)-0.4*x(1)];
% 主程序wffc1.m
t=[0 30];
x0=[0;0];
[tt,yy]=ode45(@myfun1,t,x0);
plot(tt,yy(:,1),':b',tt,yy(:,2),'-r')
hold on
plot(tt,0.2-0.2*yy(:,2)-0.4*yy(:,1),'-k')
legend('位移','速度',’加速度’)
title('微分方程的解 x(t)')
xlabel(‘时间t’)
ylabel(‘输出’)
二、方法2:
%用传递函数编程求解ksys1.m
num=1;
den=[5 1 2];
%printsys(num,den)
%t=0:0.1:10;
sys=tf(num,den);
figure(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)';
figure(6)
bode(num,den,w)
grid on
[m p]=bode(num,den,w);
[gm,pm,wpc,wgc]=margin(sys)
figure(7)
margin(sys)
figure(8)
nyquist(sys)
figure(9)
nichols(sys)
方法3:Simulink中传递函数模块求解
方法4:Simulink中积分模块求解
您可能关注的文档
- 2012年卫生保健健康教育工作计划及总结.doc
- 2012级生物技术58-60班卫生统计学考试卷.doc
- 2014年1至11月份配网跳闸分析报告.doc
- 2014年7月份浙江省学考真题.doc
- 2015产科指南解读.doc
- 2015年4月自学考试历年真题科目代码03044.doc
- 2015年下半年网站构造技术第1次作业答案.doc
- 2015年中学卫生工作计划.doc
- 2015年医院工作总结及2016年工作计划.doc
- 2015年辅导员职业能力大赛案例分析题库.doc
- 第12课 我们小点儿声 课件 二年级道德与法治上册(部编版).ppt
- 11.2我从哪里来(教学课件)二年级道德与法治下册(统编版).ppt
- 第10课 我们不乱扔 课件 二年级道德与法治上册(部编版).ppt
- 1.3过好我们的课余生活 课件五年级道德与法治上册(部编版).ppt
- 第四单元《法律保护我们健康成长》大单元整体学程设计道德与法治六年级上册统编版.pdf
- 第十一课:多姿多彩的民间艺术(分层练习)四年级道法下册 部编版.pdf
- 第八课:大家的“朋友”(分层练习)三年级道法下册 部编版.pdf
- 第5课 我爱我们班 课件 二年级道德与法治上册(部编版).ppt
- 第二单元 我们是公民 大单元整体学程设计道德与法治六年级上册统编版.pdf
- 人教部编版二年级语文下册第五单元单元教学课件.ppt
文档评论(0)