ms7r42.m

《数字信号处理实验(MATLAB版)》这

b=[0,0,5];a=[1,1,-0.6];
[r p c]=residuez(b,a)
r=real(r);p=real(p);
N=10;n=0:N-1;
x=r(1)*p(1).^n+r(2)*p(2).^n+c(1).*[n==0];
subplot(1,2,1),stem(n,x,'filled','k');
title('用部分分式法求反变换x(n)');
x2=impz(b,a,N);
subplot(1,2,2),stem(n,x2,'filled','k')
title('用impz求反变换x(n)');