1、最少拍有纹波控制系统设计最少拍有纹波控制系统设计计算机控制系统脉冲传函计算机控制系统脉冲传函D(z)Ho(s)Gc(s)e*(t)u*(t)E(z)U(z)r(t)+_R(z)(z)G(z)c(t)C(z);)(1)()()(0 sGseZsGsHZzGCTsCG(z)广义对象的脉冲传递函数：广义对象的脉冲传递函数：典型的数字反馈系统典型的数字反馈系统：sesHTso1)(D(z)Ho(s)Gc(s)e*(t)u*(t)E(z)U(z)r(t)+_R(z)(z)G(z)c(t)C(z)H0(s)零阶保持器：零阶保持器：开环冲传递函数：开环冲传递函数：()()()KzD zG z闭环脉冲传递函数
2、：闭环脉冲传递函数：()()()()1()1()()KKzD z G zzzD z G z误差脉冲传递函数：误差脉冲传递函数：()1()1()()1()()eE zzzR zD z G z D(z)Ho(s)Gc(s)e*(t)u*(t)E(z)U(z)r(t)+_R(z)(z)G(z)c(t)C(z)数字控制器输出闭环脉冲传递函数为：数字控制器输出闭环脉冲传递函数为：)()()()(1)()()()(zGzzGzDzDzRzUzU若已知脉冲传递函数，可以计算出若已知脉冲传递函数，可以计算出D(z)已知已知(z)，可计算出，可计算出D(z)：1()()()1()zD zG zzD(z)Ho(s
3、)Gc(s)e*(t)u*(t)E(z)U(z)r(t)+_R(z)(z)G(z)c(t)C(z)已知已知e(z)，可计算出，可计算出D(z)：1()1()()()eezD zG zz已知已知U(z)，可计算出，可计算出D(z)：)(1)()()(1)()(zzzGzzzDUUUD(z)必须满足以下条件：必须满足以下条件：由此而得到的由此而得到的D(z)是物理可实现的是物理可实现的 D(z)也必须是稳定的也必须是稳定的数字控制器的设计步骤：数字控制器的设计步骤：（1）根据被控对象的传递函数，求出系统的广义对象的传递）根据被控对象的传递函数，求出系统的广义对象的传递函数函数:)(1)(sGses
4、GcTs（2）求出求出G(s)所对应的广义对象脉冲传递函数所对应的广义对象脉冲传递函数G(z)：)()1()(1)()(1ssGzzsGseZsGZzGccTs（3）根据控制系统的性能指标及其输入条件，确定出整个闭环）根据控制系统的性能指标及其输入条件，确定出整个闭环系统的脉冲传递函数系统的脉冲传递函数)(z（4）确定数字控制器的脉冲）确定数字控制器的脉冲传递函数传递函数D(z):)()()(1)()()()(1)()()(zzGzzzGzzzGzzDeee（5）对最少拍无波纹系统，验证是否有波纹存在；对于带）对最少拍无波纹系统，验证是否有波纹存在；对于带纯滞后的惯性环节，还要看其是否出现振铃
5、现象。纯滞后的惯性环节，还要看其是否出现振铃现象。（6）根据采用周期、时间常数及其他条件求出相应的系数，）根据采用周期、时间常数及其他条件求出相应的系数，并将其转换成计算机能够接受的数据形式。并将其转换成计算机能够接受的数据形式。最少拍系统设计中必须满足准确性、快速性、稳定性和最少拍系统设计中必须满足准确性、快速性、稳定性和物理可实现性的要求，下面讨论最少拍有纹波控制系统的物理可实现性的要求，下面讨论最少拍有纹波控制系统的设计方法。设计方法。G(z)：G(z)中不包含单位圆上和单位圆外的零极点。中不包含单位圆上和单位圆外的零极点。u:G(z)中单位圆上和单位圆外的零点数。中单位圆上和单位圆外的
6、零点数。v:G(z)中单位圆上和单位圆外的极点数。中单位圆上和单位圆外的极点数。1.最少拍有纹波控制器设计的约束条件：最少拍有纹波控制器设计的约束条件：G(z)中有单位圆上和单位圆外的极点时：中有单位圆上和单位圆外的极点时：选择：选择：)()1()(111zFzazviieG(z)中有单位圆上和单位圆外的零点时：中有单位圆上和单位圆外的零点时：选择：选择：)()1()(211zFzbzzuiid 2.F1(z),F2(z)的确定：的确定：如果如果G(z)有有 j 个极点在单位圆上个极点在单位圆上(z=1),则则：)()1()1()(:1111zFzzazqjqjviie)()1()1()(:1
7、111zFzzazqjjjviie m=u+d n=v j+q m=u+d n=v 最少拍控制器的设计过程最少拍控制器的设计过程：根据根据R(z)、G(z)v,u,d,j,q m,n F1(z),F2(z)e(z)(z)f11,f21(z)=1-e(z)()()()(zzGzzDeU(z),Y(z),E(z)211111111)1()()1)(718.01)(68.31()56.11)(368.01(63.0)(zTzzRzzzzzzzGv=2,(圆上圆外极点数圆上圆外极点数)u=1,(圆上圆外零点数圆上圆外零点数)d=0,(z-1重数重数)j=1,(单位圆上极点数单位圆上极点数)q=2,(R
8、(z)中中(1-z-1)q)j q,m=u+d=1,n=v j+q=332322212121111)(1)(zfzfzfzFzfzF离散控制系统如图所示，被控对象为Gc(s)，采用零阶保持器，采样周期T=1s，试设计当输入信号为单位速度信号时的最少拍有纹波控制器。离散控制系统如图所示，被控对象为Gc(s)，采用零阶保持器，采样周期T=0.5s，试设计当输入信号为单位阶跃信号时的最少拍有纹波控制器。离散控制系统如图所示，被控对象为Gc(s)，采用零阶保持器，采样周期T=0.1s，试设计当输入信号为单位阶跃信号时的最少拍有纹波控制器。大林算法大林算法大林算法v在控制系统设计中，纯滞后往往是影响系统
9、动态特在控制系统设计中，纯滞后往往是影响系统动态特性的不利因素，这种系统如果控制器设计不当，常性的不利因素，这种系统如果控制器设计不当，常常会引起系统产生大的超调或振荡。常会引起系统产生大的超调或振荡。v大林算法是一种专门针对工业生产过程中含有纯滞大林算法是一种专门针对工业生产过程中含有纯滞后控制对象的直接数字设计算法。后控制对象的直接数字设计算法。大林算法大林算法v大林算法的设计目标是：大林算法的设计目标是：设计数字控制器，使系统的闭环传函为具有纯设计数字控制器，使系统的闭环传函为具有纯滞后的一阶惯性环节，而且要求闭环系统的纯滞滞后的一阶惯性环节，而且要求闭环系统的纯滞后时间等于被控对象的滞
10、后时间，即：后时间等于被控对象的滞后时间，即：期望的闭环传递函数为：期望的闭环传递函数为：NTsTess,1)(大林算法闭环传递函数为：大林算法闭环传递函数为：NTsTess,1)()(1)()()(sseZzRzYzsT采用零阶保持器，采样周期为采用零阶保持器，采样周期为T。系统闭环传递函数为。系统闭环传递函数为:1)1(111/1/zezesTesezTTNTTNTsTs控制器的传递函数为：控制器的传递函数为：)1(1)1()(1)(-1)()(1)(111/NNzezeezzGzzzGzDTTTTTT当被控对象为带有纯滞后的一阶惯性环节时：当被控对象为带有纯滞后的一阶惯性环节时：)1(1
11、)1()1)(1()(111/1/1/NzezeeKzeezDTTTTTTTTTT-11z 11)(11-N-11/1/zeeKsTKeseZzGTTTTsNTsT代入代入D(z)中中，得：，得：当被控对象为带有纯滞后的二阶惯性环节时：当被控对象为带有纯滞后的二阶惯性环节时：121112/11121()(1)(1)(1)(1)TssNT TT TCC zeKeG zZKzsTsT sezez12/11/11(1)12(1)(1)(1)()()1(1)T TT TT TT TT TNeezezD zK CC zezez12/1122111()T TT TCTeT eTT 1221(1/1/)/2
12、12211()TTTT TT TCeTeT eTT代入代入D(z)中中，得：，得：大林算法 已知被控对象传递函数为：，期望闭环传递函数的惯性时间常数 ，试根据大林算法，求数字控制器D(z)。2()(1)seG ss s1Ts2Ts大林算法 已知被控对象传递函数为：，期望闭环传递函数的惯性时间常数 ，试根据大林算法，求数字控制器D(z)。解：系统的广义对象脉冲传函为：2()(1)seG ss s2131110.368(1 0.718)()(1)(1)(1 0.368)sseezG zZzss szz1Ts2Ts 根据大林算法：期望闭环传递函数的惯性时间常数 ，因此，系统的闭环脉冲传递函数为：2T
13、s1-3-20.607z-10.393z12111Z )(1)()()(seseZsTesesseZzRzYzssTssTsT大林算法 数字控制器D(z)为单位阶跃输入下闭环系统的输出为：控制量的输出为：可以看出，控制量出现大幅度衰减振荡。111131.068(1)(1 0.368)()(10.718)(1 0.6070.393)zzD zzzz345()()()0.3930.6320.775Y zR zzzzz123()()1.0680.5120.5230.281()Y zU zzzzG z大林算法 振铃现象产生的规律振铃现象产生的规律：极点距离 z=-1越近，振铃现象越严重；单位圆内右半平
14、面的零点会加剧振铃现象；单位圆内右半平面的极点会减弱振铃现象。振铃现象的消除：振铃现象的消除：方法方法：找出D(z)中引起振铃的因子（z=-1附近的极点），令其中的z=1。大林算法 对前例 显然 z=-0.718是一个接近 z=-1的极点，它是引起振铃现象的主要原因。在因子(1+0.718z-1)中令 z=1，得到新的D(z)为：因此 111131.068(1)(1 0.368)()(10.718)(1 0.6070.393)zzD zzzz11131.068(1)(1 0.368)()1.718(1 0.6070.393)zzD zzz345()()()0.2290.5320.716Y zR
15、 zzzzz123()()0.6220.1490.0900.158()Y zU zzzzG z大林算法 例：对被控对象 采样周期T=0.5s，惯性时间常数T=0.5s，试用大林算法设计数字控制器的D(z)。)1()(1ssesGsc大林算法 例：对被控对象 采样周期T=0.5s，惯性时间常数T=0.5s，试用大林算法设计数字控制器的D(z)。)1()(1ssesGsc)6065.01)(1()8474.01(1065.01111)1(5.0)1(1111)1()1(1)(1)(111315.01211122121zzzzzezzzzzsssZzzsseseZsGsezzGsTscTs解：系统的
16、广义对象脉冲传函为：1/11111)(1)(zeeZsTeseZssezzTTTTNNTsTsTs15.0)(1sess根据大林算法，设 T=0.5s，则：因为T=0.5s，T=0.5s，N=2，所以：13368.0-1632.0)(zzz大林算法 所求数字控制器D(z)为：)632.0368.01)(8474.01()6065.01)(1(934.5)368.01632.01()6065.01()8474.01()1(1065.0368.0-1632.0)(1)()()(3111113111313zzzzzzzzzzzzzzzGzzD大林算法 3111632.0368.01)6065.01)(1(212.3)(zzzzzDD(z)中，有一个极点z=-0.8474，非常靠近z=-1，如果不采用大林算法进行修正系统将会出现振铃现象，因此令：（1+0.8474z-1）中z=1，则：大林算法 大林算法
