连续系统分析
一、已知描述连续系统的微分方程,计算该系统的响应并与理论结果比较
an=[1,10];
bn=[2];
sys=tf(bn,an);
t=0:0.01:1.5;
xt=stepfun(t,0);
yzs=lsim(sys,xt,t)';
figure(1);
subplot(2,2,1);
plot(t,yzs);
title('yzs')
yzi=exp(-10*t);
subplot(2,2,2);
plot(t,yzi);
title('yzi');
ycom=yzs+yzi;
subplot(2,1,2);
plot(t,ycom);
hold on;
yzstheory=1/5*(xt-exp(-10*t));
ytheory=yzstheory+yzi;
plot(t,ytheory,'r.');
hold off;
title('y');
axis([0,1.5,0,1]);
legend('计算值','理论值');
二、研究具有以下零极点的连续系统
(1) 利用zpk和tf命令建立系统的系统函数,画出系统的零极点图。
(2) 分析系统是否稳定。若稳定,画出系统的幅频特性曲线。
(3) 画出系统的冲激响应波形。
(4) 详细列出根据零极点分析系统特性的过程。
2.1 1个极点s=-0.1,增益k=1
a=[1,0.1];
b=[1];
sys=tf(b,a);
figure(1);
subplot(2,1,1);
pzmap(sys);
title('(a)')
w=0:0.001:2*2*pi;
H=freqs(b,a, w);
subplot(2,2,3);
plot(w/2/pi, abs(H));
title