02.10.2019 Views

UploadFile_6417

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Frequency-band Transformations 455<br />

Now in practice we will know the desired highpass frequencies ω s and<br />

ω p , and we are required to find the prototype lowpass cutoff frequencies<br />

ω ′ s and ω ′ p.Wecan choose the passband frequency ω ′ p with a reasonable<br />

value, say ω ′ p =0.2π, and determine α from ω p using the formula from<br />

Table 8.2. Now ω ′ s can be determined (for our highpass filter example)<br />

from α and<br />

Z = − z−1 + α<br />

1+αz −1<br />

where Z = e jω′ s<br />

and z = e<br />

jω s<br />

,or<br />

ω ′ s = ̸<br />

(<br />

)<br />

− e−jωs + α<br />

1+αe −jωs<br />

(8.71)<br />

Continuing our highpass filter example, let ω p =0.6π and ω s =0.4586π be<br />

the band-edge frequencies. Let us choose ω ′ p =0.2π. Then α = −0.38197<br />

from (8.70), and from (8.71)<br />

ω ′ s = ̸<br />

(<br />

)<br />

− e−j0.4586π − 0.38197<br />

1 − 0.38197e −j−0.38197 =0.3π<br />

as expected. Now we can design a digital lowpass filter and transform<br />

it into a highpass filter using the zmapping function to complete our<br />

design procedure. For designing a highpass Chebyshev-I digital filter, the<br />

above procedure can be incorporated into a MATLAB function called the<br />

cheb1hpf function shown here.<br />

function [b,a] = cheb1hpf(wp,ws,Rp,As)<br />

% IIR Highpass filter design using Chebyshev-1 prototype<br />

% function [b,a] = cheb1hpf(wp,ws,Rp,As)<br />

% b = Numerator polynomial of the highpass filter<br />

% a = Denominator polynomial of the highpass filter<br />

% wp = Passband frequency in radians<br />

% ws = Stopband frequency in radians<br />

% Rp = Passband ripple in dB<br />

% As = Stopband attenuation in dB<br />

%<br />

% Determine the digital lowpass cutoff frequencies:<br />

wplp = 0.2*pi;<br />

alpha = -(cos((wplp+wp)/2))/(cos((wplp-wp)/2));<br />

wslp = angle(-(exp(-j*ws)+alpha)/(1+alpha*exp(-j*ws)));<br />

%<br />

Copyright 2010 Cengage Learning. All Rights Reserved. May not be copied, scanned, or duplicated, in whole or in part. Due to electronic rights, some third party content may be suppressed from the eBook and/or eChapter(s).<br />

Editorial review has deemed that any suppressed content does not materially affect the overall learning experience. Cengage Learning reserves the right to remove additional content at any time if subsequent rights restrictions require it.

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

Saved successfully!

Ooh no, something went wrong!