连续系统的复频分析

实验题目：连续系统的复频分析 实验者 ：霍可栋
实验日期：2011-6-20
1、若某实验的传递函数为 H )(s =
3422+++s s s
试利用拉普拉斯变换来确定：
1）该系统的冲击响应；
2）该系统的阶跃响应； 3）该系统对于输入为)()20cos(t u t u g =的零状态响应；
4）该系统对于输入为)(t u e u t g -=的零状态响应；
syms s t;
Hs=sym('(s+2)/(s^2+4*s+3)');
Vs3=laplace(cos(20*t)); Vs4=laplace(exp(-t));
V os1=Hs; V os2=Hs*(1/s); V os3=Hs*Vs3; V os4=Hs*Vs4; V o1=ilaplace(V os1);
V o2=ilaplace(V os2);
V o3=ilaplace(V os3);
V o4=ilaplace(V os4);
subplot(2,2,1); ezplot(V o1,[0,10]);
subplot(2,2,2); ezplot(V o2,[0,10]);
subplot(2,2,3); ezplot(V o3,[0,10]);
subplot(2,2,4); ezplot(V o4,[0,10]);
title(' 霍可栋 09510441 ');
2、若某实验的传递函数为124253254
32)(12345678235-+-++--+++-+=s s s s s s s s s s s s s H
试确定其零极点，画出零极点点分布图，确定其冲击响应。

clear;
b=[2,0,1,-3,1,4]; a=[5,2,-1,-3,5,2,-4,2,-1];
t=0:0.01:10
sys=tf(b,a);
ps=roots(a); zs=roots(b);
f=step(sys,t);
subplot(2,1,1);
plot(real(zs),imag(zs),'go',real(ps),imag(ps),'mx','markersize',10); grid;
legend('零点','极点');
subplot(2,1,2);
plot(t,f);
title('霍可栋 09510441');。