ms4_6.m

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

a=[1, 0, 0.34319, 0, 0.60439, 0, 0.20407];
b=[1, 0, -3, 0, 3, 0, - 1];
[r p k]=residue(b,a);        %求时域冲激响应
nr=length(r);
for n=1:0.01:nr;
   hn(n)=r(n)*(p(n)^n);
end
stem(hn,'k');                           %显示冲激响应曲线
title('系统的单位冲激响应');
ylabel('h(n)');xlabel('n');