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

题目:P、PI和PID控制器性能分析初始条件：一二阶系统结构如图所示，其中系统对象模型为））（（）（15s1s1sG，控制器传递函数为Pk)s(D1（比例P控制），/skkIP)s(D2（比例积分PI控制），sk/skkDIP)s(D3（比例积分微分PID控制），令19Pk，5.0Ik，19/4Dk，Di(s)为上述三种控制律之一。RYe+-+W-要求完成的主要任务:（包括课程设计工作量及其技术要求，以及说明书撰写等具体要求）（1）分析系统分别在P、PI、PID控制器作用下的，由参考输入决定的系统类型及误差常数；（2）根据（1）中的条件求系统分别在P、PI、PID控制器作用下的、由扰动w(t)决定的系统类型与误差常数；（3）分析该系统的跟踪性能和扰动性能；（4）在Matlab中画出（1）和（2）中的系统响应，并以此证明（3）结论；（5）对上述任务写出完整的课程设计说明书，说明书中必须写清楚计算分析的过程，其中应包括Matlab源程序或Simulink仿真模型，并注释。说明书的格式按照教务处标准书写。时间安排：任务时间（天）指导老师下达任务书，审题、查阅相关资料2分析、计算2编写程序1撰写报告2论文答辩1指导教师签名：年月日系主任（或责任教师）签名：年月日目录摘要1设计任务及要求.................................................22设计分析与计算.................................................32.1参考输入决定的系统类型及误差常数....................................32.1.1系统类型....................................................32.1.2误差常数....................................................42.2.扰动W（t）决定的系统类型与误差常数.................................92.2.1系统类型.....................................................102.2.2误差常数.....................................................112.3系统的跟踪性能和扰动性能...........................................122.3.1跟踪性能.....................................................122.3.2扰动性能.....................................................133仿真程序及波形................................................133.1由输入决定的系统响应...............................................133.2由扰动决定的系统响应...............................................21总结参考文献摘要比例（P）控制：单独的比例控制也称“有差控制”，输出的变化与输入控制器的偏差成比例关系，偏差越大输出越大。比例控制器对偏差反应较快。实际应用中，比例度的大小应视具体情况而定。比例积分（PI）控制：比例控制不能最终消除余差的缺点限制了它的单独使用。克服余差的办法是在比例控制的基础上加上积分控制作用。积分控制器的输出与输入偏差对时间的积分成正比。只要偏差存在，输出就会不断累积（输出值越来越大或越来越小），一直到偏差为零，累积才会停止。所以，积分控制可以消除余差。两者结合，既有比例控制作用的迅速及时，又有积分控制作用消除余差的能力。比例积分微分（PID）控制：最为理想的控制当属比例-积分-微分控制规律。它集三者之长：既有比例作用的及时迅速，又有积分作用的消除余差的能力，还有微分作用的超前控制功能。本次课程设计分别叙述了系统在系统输入和扰动输入作用下的传递函数，通过传递函数反映系统类型和系统的稳态误差。本次课程设计也充分利用了MATLAB的强大功能，分别对系统的跟踪性能和扰动性能进行了曲线描述。本次课程设计最大的收获就是运用MATLAB直观地对P、PI、PID三种控制器的控制性能进行了研究。【关键词】比例积分微分稳态性能系统类型扰动跟踪1设计任务及要求本次课程设计要求对3种控制器在参考输入和扰动输入两种情况下的性能进行研究，并通过MATLAB仿真功能进行验证，主要完成以下3个要求：（1）在参考输入的情况下，分别输出在3种控制器作用下的系统在阶跃、斜坡和加速度3种输入信号作用下的不同响应曲线，同时分析系统类型，求解各个误差系数。（2）在扰动输入的情况下，分别输出在3种控制器作用下的系统在阶跃、斜坡和加速度3种输入信号作用下的不同响应曲线，同时分析系统类型，求解各个误差系数。（3）通过不同响应曲线分析系统的跟踪性能和扰动性能。本次课程设计先确定每个控制器下的传递函数，再通过软件仿真得到响应曲线并进行分析，从而得到系统的性能。2设计分析与计算2.1参考输入决定的系统类型及误差常数如图2-1所示，当参考输入决定系统类型及误差常数时，W(s)=0。图2-1二阶系统方框图由系统方框图可以知道，此二阶系统的误差信号E(s)=R(s)-Y(s)误差传递函数为：Φe(s)=E(s)R(s)=11+D(s)G(s)2.1.1系统类型（1）比例(P)控制器当传递函数D1=KP=19时，系统的开环传递函数为：D1(s)G(s)=19(s+1)(5s+1)故此时系统类型为0型系统。(2)比例积分(PI)控制器当传递函数D2(s)=KP+K1s=19+0.5s时，系统的开环传递函数为：D2(s)G(s)=19s+0.5s(s+1)(5s+1)故此时系统类型为1型系统。RYe+-+W-(2-1)(2-2)(2-3)(3)比例积分微分(PID)控制器当传递函数D3(s)=KP+KIs+KDs=19+0.5s+419s时，系统的开环传递函数为：D3(s)G(s)=4s2+361s+9.519s(s+1)(5s+1)故此时系统类型为1型系统。2.1.2误差常数下面求解在3种控制器作用下的系统在阶跃信号、斜坡信号和加速度信号三种常见的输入信号函数下的误差常数。2.1.2.1阶跃信号输入阶跃信号为r（t）=R·1（t）,其中R为输入阶跃函数的幅值，取拉氏变换可得R(s)=R/s，则有E(s)=11+D(s)G(s)Rs所以ess(∞)=lims→0sE(s)=lims→0R1+D(s)G(s)=R1+lims→0D(s)G(s)=R1+KP其中KP=lims→0D(s)G(s)对0型系统KP=lims→0K(1+T1s)(1+T2s)⋯(1+Tas)(1+Tbs)⋯=Kess(∞)=R1+KP对1型或高于1型的系统(2-4)(2-5)(2-6)(2-7)KP=lims→OK(1+T1s)(1+T2s)⋯sv(1+Tas)(1+Tbs)⋯=∞v≫1ess(∞)=R1+KP=02.1.2.2斜坡信号输入斜坡信号为r(t)=Rt,其中R为速度输入函数的斜率，取拉氏变换可得R(s)=R/s2,则有E(S)=11+D(s)G(s)Rs2所以ess(∞)=lims→0sE(s)=lims→0Rs+sD(s)G(s)=Rlims→0sD(s)G(s)=KKV其中KV=lims→0sD(s)G(s)对0型系统KV=limS→0sK(1+T1s)(1+T2s)⋯(1+Tas)(1+Tbs)⋯=0ess(∞)=RKV=∞对1型系统KV=lims→OsK(1+T1s)(1+T2s)⋯s(1+Tas)(1+Tbs)⋯=Kess(∞)=RKV=RK对2型或高于2型的系统(2-8)(2-9)(2-10)(2-11)(2-12)(2-13)(2-14)KV=lims→OsK(1+T1s)(1+T2s)⋯sN(1+Tas)(1+Tbs)⋯=∞(N≫2)ess(∞)=RKV=02.1.2.3加速度信号输入加速度信号为r(t)=Rt2/2,其中R为加速度信号输入函数的速度变换率，取拉氏变换可得R(s)=R/s3,则有E(s)=11+D(s)G(s)Rs3所以ess(∞)=lims→0sE(s)=lims→0Rs2+s2D(s)G(s)=Rlims→0s2D(s)G(s)=RKa其中Ka=lims→0s2D(s)G(s)对0型系统Ka=lims→0s2K(1+T1s)(1+T2s)⋯(1+Tas)(1+Tbs)⋯=0ess(∞)=RKa=∞对1型系统Ka=lims→0s2K(1+T1s)(1+T2s)⋯s(1+Tas)(1+Tbs)⋯=0ess(∞)=RKa=∞对2型系统(2-15)(2-16)(2-17)(2-18)(2-19)(2-20)Ka=lims→0s2K(1+T1s)(1+T2s)⋯s2(1+Tas)(1+Tbs)⋯=Kess(∞)=RKa=RK对3型或高于3型的系统Ka=lims→0s2K(1+T1s)(1+T2s)⋯sN(1+Tas)(1+Tbs)⋯=∞(N≫3)ess(∞)=RKa=0综上所述：(1)当控制器传递函数()19==1pKsD时，系统闭环特征方程()0=20+6+5=2sssD该系统劳斯表为：s2520s16s020根据劳斯判据第一列全为正，所以此系统稳定。因为系统是0型系统，开环增益K=19，因此，系统的稳态误差为：ess(∞)={11+KP=R20,阶跃输入r(t)=R×1(t)∞,斜坡输入r(t)=Rt∞,加速度输入r(t)=Rt22(2)当控制器传递函数D2(s)=kP+kIs=19+12s(比例积分PI控制)时，系统特征方程为：D(S)=10s3+12s2+40s+1=0该系统劳斯表为：s31040s2121s1236s01(2-21)(2-22)(2-23)根据劳斯判据第一列全为正，所以此系统稳定。因为系统是1型系统，系统开环增益K=12，所以系统稳态误差为：ess(∞)={0，阶跃输入r(t)=R*1(t)RＫ=2R，斜坡输入r(t)=Rt∞，加速度输入r(t)=Rt22(3)当控制器传递函数为D3(s)=kP+kIs+kDs=19+12s+4s19(比例积分微分PID控制)时，系统特征方程为：D(S)=190s3+236s2+60s+19=0该系统劳斯表为：s3190760s223619s187875118s019根据劳斯判据第一列全为正，所以此系统稳定。因为系统为1型系统,系统开环增益K=12，所以系统稳态误差为：ess(∞)={0,阶跃输入r(t)=R*1(t)Rk2R,斜坡输入r(t)=Rt∞,加速度输入r(t)=Rt22(2-24)sBsBsBKKKsAsAsAsssBsBsAsKAKsNsGsDsGsEvvvlviN32132132113213232111综合上述可得到下表。表2-1输入信号作用下的稳态误差控制器系统型别阶跃输入r(t)=R*1(t)斜坡输入r(t)=Rt加速度输入r(t)=Rt22P0R20∞∞PII02R∞PIDI02R∞2.2.扰动W（t）决定的系统类型与误差常数由于输入信号和扰动信号作用于系统的不同位置，因此即使系统对于某种形式输入信号作用的稳态误差为零，但对于同一形式的扰动作用，其稳态误差未必为零。控制系统如图2-1所示，其中W（t）代表扰动信号的拉式变换式。由于在扰动信号W（s）作用下系统的理想输出应为零，故该系统响应扰动W（t）的输出端误差信号为：设sE（s）满足终值定理条件，则()()()()()()()()()()()[]sBsBsBKKKsAsAsAsssBsBs