ms7_3.m

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

b=[1,0,0];a=[1,-1.5,0.5];
[r p c]=residuez(b,a)
N=20;n=0:N-1;
x=r(1)*p(1).^n +r(2)*p(2).^n
stem(n,x,'k');
title('用部分分式法求反变换x(n)');
axis([0,N-1,1.1*min(x),1.1*max(x)]);