计算机控制系统设计方法.ppt

《计算机控制系统设计方法.ppt》由会员分享，可在线阅读，更多相关《计算机控制系统设计方法.ppt（65页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、第三章第三章 计算机控制系统设计方法第三章 1-1计算机控制系统属于数字控制系统。所以，计算机控制系统控制器的设计属于数字控制器的设计。关于数字控制器的设计，有直接法和模拟法两种，这两种方法各有其特色，本章将分别给予介绍。第一节第一节连续域连续域离散化设计离散化设计一、设计的基本原理和步骤一、设计的基本原理和步骤连续域离散化设计分以下五个步骤完成：第1步：根据系统的性能和要求，选择采样频率。第2步：教材3-2所示，由于保持器会引入延迟，根据系统预定的性能指标，采用连续域的设计方法，设计出数字控制器的等效传递函数。保持器常采用零阶保持器，其一阶和二阶近似式表示如下：第三章 1-2（3-1）（3-

2、2）图3-2 连续域内等效设计框图第3步：选择适当的离散化方法，将D(s)离散化获得性能尽量等效的脉冲传递函数D(z)。第4步：针对由D(z)构成的离散闭环控制系统，检验其闭环性能。如图3-3所示。第5步：将D(z)编制成数字算法，在计算机上编程实现。二、前向差分法二、前向差分法 设模拟控制器传递函数为设模拟控制器传递函数为(3-3)转换成微分方程为(3-4)以一阶前向差分近似该微分，并代入（3-4）式，得(3-5)令k+1=n,上式的z变换为或（3-6）比较式（3-6）和式（3-3）可见，连续传递函数中的s在离散传递函数中的置换公式为推而广之，即给定模拟控制器传递函数D(s)，其等效离散传递

3、函数D(z)为：(3-7)下面讨论s平面和z平面之间的映射关系。因为平面上的虚轴（轴）是稳定与不稳定区域的分界线，所以应着重研究 轴在z平面内的映象。由 得：Z=1+TSZ=1+TS令 代入后可见，s平面上 轴映射在z平面上将右移1个单位，所以，采用前向差分法离散化，D(s)稳定，D(z)不一定稳定。前向差分法的特点总结如下：1、直接代换，具有串联性，变换方便；2、整个s左半平面映射到z平面z=1以左的区域，故D(s)与D(z)不具有相同的稳定性；3、因为D(s)|s=0=D(z)|z=1，故稳态增益维持不变；4、当采样周期T较小时，等效精度较好。三、后向差分法三、后向差分法设模拟控制器传递函

4、数为图3-5 前向差分法s平面稳定域在z平面内的映象(3-8)转换成微分方程(3-9)以一阶后向差分近似微分，得(3-10)代入式（3-9）得 对上式进行z变换，经整理为(3-11)比较式（3-11）与式（3-8），得s和z的置换公式为(3-12)推广而言，后向差分的离散化公式为(3-13)因为则当 时，可得 四、双线性变换法（四、双线性变换法（Tustin变换法）变换法）1、离散化公式图3-7中曲线r(t)以下的积分面积可采用个梯形面积之和来近似表示(3-15)其中前个梯形面积之和表示为C(k-1)。对式（3-15）两端求z变换，经整理，可得显然，s平面的虚轴以左映射为z平面单位圆之内的一个

5、小圆，如图3-6所示。所以，稳定的D(s)对应的D(z)也必定稳定。图3-6 后向差分法s平面稳定域 在z平面内的映象图3-7 梯形规则数值积分对式（3-15）两端求z变换，经整理，可得对比模拟控制器传递函数显然此时的转换关系为（3-16）由此可得双线性变换（梯形积分规则或Tustin变换）的离散化公式为（3-17）2、映射由式（3-16）可得 令s=j，可得|z|=1，相位随而变化，此即z平面的单位圆，如图3-8所示。可见，若D(s)稳定，D(z)也一定稳定。图3-8 双线性变换法S平面稳定域 在z平面内的映象 3、频率畸变与预修正令s=j，z=ejdt。这里，表示s域的角频率，d表示z 域

6、的角频率。根据置换公式，可得则即 （3-18）或 式（3-18）表明，s域的角频率与z域的角频率d是一个非线性关系，双线性变换产生了频率畸变。当 时，有（3-19）式（3-19）表明，D(s)在频率处的幅值D(z)等于 在的幅值，如图3-9所示。图3-9 双线性变换的频率畸变 由图3-9可见，双线性变换将s域0频段均压缩到z域的有限区间0，在系统校正装置的设计过程中，当D(z)用取代D(s)时，如果要保证在某个特征频率1处离散后，则必需进行频率预修正。进行预修正的步骤如下：（1）在特征频率1（如转折频率，自然振荡频率等）处，计算s域预修正频率（2）将模拟控制器传递函数修正为（3）将 转换成D(

7、z)，令（4）根据稳态增益相等的原则，按 进行增益匹配。总结以上各步骤，也可一步写出预修正双线性变换公式如下：例3-1 已知连续传递函数T=1s时，要 求 对 其 进 行 双 线 性 变 换，且 在 自 然 振 荡 频 率n=1rad/s处离散前后的幅值相等。解：第1步：取自然振荡频率n为特征频率1，并计算s域预修正频率1第2步：修正原传递函数D(s)为由D(s)可知 第3步：代入 进行双线性变换第4步：确定的增益根据稳态增益相等的原则：可得 解得 k=1 所以 4、双线性变换的特点（1）s域左半平面映射到z平面的单位圆内（一一对应）；（2）变换具有串联性，即针对相互串联的连续环节，可分别对各

8、个环节作双线性变换，然后相乘；（3）D(z)和D(s)具有相同稳定性；（4）频率特性发生畸变；（5）变换后稳态增益不变。五、零极点匹配法五、零极点匹配法1、转换规则令 Z=esT，将D(s)离散化为D(z)（3-20）式(3-20)表明，在z=esT时，D(s)的零点和极点均一一对应地映射到z平面，若时，则D(z)的分子上应添加(z+1)n-m因子。为保证变换前后增益不变，需进行增益匹配。D(z)的增益一般可按下式匹配：例3-2 已知连续传递函数为试采用零极点匹配法离散化，设采样周期T=1s。解：先将分解为零极点形式因T=1s，按转换公式 根据公式（3-20），得 由 可得 所以 2、特点（1

9、）D(z)和D(s)有相同稳定性；（2）D(s)的零、极点均按照Z=esT的关系与平面的零、极点一一对应；（3）稳态增益匹配，一般按 关系匹配。六、几种离散化公式的比较六、几种离散化公式的比较 以脉冲传递函数对原连续传递函数的保真度作为衡量标准，表3-1对上述几种离散化方法的主要特点进行了比较。第二章 1-30表表3-1几种离散化方法的主要特性比较几种离散化方法的主要特性比较第二节第二节 最少拍数字控制系统的设计最少拍数字控制系统的设计 在离散系统中，通常把一个采样周期称作一拍。最少拍系统，也称为最小调整时间系统或最快响应系统。一、闭环脉冲传递函数的选择一、闭环脉冲传递函数的选择 图3-10是

10、一个典型的计算机反馈控制系统。为了讨论问题方便，首先假设被控对象是稳定的。系统的闭环脉冲传递函数为由式（3-21）可以导出数字控制器的脉冲传递函数为（3-21)（3-22)图3-10 典型的计算机控制系统从式（3-22）可以看出，G(z)是零阶保持器和被控对象所固有的，不能改变。现在只需确定满足系统性能指标要求(z)的，就可以求得满足要求的数字控制器的脉冲传递函数。下面就讨论怎样确定(z)和D(z)。首先，最少拍系统要求稳态偏差e()=0。由图3-10可知，偏差的z变换为（3-23）根据终值定理，系统的稳态偏差为（3-24）设典型输入信号的一般形式为（3-25）式中,A(z)为z-1的多项式，

11、它不包含(1-z-1)因子，当k=1,2,3时，分别对应单位阶跃、单位速度和单位加速度输入。将式（3-25）代入式（3-24）则有:由于A(z)不含有(1-z-1)因子，所以1-(z)中必定含有(1-z-1)的至少k次的因式，才能使e()=0，即（(3-26）式中，F(z)是z-1和n次多项式。另一方面要使误差尽快为零，则上式右端应该是z-1最少多项式，因此应使1-(z)中(1-z-1)与分母中的(1-z-1)完全相约，即p=k，于是(3-27)不难看出，(z)具有z-1的最高幂次为k+n，其中n为F(z)的z-1的最高幂次，则可写成(3-28)式（3-28）表明，闭环系统在单位脉冲作用下，其

12、输出响应将在k+n个采样周期后变为零，或者说，在典型输入作用下，系统将经过k+n个采样周期达到稳态并实现跟踪。当n=0，即F(z)=1时，系统可经过最短时间（k个采样周期）达到稳态。(3-29)二、一般对象的最少拍系统的设计下面我们讨论在被控对象的脉冲传递函数稳定且没有纯滞后的情况下，如何根据不同的典型输入R(z)，确定数字控制器的脉冲传递函数D(z)。1、单位阶跃输入当r(t)=1时，这时 ，可选择(3-30)则数字控制器的脉冲传递函数为 （3-31）系统输出为即输出序列为C(0)=0,C(1)=1,C(2)=1,其波形如图3-11所示，由此可知离散系统的调节时间为1拍。图3-11 单位阶跃

13、输入时的输出序列2、单位速度输入当r(t)=t时，这时A(z)=Tz-1,k=2，可选择(3-32)则数字控制器的脉冲传递函数为(3-33)系统输出为 系统输出为即 C(z)=0+z-1+z-2+输出序列为C(0)=0,C(1)=1,C(2)=1，其波形如图3-11所示，由此可知离散系统的调节时间为1拍。图3-11 单位阶跃输入时的输出序列 2、单位速度输入 当r(t)=t时，这时A(z)=Tz-1，k=2，可选择（3-32）则数字控制器的脉冲传递函数为 （3-33）系统输出为即 输出序列为 其输出序列如图3-12所示。从图3-12可以看出，单位速度输入时，离散系统的调节时间为2拍。图3-12

14、 单位速度输入时的输出序列3、单位加速度输入当 时，这时 ，可选择则数字控制器的脉冲传递函数为 （3-35）系统输出为（3-34）即 输出序列为 其波形如图3-13所示。由图3-13可能看出，单位加速度输入时，离散系统的调节时间为3拍。上面的讨论结果汇总于表3-2。表3-2 三种典型输入时的最少拍系统 例3-3 设单位速度反馈线性离散系统如图3-10所示，设被控对象的传递函数 ，采样周期T=0.1s，试设计单位速度输入时的最少拍系统的数字控制器D(z)。图3-12 单位速度输入时的输出序列 解:系统广义被控对象的脉冲传递函数为将T=0.1s代入上式，经过整理后得 因为是单位速度输入，所以选择则

15、数字控制器的传递函数D(z)为系统输出为 其输出序列的波形与图3-12相同，离散系统经过两个采样周期T以后（即k 2），C(kT)=r(kT)，这就是按单位速度输入设计的最少拍系统。三、特殊情况下最少拍系统的设计现在研究考虑零点和极点分布时，最少拍数字控制系统如何设计。根据图3-10，由式(3-23)得（3-36）（3-37）由式（3-37）可以得出（3-38）闭环系统的稳定性是由(z)的极点在z平面的分布决定的，(z)的零点对系统的快速性也将产生一定的影响。由式（3-38）可以看出，被控对象的零极点对(z)有直接影响，特别是那些单位圆上或单位圆外的零极点更应引起我们的注意。（1）G(z)含单

16、位圆上或单位圆外的极点 当被控对象G(z)中含有单位圆上或单位圆外的极点时，这些极点也将作为(z)的极点，此时系统是不稳定的。为了保证系统的稳定性，有两个办法可以设法消除这些极点。一是让包含与这些极点相对应的零点通过零极相消，去掉G(z)中的不稳定极点。二是将这些极点作为的零点。(2)含单位圆上或单位圆外的零点 当G(z)含有单位圆上或单位因外的零点时，这些零点也将作为(z)的零点，此时系统的调节时间将增加。若要不增加调节时间，必然要消除这些零点。(3)G(z)含有纯滞后环节 如果G(z)中含有纯滞后环节，则G(z)可以表示为(3-39)式中,为G（z）的不含滞后因子的部分。由式（3-22）及

17、式（3-39）可以得出(3-40)式中，z因子表示D(z)具有超前环节，即在上次输入信号到来之前就有对应的输出，这在物理上是无法实现的。如果让(z)的分子中含有z-1因子，从而抵消掉G(z)中的z-r因子，这样，数字控制器的D(z)中就不会出现超前因子了。通过以上讨论，可以得出如下一些选择闭环脉冲传递函数(z)的限制条件：数字控制器在物理上应是可实现的有理多项式，即 选择(z)时，应包含G(z)中的滞后因子，G(z)的单位圆上或单位圆外的零点应保留，并作为(z)的零点。G(z)的单位圆上或单位圆外的极点，应作为e(z)=1-(z)的零点。第三节第三节 最少拍无波纹数字控制系统设计最少拍无波纹数

18、字控制系统设计一、波纹产生的原因一、波纹产生的原因 设在图3-10中，G0(s)=10/s(s+1)，对于阶跃输入，可设计出最少拍有波纹数字控制器即从第一个采样节拍起，偏差巳达到并保持为零了。如果从第一个节拍起控制量u(k)也能保持恒定(常数或零)，则输出就可保持不变。但是可见控制量实际上并不恒定，而在波动。这样一个波动的控制量作用在保持器的输入端，保持器输出也必然会波动，系统输出中也必然会出现波纹。二、消除波纹的条件二、消除波纹的条件 从以上讨论可知，波纹是由于u(k)波动引起的，要消除波纹，必须使u(k)为一恒定值。下面来讨论u(k)为一恒定值的条件。由于所以（3-41）由上一节的讨论可知

19、，最少拍系统要求闭环系统的脉冲传递函数(z)=1-(1-z-1)kF(z)，F(z)为z-1的有限项多项式。这就是说，当(z)为z-1的有限项多项式时，输出C(k)也能在若干拍后达到稳定。同样，由式（3-41）可以看出，如果(z)/G(z)可以表示成z-1的有限项多项式时，输出c(k)也能在若干拍后达到稳定。所以，系统无波纹的条件是(z)与G(z)之比可以用z-1的有限项多项式表示。由于所以，为了使(z)与G(z)之比为z-1的有限多项式，应该让(z)包含G(z)的全部零点z=zi。则系统就会是无波纹的。例 3-5 设 离 散 数 字 控 制 系 统 如 图 3-10所 述，采 样 周 期 T

20、=1s，G0(s)=10/s(s+1)，试设计单位阶跃输入的最少拍无波纹数字控制器D(z)。解 采用零阶保持器时，广义被控对象脉冲传递函数为在有波纹时可选 考虑无波纹要求，应使(z)包含G(z)的全部零点。所以应选择的零点多了一个，必然会使 的阶次增加一阶，即式中(1-z-1)是由输入的形式决定的，(1+bz-1)是由于的零点增加而相应产生的。因为 所以有比较等式两边的系数有由此解得所以最少拍无波纹数字控制器的D（z）为现在检验一下是否符合无波纹的要求两拍以后，数字控制器的输出即为零，所以能保证系统无波纹。第四节第四节 大林算法大林算法 美国IBM公司的大林，在1968年提出了一种针对工业生产

21、过程控制中的含有纯滞后对象的控制算法，具有较好的效果。一、大林算法的基本形式一、大林算法的基本形式设被控对象具有纯滞后的阶或二阶惯性环节，它们的传递函数为式中，T1和T2表示被控对象的时间常数；为被控对象的纯滞后时间；K为放大系数。假定为T的整数倍，即=LT(L=1,2,3）。大林算法的目标是设计一个数字控制器，使整个闭环系统的传递函数(s)相当于一个带纯滞后的一阶惯性环节，而且要求期望的闭环系统的纯滞后时间与被控对象的纯滞后时间相等，即式中，T0为校正后闭环系统的时间常数。由于是在平面上讨论数字控制器的设计，因此假设离散系统中采用的是零阶保持器，采样周期为T，那么闭环控制系统的脉冲传递函数为

22、 上式就是按大林算法而导出的数字控制器D(z)，可见它由被控对象唯一地确定。例3-6 设 ，试用大林算法设计阶跃输入时的数字控制器D(z)。解 包括采样保持器在内的广大对象脉冲传递函数为现构成T0=2s的一阶闭环系统，纯滞后也为2s。当输入单位阶跃时，闭环输出的z变换为 可见，系统输出是以指数形式变化的。二、振铃现象及消除方法 在设计具有纯滞后被控对象的计算机控制系统时应当注意可能会出现振铃现象。继续研究例3-6。其控制量为 很明显，数字控制器的输出控制量u(k)以2T为周期上下大幅度摆动，大林把这种现象称为振铃现象。振铃现象会导致系统执行机构的磨损增加。观察U(z)的表达式，其中包含z=-0.718的极点，该极点会引起交替变号的衰减序列。大林提出了一个消除振铃的方法是在引起振铃的D(z)的极点因子中，令z=1，从而消除振铃。如上例中，只要将D(z)中1+0.718z-1因子改为了1+0.718即可。修改后的D(z)为