ms7_4.m

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

b=[0,1];a=[1,-12,36];
[r p c]=residuez(b,a)
r=real(r);p=real(p);
N=8;n=0:N-1;
%h=(r(1)+r(2))*p(1).^n +r(2)*n.*p(2).^n;
h=r(1)*p(1).^n.*[n>=0]+r(2).*(n+1).*p(2).^n.*[n-1>=0];
subplot(1,2,1),stem(n,h,'filled','k');
axis([0,N,0,3.5e5]);
title('用部分分式法求反变换h(n)');
h2=impz(b,a,N);
subplot(1,2,2),stem(n,h2,'filled','k')
title('用impz求反变换h(n)');