东南大学计算机控制技术实验报告四

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

2、C+ ”及“ THJK_Server ”)三、实验原理1)最小拍系统在采样控制系统中，通常把一个采样周期称作一拍。在典型输入信号作用下，经过最少跟踪误差为零的系统称为最少拍，使输出量采样时刻的数值能完全跟踪参考输入量的数值, 拍系统。计算机控制系统的方框图为:tk图4-1最少拍计算机控制原理方框图根据上述方框图可知，有限拍系统的闭环脉冲传递函数为:(4-1)H 込R(z) 1+ D(z)G(z)5)燈1 + D G(z)(4-2)由(4-1)、(4-2)解得:D(z)=二G(z) 1H (z)H(z)(4-3)首先要使系统的过渡过程在有限拍内结束，显然，这样对系统的闭环脉冲传递函数H (z)提

3、出了较为苛刻的要求，即其极点应位于z平面的坐标原点处。亦即希望系统的脉冲传递函数为F(z)kz式中：F(z)为H(z)的分子多项式，k为某一整数。式(4-4)表明H(z)的极点都在z平面的原1kH(Z)=a0+6Z屮+akz-(4-4)点，系统的脉冲响应在经过了有限数 k拍以后就变为零，过渡过程结束。式(4-4)表明了离散系统中，为了使过渡过程较快地结束应符合的条件。K是个有限值，它至少应该是什么数值呢？可以分析一下闭环传递函数H(z)。将式(4-4)代入D(z)表示式，得D-胡tIzt胡凡(4-5)如果m和n分别为对象和保持器的组合脉冲传递函数G(z)的分子和分母的阶次，丨为式(4-5)中F

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

5、-H(Z)满足的条件。另外，为了使系统的暂态过程在有限时间内结束,H(z)必须是z，的有限多项式。这两者都要满足，因而应使F/z)为z，的有限多项式。最简单的情况是Fi (z) =1，这时(1)对阶跃输入1-H =1-z，1Z1从而有D(z)2Z j =丄G(z) (1-z ) G(z)2z-1因而 D(z)= 1 . H(z) = 1十 1G(z) 1-H(z)G(z) 1z 丄 G(z)(2 )对斜坡函数输入1-H(z) =(1-zj2H(z) =2z-z，(3)对加速度函数输入1-H(z) =(1-z)3H(z)=3z-3z/ +Z从而有D(z)=G(z)3z -3z +z7a(1z )

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

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

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

9、统，其传递函数为G(s) =K一朴+1)对于二阶被控对象加零阶保持器后对象的传递函数为:G3 ss(T1s + 1)选择采样周期T,将上述传递函数离散后得G(z)=K 仃一 T1+T1eT1)宀(TreT亍1_T(1-zr(1e T1Z 冷因为输入是单位斜坡信号，所以选择:1-H(z)=(1-z/)2H (z) =2zz工U (z)1D(z) = I)T(2-z)(1-eFz)甌=丽 1-H(z)T (2上)(1仝9)1匚匚_ KA(1-z)(1 + Bz)K(1 -z)(T+6 T1) +(一e T1 -Te T1 )zTT_12-(1 +2eF)z 丄KA 1+ (B-1)z4-Bz/TT

10、T芍其中 A=T+t；-T1，B1T1e1Te1由此可得斜坡输入下最少拍算法的数字控制器算法为u(k)=(1-B)u(k-1)+Bu(k-2)+ 厶(k) -1 +2e e(k -1) + kakae e(k -2) ka(2)斜坡输入下无纹波、无稳态误差最少拍系统设计TTTG(z) = K仃眄宀仃宀訂一荷). KA二(1-zJ(1-eF)z(1 + BzV)T(1-z)(1-eJ)保证系统在斜坡输入下的稳态误差为零，则系统的闭环传递函数应符合1-H(z)=(1-z)2F1(z)稳定性原则1 -H (z)应包含G(z)的不稳定极点，其形式为1 -H (z) =(1-z)2F1 (z)，其中FJ

11、z)为不包含G(z)的零点和极点的多项式。上式和保证稳态误差为零的要求相同。H (z)应包含G(z)的不稳定零点，其形式为H(z)=(1 +Bz1)F2(z)，其中F2(z)为不包含G(z)的零点和极点的多项式。无纹波有限拍暂态过程为了实现无纹波，F2(z)中必须包含G(z)的分子上的另一因子 z1。要寻找H(z)、F1(z)和F2(z)的形式，使它们的项数最少(即最小实现)而又满足上述条件。显然应为F1(z)=b0Z+1F2(z) =(C0 +卯二)zH( z) = (1 + Bz)F2 (z) = (1 + Bz)z(c。+ qz)联立求解得3B +2C0 B2 +2B +1Ci-(2B+

12、1)，B2 +2B +1b。_ B(B +1)=B2 +2B +1所以有D(zHEg=G1Z)(1 -e T1z)(C0 +C1Z)1-H(z) _ KA(1-z)(1+b0z)H(z)TT1 C0+(G-Ge T1)zGe T1z_KA1+(b0 -1)z-b0Z由此可得斜坡输入下最少拍无纹波的算法:U(k) =(1-bo)u(k -1) +bou(k -2)T+Se(k)+HSe(k-TC1e 丁1e(k -2)KA四、实验步骤1、仔细阅读PCI-1711数据采集卡驱动函数说明.doc 和THJK-Server软件使用说明.doc ”文档，掌握PCI-1711数据采集卡的数据输入输出方法和

13、THJK-Server软件（及相关函数）的使用方法。2、模拟对象的模拟电路图如图4-2下所示：图4-2模拟对象电路其中：Ri=100K，R2=100K，R3=100K，Ci=10uF，C2=10uF按上图连接一个积分环节和一个惯性环节组成的二阶被控对象的模拟电路;3、系统接线图如图 4-3所示:按照上图接线。4、设计在斜坡输入下的最少拍控制器的控制算法（有纹波、无纹波） 5、设计在单位阶跃输入下的最少拍控制器的控制算法（有纹波、无纹波） 6、打开最少拍实验文件夹下 .dsw工程文件，源程序中缺少最少拍数字控制器算法程序, 请同学用设计好的控制器算法编写程序。7、源程序编译通过后，先启动“THJ

14、K_Server ”图形显示软件，再执行程序代码，在显示界面出现的曲线并稳定后（初始化后），把+5V电源打到“开”状态，观测并记录系统的响应曲线。在实验结束后，在键盘上先按下“e”，再按下“ Enter （回车键）”键，程序退出。五、实验结果1编程实现在斜坡输入下的最少拍控制器的控制算法（有纹波、无纹波）（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 +

15、 2 ? ep) ? eix/(k ?A)+ep ? eixx/(k ? A);（2 ）无纹波算法最小拍控制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);bO = B ? (2 ? B + 1)/(B ?B+ 2 ?B+ 1);op = (1 - b0) ? opx + b0 ? opxx + c0 ? ei/(k ?A) + (c1 - c

16、0 ? ep) ? eix/(k ? A)-c1?ep? eixx/(k ?A);2、编程实现在单位阶跃输入下的最少拍控制器的控制算法(有纹波、无纹波)(1)有纹波最小拍控制 ei = SV - fVoltage;ep = exp(- T/Ti);A= T+ Ti ?ep- Ti;B = (Ti - Ti ?ep- T?ep)/A;op = (- B) ? opx + ei/(k ?A) - ep? eix/(k ?A);(2)无纹波最小拍控制 ei = SV - fVoltage;ep = exp(- T/Ti);A = T + Ti ? ep- Ti;B = (Ti - Ti ?ep- T?ep)/A;c0 = 1/(B + 1);b0 = B/(B + 1);op = (- b0) ?opx + (c0/(k ?A) ?(ei - ep ?eix);eixx = ei x;eix = ei;opxx = opx;opx = op;3、绘制在斜坡输入下的最少拍控制的响应曲线(有纹波与无纹波)并比较分析。斜坡输入，输出有纹波图3.113&-42-4-Ffc -m-LO12.2L.4.4b5.ftL3.325.427.fl20. B32图3.2斜坡输入，输出无纹波分析：通过对比输出波形可以看出斜坡输入下，有纹波最小拍控制系统在采样时刻之间存在 误差，有明