Next Article in Journal
Numerical Solution of Variable-Order Fractional Differential Equations Using Bernoulli Polynomials
Next Article in Special Issue
Rational Approximations of Arbitrary Order: A Survey
Previous Article in Journal
Existence and Uniqueness Results of Coupled Fractional-Order Differential Systems Involving Riemann–Liouville Derivative in the Space Wa+γ1,1(a,b)×Wa+γ2,1(a,b) with Perov’s Fixed Point Theorem
Previous Article in Special Issue
CMOS OTA-Based Filters for Designing Fractional-Order Chaotic Oscillators
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

FPAA-Based Realization of Filters with Fractional Laplace Operators of Different Orders

by
Stavroula Kapoulea
1,
Costas Psychalinos
1,* and
Ahmed S. Elwakil
2,3,4
1
Department of Physics, Electronics Laboratory, University of Patras, GR-26504 Rio Patras, Greece
2
Department of Electrical and Computer Engineering, University of Sharjah, Sharjah 27272, United Arab Emirates
3
Nanoelectronics Integrated Systems Center (NISC), Nile University, Giza 16453, Egypt
4
Department of Electrical and Computer Engineering, University of Calgary, Calgary, AB T2N 1N4, Canada
*
Author to whom correspondence should be addressed.
Fractal Fract. 2021, 5(4), 218; https://doi.org/10.3390/fractalfract5040218
Submission received: 10 September 2021 / Revised: 6 November 2021 / Accepted: 11 November 2021 / Published: 13 November 2021
(This article belongs to the Special Issue Fractional-Order Circuit Theory and Applications)

Abstract

:
A simple and direct procedure for implementing fractional-order filters with transfer functions that contain Laplace operators of different fractional orders is presented in this work. Based on a general fractional-order transfer function that describes fractional-order low-pass, high-pass, band-pass, band-stop and all-pass filters, the introduced concept deals with the consideration of this function as a whole, with its approximation being performed using a curve-fitting-based technique. Compared to the conventional procedure, where each fractional-order Laplace operator of the transfer function is individually approximated, the main offered benefit is the significant reduction in the order of the resulting rational function. Experimental results, obtained using a field-programmable analog array device, verify the validity of this concept.

1. Introduction

Fractional-order filters have thrust the scientific field of filtering into a new era, making it applicable in a wide range of interdisciplinary fields, such as control systems, bio-impedance spectroscopy, acoustic applications and filter applications [1]. The insertion of one or more fractional exponents in a filter function provides additional degrees of freedom that allow precise adjustment of the filter’s characteristics, including the cutoff frequency and the slope of the transition from the pass band to stop band.
A general transfer function, able to describe fractional-order low-pass (LP), high-pass (HP), band-pass (BP), band-stop (BS) and all-pass (AP) filters, is given in (1)
H ( s )   =   K H P τ s α + β + K B P τ s β + K L P τ s α + β + τ s β + 1 ,
where α , β 0 , 1 are fractional orders; s is the Laplace operator; τ is a characteristic time constant; and K L P , K H P and K B P are real coefficients denoting the gain factors of the LP, HP and BP filters, respectively. The BS filter is obtained by the sum of the LP and HP filter functions, and the AP filter by the subtraction of the LP from the HP filter function [2].
The implementation of the filter functions, generated from (1), cannot be directly performed using conventional elements. Therefore, the following alternative methods are exploited:
(a)
Employment of fractional-order capacitors, which substitute the conventional (i.e., integer-order) capacitors in the topologies of integer-order standard filters. Due to the absence of commercial availability of such elements, they are approximated by RC schemes, such as the Foster and Cauer networks. The price paid for the offered quick design procedure is the absence of on-the-fly tuning of the filter’s characteristics, making it suitable only for cases of filters with pre-defined type and frequency characteristics [3,4,5].
(b)
Approximation of the fractional-order Laplace operators in (1) using appropriate tools, such as Continued Fraction Expansion, Oustaloup, Matsuda and Carlson [6], and then substitution of the resulting rational function approximations of these operators. The derived filter function is an integer-order polynomial ratio that can be implemented utilizing classical filter design techniques, such as cascade connection of elementary filter stages and multi-feedback structures. This design procedure is more complicated than the previous one, but it offers the advantage of electronic tuning capability of the derived filter structures, considering that active elements with electronically controlled characteristics (e.g., transconductance) are utilized [2,7,8,9,10,11]. However, this approach is efficient only for approximations of filter functions with a single fractional order. In the case of functions, such as the one in (1), the order of the derived rational function will be at least double the applied order of approximation [2].
The contribution made in this work is summarized as follows:
(a)
Instead of individually approximating the fractional-order Laplace operators, the magnitude and phase frequency responses of the whole filter function in (1) are fitted by a suitable integer-order transfer function derived through the employment of specialized Symbolic Math Toolbox™ built-in functions.
(b)
The order of the derived rational function is always equal to the order of the employed approximation, and thus, this approach offers the advantage of minimizing the circuit complexity, compared to the aforementioned conventional approach.
The paper is organized as follows: an investigation of the problem regarding the approximation of the general filter function following the conventional and the proposed procedures is discussed in Section 2, and a design example for comparison purposes is given in Section 3. The validity of the suggested procedure was experimentally evaluated using a field-programmable analog array (FPAA) device [12] and the results are presented in Section 4.

2. Proposed Design Procedure

Using any of the continued fraction expansion, Oustaloup, Matsuda or Carlson methods, for the approximation of the Laplace operators in (1), the transfer function in (2) is derived this way:
( τ s ) r A n s n + A n 1 s n 1 + + A 1 s + A 0 s n + B n 1 s n 1 + + B 1 s + B 0 ,
where r   =   { β , α + β } is the order of the operator; n is the order of the applied approximation; and A i ( i   =   0 , 1 , n ) , B j ( j   =   0 , 1 , n 1 ) are real, positive coefficients.
By approximating each one of the fractional terms τ s β and τ s α + β of (1) using (2), the resulting approximated filter function will have the form of a rational, integer-order polynomial ratio, is that given in (3):
H a p p r o x ( s ) C p s p + C p 1 s p 1 + + C 1 s + C 0 s q + D q 1 s q 1 + + D 1 s + D 0 ,
where p , q is the order of the ratio’s polynomials and C k , ( k   =   0 , 1 , p ) , D l , ( l   =   0 , 1 , q 1 ) are real coefficients.
In the case that α + β 0 , 1 , the nth-order approximation of the fractional terms τ s β and τ s α + β is performed by transfer functions of the form of (2) with order equal to n. Therefore, the transfer function in (3) will be of order 2 n , due to the different values of the coefficients of the corresponding intermediate transfer functions used for approximating the Laplace operators. In the case that α + β 1 , 2 , the approximation of the term τ s α + β is performed through the utilization of the product τ s ( α + β ) 1 · τ s , and the resulting orders of the transfer function in (3) are summarized in Table 1.
The main derivation is that, for a given order of approximation n , the order of the final filter function will be at least the double, which means increased circuit complexity. The proposed concept suggests the exploitation of a curve-fitting-based approximation technique, which is based on the Sanathanan–Koerner (SK) least square iterative method, for approaching the filter function in (1) [8,9]. The application of this technique is simple, considering that it is supported by MATLAB software and specialized Symbolic Math Toolbox™ built-in functions. The main advantage of this method, compared to conventional approximation tools, is the capability of approaching a whole transfer function, leading to an expression that always has the simple form of (3). In addition, this method has already been compared with other relevant frequency response-based approximation algorithms, and it has proven to be more efficient in terms of accuracy [8,13].
The steps of the algorithm are:
step #1:
Extraction of the frequency response data of (1) within a specific frequency range of interest, a process that is performed using the frd built-in function.
step #2:
Approximation of the obtained data, based on the fitfrd built-in function, which forms the state-space model of the data for a given order of approximation n .
step #3:
Conversion of the model to an integer-order transfer function of the form in (3) using the ss2tf built-in function.
Therefore, the final approximated filter function, for an nth-order approximation, will be of order n , n independent of the type of the filter, and this is the most attractive benefit offered by the proposed procedure. As the transfer function in (3) will have a constant form, the selection of the type of filter, as well as the tuning of its characteristics, will be performed through the adjustment of the values of the coefficients C k and D l . This results in a versatile structure whose behavior can be externally determined without any change in its core. In other words, following the proposed procedure, an electronically adjustable or programmable structure can be employed for simultaneously offering the elementary LP, HP, BP, BS and AP filter functions with minimized active component count, and therefore, reduced circuit complexity, compared to the structures presented in [2].

3. Design Example and Comparative Results

The initial function in (1) was approximated using the proposed procedure, based on the curve-fitting technique, for order n   =   5 within the frequency range f   =   10 , 100 k Hz for the filter cases given in Table 2. In order to prove the efficiency of the proposed concept, a comparison with the Oustaloup approximation method of the same order n   =   5 was also performed.
The final filter function, derived following the proposed approximation process, is that given in (3) with p , q   =   5 , 5 for all types of filters. This is not the case following the Oustaloup method, where, by consulting the Table 1, the order of the obtained function is p , q   =   10 , 11 for the LP and BP types, and p , q   =   11 , 11 for the HP, BS and AP types. Therefore, the Oustaloup method leads to double the hardware compared to that required for the curve-fitting-based method. In addition, the fact that the derived filter function in the case of the proposed procedure has a common form in all filter types makes it possible to implement different types of filters with variable order using the same controllable core.
The efficiency of each approximation method was evaluated through the gain and phase frequency responses of Figure 1, obtained using the MATLAB software, where the Oustaloup results are indicated by dashed-dotted lines and the curve-fitting-based results by solid lines. Considering the associated ideal responses, which are indicated by dashed lines, the difference between the accuracy of the two methods is evident. In particular, the maximum integral absolute error (IAE) on the gain for the LP, HP, BP, BS and AP filter types was about 1.94 × 10 4 dB·Hz and in the phase 1.01 × 10 7 °·Hz in the case of the Oustaloup method. The corresponding values, in the case of the proposed approximation procedure, were 4.22 × 10 3 dB·Hz and 1.43 × 10 6 °·Hz. Consequently, the proposed approximation process offers not only a considerably compact transfer function that can be implemented with the requirement of minimum circuitry, but also a sufficiently accurate approximation over the entire frequency range of interest.

4. Experimental Results

The proposed concept was experimentally evaluated using FPAA hardware, and in particular, an Anadigm QUAD AN231K04 development kit, which is controlled through the AnadigmDesigner 2 EDA software [14,15,16,17,18]. Employing an integrator-based follow-the-leader feedback (FLF) structure, which is depicted in Figure 2, considering a 5th-order approximation, the transfer function that describes this functional block diagram is the following:
H F L F ( s )   =   K 5 s 5 + K 4 s 4 + + K 1 s + K 0 s 5 + 1 τ 1 s 4 + + 1 τ 1 τ 2 τ 4 s + 1 τ 1 τ 2 τ 5 ,
and the time constants ( τ ) and scaling factors ( K ) are calculated by equating the coefficients of (3) and (4). Therefore, the design equations, considering that D 5   =   1 , becomes as follows:
K i   =   C i D i , i   =   0 , 1 , 5 τ i + 1   =   D 5 i D 4 i , i   =   0 , 1 , 4 .
For demonstration purposes, six of the cases, given in Table 2, were examined: (i) LP filter with α , β   =   0.3 , 0.5 , (ii) HP filter with α , β   =   0.5 , 0.3 , (iii) BP filter with α , β   =   0.7 , 0.4 , (iv) BS filter with α , β   =   0.7 , 0.4 , (v) AP filter with α , β   =   0.3 , 0.5 and (vi) AP filter with α , β   =   0.6 , 0.8 . The coefficients of the derived approximated filter functions are provided in Table 3.
Using the values of Table 3 and the design equations in (5), the calculated values of time constants and scaling factors are summarized in Table 4.
The implementation of the FLF structure in Figure 2, using the FPAA hardware, was based on the configuration that is demonstrated in the Anadigm design of Figure 3. Two FPAA chips were required when Integrator, GainHold and SumDiff Configurable Analogue Blocks (CABs) were utilized to implement the integration and scaling stages, and the chip clock was set at 2 MHz. The experimental setup that was used is shown in Figure 4.
The experimental gain and phase were measured using an Agilent 4395A spectrum analyzer within the frequency range f   =   10 , 100 k Hz, and the resulting responses are given in Figure 5a for the LP and HP filter types; in Figure 5b for the BP and BS types; and in Figure 5c for the AP filter type, indicated by triangle symbols, along with the respective results obtained from the curve-fitting-based approximation (solid lines) and the theory (dashed lines).
Indicatively, a screenshot of the analyzer, showing the gain frequency response during the measurement of the BP type of filter for orders α , β   =   0.7 , 0.4 , is demonstrated in Figure 6. From these plots, the versatile features of the implemented FPAA-based structure were verified, making the presented procedure attractive from the design flexibility point of view.
The experimental values of the critical frequencies for each case are tabulated in Table 5, with the respective theoretically calculated values being given between parentheses, verifying the accurate behavior of the implemented filter functions.

5. Conclusions

The proposed concept for the implementation of fractional-order filter transfer functions with Laplace operators of multiple fractional orders, based on the exploitation of a curve-fitting-based approximation process, provides simultaneously an efficient and compact realization. The reduction in the order of the derived approximated rational transfer function was significant, compared to the case where conventional approximation methods were followed. The provided design example shows that, according to the presented procedure, the derived order of the transfer function for a 5th-order approximation was p , q   =   5 , 5 for all filter types with a maximum error of about 10 % in the frequency range 10 , 100 k Hz. When employing the Oustaloup approximation, the order of the function was p , q   =   10 , 11 for the LP and BP types and 11 , 11 for the HP, BS and AP types. The maximum observed error in the aforementioned frequency range was about 50 % . Another resulting benefit is that all types of filters with variable orders can be implemented by the same core, making the presented structure a generalized filter scheme. This was experimentally verified through the utilization of the Anadigm QUAD AN231K04 development kit.
The employed algorithm offers the advantage of quick design procedure in the sense that the whole transfer function is considered, instead of the intermediate fractional-order Laplace terms. In addition, it approximates the fractional-order function within a frequency range defined by the user, as in the case of the Oustaloup method. The same also holds for similar MATLAB software Symbolic Math Toolbox™ built-in functions, such as tfest. On the other hand, a drawback is that these curve-fitting approximation methods are only available in MATLAB software. Another option could be employing the Padé tool, which performs approximation of the transfer function around a center frequency specified by the user, as in the case of the continued fraction expansion method. A comparison among all the aforementioned methods showed that the Sanathanan–Koerner (SK) least square iterative method employed in this work outperforms any of the other approximation tools [8,13].
The presented design procedure is versatile and can be also applied for the approximation of various types of fractional-order transfer functions, including control systems [19,20] and biological/biomedical systems [21], and these are undergoing research exploitation steps.

Author Contributions

Conceptualization, C.P. and A.S.E.; methodology, S.K.; software, S.K.; validation, S.K.; formal analysis, S.K.; investigation, S.K. and C.P.; writing—original draft preparation, S.K. and C.P.; writing—review and editing, C.P., S.K. and A.S.E.; project administration, C.P. 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.

Acknowledgments

This research is co-financed by Greece and the European Union (European Social Fund-ESF) through the Operational Programme “Human Resources Development, Education and Lifelong Learning” in the context of the project “Strengthening Human Resources Research Potential via Doctorate Research-2nd Cycle” (MIS-5000432), implemented by the State Scholarships Foundation (IKY).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
APAll pass
BPBand pass
BSBand stop
CABConfigurable analogue block
FOFractional order
FBDFunctional block diagram
FLFFollow the leader feedback
FPAAField-programmable analog array
HPHigh pass
IAEIntegral absolute of error
LPLow pass

References

  1. 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]
  2. 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]
  3. Sladok, O.; Koton, J.; Kubanek, D.; Dvorak, J.; Psychalinos, C. Pseudo-differential (2 + α)-order Butterworth frequency filter. IEEE Access 2021, 1. [Google Scholar] [CrossRef]
  4. Biswal, K.; Tripathy, M.C.; Kar, S. Performance analysis of fractional order low-pass filter. In Advances in Intelligent Computing and Communication; Springer: Berlin/Heidelberg, Germany, 2020; pp. 224–231. [Google Scholar]
  5. Mijat, N.; Jurisic, D.; Moschytz, G.S. Analog Modeling of Fractional-Order Elements: A Classical Circuit Theory Approach. IEEE Access 2021, 9, 110309–110331. [Google Scholar] [CrossRef]
  6. Hamed, E.M.; AbdelAty, A.M.; Said, L.A.; Radwan, A.G. Effect of different approximation techniques on fractional-order KHN filter design. Circuits Syst. Signal Process. 2018, 37, 5222–5252. [Google Scholar] [CrossRef]
  7. 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]
  8. Kapoulea, S.; Psychalinos, C.; Elwakil, A.S. Double exponent fractional-order filters: Approximation methods and realization. Circuits Syst. Signal Process. 2021, 40, 993–1004. [Google Scholar] [CrossRef]
  9. Kapoulea, S.; Psychalinos, C.; Elwakil, A.S. Power law filters: A new class of fractional-order filters without a fractional-order Laplacian operator. AEU-Int. J. Electron. Commun. 2021, 129, 153537. [Google Scholar] [CrossRef]
  10. 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]
  11. Mahata, S.; Herencsar, N.; Kubanek, D. Optimal Approximation of Fractional-Order Butterworth Filter Based on Weighted Sum of Classical Butterworth Filters. IEEE Access 2021, 9, 81097–81114. [Google Scholar] [CrossRef]
  12. Hasler, J.; Shah, S. An SoC FPAA based programmable, ladder-filter based, linear-phase analog filter. IEEE Trans. Circuits Syst. Regul. Pap. 2020, 68, 592–602. [Google Scholar] [CrossRef]
  13. Xue, D. Fractional-Order Control Systems; de Gruyter: Berlin, Germany, 2017. [Google Scholar]
  14. Anadigm. AN231E04 dpASP: The AN231E04 dpASP Dynamically Reconfigurable Analog Signal Processor. Available online: https://anadigm.com/an231e04.asp (accessed on 5 November 2021).
  15. Freeborn, T.J.; Maundy, B.; Elwakil, A.S. Field programmable analogue array implementation of fractional step filters. IET Circuits Devices Syst. 2010, 4, 514–524. [Google Scholar] [CrossRef]
  16. Freeborn, T.J.; Maundy, B. Incorporating FPAAs into laboratory exercises for analogue filter design. Int. J. Electr. Eng. Educ. 2013, 50, 188–200. [Google Scholar] [CrossRef]
  17. Muñiz-Montero, C.; García-Jiménez, L.V.; Sánchez-Gaspariano, L.A.; Sánchez-López, C.; González-Díaz, V.R.; Tlelo-Cuautle, E. New alternatives for analog implementation of fractional-order integrators, differentiators and PID controllers based on integer-order integrators. Nonlinear Dyn. 2017, 90, 241–256. [Google Scholar] [CrossRef]
  18. 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]
  19. Lino, P.; Maione, G.; Stasi, S.; Padula, F.; Visioli, A. Synthesis of fractional-order PI controllers and fractional-order filters for industrial electrical drives. IEEE/CAA J. Autom. Sin. 2017, 4, 58–69. [Google Scholar] [CrossRef] [Green Version]
  20. Azar, A.T.; Radwan, A.G.; Vaidyanathan, S. Fractional Order Systems: Optimization, Control, Circuit Realizations and Applications; Academic Press: Cambridge, MA, USA, 2018. [Google Scholar]
  21. Mohsen, M.; Said, L.A.; Madian, A.H.; Radwan, A.G.; Elwakil, A.S. Fractional-Order Bio-Impedance Modeling for Interdisciplinary Applications: A Review. IEEE Access 2021, 9, 33158–33168. [Google Scholar] [CrossRef]
Figure 1. Gain and phase frequency responses, derived from 5th-order Oustaloup (dashed-dotted lines) and 5th-order curve-fitting-based (solid lines) approximation methods, for the (a) LP and HP, (b) BP and BS and (c) AP filter cases given in Table 2, along with the corresponding ideal responses (dashed lines).
Figure 1. Gain and phase frequency responses, derived from 5th-order Oustaloup (dashed-dotted lines) and 5th-order curve-fitting-based (solid lines) approximation methods, for the (a) LP and HP, (b) BP and BS and (c) AP filter cases given in Table 2, along with the corresponding ideal responses (dashed lines).
Fractalfract 05 00218 g001
Figure 2. FLF block diagram for implementing LP, HP, BP, BS and AP filter functions.
Figure 2. FLF block diagram for implementing LP, HP, BP, BS and AP filter functions.
Fractalfract 05 00218 g002
Figure 3. FPAA-based realization of the functional block diagram in Figure 2.
Figure 3. FPAA-based realization of the functional block diagram in Figure 2.
Fractalfract 05 00218 g003
Figure 4. Experimental setup for evaluating the performances of the filters.
Figure 4. Experimental setup for evaluating the performances of the filters.
Fractalfract 05 00218 g004
Figure 5. Experimental gain and phase frequency responses (triangle symbols), along with the corresponding approximated (solid lines) and ideal (dashed lines) responses, for (a) LP and HP, (b) BP and BS and (c) AP filter cases.
Figure 5. Experimental gain and phase frequency responses (triangle symbols), along with the corresponding approximated (solid lines) and ideal (dashed lines) responses, for (a) LP and HP, (b) BP and BS and (c) AP filter cases.
Fractalfract 05 00218 g005
Figure 6. Screenshot of the Agilent 4395A spectrum analyzer, showing the gain frequency response during the measurement of the BP type of filter for orders α , β   =   0.7 , 0.4 .
Figure 6. Screenshot of the Agilent 4395A spectrum analyzer, showing the gain frequency response during the measurement of the BP type of filter for orders α , β   =   0.7 , 0.4 .
Fractalfract 05 00218 g006
Table 1. Order of the functions in (3) for different types of filters using conventional approximation methods.
Table 1. Order of the functions in (3) for different types of filters using conventional approximation methods.
FilterOrder
α + β 0 , 1 α + β 1 , 2
LP p , q   =   2 n , 2 n p , q   =   2 n , 2 n + 1
HP p , q   =   2 n + 1 , 2 n + 1
BP p , q   =   2 n , 2 n + 1
BS p , q   =   2 n + 1 , 2 n + 1
AP p , q   =   2 n + 1 , 2 n + 1
Table 2. Values of parameters in (1) for realizing various filter cases.
Table 2. Values of parameters in (1) for realizing various filter cases.
FilterParameter
α , β K LP , K HP , K BP τ
LP 0.3 , 0.5 1 , 0 , 0 22.91 μ s
0.3 , 0.8 69.97 μ s
HP 0.5 , 0.3 0 , 1 , 0 1.106 m s
0.8 , 0.3 362 μ s
BP 0.4 , 0.7 0 , 0 , 1 186.7 μ s
0.7 , 0.4 135.7 μ s
BS 0.4 , 0.7 1 , 1 , 0 205 μ s
0.7 , 0.4 123.5 μ s
AP 0.3 , 0.5 1 , 1 , 0 184.8 μ s
0.6 , 0.8 547.8 μ s
Table 3. Values of the coefficients in (3) derived from the 5th-order approximation of the initial function in (1) for various filter types.
Table 3. Values of the coefficients in (3) derived from the 5th-order approximation of the initial function in (1) for various filter types.
CoefficientFilter
LPHPBPBSAP
0 . 3 , 0 . 5 0 . 5 , 0 . 3 0 . 7 , 0 . 4 0 . 7 , 0 . 4 0 . 3 , 0 . 5 0 . 6 , 0 . 8
C 0 1.812 × 10 20 2.383 × 10 15 4.68 × 10 16 6.689 × 10 17 1.026 × 10 19 1.146 × 10 19
C 1 4.361 × 10 17 1.811 × 10 14 1.101 × 10 15 5.767 × 10 15 3.783 × 10 16 6.374 × 10 15
C 2 1.317 × 10 14 8.018 × 10 11 2.705 × 10 12 7.219 × 10 12 8.754 × 10 12 2.488 × 10 12
C 3 6.778 × 10 9 5.37 × 10 8 1.139 × 10 9 1.886 × 10 9 2.109 × 10 9 2.18 × 10 9
C 4 5.138 × 10 4 6.26 × 10 4 2.169 × 10 4 1.237 × 10 5 1.555 × 10 5 1.386 × 10 5
C 5 0.02145 0.9617 0.01929 0.9802 0.8324 0.9842
D 0 1.871 × 10 20 1.296 × 10 17 5.653 × 10 17 7.264 × 10 17 1.11 × 10 19 1.165 × 10 19
D 1 5.069 × 10 17 1.296 × 10 15 5.51 × 10 15 7.152 × 10 15 5.894 × 10 16 1.961 × 10 16
D 2 1.966 × 10 14 2.196 × 10 12 8.418 × 10 12 1.059 × 10 13 4.652 × 10 13 1.243 × 10 13
D 3 1.65 × 10 10 8.576 × 10 8 2.753 × 10 9 3.223 × 10 9 7.759 × 10 9 3.067 × 10 9
D 4 3.144 × 10 5 7.504 × 10 4 1.368 × 10 5 1.473 × 10 5 2.532 × 10 5 1.528 × 10 5
Table 4. Values of time constants and scaling factors in (4), calculated using the values of Table 3 and the design equations in (5).
Table 4. Values of time constants and scaling factors in (4), calculated using the values of Table 3 and the design equations in (5).
CoefficientFilter
LPPBPBSAP
0 . 3 , 0 . 5 0 . 5 , 0 . 3 0 . 7 , 0 . 4 0 . 7 , 0 . 4 0 . 3 , 0 . 5 0 . 6 , 0 . 8
K 0 0.969 0.018 0.083 0.921 0.924 0.984
K 1 0.86 0.14 0.2 0.806 0.642 0.325
K 2 0.67 0.365 0.321 0.682 0.188 0.2
K 3 0.41 0.626 0.414 0.585 0.272 0.711
K 4 0.163 0.834 0.158 0.84 0.614 0.907
K 5 0.021 0.962 0.019 0.98 0.832 0.984
τ 1 3.18 ( μ s)13.32 ( μ s)7.31 ( μ s)6.79 ( μ s)3.95 ( μ s)6.546 ( μ s)
τ 2 19.05 ( μ s)87.5  ( μ s)49.7 ( μ s)45.71 ( μ s)32.63 ( μ s)49.81 ( μ s)
τ 3 83.92 ( μ s)390.5 ( μ s)327 ( μ s)304.2 ( μ s)166.8 ( μ s)246.6 ( μ s)
τ 4 387.9 ( μ s)1.69 (ms)1.53 (ms)1.78 (ms)789.25 ( μ s)634.1 ( μ s)
τ 5 2.71 (ms)10 (ms)9.75 (ms)9.85 (ms)5.31 (ms)1.68 (ms)
Table 5. Critical frequencies for the filter cases in Figure 5, along with the theoretical values given between parentheses.
Table 5. Critical frequencies for the filter cases in Figure 5, along with the theoretical values given between parentheses.
FrequencyFilter
LPHPBPBSAP
f c / f p e a k (kHz)1 (1)1 (1)1 (1)1 (1)1 (1)
f l o w (Hz)142.5 (140.7)173.8 (174.1)
f h i g h (kHz)4.36 (4.32)3.48 (3.37)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kapoulea, S.; Psychalinos, C.; Elwakil, A.S. FPAA-Based Realization of Filters with Fractional Laplace Operators of Different Orders. Fractal Fract. 2021, 5, 218. https://doi.org/10.3390/fractalfract5040218

AMA Style

Kapoulea S, Psychalinos C, Elwakil AS. FPAA-Based Realization of Filters with Fractional Laplace Operators of Different Orders. Fractal and Fractional. 2021; 5(4):218. https://doi.org/10.3390/fractalfract5040218

Chicago/Turabian Style

Kapoulea, Stavroula, Costas Psychalinos, and Ahmed S. Elwakil. 2021. "FPAA-Based Realization of Filters with Fractional Laplace Operators of Different Orders" Fractal and Fractional 5, no. 4: 218. https://doi.org/10.3390/fractalfract5040218

APA Style

Kapoulea, S., Psychalinos, C., & Elwakil, A. S. (2021). FPAA-Based Realization of Filters with Fractional Laplace Operators of Different Orders. Fractal and Fractional, 5(4), 218. https://doi.org/10.3390/fractalfract5040218

Article Metrics

Back to TopTop