控制系统CAD作业.docVIP

《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.339092, Risetime 0.508638, Setllingtime 1.271596. 同时自动绘出阶跃响应曲线如下： 其中dynamicperformance T 是我自己编写的求解动态性能指标的通用函数，详细代码 如下： function dynamicperformance T [y,t] step T ; plot t,y grid on; Stablevalue dcgain T ; [Ymax,k] max y ; if Ymax Stablevalue Peaktime t k ; Overshoot 100* Ymax-Stablevalue /Stablevalue; else %空语句，无操作； end i 1; while y i 0.9*Stablevalue i i+1; end Risetime t i ; i length t ; while abs y i -Stablevalue 0.02*Stablevalue i i-1; end Setllingtime t i ; i 1; while y i 0.5*Stablevalue i i+1; end Delaytime t i ; if Ymax Stablevalue sprintf Stablevalue %f,\nPeaktime %f,\nOvershoot %f, \nDelaytime %f,\nRisetime %f,\nSetllingtime %f.\n, Stablevalue,Peaktime,Overshoot,D