【MATLAB】【数字信号处理】产生系统的单位冲激响应h(t)与H(z)零极点分布

一、实验目的与要求

产生h(t) 与H(z) 零极点分布

二、实验仪器

微机，仿真软件MATLAB 2022a

三、实验内容与测试结果

1.已知描述连续系统的微分方程为y''(t)+5y'(t)+6y(t)=2x'(t)+8x(t) ，计算系统的单位冲激响应h(t)

程序如下：

clear all;
ts=0;te=5;dt=0.1;
num=[0 2 8];
den=[1 5 6];
sys=tf(num,den);
t=ts:dt:te;
y=impulse(sys,t);
plot(t,y);grid;xlabel('time(sec)');ylabel('h(t)');
hold on;
stem(t,y); %画出离散的和连续的

运行结果如下：

2.已知连续系统的系统函数为Hs=2s2+3s+1 / s3+2s2+2s+1 ，计算其零极点作图.

代码如下：

b=[2,3,2];
a=[1,2,2,1];

ps=roots(a);
zs=roots(b);

% subplot(121);
plot(real(zs),imag(zs),'o',real(ps),imag(ps),'x','markersize',12);
axis([-2,2,-2,2]);
grid;%绘制网格线
legend('零点','极点');
% subplot(122);
% grid;
% pzmap(b,a);
% axis([-2,2,-2,2]);

b=[2,3,2];
a=[1,2,2,1];
ps=roots(a);
zs=roots(b);
plot(real(zs),imag(zs),'o',real(ps),imag(ps),'x','markersize',12);
axis([-2,2,-2,2]);
grid;%绘制网格线
legend('零点','极点');

运行结果如下：