自动控制原理实验7

实验1：已知系统的闭环传递函数为G（s）=1/(s2+0.4s+1),试求其单位阶跃和单位斜坡响应曲线。

num1=[1];

den1=[1 0.4 1];

sys=tf(num1,den1);

t=0:0.01:4;

step(sys);

grid;

figure;

u=t;

lsim(sys,u,t,0);

grid;

实验2：

用游动鼠标法性能指标。

sys=zpk([],[-1+3*i,-1-3*i],3); figure;

[y,t]=step(sys);

plot(t,y);

grid;

c=dcgain(sys)

编程法求取系统的性能指标

sys=zpk([],[-1+3*i,-1-3*i],3); [y,t]=step(sys);

[Y,k]=max(y);

timeopeak=t(k)

c=dcgain(sys)

w=(Y-c)/c

n=1;

while y(n)

n=n+1;

end

risetime=t(n)

实验3：已知单位负反馈的二阶系统的开环传递函数G（s）=wn2/s(s+2wn§),试绘制§分别为0,0.2,0.4,0.6,0.9,1.2,1.5时，其单位负反馈系统的单位阶跃响应曲线（绘制在同一张图上）。并分析§（sigma）为多少的时候系统响应为最佳。

sigma=[0,0.2,0.4,0.6,0.9,1.2,1.5];

t=0:0.2:20

for j=1:7

num=[1];

den=[1 2*sigma 0];

s1=tf(num,den);

sys=feedback(s1,1);

Y( : ,j)=step(sys,t);

end

plot (t ,Y( : ,1:7));

grid;

gtext(‘sigma=0’);

gtext(‘sigma=0.2’)

gtext(‘sigma=0.4’)

gtext(‘sigma=0.6’)

gtext(‘sigma=0.9’)

gtext(‘sigma=1.2’)

gtext(‘sigma=1.5’)

由图所知，当§=0.6时候，系统响应为最佳。