m0536 FIR Filter Example 2.7 2000/08/10 2009/06/09 17:41:52.152 GMT-5 Don Johnson Don Johnson dhj@rice.edu Don Johnson Don Johnson dhj@rice.edu Melissa Selik Melissa Selik mselik@alumni.rice.edu John Austin Cottrell John Cottrell jac3@rice.edu Don Johnson Don Johnson dhj@rice.edu analog analog signals buffering digital signal processing discrete-time discrete-time filtering DSP example examples fast Fourier transform FFT filtering finite impulse response FIR frequency domain Science and Technology An example of using a Finite Impulse Response filter. en We want to lowpass filter a signal that contains a sinusoid and a significant amount of noise. shows a portion of this signal's waveform. If it weren't for the overlaid sinusoid, discerning the sine wave in the signal is virtually impossible. One of the primary applications of linear filters is noise removal: preserve the signal by matching filter's passband with the signal's spectrum and greatly reduce all other frequency components that may be present in the noisy signal.

The noisy input signal is sectioned into length-48 frames, each of which is filtered using frequency-domain techniques. Each filtered section is added to other outputs that overlap to create the signal equivalent to having filtered the entire input. The sinusoidal component of the signal is shown as the red dashed line.

A smart Rice engineer has selected a FIR filter having a unit-sample response corresponding a period-17 sinusoid: h n 1 2 n 1 17 18 , n 0 … 16 , which makes q 16 . Its frequency response (determined by computing the discrete Fourier transform) is shown in . To apply, we can select the length of each section so that the frequency-domain filtering approach is maximally efficient: Choose the section length N x so that N x q is a power of two. To use a length-64 FFT, each section must be 48 samples long. Filtering with the difference equation would require 33 computations per output while the frequency domain requires a little over 16; this frequency-domain implementation is over twice as fast! shows how frequency-domain filtering works.

The figure shows the unit-sample response of a length-17 Hanning filter on the left and the frequency response on the right. This filter functions as a lowpass filter having a cutoff frequency of about 0.1.

We note that the noise has been dramatically reduced, with a sinusoid now clearly visible in the filtered output. Some residual noise remains because noise components within the filter's passband appear in the output as well as the signal. Note that when compared to the input signal's sinusoidal component, the output's sinusoidal component seems to be delayed. What is the source of this delay? Can it be removed? The delay is not computational delay here—the plot shows the first output value is aligned with the filter's first input—although in real systems this is an important consideration. Rather, the delay is due to the filter's phase shift: A phase-shifted sinusoid is equivalent to a time-delayed one: 2 f n φ 2 f n φ 2 f . All filters have phase shifts. This delay could be removed if the filter introduced no phase shift. Such filters do not exist in analog form, but digital ones can be programmed, but not in real time. Doing so would require the output to emerge before the input arrives!