《2022年计控实验四-最少拍控制算法研究》由会员分享,可在线阅读,更多相关《2022年计控实验四-最少拍控制算法研究(10页珍藏版)》请在金锄头文库上搜索。
1、东南大学自动化学院实 验 报 告课程名称:电脑控制技术第4 次实验实验名称:实验四最少拍控制算法研究院 系 :自动化学院专业:自动化姓名:学号:实 验 室:416 实验组别:同组人员:实验时间:2014年 4 月 24 日评定成绩:审阅教师:精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 1 页,共 10 页一、实验目的1学习并熟悉最少拍控制器的设计和算法;2研究最少拍控制系统输出采样点间纹波的形成;3熟悉最少拍无纹波控制系统控制器的设计和实现方法。二、实验设备1THBDC-1 型 控制理论电脑控制技术实验平台2PCI-1711 数据采集卡一块3P
2、C 机 1 台(安装软件“ VC+ ”及“ THJK_Server ”) 三、实验原理1最小拍系统在采样控制系统中,通常把一个采样周期称作一拍。在典型输入信号作用下,经过最少拍,使输出量采样时刻的数值能完全跟踪参考输入量的数值,跟踪误差为零的系统称为最少拍系统。电脑控制系统的方框图为:图 4-1 最少拍电脑控制原理方框图根据上述方框图可知,有限拍系统的闭环脉冲传递函数为:)()(1)()()()()(zGzDzGzDzRzCzH(4-1) )()(11)()()(1zGzDzRzEzH(4-2) 由(4-1) 、(4-2)解得:)(1)()(1)(zHzHzGzD4-3首先要使系统的过渡过程在
3、有限拍内结束,显然,这样对系统的闭环脉冲传递函数)(zH提出了较为苛刻的要求,即其极点应位于z 平面的坐标原点处。亦即希望系统的脉冲传递函数为101( )( )kkkF zH zaa za zz4-4式中: F(z)为 H(z)的分子多项式,k 为某一整数。式4-4说明 H(z)的极点都在z 平面的原点,系统的脉冲响应在经过了有限数k 拍以后就变为零,过渡过程结束。式4-4说明了离散系统中,为了使过渡过程较快地结束应符合的条件。K 是个有限值, 它至少应该是什么数值呢?可以分析一下闭环传递函数H(z)。 将式 4-4精选学习资料 - - - - - - - - - 名师归纳总结 - - - -
4、 - - -第 2 页,共 10 页代入 D(z)表示式,得)()()(1)(1)()(1)(zFzzFzGzHzHzGzDk4-5如果 m 和 n 分别为对象和保持器的组合脉冲传递函数G(z)的分子和分母的阶次,l为式 4-5中 F(z)的阶次,要使D(z)能实现,就应使分母的阶次大于分子的阶次knml4-6由式 4-6可见,当0l时, H(z) 的分子0( )F zF常数,暂态响应的持续节拍数最少。式中n 和 m 是由对象、保持器决定的,是不可变部分。这时应有knm4-7这是过渡过程所能到达的最低极限节拍数,它规定了“最少拍”的极限数。2无稳态误差的最小拍系统由王勤主编教材P89P90 的
5、理论推导,可以知道,为保证系统稳态误差为零且拍数最少,应取111( )(1)( )kH zzF z4-8其中1( )F z为不包含G(z)的零点和极点的多项式。为式4-8表示了无稳态误差的最少拍系统,其1( )H z满足的条件。另外,为了使系统的暂态过程在有限时间内结束,H(z)必须是1z的有限多项式。 这两者都要满足,因而应使1( )F z为1z的有限多项式。 最简单的情况是1( )1F z,这时1对阶跃输入11( )1H zz1( )H zz因而11)(11)(1)(1)()(1)(11zzGzzzGzHzHzGzD2对斜坡函数输入1 21( )(1)H zz或12( )2H zzz从而有
6、22121) 1(12)(1)1(2)(1)(zzzGzzzzGzD3对加速度函数输入1 31( )(1)H zz精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 3 页,共 10 页或123( )33H zzzz从而有3231321) 1(133)(1)1 (33)(1)(zzzzGzzzzzGzD3无纹波,无稳态误差的最少拍系统用前述方法设计的最少拍控制系统,对于符合原设计的输入信号能很快地跟踪。然而,如果进一步用改良的z 变换法来研究所设计的系统,就会发现问题。 这种改良的z 变换不仅能求出采样时刻的系统输出,而且可以研究采样间隔中,输出的变化
7、情况。用这种 z 变换将发现用前述方法设计的系统,在采样时刻之间存在着波动。有纹波的系统,在采样时刻之间存在误差,而且功率损耗、振动等也很大,它将加快执行机构等可动部件的磨损。为此,必须改良设计方法,使设计出的系统满足无纹波的条件。1最少拍系统产生纹波的原因经分析可知, 最少拍系统虽然经过有限拍后能使采样时刻的稳态误差为零,从而使数字控制器的离散输入量E(z)为零。 但控制器的输出并没有到达稳态值,仍然是上下波动的。亦即控制器的输出U(z)不能在有限拍内变为零。如果整个系统以U(z)为输出量,设这时的闭环传递函数为( )DHz。同样,如果这一闭环传递函数也能表示成极点都在z 平面原点的形式,则
8、过渡过程也能在有限拍内结束。2无纹波最少拍系统的设计根据王勤主编教材P93 的理论推导可知,无纹波最少拍系统的闭环传递函数应分别为0( ) ( )( )( )kkFz P zF zH zzz4-90( )( )( )DkFz Q zHzz4-10式中:( )( )( )P zG zQ z,0( )Fz为 z 的多项式。上述传递函数能保证系统的输出Y(z) 和控制器输出U(z)的暂态过程均能在有限拍内结束。式4-9说明, 无纹波最少拍系统的闭环传递函数H(z)不仅应为1z的多项式, 而且应包含 G(z)的全部零点。由式 4-5可得)()()()()(1)()(1)(00zPzFzzQzFzHzH
9、zGzDk在最简单的情况下,0( )F z为常数。 为了保证D(z)是可实现的, 至少要使 k 大于或等于Q(z)的阶次,即kn4-11将式 4-7与式 4-11相比,发现由于要求无纹波,系统的最少拍增加了m 拍,响应的暂态过程也延长m 拍。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 4 页,共 10 页4斜坡输入下最少拍系统设计举例1斜坡输入下无稳态误差最少拍系统设计设被控对象为二阶系统,其传递函数为)1()(1sTsKsG对于二阶被控对象加零阶保持器后对象的传递函数为:)1(1)(1sTsKsesGTS选择采样周期T,将上述传递函数离散后得
10、)1)(1()()()(112111111111zezzTeeTTzeTTTKzGTTTTTTTT因为输入是单位斜坡信号,所以选择:21)1 ()(1zzH212)(zzzH)1)(1()1)(2()()(1()1)(2()(1)()(1)()()(11111111111111111BzzKAzezzTeeTTeTTTzKzezzHzHzGzEzUzDTTTTTTTTTT2121) 1(1)21(2111BzzBzezeKATTTT其中111TeTTATT,ATeeTTBTTTT1111由此可得斜坡输入下最少拍算法的数字控制器算法为u(k)=(1-B)u(k-1)+Bu(k-2)+)2() 1
11、(21)(211/keKAekeKAekeKATTTT2斜坡输入下无纹波、无稳态误差最少拍系统设计)1)(1()1()1)(1()()()(11111121111111111zezBzzKAzezzTeeTTzeTTTKzGTTTTTTTTTT保证系统在斜坡输入下的稳态误差为零,则系统的闭环传递函数应符合1211( )(1)( )H zzF z稳定性原则精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 5 页,共 10 页1( )H z应包含 G(z)的不稳定极点,其形式为1211( )(1)( )H zzF z,其中1( )F z为不包含G(z)的
12、零点和极点的多项式。上式和保证稳态误差为零的要求相同。( )H z应包含 G(z)的不稳定零点,其形式为12( )(1)( )H zBzFz,其中2( )Fz为不包含G(z)的零点和极点的多项式。无纹波有限拍暂态过程为了实现无纹波,2( )Fz中必须包含G(z)的分子上的另一因子1z。 要寻找( )H z、1( )F z和2( )F z的形式,使它们的项数最少即最小实现而又满足上述条件。显然应为110( )1F zb z11201( )()F zcc zz1111201( )(1)( )(1)()H zBzFzBzzcc z联立求解得122320BBBc,12) 12(21BBBc,12) 1
13、(20BBBBb所以有)1)(1()(1()(1)()(1)()()(10111011zbzKAzcczezHzHzGzEzUzDTT2010211010)1(1)(111zbzbzeczecccKATTTT由此可得斜坡输入下最少拍无纹波的算法:)2() 1()()2()1()1()(11101000keKAeckeKAecckeKAckubkubkuTTTT四、实验步骤1、仔细阅读 “”和“”文档, 掌握 PCI-1711 数据采集卡的数据输入输出方法和THJK-Server软件及相关函数的使用方法。2、模拟对象的模拟电路图如图4-2 下所示:图 4-2 模拟对象电路精选学习资料 - - -
14、 - - - - - - 名师归纳总结 - - - - - - -第 6 页,共 10 页其中: R1=100K , R2=100K ,R3=100K , C1=10uF,C2=10uF 按上图连接一个积分环节和一个惯性环节组成的二阶被控对象的模拟电路;3、系统接线图如图4-3 所示:图 4-3 系统接线图按照上图接线。4、设计在斜坡输入下的最少拍控制器的控制算法有纹波、无纹波。5、设计在单位阶跃输入下的最少拍控制器的控制算法有纹波、无纹波。6、打开最少拍实验文件夹下dsw 工程文件,源程序中缺少最少拍数字控制器算法程序,请同学用设计好的控制器算法编写程序。7、源程序编译通过后,先启动“THJ
15、K_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*
16、A); 斜坡输入下无纹波最小拍控制ei=sv-fVoltage; ep=exp(-T/Ti); A=T+Ti*ep-Ti; B=(Ti-Ti*ep-T*ep)/A; 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 7 页,共 10 页c0=(3*B+2)/(B*B+2*B+1); c1=-(2*B+1)/(B*B+2*B+1); b0=B*(2*B+1)/(B*B+2*B+1); op=(1-b0)*opx+b0*opxx+c0*ei/(k*A)+(c1-c0*ep)*eix/(k*A)-c1*ep*eixx/(k*A); 2、编程实现在单位阶跃输
17、入下的最少拍控制器的控制算法有纹波、无纹波。/*控制器的编程单位阶跃输入下有纹波最小拍控制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); 单位阶跃输入下无纹波最小拍控制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=eix; eix
18、=ei; opxx=opx; opx=op; 3、绘制在斜坡输入下的最少拍控制的响应曲线有纹波与无纹波并比较分析。答:有纹波图 1 无纹波精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 8 页,共 10 页图 2 分析:由图 1 可看出有纹波系统的输出围绕给定输入上下波动,对应的控制器输出值则出现正负交替的波形, 通过上面两图的比较,我们可以看出斜坡输入下有纹波最小拍控制的稳态误差更大, 在采样时刻之间存在误差,而且功率损耗、振动等也很大,它将加快执行机构等可动部件的磨损, 这显然不是我们希望的控制效果,而无纹波最小拍控制的稳态误差几乎没有。4、绘
19、制在单位阶跃输入下的最少拍控制的响应曲线有纹波与无纹波并比较分析。有纹波图 3 无纹波精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 9 页,共 10 页分析: 通过上面两图的比较,我们可以看出阶跃输入下有纹波最小拍控制的稳态误差更大,在图上可以看出误差输出曲线都是振荡的,而无纹波最小拍控制的稳态误差到第二拍结束时就使误差为零。 而且两者的动态指标也不相同:有纹波有超调量,无纹波的则没有超调,调节时间有纹波比无纹波的要大。六、实验总结本次实验研究了最小拍控制,同时又研究了在不同输入下单位阶跃输入和斜坡输入最小拍控制的响应曲线有纹波与无纹波。通过比较,我们可以看出,由于纹波的存在,在相同输入条件下,有纹波系统的响应曲线存在明显的波动,而无纹波系统则很好的解决了这一问题,改善了系统的各动态性能指标,具有很好的控制性能。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 10 页,共 10 页