Accelerating the pace of engineering and science

fir1

Window-based finite impulse response filter design

Syntax

b = fir1(n,Wn)
b = fir1(n,Wn,'ftype')
b = fir1(n,Wn,window)
b = fir1(n,Wn,'ftype',window)
b = fir1(...,'normalization')

Description

fir1 implements the classical method of windowed linear-phase FIR digital filter design [1]. It designs filters in standard lowpass, highpass, bandpass, and bandstop configurations. By default the filter is normalized so that the magnitude response of the filter at the center frequency of the passband is 0 dB.

 Note   Use fir2 for windowed filters with arbitrary frequency response.

b = fir1(n,Wn) returns row vector b containing the n+1 coefficients of an order n lowpass FIR filter. This is a Hamming-window based, linear-phase filter with normalized cutoff frequency Wn. The output filter coefficients, b, are ordered in descending powers of z.

$B\left(z\right)=b\left(1\right)+b\left(2\right){z}^{-1}+\dots +b\left(n+1\right){z}^{-N}$

Wn is a number between 0 and 1, where 1 corresponds to the Nyquist frequency.

If Wn is a two-element vector, Wn = [w1 w2], fir1 returns a bandpass filter with passband w1 < ω< w2.

If Wn is a multi-element vector, Wn = [w1 w2 w3 w4 w5 ... wn], fir1 returns an order n multiband filter with bands 0 < ω< w1, w1 < ω< w2, ..., wn < ω< 1.

By default, the filter is scaled so that the center of the first passband has a magnitude of exactly 1 after windowing.

b = fir1(n,Wn,'ftype') specifies a filter type, where 'ftype' is:

• 'high' for a highpass filter with cutoff frequency Wn.

• 'stop' for a bandstop filter, if Wn = [w1 w2]. The stopband frequency range is specified by this interval.

• 'DC-1' to make the first band of a multiband filter a passband.

• 'DC-0' to make the first band of a multiband filter a stopband.

fir1 always uses an even filter order for the highpass and bandstop configurations. This is because for odd orders, the frequency response at the Nyquist frequency is 0, which is inappropriate for highpass and bandstop filters. If you specify an odd-valued n, fir1 increments it by 1.

b = fir1(n,Wn,window) uses the window specified in column vector window for the design. The vector window must be n+1 elements long. If no window is specified, fir1 uses a Hamming window (see hamming) of length n+1.

b = fir1(n,Wn,'ftype',window) accepts both 'ftype' and window parameters.

b = fir1(...,'normalization') specifies whether or not the filter magnitude is normalized. The string 'normalization' can be the following:

• 'scale' (default): Normalize the filter so that the magnitude response of the filter at the center frequency of the passband is 0 dB.

• 'noscale': Do not normalize the filter.

The group delay of the FIR filter designed by fir1 is n/2.

Examples

expand all

FIR Bandpass Filter

Design a 48th-order FIR bandpass filter with passband .

b = fir1(48,[0.35 0.65]);
freqz(b,1,512)


FIR Highpass Filter

The file chirp.mat contains a signal, y, that has most of its power above fs/4, or half the Nyquist frequency. Load the file. Design a 34th-order FIR highpass filter to attenuate the components of the signal below fs/4. Use a cutoff frequency of 0.48 and a Chebyshev window with 30 dB of ripple.

load chirp
b = fir1(34,0.48,'high',chebwin(35,30));
freqz(b,1,512)


expand all

Algorithms

fir1 uses the window method of FIR filter design [1]. If w(n) denotes a window, where 1 ≤ nN, and the impulse response of the ideal filter is h(n), then the windowed digital filter coefficients are given by

$b\left(n\right)=w\left(n\right)h\left(n\right)\text{ }1\le n\le N$

References

[1] Digital Signal Processing Committee of the IEEE Acoustics, Speech, and Signal Processing Society, eds. Programs for Digital Signal Processing. New York: IEEE Press, 1979, algorithm 5.2.