《计控实验四-最少拍控制算法研究报告》由会员分享，可在线阅读，更多相关《计控实验四-最少拍控制算法研究报告（11页珍藏版）》请在金锄头文库上搜索。
1、东南大学自动化学院实验报告课程名称：计算机控制技术第4次实验实验名称：实验四最少拍控制算法研究院系:自动化学院专业：自动化姓名：学号：实验室：416实验组别：同组人员：实验时间：2014年4月24日评定成绩：审阅教师：、实验目的1 .学习并熟悉最少拍控制器的设计和算法；2 .研究最少拍控制系统输出采样点间纹波的形成；3 .熟悉最少拍无纹波控制系统控制器的设计和实现方法。、实验设备1. THBDC-1型控制理论计算机控制技术实验平台2. PCI-1711数据采集卡一块3. PC机1台(安装软件VC+及THJK_Server)三、实验原理1最小拍系统在采样控制系统中，通常把一个采样周期称作一拍。在
2、典型输入信号作用下，经过最少拍，使输出量采样时刻的数值能完全跟踪参考输入量的数值，跟踪误差为零的系统称为最少拍系统。计算机控制系统的方框图为：y图4-1最少拍计算机控制原理方框图根据上述方框图可知，有限拍系统的闭环脉冲传递函数为：(4-1)CD(z)G(z)R(z)1D(z)G(z)1 H(z)E(z)R(z)11 D(z)G(z)(4-2)由(4-1)、(4-2)解得:D(z)1 H(z)G(z) 1 H(z)4-3首先要使系统的过渡过程在有限拍完毕，显然，这样对系统的闭环脉冲传递函数H(z)提出了较为苛刻的要求，即其极点应位于z平面的坐标原点处。亦即希望系统的脉冲传递函数为H(z)aoz1
3、akzkF(z)4-4z式中：F(z)为H(z)的分子多项式，k为某一整数。式4-4说明H(z)的极点都在z平面的原点，系统的脉冲响应在经过了有限数k拍以后就变为零，过渡过程完毕。式4-4说明了离散系统中，为了使过渡过程较快地完毕应符合的条件。K是个有限值，它至少应该是什么数值呢？可以分析一下闭环传递函数H(z)o将式4-4代入D(z)表示式，得1 H(z)G(z) 1 H(z)1 FG(z) zk F4-5如果m和n分别为对象和保持器的组合脉冲传递函数G(z)的分子和分母的阶次，l为式4-5中F(z)的阶次，要使D(z)能实现，就应使分母的阶次大于分子的阶次knml4-6由式4-6可见，当l
4、0时，H(z)的分子F(z)F0常数，暂态响应的持续节拍数最少。式中n和m是由对象、保持器决定的，是不可变局部。这时应有knm4-7这是过渡过程所能到达的最低极限节拍数，它规定了“最少拍的极限数。2无稳态误差的最小拍系统由王勤主编教材P89P90的理论推导，可以知道，为保证系统稳态误差为零且拍数最少，应取4-81k1H(z)(1z)Fi(z)其中Fi(z)为不包含G(z)的零点和极点的多项式。为式4-8表示了无稳态误差的最少拍系统，其1H(z)满足的条件。另外，为了使系统的暂态过程在有限时间完毕，H(z)必须是z的有限多项式。这两者都要满足，因而应使Fi(z)为z 1的有限多项式。最简单的情况
5、是Fi(z)1,这时1对阶跃输入11H(z)1z1H(z)z11H(z)1z11因而D(z)1G(z)1H(z)G(z)1zG(z)z12对斜坡函数输入121H(z)(1z)H(z)2z1z2从而有D(z)121 2z z_1 2G(z) (1 z )1 2z 1_2G(z) (z 1)3对加速度函数输入131H(z)(1z)1 3z2 3z 1 -3G(z) (z 1)H(z)3z13z2z3从而有13z13z2z31-3G(z)(1z)3无纹波，无稳态误差的最少拍系统用前述方法设计的最少拍控制系统，对于符合原设计的输入信号能很快地跟踪。然而，如果进一步用改良的z变换法来研究所设计的系统，就
6、会发现问题。这种改良的z变换不仅能求出采样时刻的系统输出，而且可以研究采样间隔中，输出的变化情况。用这种z变换将发现用前述方法设计的系统，在采样时刻之间存在着波动。有纹波的系统，在采样时刻之间存在误差，而且功率损耗、振动等也很大，它将加快执行机构等可动部件的磨损。为此，必须改良设计方法，使设计出的系统满足无纹波的条件。1最少拍系统产生纹波的原因经分析可知，最少拍系统虽然经过有限拍后能使采样时刻的稳态误差为零，从而使数字控制器的离散输入量E(z)为零。但控制器的输出并没有到达稳态值，仍然是上下波动的。亦即控制器的输出U(z)不能在有限拍变为零。如果整个系统以U(z)为输出量，设这时的闭环传递函数
7、为Hd(z)。同样，如果这一闭环传递函数也能表示成极点都在z平面原点的形式，那么过渡过程也能在有限拍完毕。2无纹波最少拍系统的设计根据王勤主编教材P93的理论推导可知，无纹波最少拍系统的闭环传递函数应分别为H(z),zzF0Q4-10kz一P(z)_.一.一.式中：G(z)(-),F0(z)为z的多项式。Q(z)上述传递函数能保证系统的输出Y(z)和控制器输出U(z)的暂态过程均能在有限拍完毕。式4-9说明，无纹波最少拍系统的闭环传递函数H(z)不仅应为z1的多项式，而且应包含G(z)的全部零点。由式4-5可得Fo(z)Q(z)zkFo(z)P(z)在最简单的，情况下，Fo(z)为常数。为了保
8、证D(z)是可实现的，至少要使k大于或等于Q(z)的阶次，即kn4-11将式4-7与式4-11相比，发现由于要求无纹波，系统的最少拍增加了m拍，响应的暂态过程也延长m拍。4斜坡输入下最少拍系统设计举例1斜坡输入下无稳态误差最少拍系统设计设被控对象为二阶系统，其传递函数为G(s)s(TiS1)对于二阶被控对象加零阶保持器后对象的传递函数为:dTSGTHs(*1)选择采样周期T,将上述传递函数离散后得TG(z)K(T工T1、)z1TT(T1TgT1TeT1)z2(1z1)(1TT11e1z)因为输入是单位斜坡信号，所以选择:1H(z)(1z1)2H(z)2z1U(z)D(z)E(z)H(z)H(z
9、)(2zT1)(1eG1)T(2z1)(1e%1)K(1z)(TT1T2)TT(T1T1eT1TeT1)z111KA(1z1)(1Bz1)T12(12eT1)zKA1(B1)z1TT1z2Bz2T其中ATeT1T1eTTeT1A由此可得斜坡输入下最少拍算法的数字控制器算法为212eT/T1u(k)=(1-B)u(k-1)+Bu(k-2)+e(k)12eKAKAe(k1)T/Te1e(k2)KA2斜坡输入下无纹波、无稳态误差最少拍系统设计TTT(TT1T1eT1)z1(T1T1eT1TeT1)z2G(z)KT(1z1)(1eT1z1)KAz1(1Bz1)1(1z)(1TT1z1)保证系统在斜坡输
10、入下的稳态误差为零，那么系统的闭环传递函数应符合12一1H(z)(1z)E(z)稳定性原那么和F2(z)的形式，使它们的项数最少即最小实现而又满足上述条件。显然应为Fi(z)bozF211cz )zH(z)(1_1 一Bz )F2(z)(1 Bz1)z1(CoGz 1)1H(z)应包含G(z)的不稳定极点，其形式为12_1H(z)(1z)F1(z),其中E(z)为不包含G(z)的零点和极点的多项式。上式和保证稳态误差为零的要求一样。H(z)应包含G(z)的不稳定零点，其形式为1H(z)(1Bz)Fz(z),其中F2(z)为不包含G(z)的零点和极点的多项式。无纹波有限拍暂态过程1z 。要寻找
11、H (z)、F1 (z)为了实现无纹波，F2(z)中必须包含G(z)的分子上的另一因子联立求解得3B 2Co -2-B2 2BCi(2BB2BboB(B 1)B2 2B所以有D(z)U(z)ez1G()H(z)H(z)(1TT；z1)(CoKA(1 zqz 1)(1 boz1)LKACo(GTCoe T1)1 (bo1)z 1Ge z boz 2由此可得斜坡输入下最少拍无纹波的算法:Cou(k) (1 bo)u(k 1) bou(k 2) e(k)KAc1 coeKATT1一e(k1)TT1C1e 1KAe(k 2)四、实验步骤1、仔细阅读“PCI-1711数据采集卡驱动函数说明.doC”和“
12、THJK-Server软件使用说明.doC”文档，掌握PCI-1711数据采集卡的数据输入输出方法和THJK-Server软件及相关函数的使用方法。2、模拟对象的模拟电路图如图4-2下所示：七2rJ图4-2模拟对象电路其中：Ri=100K,R2=100K,R3=100K,Ci=10uF,Q=10uF按上图连接一个积分环节和一个惯性环节组成的二阶被控对象的模拟电路;3、系统接线图如图4-3所示：图4-3系统接线图按照上图接线。4、设计在斜坡输入下的最少拍控制器的控制算法有纹波、无纹波。5、设计在单位阶跃输入下的最少拍控制器的控制算法有纹波、无纹波。6、翻开最少拍实验文件夹下.dsw工程文件，源程
13、序中缺少最少拍数字控制器算法程序,请同学用设计好的控制器算法编写程序。7、源程序编译通过后，先启动THJK_Server图形显示软件，再执行程序代码，在显示界面出现的曲线并稳定后初始化后，把+5V电源打到“开状态，观测并记录系统的响应曲线。在实验完毕后，在键盘上先按下“e，再按下“Enter回车键键，程序退出。五、实验记录1、编程实现在斜坡输入下的最少拍控制器的控制算法有纹波、无纹波/*控制器的编程斜坡输入下有纹波最小拍控制ei=sv-fVoltage;ep=exp(-T/Ti);A=T+Ti*ep-Ti;B=(Ti-Ti*ep-T*ep)/A;op=(1-B)*opx+B*opxx+2*ei/(k*A)-(1+2*ep)*eix/(k*A)+ep*eixx/(k*A);斜坡输入下无纹波最小拍控制ei=sv-fVoltage;ep=exp(-T/Ti);A=T+Ti*ep-Ti;B=(Ti-Ti*ep-T*ep)/A;c0=(3*B+2)/(B*B+2*B+1);c1=-(2*B+1)/(B*B+2*B+1);
