program_8_1.m

来自「这是国外一本经典教材附带的程序」· M 代码 · 共 14 行

% Program 8_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(n-1,y)
xlabel('Time index n'); ylabel('Amplitude');
title('Impulse response samples');