2.2 Least squared error design of fir filters  (Page 5/13)

where $A_d\left(\omega \right)$ is the desired ideal amplitude response, $A\left(\omega \right)={\sum }_{n}a\left(n\right)cos\left(\omega (M-N)n\right)$ is the achieved amplitude response with the length $h\left(n\right)$ related to $h\left(n\right)$ by Equation 29 from FIR Digital Filters . This integral squared error is approximated by the discrete squared errordefined in [link] for $L>>N$ which in some cases is much easier to minimize. However for some very useful cases, formulas can be foundfor $h\left(n\right)$ that minimize [link] and that is what we will be considering in this section.

The unweighted least integral squared error approximation

If the error measure is the unweighted integral squared error defined in [link] , Parseval's theorem gives the equivalent time-domain formulation for the error to be

In general, this ideal response is infinite in duration and, therefore, cannot be realized exactly byan actual FIR filter.

As was done in the case of the discrete error measure, we break the infinite sum in [link] into two parts, one of which depends on $h\left(n\right)$ and the other does not.

Again, we see that the minimum $q$ is achieved by using $h\left(n\right)=h_d\left(n\right)$ for $-M\le n\le M$ . In other words, the infinitely long $h_d\left(n\right)$ is symmetrically truncated to give the optimal least integral squared errorapproximation. The problem then becomes one of finding the $h_d\left(n\right)$ to truncate.

Here the integral definition of approximation error is used. This is usually what we really want, but in some cases the integrals can not becarried out and the sampled method above must be used.

Ideal constant gain passband lowpass filter

Here we assume the simplest ideal lowpass single band FIR filter to have unity passband gain for $0<\omega <{\omega }_0$ and zero stopband gain for ${\omega }_0<\omega <\pi$ similar to those in Figure 8a from FIR Digital Filters and [link] . This gives

as the ideal desired amplitude response. The ideal shifted filter coefficients are the inverse DTFT from Equation 15 from Chebyshev or Equal Ripple Error Approximation Filters of this amplitude which for $N$ odd are given by

which is sometimes called a “sinc" function. Note ${\widehat{h}}_d\left(n\right)$ is generally infinite in length. This is now symmetricallytruncated and shifted by $M=(N-1)/2$ to give the optimal, causal length- $N$ FIR filter coefficients as

and $h\left(n\right)=0$ otherwise. The corresponding derivation for an even length starts with the inverseDTFT in Equation 5 from Constrained Approximation and Mixed Criteria for a shifted even length filter is

which when truncated and shifted by $N/2$ gives the same formula as for the odd length design in [link] but one should note that $M=(N-1)/2$ is not an integer for an even $N$ .

Ideal linearly increasing gain passband lowpass filter

We now derive the design formula for a filter with an ideal amplitude response that is a linearly increasing function in the passband ratherthan a constant as was assumed above. This ideal amplitude response is given byand illustrated in [link] For $N$ odd, the ideal infinitely long shifted filter coefficients are the inverse DTFT of this amplitude given by