% Chapter 8: Example 8.12%            Impulse Invariance Transformation%                  Chebyshev-1 Lowpass Filter Design%% Digital Filter Specifications:wp = 0.2*pi;                         % digital Passband freq in Hzws = 0.3*pi;                         % digital Stopband freq in HzRp = 1;                              % Passband ripple in dBAs = 15;                             % Stopband attenuation in dB% Analog Prototype Specifications: Inverse mapping for frequenciesT = 1;                               % Set T=1OmegaP = wp * T;                     % Prototype Passband freqOmegaS = ws * T;                     % Prototype Stopband freqep = sqrt(10^(Rp/10)-1);             % Passband Ripple parameterRipple = sqrt(1/(1+ep*ep));          % Passband RippleAttn = 1/(10^(As/20));               % Stopband Attenuation% Analog Chebyshev-1 Prototype Filter Calculation:[cs,ds] = afd_chb1(OmegaP,OmegaS,Rp,As);%%*** Chebyshev-1 Filter Order =  4 % Impulse Invariance transformation:[b,a] = imp_invr(cs,ds,T);[C,B,A] = dir2par(b,a)%%C =  []%%B =-0.0833   -0.0246%%    0.0833    0.0239%%A = 1.0000   -1.4934    0.8392%%    1.0000   -1.5658    0.6549% Plottingfigure(1); subplot(1,1,1)[db,mag,pha,grd,w] = freqz_m(b,a);subplot(2,2,1); plot(w/pi,mag); title('Magnitude Response')xlabel('frequency in pi units'); ylabel('|H|'); axis([0,1,0,1.1])set(gca,'XTickMode','manual','XTick',[0,0.2,0.3,1]);set(gca,'YTickmode','manual','YTick',[0,Attn,Ripple,1]); gridsubplot(2,2,3); plot(w/pi,db); title('Magnitude in dB');xlabel('frequency in pi units'); ylabel('decibels'); axis([0,1,-40,5]);set(gca,'XTickMode','manual','XTick',[0,0.2,0.3,1]);set(gca,'YTickmode','manual','YTick',[-50,-15,-1,0]); gridset(gca,'YTickLabelMode','manual','YTickLabels',['50';'15';' 1';' 0'])subplot(2,2,2); plot(w/pi,pha/pi); title('Phase Response')xlabel('frequency in pi units'); ylabel('pi units'); axis([0,1,-1,1]);set(gca,'XTickMode','manual','XTick',[0,0.2,0.3,1]);set(gca,'YTickmode','manual','YTick',[-1,0,1]); gridsubplot(2,2,4); plot(w/pi,grd); title('Group Delay')xlabel('frequency in pi units'); ylabel('Samples'); axis([0,1,0,15])set(gca,'XTickMode','manual','XTick',[0,0.2,0.3,1]);set(gca,'YTickmode','manual','YTick',[0:5:15]); grid