Next Article in Journal
Attention-Oriented Deep Multi-Task Hash Learning
Next Article in Special Issue
FPGA-Based Chaotic Image Encryption Using Systolic Arrays
Previous Article in Journal
Private and Public Digitalization and Economic Resilience during COVID-19 Pandemic
Previous Article in Special Issue
Multi-Stage Ensemble-Based System for Glaucomatous Optic Neuropathy Diagnosis in Fundus Images
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

A 1 + α Order Generalized Butterworth Filter Structure and Its Field Programmable Analog Array Implementation

1
Electronics Laboratory, Department of Physics, University of Patras, GR-26504 Rio Patras, Greece
2
Department of Electrical and Computer Engineering, University of Sharjah, Sharjah P.O. Box 27272, United Arab Emirates
3
Nanoelectronics Integrated Systems Center (NISC), Nile University, Giza 12677, Egypt
4
Department of Electrical and Software Engineering, University of Calgary, Calgary, AB T2N 1N4, Canada
*
Author to whom correspondence should be addressed.
Electronics 2023, 12(5), 1225; https://doi.org/10.3390/electronics12051225
Submission received: 30 January 2023 / Revised: 22 February 2023 / Accepted: 27 February 2023 / Published: 3 March 2023

Abstract

:
Fractional-order Butterworth filters of order 1 + α (0 < α < 1) can be implemented by a unified structure, using the method presented in this paper. The main offered benefit is that the cutoff frequencies of the filters are fully controllable using a very simple method and, also, various types of filters (e.g., low-pass, high-pass, band-pass, and band-stop) could be realized. Thanks to the employment of a Field Programmable Analog Array device, the implementation of the introduced method is fully reconfigurable, in the sense that various types of filter functions as well as their order are both programmable.

1. Introduction

Fractional-order filters offer a fine adjustment of the slope of the attenuation from the pass-band to the stop-band, and this is the result of the fact that it is equal to ± 20 · ( n + α ) dB/dec, where n is the integer part of the order and 0 < α < 1 is its non-integer counterpart [1,2,3,4]. In addition, they offer the scaling of the cutoff frequencies due to their dependence on both the pole frequency and the order of the filter. In the case of filters of order 0 < α < 1, there are simple expressions describing the dependence of the cutoff frequencies on the aforementioned parameters. This is not the case for filters with the order 1 + α , where complicated non-linear equations must be solved to find the value of the cutoff frequency. In addition, the realization of well-known filter functions such as the Butterworth transfer functions is an even more difficult task [5,6].
A significant research effort has been dedicated to overcoming this obstacle, where appropriate cost functions error minimization metaheuristic and genetic algorithms were reported. A first attempt for approximating fractional-order Butterworth filters has been performed in [7], where a polynomial fitting-based method was used for expressing the coefficients of the fractional-order transfer function as a function of the order of the filter. This method is oriented to the realization of low-pass filters, and the high-pass filter functions have been derived through a frequency transformation. In [8], the approximation of the fractional-order Butterworth low-pass filters is performed using the pole placement in the W-plane. In [9], the Flower Pollination Algorithm (FPA) with an appropriate cost function is employed for approximating Butterworth low-pass filters, while in [10], a nature-inspired optimization technique called the Gravitational Search Algorithm (GSA) has been employed for the same purpose. In [11], the Symbiotic Organisms Search (SOS) algorithm is presented for approximating Butterworth high-pass filters. In [12], a genetic algorithm is introduced for approximating Butterworth low-pass filters.
From the above, it is readily obtained that the solutions which are available in the literature are based on the employment of relatively complicated optimization algorithms. In addition, they do not offer a simultaneous approximation of different types of Butterworth filter functions, such as the low-pass, high-pass, band-pass, and band-stop.
In the present work, a simple procedure for realizing various types of 1 + α order Butterworth transfer functions, with fully controllable cutoff frequencies, is introduced. It is based on the fit of the frequency response magnitude data by a minimum-phase state-space model, using a Log-Chebyshev magnitude design. The main important benefits are as follows: (a) The procedure is independent of the type of the filter; low-pass, high-pass, band-pass, and band-stop filter functions can be implemented. (b) The resulting structure has the capability of implementing all the aforementioned filter functions without altering its core.
This paper is organized as follows: The description of the basic characteristics of the fractional filters of order 1 + α is performed in Section 2. The proposed procedure is introduced in Section 3, where a possible implementation of the resulting approximation functions using a Field Programmable Array (FPAA) device is also demonstrated. The performance of the designed filters is experimentally verified in Section 4.

2. Fractional-Order Filters of Order 1 + α (0 < α < 1)

Let us consider the transfer function of a low-pass filter of order 0 < α < 1 and pole frequency ω 0 = 1 / τ , with the variable τ being a time constant
H ( s ) = 1 τ s α + 1 .
Setting s α = ω α · [ cos ( 0.5 α π ) + j sin ( 0.5 α π ) ] in (1), the following expression of the magnitude response is derived:
H ( ω ) = 1 ω ω 0 2 α + 2 ω ω 0 α cos α π 2 + 1 1 / 2 .
The cutoff or half-power frequency of the filter ( ω h ), defined as the frequency where there is a 3 dB drop in the magnitude from its maximum value, is determined using (2) and setting H ( ω ) = 0.707 . The resulting expression is given by (3)
ω h = ω 0 1 + cos 2 α π 2 cos α π 2 M M 1 / α .
In a similar way, the behavior of a high-pass filter is described by (4) and (5) and the associated cutoff frequency is given by the expression in (6)
H ( s ) = τ s α τ s α + 1 ,
H ( ω ) = ω ω 0 α ω ω 0 2 α + 2 ω ω 0 α cos α π 2 + 1 1 / 2 ,
ω h = ω 0 1 + cos 2 α π 2 + cos α π 2 1 / α .
Inspecting (3) and (6), it is readily obtained that, having available the pole frequency and the order of the filter, the cutoff frequency can be fully determined.
In the case of a 1 + α order low-pass filter, its transfer function is given by (7)
H ( s ) = 1 τ s 1 + α + 1 Q τ s α + 1 ,
with Q being a factor equal to the quality factor of the corresponding second-order filter. The half-power frequency is the solution of the following non-linear equation:
ω h ω 0 2 1 + α 2 ω h ω 0 1 + α sin α π 2 + 1 Q 2 ω h ω 0 2 α + 2 Q ω h ω 0 α cos α π 2 1 = 0 .
The corresponding non-linear equation for obtaining the half-power frequency of a high-pass filter, described by (9),
H ( s ) = τ s 1 + α τ s 1 + α + 1 Q τ s α + 1 ,
is provided by (10)
ω h ω 0 2 1 + α + 2 ω h ω 0 1 + α sin α π 2 1 Q 2 ω h ω 0 2 α 2 Q ω h ω 0 α cos α π 2 1 = 0 .
For band-pass and band-stop filters described by (11) and (12), respectively,
H ( s ) = 1 Q τ s α τ s 1 + α + 1 Q τ s α + 1 ,
H ( s ) = τ s 1 + α + 1 τ s 1 + α + 1 Q τ s α + 1 .
The peak (center) frequency is calculated by solving the differential equation:
ω H ( ω ) ω = ω p e a k = 0 , while the associated lower ( ω l o w ) and upper ( ω h i g h ) cutoff frequencies are calculated by setting H ( ω ) ω = ω l o w , h i g h = 0.707 · H ( ω ) ω = ω p e a k .
From the above, it is clear that the determination of the characteristic frequencies in the case of 1 + α order filters is a difficult procedure, which is the reason for using the aforementioned algorithmic methods. This obstacle can be overcome through the utilization of the curve fitting-based procedure which will be described in the next section.

3. Proposed Procedure for Approximating 1 + α Order Filters

3.1. Derivation of the Rational Integer-Order Approximation Function

According to the proposed method, the starting point is the consideration of the gain responses of the filters. For example, the corresponding expressions of low-pass, high-pass, band-pass, and band-stop Butterworth-type filters are those in (13)–(16).
H L P , 1 + α ( ω ) = 1 1 + ω ω 0 2 ( 1 + α ) ,
H H P , 1 + α ( ω ) = 1 1 + ω 0 ω 2 ( 1 + α ) ,
H B P , 1 + α ( ω ) = 1 1 + ω 0 2 ω 2 ω · ( ω h i g h ω l o w ) 2 ( 1 + α ) ,
H B S , 1 + α ( ω ) = 1 1 + ω · ( ω h i g h ω l o w ) ω 0 2 ω 2 2 ( 1 + α ) .
The next step is to collect the values of the gain of the filter transfer function of interest within the desired range ( ω m i n , ω m a x ). This can be performed through the utilization of the MATLAB frd built-in function. The last steps include the fit of the collected gain data with a minimum-phase state-space model using the Log-Chebyshev magnitude design, realized using the MATLAB built-in command fitmagfrd, and the derivation of the corresponding rational integer-order transfer function using the tf MATLAB command.
The resulting transfer function will have the form of (17)
H ( s ) = B n s n + B n 1 s n 1 + + B 1 s + B 0 s n + A n 1 s n 1 + + A 1 s + A 0 ,
with A i and B j i = 0 , 1 . . . n 1 , j = 0 , 1 . . . n being positive and real coefficients, and n being the order of the approximation.
Considering Butterworth low-pass and high-pass filters of order 1.2, 1.4, and 1.6, with cutoff frequency ω h =1 rad/s and approximated by a 4th-order approximation within the range of [10 1 , 101] rad/s, the resulting values of the coefficients A i and B j are summarized in Table 1 and Table 2. In the case of a band-pass filter, the same order of approximation is considered. While the peak and cutoff frequencies are ω p e a k = 1 rad/s, ω l o w = 0.5 rad/s, and ω h i g h = 2 rad/s, for the band-stop filter, the corresponding values are ω p e a k = 1 rad/s, ω l o w = 0.2 rad/s, and ω h i g h = 5 rad/s. The values of the coefficients are provided in Table 3 and Table 4.
The derived plots are demonstrated in Figure 1 and Figure 2, accompanied by their associated error plots.

3.2. A Possible Implementation

A possible realization of the transfer function in (17) is using the Follow-the-Leader Feedback (FLF) structure in Figure 3, where the realized transfer function is given by (18) [13].
H ( s ) = G 4 s 4 + G 3 τ 1 s 3 + G 2 τ 1 τ 2 s 2 + G 1 τ 1 τ 2 τ 3 s + G 0 τ 1 τ 2 τ 3 τ 4 s 4 + 1 τ 1 s 3 + 1 τ 1 τ 2 s 2 + 1 τ 1 τ 2 τ 3 s + 1 τ 1 τ 2 τ 3 τ 4 .
Equalizing the coefficients of (17) and (18), the resulting design equations are given by (19)
τ i + 1 = A 4 i A 3 i ( i = 0 . . . 3 ) G j = B j A j ( j = 0 . . . 4 ) .
The implementation of the integration blocks in Figure 3 can be performed using a variety of active blocks, including Operational Amplifiers (op-amps), Second-Generation Current Conveyors (CCIIs), and Current Feedback Operational Amplifiers (CFOAs). The derived implementations are active-RC structures which do not offer the electronic tunability of the characteristics of the filters. The employment of the Operational Transconductance Amplifiers (OTAs) as active elements is a solution for overcoming this obstacle due to the fact that the transconductance parameter of the OTA can be electronically adjusted by an appropriate dc current/voltage. The price paid for this achievement is the reduced linearity of the resulting filter structures, and this is caused by the small-signal nature of the transconductance [14]. Another option could be the utilization of a digitally programmable device [15], such as the FPAA AN231E04 device provided by Anadigm [16].
Denormalizing (17) at ω h = 150 krad/s and taking into account that the integration constants ( M i ) of the FPAA device are formed as M i = 10 6 / τ i ( i = 1 . . . 4 ), (having units of μ s 1 ), their values are summarized in Table 5, Table 6, Table 7 and Table 8, where the values of the scaling factors G j (j = 0 ...4) are also given. The clock frequency is equal to f c l k = 2 MHz. In addition, according to the specifications of the FPAA, the realizable values of the integration constants are in the range [0.02, 17.7] μ s 1 , while the weight factors of the summation stages for realizing the scaling factors G j are in the range [0.01, 8.83]. Inspecting Table 5, Table 6, Table 7 and Table 8, it is derived that there are values of M i and G j that are not realizable. As the range of the realizable values of the gain stages is [0.01, 100], the problem can be solved by following the technique presented in Figure 4. According to this, the desired value of the gain scaling factor G j is implemented as G j = G j · G j , making the values of G j to be in the realizable range. In the FPAA implementation presented in the next section, this will be performed by adding an extra gain stage to each one of the inputs of the summation stages, and both of them will have a value equal to G j , establishing a total gain equal to G j in the path of the associated signal. A similar approach will be followed for the integration constants.

4. Experimental Results

Using the Anadigm Designer® 2 EDA software, the resulting design is depicted in Figure 5, while the experimental setup is demonstrated in Figure 6a. The FPAA device with the accompanied printed circuit board (PCB) for implementing the required single-to-differential and differential-to-single voltage conversion is demonstrated in Figure 6b. The corresponding circuitry of the interface for performing a single-to-differential conversion of the input signal is depicted in Figure 7a, where the corresponding one for performing the differential-to-single conversion of the output signal is demonstrated in Figure 7b. Taking into account the terminal properties of the CFOA, described by the following formulas, υ X = υ Y , i X = i Z , υ Z = υ O , the expressions in (20) and (21) are readily obtained for the topology in Figure 7a.
υ i n , p = υ i n 2 + E 2 ,
υ i n , n = υ i n 2 + E 2 .
The value of the dc voltage source is E = 3 V, and this is originated from the fact that the common-mode voltage of the FPAA is equal to 1.5 V.
The topology in Figure 7b implements the following input–output relationship
υ o u t = υ o u t , p υ o u t , n ,
The υ o u t , p and υ o u t , n are given by the formulas υ o u t , p = υ o u t / 2 + E / 2 and υ o u t , n = υ o u t / 2 + E / 2 . The AD844 discrete component is chosen for implementing the required CFOAs [17] while the value of the basic resistor is R = 10 k Ω .
The frequency responses of the Butterworth low-pass and high-pass filters, which are obtained using the HP4395A Network/Impedance analyzer, are demonstrated in Figure 8a,b, respectively. The measured values of the cutoff frequencies were 151.1 rad/s for the low-pass filter and 147.8 rad/s for the high-pass filter, close to the theoretical value of 150 krad/s. The experimental results of the band-pass and band-stop filters are provided in Figure 9a,b. The measured value of the lower cutoff frequency of the band-pass filters was 74.7 krad/s, with the theoretical value being 75 krad/s. The corresponding values of the band-stop filters were 11.8 krad/s and 298.1 krad/s, with the associated theoretical values being 12 krad/s and 300 krad/s, respectively.
These slight deviations from the theory are mainly caused by the tolerances of the resistors employed in the interfacing stages as well as by the truncation of the values of the time constants and the scaling factors performed by the FPAA device.
The time-domain behavior of the filters has been evaluated by stimulating the low-pass and high-pass filters (order equal to 1.4) by a 1Vp-p sinusoidal signal with a frequency equal to their cutoff frequency. The waveforms, obtained using a DSO6034A digital oscilloscope, are demonstrated in Figure 10a,b. The values of the corresponding gains were −3.06 dB and −3.4 dB, confirming their time-domain performance. In the case of the band-pass filters, the stimulation of a filter of an order equal to 1.4 by sinusoidal inputs of the same amplitude as in the previous case and of a frequency equal to the peak frequency results in the waveforms in Figure 10c, where the expected gain of 0 dB is confirmed.

5. Conclusions

A versatile procedure for approximating various types of Butterworth filters with fully controllable frequency characteristics is introduced in this work. The algorithm is based on the built-in commands frd, fitmagfrd, and tf of the MATLAB software, leading to a rational integer-order transfer function. This can be implemented using the well-known classical filter design techniques, such as multi-feedback or cascaded structures. The choice of the technique, as well as of the utilized active element, is up to the designer, performing a trade-off between the design complexity and electronic tunability.
The experimental results derived using an FPAA device confirm the validity of the proposed concept. Future research steps include the application of the presented concept for implementing other types of filter functions, such as Chebyshev filters [18,19,20,21], Bessel filters [22,23,24], and elliptic filters [25].

Author Contributions

Conceptualization, C.P., J.N. and A.S.E.; methodology, J.N. and C.P.; software, J.N.; validation, J.N; formal analysis, J.N.; investigation, J.N. and C.P.; writing—original draft preparation, C.P.; writing—review and editing, C.P. and A.S.E.; project administration, C.P. and A.S.E. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
AMSAustria Mikro Systeme
CCIISecond-Generation Current Conveyor
CFOACurrent Feedback Operational Amplifier
CMOSComplimentary Metal-Oxide Semiconductor
EDAElectronic Design Automation
FLFFollow-the-Leader Feedback
FPAFlower Pollination Algorithm
FPAAField Programmable Analog Array
GSAGravitational Search Algorithm
MOSMetal-Oxide Semiconductor
OP-AMPOperational Amplifier
OTAOperational Transconductance Amplifier
PCBPrinted Circuit Board
RCResistor Capacitor
SOSSymbiotic Organisms Search

References

  1. Bertsias, P.; Khateb, F.; Kubanek, D.; Khanday, F.A.; Psychalinos, C. Capacitorless digitally programmable fractional-order filters. AEU-Int. J. Electron. Commun. 2017, 78, 228–237. [Google Scholar] [CrossRef]
  2. Jerabek, J.; Sotner, R.; Dvorak, J.; Polak, J.; Kubanek, D.; Herencsar, N.; Koton, J. Reconfigurable fractional-order filter with electronically controllable slope of attenuation, pole frequency and type of approximation. J. Circuits Syst. Comput. 2017, 26, 1750157. [Google Scholar] [CrossRef] [Green Version]
  3. Dvorak, J.; Langhammer, L.; Jerabek, J.; Koton, J.; Sotner, R.; Polak, J. Synthesis and analysis of electronically adjustable fractional-order low-pass filter. J. Circuits Syst. Comput. 2018, 27, 1850032. [Google Scholar] [CrossRef]
  4. Langhammer, L.; Dvorak, J.; Sotner, R.; Jerabek, J.; Bertsias, P. Reconnection–Less reconfigurable low–Pass filtering topology suitable for higher–Order fractional–Order design. J. Adv. Res. 2020, 25, 257–274. [Google Scholar] [CrossRef] [PubMed]
  5. Freeborn, T.J. Comparison of (1 + α) Fractional-Order Transfer Functions to Approximate Lowpass Butterworth Magnitude Responses. Circuits, Syst. Signal Process. 2016, 35, 1983–2002. [Google Scholar] [CrossRef]
  6. Tsirimokou, G.; Psychalinos, C.; Elwakil, A. Design of CMOS Analog Integrated Fractional-Order Circuits: Applications in Medicine and Biology; Springer: Berlin/Heidelberg, Germany, 2017. [Google Scholar]
  7. Freeborn, T.J.; Maundy, B.; Elwakil, A.S. Field programmable analogue array implementation of fractional step filters. IET Circ. Devices Syst. 2010, 4, 514–524. [Google Scholar] [CrossRef]
  8. Mishra, S.K.; Upadhyay, D.K.; Gupta, M. Approximation of Fractional-Order Butterworth Filter Using Pole-Placement in W-Plane. IEEE Trans. Circ. Syst. II Express Briefs 2021, 68, 3229–3233. [Google Scholar] [CrossRef]
  9. Mahata, S.; Kar, R.; Mandal, D. Optimal design of fractional-order Butterworth filter with improved accuracy and stability margin. In Fractional-Order Modeling of Dynamic Systems with Applications in Optimization, Signal Processing and Control; Elsevier: Amsterdam, The Netherlands, 2022; pp. 293–321. [Google Scholar]
  10. Mahata, S.; Saha, S.K.; Kar, R.; Mandal, D. Optimal design of fractional order low pass Butterworth filter with accurate magnitude response. Digit. Signal Process. 2018, 72, 96–114. [Google Scholar] [CrossRef]
  11. Mahata, S.; Kar, R.; Mandal, D. Optimal fractional-order highpass Butterworth magnitude characteristics realization using current-mode filter. AEU-Int. J. Electron. Commun. 2019, 102, 78–89. [Google Scholar] [CrossRef]
  12. Mahata, S.; Banerjee, S.; Kar, R.; Mandal, D. Revisiting the use of squared magnitude function for the optimal approximation of (1 + α)-order Butterworth filter. AEU-Int. J. Electron. Commun. 2019, 110, 152826. [Google Scholar] [CrossRef]
  13. Tsirimokou, G.; Psychalinos, C.; Elwakil, A.S.; Salama, K.N. Electronically tunable fully integrated fractional-order resonator. IEEE Trans. Circ. Syst. II Express Briefs 2017, 65, 166–170. [Google Scholar] [CrossRef] [Green Version]
  14. Tsirimokou, G.; Psychalinos, C.; Elwakil, A.S. Fractional-order electronically controlled generalized filters. Int. J. Circuit Theory Appl. 2017, 45, 595–612. [Google Scholar] [CrossRef]
  15. Tlelo-Cuautle, E.; Pano-Azucena, A.D.; Guillén-Fernández, O.; Silva-Juárez, A. Analog/Digital Implementation of Fractional Order Chaotic Circuits and Applications; Springer: Berlin/Heidelberg, Germany, 2020. [Google Scholar]
  16. Anadigm. AN231E04 dpASP: The AN231E04 dpASP Dynamically Reconfigurable Analog Signal Processor. Available online: https://anadigm.com/an231e04.asp (accessed on 3 January 2023).
  17. Devices, A. MHz 2000 V/μs Monolithic Op Amp with Quad Low Noise AD844 Data Sheet, Revision G; One Technology Way: Noorwood, MA, USA, 2009; p. 60. [Google Scholar]
  18. Freeborn, T.; Maundy, B.; Elwakil, A.S. Approximated fractional order Chebyshev lowpass filters. Math. Probl. Eng. 2015, 2015, 832468. [Google Scholar] [CrossRef]
  19. Freeborn, T.J.; Elwakil, A.S.; Maundy, B. Approximated fractional-order inverse Chebyshev lowpass filters. Circuits Syst. Signal Process. 2016, 35, 1973–1982. [Google Scholar] [CrossRef]
  20. AbdelAty, A.M.; Soltan, A.; Ahmed, W.A.; Radwan, A.G. Fractional order Chebyshev-like low-pass filters based on integer order poles. Microelectron. J. 2019, 90, 72–81. [Google Scholar] [CrossRef]
  21. Daryani, R.; Aggarwal, B.; Gupta, M. Design of Fractional-Order Chebyshev Low-Pass Filter for Optimized Magnitude Response Using Metaheuristic Evolutionary Algorithms. Circuits Syst. Signal Process. 2022, 1–31. [Google Scholar] [CrossRef]
  22. Ćoza, A.; Županović, V.; Vlah, D.; Jurišić, D. Group delay of fractional n+ α-order Bessel filters. In Proceedings of the IEEE 2020 43rd International Convention on Information, Communication and Electronic Technology (MIPRO), Opatija, Croatia, 28 September–2 October 2020; pp. 163–168. [Google Scholar]
  23. Soni, A.; Gupta, M. Analysis and design of optimized fractional order low-pass Bessel filter. J. Circuits Syst. Comput. 2021, 30, 2150035. [Google Scholar] [CrossRef]
  24. Soni, A.; Gupta, M. Designing of Fractional Order Bessel Filter using Optimization Techniques. Int. J. Electron. Lett. 2022, 10, 71–86. [Google Scholar] [CrossRef]
  25. Kubanek, D.; Freeborn, T.J.; Koton, J.; Dvorak, J. Validation of fractional-order lowpass elliptic responses of (1 + α)-order analog filters. Appl. Sci. 2018, 8, 2603. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Frequency responses of the approximated Butterworth (a) low-pass filters and (b) high-pass filters of orders 1.2, 1.4, and 1.6.
Figure 1. Frequency responses of the approximated Butterworth (a) low-pass filters and (b) high-pass filters of orders 1.2, 1.4, and 1.6.
Electronics 12 01225 g001
Figure 2. Frequency responses of the approximated Butterworth (a) band-pass filters and (b) band-stop filters of orders 1.2, 1.4, and 1.6.
Figure 2. Frequency responses of the approximated Butterworth (a) band-pass filters and (b) band-stop filters of orders 1.2, 1.4, and 1.6.
Electronics 12 01225 g002
Figure 3. Follow-the-Leader Feedback block diagram for implementing the transfer function in (17).
Figure 3. Follow-the-Leader Feedback block diagram for implementing the transfer function in (17).
Electronics 12 01225 g003
Figure 4. Concept for enhancing the range of the implementable scaling factors.
Figure 4. Concept for enhancing the range of the implementable scaling factors.
Electronics 12 01225 g004
Figure 5. FPAA configuration for realizing Butterworth low-pass, high-pass, band-pass, and band-stop filters.
Figure 5. FPAA configuration for realizing Butterworth low-pass, high-pass, band-pass, and band-stop filters.
Electronics 12 01225 g005
Figure 6. (a) Experimental setup for evaluating the performance of the implemented filters, (b) FPAA device along with the associated interface.
Figure 6. (a) Experimental setup for evaluating the performance of the implemented filters, (b) FPAA device along with the associated interface.
Electronics 12 01225 g006
Figure 7. Circuitry for performing (a) single-to-differential conversion of the input signal and (b) circuitry for performing differential-to-single conversion of the output signal.
Figure 7. Circuitry for performing (a) single-to-differential conversion of the input signal and (b) circuitry for performing differential-to-single conversion of the output signal.
Electronics 12 01225 g007
Figure 8. Experimental frequency responses of the Butterworth (a) low-pass filters and (b) high-pass filters of orders 1.2, 1.4, 1.6, realized using the FPAA configuration in Figure 5.
Figure 8. Experimental frequency responses of the Butterworth (a) low-pass filters and (b) high-pass filters of orders 1.2, 1.4, 1.6, realized using the FPAA configuration in Figure 5.
Electronics 12 01225 g008
Figure 9. Experimental frequency responses of the Butterworth (a) band-pass filters and (b) band-stop filters of orders 1.2, 1.4, 1.6, realized using the FPAA configuration in Figure 5.
Figure 9. Experimental frequency responses of the Butterworth (a) band-pass filters and (b) band-stop filters of orders 1.2, 1.4, 1.6, realized using the FPAA configuration in Figure 5.
Electronics 12 01225 g009
Figure 10. Experimental input and output waveforms of (a) low-pass filter, (b) high-pass filter, and (c) band-pass filter of order equal to 1.4.
Figure 10. Experimental input and output waveforms of (a) low-pass filter, (b) high-pass filter, and (c) band-pass filter of order equal to 1.4.
Electronics 12 01225 g010aElectronics 12 01225 g010b
Table 1. Values of coefficients for approximating Butterworth low-pass filters of order 1.2, 1.4, 1.6 ( ω h = 1 rad/s).
Table 1. Values of coefficients for approximating Butterworth low-pass filters of order 1.2, 1.4, 1.6 ( ω h = 1 rad/s).
Coefficient Order
1.21.41.6
B 4 0.0058670.002280.001421
B 3 0.61670.33360.1768
B 2 5.1493.8412.674
B 1 6.856.966.071
B 0 1.7041.7691.64
A 3 7.1967.0146.537
A 2 12.6612.410.97
A 1 8.6068.8777.976
A 0 1.7011.771.644
Table 2. Values of coefficients for approximating Butterworth high-pass filters of order 1.2, 1.4, 1.6 ( ω h = 1 rad/s).
Table 2. Values of coefficients for approximating Butterworth high-pass filters of order 1.2, 1.4, 1.6 ( ω h = 1 rad/s).
Coefficient Order
1.21.41.6
B 4 1.0020.99910.9974
B 3 4.0283.9323.693
B 2 3.0282.171.627
B 1 0.36260.18840.1075
B 0 0.003450.0012880.0008641
A 3 5.065.0154.852
A 2 7.4467.0036.676
A 1 4.2313.9623.976
A 0 0.5880.5650.6083
Table 3. Values of coefficients for approximating Butterworth band-pass filters of order 1.2, 1.4, 1.6 ( ω p e a k = 1 rad/s, ω l o w = 0.5 rad/s, ω h i g h = 2 rad/s).
Table 3. Values of coefficients for approximating Butterworth band-pass filters of order 1.2, 1.4, 1.6 ( ω p e a k = 1 rad/s, ω l o w = 0.5 rad/s, ω h i g h = 2 rad/s).
Coefficient Order
1.21.41.6
B 4 0.0027950.0023580.004134
B 3 0.98570.6320.4071
B 2 4.8284.0163.325
B 1 0.98570.6320.4071
B 0 0.0027950.0023580.004134
A 3 4.4483.733.118
A 2 6.7725.9395.246
A 1 4.4483.733.118
A 0 111
Table 4. Values of coefficients for approximating Butterworth band-stop filters of order 1.2, 1.4, 1.6 ( ω p e a k = 1 rad/s, ω l o w = 0.2 rad/s, ω h i g h = 5 rad/s).
Table 4. Values of coefficients for approximating Butterworth band-stop filters of order 1.2, 1.4, 1.6 ( ω p e a k = 1 rad/s, ω l o w = 0.2 rad/s, ω h i g h = 5 rad/s).
Coefficient Order
1.21.41.6
B 4 1.2711.8152.373
B 3 0.19480.19590.1668
B 2 2.5423.634.747
B 1 0.19480.19590.1668
B 0 1.2711.8152.373
A 3 8.07413.5918.94
A 2 4.528.95719.05
A 1 8.07413.5918.94
A 0 111
Table 5. Values of the integration constants and scaling factors of the FPAA device for approximating Butterworth low-pass filters of order 1.2, 1.4, 1.6 ( ω h = 150 krad/s).
Table 5. Values of the integration constants and scaling factors of the FPAA device for approximating Butterworth low-pass filters of order 1.2, 1.4, 1.6 ( ω h = 150 krad/s).
Coefficient Order
1.21.41.6
M 1 1.0791.0520.980
M 2 0.2640.2650.252
M 3 0.1020.1070.109
M 4 0.0300.0300.031
G 0 1.0020.9990.997
G 3 0.7960.7840.761
G 2 0.4010.3100.244
G 3 0.0860.0480.027
G 4 0.0060.0020.001
Table 6. Values of the integration constants and scaling factors of the FPAA device for approximating Butterworth high-pass filters of order 1.2, 1.4, 1.6 ( ω h = 150 krad/s).
Table 6. Values of the integration constants and scaling factors of the FPAA device for approximating Butterworth high-pass filters of order 1.2, 1.4, 1.6 ( ω h = 150 krad/s).
Coefficient Order
1.21.41.6
M 1 0.7590.7520.728
M 2 0.2200.2090.206
M 3 0.0850.0850.089
M 4 0.0210.0210.023
G 0 0.0060.0020.001
G 3 0.0860.0470.027
G 2 0.4070.3100.244
G 3 0.7960.7840.761
G 4 1.0020.9990.997
Table 7. Values of the integration constants and scaling factors of the FPAA device for approximating Butterworth band-pass filters of order 1.2, 1.4, 1.6 ( ω p e a k = 150 krad/s, ω l o w = 75 krad/s, ω h i g h = 300 krad/s).
Table 7. Values of the integration constants and scaling factors of the FPAA device for approximating Butterworth band-pass filters of order 1.2, 1.4, 1.6 ( ω p e a k = 150 krad/s, ω l o w = 75 krad/s, ω h i g h = 300 krad/s).
Coefficient Order
1.21.41.6
M 1 0.6670.5590.468
M 2 0.2280.2390.252
M 3 0.0980.0940.089
M 4 0.0340.0400.048
G 0 0.0030.0020.004
G 3 0.2210.1690.131
G 2 0.7130.6760.634
G 3 0.2230.1690.131
G 4 0.0030.0020.004
Table 8. Values of the integration constants and scaling factors of the FPAA device for approximating Butterworth band-stop filters of order 1.2, 1.4, 1.6 ( ω p e a k = 60 krad/s, ω l o w = 12 krad/s, ω h i g h = 300 krad/s).
Table 8. Values of the integration constants and scaling factors of the FPAA device for approximating Butterworth band-stop filters of order 1.2, 1.4, 1.6 ( ω p e a k = 60 krad/s, ω l o w = 12 krad/s, ω h i g h = 300 krad/s).
Coefficient Order
1.21.41.6
M 1 0.4840.8161.1362
M 2 0.0330.0390.060
M 3 0.2070.0920.060
M 4 0.0070.0040.003
G 0 1.2711.8152.373
G 3 0.0240.0140.009
G 2 0.5620.4050.2492
G 3 0.0240.0140.009
G 4 1.2711.8152.373
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Nako, J.; Psychalinos, C.; Elwakil, A.S. A 1 + α Order Generalized Butterworth Filter Structure and Its Field Programmable Analog Array Implementation. Electronics 2023, 12, 1225. https://doi.org/10.3390/electronics12051225

AMA Style

Nako J, Psychalinos C, Elwakil AS. A 1 + α Order Generalized Butterworth Filter Structure and Its Field Programmable Analog Array Implementation. Electronics. 2023; 12(5):1225. https://doi.org/10.3390/electronics12051225

Chicago/Turabian Style

Nako, Julia, Costas Psychalinos, and Ahmed S. Elwakil. 2023. "A 1 + α Order Generalized Butterworth Filter Structure and Its Field Programmable Analog Array Implementation" Electronics 12, no. 5: 1225. https://doi.org/10.3390/electronics12051225

APA Style

Nako, J., Psychalinos, C., & Elwakil, A. S. (2023). A 1 + α Order Generalized Butterworth Filter Structure and Its Field Programmable Analog Array Implementation. Electronics, 12(5), 1225. https://doi.org/10.3390/electronics12051225

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop