原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
控制系统AD作业
《MATLAB基础与控制系统CAD》仿真作业
学生学号： 学生姓名：
一、 求如下非线性二阶系统的时间响应（20分）
其中，要求绘出的系统状态响应曲线。
解：在命令行中输入：
[t,y]=ode45(final_01,[0,3],[0,2]);
plot(t,y)
得到如下结果：
注：曲线中的坐标说明，网格等都是用编辑器而非命令行添加的。
其中，final_01(t,y)是自己编写的M文件，代码如下：
function dy=final_01(t,y)
dy=[-y(1)*exp(1-t)+0.8*y(2);y(1)-y(2)^3];
二、 已知系统的开环传递函数如下 （20分）
把G(s)转换成零极点形式的传递函数，判断开环系统稳定性。
判别系统在单位负反馈下的稳定性，并求出闭环系统在0～10秒内的脉冲响应和单位阶跃响应，分别绘出响应曲线。
解：(1)在命令行中输入：
p=roots([1 5 25])
得极点为p=-2.5000 + 4.3301i
-2.5000 - 4.3301i
所以零极点形式的传递函数为：
因为极点都在左半平面，所以系统稳定。
(2)在命令行中输入：
G=tf([10],[1 5 25]);
H=tf([1],[1]);
T=feedback(G,H)
得到闭环传递函数为：
10
T= ------------------
s^2 + 5 s + 35
由roots([1 5 35])得其极点为：p= -2.5000 + 5.3619i
-2.5000 - 5.3619i
因为极点均在左半平面，所以系统稳定。
(3)在命令行中输入：
G=tf([10],[1 5 25]);
H=tf([1],[1]);
T=feedback(G,H);
impulse_and_step(T,10)
得到如下结果：
其中impulse_and_step(T,t)是自己编写的M文件，代码如下：
function impulse_and_step(T,t)
impulse(T,t);
hold on;
step(T,t);
grid on;
title(Impulse and Step Response);
三、 某单位负反馈系统如下图所示，（20分）
当比例控制器增益K＝1时，在Simulink中搭建系统，当输入为阶跃函数时，用示波器观察系统的输出，绘出输出曲线。
把(1)中的对象输出和时钟输出输入Workspace中，通过在命令窗口中执行M文件求出系统在阶跃输入下的超调量()和峰值时间()，写出源程序。
调节控制器增益，使超调量且稳态误差，给出此时K值的范围。
解：(1)在simulink中搭建的系统如下图：
阶跃响应为：
(2)在simulink中搭建的框图如下：
在命令行中输入：
G=tf([10],[1 5 10]);
H=tf([1],[1]);
T=feedback(G,H);
dynamicperformance(T)
得到如下结果：Stablevalue=0.500000,
Peaktime=0.847731, //峰值时间；
Overshoot=12.026425, //超调量；
Delaytime=0.339
文档评论（0）