计控实验四-最少拍控制算法研究

1、东南大学自动化学院实验报告课程名称:计算机控制技术评定成绩:审阅教师:第4次实验实验名称:实验四最少拍控制算法研究院（系）:自动化学院专 业:自动化姓 名:学 号:实验室:416实验组别:同组人员:实验时间:2014年4月24日、实验目的1. 学习并熟悉最少拍控制器的设计和算法；2. 研究最少拍控制系统输出采样点间纹波的形成；3. 熟悉最少拍无纹波控制系统控制器的设计和实现方法。二、实验设备1. THBDC-1型控制理论计算机控制技术实验平台2. PCI-1711数据采集卡一块3. PC 机 1 台(安装软件“ VC+ ”及“ THJK_Server ”)三、实验原理1)最小拍系统在采样控制系

2、统中，通常把一个采样周期称作一拍。在典型输入信号作用下，经过最少 拍，使输出量采样时刻的数值能完全跟踪参考输入量的数值，跟踪误差为零的系统称为最少拍系统。计算机控制系统的方框图为：CW1k图4-1最少拍计算机控制原理方框图根据上述方框图可知，有限拍系统的闭环脉冲传递函数为：H(z)他 _D( z)G(z) )R(z) 1D(z)G(z)(4-1)1 H(z)Rz1 D(z)G(z)(4-2)由(4-1)、(4-2)解得:(4-3)首先要使系统的过渡过程在有限拍内结束，显然，这样对系统的闭环脉冲传递函数H(z)提出了较为苛刻的要求，即其极点应位于z平面的坐标原点处。亦即希望系统的脉冲传递函数为F

3、(z) zk 式中：F(z)为H(z)的分子多项式，k为某一整数。式(4-4)表明H(z)的极点都在z平面的原 点，系统的脉冲响应在经过了有限数k拍以后就变为零，过渡过程结束。式(4-4)表明了离散系统中，为了使过渡过程较快地结束应符合的条件。K是个有限值，它至少应该是什么数值呢？可以分析一下闭环传递函数H (z) a。a,z 1 LakZ k(4-4)H(z)。将式(4-4)代入D(z)表示式，得D(z)丄止丄kF(z)G(z) 1 H(z) G(z) zk F(z)(4-5)如果m和n分别为对象和保持器的组合脉冲传递函数G(z)的分子和分母的阶次，I为式(4-5)中F(z)的阶次，要使D(

4、z)能实现，就应使分母的阶次大于分子的阶次k n m I(4-6)由式(4-6)可见，当l 0时，H(z)的分子F(z)F0常数，暂态响应的持续节拍数最少。式中n和m是由对象、保持器决定的，是不可变部分。这时应有k n m这是过渡过程所能达到的最低极限节拍数，它规定了“最少拍”的极限数。(4-7)2)无稳态误差的最小拍系统由王勤主编教材 P89P90的理论推导，可以知道，为保证系统稳态误差为零且拍数最 少，应取1 k1 H(z) (1 z ) F1(z)(4-8)其中F1(z)为不包含 G(z)的零点和极点的多项式。为式(4-8)表示了无稳态误差的最少拍系统，其1 H (z)满足的条件。另外，

5、为了使系统的暂态过程在有限时间内结束,H(z)必须是z 1的有限多项式。这两者都要满足，因而应使Fi(z)为z 1的有限多项式。最简单的情况是Fi (z)1，这时(1)对阶跃输入21 H (z)1H(z) z因而1 D(z)H(z)G(z) 1 H(z)1 z 1G(z) 1 z11G(z)对斜坡函数输入1 H (z)(11)2H(z) 2z1从而有D(z)1 2z 1 z 21 2G(z) (1 z )2zG(z) (z 1)(3)对加速度函数输入1 H (z)(1 z1)3H(z)3z 1 3z 2 z 3从而有D(z)G(z)3Z1 3z2 z313z2 3z 11 33(1 Z ) G

6、(z) (z 1)3)无纹波，无稳态误差的最少拍系统对于符合原设计的输入信号能很快地跟踪。然而，就会发现问题。这种改进的Z变换不仅 输出的变化情况。用这种Z变换将用前述方法设计的最少拍控制系统， 如果进一步用改进的 Z变换法来研究所设计的系统， 能求出采样时刻的系统输出，而且可以研究采样间隔中，发现用前述方法设计的系统，在采样时刻之间存在着波动。有纹波的系统，在采样时刻之间存在误差，而且功率损耗、振动等也很大，它将加快执行机构等可动部件的磨损。为此，必须改进设计方法，使设计出的系统满足无纹波的条件。(1)最少拍系统产生纹波的原因经分析可知，最少拍系统虽然经过有限拍后能使采样时刻的稳态误差为零，

7、从而使数字控制器的离散输入量 E(z)为零。但控制器的输出并没有达到稳态值， 仍然是上下波动的。亦 即控制器的输出 U(z)不能在有限拍内变为零。如果整个系统以U(z)为输出量，设这时的闭z平面原点的形环传递函数为Hd(z)。同样，如果这一闭环传递函数也能表示成极点都在 式，则过渡过程也能在有限拍内结束。(2)无纹波最少拍系统的设计根据王勤主编教材 P93的理论推导可知，无纹波最少拍系统的闭环传递函数应分别为F(z) F0(z) P(z)H(z) kkZZFo(z)Q(z)(4-9)Hd(z)(4-10)式中：G0(1，fo为z的多项式。U(z)的暂态过程均能在有限拍内结上述传递函数能保证系统

8、的输出Y(z)和控制器输出束。式(4-9)说明，无纹波最少拍系统的闭环传递函数H(z)不仅应为Z 1的多项式，而且应包含G(z)的全部零点。由式(4-5)可得D(z)丄皿Fo(z)Q(z)k 一 ， 、_ ,G(z) 1 H(z) zFo(z) P(z)在最简单的情况下,Fo (z)为常数。为了保证D(z)是可实现的，至少要使k大于或等于Q(z)的阶次，即k n(4-11)m拍，响将式(4-7)与式(4-11)相比，发现由于要求无纹波，系统的最少拍增加了 应的暂态过程也延长m拍。4)斜坡输入下最少拍系统设计举例(1)斜坡输入下无稳态误差最少拍系统设计设被控对象为二阶系统，其传递函数为G(s)s

9、(T1S 1)对于二阶被控对象加零阶保持器后对象的传递函数为:. TSG(s)ss(T1s 1)选择采样周期T,将上述传递函数离散后得TTTk(T T1 T1eT1)z1 (T1e Te1)z2(1 zT1)(1 eT1z1)因为输入是单位斜坡信号，所以选择:1 H(Z)(1z1)2H(z) 2z 1D(z) EzG(z)H(z)1 H(z)(2T1 s1z )(1e z )K(1)(T T1T(T1廿TTeT1)z1T(2 z1)。eJ1)KA(1 z 1)(1 Bz 1)2 (1KA 1 (BT1)z 1 Bz 2T2e T1)zT其中 A T l；eT1，BTT；T1e EA由此可得斜坡

10、输入下最少拍算法的数字控制器算法为u(k)=(1-B)u(k-1)+Bu(k-2)+e(k) 2eKAKAT/T1e(k1)匚(kKA2)(2)斜坡输入下无纹波、无稳态误差最少拍系统设计TTTKA(1 z1)(11 1z (1 Bz )G(z)k(T T1 。丁11 (T1 严TeT1)z2 (1 z1)(1 eJ1) 保证系统在斜坡输入下的稳态误差为零，则系统的闭环传递函数应符合1 H(z) (1 z1)2F1 (z) 稳定性原则1 H(z)应包含G(z)的不稳定极点，其形式为1 21 H (z) (1 z ) F1(z)，其中F1(z)为不包含G(z)的零点和极点的多项式。上式和保证稳态误

11、差为零的要求相同。H(z)应包含G(z)的不稳定零点，其形式为H(z) (1 Bz1)F2(z)，其中F2(z)为不包含G(z)的零点和极点的多项式。z 1。要寻找 H (z)、F1 (z) 无纹波有限拍暂态过程为了实现无纹波，F2(z)中必须包含G(z)的分子上的另一因子和F2(Z)的形式，使它们的项数最少(即最小实现)而又满足上述条件。显然应为Fi(z)boZF2(z)(CoC1z 1)z 1H(z)(1Bz1)F2(z)(1 Bz1)Z 1(CoGZ 1)联立求解得2B 1C1(2BB_ B 11)，boB(B 1)B2 2B 1所以有D(z)U(z)E(z)H(z)G(z) 1 H(z

12、)(1TT1 11e z )(Co cz )KA(1 z 1)(1 boZ 1)T1 Co (c Coe T1 )z 1 GeKA1 (bo 1)z 1 boz 2由此可得斜坡输入下最少拍无纹波的算法:u(k) (1 bo)u(k 1) bou(k 2)CoKXe(k)Tt7C1 Coe 1e(kKA1)Tc1e T1 -一e(k 2)KA软件使用说明.doc ”THJK-Server软件(及相关函数)四、实验步骤1、仔细阅读“PCI-1711数据采集卡驱动函数说明.doc 和THJK-Server 文档，掌握P CI-1711数据采集卡的数据输入输出方法和的使用方法。2、模拟对象的模拟电路图

13、如图图4-2模拟对象电路其中：R1=1OOK，R2=100K，R3=1OOK，C1=10uF，C2=10uF按上图连接一个积分环节和一个惯性环节组成的二阶被控对象的模拟电路;3、系统接线图如图 4-3所示：图4-3系统接线图按照上图接线。4、设计在斜坡输入下的最少拍控制器的控制算法（有纹波、无纹波）5、 设计在单位阶跃输入下的最少拍控制器的控制算法（有纹波、无纹波）。6、打开最少拍实验文件夹下 .dsw工程文件，源程序中缺少最少拍数字控制器算法程序, 请同学用设计好的控制器算法编写程序。7、源程序编译通过后，先启动“ THJK_Server ”图形显示软件，再执行程序代码，在显示 界面出现的曲

14、线并稳定后（初始化后） ，把+5V电源打到“开”状态，观测并记录系统的响应曲线。在实验结束后，在键盘上先按下“e”，再按下“ Enter （回车键）”键，程序退出。五、实验记录1、编程实现在斜坡输入下的最少拍控制器的控制算法（有纹波、无纹波）/* 控制器的编程斜坡输入下有纹波最小拍控制ei=sv-fVoltage;ep=ex p(-T/Ti);A=T+Ti*e p-Ti;B=(Ti-Ti*e p-T*e p)/A;op=(1-B)* op x+B*o pxx+2*ei/(k*A)-(1+2*e p)*eix/(k*A)+e p*eixx/(k*A); 斜坡输入下无纹波最小拍控制ei=sv-fV

15、oltage;ep=ex p(-T/Ti);A=T+Ti*e p-Ti;B=(Ti-Ti*e p-T*e p)/A;cO=(3*B+2)/(B*B+2*B+1);c1=-(2*B+1)/(B*B+2*B+1);bO=B*(2*B+1)/(B*B+2*B+1);op=(1-b0)* op x+b0*o pxx+cO*ei/(k*A)+(c1-cO*e p)*eix/(k*A)-c1*e p*eixx/(k*A);2、编程实现在单位阶跃输入下的最少拍控制器的控制算法（有纹波、无纹波）/* 控制器的编程单位阶跃输入下有纹波最小拍控制ei=sv-fVoltage;ep=ex p(-T/Ti);A=T+

16、Ti*e p-Ti;B=(Ti-Ti*e p-T*e p)/A;op=(-B)*o px+ei/(k*A)-e p*eix/(k*A);单位阶跃输入下无纹波最小拍控制ei=sv-fVoltage;ep=ex p(-T/Ti);A=T+Ti*e p-Ti;B=(Ti-Ti*e p-T*e p)/A; c0=1/(B+1);b0=B/(B+1);op=(-b0)* op x+(c0/(k*A)*(ei-e p*eix); eixx=eix;eix=ei;opxx=op兀opx=op;3、绘制在斜坡输入下的最少拍控制的响应曲线（有纹波与无纹波）并比较分析。答：有纹波-S-6 =2-2-4-S 3-1

17、0-1tiJ.4Lkti匕/吋间/Id73.075.B784-无纹波ia =A2- 帛0-2-T-6 a-1D-|27阿冋图2分析：由图1可看出有纹波系统的输出围绕给定输入上下波动，:6230.422 6如弓31 433 .S 35 3对应的控制器输出值则出现30正负交替的波形，通过上面两图的比较， 我们可以看出斜坡输入下有纹波最小拍控制的稳态 误差更大，在采样时刻之间存在误差，而且功率损耗、振动等也很大，它将加快执行机构等可动部件的磨损，这显然不是我们希望的控制效果，而无纹波最小拍控制的稳态误差几乎没有。4、绘制在单位阶跃输入下的最少拍控制的响应曲线（有纹波与无纹波）并比较分析。有纹波3-評一2-3-3D32.234.43G.04L图343.2胡7.643.9无纹波示加2-2 Ib.ty?e.juby.CM2 y.a?u时冋JU.Ll.=b12.JUx14.U59分析：通过上面两图的比较，我们可