ms7_5.m

来自「《数字信号处理实验(MATLAB版)》这」· M 代码 · 共 10 行

a=[1, 0, 0.34319, 0, 0.60439, 0, 0.20407];
b=[0.1321, 0, -0.3963, 0, 0.3963, 0, -0.1321];
[r p c]=residuez(b,a)
N=40;n=0:N-1;
h=r(1)*p(1).^n+r(2)*p(2).^n+r(3)*p(3).^n+r(4)*p(4).^n+r(5)*p(5).^n+r(6)*p(6).^n+c(1).*[n==0];
subplot(1,2,1),stem(n,real(h),'k');
title('用部分分式法求反变换h(n)');
h2=impz(b,a,N);
subplot(1,2,2),stem(n,h2,'k');
title('用impz求反变换h(n)');