ms7r41.m

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

b=[10,20,0,0];a=[1,8,19,12];
[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+r(3)*p(3).^n;
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)');