program_11_1.m

《数字信号处理：基于计算机的方法（第三版）》资源下载第三部分。 著者： （美）米特拉著 作译者： 孙洪等译 ISBN号： 7-121-02652-X 出版日期： 2006-06 出版

% Program 11_1
% Impulse Response Computation
%
n = input('Impulse response length desired = ');
num = input('Numerator coefficients = ');
den = input('Denominator coefficients = ');
x = [1 zeros(1,n-1)];
order = max(length(num),length(den))-1;
si = [zeros(1,order)];
y = filter(num,den,x,si);
stem(y)
xlabel('Time index n'); ylabel('Amplitude');
title('Impulse response samples')