<< Chapter < Page Chapter >> Page >

Given a desired frequency response, the frequency sampling design method designs a filter with a frequency response exactly equal to the desired response at a particular set of frequencies k .

Procedure

k k
    o 1 N 1
H d k n M 1 0 h n k n
Desired Response must incluce linear phase shift (if linear phase is desired)

What is H d for an ideal lowpass filter, cotoff at c ?

M 1 2 c c 0 c c

Got questions? Get instant answers now!
This set of linear equations can be written in matrix form
H d k n M 1 0 h n k n
H d 0 H d 1 H d N - 1 0 0 0 1 0 M 1 1 0 1 1 1 M 1 M - 1 0 M - 1 1 M - 1 M 1 h 0 h 1 h M 1

or H d W h So

h W H d
W is a square matrix for N M , and invertible as long as i j 2 l , i j

Important special case

What if the frequencies are equally spaced between 0 and 2 , i.e. k 2 k M

Then H d k n M 1 0 h n 2 k n M n n M 1 0 h n n 2 k n M DFT! so h n n 1 M k M 1 0 H d k 2 n k M or h n n M k M 1 0 H d k 2 n k M n IDFT H d k

Important special case #2

h n symmetric, linear phase, and has real coefficients. Since h n h M n 1 , there are only M 2 degrees of freedom, and only M 2 linear equations are required.

H k n M 1 0 h n k n n M 2 1 0 h n k n k M n 1 M even n M 3 2 0 h n k n k M n 1 h M 1 2 k M 1 2 M odd k M 1 2 2 n M 2 1 0 h n k M 1 2 n M even k M 1 2 2 n M 3 2 0 h n k M 1 2 n h M 1 2 M odd

Removing linear phase from both sides yields A k 2 n M 2 1 0 h n k M 1 2 n M even 2 n M 3 2 0 h n k M 1 2 n h M 1 2 M odd Due to symmetry of response for real coefficients, only M 2 k on 0 need be specified, with the frequencies k thereby being implicitly defined also. Thus we have M 2 real-valued simultaneous linear equations to solve for h n .

Special case 2a

h n symmetric, odd length, linear phase, real coefficients, and k equally spaced: k 0 k M 1 k n k M

h n IDFT H d k 1 M k M 1 0 A k 2 k M M 1 2 2 n k M 1 M k M 1 0 A k 2 k M n M 1 2

To yield real coefficients, A mus be symmetric A A A k A M k

h n 1 M A 0 k M 1 2 1 A k 2 k M n M 1 2 2 k n M 1 2 1 M A 0 2 k M 1 2 1 A k 2 k M n M 1 2 1 M A 0 2 k M 1 2 1 A k 1 k 2 k M n 1 2

Simlar equations exist for even lengths, anti-symmetric, and 1 2 filter forms.

Comments on frequency-sampled design

This method is simple conceptually and very efficient for equally spaced samples, since h n can be computed using the IDFT.

H for a frequency sampled design goes exactly through the sample points, but it may be very far off from the desired response for k . This is the main problem with frequency sampled design.

Possible solution to this problem: specify more frequency samples than degrees of freedom, and minimize the total errorin the frequency response at all of these samples.

Extended frequency sample design

For the samples H k where 0 k M 1 and N M , find h n , where 0 n M 1 minimizing H d k H k

For l norm, this becomes a linear programming problem (standard packages availble!)

Here we will consider the 2 l norm.

To minimize the 2 l norm; that is, n N 1 0 H d k H k , we have an overdetermined set of linear equations: 0 0 0 M 1 N - 1 0 N - 1 M 1 h H d 0 H d 1 H d N - 1 or W h H d

The minimum error norm solution is well known to be h W W W H d ; W W W is well known as the pseudo-inverse matrix.

Extended frequency sampled design discourages radical behavior of the frequency response between samples forsufficiently closely spaced samples. However, the actual frequency response may no longer pass exactly through any of the H d k .

Get Jobilize Job Search Mobile App in your pocket Now!

Get it on Google Play Download on the App Store Now




Source:  OpenStax, Digital filter design. OpenStax CNX. Jun 09, 2005 Download for free at http://cnx.org/content/col10285/1.1
Google Play and the Google Play logo are trademarks of Google Inc.

Notification Switch

Would you like to follow the 'Digital filter design' conversation and receive update notifications?

Ask