02.10.2019 Views

UploadFile_6417

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Characteristics of Prototype Analog Filters 417<br />

value k. Weneed an unnormalized Chebyshev-I filter with arbitrary Ω c .<br />

This is achieved by scaling the array p of the normalized filter by Ω c . Since<br />

this filter has zeros, we also have to scale the array z by Ω c . The new gain<br />

k is determined using (8.56), which is achieved by scaling the old k by the<br />

ratio of the unnormalized to the normalized rational functions evaluated<br />

at s =0.Inthe following function, called U chb2ap(N,As,Omegac), we<br />

design an unnormalized Chebyshev-II analog prototype filter that returns<br />

H a (s) inthe direct form.<br />

function [b,a] = u_chb2ap(N,As,Omegac);<br />

% Unnormalized Chebyshev-2 Analog Lowpass Filter Prototype<br />

% --------------------------------------------------------<br />

% [b,a] = u_chb2ap(N,As,Omegac);<br />

% b = numerator polynomial coefficients<br />

% a = denominator polynomial coefficients<br />

% N = Order of the Elliptic Filter<br />

% As = Stopband Ripple in dB; As > 0<br />

% Omegac = Cutoff frequency in radians/sec<br />

%<br />

[z,p,k] = cheb2ap(N,As);<br />

a = real(poly(p)); aNn = a(N+1);<br />

p = p*Omegac; a = real(poly(p)); aNu = a(N+1);<br />

b = real(poly(z)); M = length(b); bNn = b(M);<br />

z = z*Omegac; b = real(poly(z)); bNu = b(M);<br />

k = k*(aNu*bNn)/(aNn*bNu);<br />

b0 = k; b = k*b;<br />

The design equations for the Chebyshev-II prototype are similar to<br />

those of the Chebyshev-I except that Ω c =Ω s since the ripples are in the<br />

stopband. Therefore we can develop a MATLAB function similar to the<br />

afd chb1 function for the Chebyshev-II prototype.<br />

function [b,a] = afd_chb2(Wp,Ws,Rp,As);<br />

% Analog Lowpass Filter Design: Chebyshev-2<br />

% -----------------------------------------<br />

% [b,a] = afd_chb2(Wp,Ws,Rp,As);<br />

% b = Numerator coefficients of Ha(s)<br />

% a = Denominator coefficients of Ha(s)<br />

% Wp = Passband edge frequency in rad/sec; Wp > 0<br />

% Ws = Stopband edge frequency in rad/sec; Ws > Wp > 0<br />

% Rp = Passband ripple in +dB; (Rp > 0)<br />

% As = Stopband attenuation in +dB; (As > 0)<br />

%<br />

if Wp

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!