一、实验目的与要求
产生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('零点','极点');
运行结果如下: