<< Chapter < Page Chapter >> Page >

As mentioned earlier, this design problem might be addressed by general constrained quadratic programming methods [link] , [link] , [link] , [link] , [link] , [link] , [link] , [link] , [link] , [link] .

Examples and observations on cls designs

Here we show that the CLS FIR filter design approach is probably the best general FIR filter design method. For example, a length-31 linear phase lowpass FIRfilter is designed for a band edge of 0.3 and the constraint that the response in the stop cannot be greater than 0.03 is illustrated in [link] .

Figure two is a graph titled Constrained Least Squared Error Design Filter Response. Its horizontal axis is labeled, Normalized Frequency, and it ranges in value from 0 to 2000 in increments of 200. Its vertical axis is labeled Magnitude Response, Log Scale, and ranges in value from 10^-2 to 10^0. There is a horizontal line in the graph, labeled constraint, that runs approximately along the vertical value that is halfway between 10^-1 and 10^-2. A single curve begins at (0, 10^0) and, while wavering, runs horizontally to the right to approximately (500, 10^0), where it sharply decreases with a strengthening slope to (700, 10^-2). After this point, the curve has 11 rounded peaks that are approximately 100 horizontal units in width. The first couple peaks have a height that reaches the horizontal line labeled constraint, and afterwards the heights become smaller as they move to the right.
Response of a Constrained Least Squared Error Filter Design

This filter was designed using the Matlab command: ' fircls1() ' function.

L^p approximation and the iterative reweighted least squares method

We now consider the general L p approximation which contains the least squares L 2 and the Chebyshev L cases. This approach is described in [link] , [link] , [link] using the iterative reweighted least squared (IRLS) alorithm and looks attractive in that it can use different p in different frequency bands. This would allow, for example, a least squared error approximation in thepassband and a Chebyshev approximation in the stopband. The IRLS method can also be used for complex Chebyshev approximations [link] and constrained L 2 approximatin.

Iterative reweighted least squares filter design methods

There are cases where it is desirable to design an FIR filter that will minimize the L p error norm. The error is defined by

q = Ω | A ( ω ) - A d ( ω ) | p d ω

but we usually work with Q 2 . For large p , the results are essentially the same as the Chebyshev filter and this gives a continuum of design between L 2 and L . It also allows the very interesting important possibility of allowing p ( ω ) to be a function of frequency. This means one could have different error criteria in different frequency bands. It can be modifiedto give the same effects as a constraint. This approach is discussed in [link] . It can be applied to complex approximation and to two-dimensional filter design [link] , [link] .

The least squared error and the minimum Chebyshev error criteria are thetwo most commonly used linear-phase FIR filter design methods [link] . There are many situations where better total performance would be obtainedwith a mixture of these two error measures or some compromise design that would give a trade-off between the two. We show how to design a filterwith an L 2 approximation in the passband and a Chebyshev approximation in the stopband. We also show that by formulating the L p problem we can solve the constrained L 2 approximation problem [link] .

This section first explores the minimization of the p t h power of the error as a general approximation method and then shows how this allows L 2 and L approximations to be used simultaneous in different frequency bands of one filter and how the method can be used to imposeconstraints on the frequency response. There are no analytical methods to find this approximation, therefore, an iterative method is used oversamples of the error in the frequency domain. The methods developed here [link] , [link] are based on what is called an iterative reweighted least squared (IRLS) error algorithm [link] , [link] , [link] and they can solve certain FIR filter design problems that neither the Remezexchange algorithm nor analytical L 2 methods can.

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