program_2_6.m

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

% 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');