example6_11.m

西交大刘树棠写的《数字信号处理》那本书的所有源代码

Wp=0.2*pi;Ws=0.3*pi;Rp=1;As=16;
Ripple=10^(-Rp/20);Attn=10^(-As/20);
%Analog filter design:
[b,a]=afd_chb1(Wp,Ws,Rp,As)
%Chebyshev-1 Filter Order=4
% Calculation of second-order sections:
[C,B,A]=sdir2cas(b,a)
%Calculation of Frequency Response:
[db,mag,pha,w]=freqs_m(b,a,0.5*pi);
%Calculation of Impulse response:
%[ha,x,t]=impulse(b,a);