第四章 计算机控制系统的控制算法 (第八讲)

1、计算机控制系统计算机控制系统电气与自动化工程学院电气与自动化工程学院自动化系计算机控制系统计算机控制系统1计算机控制系统计算机控制系统电气与自动化工程学院电气与自动化工程学院2第四章第四章 计算机控制系统的计算机控制系统的 控制算法控制算法第八讲第八讲- -第十二讲第十二讲第八讲第八讲计算机控制系统计算机控制系统电气与自动化工程学院电气与自动化工程学院3本章知识结构本章知识结构计算机控制系统计算机控制系统电气与自动化工程学院电气与自动化工程学院 计算机控制系统的设计，是指在给定系统性能计算机控制系统的设计，是指在给定系统性能指标的条件下，设计指标的条件下，设计控制器的控制规律和相应的数控制器的

2、控制规律和相应的数字控制算法字控制算法。 概概 述述数字控制器的设计方法数字控制器的设计方法按其设计特点分为三大类：按其设计特点分为三大类：4计算机控制系统计算机控制系统电气与自动化工程学院电气与自动化工程学院1. 模拟化设计方法模拟化设计方法(间接法间接法) 先设计连续系统的校正环节的传递函数先设计连续系统的校正环节的传递函数D(s)，然后将其离散化为然后将其离散化为D(z) ，进而设计数字控制器。，进而设计数字控制器。 典型代表为典型代表为PID算法算法。 (结构简单，各环节并联互相独立，调参方便结构简单，各环节并联互相独立，调参方便)5计算机控制系统计算机控制系统电气与自动化工程学院电气

3、与自动化工程学院2. 离散化设计方法离散化设计方法(直接法直接法) 已知被控对象的传递函数或特性已知被控对象的传递函数或特性G(z)，根据，根据所要求的性能指标，直接设计数字控制器所要求的性能指标，直接设计数字控制器D(z)。 常用方法有常用方法有最少拍有纹波设计最少拍有纹波设计、最少拍无最少拍无纹波设计纹波设计和和达林算法达林算法。3. 状态空间设计法状态空间设计法 能处理能处理多输入多输入-多输出多输出(MIMO)系统系统6计算机控制系统计算机控制系统电气与自动化工程学院电气与自动化工程学院 模拟化设计模拟化设计：适于采样周期短、控制算法简：适于采样周期短、控制算法简单的系统。忽略零阶保持

4、器和采样器，求出系统单的系统。忽略零阶保持器和采样器，求出系统的连续控制器，以近似方式离散化为数字控制器。的连续控制器，以近似方式离散化为数字控制器。 离散化设计离散化设计：适于采样周期长或控制复杂的：适于采样周期长或控制复杂的系统。直接使用采样控制理论设计数字控制器。系统。直接使用采样控制理论设计数字控制器。 7计算机控制系统计算机控制系统电气与自动化工程学院电气与自动化工程学院 数字控制理论基础数字控制理论基础 1. 计算机只能接受和处理二进制代码计算机只能接受和处理二进制代码0和和1及其组合，及其组合，这些二进制数可以表示某一物理量的大小，称之为离散量这些二进制数可以表示某一物理量的大小

5、，称之为离散量或数字量。但实际系统中的被控量是在时间上连续的信号，或数字量。但实际系统中的被控量是在时间上连续的信号，一般称之为连续量或模拟量。模拟量进行离散化并转换成一般称之为连续量或模拟量。模拟量进行离散化并转换成数字量后，才能由计算机处理。数字量后，才能由计算机处理。 因此，计算机控制系统也可以称为数字控制系统、离因此，计算机控制系统也可以称为数字控制系统、离散控制系统或采样控制系统。散控制系统或采样控制系统。 模拟控制系统称为连续控制系统。模拟控制系统称为连续控制系统。8计算机控制系统计算机控制系统电气与自动化工程学院电气与自动化工程学院 2. 离散离散(数字数字)控制系统与连续控制系

6、统与连续(模拟模拟)控制系统的本质控制系统的本质区别在于：模拟系统中的给定量、反馈量和被控量都是连区别在于：模拟系统中的给定量、反馈量和被控量都是连续型的时间函数，而在离散系统中，通过计算机处理得给续型的时间函数，而在离散系统中，通过计算机处理得给定量、反馈量和被控量是在时间上离散的数字信号。定量、反馈量和被控量是在时间上离散的数字信号。 把计算机引入连续控制系统中作为控制器使用，便把计算机引入连续控制系统中作为控制器使用，便构成了计算机控制系统。构成了计算机控制系统。 由计算机构成的控制系统，在本质上是一个离散系统。由计算机构成的控制系统，在本质上是一个离散系统。9计算机控制系统计算机控制系

7、统电气与自动化工程学院电气与自动化工程学院 传递函数的定义传递函数的定义 1. 线性连续系统中传递函数的定义是：线性连续系统中传递函数的定义是：零初始条件零初始条件下下，一个，一个线性线性环节环节(系统系统)的输出量的拉氏变换与输入量的输出量的拉氏变换与输入量的拉氏变换之比。的拉氏变换之比。 2. 线性离散线性离散(数字数字)控制系统中脉冲传递函数的定义：控制系统中脉冲传递函数的定义：零初始条件下零初始条件下，一个，一个线性线性环节环节(系统系统)的输出脉冲序列的的输出脉冲序列的Z变换与输入脉冲序列的变换与输入脉冲序列的Z变换之比。变换之比。10计算机控制系统计算机控制系统电气与自动化工程学院

8、电气与自动化工程学院4.1 数字控制器的间接设计方法数字控制器的间接设计方法4.2 数字数字PID控制算法控制算法4.3 数字控制器的直接设计方法数字控制器的直接设计方法4.4 纯滞后控制技术纯滞后控制技术4.5 控制算法控制算法Matlab仿真举例仿真举例4.6 小小 结结11计算机控制系统计算机控制系统电气与自动化工程学院电气与自动化工程学院4.1 数字控制器的间接设计方法数字控制器的间接设计方法 数字控制器的间接设计技术，是立足于连续控制数字控制器的间接设计技术，是立足于连续控制系统控制器的设计，然后在计算机上进行数字模拟来系统控制器的设计，然后在计算机上进行数字模拟来实现的。实现的。

9、实质：实质：在采样周期足够短的情况下，把数字控制在采样周期足够短的情况下，把数字控制器器(A/D、计算机、计算机、D/A)看作一个整体，其输入和输看作一个整体，其输入和输出为模拟量，将其等效为连续传递函数。出为模拟量，将其等效为连续传递函数。 12计算机控制系统计算机控制系统电气与自动化工程学院电气与自动化工程学院图图4.1 计算机控制系统的典型结构图计算机控制系统的典型结构图 工作过程工作过程：e(t)经经A/D转换器后变成数字量转换器后变成数字量e(k)，计算，计算机机(数字调节器数字调节器)按照一定的控制规律对按照一定的控制规律对e(k)进行计算处进行计算处理，将所得结果理，将所得结果u

10、(k)经经D/A转换器转换得到连续的控制转换器转换得到连续的控制量量u(t)作用到被控对象上，实现对被控参数作用到被控对象上，实现对被控参数y(t)的调节。的调节。信号采样信号采样+A/D转换转换D/A转换转换+信号保持信号保持13计算机控制系统计算机控制系统电气与自动化工程学院电气与自动化工程学院 4.1.1 采样周期与模拟化设计采样周期与模拟化设计 香农采样定理是使数字控制器进行模拟化设香农采样定理是使数字控制器进行模拟化设计的最基本的前提条件计的最基本的前提条件。 零阶保持器是零阶保持器是由由D/A转换器的转换器的输出保持来实输出保持来实现的。现的。频率特性：频率特性：max2s14计算

11、机控制系统计算机控制系统电气与自动化工程学院电气与自动化工程学院 信号通过零阶保持器后存在信号通过零阶保持器后存在幅值衰减幅值衰减和和相位滞后相位滞后。但如果采样周期但如果采样周期T足够小，即采样频率足够高时，可以足够小，即采样频率足够高时，可以忽略这一影响。对于小的采样周期，用幂级数展开，用忽略这一影响。对于小的采样周期，用幂级数展开，用T/2的时间滞后环节来近似：的时间滞后环节来近似： 设相位裕量减少设相位裕量减少5-155-15度，则采样周期应选为：度，则采样周期应选为：151(0.15 0.5)cT计算机控制系统计算机控制系统电气与自动化工程学院电气与自动化工程学院间接设计方法得以实现

12、的重要依据是：间接设计方法得以实现的重要依据是： (1) 采样周期要满足香农采样定理；采样周期要满足香农采样定理； (2) 采样周期足够小，达到零阶保持器的相位采样周期足够小，达到零阶保持器的相位滞后可以忽略不计的程度。滞后可以忽略不计的程度。16计算机控制系统计算机控制系统电气与自动化工程学院电气与自动化工程学院4.1.2 模拟化设计步骤模拟化设计步骤1. 设计假想的模拟控制器设计假想的模拟控制器D(s) 2. 正确地选择采样周期正确地选择采样周期 3. 将模拟控制器将模拟控制器D(s)离散化为数字控制器离散化为数字控制器D(z)4. 求出与求出与D(z)对应的差分方程对应的差分方程5. 编

13、制计算机程序，实现计算机控制编制计算机程序，实现计算机控制17计算机控制系统计算机控制系统电气与自动化工程学院电气与自动化工程学院1. 设计假想的模拟控制器设计假想的模拟控制器 解决方案：首先将计算机与零阶保持器合在一解决方案：首先将计算机与零阶保持器合在一起看作一个模拟环节，输入为起看作一个模拟环节，输入为e(t)，输出为，输出为u(t)，等，等效传递函数为效传递函数为D(s)；然后根据自控原理中的连续系；然后根据自控原理中的连续系统的统的频域设计法频域设计法、根轨迹法根轨迹法等校正方法设计等校正方法设计D(s)。零阶保持零阶保持器器18计算机控制系统计算机控制系统电气与自动化工程学院电气与

14、自动化工程学院2. 正确地选择采样周期正确地选择采样周期 实际应用中工程技术人员通常综合考虑来选择实际应用中工程技术人员通常综合考虑来选择采样周期采样周期T：(1) 从调节品质上考虑，希望从调节品质上考虑，希望T短以减少系统纯滞后短以减少系统纯滞后的影响，提高控制精度。保证在的影响，提高控制精度。保证在95%的系统的过渡的系统的过渡过程内采样过程内采样615次。次。(2) 从快速性和抗干扰性方面考虑，希望从快速性和抗干扰性方面考虑，希望T短使系统短使系统快速跟随输入的变化。快速跟随输入的变化。191(0.15 0.5)cT计算机控制系统计算机控制系统电气与自动化工程学院电气与自动化工程学院(3

15、) 从计算机的工作量和回路成本考虑，希望从计算机的工作量和回路成本考虑，希望T长使多长使多回路的每回路都有足够的计算时间。当被控对象的纯回路的每回路都有足够的计算时间。当被控对象的纯滞后时间滞后时间 较大时，常选较大时，常选 。(4) 从计算精度方面考虑，采样周期不应过短，应大从计算精度方面考虑，采样周期不应过短，应大于执行器的调节时间。于执行器的调节时间。 此外，被控量不同，采样周期也不同。此外，被控量不同，采样周期也不同。T20计算机控制系统计算机控制系统电气与自动化工程学院电气与自动化工程学院被被 控控 量量采采 样样 周周 期期/s注注流量流量15优选优选1s压力压力310优选优选5s

16、液位液位68优选优选7s温度温度1520优选纯滞后时间优选纯滞后时间成分成分1520优选优选18s表表4-1 4-1 采样周期的经验数据表采样周期的经验数据表21计算机控制系统计算机控制系统电气与自动化工程学院电气与自动化工程学院3. 将模拟控制器将模拟控制器D(s)离散化为数字控制器离散化为数字控制器D(z) 将连续控制器将连续控制器D(s)离散化为数字控制器离散化为数字控制器D(z)的方的方法有很多，如法有很多，如双线性变换法双线性变换法、前向差分法前向差分法、后向差后向差分法分法、冲击响应不变法、零极点匹配法、零阶保持、冲击响应不变法、零极点匹配法、零阶保持法等法等 。 通过近似方法，把

17、连续控制器离散化为数字通过近似方法，把连续控制器离散化为数字控制器。控制器。 22计算机控制系统计算机控制系统电气与自动化工程学院电气与自动化工程学院4. 求出与求出与D(s)对应的差分方程对应的差分方程 为了用计算机实现数字控制器，必须求出相应为了用计算机实现数字控制器，必须求出相应的差分方程，实现的方法有两种：的差分方程，实现的方法有两种：l 由数字控制器写出系统的由数字控制器写出系统的微分方程微分方程，然后进行，然后进行差差分处理分处理得到相应的差分方程，如得到相应的差分方程，如4.2节的数字节的数字PID控控制算法就是由此推导出来的。制算法就是由此推导出来的。l 根据数字控制器用根据数

18、字控制器用直接程序设计法直接程序设计法、串联实现法串联实现法等将其变为差分方程。等将其变为差分方程。23计算机控制系统计算机控制系统电气与自动化工程学院电气与自动化工程学院数字控制器数字控制器D(z)的一般形式：的一般形式：m个零点和个零点和n个极点，且个极点，且n=m，可写为：，可写为：对应的差分方程为：对应的差分方程为：上式即为直接程序设计法将上式即为直接程序设计法将D(z)转为差分方程的通式。转为差分方程的通式。120121212( )1mnnnbb zb zb zD za za za z121212012( ) () ( ) () ( )nmnnU zaza za zU zbbzbzb

19、zE z 为何为何n=m?24计算机控制系统计算机控制系统电气与自动化工程学院电气与自动化工程学院5. 根据差分方程编制相应程序，以实现计算机根据差分方程编制相应程序，以实现计算机控制控制 控制器控制器D(z)设计完成并求出控制算法后，按照图设计完成并求出控制算法后，按照图4-2所示的计算机控制系统检验其闭环特性是否符合设计所示的计算机控制系统检验其闭环特性是否符合设计要求，这一步可由计算机控制系统的数字仿真计算来要求，这一步可由计算机控制系统的数字仿真计算来验证，若满足要求则设计结束，否则应修改设计。验证，若满足要求则设计结束，否则应修改设计。25计算机控制系统计算机控制系统电气与自动化工程

20、学院电气与自动化工程学院1) 双线性变换法双线性变换法(梯形法或图斯汀法梯形法或图斯汀法)推导方法推导方法1：利用泰勒级数展开利用泰勒级数展开 得：得：则则D(s)离散化后的脉冲传递函数为：离散化后的脉冲传递函数为：2211e22e11e22sTsTsTsTsTzsTsT211zsTz211( )( )zsT zD zD s26计算机控制系统计算机控制系统电气与自动化工程学院电气与自动化工程学院推导方法推导方法2：梯形法数值积分梯形法数值积分 积分控制器积分控制器 用梯形法求积分运算用梯形法求积分运算 两边求两边求Z变换变换 则则D(s)离散化后的脉冲传递函数为：离散化后的脉冲传递函数为：0(

21、 )( )tu te t dt( )1( )( )U sD sE ss( )(1) ( )(1)2Tu ku ke ke k11( )( ) ( )( )2TU zz U zE zz E z12111( )11( )( )21( )2 11zsT zU zTzD zD szE zzT z27计算机控制系统计算机控制系统电气与自动化工程学院电气与自动化工程学院 映射关系：映射关系： 双线性变换置换公式双线性变换置换公式 把把 代入，代入， 取模的平方有：取模的平方有：sj22222122122TTzTT211zsTz28计算机控制系统计算机控制系统电气与自动化工程学院电气与自动化工程学院结论结论

22、：(1) 双线性变换法将双线性变换法将s的左半平面映射到的左半平面映射到z平面平面内的单位圆内，因此如果内的单位圆内，因此如果D(s)稳定，则稳定，则D(z)稳定。稳定。 (2) 频率畸变：双线性变换的一对一映射，保频率畸变：双线性变换的一对一映射，保证了离散频率特性不产生频率混叠现象，但产生了证了离散频率特性不产生频率混叠现象，但产生了频率畸变。频率畸变。 (3) 主要用于低通环节的离散化，此时频率主要用于低通环节的离散化，此时频率畸变程度较小。畸变程度较小。2tan2DATT29计算机控制系统计算机控制系统电气与自动化工程学院电气与自动化工程学院双线性变换的特点：双线性变换的特点： (1)

23、 变换关系简单，应用方便。可用计算机算出变换关系简单，应用方便。可用计算机算出D(z)的系数。的系数。(2) Z变换是将变换是将S左半平面变换到左半平面变换到Z平面的单位圆内，平面的单位圆内，所以如果所以如果D(s)稳定，则稳定，则D(z)亦稳定。亦稳定。 (3) 双线性变换不会引起高频混迭现象。双线性变换不会引起高频混迭现象。 (4) 它不能保持它不能保持D(s)的脉冲响应和频率响应。的脉冲响应和频率响应。(5) 变换前后，稳态增益保持不变。变换前后，稳态增益保持不变。3001( )( )szD sD z计算机控制系统计算机控制系统电气与自动化工程学院电气与自动化工程学院2) 前向差分法前向

24、差分法推导方法推导方法1：泰勒级数展开泰勒级数展开 , 有有 得到：得到： 则则D(s)离散化后的脉冲传递函数为：离散化后的脉冲传递函数为：sTze11sTzesTsT 1zsT1( )( )zsTD zD s31计算机控制系统计算机控制系统电气与自动化工程学院电气与自动化工程学院推导方法推导方法2：用一阶前向差分近似代替微分用一阶前向差分近似代替微分 微分控制器微分控制器 用前向差分近似代替用前向差分近似代替 令令n=k+1，并对两边作，并对两边作Z变换有：变换有： 得得D(s)离散化后的脉冲传递函数为：离散化后的脉冲传递函数为：( )( )( )U sD ssE s(1)( )( )e k

25、e ku kT1( )1( )( )( )zsTU zzD zD sE zT( )( )de tu tdt11( )(1)( )(1)( )(1)( )( )( )Tu ke ke kTu ne ne nTz U zE zz E z32计算机控制系统计算机控制系统电气与自动化工程学院电气与自动化工程学院 映射关系：映射关系： 前向差分法前向差分法(泰勒级数泰勒级数)置换公式置换公式 把把 代入，代入， 取模的平方有：取模的平方有： 令令|z|=1，则对应到，则对应到s平面上是一个圆，有：平面上是一个圆，有： 即当即当D(s)的极点位于左半平面以的极点位于左半平面以(-1/T,0)为圆心，为圆心

26、，1/T为半径的圆内，为半径的圆内，D(z)才在单位圆内，稳定。才在单位圆内，稳定。 1zsTsj33计算机控制系统计算机控制系统电气与自动化工程学院电气与自动化工程学院结论：稳定的系统经前向差分法转换后可能不稳定。结论：稳定的系统经前向差分法转换后可能不稳定。34计算机控制系统计算机控制系统电气与自动化工程学院电气与自动化工程学院(1)( )( )e ke ku kT前向差分法前向差分法(一阶前向差分一阶前向差分)置换公式置换公式 用用D(z)计算控制量计算控制量U(z)时，需要在时，需要在k时刻知道时刻知道k+1时时刻的刻的e(k+1)，这在物理上难以实现。，这在物理上难以实现。35前向差

27、分的特点：前向差分的特点： 计算机控制系统计算机控制系统电气与自动化工程学院电气与自动化工程学院3) 后向差分法后向差分法推导方法推导方法1：泰勒级数展开泰勒级数展开 , 有有 得到：得到： 则则D(s)离散化后的脉冲传递函数为：离散化后的脉冲传递函数为：sTeezTsTs1111zsTz1( )( )zsTzD zD ssTze36计算机控制系统计算机控制系统电气与自动化工程学院电气与自动化工程学院推导方法推导方法2：用一阶后向差分近似代替微分用一阶后向差分近似代替微分 用后向差分近似代替用后向差分近似代替 对两边作对两边作Z变换得离散化后的脉冲传递函数为：变换得离散化后的脉冲传递函数为：( )(1)( )e ke ku kT1( )1( )( )( )zsTzU zzD zD sE zTz37计算机控制系统计算机控制系统电气与自动化工程学院电气与自动化工程学院映射关系：映射关系： 后向差分法后向差分法(泰勒级数泰勒级数)置换公式置换公式 有：有：把把 代入，代入， 取模的平方有：取模的平方有：则：则： =0（s平面虚轴），平面虚轴）， 0（s右半平面），右半平面）， 1zsTzsj2222211 (1)()24 (1)()TTzTT38计算机控制系统计算机控制系统电气与自动化工程学院电气与自动化工程学院结论结论: (1)