extra.m

来自「该程序用于讨论各种离散时间信号和系统的性质及特性。」· M 代码 · 共 14 行

% Program P2_5
% Compute the impulse response y
clf;
x = [1 ones(1,40)]; % Generate the input
n = 0:41;
% Compute the output of 4th order system
num = [2.2403 2.4908 2.2403];
den = [1 -0.4 0.75];
y = filter(num,den,x);
% Plot the impulse response
stem(y);
xlabel('Time index n'); ylabel('Amplitude');
title('Impulse Response'); grid;