Sinc filter

Sinc-in-time is an ideal filter that removes all frequency components above a given cutoff frequency, without attenuating lower frequencies, and has linear phase response. It may thus be considered a brick-wall filter or rectangular filter.

Its impulse response is a sinc function in the time domain:

$${\displaystyle {\frac {\sin(\pi t)}{\pi t}}}$$ 

while its frequency response is a 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|<B,\end{cases}}}$ 

where ${\displaystyle B}$  (representing its bandwidth) is an arbitrary cutoff frequency.

Its impulse response is given by the inverse Fourier transform of its frequency response:

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

where sinc is the normalized sinc function.

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-in-time 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}t\right)}$ 

${\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-in-time filters (since the filters are zero phase, their magnitude responses subtract directly):^[1]

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

${\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-in-time filter, which makes it clear that the Dirac delta function is the limit of a narrow-in-time sinc-in-time filter:

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

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

As the sinc-in-time filter has infinite impulse response in both positive and negative time directions, it is non-causal and has an infinite delay (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). However, it is used in conceptual demonstrations or proofs, such as the sampling theorem and the Whittaker–Shannon interpolation formula.

Sinc-in-time filters must be approximated for real-world (non-abstract) applications, typically by windowing and truncating an ideal sinc-in-time filter kernel, but doing so reduces its ideal properties. This applies to other brick-wall filters built using sinc-in-time filters.

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.

Transmission plots for group averaging filters using 1000 Hz sampling rate:

 

4-sample averaging (above)

 

32-sample averaging (above)

 

16-sample averaging (above), extended to 4x the Nyquist frequency. Because the transfer function is periodic, this repeated pattern continues forever.

The simplest implementation of a sinc-in-frequency filter uses a boxcar impulse response to produce a simple moving average (specifically if divide by the number of samples), also known as accumulate-and-dump filter (specifically if simply sum without a division). It can be modeled as a FIR filter with all ${\displaystyle N}$  coefficients equal. It is sometimes cascaded to produce higher-order moving averages (see Finite impulse response § Moving average example and cascaded integrator–comb filter).

This filter can be used for crude but fast and easy downsampling (a.k.a. decimation) by a factor of ${\displaystyle N.}$  The simplicity of the filter (accumulate ${\displaystyle N}$  data samples, output the accumulator result, zero the accumulator, and repeat) is foiled by its mediocre low-pass capabilities. Its poorest attenuation in the stop-band is -13.3 dB^[2] and most high frequency components are only slightly more attenuated than that. An ${\displaystyle N}$ -sample filter sampled at ${\displaystyle f_{S}}$  will alias all non-fully attenuated signal components lying above ${\textstyle {\frac {f_{S}}{2N}}}$  to the baseband ranging from DC to ${\textstyle {\frac {f_{S}}{2N}}.}$ 

A group averaging filter processing ${\displaystyle N}$  samples has ${\displaystyle {\tfrac {N}{2}}}$  transmission zeroes evenly-spaced by ${\displaystyle {\tfrac {f_{S}}{N}},}$  with the lowest zero at ${\displaystyle {\tfrac {f_{S}}{N}}}$  and the highest zero at ${\displaystyle {\tfrac {f_{S}}{2}}}$  (the Nyquist frequency). Above the Nyquist frequency, the frequency response is mirrored and then is repeated periodically above ${\displaystyle f_{S}}$  forever.

The magnitude of the frequency response (plotted in these graphs) is useful when one wants to know how much frequencies are attenuated. Though the sinc function really oscillates between negative and positive values, negative values of the frequency response simply correspond to a 180-degree phase shift.

An inverse sinc filter may be used for equalization in the digital domain (e.g. a FIR filter) or analog domain (e.g. opamp filter) to counteract undesired attenuation in the frequency band of interest to provide a flat frequency response.^[3]

See Window function § Rectangular window for application of the sinc kernel as the simplest windowing function.

• Lanczos resampling
• Aliasing
• Anti-aliasing filter

• Brick Wall Digital Filters and Phase Deviations
• Brick-wall filters