《DSP using MATLAB》示例Example 8.16

%% ------------------------------------------------------------------------
%%            Output Info about this m-file
fprintf(‘\n***********************************************************\n‘);
fprintf(‘        <DSP using MATLAB> Exameple 8.16 \n\n‘);

time_stamp = datestr(now, 31);
[wkd1, wkd2] = weekday(today, ‘long‘);
fprintf(‘      Now is %20s, and it is %8s  \n\n‘, time_stamp, wkd2);
%% ------------------------------------------------------------------------

c = [1, 1];               % Numerator coefficient of Ha(s)
d = [1, 5, 6];            % Denominator coefficient of Ha(s)
T = 1; Fs = 1/T;

[b, a] = bilinear(c, d, Fs)