《第5章数字控制器离散直接设计方法课件.ppt》由会员分享，可在线阅读，更多相关《第5章数字控制器离散直接设计方法课件.ppt（67页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、第五章第五章 数字控制器离散直接设计方法数字控制器离散直接设计方法v第一节 脉冲传递函数解析设计原理v第二节 最少拍控制系统设计v第三节 纯滞后对象的控制算法史密斯预估器v 第四节 纯滞后对象的控制算法大林算法v第五节 数字控制器的程序实现 解析设计方法是一种直接在z 域中设计方法，其基本思想是依据给定的控制环的系统结构，由系统的指标要求及实现的约束条件确定期望闭环z传递函数，通过代数方法求出所设计控制器的z传递函数第一节第一节 脉冲传递函数解析设计原理图5.2 典型计算机控制系统结构图一、数字控制器一、数字控制器D(z)的一般形式的一般形式零阶保持器的传递函数为：sesHTs1)(广义被控对
2、象的脉冲传递函数为：)()()(0sGsHZzG求出开环系统的脉冲传递函数为)()()()()(zGzDzEzCzW闭环系统的脉冲传递函数为：( )( )( )( )( )1( )( )C zD z G zzR zD z G z误差的脉冲传递函数为：( )1( )( )1( ) ( )eE zzR zD z G z又)(1)(zze求出数字控制器的脉冲传递函数为：)(1)()()(zzGzzD二二 在物理上的可实现性在物理上的可实现性 所谓数字控制器 的物理可实现，即要求数字控制器算法中不允许出现对未来时刻的信息的要求，反映在z传递函数上，即 的无穷级数展开式不能出现z的正幂次项，表现为分子阶
3、次必然要低于或等于分母阶次。)(zD设被控对象具有 个采样周期纯滞后，纯滞后脉冲传递函数为 ，其脉冲传递函数为： 111101211(1)()1tKzb zG zzgg zga z 12012( )( )( )1( )() (1( )zzzD zG zzgg zg zLz)(zD设期望闭环传递函数为：123123( ) zzzzL -1-2-3123-1-2-1-201212+( )( )=( )1-( )+1-zzzzD zG zzgzzLzzgg式(5-10)代入(5-9)得： (5-9) (5-10)ttz将式(5-8)代入式(5-7)( )D zz若 物理可实现， 中不能包含超前因子
4、，所以应满足下式：12-1 = = =0因此，期望闭环 z的传递函数为：-t-+1-t-1+1+1( )=+ = +zzzL zzL（t）（）即期望闭环传递函数 与广义对象 具有同样的滞后。（5-12）( )D z)(z)(zG三三 闭环稳定性要求闭环稳定性要求)()()()(zRzGzzU为保证闭环系统稳定， 的零点应包括 所有不稳定的极点，而 应包括 的所有不稳定的零点。因此能够设 为下式jjiizpzbzqzazG)()1 ()()1 ()(11为使闭环系统稳定，可建立(5-15)式，并能推导出(5-16)式。iizFzaz)()1 ()(11jjzFzbz)()1 ()(121)(1z
5、(5-15)(5-16)(zG)(z四四 离散直接设计一般步骤离散直接设计一般步骤 离散化设计是把计算机控制系统近似看作离散系统，所用的数学工具是差分方程和z变换，完全采用离散控制系统理论进行分析，直接设计数字控制器。(1) 根据式(5-2)求出带零阶保持器的广义被控对象的脉冲传递函数 。(2) 根据系统的性能指标要求和其它约束条件，确定闭环系统的脉冲传递函数 。(3) 根据式(5-7)求出数字控制器的脉冲传递函数 。(4) 利用仿真软件，对求出的数字控制器 进行校验，若达到设计要求，进行下一个步骤，否则进行再设计。(5) 根据数字控制器的脉冲传递函数 ，求出差分方程，编写控制程序。(6) 接
6、入硬件，进行系统调试。（1）单位阶跃输入：（2）单位速度输入：（3）单位加速度输入：它们都可以表示为：(5-17)式中， 是不包括 的 多项式。 为正整数，对于不同的输入，只是 不同而已，一般只讨论 的情况。)(zA)1 (1 z1z3 , 2 , 1qqzzAzR)1 ()()(1qq111)(,1)(),( 1)(zzRssRttR2112)1 ()(,1)(,)(zTzzRssRttR3111232)1 (2)1 ()(,1)(,21)(zzzTzRssRttR自动控制系统中，有三种典型的输入形式，其表示形式为：第二节 最少拍控制系统设计 所谓最少拍控制，就是要求设计的数字调节器能使闭环
7、系统在典型输入作用下，具有最快的响应速度，能在有限采样周期内达到采样点上无稳态误差或无静态误差。最少拍系统是建立在时间最优控制的基础上，因此也叫最快响应系统一一 最少拍无差系统最少拍无差系统 最少拍无差控制器的设计任务就是根据式(5-7)求出数字控制器的脉冲传递函数 ，使闭环系统在特定的典型输入作用下，以最少拍结束响应过程，并在采样时刻系统不存在稳态误差，输出能够准确地跟踪输入。因此最少拍无差系统的闭环脉冲传递函数形式为：)(zDnnzazazaz2211)(（1）调节时间最短，即系统跟踪输入信号所需的采样周期数最少；（2）在采样点处无静差，即对特定的参考输入信号，在达到稳态后，系统在采样点能
8、精确实现对输入信号的跟踪；（3）设计出来的数字控制器必须是物理上可以实现的；（4）闭环系统必须是稳定的。对最小拍控制系统设计的要求是：一一 最少拍无差系统最少拍无差系统)(ze二二 最小拍闭环脉冲传递函数的确定最小拍闭环脉冲传递函数的确定 最小拍控制系统的设计要求是对特定的参考输入信号，在系统达到稳态后，系统在采样点处静差为零。根据此约束条件可以构造出系统的误差脉冲传递函数 。典型计算机控制系统结构图如图5.2所示。图5.2 典型计算机控制系统结构图由式（5-5）和式（5-6），得(5-19)利用Z变换的终值定理可以求出稳态误差为mezzkzzAzzzEzke)1 ()()()1 (lim)(
9、)1 (lim)(lim11111(5-20)由于 不包括 的因子，因此稳态误差为零的条件是 含有 ，则可为下列形式)(zA)1 (1z)(zemz)1 (1（5-21）式中 为 的有限多项式，即)(zF1znnzfzfzfzF22111)()()(1)()()(zRzzRzzEe)()1 ()(1)(1zFzzzpe)(zF)(zF 由最小拍控制系统的时间最短约束条件来确定 的形式。当取 1，不仅可以简化数字控制器，降低控制器阶数,而且还可以使 的项数最少，调节时间最短。由式（5-21）得 为)(zE)(zemezz)1 ()(1（5-23）mezzz)1 (1)(1)(1 对于三种典型输入
10、信号下，最小拍控制系统的 和 汇总于表5-1中。)(ze)(z三三 最小拍控制器的设计最小拍控制器的设计 的确定的确定00( )( )( ) ( )( )1 ( )CCTsCG zZ Hs GsZ Hs GseZGss 1、广义对象的脉冲传递函数广义对象的脉冲传递函数:数字控制器数字控制器零阶保持器零阶保持器被控对象被控对象)(zD2、系统的闭环脉冲传递函数为：( ) ( )( )1( ) ( )D z G zzD z G z ( )1( )1( )( )1( ) ( )eE zzzR zD z G z 3、误差E(z)的脉冲传递函数为： ( )( )( )( )( )( ) 1( )ezzD
11、 zG zzG zz 4、数字控制器的脉冲传递函数为： 时间序列脉冲传递函数 单位阶跃输入 单位速度输入单位加速度输入 通式 ()()R nTu nT 11( )1R zz ()R nTnT 112( )(1)TzR zz 21()()2R nTnT 21113(1)( )2(1)T zzR zz 1( )( )(1)mA zR zz 典型的输入形式:A(z)是不包括(1-z-1)因子的关于z-1的多项式 1( ) ( )( )( ) ( )(1)eemz A zE zz R zz 11111 ( )lim (1) ( )( ) lim (1)( )(1) 0zemzezE zA zzzZ 1
12、( )(1)( )MezzF zMm F(z)是不包含零点z=1的z-1的多项式根据z变换的终值定理，系统的稳态误差e()取 F(z) = 1，M = m讨论：不同的输入最少拍系统的调整时间？讨论：不同的输入最少拍系统的调整时间？则有：1( )(1)mezz 单位阶跃输入：11( )(1( )ezzzz 单位速度输入：1221( )(1( )2)ezzzzz 单位加速度输入：11233( ) (1) 33ezzzzzz 单位阶跃输入时：111( )( ) ( )(1)()11eE zz R zzz (0)1 ( )(2 )()0ee Te Te kT，系统经过系统经过T，系统稳定，系统稳定11
13、2112( )( ) ( )(1)(1)eTzE zz R zzTzz (0)0,( ),(2 )(3 )()0ee TTeTeTe kT系统经过系统经过2T，系统稳定，系统稳定单位速度输入时：单位加速度输入时：211221 3121 3(1)( )( ) ( )(1)2(1)22eT zzTTE zz R zzzzz 22(0)0( )(2 )(3 )(4 )()022TTee TeTeTe Te kT，系统经过系统经过3T，系统稳定，系统稳定1、对应于三种不同典型输入，系统分别经过、对应于三种不同典型输入，系统分别经过T，2T，3T 系统达到稳定，系统的稳态误差为系统达到稳定，系统的稳态误
14、差为0。2、对应于不同的典型输入，为了得到最少拍响应，应选、对应于不同的典型输入，为了得到最少拍响应，应选 择合适的择合适的e(z)。3、对应于典型输入，选定、对应于典型输入，选定e(z)后，可根据后，可根据G(z)得到得到D(z)。 1( )( )( )( )( )( )( )eeezzD zG zzG zz 结论表5.1 三种典型输入下的理想最少拍无差系统 例5.1 设最小拍控制系统如图5.2所示，被控对象的传递函数 ,采样周期 ，采用零阶保持器，试设计在单位速度输入时的最小拍控制器。 解：根据图5.2可求出系统广义被控对象脉冲传递函数 )15 . 0(2)(sssGs5 . 0T)1)(
15、1 ()1 ( )2112(2112 )2(4)2(4 )2(4)1( ) 15 . 0(21)(1211211222222zezzezzessseZsssZsseZssZsseZssseZzGTTTTsTsTsTs将 代入，有 根据题意，输入信号为单位速度输入，即 ，则有： 代入式（5-7）求出最小拍控制器为下面对设计出来的最小拍控制器进行分析与校验。系统闭环脉冲传递函数为当输入为单位速度信号时，系统输出序列的 变换为 s 5 . 0T)368. 01)(1 ()718. 01 (368. 0)(1111zzzzzGttr)(21)1 ()(zze)718. 01)(1 ()368. 01)
16、(5 . 01 (435. 5)(1111zzzzzD212)(zzz5432211215432 )1 ()2()()()(TzTzTzTzzTzzzzzRzYZ 即 输出响应如图5.3所示。从图中可以看出，当系统为单位速度输入时，经过两拍以后，输出量完全等于输入采样值，即 。但在各采样点之间还存在着一定的误差，即存在着一定的波纹。 图5.3 单位速度输入时最小拍控制系统输出响应曲线图,3)3(,2)2(, 0) 1 (, 0)0(TyTyyy)()(krkyTtY(kT) 2T3T4T5Tx(t) y(kT) 123450第三节第三节 纯滞后补偿控制算法纯滞后补偿控制算法-史密斯预估器史密斯
17、预估器 具有纯滞后特性的对象属于比较难于控制的一类对象，对其控制需采用特殊处理方法。一般来说，对于控制这类对象，快速性要求是次要的，调整时间允许较多的采样周期，而对稳定性、不产生超调的要求是主要的。 基于此，人们提出了许多设计方法，比较有代表性的方法有纯滞后补偿控制史密斯(Smith)预估和大林(Dahlin)算法。spesGsUsYsG)()()()(一、纯滞后对系统控制品质的影响一、纯滞后对系统控制品质的影响 常规控制系统的结构框图如图5.4所示。被控对象含有纯滞后特性，其传递函数为式中， 为被控对象不含纯滞后特性的传递函数。)(sGp图5.4 有纯滞后的常规反馈控制结构图一、纯滞后对系统
18、品质的影响一、纯滞后对系统品质的影响 系统的闭环传递函数（不考虑扰动时）为 （5-27） 系统的特征方程为 （5-28） 这是一个复变数 的超越方程，方程的根也就是系统闭环特征根，将受到纯滞后时间 的影响。通过对系统的频域分析可知， 的增加不利于闭环系统的稳定性，使闭环系统的控制品质下降。因此，在进行控制系统设计时，为了提高系统的控制品质，应设法努力减小处于闭环回路中的纯滞后。除了选择合适的被控变量来减小对象的纯滞后外，在控制方案上，也应该采用各种补偿的方法来减小或补偿纯滞后造成的不利影响。spspesGsDesGsDsRsYs)()(1)()()()()(0)()(1 spesGsDs二、史
19、密斯补偿控制原理二、史密斯补偿控制原理 Smith预估补偿控制是克服纯滞后的一个有效的控制方法，其思想是根据系统的动态特性建立一个模型加入到反馈控制系统中，使延迟 时间的被控量提前反映到控制器，从而减少超调量和加快控制过程。根据这个控制思想，控制器D(s)联接一个补偿环节，用来补偿被控对象中的纯滞后部分，这个补偿环节称为预估器，其传递函数为 ，由Smith预估器和控制器组成的补偿回路称为纯滞后补偿器，其传递函数为 ( )(1)spGse 由Smith预估器和控制器组成的补偿回路称为纯滞后补偿器，其传递函数为 。补偿后系统框图如图5.5所示，图5.6为其转换后的等效形式。 图5.5 Smith预
20、估器控制系统结构图)(sD 实际工程上设计Smith预估器时，将其并联在控制器 上，对图5.5作方框图等效变换，得到图5.6所示的形式。 图5.6 Smith预估器控制系统等效图 图中虚线部分是带纯滞后补偿控制的控制器，其传递函数为 (5-29) 经过纯滞后补偿控制后系统的闭环传递函数为 -( )( )( )=( )+ ( )( )(1-)spU sD sD sE sD s Gse1-( )( )( )( )+ ( )( )(1-)( )( )=( )( )( )+( )( )1+ ( )( )(1-)spssppssppspD s Gs eD s Gs eD s GseY ssD s Gs
21、eU sD s Gs eD s Gse111 (5-30) 由式(5-30)可见，带纯滞后补偿的闭环系统与图5.4所示的理想结构是一致的，其特征方程为： 。纯滞后环节 已经不出现在特征方程中，故不再影响闭环系统的稳定性。分子中的 并不影响系统输出量 的响应曲线和系统的其他性能指标，只是把控制过程推迟了时间 。换句话说，纯滞后补偿控制系统在单位阶跃输入时，输出量 的响应曲线和系统的其他性能指标与控制对象不含纯滞后特性时完全相同，只是在时间轴上滞后 ，闭环系统输出特性如图5.7所示。)()(1)()(s) sGsDesGsDpsp0)()(1sGsDpsese)(ty)(ty 图5.7 闭环系统输
22、出特性示意图0) (ty) (ty) (ty0t三、三、 史密斯补偿器的计算机实现史密斯补偿器的计算机实现 带有纯滞后Smith补偿器的计算机控制系统如图5.8所示。图5.8 纯滞后补偿计算机控制系统结构图 图中 为数字PID控制器；Smith补偿器 与对象特性有关； 为被控对象传递函数中不包含纯滞后环节的部分。)(zD)(sGp-( )=( )(1-)spD sGsesppsesGsTKesG)(1)(1)(sTKsGpp 下面以一阶惯性纯滞后对象为例，说明Smith纯滞后补偿器的计算机实现过程。设被控对象的传递函数为式中Smith补偿器为：（5-32）离散化处理为：（5-33）式中， ，
23、， （取整数）。pTTea1)1 (1pTTeKbTN1)1 ()(sTeKsDps-1-N1-11(1-)(1-)( )=Z=(1-)+11-sspb zeKeD szsT sa z 为了便于说明Smith补偿器的计算机实现过程，将图中5.8的虚框部分变换为图5.9所示形式。图5.9 Smith补偿器计算机实现结构图由图5.9有)()()()()()()(zUzPzPzYzUzYzD（5-34）为了便于计算机实现，由式（5-33），令)1 ()()(NzzPzY11111)()(zazbzUzP可得到Smith补偿器的差分方程为)()()() 1() 1()(11Nkpkpkykubkpak
24、p（5-35）)(Nkp由式（5-35）可见，Smith补偿器的差分方程中有 项。 那么如何用计算机产生该纯滞后信号，对纯滞后补偿控制的计算机实现是至关重要。 下面介绍一种在计算机控制系统中常用的产生纯滞后信号的方法，即存储单元法。 为了形成纯滞后 步的信号，需在内存中开辟 个存储单元，用来存储 的历史数据N1N)( kp 用上述方法产生纯滞后信号后，由式（5-35）即可求出 。Smith补偿控制算法的实现步骤为：)(ky（1）计算偏差（2）计算控制器输出)(2ke) 1() 1()(11kubkpakp)()()(Nkpkpky)()()()(2kykykrke式中， ， ， （取整数）。p
25、TTea1)1 (1pTTeKbTN)(ku)2() 1(2)()( )1()() 1( )() 1()(222222kekekekkekkekekkukukukudip式中， 为比例系数； 为积分系数； 为微分系数。pkikdk第四节第四节 纯滞后对象的控制算法纯滞后对象的控制算法大林算法大林算法 大林算法的目标是设计一个合适的数字调节器D(z)，使整个系统的闭环传递函数相当于一个带有纯滞后的一阶惯性环节，而且要求闭环系统的纯滞后时间等于被控对象的纯滞后时间。大林算法方法比较简单，只要能设计出合适的且可以物理实现的数字调节器D(z)，就能够有效地克服纯滞后的不利影响，因而在工业生产中得到了广
26、泛应用。但它的缺点是设计中存在振铃现象，且与Smith算法一样，需要一个准确的过程数字模型，当模型误差较大时，控制质量将大大恶化，甚至系统会变得不稳定。第三节第三节 大林控制算法大林控制算法1)()()(0sTesRsYss 按照计算机控制系统直接化设计方法，Dahlin算法根据纯滞后系统的主要控制要求，将期望的闭环脉冲传递函数设计为一个带有纯滞后的一阶惯性环节，且纯滞后时间与被控对象的纯滞后时间相同。Dahlin算法的设计目标是设计一个合适的数字控制器，使整个闭环系统的传递函数相当于一个一阶惯性纯滞后环节，即式中， 为被控对象的纯滞后时间， 。为简单起见，设 为采样周期的整数倍，即 为正整数
27、。 为期望闭环传递函数的时间常数，其值由设计者用试凑法给出。NTN0T一一 大林算法基本原理大林算法基本原理大林算法的设计目标是将整个系统的传递函数设计为具有一阶惯性加纯滞后环节形式，其中纯滞后时间与被控过程的纯滞后时间相等或近似相等，即将期望闭环传递函数设计为1)(0sTess如果采用计算机控制，需要加入零阶保持器，此时期望闭环传递函数可离散化为：1) 11/) 1(/01)1 (1)-1 (11)(00zzzezesesTeZzNTTTTsTs（ 则可导出数字控制器的z传递函数如式(5-40)所示)1 (1)()1 ()(1)()()()1(1)1(NNzzzGzzzGzzD(5-40)(
28、5-38)(5-39)NTsTKesGs ,1)(11/)1(1111_111)1 ( 11)1 ( 11)(zeeKzsTTsZzzKsTKeseZzGTTTTNNNTssT1含纯滞后的一阶惯性环节的控制形式带纯滞后的一阶被控对象的传递函数为广义被控对象的脉冲传递函数为将式(5-41)代入式(5-38)，得(5-36)(5-41)1 (1)1 ()1)(1 ()()1(/1/1/00110NTTTTTTTTTTzezeeKzeezD(5-42)二二 大林算法数字控制器的基本形式大林算法数字控制器的基本形式2含纯滞后的二阶惯性环节的控制器形式带有纯滞后特性的二阶被控对象的传递函数为NTsTsT
29、KesGs ,) 1)(1()(21广义被控对象的脉冲传递函数为12121121(1)12/111( )(1)(1)1 (1)(1)(1)() (1)(1)sTNTsNNT TT TeKeG zZsTsT sKzzZs TsT sK CC zzezez式中，)(112121121TTTTeTeTTTC)(112212112)11(2TTTTTTTeTeTTTeC(5-43)得)1 (1)()1)(1)(1 ()()1(11211100210NTTTTTTTTTTzezezCCKzezeezD(5-44)132)(4 . 2sesGscsT1例例5.2已知被控对象的传递函数采样周期试用大林算法设
30、计数字控制器，并分析其输出响应和控制器的输出序列采样周期s2sT1解：根据系统的性能指标选择闭环系统的时间常数被控对象的滞后时间为2.4s，它不是T采样周期的整数倍，必须利用扩展z变换求广义被控对象的脉冲传递函数。) 13()1 (21321)(6 . 0314 . 2sseZzseseZzGsssTs11313/13/6 . 0113718. 0155. 01364. 0111)1 (2zzzzeezzz据大林算法，闭环系统的滞后时间也应该为2.4s，但为了简化运算，取纯滞后时间为2s，则闭环系统的传递函数为：sess2121)( 13605. 01395. 0)(zzz)395. 0605
31、. 01)(55. 01 ()718. 01 (09. 1)(3111zzzzzD213605. 0605. 11395. 0)()()(zzzzRzzC3211333. 0278. 0055. 11779. 0085. 1)()()(zzzzzGzCzU其脉冲传递函数为：数字控制器的脉冲传递函数为：当输入信号为单位阶跃时，系统的输出为：876543951. 0919. 0866. 0779. 0634. 0395. 0zzzzzz数字控制器的控制量输出为4321562. 0467. 0688. 0366. 0085. 1zzzzNTC(k) 0.51012345678NTu(k) 0.510
32、12345678 图5.10 系统输出波形图5.11控制量输出波形 由以上输出序列与控制序列可见，输出采样值是呈惯性平滑上升的，但控制量出现大幅振荡，如5.10所示，这种控制量以二分之一采样频率大幅度衰减的振荡现象称为振铃。)()()(zUzGzY)()()(zRzzY三、振铃现象及消除方法三、振铃现象及消除方法 数字控制器的输出以接近二分之一的采样频率大幅度上下摆动，这称为振铃现象。它对系统的输出几乎是没有影响的，但会使执行机构因磨损而造成损坏。在有交互作用的多参数控制系统中，振铃现象还有可能影响到系统的稳定性，所以在系统设计中，应该设法消除振铃现象。1振铃现象分析 振铃现象与被控对象的特性
33、、闭环时间常数、采样周期、纯滞后时间的大小等有关。 在计算机控制系统中，系统的输出和数字控制器的输出之间的关系为系统输出与闭环系统的输入的关系为)()()()(zGzzRzU)()()(zRzGzUu)()()(1)()(zRzGzDzDzU可以得出 它描述了数字控制器的输出与闭环系统的输入之间的关系，可进一步写作 表示的是数字控制器的输出与闭环系统输入之间的关系，它是分析振铃的基础。有 由分析可知，产生振铃现象的原因是数字控制器 在Z平面上 附近有极点。当 时，振铃现象最严重，在单位圆内离 越远，振铃现象越弱。)(zD1z1z1z)(zD) 1 ()0(uuRA 用振铃幅度来衡量振铃程度的强
34、弱。它的定义是，在单位阶跃输入作用下，数字控制器 的第0次输出减去第1次输出所得的差值，即3振铃现象的消除假设)(zKu写成如下形式：2211221111)(zazazbzbzKu在单位阶跃输入时，数字控制器的输出为)()()(zRzKzUu1221122111111zzazazbzb212112211)() 1(11zaazazbzb2122111)() 1(1zabazab1)0(u1) 1 (11abu由上式可见根据振铃幅度定义，有1111) 1(1) 1 ()0(baabuuRA 当 时，系统不会产生振铃现象； 时，系统产生振铃现象，其值越大，振铃幅度就越大。 由以上分析可得到一种消除
35、振铃极点的可能方法，即在控制器设计时，通过合理选取采样周期，避免在 中出现可能引起振铃的零点。大林提出一种更为简单的修正设计方法。先找出造成振铃现象的极点因子，令其中z =1，这样便消除了这个造成振铃的极点。根据终值定理，这样的处理不会影响输出的稳态值，却可改善系统动态性能，即消除了振铃现象。0RA0RA例5-2中，按照大林算法设计的数字控制器的脉冲传递函数为：)395. 0605. 01)(55. 01 ()718. 01 (09. 1)(3111zzzzzD其中极点 在负实轴上，会产生振铃现象，可用上面的方法消除振铃现象。55. 0z令因子 中的 ，代入 ，得155. 01z1z)(zD)
36、395. 0605. 01 ()718. 01 (703. 0)395. 0605. 01 (55. 1)718. 01 (09. 1)(311311zzzzzzzD这样数字控制器的控制量输出为：211605. 0605. 11457. 0637. 0)()()(zzzzGzCzU4321-48. 0496. 0522. 0565. 0637. 0zzzz可以看出，控制量序列朝一个方向衰减，不再以二分之一的采样频率振荡，因此消除了振铃现象。四、大林算法的设计步骤四、大林算法的设计步骤 （1）根据系统性能要求，确定期望闭环系统的参数，给出振铃幅度的指标。 （2）根据振铃幅度 的要求，确定采样周期
37、 ，如果有多解，则选择较大的 。RATT（3）确定整数 。TN （4）求广义对象的脉冲传递函数 及期望闭环系统的脉冲传递函数 。)(zG)(z （5）求数字控制器的脉冲传递函数 。)(zD （6）将 变换为差分方程，以便于计算机编写相应算法程序。)(zD例 5.3 已知某控制系统被控对象的传递函数为 。 试用大林算法设计数字控制器D(z)。设采样周期为T=0.5s，并讨论该系统是否会发生振铃现象。如果振铃现象出现，如何消除。1)(sesGsc211TNK，11T解：由题可知, 当被控对象与零阶保持器相连时，系统的广义对象的传递函数为 ) 1()1 ()(1)(5 . 0sseesGsesGss
38、cTs于是，可求出广义对象的数字脉冲传递函数为1315 . 05 . 031/16065. 013935. 01111)(11zzzeezzeeKzzGTTTTN大林算法的设计目标是使整个闭环系统的脉冲传递函数相当于一个带有纯滞后的一阶惯性环节。据此可设 ，可得 sT1 . 00)9933. 09933. 01)(1 ()6065. 01 (2524. 2)1 (1 )1 (3935. 06065. 01)1 (1 )1 ()(1)(2111351535311/1/1000zzzzzezezezzzezeezzGzDNTTTTTTN由上式可知， 有三个极点： ，根据前边的讨论z=1处的极点不会
39、引起振铃现象，引起振铃现象的极点为 jzzz864. 04967. 01321 ，19966. 011|5/32eezzxT依据前述讨论，要想消除振铃现象，应去掉分母中的因子 即令 ， 代入上式即可消除振铃现象。 这样，无振铃时，数字控制器的脉冲传递函数 为 )9933. 09933. 01 (21zz1z19966. 0|32 zz)(zD11111)6065. 01 (8451. 0)9933. 09933. 01)(1 ()6065. 01 (524. 2)(zzzzzD第五节、数字控制器的程序实现第五节、数字控制器的程序实现 数字控制器的设计方法要变成在计算机上实现则需要编制算法程序。
40、若状态空间形式表示 对应的差分方程时，数字控制器的状态空间方程可直接在计算机上编程。若以Z的脉冲传递函数表示 时，数字控制器 都是由软件来实现的。软件实现包括三种方法：直接程序法、串行程序法和并行程序法。)(zD)(zD)(zD一.直接程序设计法所谓直接程序法，是指 将 离散化的差分方程不做任何变化，直接编制软件的方法。直接编排结构就是按高阶传递函数分子、分母多项式系数进行编程实现。数字控制器 通常可表示为)(zD)(zDmjjjniiimmnnzbzazbzbzbzazazaazEzUzD1022112211011)()()( (5-52)01( )E()U()nmijiju kakibkj
41、 (5-53)为使计算机实现方便，把式(5-53)进行z反变换，写成差分方程的形式。mjjniijkubikeaku10)()()( (5-54)直接程序设计法的算法比较简单，可大大减少计算机运算的延时，提高系统的动态性能。但是其运算量比较大，程序占用的内存容量也比较大，程序的可读性较差，调试不方便。例例5-4 已知数字控制器脉冲传递函数D(z)为2221( )56zzD zzz试用直接程序设计法写出实现D(z)的表达式，求出D(z)的差分方程后，画出相应的程序流程图。解：根据直接程序设计法知：对给定的数字控制器的D(z)的分子、分母都乘以z22122212(21)12( )(56)156zz
42、zzzD zzzzzz用直接程序法去实现D(z)的表达式：1212( )( )2 ( )( )5 ( )6( )U zE zE z zE z zU z zU z z根据上式可得结果知： 再进行逆z变换，便可求得数字控制器的差分方程为6, 5, 1, 2, 121210bbaaamn，)2(6) 1(5)2() 1(2)()(kukukekekeku二 串行程序设计法串行程序设计法也叫迭代程序设计法。如果数字控制器的脉冲传递函数 中的零点、极点均已知或者其多项式形式可以分解为一些简单的一阶或二阶环节的串联时，就可以采用串行程序设计法来实现。)()()()(10zDdzEzUzDiki)().()
43、(210zDzDzDdknk 1= =即：1111)()()(zzzEzUzDiiiii或二阶形式2211221111)()()(zzzzzEzUzDiiiiiii（5-56）（5-57）串行程序设计框图如图5.11E(z). .D (z2)D (zk)d0D (z1)e(k)E(z)E(z)e (k1)u (k1)u (k2)U(z)u(k)图5.11 串行程序设计框图为计算 ，可先求出 ，再算出 ， ，最后算出 。分别对各个串联脉冲传递函数作z反变换，求得各个环节的输出 、 、 、 ，最后求出整个数字控制器的输出 ，即)(ku)(1ku)(2ku)(3ku)(ku1( )u k2( )u
44、k( )nu k 111 1( )( )(1)(1)u ke kz e kpu k=+-212122( )( )(1)(1)uku kz u kp uk=+-11( )( )(1)(1)mmmmmmukukz ukp uk-=+-111( )(1)(1)mmmmukukpuk+=-1( )( )(1)nuu kKukp u k-=-(5-58）例例5.5 设数字控制器的脉冲传递函数 ，试用串行程序设计法写出 的迭代表达式。2256( )920zzD zzz)(zD解：2256(2)(3)( )920(4)(5)zzzzD zzzzz根据(5-56)式，得：11111211( )212( )41
45、4( )31 3( )( )51 5( )U zzzD zzzE zzzU zD zzzU z由式(5-58)，得到控制量的输出序列为：1111( )( )2 (1)4 (1)( )( )3 (1)5 (1)u ke ke ku ku ku ku ku k三 并行程序设计法如果 可以写成部分分式的形式，数字控制器就可以用并行程序设计法来实现。此时 的一般表达式为)(zD)(zD1112121111111)(zbzazbzazbzazDnn可以把 表示成多个控制器并联的形式，即)(zD)()()()()(321zDzDzDzDzDn11111111( )( )( )1( )( )( )1nnnn
46、U za zD zE zb zUza zD zE zb z则： 由各个子函数并联组成，如图5.12所示。( )D zD (z1)D (zk)d0E(z)e (z1)u (z1)u (zk)U(z). 图5.12 并行程序设计法框图同样可以求出各个控制器输出的控制量序列为：111 1( )(1)(1)( )(1)(1)nnnnu ka e kbu ku ka e kb u k所以数字控制器的输出序列为：)()()()(21kukukukun例例5.6设数字控制器的脉冲传递函数 试用并行程序设计法写出 的表达式。127175)(2zzzzD)(zD解：首先将数字控制器的脉冲传递函数写成部分分式的形式：111124133124332127175)(zzzzzzzzzzD)()(413)()()(312)(21121111zEzUzzzDzEzUzzzD由式(5-62)，得到控制量的输出序列为：) 1(4) 1(3)() 1(3) 1(2)(2211kukekukukeku所以并行程序设计法的 表达式为：) 1(4) 1(3) 1(5)()()(2121kukukekukuku