<< Chapter < Page Chapter >> Page >

To develop the theory for linear phase FIR filters, a careful definition of phase shift is necessary. If the real and imaginary parts of H ( ω ) are given by

H ( ω ) = R ( ω ) + j I ( ω )

where j = - 1 and the magnitude is defined by

| H ( ω ) | = R 2 + I 2

and the phase by

Φ ( ω ) = arctan ( I / R )

which gives

H ( ω ) = | H ( ω ) | e j Φ ( ω )

in terms of the magnitude and phase. Using the real and imaginary parts is using a rectangular coordinate system and using the magnitude and phase isusing a polar coordinate system. Often, the polar system is easier to interpret.

Mathematical problems arise from using | H ( ω ) | and Φ ( ω ) , because | H ( ω ) | is not analytic and Φ ( ω ) not continuous. This problem is solved by introducing an amplitude function A ( ω ) that is real valued and may be positive or negative. The frequency response is written as

H ( ω ) = A ( ω ) e j Θ ( ω )

where A ( ω ) is called the amplitude in order to distinguish it from the magnitude | H ( ω ) | , and Θ ( ω ) is the continuous version of Φ ( ω ) . A ( ω ) is a real, analytic function that is related to the magnitude by

A ( ω ) = ± | H ( ω ) |

or

| A ( ω ) | = | H ( ω ) |

With this definition, A ( ω ) can be made analytic and Θ ( ω ) continuous. These are much easier to work with than | H ( ω ) | and Φ ( ω ) . The relationship of A ( ω ) and | H ( ω ) | , and of Θ ( ω ) and Φ ( ω ) are shown in [link] .

This image consist of four graphs. Above each graph is the label FIR Filter Frequency Response. The top left graph compares the normalized frequency on the x axis with the magnitude response, M on the y axis. The line for this graph follows a negative slope from point (0,14) to about (.5,0). The line then makes a shallow arch from that point to about (.75.0). The line then follows a positive slope as it runs off the graph. The top right graph is similar except that y axis is labeled Amplitude response, A. The line starts at (14,0) and follows a negative slope falls below the x axis and the then the slope turns positive and proceeds off the graph. The bottom left graph has a y axis labeled Phase in Radians, F. This line is very jagged and consist only of straight lines. The line begins with a negative slope from the origin to (.5,-3). Then the line proceeds vertically to point(.5,0). Then the line takes on the same negative slope as the first line and then goes vertical again crossing over the x axis to a point around(.75,1), and then the line takes on the same negative slope again and ends at (1,0). The bottom right graph is simpler. A single straight line follows a negative slope from the origin to about (1,-6) where the line ends. The y axis in this graph is labeled Phase in Radians Q.
Magnitude and Amplitude Frequency Responses and Corresponding Phase Frequency Response of Example Filter

To develop the characteristics and properties of linear-phase filters, assume a general linear plus constant form for the phasefunction as

Θ ( ω ) = K 1 + K 2 ω

This gives the frequency response function of a length-N FIR filter as

H ( ω ) = n = 0 N - 1 h ( n ) e - j ω n = e - j ω M n = 0 N - 1 h ( n ) e j ω ( M - n )

and

H ( ω ) = e - j ω M [ h 0 e j ω M + h 1 e j ω ( M - 1 ) + + h N - 1 e j ω ( M - N + 1 ) ]

[link] can be put in the form of

H ( ω ) = A ( ω ) e j ( K 1 + K 2 ω )

if M (not necessarily an integer) is defined by

M = N - 1 2

or equivalently,

M = N - M - 1

[link] then becomes

H ( ω ) = e - j ω M [ ( h 0 + h N - 1 ) cos ( ω M ) + j ( h 0 - h N - 1 ) sin ( ω M ) + ( h 1 + h N - 2 ) cos ( ω ( M - 1 ) ) + j ( h 1 - h N - 2 ) sin ( w ( M - 1 ) ) + ]

There are two possibilities for putting this in the form of [link] where A ( ω ) is real: K 1 = 0 or K 1 = π / 2 . The first case requires a special even symmetry in h ( n ) of the form

h ( n ) = h ( N - n - 1 )

which gives

H ( ω ) = A ( ω ) e - j M ω

where A ( ω ) is the amplitude, a real-valued function of ω and e - j M ω gives the linear phase with M being the group delay. For the case where N is odd, using [link] , [link] , and [link] , we have

A ( ω ) = n = 0 M - 1 2 h ( n ) cos ω ( M - n ) + h ( M )

or with a change of variables,

A ( ω ) = n = 1 M 2 h ( M - n ) cos ( ω n ) + h ( M )

which becomes

A ( ω ) = n = 1 M 2 h ^ ( n ) cos ( ω n ) + h ( M )

where h ^ ( n ) = h ( M - n ) is a shifted h ( n ) . These formulas can be made simpler by defining new coefficients so that [link] becomes

A ( ω ) = n = 0 M a ( n ) cos ( ω ( M - n ) )

where

a ( n ) = { 2 h ( n ) for 0 n M - 1 h ( M ) for n = M 0 otherwise

and [link] becomes

A ( ω ) = n = 0 M a ( n ) cos ( ω n )

with

a ( n ) = { h ( M ) for n = 0 2 h ( M + n ) for 1 n M 0 otherwise .

Notice from [link] for N odd, A ( ω ) is an even function around ω = 0 and ω = π , and is periodic with period 2 π .

For the case where N is even,

A ( ω ) = n = 0 N / 2 - 1 2 h ( n ) cos ω ( M - n )

or with a change of variables,

A ( ω ) = n = 1 N / 2 2 h ( N / 2 - n ) cos ω ( n - 1 / 2 )

These formulas can also be made simpler by defining new coefficients so that [link] becomes

A ( ω ) = n = 0 N / 2 - 1 a ( n ) cos ( ω ( M - n ) )

where

a ( n ) = { 2 h ( n ) for 0 n N / 2 - 1 0 otherwise

and [link] becomes

A ( ω ) = n = 1 N / 2 a ( n ) cos ( ω ( n - 1 / 2 ) )

with

a ( n ) = { 2 h ( N / 2 - n ) for 1 n N / 2 0 otherwise

Get Jobilize Job Search Mobile App in your pocket Now!

Get it on Google Play Download on the App Store Now




Source:  OpenStax, Digital signal processing and digital filter design (draft). OpenStax CNX. Nov 17, 2012 Download for free at http://cnx.org/content/col10598/1.6
Google Play and the Google Play logo are trademarks of Google Inc.

Notification Switch

Would you like to follow the 'Digital signal processing and digital filter design (draft)' conversation and receive update notifications?

Ask