program_2_6.m

这是国外一本经典教材附带的程序

% Program 2_6
% Illustration of Impulse Response Computation
%
N = input('Desired impulse response length = ');
p = input('Type in the vector p =');
d = input('Type in the vector d =');
x = [1 zeros(1,N-1)];
y = filter(p,d,x);
k = 0:1:N-1;
stem(k,y)
xlabel('Time index n'); ylabel('Amplitude');