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.

466 Chapter 8 IIR FILTER DESIGN<br />

P8.24 In this problem we will develop a technique called the step invariance transformation. In<br />

this technique, the step response of an analog prototype filter is preserved in the resulting<br />

digital filter; i.e., if v a(t) isthe step response of the prototype and if v(n) isthe step<br />

response of the digital filter, then<br />

v(n) = v a(t)| t=nT<br />

,<br />

T : sampling interval<br />

Note that the frequency-domain quantities are related by<br />

V a(s) △ = L [v a(t)] = H a(s)/s<br />

and<br />

V (z) = △ 1<br />

Z [v(n)] = H(z)<br />

1 − z −1<br />

Hence the step invariance transformation steps are: Given H a(s)<br />

• divide H a(s) bys to obtain V a(s),<br />

• find residues {R k } and poles {p k } of V a(s),<br />

• transform analog poles {p k } into digital poles { e p kT } where T is arbitrary,<br />

• determine V (z) from residues {R k } and poles { e p kT } , and finally<br />

• determine H(z) bymultiplying V (z) by ( 1 − z −1) .<br />

Use the above procedure to develop a MATLAB function to implement the step invariance<br />

transformation. The format of this function should be<br />

function [b,a] =stp_invr(c,d,T)<br />

% Step Invariance Transformation from Analog to Digital Filter<br />

% [b,a] =stp_invr(c,d,T)<br />

% b = Numerator polynomial in z^(-1) of the digital filter<br />

% a = Denominator polynomial in z^(-1) of the digital filter<br />

% c = Numerator polynomial in s of the analog filter<br />

% d = Denominator polynomial in s of the analog filter<br />

% T = Sampling (transformation) parameter<br />

P8.25 Design the lowpass Butterworth digital filter of Problem P8.22 using the step invariance<br />

method. Plot the log-magnitude response in dB and compare it with that in Problem P8.22.<br />

Plot the step response v(n) and the impulse response v a(t) ofthe analog prototype and<br />

compare their shapes.<br />

P8.26 In this chapter we discussed a filter transformation technique called the matched-z<br />

transformation. Using (8.69) write a MATLAB function called mzt that maps the analog<br />

system function H a(s) into the digital system function H(z). The format of the function<br />

should be<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!