自动控制原理校正课程设计单位负反馈系统设计校正

1、课程设计报告题 目 单位负反馈系统设计校正课 程 名 称 自动控制原理课程设计院 部 名 称 机电工程学院 专 业电气工程及其自动化班 级 学 生 姓 名 学 号 课程设计地点 C214 课程设计学时1周指 导 教 师 金陵科技学院教务处制目录1. 绪论·····························

2、;··························31.1相关背景知识······················&

3、#183;····················31.2课程设计目的···························

4、83;···············31.3课程设计任务·································

5、;··········32设计过程·······································

6、;·············42.1确定及校正传递函数··································4 2.

7、1.1确定及校正数····························42.1.2.校正前系统bode图·················

8、3;·············6 2.1.3.校正后系统bode图·······························72.2校正前后分析对比

9、3;·····································8 2.2.1特征根···········

10、·······························9 2.2.2三种响应曲线·················

11、;···················102.2.3动态性能指标····························

12、3;······13根轨迹图·······································152.2.5奈奎斯特图·

13、3;···································203 课程设计总结·············

14、·································214. 参考文献···············

15、83;·································221.绪论1.1相关背景知识自动控制技术已经广泛应用到工业、农业生产，交通运输和国防建设。指导自动控制系统的分析和设计的控制理论也有了很大的发展，它的概念、方法、和体系已经渗透到许多科学领域。在2

16、0世纪40和50年代中发展起来的经典控制理论至今仍然被成功的应用于单变量定常系统的分析和设计中。在20世纪50、60年代初发展起来的状态空间方法更具有广泛的适应性。MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。1.2课程设

17、计目的1.掌握自动控制原理的时域分析法，根轨迹法，频域分析法，以及各种补偿（校正）装置的作用及用法，能够利用不同的分析法对给定系统进行性能分析，能根据不同的系统性能指标要求进行合理的系统设计，并调试满足系统的指标。 2.学会使用MATLAB语言及Simulink动态仿真工具进行系统仿真与调试。1.3课程设计任务 题目：已知单位负反馈系统的开环传递函数，试用频率法设计串联滞后超前校正装置，使之满足在单位斜坡作用下，系统的速度误差系数，系统的相角裕量，校正后的剪切频率 要求：根据给定的性能指标选择合适的校正方式对原系统进行校正（须写清楚校正过程），使其满足工作要求。然后利用MATLAB对未校正系统

18、和校正后系统的性能进行比较分析，针对每一问题分析时应写出程序，输出结果图和结论。最后还应写出心得体会与参考文献等。具体要求如下：1.首先，根据给定的性能指标选择合适的校正方式对原系统进行校正，使其满足工作要求。要求程序执行的结果中有校正装置传递函数和校正后系统开环传递函数，校正装置的参数T，等的值。2.利用MATLAB函数求出校正前与校正后系统的特征根，并判断其系统是否稳定，为什么? 3.利用MATLAB作出系统校正前与校正后的单位脉冲响应曲线，单位阶跃响应曲线，单位斜坡响应曲线，分析这三种曲线的关系？求出系统校正前与校正后的动态性能指标%、tr、tp、ts以及稳态误差的值，并分析 其有何变化

19、？4.绘制系统校正前与校正后的根轨迹图，并求其分离点、汇合点及与虚轴交点的坐标和相应点的增益值，得出系统稳定时增益的变化范围。绘制系统校正前与校正后的Nyquist图，判断系统的稳定性，并说明理由？ 5.绘制系统校正前与校正后的Bode图，计算系统的幅值裕量，相位裕 量，幅值穿越频率和相位穿越频率。判断系统的稳定性，并说明理由？2.设计过程 2.1确定及校正传递函数： 2.1.1确定：由于系统的速度误差系数 ,且 所以，即k0=20 利用MATLAB计算出校正函数：程序如下：（1）确定滞后环节的校正函数（取Wc=2.0rad/s），一般在810之间，取=9>> wc=2.0;k0=

20、20;n=1;>> den=1 3 2 20;>> be=9;>> T=1/(0.1*wc);>> bt=be*T;>> gc1=tf(T 1,bt 1)Transfer function: 5 s + 1-45 s + 1（2）确定超前环节的校正函数：n=conv(0 20,5 1);d=conv(conv(conv(1 0,1 1),1 2),45 1);sope=tf(n,d);wc=2.0;num=sope.num1;den=sope.den1;na=polyval(num,j*wc);da=polyval(den,j*wc)

21、;g=na/da;g1=abs(g);h=20*log10(g1);a=10(h/10);wm=wc;T=1/(wm*(a)(1/2);al=a*T;gc=tf(T 1,al 1)Transfer function: 2.832 s + 1-0.08827 s + 1所以得到滞后超前校正后的系统函数为：所以：T1=45 b=0.1111 T2=0.08827 a=32.082.1.2利用MATLAB绘画校正前系统的bode图：>> num=20;>> den=1 3 2 0;>> sys=tf(num,den);>> margin(sys) %绘

22、制简略波德图 >> hold on>> grid>> hold off 由图中可知：模稳定裕量：Gm=-10.5dB； -穿越频率：Wcg= 1.43rad/sec 相相稳定裕量Pm=28.1deg剪切频率 Wc=2.44rad/sec判断稳定性：因为该系统含有一个积分环节，所以对数相频特性曲线向上补画90，在的范围内对有一次负穿越，又系统无开环极点，根据对数频率稳定判据 ，Z=P-R=1>0，可知该系统不稳定。2.1.3利用MATLAB绘画校正后系统的bode图>> n=20;>> d=conv(conv(1 0,1 1),1

23、 2);>> s=tf(n,d);>> s1=tf(5 1,45 1);>> s2=tf(2.832 1,0.08827 1);>> sope=s*s1*s2;>> mag,pha,w=bode(sope);>> margin(mag,pha,w) 由图中可得：模温度裕量：Gm=14.8dB ； -穿越频率;cg=5.29rad/sec相稳定裕量：Pm=46.6deg ； 剪切频率：cp=2.02ad/sec判断稳定性：含有一个积分环节，向上补画，图中在的范围内对无穿越，又系统无开环不稳定极点，根据对数频率稳定判据 Z=P

24、-R=0，可知该系统稳定。由上可知校正后系统的相角裕度=Pm=45，Wc=2.0rad/s>1.5rad/s 满足题目要求，所以校正装置的传递函数确定为校正后系统开环传递函数为2.2校正前后分析对比 2.2.1特征根 1.利用MATLAB求校正前系统特征根：>> num=20;den=1 3 2 0;g=tf(num,den);sys=feedback(g,1);pzmap(g); %绘制系统零极点>> den=1 3 2 0;>> t=tf(num,den);>> pzmap(t);>> p,z=pzmap(g); %计算并返

25、回系统零点Z和极点P>> den=sys.den1;>> r=roots(den);>> disp(r) -3.8371 0.4186 + 2.2443i 0.4186 - 2.2443i判断稳定性：由于闭环特征根有两个右半平面的根，故系统不稳定2.利用MATLAB求校正后系统特征根： >> num=283.2 156.64 20;den=3.97 57 144.21 376.377 158.64 20;g=tf(num,den);sys=feedback(g,1);pzmap(g);den=1 3 2 0;t=tf(num,den);pzmap

26、(t);p,z=pzmap(g);den=sys.den1;r=roots(den);disp(r) -12.4728 -0.6846 + 3.4019i -0.6846 - 3.4019i -0.2579 + 0.0242i -0.2579 - 0.0242i判断稳定性：由于闭环特征根均在左半平面，故系统稳定。2.2.2响应曲线： （1.1）.校正前单位阶跃响应曲线：>> num=20;den=1 3 2 0;sys=tf(num,den);sys1=feedback(sys,1);>> t=0:0.1:15;>> step(sys1,t) %绘制单位阶跃

27、响应曲线>> hold on>> grid>> hold off（1.2）校正后单位阶跃响应曲线：>> n1=20;d1=conv(conv(1 0,1 1),1 2);s1=tf(n1,d1);s2=tf(5 1,45 1);>> s3=tf(2.832 1,0.08827 1);>> sope=s1*s2*s3;>> sys=feedback(sope,1);>> step(sys)由图知：超调量： 上升时间：tr=0.897-0.25=0.627s峰值时间：tp=1.37s调节时间（5%）：t

28、s=3.2s （2.1）、 校正前单位脉冲响应曲线：num=20;>> den=1 3 2 0;>> sys=tf(num,den);>> closys=feedback(sys,1);>> figure(1);>> impulse(closys); hold on（2.2）校正后单位脉冲响应曲线： k=20;num=conv(5 1,2.832 1);den=conv(conv(conv(conv(1 0,1 1),1 2),45 1),0.08827 1);sys=tf(k*num,den);Lsys=feedback(sys,1

29、,-1);y,t,x=impulse(Lsys);plot(t,y);（3.1）、校正前的斜坡响应曲线：>> t=0:0.02:10;u=t;num=20;den=1 3 2 0;sys=tf(num,den);lsim(sys,u,t);t=0:0.01:10;(3.2)校正后斜坡响应曲线：>> num=conv(5 1,2.932 1);den=conv(conv(conv(conv(1 0,1 1),1 2),45 1),0.08827 1);sys=tf(num,den);lsim(sys,0:0.01:40,0:0.01:40) 三条曲线关系：单位斜坡响应曲线

30、的导数是阶跃信号，阶跃响应曲线的导数是脉冲响应曲线2.2.3动态性能指标： 1. 校正前： 求峰值时间tp>> num=20;>> den=1 3 2 0;>> G1=tf(num,den);>> G2=feedback(G1,1,-1);>> y,t=step(G2);C=dcgain(G2);max_y,k=max(y);tp=t(k) tp = 4.4009>> 求超调量 max_overshoot=100*(max_y-C)/Cmax_overshoot = 506.0398>> 上升时间tr (以从稳

31、态值的10%上升到90%的定义)r1=1;while(y(r1)<0.1*C)r1=r1+1;endr2=1;while(y(r2)<0.9*C)r2=r2+1;endtr=t(r2)-t(r1)tr = 0.4401>> 求调节时间ts(取2%误差带) s=length(t);while y(s)>0.98*C&&y(s)<1.02*Cs=s-1;end>> ts=t(s)ts = 5.9412求稳态误差esss=tf('s') >> g=20/(s*s*s+3*s*s+2*s);>> G

32、c=feedback(g,1) Transfer function: 20-s3 + 3 s2 + 2 s + 20 >> ess=1-dcgain(Gc)ess = 0分析：由上述程序可知： tr = 0.4401，tp =4.4009s，ts =5.9412s，= 506，0398%，由校正前单位阶跃响应曲线图和单位斜坡响应曲线图，可知系统在校正之前是不稳定的。故不存在稳态误差。校正后求峰值时间>> num=283.2 156.64 20;den=conv(conv(conv(conv(1 0,1 1),1 2),45 1),0.08827 1);>>

33、G1=tf(num,den);G2=feedback(G1,1,-1);y,t=step(G2);C=dcgain(G2);max_y,k=max(y);tp=t(k)tp = 1.3806求超调量>> max_overshoot=100*(max_y-C)/Cmax_overshoot =19.5156求上升时间>> r1=1;while(y(r1)<0.1*C)r1=r1+1;endr2=1;while(y(r2)<0.9*C)r2=r2+1;endtr=t(r2)-t(r1)tr = 0.5917求调节时间>> s=length(t);wh

34、ile y(s)>0.98*C&&y(s)<1.02*Cs=s-1;endts=t(s)ts =12.6226求稳态误差>> num=283.2 156.64 20;den=conv(conv(conv(conv(1 0,1 1),1 2),45 1),0.08827 1);G1=tf(num,den);>> Gc=feedback(G1,1)Transfer function: 283.2 s2 + 156.6 s -3.972 s5 + 57 s4 + 144.2 s3 + 376.4 s2 + 158.6 s + 20>>

35、ess=1-dcgain(Gc)ess = 0分析：前经过校正后系统的减小，误差有明显的减小，动态性能变好了。后对比可知校正后tr，tp，变小，tp变大，可知动态性能变好。校正后系统的稳态误差也变小了，比校正前性能好了很多，有更好的适用性。2.2.4根轨迹图：校正前根轨迹：num=20;>> den=1 3 2 0;>> sys=tf(num,den);>> rlocus(sys) %绘制根轨迹（开环增益K从零到无穷大）由上可知分离点处：坐标point = -0.423 +3.37e-0.009i，相应点的增益= 0.0192，与虚轴交点：坐标point =

36、-0.0039+ 1.4i，相应点的增益=0.294所以可知系统稳定的值范围为0.0192<<0.294。2. 校正后： num=283.2 156.64 20;den=conv(conv(conv(conv(1 0,1 1),1 2),45 1),0.08827 1);rlocus(num,den) %计算并返回根r和开环增益k，不绘图。k,p=rlocfind(num,den) %用鼠标选择分离点坐标Select a point in the graphics window -1.4221 - 0.0967ik =0.0536p = -11.3664 -1.4458 + 0.0

37、949i -1.4458 - 0.0949i -0.0465 + 0.0956i -0.0465 - 0.0956i>> k,p=rlocfind(num,den) %与虚轴交点>Select a point in the graphics windowselected_point = -0.0024 + 5.4159ik = 5.5824p = -13.8452 0.0170 + 5.4083i 0.0170 - 5.4083i -0.3288 -0.2111 由上可知分离点处：坐标point = -0.0024 + 5.4159i相应点的增益=0.0536，与虚轴交点：坐标p