%% Frequency Response of a Digital Bessel Filter
% Design an analog Bessel filter of order 5. Convert it to a digital IIR filter using bilinear. Display its frequency response.

Fs = 100;                           % Sampling Frequency
[z,p,k] = besself(5,1000);          % Bessel analog filter design
[zd,pd,kd] = bilinear(z,p,k,Fs);    % Analog to digital mapping
sos = zp2sos(zd,pd,kd);             % Convert to SOS form
fvtool(sos)                         % Visualize the digital filter

% Copyright 2015 The MathWorks, Inc.