自控-二阶系统Matlab仿真.docxVIP

PAGE PAGE 2 自动控制原理 二阶系统性能分析 Matlab 仿真大作业附题目+ 完整报告内容 设二阶控制系统如图1所示，其中开环传递函数 图1 图2 图3 要求： 1、分别用如图2和图3所示的测速反馈控制和比例微分控制两种方式改善系统的性能，如果要求改善后系统的阻尼比 =0.707，则和分别取多少？ 解： 由得 对于测速反馈控制，其开环传递函数为：； 闭环传递函数为：； 所以当=0.707时，； 对于比例微分控制，其开环传递函数为：； 闭环传递函数为：； 所以当=0.707时，； 2、请用MATLAB分别画出第1小题中的3个系统对单位阶跃输入的响应图； 解： ①图一的闭环传递函数为： ， Matlab代码如下： clc clear wn=sqrt(10); zeta=1/(2*sqrt(10)); t=0:0.1:12; Gs=tf(wn^2,[1,2*zeta*wn,wn^2]); step(Gs,t) title(图一单位阶跃响应曲线); xlabel(t/s);ylabel(c(t)); 响应图如下： ②图二的闭环传递函数为： ， Matlab代码如下： clc clear wn=sqrt(10); zeta=0.707; t=0:0.1:12; Gs=tf(wn^2,[1,2*zeta*wn,wn^2]); step(Gs,t) title(图二单位阶跃响应曲线); xlabel(t/s);ylabel(c(t)); 响应图如下： ③图三的闭环传递函数为： ， Matlab代码如下： clc clear wn=sqrt(10); zeta=0.707; t=0:0.1:12; Gs=tf([0.347*wn^2,wn^2],[1,2*zeta*wn,wn^2]); step(Gs,t) title(图三单位阶跃响应曲线); xlabel(t/s);ylabel(c(t)); 响应图如下： 3、分别求出在单位斜坡输入下，3个系统的稳态误差； 解： ①当时，图一的开环传递函数为： 是I型系统 其中K=10，所以 ②当时，图二的开环传递函数为： 是I型系统 其中K=2.237，所以 ③当时，图三的开环传递函数为： 是I型系统 其中K=10，所以 4、列表比较3个系统的动态性能和稳态性能，并比较分析测速反馈控制和比例微分控制对改善系统性能的不同之处； 解： 可以利用Matlab求峰值时间、超调量、上升时间、调节时间，代码以系统一为例： clc clear wn=sqrt(10); zeta=1/(2*sqrt(10)); t=0:0.1:12; G=tf(wn^2,[1,2*zeta*wn,wn^2]); C=dcgain(G); [y,t]=step(G); plot(t,y); [Y,k]=max(y); timetopeak=t(k) percentovershoot=100*(Y-C)/C n=1; while y(n)C n=n+1; end risetime=t(n) i=length(t) while(y(i)0.98*C)(y(i)1.02*C) i=i-1; end settingtime=t(i) 得到结果如下： 动态性能比较 峰值时间(s) 超调量(﹪) 上升时间(s) 调节时间(s) 系统一 1.0154 60.4417 0.5712 7.2985 系统二 1.4077 4.3253 1.0619 1.8769 系统三 0.8397 12.6740 0.4939 1.5806 稳态性能比较 单位阶跃输入下的稳态误差 系统一 0 系统二 0 系统三 0 由上述数据可以看出，测速反馈控制着重改善系统的平稳性（超调量明显降低），而比例微分控制着重改善系统的快速性（峰值时间、上升时间、调节时间降低）。 5、试用绘制图3对应的系统中参数变化时的根轨迹图，分析变化对系统性能的影响；用MATLAB画出分别为0，0.1，0.2，0.5和1时的系统单位阶跃响应图，比较其动态性能。 解： ①，由特征方程得： ， 此时可利用Matlab编程得到根轨迹 Matlab代码如下： clc clear num=[10 0]; den=[1 1 10]; G=tf(num,den); rlocus(G); title(Td变化的参数根轨迹); xlabel(实轴);ylabel(虚轴); 根轨迹图如下： ②图三的闭环传递函数为： ，，在Td分别取0，0.1，0.2，0.5和1时，可以用for语句实现 Matlab代码如下： clc clear wn=3.1623; zeta=0.1581; t=0:0.1:12; Td=[0,0.1,0.2,0.5,1]; hold on; for i=1:le