实验三 连续系统分析

实验三 连续系统分析

实验目的

 深刻理解连续时间系统的系统函数在分析连续系统的时域特性、频域特性及稳定性中的重要作用及意义，掌握根据系统函数的零极点设计简单的滤波器的方法。掌握利用MATLAB分析连续系统的时域响应、频响特性和零极点的基本方法。

实验内容

已知连续系统的系统函数为$y^{\prime \prime }(t)+3y^{\prime }(t)+2y(t)=4x^{\prime }(t)+x(t)$

1. 计算并画出该系统的单位冲激响应h(t),单位阶跃响应g(t)；

求LTI系统的冲激响应的函数为：impulse(b,a)
求LTI系统的阶跃响应的函数为：step(b,a)

t=0:0.01:5;
sys=tf([4 1],[1 3 2]);
y1=impulse(sys,t);
figure(1);plot(t,y1);
xlabel('Time(sec)');
ylabel('h(t)');
y2=step(sys,t);
figure(2);plot(t,y2);
xlabel('Time(sec)');
ylabel('g(t)');

另一种方式
a=[1 3 2];%构造系数向量
b=[4 1];
subplot(2,1,1)
impulse(b,a);
subplot(2,1,2)
step(b,a);

2. 求输入分别为u(t),sintu(t),e^(-t)u(t) 时，计算并画出系统的零状态输出响应y(t)；

lsim()表示输入函数的响应

sys=tf([4 1],[1 3 2]);
t1=0:0.01:20;
t2=0:0.01:50;
t0=0;
f1=stepfun(t1,t0);
y1=lsim(sys,f1,t1);
figure(1);plot(t1,y1);
xlabel('Time(sec)');
ylabel('y1(t)');
f2=sin(t2.*stepfun(t2,t0));
y2=lsim(sys,f2,t2);
figure(2);plot(t2,y2);
xlabel('Time(sec)');
ylabel('y2(t)');
f3=exp(-t1).*stepfun(t1,t0);
y3=lsim(sys,f3,t1);
figure(3);plot(t1,y3);
xlabel('Time(sec)');
ylabel('y3(t)');

3. 计算该系统的零极点，并作出零极点分布图；

num=[4 1];
den=[1 3 2];
sys=tf(num,den);
poles=roots(den);
pzmap(sys);

4. 分析系统的幅频特性

num=[4 1];
den=[1 3 2];
sys=tf(num,den);
t=0:0.02:10;
[H,w]=freqs(num,den);
figure(3);plot(w,abs(H));