Advanced Low-Pass Filter

Component Icon

Description

The Low-Pass Filter component implements a low-pass filter on the input signal. The cutoff frequency f_c can be defined through the Component Properties menu or connected to the Cutoff frequency input of the component.

All the filters are designed to have a maximum magnitude of 1 (0 dB) in the passband.

Passband is the range of frequencies over which the magnitude response does not fall below ( ${\mathrm{-A}}_p$ ). For a low-pass filter it is the range between 0 and ${\omega }_p$ . ${\omega }_p$ is called the passband frequency.

Stopband is the range of frequencies over which the magnitude response does not rise above ( ${\mathrm{-A}}_s$ ). For a low-pass filter it is the range between 0 and ${\omega }_s$ . ${\omega }_s$ is called the stopband frequency.

Transition band is the range of frequencies between the passband and the stopband. Roll-off is the term used for the decrease of the magnitude after the passband i.e. it represents the transition band. When comparing roll-offs between filters, it is assumed that the comparison is done for filters with same number of poles.

Cutoff frequency is by definition the frequency in which the magnitude plot has has an attenuation of 3 dB. However, this component defines the cutoff frequency for the Chebyshev1 and Elliptic filters as the frequency in which the passband regime is over and after which the magnitude infinitely drops (the official term for this frequency is the passband frequency). Also, the Basic filter represents a filter whose transfer function has all poles within a defined cut-off frequency. Therefore, cutoff frequency is in the -3 dB spot only for the first order of the Basic filter.

Passband attenuation is the width of the passband ripple for the Chebyshev1 and Elliptic filters ( $A_p$ ). Stopband attenuation is the minimum attenuation in the stopband for the Chebyshev2 and Elliptic filters ( $A_s$ ). Basic and Butterworth filters do not have ripples, so they do not have these properties.

In order to implement discrete filters, a Bilinear discretization is applied to all filters.

Figure 3 shows the 7th Order Basic filter with 100 Hz cutoff frequency obtained from the component. It also shows the cutoff frequency obtained from the component, as well as the real 3 dB cutoff frequency. The link between those two frequencies is ${\omega }_{\mathrm{3\; dB}}={\omega }_c\sqrt{{10}^{\frac{3}{10N}}-1}$ ,

where ${\omega }_c=2\pi f_c$ is the radian cutoff frequency obtained from the component, and $N$ is the filter order.

The transfer function of the basic filter is:

$\frac{Y\left(s\right)}{U\left(s\right)}=\frac{{\omega }_c^N}{{(s+{\omega }_c)}^N}$

Figure 4 shows the 4th Order Butterworth filter with 100 Hz cutoff frequency obtained from the component. It also shows cutoff frequency obtained from the component, which is also the 3 dB cutoff frequency. Butterworth filter has a maximally flat passband, but a slower roll-off compared to Chebyshev1, Chebyshev2 and Elliptic filters.

The transfer function of the Butterworth filter is:

$\frac{Y\left(s\right)}{U\left(s\right)}={\displaystyle \prod _{k=1}^N}\frac{{\omega }_c}{(s-s_{pk})}$ ,

where $s_k$ are the poles obtained from:

$s_{pk}={\omega }_c{e}^{\frac{\mathrm{j(2k+N-1)\pi }}{\mathrm{2N}}}$ ,

for $k=1,2,\mathrm{...}N$ , $N$ is the filter order.

Figure 5 shows the 3rd Order Chebyshev1 filter with a 5 dB passband ripple ( $A_p$ ) and with 100 Hz cutoff frequency obtained from the component. As mentioned earlier, that cutoff frequency obtained from the component is actually the passband frequency of the filter. Chebyshev1 filter has a ripple in the passband, but also a steeper roll-off than Basic, Butterworth and Chebyshev2 filters. The link between the cutoff frequency obtained from the component and the 3 dB cutoff frequency is:

${\omega }_{\mathrm{3\; dB}}={\omega }_c\cos \left(\frac{\mathrm{arccosh}\left(\frac{1}{{\epsilon }_p}\right)}{N}\right)$ , for $A_p<3dB$

where ${\omega }_c=2\pi f_c$ is the radian cutoff frequency obtained from the component, and $N$ is the filter order. ${\epsilon }_p$ is the passband ripple factor obtained from:

${\epsilon }_p=\sqrt{{10}^{\frac{Ap}{10}}-1}$ .

The transfer function of the Chebyshev1 filter is:

$\frac{Y\left(s\right)}{U\left(s\right)}=\frac{1}{2^{\mathrm{N-1}}\varepsilon }{\displaystyle \prod _{k=1}^N}\frac{{\omega }_c}{(s-s_{pk})}$ ,

where $s_{pk}$ are the poles obtained from:

$s_{pk}={\omega }_c\left(-sinh\left(\frac{arsinh\left(\frac{1}{{\epsilon }_p}\right)}{N}\right)sin\left({\theta }_k\right)+jcosh\left(\frac{arsinh\left(\frac{1}{{\epsilon }_p}\right)}{N}\right)cos\left({\theta }_k\right)\right)$ , ${\theta }_k=\frac{\pi }{2}\frac{2k-1}{N}$ ,

for $k=1,2,\mathrm{...}N$

Figure 6 shows the 4th Order Chebyshev2 Filter with a 40 dB attenuation in the stopband ( $A_s$ ) and with 100 Hz cutoff frequency obtained from the component, which is also the 3 dB cutoff frequency. Chebyshev2 filter has a ripple in the stopband, but also a steeper roll-off than Basic and Butterworth filters. The link between cutoff and stopband frequency is

${\omega }_c={\omega }_s/\cos \left(\frac{\mathrm{arccosh}\left({\left({\epsilon }_s\sqrt{{10}^{0.3}-1}\right)}^{-1}\right)}{N}\right)$ ,

where ${\omega }_c=2\pi f_c$ is the radian cutoff frequency obtained from the component, ${\omega }_s$ is the stopband frequency, and $N$ is the filter order. ${\epsilon }_s$ is the stopband ripple factor obtained from

${\epsilon }_s=\mathrm{1/}\sqrt{{10}^{\frac{As}{10}}-1}$ ,

The transfer function of the Chebyshev2 filter is:

$\frac{Y\left(s\right)}{U\left(s\right)}=K{\displaystyle \prod _{k=1}^N}\frac{(s-s_{zk})}{(s-s_{pk})}$ ,

$s_{pk}$ are the poles obtained from:

$s_{pk}=\frac{-{\omega }_s}{D_k}\left(sinh\left(\frac{arsinh\left(\frac{1}{{\epsilon }_p}\right)}{N}\right)sin\left({\theta }_k\right)+jcosh\left(\frac{arcosh\left(\frac{1}{{\epsilon }_p}\right)}{N}\right)cos\left({\theta }_k\right)\right)$ ,

where ${\theta }_k=\frac{\pi }{2}\frac{2k-1}{n}$ , $D_k={\sinh }^2\left(\frac{arsinh\left(\frac{1}{{\epsilon }_s}\right)}{N}\right){\sin }^2\left({\theta }_k\right)+{\cosh }^2\left(\frac{arsinh\left(\frac{1}{{\epsilon }_s}\right)}{N}\right){\cos }^2\left({\theta }_k\right)$ , for $k=1,2,\mathrm{...}N$

$s_{zk}$ are the zeros obtained from:

$s_{zk}=\frac{\pm j{\omega }_s}{cos([2k-1]\pi /(2N))}$ , for $k=1,2,...[N/2]$

If the order of the filter is odd, the last zero (zero of the highest frequency) is infinity.

$K$ is the gain of the transfer function, equal to ${\epsilon }_s$ if the order is even, or ${\epsilon }_s{\omega }_{s}N$ if the order is odd.

Figure 7 shows the 4th Order Elliptic Filter with a 3 dB passband ripple ( $A_p$ ), a 40 dB attenuation in the stopband ( $A_s$ ) and with 500Hz cutoff frequency obtained from the component. Elliptic filter has ripples both in the passband and in the stopband, but it also has the steepest roll-off.

The magnitude response of the Elliptic filter is:

$R_n\left(\omega ,{\omega }_s,{\omega }_c,{\epsilon }_s,{\epsilon }_p\right)=sn[b{\mathrm{sn}}^{\mathrm{-1}}(\omega /{\omega }_c,{\omega }_c/{\omega }_s)+q{X}_1,{\epsilon }_s{\epsilon }_p]$

where $sn\left(x,\tau \right)$ is the Jacobian elliptic sine of $x$ with modulus $\tau$ , where ${sn}^{\mathrm{-1}}\left(x,\tau \right)$ is the inverse Jacobian elliptic sine of $x$ with modulus $\tau$ , $q$ is zero if the order is odd and unity if order is even, $X_1$ is the complete elliptic integral of first kind of ${\tau }_1={\epsilon }_s{\epsilon }_p$ , ${\epsilon }_s$ and ${\epsilon }_p$ are ripple parameters as same as for the Chebyshev1 and Chebyshev2 filters, and $b$ can be obtained from:

$b=\frac{N{X}_1}{X_2}$ ,

where $N$ is the filter order and $X_2$ is the complete elliptic integral of first kind of ${\tau }_2={\omega }_p/{\omega }_s$ .

Ports

• In (in)
  • The signal that the filtering should be applied to.
    • Supported types: real.
    • Vector support: no.

• Out (out)
  • Filtered input signal.
    • Supported types: real.
    • Vector support: no.
• Cutoff frequency (fc_input)
  • Port enabled when the frequency source of the component is external (except for the Elliptic filter).
    • Supported types: real.
    • Vector support: no.

Properties

• Filter type
  • Choose the type of the low-pass filter.

• Filter order
  • Type in the order of the low-pass filter.

• Cutoff frequency source
  • Choose the source of the cutoff frequency for the low-pass filter. Option "Internal" displays Cutoff frequency property. Option "External" creates a new input port on the component (except for the Elliptic filter).

• Cutoff frequency
  • Type in the cutoff frequency of the low-pass filter. Property enabled when the frequency source of the component is internal.

• Passband attenuation
  • Type in the passband attenuation of the low-pass filter. It is the attenuation at cut-off frequency for both Chebyshev1 and Elliptic filters.

• Stopband attenuation
  • Type in the stopband attenuation of the low-pass filter. It is the attenuation at stopband frequency for both Chebyshev2 and Elliptic filters.

• Execution rate
  • Type in the desired signal processing execution rate. This value must be compatible with other signal processing components of the same circuit: the value must be a multiple of the fastest execution rate in the circuit. There can be up to four different execution rates, but they must all be multiple of the basic simulation timestep. To specify the execution rate, you can use either decimal (e.g. 0.001) or exponential values (e.g. 1e-3) in seconds. Alternatively, you can type in ‘inherit’ in which case the component will be assigned execution rate based on the execution rate of the components it is receiving input from.

• Bode plot
  • Button which displays the magnitude Bode plot of the created filter. If any property is defined through the namespace, or the cutoff frequency is defined externally, a Generic Bode plot with valid parameters will be shown (properties which cannot be evaluated will be substituted with unity or other generic values). Also, a warning will occur in the output console.