3차 저역 통과 필터가 다음 차분 방정식으로 주어질 때,

이 필터의 크기와 응답을 그리고, 저역 통과 필터임을 확인해보자.

ex3_16.cem

b = [0.0181,  0.0543, 0.0543,  0.0181];
a = [1.0000, -1.7600, 1.1829, -0.2781];
m = 0:length(b)-1; l = 0:length(a)-1;
K = 500; k = 1:K:1;
w = pi*k/K;  
% [0, pi] axis divided into 501 points.
num = b * exp(-j*m'*w); 
% Numerator calculations
den = a * exp(-j*l'*w); 
% Denominator calculations
H = num ./ den;

magH = abs(H); angH = angle(H);
subplot(1,1,1);
subplot(2,1,1); plot(w/pi,magH); 
xtitle("frequency in pi units"); ytitle("|H|");
title("Magnitude Response");
subplot(2,1,2); plot(w/pi,angH/pi); 
xtitle("frequency in pi units"); ytitle("Phase in pi Radians");
title("Phase Response");