数字信号处理（2）

由于时间太久，我找不到实验报告的要求，而且前几次实验比较简单，直接发出来大家对比一下

%1.求单位取样响应并判断稳定性 参数a,b,a1,b1
a=[3 4 1];
b=[1 1];
subplot(4,1,1)
impz(b,a,100)
grid on;
title('单位取样函数1稳定')

a1=[2.5 6 10];
b1=[1];
subplot(4,1,2)
impz(b1,a1,100)
grid on;
title('单位取样函数2不稳定')

运行结果：

%2.求零状态响应 参数nx nh x h y ny1 ny2 ny
nx=-1:5;
nh=-1:10;
x=uDT(nx)-uDT(nx-5);
h=(7/8).^nh.*(uDT(nh)-uDT(nh-10));
ny1=nx(1)+nh(1);
ny2=nx(end)+nh(end);
ny=ny1:ny2;
subplot(4,1,3)
y=conv(h,x);
stem(ny,y);
grid on;
title('零状态响应')

运行结果：