Sinc filter – Wikipedia

Ideal low-pass filter

In signal processing, a sinc filter is an idealized filter that removes all frequency components above a given cutoff frequency, without affecting lower frequencies, and has linear phase response. The filter’s impulse response is a sinc function in the time domain

and its frequency response is a rectangular function. It is an “ideal” low-pass filter in the frequency sense, perfectly passing low frequencies, perfectly cutting high frequencies; and thus may be considered to be a brick-wall filter.

Real-time filters can only approximate this ideal, since an ideal sinc filter (a.k.a. rectangular filter) is non-causal and has an infinite delay, but it is commonly found in conceptual demonstrations or proofs, such as the sampling theorem and the Whittaker–Shannon interpolation formula.

In mathematical terms, the desired frequency response is the rectangular function:

${displaystyle H(f)=operatorname {rect} left({frac {f}{2B}}right)={begin{cases}0,&{text{if }}|f|>B,\{frac {1}{2}},&{text{if }}|f|=B,\1,&{text{if }}|f|$

where B is an arbitrary cutoff frequency (a.k.a. bandwidth). The impulse response of such a filter is given by the inverse Fourier transform of the frequency response:

${displaystyle {begin{aligned}h(t)={mathcal {F}}^{-1}{H(f)}&=int _{-B}^{B}exp(2pi ift),df\&=2Boperatorname {sinc} (2Bt)end{aligned}}}$

where sinc is the normalized sinc function.

As the sinc filter has infinite impulse response in both positive and negative time directions, it must be approximated for real-world (non-abstract) applications; a windowed sinc filter is often used instead. Windowing and truncating a sinc filter kernel in order to use it on any practical real world data set reduces its ideal properties.

Brick-wall filters[edit]

An idealized electronic filter with full transmission in the pass band, complete attenuation in the stop band, and abrupt transitions is known colloquially as a “brick-wall filter” (in reference to the shape of the transfer function). The sinc filter is a brick-wall low-pass filter, from which brick-wall band-pass filters and high-pass filters are easily constructed.

The lowpass filter with brick-wall cutoff at frequency B_L has impulse response and transfer function given by:

${displaystyle h_{LPF}(t)=2B_{L}operatorname {sinc} left(2B_{L}tright)}$

${displaystyle H_{LPF}(f)=operatorname {rect} left({frac {f}{2B_{L}}}right).}$

The band-pass filter with lower band edge B_L and upper band edge B_H is just the difference of two such sinc filters (since the filters are zero phase, their magnitude responses subtract directly):^[1]

${displaystyle h_{BPF}(t)=2B_{H}operatorname {sinc} left(2B_{H}tright)-2B_{L}operatorname {sinc} left(2B_{L}tright)}$

${displaystyle H_{BPF}(f)=operatorname {rect} left({frac {f}{2B_{H}}}right)-operatorname {rect} left({frac {f}{2B_{L}}}right).}$

The high-pass filter with lower band edge B_H is just a transparent filter minus a sinc filter, which makes it clear that the Dirac delta function is the limit of a narrow-in-time sinc filter:

${displaystyle h_{HPF}(t)=delta (t)-2B_{H}operatorname {sinc} left(2B_{H}tright)}$

${displaystyle H_{HPF}(f)=1-operatorname {rect} left({frac {f}{2B_{H}}}right).}$

Brick-wall filters that run in realtime are not physically realizable as they have infinite latency (i.e., its compact support in the frequency domain forces its time response not to have compact support meaning that it is ever-lasting) and infinite order (i.e., the response cannot be expressed as a linear differential equation with a finite sum), but approximate implementations are sometimes used and they are frequently called brick-wall filters.^{[citation needed]}

Frequency-domain sinc[edit]

The name “sinc filter” is applied also to the filter shape that is rectangular in time and a sinc function in frequency, as opposed to the ideal low-pass sinc filter, which is sinc in time and rectangular in frequency. In case of confusion, one may refer to these as sinc-in-frequency and sinc-in-time, according to which domain the filter is sinc in.

Sinc-in-frequency Cascaded integrator–comb (CIC) filters, among many other applications, are almost universally used for decimating delta-sigma ADCs, as they are easy to implement and nearly optimal for this use.^[2]

The simplest implementation of a Sinc-in-frequency filter is a group-averaging filter, also known as accumulate-and-dump filter. This filter also performs a data rate reduction.

Filter transmission function of a 4sample group averaging filter

It collects N data samples, accumulates them and provides the accumulator value as output. Thus, the decimation factor of this filter is N.
It can be modeled as a FIR filter with all N coefficients equal, followed by a N-time downsampling block. The simplicity of the filter, requiring just an accumulator as central data processing block, is foiled with strong aliasing effects: an N sample filter aliases all non-fully attenuated signal components lying above

to the baseband ranging from 0 to

(f_S is the input sample rate).

Filter transmission function of a 32sample group averaging filter

A group averaging filter processing N samples has N/2 transmission zeroes.
The picture “transmission function of a 16sample group averaging filter” shows how the transmission function looks above the Nyquist frequency.

transmission function of a 16sample group averaging filter working on a 1kHz input data rate, extended to 4x the Nyquist frequency

Stability[edit]

The sinc filter is not bounded-input–bounded-output (BIBO) stable. That is, a bounded input can produce an unbounded output, because the integral of the absolute value of the sinc function is infinite. A bounded input that produces an unbounded output is sgn(sinc(t)). Another is sin(2πBt)u(t), a sine wave starting at time 0, at the cutoff frequency.

See also[edit]

References[edit]

External links[edit]