实验三、最少拍控制算法matlab仿真实验2017

实验三、最少拍控制算法 matlab 仿真实验

实验目的： 1．掌握最少拍有纹波、无纹波系统的设计方法； 2．学会最少拍控制系统的分析方法； 3．了解输入信号对最少拍控制系统的影响及改进措施

二、实验内容： 设单位反馈线性定常系统的连续部分和零阶保持器的传递函数分别为：

1 1 e

G 0(s)

，

G h (s)

s(s 1) s

，采样周期 T=1 秒。要求系统在单位阶跃输入时实现最少拍控制，求数字控制器的脉冲传递函数 D (z)。讨 论加上控制器后的系统输出在采样点之间是否存在纹波，若存在纹波，请设计最少拍无纹波控制器。系统 结构如图：

1、对未加控制器时的闭环系统分析其是否满足最少拍要求

clear all ; clc;

G=zpk([],[0,-1],1); T=1;

Gd=c2d(G,T, 'zoh' ) [num1,den1]=tfdata(Gd, GG=feedback(Gd,1) yd_1=0;yd_2=0; ed_1=0;ed_2=0; for k=1:1:35 time(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); %求

偏差 %数据更新

ed_2=ed_1;ed_1=ed(k); yd_2=yd_1;yd_1=yd(k); end

yd %原闭环系统的输出序列 ed %原闭环系统的偏差序列 figure(1);

plot(time,rin, 'r' ,time,yd, 'g' ,time,ed, 'b' ); xlabel( 'time' ); ylabel( 'rin,yd' );

title( ' 闭环系统在单位速度信号下的响应 ' );

%零极点形式的被控

对象

%采样周期

%广义被控对象脉冲传函

'v' );

text(3,23, ' 红线：输入信号；绿线：系统输出响应；蓝线：偏差' );

1)原系统闭环脉冲传递函数为：

GG =

0.36788 (z+0.7183) (z^2 - z + 0.6321)

得系统特征根为z 0.5 j0.618 ，

特征根z 20.520.61820.25 0.49 0.74 1，即闭环系统稳定。

2)广义对象的开环脉冲传递函数为：

Gd =

0.36788 (z+0.7183)

(z-1) (z-0.3679) 可知，开环脉冲传递函数有z 1 的极点个数为 1，系统为 I 型系统。对于单位速度输入，稳态误差不为零。即原来的闭环系统达不到最少拍系统中对某个输入稳定无静差的要求。

3)原闭环系统输出序列为：

yd =

Columns 1 through 9

0 0.3679 1.0000 1.3996 1.3996 1.1470 0.8944 0.8015 0.8682

Columns 10 through 18

0.9937 1.0770 1.0810 1.0323 0.9811 0.9607 0.9726 0.9975 1.0148

Columns 19 through 27

1.0164 1.0070 0.9967 0.9922 0.9943 0.9992 1.0028 1.0033 1.0015

Columns 28 through 35

0.9994 0.9985 0.9988 0.9998 1.0005 1.0007 1.0003 0.9999

可知，系统输出在第十三拍以后，基本稳定为1，满足稳定性要求。调节时间为13s。动态过程为 13 拍，达不到最少怕的要求。

原闭环系统偏差序列为：

ed =

Columns 1 through 9

1.0000 0.6321 0.0000 -0.3996 -0.3996 -0.1470 0.1056 0.1985 0.1318

Columns 10 through 18

0.0063 -0.0770 -0.0810 -0.0323 0.0189 0.0393 0.0274 0.0025 -0.0148

Columns 19 through 27

-0.0164 -0.0070 0.0033 0.0078

0.0057 0.0008 -0.0028 -0.0033 -0.0015

Columns 28 through 35

0.0006 0.0015 0.0012 0.0002 -0.0005 -0.0007 -0.0003 0.0001 可知，系统偏差输出在第十三拍(系统稳定)以后，基本保持 0 大小，无静差。

4) 原闭环系统输出曲线

1.4 1.2 1 0.8 y

d 0.6 ,y i ,r n

0.4

0.2

0 -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 (z+0.7183)

%零极点形式的被控对象 %采样周期

%广义被控对象脉冲

time

闭环系统在单位速度信号下的响应