Low-pass filter

A low-pass filter is a filter that passes signals with a frequency lower than a selected cutoff frequency and attenuates signals with frequencies higher than the cutoff frequency. The exact frequency response of the filter depends on the filter design. The filter is sometimes called a high-cut filter, or treble-cut filter in audio applications. A low-pass filter is the complement of a high-pass filter.

In optics, high-pass and low-pass may have different meanings, depending on whether referring to the frequency or wavelength of light, since these variables are inversely related. High-pass frequency filters would act as low-pass wavelength filters, and vice versa. For this reason, it is a good practice to refer to wavelength filters as short-pass and long-pass to avoid confusion, which would correspond to high-pass and low-pass frequencies. ^[1]

Low-pass filters exist in many different forms, including electronic circuits such as a hiss filter used in audio, anti-aliasing filters for conditioning signals before analog-to-digital conversion, digital filters for smoothing sets of data, acoustic barriers, blurring of images, and so on. The moving average operation used in fields such as finance is a particular kind of low-pass filter and can be analyzed with the same signal processing techniques as are used for other low-pass filters. Low-pass filters provide a smoother form of a signal, removing the short-term fluctuations and leaving the longer-term trend.

Filter designers will often use the low-pass form as a prototype filter. That is a filter with unity bandwidth and impedance. The desired filter is obtained from the prototype by scaling for the desired bandwidth and impedance and transforming into the desired bandform (that is, low-pass, high-pass, band-pass or band-stop).

Examples

Examples of low-pass filters occur in acoustics, optics and electronics.

A stiff physical barrier tends to reflect higher sound frequencies, acting as an acoustic low-pass filter for transmitting sound. When music is playing in another room, the low notes are easily heard, while the high notes are attenuated.

An optical filter with the same function can correctly be called a low-pass filter, but conventionally is called a longpass filter (low frequency is long wavelength), to avoid confusion. ^[1]

In an electronic low-pass RC filter for voltage signals, high frequencies in the input signal are attenuated, but the filter has little attenuation below the cutoff frequency determined by its RC time constant. For current signals, a similar circuit, using a resistor and capacitor in parallel, works in a similar manner. (See current divider discussed in more detail below.)

Electronic low-pass filters are used on inputs to subwoofers and other types of loudspeakers, to block high pitches that they cannot efficiently reproduce. Radio transmitters use low-pass filters to block harmonic emissions that might interfere with other communications. The tone knob on many electric guitars is a low-pass filter used to reduce the amount of treble in the sound. An integrator is another time constant low-pass filter. ^[2]

Telephone lines fitted with DSL splitters use low-pass filters to separate DSL from POTS signals (and high-pass vice versa), which share the same pair of wires (transmission channel). ^{[3] [4]}

Low-pass filters also play a significant role in the sculpting of sound created by analogue and virtual analogue synthesisers. See subtractive synthesis.

A low-pass filter is used as an anti-aliasing filter before sampling and for reconstruction in digital-to-analog conversion.

Ideal and real filters

The sinc function, the time-domain impulse response of an ideal low-pass filter. The ripples of a true sinc extend infinitely to the left and right while getting smaller and smaller, but this particular graph is truncated.

The gain-magnitude frequency response of a first-order (one-pole) low-pass filter. Power gain is shown in decibels (i.e., a 3 dB decline reflects an additional half-power attenuation). Angular frequency is shown on a logarithmic scale in units of radians per second.

An ideal low-pass filter completely eliminates all frequencies above the cutoff frequency while passing those below unchanged; its frequency response is a rectangular function and is a brick-wall filter. The transition region present in practical filters does not exist in an ideal filter. An ideal low-pass filter can be realized mathematically (theoretically) by multiplying a signal by the rectangular function in the frequency domain or, equivalently, convolution with its impulse response, a sinc function, in the time domain.

However, the ideal filter is impossible to realize without also having signals of infinite extent in time, and so generally needs to be approximated for real ongoing signals, because the sinc function's support region extends to all past and future times. The filter would therefore need to have infinite delay, or knowledge of the infinite future and past, to perform the convolution. It is effectively realizable for pre-recorded digital signals by assuming extensions of zero into the past and future, or, more typically, by making the signal repetitive and using Fourier analysis.

Real filters for real-time applications approximate the ideal filter by truncating and windowing the infinite impulse response to make a finite impulse response; applying that filter requires delaying the signal for a moderate period of time, allowing the computation to "see" a little bit into the future. This delay is manifested as phase shift. Greater accuracy in approximation requires a longer delay.

Truncating an ideal low-pass filter result in ringing artifacts via the Gibbs phenomenon, which can be reduced or worsened by the choice of windowing function. Design and choice of real filters involves understanding and minimizing these artifacts. For example, simple truncation of the sinc function will create severe ringing artifacts, which can be reduced using window functions that drop off more smoothly at the edges. ^[5]

The Whittaker–Shannon interpolation formula describes how to use a perfect low-pass filter to reconstruct a continuous signal from a sampled digital signal. Real digital-to-analog converters uses real filter approximations.

Time response

The time response of a low-pass filter is found by solving the response to the simple low-pass RC filter.

A simple low-pass RC filter

Using Kirchhoff's Laws we arrive at the differential equation ^[6]

${\displaystyle v_{\text{out}}(t)=v_{\text{in}}(t)-RC{\frac {\operatorname {d} v_{\text{out}}}{\operatorname {d} t}}}$

Step input response example

If we let ${\displaystyle v_{\text{in}}(t)}$ be a step function of magnitude ${\displaystyle V_{i}}$ then the differential equation has the solution ^[7]

${\displaystyle v_{\text{out}}(t)=V_{i}(1-e^{-\omega _{0}t}),}$

where ${\displaystyle \omega _{0}={1 \over RC}}$ is the cutoff frequency of the filter.

Frequency response

The most common way to characterize the frequency response of a circuit is to find its Laplace transform ^[6] transfer function, ${\displaystyle H(s)={V_{\rm {out}}(s) \over V_{\rm {in}}(s)}}$. Taking the Laplace transform of our differential equation and solving for ${\displaystyle H(s)}$ we get

${\displaystyle H(s)={V_{\rm {out}}(s) \over V_{\rm {in}}(s)}={\omega _{0} \over s+\omega _{0}}}$

Difference equation through discrete time sampling

A discrete difference equation is easily obtained by sampling the step input response above at regular intervals of ${\displaystyle nT}$ where ${\displaystyle n=0,1,...}$ and ${\displaystyle T}$ is the time between samples. Taking the difference between two consecutive samples we have

${\displaystyle v_{\rm {out}}(nT)-v_{\rm {out}}((n-1)T)=V_{i}(1-e^{-\omega _{0}nT})-V_{i}(1-e^{-\omega _{0}((n-1)T)})}$

Solving for ${\displaystyle v_{\rm {out}}(nT)}$ we get

${\displaystyle v_{\rm {out}}(nT)=\beta v_{\rm {out}}((n-1)T)+(1-\beta )V_{i}}$

Where ${\displaystyle \beta =e^{-\omega _{0}T}}$

Using the notation ${\displaystyle V_{n}=v_{\rm {out}}(nT)}$ and ${\displaystyle v_{n}=v_{\rm {in}}(nT)}$, and substituting our sampled value, ${\displaystyle v_{n}=V_{i}}$, we get the difference equation

${\displaystyle V_{n}=\beta V_{n-1}+(1-\beta )v_{n}}$

Error analysis

Comparing the reconstructed output signal from the difference equation, ${\displaystyle V_{n}=\beta V_{n-1}+(1-\beta )v_{n}}$, to the step input response, ${\displaystyle v_{\text{out}}(t)=V_{i}(1-e^{-\omega _{0}t})}$, we find that there is an exact reconstruction (0% error). This is the reconstructed output for a time-invariant input. However, if the input is time variant, such as ${\displaystyle v_{\text{in}}(t)=V_{i}\sin(\omega t)}$, this model approximates the input signal as a series of step functions with duration ${\displaystyle T}$ producing an error in the reconstructed output signal. The error produced from time variant inputs is difficult to quantify^{[ citation needed ]} but decreases as ${\displaystyle T\rightarrow 0}$.

Discrete-time realization

For another method of conversion from continuous- to discrete-time, see Bilinear transform.

Many digital filters are designed to give low-pass characteristics. Both infinite impulse response and finite impulse response low pass filters, as well as filters using Fourier transforms, are widely used.

Simple infinite impulse response filter

The effect of an infinite impulse response low-pass filter can be simulated on a computer by analyzing an RC filter's behavior in the time domain, and then discretizing the model.

A simple low-pass RC filter

From the circuit diagram to the right, according to Kirchhoff's Laws and the definition of capacitance:

${\displaystyle v_{\text{in}}(t)-v_{\text{out}}(t)=R\;i(t)}$

(V)

${\displaystyle Q_{c}(t)=C\,v_{\text{out}}(t)}$

(Q)

${\displaystyle i(t)={\frac {\operatorname {d} Q_{c}}{\operatorname {d} t}}}$

(I)

where ${\displaystyle Q_{c}(t)}$ is the charge stored in the capacitor at time t. Substituting equation Q into equation I gives ${\displaystyle i(t)\;=\;C{\frac {\operatorname {d} v_{\text{out}}}{\operatorname {d} t}}}$, which can be substituted into equation V so that

${\displaystyle v_{\text{in}}(t)-v_{\text{out}}(t)=RC{\frac {\operatorname {d} v_{\text{out}}}{\operatorname {d} t}}.}$

This equation can be discretized. For simplicity, assume that samples of the input and output are taken at evenly spaced points in time separated by ${\displaystyle \Delta _{T}}$ time. Let the samples of ${\displaystyle v_{\text{in}}}$ be represented by the sequence ${\displaystyle (x_{1},\,x_{2},\,\ldots ,\,x_{n})}$, and let ${\displaystyle v_{\text{out}}}$ be represented by the sequence ${\displaystyle (y_{1},\,y_{2},\,\ldots ,\,y_{n})}$, which correspond to the same points in time. Making these substitutions,

${\displaystyle x_{i}-y_{i}=RC\,{\frac {y_{i}-y_{i-1}}{\Delta _{T}}}.}$

Rearranging terms gives the recurrence relation

${\displaystyle y_{i}=\overbrace {x_{i}\left({\frac {\Delta _{T}}{RC+\Delta _{T}}}\right)} ^{\text{Input contribution}}+\overbrace {y_{i-1}\left({\frac {RC}{RC+\Delta _{T}}}\right)} ^{\text{Inertia from previous output}}.}$

That is, this discrete-time implementation of a simple RC low-pass filter is the exponentially weighted moving average

${\displaystyle y_{i}=\alpha x_{i}+(1-\alpha )y_{i-1}\qquad {\text{where}}\qquad \alpha$ :={\frac {\Delta _{T}}{RC+\Delta _{T}}}.}

By definition, the smoothing factor is within the range ${\displaystyle 0\;\leq \;\alpha \;\leq \;1}$. The expression for α yields the equivalent time constant RC in terms of the sampling period ${\displaystyle \Delta _{T}}$ and smoothing factor α,

${\displaystyle RC=\Delta _{T}\left({\frac {1-\alpha }{\alpha }}\right).}$

Recalling that

${\displaystyle f_{c}={\frac {1}{2\pi RC}}}$ so ${\displaystyle RC={\frac {1}{2\pi f_{c}}},}$

note α and ${\displaystyle f_{c}}$ are related by,

${\displaystyle \alpha ={\frac {2\pi \Delta _{T}f_{c}}{2\pi \Delta _{T}f_{c}+1}}}$

and

${\displaystyle f_{c}={\frac {\alpha }{(1-\alpha )2\pi \Delta _{T}}}.}$

If α=0.5, then the RC time constant equals the sampling period. If ${\displaystyle \alpha \;\ll \;0.5}$, then RC is significantly larger than the sampling interval, and ${\displaystyle \Delta _{T}\;\approx \;\alpha RC}$.

The filter recurrence relation provides a way to determine the output samples in terms of the input samples and the preceding output. The following pseudocode algorithm simulates the effect of a low-pass filter on a series of digital samples:

// Return RC low-pass filter output samples, given input samples, // time interval dt, and time constant RCfunction lowpass(real[1..n] x, real dt, real RC)     varreal[1..n] y     varreal α := dt / (RC + dt)     y[1] := α * x[1]     for i from 2 to n         y[i] := α * x[i] + (1-α) * y[i-1]     return y

The loop that calculates each of the n outputs can be refactored into the equivalent:

for i from 2 to n         y[i] := y[i-1] + α * (x[i] - y[i-1])

That is, the change from one filter output to the next is proportional to the difference between the previous output and the next input. This exponential smoothing property matches the exponential decay seen in the continuous-time system. As expected, as the time constant RC increases, the discrete-time smoothing parameter ${\displaystyle \alpha }$ decreases, and the output samples ${\displaystyle (y_{1},\,y_{2},\,\ldots ,\,y_{n})}$ respond more slowly to a change in the input samples ${\displaystyle (x_{1},\,x_{2},\,\ldots ,\,x_{n})}$; the system has more inertia . This filter is an infinite-impulse-response (IIR) single-pole low-pass filter.

Finite impulse response

Finite-impulse-response filters can be built that approximate the sinc function time-domain response of an ideal sharp-cutoff low-pass filter. For minimum distortion, the finite impulse response filter has an unbounded number of coefficients operating on an unbounded signal. In practice, the time-domain response must be time truncated and is often of a simplified shape; in the simplest case, a running average can be used, giving a square time response. ^[8]

Fourier transform

For non-realtime filtering, to achieve a low pass filter, the entire signal is usually taken as a looped signal, the Fourier transform is taken, filtered in the frequency domain, followed by an inverse Fourier transform. Only O(n log(n)) operations are required compared to O(n²) for the time domain filtering algorithm.

This can also sometimes be done in real time, where the signal is delayed long enough to perform the Fourier transformation on shorter, overlapping blocks.

Continuous-time realization

Plot of the gain of Butterworth low-pass filters of orders 1 through 5, with cutoff frequency ${\displaystyle \omega _{0}=1}$. Note that the slope is 20n dB/decade where n is the filter order.

There are many different types of filter circuits, with different responses to changing frequency. The frequency response of a filter is generally represented using a Bode plot, and the filter is characterized by its cutoff frequency and rate of frequency rolloff. In all cases, at the cutoff frequency, the filter attenuates the input power by half or 3 dB. So the order of the filter determines the amount of additional attenuation for frequencies higher than the cutoff frequency.

• A first-order filter, for example, reduces the signal amplitude by half (so power reduces by a factor of 4, or 6 dB), every time the frequency doubles (goes up one octave); more precisely, the power rolloff approaches 20 dB per decade in the limit of high frequency. The magnitude Bode plot for a first-order filter looks like a horizontal line below the cutoff frequency, and a diagonal line above the cutoff frequency. There is also a "knee curve" at the boundary between the two, smoothly transitioning between the two straight-line regions. If the transfer function of a first-order low-pass filter has a zero as well as a pole, the Bode plot flattens out again, at some maximum attenuation of high frequencies; such an effect is caused for example by a little bit of the input leaking around the one-pole filter; this one-pole–one-zero filter is still a first-order low-pass. See Pole–zero plot and RC circuit.
• A second-order filter attenuates high frequencies more steeply. The Bode plot for this type of filter resembles that of a first-order filter, except that it falls off more quickly. For example, a second-order Butterworth filter reduces the signal amplitude to one-fourth of its original level every time the frequency doubles (so power decreases by 12 dB per octave, or 40 dB per decade). Other all-pole second-order filters may roll off at different rates initially depending on their Q factor, but approach the same final rate of 12 dB per octave; as with the first-order filters, zeroes in the transfer function can change the high-frequency asymptote. See RLC circuit.
• Third- and higher-order filters are defined similarly. In general, the final rate of power rolloff for an order- n all-pole filter is 6n dB per octave (20n dB per decade).

On any Butterworth filter, if one extends the horizontal line to the right and the diagonal line to the upper-left (the asymptotes of the function), they intersect at exactly the cutoff frequency, 3 dB below the horizontal line. The various types of filters (Butterworth filter, Chebyshev filter, Bessel filter, etc.) all have different-looking knee curves. Many second-order filters have "peaking" or resonance that puts their frequency response above the horizontal line at this peak.

The meanings of 'low' and 'high'—that is, the cutoff frequency—depend on the characteristics of the filter. The term "low-pass filter" merely refers to the shape of the filter's response; a high-pass filter could be built that cuts off at a lower frequency than any low-pass filter—it is their responses that set them apart. Electronic circuits can be devised for any desired frequency range, right up through microwave frequencies (above 1 GHz) and higher.

Laplace notation

Continuous-time filters can also be described in terms of the Laplace transform of their impulse response, in a way that lets all characteristics of the filter be easily analyzed by considering the pattern of poles and zeros of the Laplace transform in the complex plane. (In discrete time, one can similarly consider the Z-transform of the impulse response.)

For example, a first-order low-pass filter can be described in Laplace notation as:

${\displaystyle {\frac {\text{Output}}{\text{Input}}}=K{\frac {1}{\tau s+1}}}$

where s is the Laplace transform variable, τ is the filter time constant, and K is the gain of the filter in the passband.

Electronic low-pass filters

First order

RC filter

Main article: RC circuit § Series circuit

Passive, first order low-pass RC filter

One simple low-pass filter circuit consists of a resistor in series with a load, and a capacitor in parallel with the load. The capacitor exhibits reactance, and blocks low-frequency signals, forcing them through the load instead. At higher frequencies, the reactance drops, and the capacitor effectively functions as a short circuit. The combination of resistance and capacitance gives the time constant of the filter ${\displaystyle \tau \;=\;RC}$ (represented by the Greek letter tau). The break frequency, also called the turnover frequency, corner frequency, or cutoff frequency (in hertz), is determined by the time constant:

${\displaystyle f_{\mathrm {c} }={1 \over 2\pi \tau }={1 \over 2\pi RC}}$

or equivalently (in radians per second):

${\displaystyle \omega _{\mathrm {c} }={1 \over \tau }={1 \over RC}}$

This circuit may be understood by considering the time the capacitor needs to charge or discharge through the resistor:

• At low frequencies, there is plenty of time for the capacitor to charge up to practically the same voltage as the input voltage.
• At high frequencies, the capacitor only has time to charge up a small amount before the input switches direction. The output goes up and down only a small fraction of the amount the input goes up and down. At double the frequency, there's only time for it to charge up half the amount.

Another way to understand this circuit is through the concept of reactance at a particular frequency:

• Since direct current (DC) cannot flow through the capacitor, DC input must flow out the path marked ${\displaystyle V_{\mathrm {out} }}$ (analogous to removing the capacitor).
• Since alternating current (AC) flows very well through the capacitor, almost as well as it flows through a solid wire, AC input flows out through the capacitor, effectively short circuiting to the ground (analogous to replacing the capacitor with just a wire).

The capacitor is not an "on/off" object (like the block or pass fluidic explanation above). The capacitor variably acts between these two extremes. It is the Bode plot and frequency response that show this variability.

RL filter

Main article: RL circuit § Series circuit

A resistor–inductor circuit or RL filter is an electric circuit composed of resistors and inductors driven by a voltage or current source. A first-order RL circuit is composed of one resistor and one inductor and is the simplest type of RL circuit.

A first-order RL circuit is one of the simplest analogue infinite impulse response electronic filters. It consists of a resistor and an inductor, either in series driven by a voltage source or in parallel driven by a current source.

Second order

RLC filter

RLC circuit as a low-pass filter

An RLC circuit (the letters R, L, and C can be in a different sequence) is an electrical circuit consisting of a resistor, an inductor, and a capacitor, connected in series or in parallel. The RLC part of the name is due to those letters being the usual electrical symbols for resistance, inductance, and capacitance, respectively. The circuit forms a harmonic oscillator for current and will resonate in a similar way as an LC circuit will. The main difference that the presence of the resistor makes is that any oscillation induced in the circuit will die away over time if it is not kept going by a source. This effect of the resistor is called damping. The presence of the resistance also reduces the peak resonant frequency somewhat. Some resistance is unavoidable in real circuits, even if a resistor is not specifically included as a component. An ideal, pure LC circuit is an abstraction for the purpose of theory.

There are many applications for this circuit. They are used in many different types of oscillator circuits. Another important application is for tuning, such as in radio receivers or television sets, where they are used to select a narrow range of frequencies from the ambient radio waves. In this role, the circuit is often called a tuned circuit. An RLC circuit can be used as a band-pass filter, band-stop filter, low-pass filter, or high-pass filter. The RLC filter is described as a second-order circuit, meaning that any voltage or current in the circuit can be described by a second-order differential equation in circuit analysis.

Second-Order Low-Pass Filter Standard Form

The transfer function ${\displaystyle H_{LP}(f)}$ of a second-order low-pass filter can be expressed as a function of frequency ${\displaystyle f}$ as shown in Equation 1, the Second-Order Low-Pass Filter Standard Form.

${\displaystyle H_{LP}(f)=-{\frac {K}{f_{FSF}\cdot f_{c}^{2}+{\frac {1}{Q}}\cdot jf_{FSF}\cdot f_{c}+1}}\quad (1)}$

In this equation, ${\displaystyle f}$ is the frequency variable, ${\displaystyle f_{c}}$ is the cutoff frequency, ${\displaystyle f_{FSF}}$ is the frequency scaling factor, and ${\displaystyle Q}$ is the quality factor. Equation 1 describes three regions of operation: below cutoff, in the area of cutoff, and above cutoff. For each area, Equation 1 reduces to:

• ${\displaystyle f\ll f_{c}}$: ${\displaystyle H_{LP}(f)\approx K}$ - The circuit passes signals multiplied by the gain factor ${\displaystyle K}$.
• ${\displaystyle {\frac {f}{f_{c}}}=f_{FSF}}$: ${\displaystyle H_{LP}(f)=jKQ}$ - Signals are phase-shifted 90° and modified by the quality factor ${\displaystyle Q}$.
• ${\displaystyle f\gg f_{c}}$: ${\displaystyle H_{LP}(f)\approx -{\frac {K}{f_{FSF}\cdot f^{2}}}}$ - Signals are phase-shifted 180° and attenuated by the square of the frequency ratio. This behavior is detailed by Jim Karki in "Active Low-Pass Filter Design" (Texas Instruments, 2023). ^[9]

With attenuation at frequencies above ${\displaystyle f_{c}}$ increasing by a power of two, the last formula describes a second-order low-pass filter. The frequency scaling factor ${\displaystyle f_{FSF}}$ is used to scale the cutoff frequency of the filter so that it follows the definitions given before.

Higher order passive filters

Higher-order passive filters can also be constructed (see diagram for a third-order example).

A third-order low-pass filter (Cauer topology). The filter becomes a Butterworth filter with cutoff frequency ω_c=1 when (for example) C₂=4/3 farad, R₄=1 ohm, L₁=3/2 henry and L₃=1/2 henry.

Active electronic realization

An active low-pass filter

See also: operational amplifier applications § Inverting integrator, and Op amp integrator

An active low-pass filter adds an active device to create an active filter that allows for gain in the passband.

In the operational amplifier circuit shown in the figure, the cutoff frequency (in hertz) is defined as:

${\displaystyle f_{\text{c}}={\frac {1}{2\pi R_{2}C}}}$

or equivalently (in radians per second):

${\displaystyle \omega _{\text{c}}={\frac {1}{R_{2}C}}}$

The gain in the passband is −R₂/R₁, and the stopband drops off at −6 dB per octave (that is −20 dB per decade) as it is a first-order filter.

See also

• Electronics portal

• Baseband

References

1. Long Pass Filters and Short Pass Filters Information , retrieved 2017-10-04
2. Sedra, Adel; Smith, Kenneth C. (1991). Microelectronic Circuits, 3 ed . Saunders College Publishing. p.  60. ISBN   0-03-051648-X.
3. "ADSL filters explained". Epanorama.net. Retrieved 2013-09-24.
4. "Home Networking – Local Area Network". Pcweenie.com. 2009-04-12. Archived from the original on 2013-09-27. Retrieved 2013-09-24.
5. Mastering Windows: Improving Reconstruction
6. Hayt, William H. Jr. and Kemmerly, Jack E. (1978). Engineering Circuit Analysis. New York: McGRAW-HILL BOOK COMPANY. pp. 211–224, 684–729.
7. Boyce, William and DiPrima, Richard (1965). Elementary Differential Equations and Boundary Value Problems. New York: JOHN WILEY & SONS. pp. 11–24.
8. Whilmshurst, T H (1990) Signal recovery from noise in electronic instrumentation. ISBN   9780750300582
9. Active Low-Pass Filter Design" (Texas Instruments, 2023)

External links

• Low Pass Filter java simulator
• ECE 209: Review of Circuits as LTI Systems, a short primer on the mathematical analysis of (electrical) LTI systems.
• ECE 209: Sources of Phase Shift, an intuitive explanation of the source of phase shift in a low-pass filter. Also verifies simple passive LPF transfer function by means of trigonometric identity.
• C code generator for digital implementation of Butterworth, Bessel, and Chebyshev filters created by the late Dr. Tony Fisher of the University of York (York, England).