第七章控制系统设计方法ppt课件.ppt

《第七章控制系统设计方法ppt课件.ppt》由会员分享，可在线阅读，更多相关《第七章控制系统设计方法ppt课件.ppt（204页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、控制系统的设计过程可以在时域进行，也可以控制系统的设计过程可以在时域进行，也可以在频域进行。在频域进行。如果对象模型是以传递函数的形式给出，通常如果对象模型是以传递函数的形式给出，通常采用经典控制理论中的频率特性法或根轨迹法采用经典控制理论中的频率特性法或根轨迹法完成控制器的设计，即在原有系统中引入适当完成控制器的设计，即在原有系统中引入适当的环节，用以对原有系统的某些性能（如相角的环节，用以对原有系统的某些性能（如相角裕度、剪切频率、误差系数等）进行校正，使裕度、剪切频率、误差系数等）进行校正，使校正后的系统达到期望的性能要求。校正后的系统达到期望的性能要求。如果对象模型是在状态空间以状态方

2、程形如果对象模型是在状态空间以状态方程形式描述的，则系统的设计过程是在时域进行的，式描述的，则系统的设计过程是在时域进行的，通常是采用状态反馈和极点配置的方法得到控通常是采用状态反馈和极点配置的方法得到控制策略，其中包括状态观测器的设计以及最优制策略，其中包括状态观测器的设计以及最优控制系统的设计等，其研究内容习惯上称为现控制系统的设计等，其研究内容习惯上称为现代控制理论。代控制理论。本章主要以线性时不变系统为对象讨论几种本章主要以线性时不变系统为对象讨论几种常用的设计方法，包括常用的设计方法，包括串联校正、串联校正、PIDPID控制器设计、控制器设计、极点配置与状态观测器设计、极点配置与状态

3、观测器设计、线性二次型最优控制系统设计等。线性二次型最优控制系统设计等。7.1 7.1 串联校正串联校正 本节讨论三种串联校正装置的频域设计方法，即本节讨论三种串联校正装置的频域设计方法，即相位超前、相位滞后、相位滞后超前校正装置设计。相位超前、相位滞后、相位滞后超前校正装置设计。相位超前校正相位超前校正主要用于改善闭环系统的动态特性，主要用于改善闭环系统的动态特性，对于系统的稳态精度影响较小；对于系统的稳态精度影响较小；相位滞后校正相位滞后校正可以明显地改善系统的稳态性能，但可以明显地改善系统的稳态性能，但会使动态响应过程变缓；会使动态响应过程变缓；相位滞后超前校正相位滞后超前校正则把两者的

4、校正特性结合起来，则把两者的校正特性结合起来，用于动态、静态特性均要求较高的系统。用于动态、静态特性均要求较高的系统。下面具体讨论三种校正装置的设计与实现问题下面具体讨论三种校正装置的设计与实现问题7.1.1 7.1.1 相位超前校正相位超前校正相位超前校正环节可以等效相位超前校正环节可以等效地由电阻电容构成的地由电阻电容构成的RCRC网络网络来表示。来表示。其网络传递函数可以写为其网络传递函数可以写为其中其中z=1/T 说明说明 超前校正环节具有极点超前校正环节具有极点 ,零点零点 ，由于，由于 ，因此在，因此在s s平面极点位于零点平面极点位于零点的左侧。的左侧。例例7-17-1使用使用M

5、ATLABMATLAB绘制当绘制当 时时的的bodebode图和图和NyquistNyquist图图T=1，MATLAB 程序如下：程序如下：alpha0=0.1T=1for i=1:5 alpha(i)=i*alpha0 G(i)=tf(alpha(i)*T alpha(i),alpha(i)*T 1)endbode(G(1),G(2),G(3),G(4),G(5)figurenyquist(G(1),G(2),G(3),G(4),G(5)运行程序，得到结果：运行程序，得到结果：由图可知：由图可知：1 1）最大超前角）最大超前角 与所对应的频率与所对应的频率 随随 的减的减 小而升高，并有关

6、系式小而升高，并有关系式2 2）处于两个转折频率处于两个转折频率 和和1/T1/T的几何中心，的几何中心，即即3 3）超前校正环节提供的最大相位超前角约在）超前校正环节提供的最大相位超前角约在 之间。若需要更大的超前角，可以采用多个超前之间。若需要更大的超前角，可以采用多个超前环节的串联。环节的串联。或或说明：由于此网络的增益为，所以在实际说明：由于此网络的增益为，所以在实际应用时，为保证系统的稳态性能，必须增应用时，为保证系统的稳态性能，必须增加一个增益为的放大器，即校正网络为加一个增益为的放大器，即校正网络为例例7-2 7-2 已知原系统开环传递函数为已知原系统开环传递函数为试（试（1 1

7、）用）用bodebode图设计超前校正装置，设计指标为：图设计超前校正装置，设计指标为：静态速度误差系数静态速度误差系数 相角裕度相角裕度（2 2）用）用MATLABMATLAB语言绘制校正前后的语言绘制校正前后的bodebode图及单图及单位阶跃响应。位阶跃响应。理论分析理论分析根据根据 可以求得校正环节的增益可以求得校正环节的增益MATLAB 程序如下：程序如下：ng=400dg=1 30 200 0G0=tf(ng,dg)kc=5dPm=40+10mag,phase,w=bode(G0*kc)Mag=20*log10(mag)Gm,Pm,Wcg,Wcp=margin(G0*kc)phi=

8、(dPm-Pm)*pi/180alpha=(1+sin(phi)/(1-sin(phi)Mn=-10*log10(alpha)Wcgn=spline(Mag,w,Mn)T=1/Wcgn/sqrt(alpha)Tz=alpha*TGc=tf(Tz,1,T,1)figure(1)bode(G0*kc,G0*kc*Gc)；F0=feedback(G0*kc,1)F=feedback(G0*kc*Gc,1)figure(2)step(F0,F)运行程序，得到结果：运行程序，得到结果：说明：运行以下语句，可以验证性能指标说明：运行以下语句，可以验证性能指标 margin(G0*kc*Gc)经比较看出：校

9、正后系统的快速性得到提高，经比较看出：校正后系统的快速性得到提高，超调量有所降低，大约由原来的超调量有所降低，大约由原来的40%40%下降到下降到29%29%。例例7-3 7-3 已知原系统开环传递函数为已知原系统开环传递函数为（1 1）绘制原系统的）绘制原系统的bodebode图，标出相角裕度和幅值图，标出相角裕度和幅值裕度；裕度；（2 2）现引入超前校正装置，）现引入超前校正装置，绘制校正后系统的绘制校正后系统的bodebode图，并与原系统的图，并与原系统的bodebode图进图进行比较；行比较；（3 3）绘制校正前后的阶跃响应曲线并进行比较。）绘制校正前后的阶跃响应曲线并进行比较。要求

10、：要求：MATLAB 程序如下：程序如下：G0=tf(100,0.04,1,0)Gm0,Pm0,Wcg0,Wcp0=margin(G0)Gc=tf(0.0262,1,0.0106,1)G=Gc*G0Gm,Pm,Wcg,Wcp=margin(G)bode(G0,G)T0=feedback(G0,1)T=feedback(G,1)figurestep(T0,T)运行程序，得到结果：运行程序，得到结果：Gm0=InfPm0=28.0243Wcg0=InfWcp0=46.9701Gm=InfPm=47.5917Wcg=InfWcp=60.3251%原系统原系统%校正后系统校正后系统从图中看出：从图中看

11、出：展宽了频带，展宽了频带，增加了相角裕增加了相角裕度，提高了系度，提高了系统的快速性。统的快速性。从图中看出：从图中看出：证明了不但证明了不但减小了系统减小了系统的超调量，的超调量，而且提高了而且提高了系统的快速系统的快速性。性。7.1.2 7.1.2 相位滞后校正相位滞后校正 相位滞后校正的等效相位滞后校正的等效RCRC网络如图所示网络如图所示 其网络传递函数可以写为其网络传递函数可以写为其中其中说明说明滞后校正环节具有极点滞后校正环节具有极点 ，零点，零点因此在因此在s s平面上，极点位于零点的右侧平面上，极点位于零点的右侧由于由于 ，例例7-4 7-4 设滞后校正环节的传递函数为设滞后

12、校正环节的传递函数为设：设：T=1T=1；=5=5，1010，1515，试绘制滞后校正环节的试绘制滞后校正环节的bodebode图和图和NyquistNyquist图图MATLAB 程序如下：程序如下：beta0=5T=1for i=1:3 beta(i)=i*beta0 G(i)=tf(T 1,beta(i)*T 1)endfigure(1)bode(G(1),G(2),G(3)figure(2)nyquist(G(1),G(2),G(3)运行程序，得到结果：运行程序，得到结果：说明：说明：最大相位滞后角所对应的频率最大相位滞后角所对应的频率在转折频率在转折频率 处，校正环节的幅值裕度衰处，

13、校正环节的幅值裕度衰减为减为 。滞后校正环节对高频噪声可以有效的抑制，因滞后校正环节对高频噪声可以有效的抑制，因此具有低通滤波器特性。此具有低通滤波器特性。例例7-5 7-5 已知原系统开环传递函数已知原系统开环传递函数要求：稳态误差要求：稳态误差 ,剪切频率剪切频率 ，相角裕度相角裕度 ，试用，试用MATLABMATLAB语句编制滞后校语句编制滞后校正装置的设计程序，绘制校正前后的正装置的设计程序，绘制校正前后的bodebode图和阶图和阶跃响应曲线。跃响应曲线。MATLAB 程序如下：程序如下：num=10den=1 5 0G0=tf(num,den)wcg=2kc=10mag,phase

14、,w=bode(G0*kc)magdb=20*log10(mag)Gr=-spline(w,magdb(1,:),wcg)alpha=10(Gr/20)T=10/(alpha*wcg)Gc=tf(alpha*T 1,T 1)F0=feedback(G0*kc,1)F=feedback(G0*kc*Gc,1)figure(1)bode(G0*kc,G0*kc*Gc)figure(2)step(F0,F)运行程序，得到结果：运行程序，得到结果：说明：说明：从从阶跃响应曲阶跃响应曲线看出，滞线看出，滞后校正使系后校正使系统响应速度统响应速度变慢，但平变慢，但平稳性得到改稳性得到改善，超调变善，超调变

15、小，振荡次小，振荡次数减少。数减少。校正前相校正前相角裕度为角裕度为校正后相校正后相角裕度为角裕度为例例7-6 7-6 已知原系统开环传递函数为已知原系统开环传递函数为若采用滞后校正装置若采用滞后校正装置 ，试绘制校正前后的试绘制校正前后的bodebode图和阶跃响应曲线，并图和阶跃响应曲线，并与超前校正结果进行比较。与超前校正结果进行比较。MATLAB 程序如下：程序如下：G0=tf(100,0.04,1,0)Gc=tf(0.5,1,2.5,1)G=Gc*G0Gm,Pm,Wcg,Wcp=margin(G)figure(1)bode(G0,G)T0=feedback(G0,1)T=feedba

16、ck(G,1)figure(2)step(T0,T)运行程序，得到结果：运行程序，得到结果：说明说明滞后校正后系滞后校正后系统的相角裕度统的相角裕度由原来由原来 的的增加到增加到 ，而剪切频率由而剪切频率由原来原来 的减少到的减少到 ，同时阶跃响，同时阶跃响应的平稳性得应的平稳性得到改善。到改善。MATLAB 程序如下：程序如下：超前校正结果进行比较超前校正结果进行比较G0=tf(100,0.04,1,0)Ga=tf(0.5,1,2.5,1)G1=Ga*G0Gb=tf(0.0262,1,0.0106,1)G2=Gb*G0T1=feedback(G1,1)T2=feedback(G2,1)ste

17、p(T1,r,T2,b,1.5)运行程序，得到结果：运行程序，得到结果：7 71 13 3 相位超前滞后校正相位超前滞后校正相位超前滞后校正装置的等相位超前滞后校正装置的等效效RCRC网络如图所示，其传递网络如图所示，其传递函数函数其中其中超前环节超前环节滞后环节滞后环节例例7-7 7-7 设设试绘制超前滞后校正环节的试绘制超前滞后校正环节的bodebode图和图和NyquistNyquist图图=10=10 beta=10T1=1T2=10num=conv(T1 1,T2 1)den=conv(T1/beta 1,T2*beta 1)Gc=tf(num,den)figure(1)bode(G

18、c)figure(2)nyquist(Gc)MATLAB 程序如下：程序如下：运行程序，得到结果：运行程序，得到结果：说明说明在在bodebode图幅频特性上，低频段与高频段均为图幅频特性上，低频段与高频段均为0dB0dB，这表明该环节的总增益为，这表明该环节的总增益为1 1；在相频特；在相频特性曲线上，曲线与性曲线上，曲线与 线交点处的频率刚好是线交点处的频率刚好是在幅频特性曲线的最低点，其值在幅频特性曲线的最低点，其值 ，当当 时，整个环节起滞后作用，时，整个环节起滞后作用，当当 时，起超前校正作用。时，起超前校正作用。例例7-87-8已知原系统开环传递函数已知原系统开环传递函数试设计超前

19、滞后校正装置试设计超前滞后校正装置 ，满足下列性能指，满足下列性能指标：速度误差系数标：速度误差系数 ，相角裕度，相角裕度 ，幅，幅值裕度值裕度要求：要求：（1）用）用MATLAB语句绘制校正前后系统语句绘制校正前后系统的的bode图和校正后的阶跃响应；图和校正后的阶跃响应；（2）用）用Simulink进行系统仿真，绘制仿进行系统仿真，绘制仿真结构图。真结构图。理论分析：理论分析：为了满足速度误差系数为了满足速度误差系数 ，可求得，可求得 ，得未校正系，得未校正系统的相角裕度为统的相角裕度为 ，而幅值裕度，而幅值裕度 因此系统是不稳定的因此系统是不稳定的 现采用滞后超前校正。为了保持足够的响应

20、速现采用滞后超前校正。为了保持足够的响应速度，取校正后的剪切频率度，取校正后的剪切频率 （不（不至于距离校正前的剪切太远），求得滞后校正至于距离校正前的剪切太远），求得滞后校正部分为：部分为：超前校正部分为：超前校正部分为：于是得滞后超前校正装置的传递函数为于是得滞后超前校正装置的传递函数为MATLAB 程序如下：程序如下：G0=tf(10,conv(1 0,conv(1 1,0.5 1)figure(1)margin(G0)Gc1=tf(1.43 1,0.143 1)Gc2=tf(7.14 1,71.4 1)G=Gc1*Gc2*G0figure(2)margin(G)T=feedback(G

21、,1)figure(3)step(T)运行程序，得到结果：运行程序，得到结果：说明：说明：校正后系统校正后系统的的bodebode图和图和阶跃响应曲阶跃响应曲线表明系统线表明系统具有良好的具有良好的动态特性，动态特性，比原来有了比原来有了明显的改善。明显的改善。说明：说明：校正满足设计要求校正满足设计要求用用SimulinkSimulink进行系统仿真，仿真结构图如图所示。进行系统仿真，仿真结构图如图所示。7.2 7.2 反馈校正反馈校正除了前面介绍的三种串联校正方法之外，反馈校除了前面介绍的三种串联校正方法之外，反馈校正（又称并联校正），也是广泛采用的系统设计正（又称并联校正），也是广泛采用

22、的系统设计方法之一。方法之一。设含有反馈校正的控制系统框图如图所示。设含有反馈校正的控制系统框图如图所示。待校正系统开环传递函数为待校正系统开环传递函数为 校正后系统开环传递函数为校正后系统开环传递函数为若在我们感兴趣的频段（即可接受校正的频段，一若在我们感兴趣的频段（即可接受校正的频段，一般在低、中频段）内般在低、中频段）内则有则有或写为或写为 表明表明，在，在bodebode图上只要待校正的开环幅频特性与期图上只要待校正的开环幅频特性与期望开环幅频特性相减，即可近似地获得望开环幅频特性相减，即可近似地获得 ，由于由于 已知，已知，可以立刻求出。可以立刻求出。用期望特性法设计反馈校正装置的步

23、骤如下：用期望特性法设计反馈校正装置的步骤如下：1）根据稳态指标绘制未校正系统的开环对数幅）根据稳态指标绘制未校正系统的开环对数幅频特性，即频特性，即2）根据给定性能指标绘制期望开环幅频特性，即）根据给定性能指标绘制期望开环幅频特性，即3）将以上两式相减，即可求得的对数频率特性，即）将以上两式相减，即可求得的对数频率特性，即（要求（要求 ）4）检查局部反馈回路的稳定性，并检查在期望）检查局部反馈回路的稳定性，并检查在期望剪切频率剪切频率 附近附近 的程度是否符合近似条件的程度是否符合近似条件5）由）由 求出求出6）检验校正后的系统是否满足性能指标要求。）检验校正后的系统是否满足性能指标要求。例

24、例7-97-9已知某随动系统如图所示，试设计反馈校正已知某随动系统如图所示，试设计反馈校正装置装置设计指标为：设计指标为：绘制校正前后的绘制校正前后的bodebode图及单位阶跃响应曲线，并作图及单位阶跃响应曲线，并作出仿真结构图。出仿真结构图。首先对本题做理论分析首先对本题做理论分析 1 1）由）由 可确定可确定k=100k=100，故固有部，故固有部分传递函数分传递函数为为由由BodeBode图可知系图可知系统的性能指标统的性能指标 2 2）求取期望频率特性。）求取期望频率特性。根据性能指标根据性能指标 ，可求得可求得 再由再由 ，可得，可得由经验公式由经验公式近似取近似取取中频宽取中频宽

25、 取取 取取 低频段的转折频率可由几何方法求得。低频段的转折频率可由几何方法求得。则校正后系统开环传递函数（即期望特性）则校正后系统开环传递函数（即期望特性）为为 3 3）确定校正装置）确定校正装置 的参数。的参数。将原系统化为结构图如图所示将原系统化为结构图如图所示由于低于由于低于 和高于和高于 的频率无需校正，故考虑的频率无需校正，故考虑 的校正频率段。由期望特性和原系统的频率特的校正频率段。由期望特性和原系统的频率特性可推得性可推得故反馈校正装置的传递函数为故反馈校正装置的传递函数为式中式中 是反馈校正装置的增益是反馈校正装置的增益G01=tf(100,conv(0.1 1,0.0067

26、 1)G02=tf(1,1 0)Hc=tf(0.0167 0,0.2 1)G0=G01*G02Ga=feedback(G01,Hc)G=Ga*G02figure(1)bode(G0,G)T0=feedback(G0,1)T=feedback(G,1)figure(2)step(T0,T)MATLAB 程序如下：程序如下：运行程序，得到结果：运行程序，得到结果：仿真结构图仿真结构图7.3 PID7.3 PID控制器设计控制器设计PIDPID控制是最早发展起来方法之一。由于其结构控制是最早发展起来方法之一。由于其结构简单，应用中参数整定方便，因此在工业控制简单，应用中参数整定方便，因此在工业控制中

27、得到广泛的应用。事实上，当今应用的工业中得到广泛的应用。事实上，当今应用的工业控制器中，有半数以上是采用控制器中，有半数以上是采用PIDPID或变形或变形PIDPID控控制方案的，其中包括传统的模拟式制方案的，其中包括传统的模拟式PIDPID控制和近控制和近年来微处理器实现的数字年来微处理器实现的数字PIDPID控制器。控制器。PIDPID控制器的数学表达式为控制器的数学表达式为 其传递函数为其传递函数为7.3.1 PID7.3.1 PID控制器的控制特性控制器的控制特性PIDPID控制有多种应用形式，如控制有多种应用形式，如P P、PIPI、PIDPID等，下面等，下面通过具体实例说明比例、

28、积分、微分各环节的控制通过具体实例说明比例、积分、微分各环节的控制作用。作用。例例7-107-10考虑一个三阶对象模型，考虑一个三阶对象模型，研究分别采用研究分别采用P P、PIPI、PIDPID控制策略下闭环系统的控制策略下闭环系统的阶跃响应。阶跃响应。G=tf(1,1 3 3 1)P=0.2:0.6:2.0for i=1:length(P)G_c=feedback(P(i)*G,1)step(G_c)hold onendMATLAB 程序如下程序如下（采用（采用P控制）控制）运行程序，得到结果：运行程序，得到结果：说明说明当当 值增大值增大时，系统响时，系统响应速度加快，应速度加快，幅值增

29、高幅值增高 说明说明根轨迹分根轨迹分析可知，析可知，当当 系统将不系统将不稳定稳定。G=tf(1,1 3 3 1)Kp=1Ti=0.7:0.2:1.5for i=1:length(Ti)Gc=tf(Kp*1,1/Ti(i),1,0)G_c=feedback(G*Gc,1)step(G_c)hold onendaxis(0,18,0,2)MATLAB 程序如下程序如下（采用（采用PI控制）控制）运行程序，得到结果：运行程序，得到结果：说明说明PIPI的控制作用时的控制作用时可以消除静差，可以消除静差，当当 值增大时，值增大时，系统超调变小，系统超调变小，响应速度变慢；响应速度变慢；若若 变小，则

30、变小，则超调增大，响应超调增大，响应加快。加快。G=tf(1,1 3 3 1)Kp=1Ti=1Td=0.1:0.4:2.1for i=1:length(Td)Gc=tf(Kp*Ti*Td(i),Ti,1/Ti,1,0)G_c=feedback(G*Gc,1)step(G_c)hold onendaxis(0,20,0,1.6)MATLAB 程序如下程序如下（采用（采用PID控制）控制）运行程序，得到结果：运行程序，得到结果：说明说明当当 增大时，增大时，系统的响应速度系统的响应速度加快，响应峰值加快，响应峰值提高提高 在实际应用时，由于纯微分环节无法实现，通在实际应用时，由于纯微分环节无法实现

31、，通常用带有滞后的近似一阶环节来代替。相应的常用带有滞后的近似一阶环节来代替。相应的PIDPID控制器传递函数为控制器传递函数为为了较好地近似，为了较好地近似，N N取值的可以大些。理论上取值的可以大些。理论上当当 时，近似微分环节将趋于理想微分时，近似微分环节将趋于理想微分环节。环节。例例7-11 7-11 考虑一个三阶对象模型，考虑一个三阶对象模型，采用采用PIDPID控制方法，令控制方法，令 ，研究近，研究近似微分环节对闭环系统阶跃响应的影响。似微分环节对闭环系统阶跃响应的影响。G=tf(1,1 3 3 1)Td=1;Kp=1;Ti=1N=100,1000,1:10Gc=tf(Kp*Ti

32、*Td,Ti,1/Ti,1,0)G_c=feedback(G*Gc,1)step(G_c)hold onfor i=1:length(N)nn=Kp*(Ti*Td,0,0+conv(Ti,1,Td/N(i),1)/Ti dd=Td/N(i),1,0 Gc=tf(nn,dd)G_c=feedback(G*Gc,1)step(G_c)hold onendfigurey,t=step(G_c)err=1-yplot(t,err)MATLAB 程序如下程序如下运行程序，得到结果：运行程序，得到结果：说明说明当选择当选择N=10N=10时，时，近似精度近似精度是令人满是令人满意的意的 说明说明N=10N

33、=10时误差时误差信号的变化信号的变化情况情况 7.3.2 PID7.3.2 PID控制器的参数整定控制器的参数整定实现实现PIDPID控制器的核心问题是根据给定的受控对控制器的核心问题是根据给定的受控对象，合理的选择控制参数。象，合理的选择控制参数。PIDPID控制器的三个参控制器的三个参数选择，本质上是在三维空间的搜索问题。早在数选择，本质上是在三维空间的搜索问题。早在19421942年，齐格勒年，齐格勒-尼柯尔斯（尼柯尔斯（Ziegler-NicholsZiegler-Nichols）就在大量的实验基础上，提出了一种实用的参数就在大量的实验基础上，提出了一种实用的参数整定规则，简称整定规

34、则，简称Z-NZ-N规则规则。具有具有PIDPID控制器的闭环系统框图。控制器的闭环系统框图。由图可见，由图可见，PIDPID控制器是一种串联校正装置。当被控制器是一种串联校正装置。当被控对象的数学模型已知时，可以采用各种不同的设控对象的数学模型已知时，可以采用各种不同的设计方法确定控制器的参数，包括解析方法和阶跃曲计方法确定控制器的参数，包括解析方法和阶跃曲线等。但是被控对象模型无法精确获得，则不能用线等。但是被控对象模型无法精确获得，则不能用解析方法去设计控制器。解析方法去设计控制器。在这种情况下，只能借助在这种情况下，只能借助于实验的方法来整定控制器的参数。此时，于实验的方法来整定控制器

35、的参数。此时，Z-NZ-N整整定规则更显出它的实用价值。定规则更显出它的实用价值。Z-NZ-N整定规则有两种实施的方法，它们共同的目整定规则有两种实施的方法，它们共同的目标都是使系统阶跃响应的最大超调量不超过标都是使系统阶跃响应的最大超调量不超过25%25%。第一种方法（响应曲线法第一种方法（响应曲线法 ）在被控对象的输入端加一阶跃信号，然后测出输出在被控对象的输入端加一阶跃信号，然后测出输出的响应曲线。如果被控对象既无积分环节，又无共的响应曲线。如果被控对象既无积分环节，又无共轭复数极点存在，则阶跃响应曲线呈轭复数极点存在，则阶跃响应曲线呈S S形。形。该曲线的特征可以用测得的延迟时间该曲线

36、的特征可以用测得的延迟时间 和时间常数和时间常数 T来表征来表征其相应的数学模型可以用下面的传递函数近其相应的数学模型可以用下面的传递函数近似的描述：似的描述：（如果阶跃响应的曲线不是（如果阶跃响应的曲线不是S S形，则不能应用此方法）形，则不能应用此方法）根据实验测得的根据实验测得的、T T参数，再按表参数，再按表7-17-1，即可整定，即可整定PIDPID控制器的参数。控制器的参数。表表7-1 7-1 第一种第一种Z-NZ-N整定方法整定方法 第二种方法（临界比例度法）第二种方法（临界比例度法）对于图所示闭环系统，设对于图所示闭环系统，设 ,即只采用即只采用比例控制。令比例控制。令 从零逐

37、渐增大，直至系统阶跃响应从零逐渐增大，直至系统阶跃响应呈现持续的等幅振荡。记下此时输出曲线对应的临呈现持续的等幅振荡。记下此时输出曲线对应的临界增益值界增益值 和振荡周期和振荡周期 （若无论怎样变化（若无论怎样变化 值，值，系统都不会呈现持续振荡，则不能用此方法），再系统都不会呈现持续振荡，则不能用此方法），再按表按表7-27-2给出的经验公式确定给出的经验公式确定 、和和 。这种。这种方法又称为方法又称为临界比例度法临界比例度法。表中比例度。表中比例度 ，临，临界比例度界比例度表表7-2 7-2 第二种第二种Z-NZ-N整定方法（整定方法（临界比例度法临界比例度法）例例7-12 7-12 已

38、知被控对象传递函数为已知被控对象传递函数为试用试用Z-NZ-N两种整定方法确定两种整定方法确定PIDPID控制器的参数，并控制器的参数，并绘制阶跃响应曲线。绘制阶跃响应曲线。理论分析理论分析 根据实验测得的根据实验测得的、T T参数，按表参数，按表7-17-1，即可整，即可整定定PIDPID控制器的参数。控制器的参数。MATLAB 程序如下程序如下运行程序，得到结果：运行程序，得到结果：G=tf(1,conv(1 1,conv(2 1,conv(5 1,10 1)step(G)可得可得K=1，=4.05，T=26.9-4.05=22.85理论分析理论分析 根据实验测得的根据实验测得的 、参数，

39、按表参数，按表7-27-2，即可，即可整定整定PIDPID控制器的参数。控制器的参数。MATLAB 程序如下程序如下运行程序，得到结果：运行程序，得到结果：G1=tf(1,conv(1 1,conv(2 1,conv(5 1,10 1)kc=7.698G=feedback(G1*kc,1)step(G,40)第二中方法基本参数第二中方法基本参数MATLAB 程序如下程序如下t1=0:0.5:150k=1;t=22.85;tor=4.05G=tf(1,conv(1 1,conv(2 1,conv(5 1,10 1)Gc,kc,ti,td=z_n1(k,t,tor,pid)Gcl1=feedbac

40、k(G*Gc,1)step(Gcl1,150)hold onKc=7.698;Tc=19.7Gc,kc,ti,td=z_n2(Kc,Tc,pid)Gcl2=feedback(G*Gc,1)step(Gcl2,150)z_n1自定义函数程序清单自定义函数程序清单functionGc,kc,ti,td=z_n1(k,t,tor,varargin)controllertype=varargin1switch controllertype case p delta=k*tor/t;Gc=tf(1/delta,1)kc=1/delta;ti=inf;td=0 case pi delta=1.1*k*to

41、r/t;ti=3.3*tor;kc=1/delta;Gc=tf(ti 1,ti 0)*kc;td=0 case pid delta=0.85*k*tor/t;ti=2*tor;td=0.5*tor;kc=1/delta;Gc1=tf(ti 1,ti 0)*kc;Gc2=tf(td 0,1)*kc Gc=Gc1+Gc2endz_n2自定义函数程序清单自定义函数程序清单functionGc,kc,ti,td=z_n2(Kc,Tc,varargin)controllertype=varargin1switch controllertype case p kc=Kc/2;Gc=tf(kc,1)ti=i

42、nf;td=0 case pi kc=Kc/2.2;ti=0.833*Tc;Gc=tf(ti 1,ti 0)*kc;td=0 case pid kc=Kc/1.7;ti=0.5*Tc;td=0.125*Tc Gc1=tf(ti 1,ti 0)*kc;Gc2=tf(td 0,1)*kc Gc=Gc1+Gc2end运行程序，得到结果：运行程序，得到结果：例例7-13 7-13 已知被控对象传递函数为已知被控对象传递函数为(1)(1)试用试用MATLABMATLAB语句输入语句输入 ，并绘制阶跃响应，并绘制阶跃响应曲线，估计模型参数曲线，估计模型参数K K、T值。值。(2)(2)根据模型参数，用根据

43、模型参数，用Z-NZ-N第一种方法设计第一种方法设计P P、PIPI、PIDPID控制器，并绘制阶跃响应曲线。控制器，并绘制阶跃响应曲线。理论分析理论分析 MATLAB 程序如下程序如下运行程序，得到结果：运行程序，得到结果：K=0.42G=tf(10,conv(1 1,conv(1 2,conv(1 3,1 4)step(G)k=dcgain(G)=0.61T=3.11-0.61=2.5可知可知MATLAB 程序如下程序如下k=0.42;l=0.61;t=3.11-lG=tf(10,conv(1 1,conv(1 2,conv(1 3,1 4)Gc1,Kp1=ziegler(1,k,l,t,

44、10);Kp1Gc2,Kp2,Ti2=ziegler(2,k,l,t,10);Kp2,Ti2Gc3,Kp3,Ti3,Td3=ziegler(3,k,l,t,10);Kp3,Ti3,Td3G_c1=feedback(G*Gc1,1);step(G_c1);hold onG_c2=feedback(G*Gc2,1);step(G_c2)G_c3=feedback(G*Gc3,1);step(G_c3)function Gc,Kp,Ti,Td=ziegler(key,vars)Ti=;Td=;k=vars(1);l=vars(2);T=vars(3);N=vars(4);a=k*l/T;if key

45、=1,Kp=1/a;elseif key=2,Kp=0.9/a;Ti=3.33*1;elseif key=3 Kp=1.2/a;Ti=2*1;Td=1/2;endswitch key case 1 Gc=Kp;case 2 Gc=tf(Kp*Ti,1,Ti,0);case 3 nn=Kp*Ti*Td*(N+1)/N,Kp*(Ti+Td/N),Kp;dd=Ti*Td/N,1,0;Gc=tf(nn,dd);endZiegler自定义函数程序清单自定义函数程序清单运行程序，得到结果：运行程序，得到结果：说明说明在在P控制下，控制下，系统存在稳态系统存在稳态误差。可以看误差。可以看出出PID的响应的响

46、应速度一般高于速度一般高于PI控制器控制器关于关于Z-NZ-N参数整定方法的几点说明：参数整定方法的几点说明：对象数学模型已知时，对象数学模型已知时，Z-NZ-N规则当然可以使用，但同时还规则当然可以使用，但同时还有许多其它的解析方法有许多其它的解析方法 和图解方法。实践证明，在对象和图解方法。实践证明，在对象模型不能精确确定的情况下（特别是在工业过程控制中），模型不能精确确定的情况下（特别是在工业过程控制中），Z-NZ-N规则的使用价值更加明显。规则的使用价值更加明显。注意注意Z-NZ-N规则的使用条件：第一种方法要求对象的阶跃响应规则的使用条件：第一种方法要求对象的阶跃响应曲线呈曲线呈S

47、S 形；第二种方法要求当以增益为参量变化时，闭环形；第二种方法要求当以增益为参量变化时，闭环系统的输出响应必须呈现等幅振荡状态，否则不适用。系统的输出响应必须呈现等幅振荡状态，否则不适用。使用使用Z-NZ-N整定方法，闭环系统阶跃响应呈现出整定方法，闭环系统阶跃响应呈现出10%60%10%60%的超调的超调量，经过量，经过Z-NZ-N整定，最大超调量的平均值约为整定，最大超调量的平均值约为25%25%。实际上，。实际上，Z-NZ-N整定方法对参数值提供了一种科学的估测，并为进一步整定方法对参数值提供了一种科学的估测，并为进一步精调提供了起点。精调提供了起点。PID PID控制器由于其算法简单，

48、鲁棒性好以及可靠控制器由于其算法简单，鲁棒性好以及可靠性高，被广泛应用于过程控制和运动控制。性高，被广泛应用于过程控制和运动控制。随着计算机技术与智能理论的发展，近年来出现随着计算机技术与智能理论的发展，近年来出现了许多智能了许多智能PIDPID控制器的设计和各种参数优化方法，控制器的设计和各种参数优化方法，为复杂动态不确定性系统的控制提供了新的途径，出为复杂动态不确定性系统的控制提供了新的途径，出现了现了神经网络神经网络PIDPID控制、遗传算法控制、遗传算法PIDPID参数优化、非线参数优化、非线性系统性系统PIDPID鲁棒控制、专家系统和模糊鲁棒控制、专家系统和模糊PIDPID参数整定方

49、参数整定方法法等，使传统的等，使传统的PIDPID控制方法与现代各种控制技术融控制方法与现代各种控制技术融为一体，展现出广阔的发展前景。为一体，展现出广阔的发展前景。MATLABMATLAB辅助工具的作用愈显突出。辅助工具的作用愈显突出。7.3.3 PID7.3.3 PID控制器参数的智能整定（参数优化方法）控制器参数的智能整定（参数优化方法）经典控制理论经典控制理论通常采用输出反馈，用频率法通常采用输出反馈，用频率法或根轨迹法设计控制系统，使闭环极点在期望的或根轨迹法设计控制系统，使闭环极点在期望的位置上。位置上。现代控制理论现代控制理论则更多地采用状态和极点任意则更多地采用状态和极点任意配

50、置的方法进行综合。由于状态反馈可以提供更配置的方法进行综合。由于状态反馈可以提供更丰富的状态信息和选择的自由度，因此使系统容丰富的状态信息和选择的自由度，因此使系统容易获得更优良的性能。易获得更优良的性能。7.4 7.4 状态反馈与极点配置状态反馈与极点配置设受控系统设受控系统 的状态空间表达式为的状态空间表达式为7.4.1 7.4.1 状态反馈状态反馈若用状态向量的线性反馈构成闭环系统，则称若用状态向量的线性反馈构成闭环系统，则称其为状态反馈系统其为状态反馈系统状态反馈系统的结构图如图所示状态反馈系统的结构图如图所示由图可见，状态反馈控制规律为由图可见，状态反馈控制规律为将此式代入状态方程，