IIR Digital Filter Design via the Bilinear Transform 1.1 2005/02/03 22:52:26 US/Central 2005/04/24 16:05:17 GMT-5 Douglas L. Jones dl-jones@uiuc.edu Douglas L. Jones dl-jones@uiuc.edu Jeffrey Silverman jsilv@rice.edu Kyle Clarkson kclarks@rice.edu A bilinear transform maps an analog filter Ha s to a discrete-time filter Hz of the same order. If only we could somehow map these optimal analog filter designs to the digital world while preserving the magnitude response characteristics, we could make use of the already-existing body of knowledge concerning optimal analog filter design.

Bilinear Transformation The Bilinear Transform is a nonlinear → mapping that maps a function of the complex variable s to a function of a complex variable z . This map has the property that the LHP in s ( s 0 ) maps to the interior of the unit circle in z , and the λ s axis maps to the unit circle ω in z . Bilinear transform: s α z 1 z 1 H z H a s α z 1 z 1 λ α ω 1 ω 1 α ω 1 ω 1 ω 1 ω 1 2 ω 2 2 ω α ω 2 , so λ α ω 2 , ω 2 λ α . .

The magnitude response doesn't change in the mapping from λ to ω , it is simply warped nonlinearly according to H ω H a α ω 2 , .

The first image implies the second one.

This mapping preserves L errors in (warped) frequency bands. Thus optimal Cauer ( L ) filters in the analog realm can be mapped to L optimal discrete-time IIR filters using the bilinear transform! This is how IIR filters with L optimal magnitude responses are designed. The parameter α provides one degree of freedom which can be used to map a single λ 0 to any desired ω 0 : λ 0 α ω 0 2 or α λ 0 ω 0 2 This can be used, for example, to map the pass-band edge of a lowpass analog prototype filter to any desired pass-band edge in ω . Often, analog prototype filters will be designed with λ 1 as a band edge, and α will be used to locate the band edge in ω . Thus an Mth order optimal lowpass analog filter prototype can be used to design any Mth order discrete-time lowpass IIR filter with the same ripple specifications.

Prewarping Given specifications on the frequency response of an IIR filter to be designed, map these to specifications in the analog frequency domain which are equivalent. Then a satisfactory analog prototype can be designed which, when transformed to discrete-time using the bilinear transformation, will meet the specifications. The goal is to design a high-pass filter, ω s ω s , ω p ω p , δ s δ s , δ p δ p ; pick up some α α 0 . In the δi remain the same and the band edges are mapped by λ i α 0 ω i 2 .

Where λ s α 0 ω s 2 and λ p α 0 ω p 2 .