P、PI和PID控制器性能分析

1、题 目: P、PI和PID控制器性能分析 初始条件： 一二阶系统结构如图所示，其中系统对象模型为 ， 控制器传递函数为（比例P控制），（比例积分PI控制），（比例积分微分PID控制），令，Di(s)为上述三种控制律之一。RYe+-+W-要求完成的主要任务: （包括课程设计工作量及其技术要求，以及说明书撰写等具体要求）（1） 分析系统分别在P、PI、PID控制器作用下的，由参考输入决定的系统类型及误差常数；（2） 根据（1）中的条件求系统分别在P、PI、PID控制器作用下的、由扰动w(t)决定的系统类型与误差常数； （3） 分析该系统的跟踪性能和扰动性能；（4） 在Matlab中画出（1）和（2

2、）中的系统响应，并以此证明（3）结论；（5） 对上述任务写出完整的课程设计说明书，说明书中必须写清楚计算分析的过程，其中应包括Matlab源程序或Simulink仿真模型，并注释。说明书的格式按照教务处标准书写。时间安排： 任务时间（天）指导老师下达任务书，审题、查阅相关资料2分析、计算2编写程序1撰写报告2论文答辩1指导教师签名： 年 月 日系主任（或责任教师）签名： 年 月 日目录摘要1 设计任务及要求12 设计分析与计算22.1参考输入决定的系统类型及误差常数22.1.1系统类型22.1.2误差常数32.2.扰动W（t）决定的系统类型与误差常数82.2.1 系统类型92.2.2 误差常数

3、102.3系统的跟踪性能和扰动性能112.3.1 跟踪性能112.3.2扰动性能113 仿真程序及波形123.1由输入决定的系统响应123.2由扰动决定的系统响应20总结参考文献摘要比例（P）控制：单独的比例控制也称“有差控制”，输出的变化与输入控制器的偏差成比例关系，偏差越大输出越大。比例控制器对偏差反应较快。实际应用中，比例度的大小应视具体情况而定。比例积分（PI）控制：比例控制不能最终消除余差的缺点限制了它的单独使用。克服余差的办法是在比例控制的基础上加上积分控制作用。积分控制器的输出与输入偏差对时间的积分成正比。只要偏差存在，输出就会不断累积（输出值越来越大或越来越小），一直到偏差为零

4、，累积才会停止。所以，积分控制可以消除余差。两者结合，既有比例控制作用的迅速及时，又有积分控制作用消除余差的能力。比例积分微分（PID）控制：最为理想的控制当属比例-积分-微分控制规律。它集三者之长：既有比例作用的及时迅速，又有积分作用的消除余差的能力，还有微分作用的超前控制功能。本次课程设计分别叙述了系统在系统输入和扰动输入作用下的传递函数，通过传递函数反映系统类型和系统的稳态误差。本次课程设计也充分利用了MATLAB的强大功能，分别对系统的跟踪性能和扰动性能进行了曲线描述。 本次课程设计最大的收获就是运用MATLAB直观地对P、PI、PID三种控制器的控制性能进行了研究。【关键词】 比例

5、积分 微分 稳态性能 系统类型 扰动 跟踪1 设计任务及要求本次课程设计要求对3种控制器在参考输入和扰动输入两种情况下的性能进行研究，并通过MATLAB仿真功能进行验证，主要完成以下3个要求：（1）在参考输入的情况下，分别输出在3种控制器作用下的系统在阶跃、斜坡和加速度3种输入信号作用下的不同响应曲线，同时分析系统类型，求解各个误差系数。（2）在扰动输入的情况下，分别输出在3种控制器作用下的系统在阶跃、斜坡和加速度3种输入信号作用下的不同响应曲线，同时分析系统类型，求解各个误差系数。（3）通过不同响应曲线分析系统的跟踪性能和扰动性能。本次课程设计先确定每个控制器下的传递函数，再通过软件仿真得到

6、响应曲线并进行分析，从而得到系统的性能。2 设计分析与计算2.1参考输入决定的系统类型及误差常数如图2-1所示，当参考输入决定系统类型及误差常数时，Ws=0。RYe+-+W-图2-1 二阶系统方框图由系统方框图可以知道，此二阶系统的误差信号Es=Rs-Y(s) 误差传递函数为：(2-1)es=EsRs =11+D(s)G(s)2.1.1 系统类型（1）比例(P)控制器当传递函数D1=KP=19 时，系统的开环传递函数为： (2-2)D1 sGs=19s+15s+1故此时系统类型为0型系统。(2)比例积分(PI)控制器当传递函数 D2s=KP+K1s=19+0.5s 时，系统的开环传递函数为：(

7、2-3)D2sGs=19s+0.5ss+15s+1故此时系统类型为1型系统。 (3)比例积分微分(PID)控制器当传递函数D3s=KP+KIs+KDs=19+0.5s+419 s 时，系统的开环传递函数为：D3sGs=4s2+361s+9.519ss+1(5s+1)(2-4)故此时系统类型为1型系统。2.1.2 误差常数下面求解在3种控制器作用下的系统在阶跃信号、斜坡信号和加速度信号三种常见的输入信号函数下的误差常数。2.1.2.1 阶跃信号输入阶跃信号为r（t）=R·1（t）,其中R为输入阶跃函数的幅值，取拉氏变换可得Rs=R/s，则有(2-5)Es=11+DsG(s)Rs所以(2

8、-6)ess=lims0sEs=lims0R1+DsGs=R1+lims0DsGs=R1+KP其中KP=lims0DsG(s)对0型系统(2-7)KP=lims0K1+T1s(1+T2s)1+Tas(1+Tbs)=Kess=R1+KP 对1型或高于1型的系统(2-8)KP=limsOK1+T1s(1+T2s)sv1+Tas(1+Tbs)= v1ess=R1+KP=02.1.2.2 斜坡信号输入斜坡信号为rt=Rt,其中R为速度输入函数的斜率，取拉氏变换可得Rs=R/s2,则有(2-9)ES=11+DsG(s)Rs2所以(2-10)ess=lims0sEs=lims0Rs+sDsGs=Rlims

9、0sDsG(s)=KKV其中KV=lims0sDsG(s)(2-11)对0型系统(2-12)KV=limS0sK1+T1s(1+T2s)1+Tas(1+Tbs)=0ess=RKV=对1型系统(2-13)KV=limsOsK1+T1s(1+T2s)s1+Tas(1+Tbs)=Kess=RKV=RK对2型或高于2型的系统(2-14)KV=limsOsK1+T1s(1+T2s)sN1+Tas(1+Tbs)= (N2)ess=RKV=02.1.2.3 加速度信号输入加速度信号为rt=Rt2/2,其中R为加速度信号输入函数的速度变换率，取拉氏变换可得Rs=R/s3,则有(2-15)Es=11+DsG(s

10、)Rs3所以(2-16)ess=lims0sEs=lims0Rs2+s2DsGs=Rlims0s2DsG(s)=RKa其中(2-17)Ka=lims0s2DsG(s)对0型系统(2-18)Ka=lims0s2K1+T1s(1+T2s)1+Tas(1+Tbs)=0ess=RKa=对1型系统Ka=lims0s2K1+T1s(1+T2s)s1+Tas(1+Tbs)=0(2-19)ess=RKa=对2型系统(2-20)Ka=lims0s2K1+T1s(1+T2s)s21+Tas(1+Tbs)=Kess=RKa=RK对3型或高于3型的系统Ka=lims0s2K1+T1s(1+T2s)sN1+Tas(1+

11、Tbs)= (N3)ess=RKa=0(2-21)综上所述：(1)当控制器传递函数时，系统闭环特征方程(2-22) 该系统劳斯表为： s2 5 20 s1 6 s0 20根据劳斯判据第一列全为正，所以此系统稳定。 因为系统是0型系统，开环增益K=19，因此，系统的稳态误差为：ess= 11+KP=R20 , 阶跃输入rt=R×1t,斜坡输 入rt=Rt ,加速度输入rt=Rt22 (2)当控制器传递函数D2s=kP+kIs=19+12s(比例积分PI控制)时，系统特征方程为：(2-23)DS=10s3+12s2+40s+1=0该系统劳斯表为： s3 10 40 s2 12 1 s1

12、236 s0 1根据劳斯判据第一列全为正，所以此系统稳定。因为系统是1型系统，系统开环增益K=12，所以系统稳态误差为：ess()= 0， 阶跃输入 rt=R*1(t) R=2R， 斜坡输入 rt=Rt ， 加速度输入 rt=Rt22(3)当控制器传递函数为D3s=kP+kIs+kDs=19+12s+4s19(比例积分微分PID控制)时，系统特征方程为：(2-24)DS=190s3+236s2+60s+19=0该系统劳斯表为： s3 190 760 s2 236 19 s1 87875118 s0 19根据劳斯判据第一列全为正，所以此系统稳定。因为系统为1型系统,系统开环增益K=12，所以系统

13、稳态误差为：ess()= 0, 阶跃输入 rt=R*1(t)Rk2R, 斜坡输入 rt=Rt , 加速度输入 rt=Rt22综合上述可得到下表。表2-1 输入信号作用下的稳态误差控制器系统型别阶跃输入r(t)=R*1(t)斜坡输入r(t)=Rt加速度输入r(t)=Rt22P0R20PII02RPIDI02R2.2.扰动W（t）决定的系统类型与误差常数由于输入信号和扰动信号作用于系统的不同位置，因此即使系统对于某种形式 输入信号作用的稳态误差为零，但对于同一形式的扰动作用，其稳态误差未必为零。控制系统如图2-1所示，其中W（t）代表扰动信号的拉式变换式。由于在扰动信号W（s）作用下系统的理想输出

14、应为零，故该系统响应扰动W（t）的输出端误差信号为：(2-25)设sE（s）满足终值定理条件，则(2-26)当扰动为阶跃信号时，l=0,则当扰动信号为斜坡信号时，l=1，则当扰动信号为加速度信号时，l=2，则由于误差传递函数所含s=0的零点数，等价于系统扰动作用点前向通道串联积分环节与主反馈通道串联积分环节之和，故对于响应扰动作用的系统有下列结论：（1） 扰动作用点之前的前向通道积分环节数与主反馈通道积分环节数之和决定系统响应扰动作用的型别，该型别与扰动作用点之后前向通道的积分环节数无关。（2） 如果在扰动作用点之前的前向通道或主反馈通道中设置个积分环节，必可消除系统在扰动信号作用下的稳态误差

15、。2.2.1 系统类型因为扰动作用点之前的向前通道积分环节数与主反馈通道积分环节数之和决定系统相应扰动作用的型别，该型别与扰动作用点之后向前通道的积分环节数无关，所以（1）当控制器传递函数为D1s=kP=19(比例P控制)时，系统扰动作用点前的向前通道所含的积分环节数为0，系统主反馈通道所含的积分环节数为0，所以此系统为相应扰动信号的0型系统。（2）当控制器传递函数D2s=kP+kIs=19+12s(比例积分PI控制)时，系统扰动作用点前的向前通道所含的积分环节数为1，系统主反馈通道所含的积分环节数为0，所以此系统为相应扰动信号的I型系统。（3）当控制器传递函数为D3s=kP+kIs+kDs=

16、19+12s+4s19(比例积分微分PID控制)时，系统扰动作用点前的向前通道所含的积分环节数为1，系统主反馈通道所含的积分环节数为0，所以此系统为相应扰动信号的I型系统2.2.2 误差常数控制系统如图2-1所示 ，由于在扰动信号W(s)作用下系统理想输出应为零，故该反馈系统扰动w(t)的输出端误差信号为(2-27)Ens=-Yns=GsW(s)1+DsG(s)下面考虑扰动输入分别在阶跃扰动转矩、斜坡扰动转矩和加速度扰动转矩作用下的稳态误差。（1）当控制器传递函数为D1s=kP=19(比例P控制)时，sEns在s右半平面及虚轴上解析。其中Ens=W(s)s+15s+1+19(2-28)所以系统

17、稳态误差为：ess()=lims0sEn(s)=n020,Ws=n0s,Ws=n1s2,Ws=n2s3（2）当控制器传递函数D2s=kP+kIs=19+12s(比例积分PI控制)时，sEns在s右半平面及虚轴上解析。其中(2-29)Ens=sW(s)ss+15s+1+19s+0.5所以系统稳态误差为：ess()=lims0sEn(s)=0,Ws=n0s2n1,Ws=n1s2,Ws=n2s3(2-30)（3）当控制器传递函数为D3s=kP+kIs+kDs=19+12s+4s19(比例积分微分PID控制)时sE_n (s)在s右半平面及虚轴上解析。其中Ens=sW(s)ss+15s+1+19s+0

18、.5+4s219所以系统稳态误差为：ess()=lims0sEn(s)=0,Ws=n0s2n1,Ws=n1s2,Ws=n2s3综上所述得下表。表2-2不同控制器和输入下的稳态误差和型别控制器阶跃输入W(s)=n0s斜坡输入W(s)=n1s2加速度输入W(s)=n2s3系统型别Pn0200PI02n1PID02n12.3系统的跟踪性能和扰动性能2.3.1 跟踪性能（1）在P控制器作用下能够跟踪阶跃输入，但存在一个稳态位置误差R/20,不能跟踪斜坡输入与加速度输入。（2）在PI控制器作用下系统能很好的跟踪阶跃输入，跟踪斜坡输入存在稳态位置误差2R，不能跟踪加速度输入。（3）在PID控制器作用下系统

19、能很好的跟踪阶跃输入，跟踪斜坡输入存在稳态位置误差2R，不能跟踪加速度输入。2.3.2扰动性能（1）在P控制器作用下系统对阶跃扰动的抵抗性比较好但存在一个扰动误差n019，对斜坡和加速度扰动抵抗能力很差。（2）在PI控制器作用下系统对阶跃和斜坡扰动抵抗力较强，但对加速度扰动抵抗力很差。（3）在PID控制器作用下系统对阶跃和斜坡扰动抵抗力较强，但对加速度扰动抵抗力很差。3 仿真程序及波形3.1由输入决定的系统响应(1) 当控制器传递函数时，系统的开环传递函数为：(3-1)D1 sGs=19s+15s+1系统的闭环传递函数为：(3-2)(s)=19s+15s+1单位阶跃响应的MATLAB程序命令如

20、下：num=19; %分子多项式den=5,6,20; %分母多项式t=0:0.1:10; %时间矢量y,x,t=step(num,den,t); %阶跃响应plot(t,y); %绘制曲线grid; %绘制网格xlabel('t'); %设置横坐标ylabel('y'); %设置纵坐标系统响应曲线图如图3-1所示。图3-1单位阶跃相应单位斜坡响应的MATLAB程序命令如下：num=19; %分子多项式den=5,6,20,0; %分母多项式t=0:0.1:10; %时间矢量y,x,t=step(num,den,t); %阶跃响应plot(t,y); %绘制曲线

21、grid; %绘制网格xlabel('t'); %绘制横坐标ylabel('y'); %绘制纵坐标系统响应曲线图如图3-2所示。图3-2单位斜坡响应单位加速度响应的MATLAB程序命令如下：num=19; %分子多项式den=5,6,20,0,0; %分母多项式t=0:0.1:10; %时间矢量y,x,t=step(num,den,t); %阶跃响应plot(t,y); %绘制曲线grid; %绘制网格xlabel('t'); %绘制横坐标ylabel('y'); %绘制纵坐标系统响应曲线图如图3-3所示。图3-3单位加速度响应(

22、3-3)(2)当控制器传递函数D2s=kP+kIs=19+12s(比例积分PI控制)时，系统闭环传递函数为：s=38s+110s3+12s240s+1单位阶跃响应的MATLAB程序命令如下：num=38,1; %分子多项式den=10,12,40,1; %分母多项式t=0:0.1:10; %时间矢量y,x,t=step(num,den,t); %阶跃响应plot(t,y); %绘制曲线grid; %绘制网格xlabel('t'); %设置横坐标ylabel('y'); %设置纵坐标系统响应曲线图如图3-4所示。图3-4单位阶跃响应单位斜坡响应的MATLAB程序命

23、令如下：num=38,1; %分子多项式den=10,12,40,1,0; %分母多项式t=0:0.1:10; %时间矢量y,x,t=step(num,den,t); %阶跃响应plot(t,y); %绘制曲线grid; %绘制网格xlabel('t'); %设置横坐标ylabel('y'); %设置纵坐标系统响应曲线图如图3-5所示。图3-5单位斜坡响应单位加速度响应的MATLAB程序命令如下：num=38,1; %分子多项式den=10,12,40,1,0,0; %分母多项式t=0:0.1:10; %时间矢量y,x,t=step(num,den,t); %阶

24、跃响应plot(t,y); %绘制曲线grid; %绘制网格xlabel('t'); %设置横坐标ylabel('y'); %设置纵坐标系统响应曲线图如图3-6所示。图3-6单位加速度响应(3-4)（3)当控制器传递函数为D3s=kP+kIs+kDs=19+12s+4s19(比例积分微分PID控制)时，系统闭环传递函数为：s=8s2+722s+19190s3+236s2+760s+19单位阶跃响应的MATLAB程序命令如下：num=8,342,19; %分子多项式den=190,236,760,19; %分母多项式t=0:0.1:10; %时间矢量y,x,t=s

25、tep(num,den,t); %阶跃响应plot(t,y); %绘制曲线grid; %绘制网格xlabel('t'); %设置横坐标ylabel('y'); %设置纵坐标系统响应曲线图如图3-7所示。图3-7单位阶跃响应单位斜坡响应的MATLAB程序命令如下：num=8,342,19; %分子多项式den=190,236,760,19,0; %分母多项式t=0:0.1:10; %时间矢量y,x,t=step(num,den,t); %阶跃响应plot(t,y); %绘制曲线grid; %绘制网格xlabel('t'); %设置横坐标ylabel

26、('y'); %设置纵坐标系统响应曲线图如图3-8所示。图3-8单位斜坡响应单位加速度响应的MATLAB程序命令如下：num=8,342,19; %分子多项式den=190,236,760,19,0,0; %分母多项式t=0:0.1:10; %时间矢量y,x,t=step(num,den,t); %阶跃响应plot(t,y); %绘制曲线grid; %绘制网格xlabel('t'); %设置横坐标ylabel('y'); %设置纵坐标系统响应曲线图如图3-9所示。图3-9单位加速度响应综上分析：对于比例控制系统、比例积分控制系统不能跟踪斜坡输入，

27、对于比例积分微分控制系统能够跟踪斜坡输入；对于比例、比例积分、比例积分微分控制系统都不能跟踪加速度输入。3.2由扰动决定的系统响应(3-5)（1）当控制器传递函数为D1s=kP=19(比例P控制)时，扰动系统闭环传递函数为：s=-G(s)1+GsD(s)=-15s2+6s+20单位阶跃响应的MATLAB程序命令如下：num=-1; %分子多项式den=5,6,20; %分母多项式t=0:0.1:10; %时间矢量y,x,t=step(num,den,t); %阶跃响应plot(t,y); %绘制曲线grid; %绘制网格xlabel('t'); %设置横坐标ylabel(

28、9;y'); %设置纵坐标系统响应曲线图如图3-10所示。图3-10单位阶跃响应单位斜坡响应的MATLAB程序命令如下：num=-1; %分子多项式den=5,6,20,0; %分母多项式t=0:0.1:10; %时间矢量y,x,t=step(num,den,t); %阶跃响应plot(t,y); %绘制曲线grid; %绘制网格xlabel('t'); %设置横坐标ylabel('y'); %设置纵坐标系统响应曲线图如图3-11所示。图3-11单位斜坡响应单位加速度响应的MATLAB程序命令如下：num=-1; %分子多项式den=5,6,20,0,0

29、; %分母多项式t=0:0.1:10; %时间矢量y,x,t=step(num,den,t); %阶跃响应plot(t,y); %绘制曲线grid; %绘制网格xlabel('t'); %设置横坐标ylabel('y'); %设置纵坐标系统响应曲线图如图3-12所示。图3-12单位加速度响应(3-6)(2)当控制器传递函数D2s=kP+kIs=19+12s(比例积分PI控制)时，扰动系统闭环传递函数为：s=-2s10s3+12s2+40s+1单位阶跃响应的MATLAB程序命令如下：num=-2,0; %分子多项式den=10,12,40,1; %分母多项式t=0

30、:0.1:20; %时间矢量y,x,t=step(num,den,t); %阶跃响应plot(t,y); %绘制曲线grid; %绘制网格xlabel('t'); %设置横坐标ylabel('y'); %设置纵坐标系统响应曲线图如图3-13所示。图3-13单位阶跃响应单位斜坡响应的MATLAB程序命令如下：num=-2,0; %分子多项式den=10,12,40,1,0; %分母多项式t=0:0.1:20; %时间矢量y,x,t=step(num,den,t); %阶跃响应plot(t,y); %绘制曲线grid; %绘制网格xlabel('t'); %设置横坐标ylabel(&#