试验最少拍控制算法matlab仿真试验2017

1、实验三、最少拍控制算法matlab仿真实验一、 实验目的：.掌握最少拍有纹波、无纹波系统的设计方法；.学会最少拍控制系统的分析方法；. 了解输入信号对最少拍控制系统的影响及改进措施二、实验内容：设单位反馈线性定常系统的连续部分和零阶保持器的传递函数分别为：、1、1 -*G0(s) =，Gh(s)=s(s 1)sD(z)。讨,采样周期T=1秒。要求系统在单位阶跃输入时实现最少拍控制，求数字控制器的脉冲传递函数论加上控制器后的系统输出在采样点之间是否存在纹波，若存在纹波，请设计最少拍无纹波控制器。系统 结构如图：clear all ; clc;G=zpk(,0,-1,1);T=1;Gd=c2d(G

2、,T,zoh)num1,den1=tfdata(Gd, GG=feedback(Gd,1)1、对未加控制器时的闭环系统分析其是否满足最少拍要求%零极点形式的被控对象%采样周期%广义被控对象脉冲传函v);%原系统闭环传递函数yd_1=0;yd_2=0;ed_1=0;ed_2=0;for k=1:1:35time(k)=k*T;rin(k)=1;%单位阶跃信号%闭环系统的差分方程yd(k)=num1(2)*ed_1+num1(3)*ed_2-den1(2)*yd_1-den1(3)*yd_2;ed(k)=rin(k)-yd(k);%t 偏差 %数据更新ed_2=ed_1;ed_1=ed(k);yd

3、_2=yd_1;yd_1=yd(k);endyd%原闭环系统的输出序列ed%原闭环系统的偏差序列figure(1);plot(time,rin,r ,time,yd, g ,time,ed, b);xlabel( time );ylabel( rin,yd );title(闭环系统在单位速度信号下的响应);text(3,23,红线：输入信号；绿线：系统输出响应；蓝线：偏差 );1)原系统闭环脉冲传递函数为：GG =0.36788 亿+0.7183) (zA2 - z + 0.6321)得系统特征根为 z = 0.5 j0.618,特征根 |z|2 =0.52 +0.6182 0.25 +0.4

4、9 =0.74 1, 即闭环系统稳定。2)广义对象的开环脉冲传递函数为：Gd =0.36788 (z+0.7183)(z-1) (z-0.3679)可知，开环脉冲传递函数有z=1的极点个数为1,系统为I型系统。对于单位速度输入，稳态误差不为零。即原来的闭环系统达不到最少拍系统中对某个输入稳定无静差的要 求。3)原闭环系统输出序列为：yd =Columns 1 through 90 0.3679 1.0000 1.3996 1.3996 1.1470 0.8944 0.8015 0.8682Columns 10 through 180.9937 1.0770 1.0810 1.0323 0.98

5、11Columns 19 through 271.0164 1.0070 0.9967 0.9922 0.9943Columns 28 through 350.9994 0.9985 0.9988 0.9998 1.0005 可知，系统输出在第十三拍以后，基本稳定为 达不到最少怕的要求。原闭环系统偏差序列为：ed =0.96070.99921.00070.97261.00281.00030.99751.00330.99991.01481.00151,满足稳定性要求。调节时间为13s。动态过程为13拍,0.1985 0.1318Columns 1 through 91.0000 0.6321 0

6、.0000 -0.3996 -0.3996 -0.1470 0.1056Columns 10 through 180.0063-0.0770-0.0810-0.03230.01890.03930.02740.0025-0.0148Columns 19 through 27-0.0164-0.00700.00330.00780.00570.0008-0.0028-0.0033-0.0015Columns 28 through 350.00060.00150.00120.0002-0.0005-0.0007-0.00030.0001可知，系统偏差输出在第十三拍(系统稳定)以后，基本保持0大小，无静

7、差。闭环系统在单位速度信号下的响应4)原闭环系统输出曲线1.41.210.805101520253035time红线：输入信号；绿线：系统输出响应；蓝线：偏差d 0.6 y n 0.40.20-0.2-0.4从图中也可看出，原系统不满足最少始要求。需要串入最少扣调节器。2、设计最少拍控制器：步骤：(1)确定有零阶保持器的广义对象G(z)clear all ;Close all;%零极点形式的被控对象%采样周期%广义被控对象脉冲传函clc;G=zpk(口,0,-1,1);T=1;Gd=c2d(G,T, zoh);广义被控对象的脉冲传递函数为:0.36788 亿+0.7183)(z-1) (z-0

8、.3679)可知，广义被控对象没有单位圆外的零极点，也没有纯滞后。即 i=j=l=0。(2)输入为单位阶跃时，m=1 , s=j+m=1 , t=l+i=0最少拍系统应具有的误差脉冲传递函数、闭环脉冲传递函数分别为：，(z): m1zm2z-1 - :：J(z) = 1 -z求待定系数，得:1 (z)=1-zzz -1(3)确定数字控制器 D(z)D(z)=：(z)G(zI1 -：，(z)(4)闭环系统输出: 数字控制器的输出：Y(z) = R(z) ：.：，(z) =U (z) G(z)U (z)= E(z) D(z)程序:clear all ;clc;G=zpk(,0,-1,1);%零极点

9、形式的被控对象T=1;Gd=c2d(G,T,zoh)；%采样周期%广义被控对象脉冲传函phi=tf(1,1 0,T); phie=tf(1 -1,1 0,T);D=phi/(Gd*phie); phie,phi,Gd,D num1,den1=tfdata(Gd, num2,den2=tfdata(D,%闭环脉冲传函% 误差脉冲传函%数字控制器脉冲传函v);v);D(z)ud_1=0;ud_2=0;ud_3=0;yd_1=0;yd_2=0;ed_1=0;ed_2=0;ed_3=0;for k=1:1:35time(k)=k*T;rin(k)=1;%单位阶跃信号% 加数字控制器的差分方程%求偏差y

10、d(k)=num1(2)*ud_1+num1(3)*ud_2-den1(2)*yd_1-den1(3)*yd_2;ed(k)=rin(k)-yd(k);%数字控制器的差分方程ud(k)=(num2(1)*ed(k)+num2(2)*ed_1+num2(3)*ed_2+num2(4)*ed_3-den2(2)*ud_1-den2(3)* ud_2-den2(4)*ud_3)/den2(1);%数据更新ud_3=ud_2;ud_2=ud_1;ud_1=ud(k);ed_3=ed_2;ed_2=ed_1;ed_1=ed(k);yd_2=yd_1;yd_1=yd(k);end yd,ud figure

11、(2);plot(time,rin,r ,time,yd, b ,time,ud, k);xlabel( time );););ylabel( rin,yd,ud );title(加数字控制器后单位阶跃信号的闭环响应及数字控制器的输出text(3,1.5,红线：系统输入；蓝线：系统输出；黑线：控制器输出加数字控制器后单位阶跃信号的闭环响应及数字控制器的输出-i !红线：系统输入；蓝线：系统输出；黑线：控制器输出! I： !d u d y-1-2-3 l 0101520time253035系统输出序列：yd =Columns 1 through 100 1.0000 1.0000 1.0000

12、1.0000 1.0000 1.0000 1.0000 1.00001.0000Columns 11 through 201.0000 1.0000 1.0000 1.00001.0000Columns 21 through 301.0000 1.0000 1.0000 1.00001.0000Columns 31 through 351.0000 1.0000 1.0000 1.0000控制器输出序列：ud =Columns 1 through 101.0000 1.0000 1.0000 1.0000 1.00001.0000 1.0000 1.0000 1.0000 1.00001.00

13、002.7183 -2.9525 2.1207 -1.5233 1.0941 -0.7859 0.5645 -0.4055 0.2912-0.2092Columns 11 through 200.1503 -0.1079 0.0775 -0.0557 0.0400 -0.0287 0.0206 -0.0148 0.0106-0.0076Columns 21 through 300.0055 -0.0039 0.0028 -0.0020 0.0015 -0.0011 0.0008 -0.0005 0.0004-0.0003Columns 31 through 350.0002 -0.0001 0

14、.0001 -0.0001 0.0001可知，加数字控制器后，系统输出与系统输入在1s以后相等，说明系统是稳定的且稳态无静差。调节时间为1s,但观察其数字控制器输出曲线及输出序列发现，其控制器输出在28拍(28s)后才变为恒定值,说明输出在第3拍至第28拍间，仅能保证在采样点上与输入相同，在采样点间不相同，即输出有纹波。.11此时闭环脉冲传函为：？(z)=z = zr ”、，一，1 Z - 1误差脉冲传函为：1 -：，(z) =1 -Z =Z数字控制器为：2.7183 z (z-1) (z-0.3679)z 亿+0.7183) (z-1)3、最少拍无纹波控制器：步骤：(1)确定有零阶保持器的广

15、义对象G(z),广义被控对象脉冲传函:0.36788 (z+0.7183) (z-1) (z-0.3679)广义被控对象没有单位圆外的零极点，也没有纯滞后。即 i=j=l=0。(2)由G(z)有一个单位圆内的零点-0.7183,要实现无纹波，则最少拍系统的闭环脉冲传递函数要求将单位 圆内圆外的所有零点都包含进来。此时零点个数i=1。(3)输入为单位阶跃时，m=1 , s=j+m=1 , t=l+i=0+1=1 。应具有的误差脉冲传递函数、闭环脉冲传递函数分别为：中(z) =(1 0.7183z,)mz4 =mz，0.7183mz，_ mz 0.7183m一2z1 一力(z) =(1 z4)(1

16、 fz,)=1 fz -zJ - fz z2 fz -z _ f -2z因1 6(z)十(z) =1，有mz 0.7183mzi 1 fz4-z-fz=1得方程组:m f =1f =0.7183m解方程组，得 m=1/1.7183;f=0.7183/1.7183;即得10.7183z -中 二 1.71831.7183z2210.7183z - z -1.71831.71832z(3)确定数字控制器 D(z)D(z)=中(z)G(zI1-：，(z)(4)闭环系统输出：丫(z)=R(z)(z)=U(z) G(z)数字控制器的输出：U(z)=E(z)，D(z) 程序：clearall;closea

17、ll;%婿极点形式的被控对象%采样周期%广义被控对象脉冲传函cic;G=zpk(,0,-1,1);T=1;Gd=c2d(G,T, zoh);%求误差脉冲传函及闭环脉冲传函phi=tf(1/1.7183 0.7183/1.7183,1 0 0,T);%闭环脉冲传函phie=tf(1 -1/1.7183 -0.7183/1.7183,1 0 0,T);%误差脉冲传函D=phi/(Gd*phie);瀚字控制器脉冲传函D(z)Gd,D,phie,phinum1,den1=tfdata(Gd,v ); %求脉冲传函的分子分母多项式num2,den2=tfdata(D,v);ud_1=0;ud_2=0;u

18、d_3=0;ud_4=0;ud_5=0;yd_1=0;yd_2=0;ed_1=0;ed_2=0;ed_3=0;ed_4=0;ed_5=0;for k=1:1:25time(k)=k*T;rin(k)=1;%单位阶跃信号% 加数字控制器的差分方程yd(k)=num1(2)*ud_1+num1(3)*ud_2-den1(2)*yd_1-den1(3)*yd_2;ed(k)=rin(k)-yd(k);%数字控制器的差分方程ud(k)=(num2(1)*ed(k)+num2(2)*ed_1+num2(3)*ed_2+num2(4)*ed_3+num2(5)*ed_4+num2(6)*ed_5-den2

19、*ud_1-den2(3)*ud_2-den2(4)*ud_3-den2(5)*ud_4-den2(6)*ud_5)/den2；ud 5=ud 4;ud 4=ud 3;ud 3=ud 2;ud 2=ud 1;ud 1=ud(k);ed 5=ed 4;ed 4=ed 3;ed 3=ed 2;ed 2=ed 1;ed 1=ed(k);yd_2=yd_1;yd_1=yd(k);endyd,ud figure(1);plot(time,rin,r ,time,yd, b ,time,ud, k);xlabel( time );ylabel( rin,yd,ud );title(加数字控制器后单位阶跃信

20、号的闭环响应及数字控制器的输出)；text(1.5,1.5,红线：系统输入；蓝线：系统输出；黑线：控制器输出 );系统输出序列：yd =Columns 1 through 10 00.58201.00001.00001.00001.00001.00001.00001.0000Columns 11 through 201.00001.00001.00001.00001.00001.00001.00001.00001.0000Columns 21 through 251.00001.00001.00001.00001.00001.00001.0000控制器输出序列:ud =Columns 1 th

21、rough 101.5820-0.5819-0.00000.0000-0.00000.0000-0.00000.0000-0.00000.0000Columns 11 through 20-0.00000.0000-0.00000.0000-0.00000.0000-0.00000.0000-0.00000.0000由图可以看出，此时调节时间变为2s,调节时间比原来延长1拍，但其控制器输出也在 2拍后变为恒定值,可以保证2拍后，输出不仅在采样点上与输入相同，还在采样点间保证相同，实现输出无纹波。广义被控对象脉冲传函：0.36788 (z+0.7183) (z-1) (z-0.3679) 闭环脉冲传函为:0.582 z + 0.418zA2误差脉冲传函为：zA2 - 0.582 z - 0.418zA2数字控制器为：1.582 zA2 (z+