On the Design of Power Law Filters and Their Inverse Counterparts

: This paper presents the optimal modeling of Power Law Filters (PLFs) with the low-pass (LP), high-pass (HP), band-pass (BP), and band-stop (BS) responses by means of rational approximants. The optimization is performed for three different objective functions and second-order ﬁlter mother functions. The formulated design constraints help avoid placement of the zeros and poles on the right-half s -plane, thus, yielding stable PLF and inverse PLF (IPLF) models. The performances of the approximants exhibiting the fractional-step magnitude and phase responses are evaluated using various statistical indices. At the cost of higher computational complexity, the proposed approach achieved improved accuracy with guaranteed stability when compared to the published literature. The four types of optimal PLFs and IPLFs with an exponent α of 0.5 are implemented using the follow-the-leader feedback topology employing AD844AN current feedback operational ampliﬁers. The experimental results demonstrate that the Total Harmonic Distortion achieved for all the practical PLF and IPLF circuits was equal or lower than 0.21%, whereas the Spurious-Free Dynamic Range also exceeded 57.23 and 54.72 dBc, respectively.


Introduction
The theoretical concepts of fractional calculus [1][2][3], which generalized differ-integral operators, have led to significant developments in circuit theory, signal processing, control theory, bio-impedance modeling, etc. [4][5][6][7][8].Fractional-order (FO) filters are considered as the generalization of the traditional filters [9].This is due to the ability of the FO filters to achieve any roll-off rate [10]; in contrast, an integer-order filter can only achieve a rolloff at −20 log 10 n decibels/decade (dB/dec), where n is an integer [11].FO analog filter transfer functions are generally realized from the integer-order filters by substitution of the Laplacian operator s with the non-integer Laplacian operator s α , where α ∈ (0, 1).The frequency-domain transfer function of s α is given by ( 1): where j = √ −1 and ω is the angular frequency in radians per second (rad/s).Since s α is an irrational function, various rational approximations based on series truncation, frequency-domain curve-fitting, pole-zero placement, optimization techniques, etc., have been reported [12][13][14][15].The impedance characteristics of the operator s α may be practically realized using the FO elements (also known as the fractance devices or the constant phase elements) [16][17][18].Due to the unavailability of the commercial FO device, their behavior may be emulated using the passive and active circuits [19][20][21][22].
Recent works have demonstrated the generalization of the Butterworth [23], Chebyshev [24], inverse Chebyshev [25], and elliptic filters [26] to the FO domain.Another design strategy that involves the approximation of FO filter characteristics using the integer order transfer function was also reported in the literature [27].The integer order approximant can be realized using the field programmable analog array [28], voltage mode operational amplifier [29], switched capacitor [30], operational transconductance amplifier (OTA) [31], and current feedback operational amplifier (CFOA) [32].
The application of numerical and metaheuristic optimization algorithms for the approximation of FO filters has gained traction in recent years.The modeling of FO RLC filter and low-pass filter transfer functions of the form 1/(s + 1) α using classical optimization techniques has been reported [33,34].Numerical optimization methods were employed for the approximation of low-pass [35,36] and band-pass [37,38] filters exhibiting fractional-step behavior.The Nelder-Mead simplex [39], Cuckoo Search algorithm [40], and MATLAB-based optimization function fmincon [41,42] were employed to model the magnitude-frequency characteristics of the FO low-pass Butterworth filter.The performances of several nature-inspired algorithms were compared for the optimal rational approximation of the s α operator [43].The optimal design of a grounded FO inductor based on the generalized impedance converter was reported in [44].(1 + α) and (α + β) order, where α, β ∈ (0, 1), low-pass Bessel filter characteristics were optimally approximated using the Interior Search algorithm [45].
Since inverse filters yield the reciprocal frequency characteristic of the system that causes distortion during transmission or reception, these filters are widely used in communication systems to alleviate signal distortion [46].Inverse filters are also employed in acoustic systems [47], proportional integral derivative controllers [48], and digital filtering [49].Realizations of FO inverse filters using the operational amplifiers [50], CFOAs [51], operational trans-resistance amplifiers [52], and OTAs [53] have been recently reported.
A new class of FO filters, namely the Power Law Filter (PLF), was approximated based on a frequency-domain curve fitting using the Sanathanan-Koerner (S-K) leastsquare iterative method in [54].Applications of power law compensators and filters in achieving robust frequency compensation of process plants with uncertainties [55] and in bio-impedance modeling of fruits [56] have also been exemplified.The transfer function, magnitude-frequency, and phase-frequency relationships for these filters exhibiting the low-pass (LP), high-pass (HP), band-pass (BP), and band-stop (BS) characteristics are presented in Table 1 [54].In the PLFs, a FO exponent term α is introduced in the transfer function of the second-order filter (referred to as the mother filter function).Consequently, the second-order filter characteristics may now be considered as a particular case of the PLFs when α = 1.

Type Transfer Function Magnitude Phase
Low-pass This paper presents the optimal modeling of PLFs exhibiting the LP, HP, BP, and BS characteristics.Comparative studies concerning the modeling performance based on three different objective function formulations are presented.A non-parametric statistical null hypothesis test is employed to investigate the similarity in the design performances between the objective functions.The proposed strategy incorporates constraints to ensure a stable PLF approximant.Additional constraints are also incorporated to avoid the placement of zeros in the right-half s-plane.This allows the attainment of inverse PLF (IPLF) characteristics through inversion of the optimal PLF models.Although the S-K method [54] is simpler and computationally superior compared to the proposed technique, the suggested approach demonstrates the following advantages: (i) The stability of the PLF approximant in S-K method [54] is governed by the stability of the mother filter function.While this guarantees a stable rational approximant, however, the zeros of the filter transfer function may lie on the right-half of s (RHS)plane.For example, while the zeros of the BPPLFs and BSPLFs designed using the S-K method lie on the left-half s-plane; however, the LPPLF and HPPLF models for α = 0.7 have zeros located at {+183.0053, -23.6922, -8.6177, -1.7485} and {+0.0055, -0.5718, -0.1161, -0.0422}, respectively.It may be noted that for the S-K method-based LPPLFs and HPPLFs with α ∈ [0.51, 0.99], one zero lies on the RHS-plane.Consequently, inverting such a transfer function will lead to an unstable inverse-LPPLF (ILPPLF) and inverse-HPPLF (IHPPLF) model.In contrast, the proposed approach can guarantee the generation of stable designs for both PLFs and IPLFs.Hence, this paper also presents the design of IPLFs that has not yet been reported in the literature.(ii) The modeling accuracy of the proposed optimal PLF approximants, as justified by the Mean Absolute Relative Error (MARE) metric, is significantly better (particularly for the LP, HP, and BP-types) in comparison to the S-K method.
To demonstrate the practical efficacy, discrete component-based circuit realization using CFOAs as active elements for the proposed PLFs and their inverse counterparts with α = 0.5 is conducted.The experimental results reveal excellent agreement with the theoretical magnitude-and phase-frequency behavior for the PLFs and magnitude responses for the IPLFs.The time-domain and Fast Fourier Transformation (FFT) characteristics of the practical filters are also investigated.
In this paper, Section 2 presents the proposed optimization problem formulation and PLF/IPLF design strategy.The modeling accuracy of the proposed method is investigated using MATLAB simulations in Section 3. Statistical studies concerning the hypothesis test are also presented in this section.Practical circuit implementation and experimental results (magnitude-frequency, phase-frequency, time-domain, and FFT) for the designed filters are presented in Section 4. Finally, the paper concludes in Section 5.

Proposed Technique
The proposed rational approximant of the PLF is modeled as per (2): where the coefficients of the numerator and denominator polynomials of H P (s) are denoted by a k (k = 0, 1, . . ., N) and b k (k = 0, 1, . . ., N -1), respectively; and N is the order of the filter.Figure 1 presents the flowchart of the proposed PLF design technique.The magnitude and phase characteristics of the PLFs may be approximated by determining the optimal values of the coefficients of H P (s) such that the error between the theoretical and proposed responses is minimized.For this purpose, three different objective functions, as defined by (3)- (5), are proposed: where L denotes the total number of logarithmically spaced sample points in the bandwidth of interest ω ∈ [ω min , ω max ] rad/s; the decision variables vector is represented by To achieve a stable design, the inequality constraints as given by ( 6) are incorporated in the proposed optimization method: where  6) with b k substituted with a k (k = 0, 1, • • • , N).Thus, the proposed formulation can be solved using a global search constrained optimization problem solver.The transfer function of the proposed IPLFs can be obtained as H I (s) = [H P (s)] −1 .As the value of α approaches 1, it is possible that the coefficients a N and a 0 of the optimal LPPLF and HPPLF, respectively, may attain a value of 0. The corresponding ILPPLF would possess one extra zero as compared to the number of poles, while the IHPPLF model will have a pole located at the origin of the s-plane.
To circumvent these particular issues, (i) the proposed ILPPLF transfer function may be represented as p.H I (s)/(s + p), where p is a large positive number such as 200, 500, or 1000.The accuracy of approximation for the IPLFs increases as p is increased; and (ii) for the case when a 0 = 0 arises for the HPPLF, a 0 may be replaced by a small, positive number (q), for instance, smaller than 0.005.This technique allows the pole of the IHPPLF to be shifted away from the origin towards the left-half of the s-plane, thus, avoiding potential instability issues.
A numerical optimizer requires a user-supplied initial point for the decision variables at the start of the optimization procedure.Subsequently, the optimization algorithm iteratively minimizes the objective function by varying the decision variables.For solving a constrained optimization problem, an additional task of the optimizer is to satisfy the design constraints (i.e., generate a feasible solution).At the end of a single run of the optimization routine (when the maximum number of function evaluations or iterations or function tolerance value is reached), the optimal values of decision variables are obtained.
The final solution quality of a numerical optimization algorithm may be influenced by the choice of the initial point.For solving the global search optimization problems, identifying an appropriate initial point may not be easy.To circumvent this problem, a standard technique is to independently execute the optimization routine several times (iter m ) with randomly chosen initial points in each run.Hence, iter m number of near-global optimal solution vectors can be generated in this process.The best decision variables vector (X best ) is selected as the one that achieves the smallest error (MARE in the present case) from the iter m solutions.The previously-mentioned strategy is employed in this paper for the optimal modeling of PLFs.

MATLAB Simulations and Performance Analysis
The optimization procedures to minimize the objective functions (3)-( 5) are implemented in MATLAB programming language using the function fmincon (algorithm: active-set) with the following parameter settings: maximum number of function evaluations = 50,000; maximum number of iterations = 50,000; and termination tolerance on the function value = 10 −10 .In each trial run, the initial point for the decision variables vector is randomly chosen from a uniform distribution in the interval [0, 50].For each design case, iter m = 100 independent trial runs of the optimization routine for each objective function are carried out.
Quantitative comparisons of the design accuracy are carried out based on the MARE metric as defined by ( 7): For demonstration purposes, the values of L, α, N, ω min , ω max , ω 0 , and Q are chosen as 1000, {0.3, 0.5, 0.7}, 4, 0.01 rad/s, 100 rad/s, 1 rad/s, and 1 √ 2, respectively.Detailed results for various other values of α are also available from the authors and can be shared with interested readers.

Comparisons Based on the MARE Metric
Statistical comparisons about the MARE metric were carried out to determine the average performance of the proposed optimization strategy for the design of PLFs using the three objective functions.The minimum (min), maximum (max), mean, and standard deviation (SD) indices for the various design orders are shown in Table 2. Graphical comparisons for the MARE attained for the designed PLFs are presented using boxplots in Figure 2. In the case of LPPLFs, we found that: (i) for α = 0.3, f 3 attained the best performance for all the statistical indices; (ii) the most accurate model (i.e., the design with the minimum value of MARE (MARE min )) for α = 0.5 was achieved by f 1 (1.11 × 10 −4 ), although f 2 yielded the best results for the other indices; and (iii) very similar performances for MARE min were achieved for α = 0.7, with f 3 (0.0068) attaining marginally better accuracy.
However, f 3 distinctly outperformed the other two objective functions regarding the max, mean, and SD values.For the HPPLFs, (i) f 3 achieved significantly better MARE min (0.0081) as compared to f 1 (0.0371) and f 2 (0.0380) for α = 0.3.However, the max, mean, and SD performances for f 3 were the worst among the three functions; (ii) with α = 0.5, f 1 yielded inferior performance about the min, mean, and SD indices, while the best values concerning MARE min and mean MARE were achieved by f 3 (1.20 × 10 −5 ) and f 2 (0.1516), respectively; and (iii) for α = 0.7, f 1 attained the most accurate model (MARE min = 0.0068), while f 2 yielded the best performance for all the other indices.

Filter
Index In the case of BPPLFs, it is revealed that: (i) for all the cases, the best performer about MARE min is f 3 ; (ii) regarding the mean and SD indices, f 2 attained the minimum value for α = 0.5 and f 3 for the other two orders; and (iii) the best results for the max MARE were achieved by f 2 (5.6081), f 2 (3.0704), and f 1 (3.2517) for α = 0.3, 0.5, and 0.7, respectively.Comparisons for the designed BSPLFs show that f 1 and f 2 achieved the best values for all the statistical indices with α = 0.3 and α = 0.7, respectively.In the case of α = 0.5, the most accurate model is attained by f 2 (MARE min = 0.0123), whereas f 1 yields superior performance for the other three indices.

Comparisons Based on the Wilcoxon Rank-Sum Hypothesis Test
Pair-wise comparisons based on the Wilcoxon rank-sum test for the PLFs designed using f 1 , f 2 , and f 3 are conducted to determine whether a statistically significant difference exists in the modeling performance [58].If no significant difference exists in the design accuracy, then a similar average-case modeling error performance can be expected from all the proposed objective functions.Consequently, all the three objective functions may exhibit similar robustness.
Thus, in terms of solution consistency, the objective functions demonstrate similar performance.For this purpose, the null hypothesis (H 0 : 'equality of medians for the MARE metric') is considered at a confidence level of 95%.The decision index is represented by H, where H = 0/1 indicates that H 0 cannot/can be rejected.Table 3 presents the Wilcoxon rank-sum test results along with the p-value (p-val) index.A smaller p-val indicates stronger evidence in favor of rejection of H 0 .
It is revealed that: (i) H 0 is rejected for all the pair-wise combinations concerning the LPPLF, HPPLF, and BSPLF with α = 0.3, 0.7, and 0.3, respectively; (ii) only one case (HPPLF with α = 0.3) exists where all the combinations may result in non-rejection of H 0 ; (iii) for comparisons involving all the three design orders for a particular pair, the hypothesis can be rejected for f 1 / f 2 with the BSPLF, f 2 / f 3 with the LPPLF, and f 1 / f 3 for the BPPLF and BSPLF, whereas H 0 may not be rejected for the BPPLF with f 2 / f 3 ; and (iv) across all the PLFs, 9, 8, and 5 cases out of 12 exist for f 1 / f 2 , f 1 / f 3 , and f 2 / f 3 , respectively, that lead to the rejection of H 0 .Overall, it may be concluded from the statistical analysis that, in terms of attaining a similar performance consistency for all the design cases, the three objective functions cannot be used interchangeably.
The optimal coefficients of the PLFs that achieve the smallest value of MARE for each of the three objective functions are presented in Appendix A. The MATLAB-simulated magnitude-and phase-frequency responses for the most accurate (least MARE min ) proposed LPPLFs, HPPLFs, BPPLFs, and BSPLFs are presented in Figures 3-6, respectively.The plots for the corresponding IPLFs are also illustrated in these figures.The values of p and q are considered as 1000 and 10 −6 for the ILPPLF and IHPPLF, respectively, for α = 0.5 and 0.7.We observed that all the design cases attained good agreement with the theoretical responses.The MARE min attained for the ILPPLF (α = 0.5) with p = {100, 200, 500, 1000} is {0.0919, 0.0439, 0.0164, 0.0079}.In the case of IHPPLF for α = 0.5 with q = {0.005,0.002, 0.001, 0.0001}, the MARE min is obtained as {0.0447, 0.0170, 0.0084, 0.0008}.These results confirm that the approximation accuracy of the ILPPLF and IHPPLF improves for larger values of p and smaller values of q, respectively.The MARE min attained for the proposed inverse BPPLFs (IBPPLFs) with α = {0.3,0.5, 0.7} is {0.0790, 0.0745, 0.0548}; the same for the inverse BSPLFs (IBSPLFs) is {0.0147, 0.0121, 0.0092}.

Comparison with the Literature
The proposed most accurate PLFs are compared with the published literature [54] based on the MARE metric, as shown in Table 4.The proposed designs outperform the reported models for all the considered cases by achieving the least MARE values.
The improvements in the percentage absolute relative error (100 ), where MARE Prop and MARE Rep are the MARE values achieved for the proposed model and [54], respectively, for the LPPLF, HPPLF, BPPLF, and BSPLF with α = {0.3,0.5, 0.7} are {47.40,58.58, 83.41}%, {48.73, 87.22, 83.41}%, {7.75, 16.85, 27.71}%, and {18.23, 10.87, 8.16}%, respectively.Graphical visualizations of these percentage improvements for all the proposed designs are shown using bar plots in Figure 7. Improvements over [54] for the proposed LPPLF and HPPLF are distinctly pronounced as compared to those of the BPPLF and BSPLF.Hence, the proposed technique may be considered as a more accurate modeling tool when compared against [54].Comparisons about the computational time (t c ) required by the S-K method [54] and the proposed technique for the design of PLFs with α = 0.5 are conducted in the following environment -CPU: Intel i3 @ 1.70 GHz, RAM: 2.0 GB, Operating System: Windows 7 (64 bit), and Software: MATLAB 2014a.The t c required by [54] for the design of LPPLF, HPPLF, BPPLF, and BSPLF is 1.601, 1.600, 1.665, and 1.663 s, respectively.For the proposed technique with objective functions { f 1 , f 2 , f 3 }, the t c (expressed in seconds) required to complete 100 iterations is, respectively, obtained as {254.These results show that no specific objective function can attain the smallest value of t c for all the designed PLFs.In terms of computational efficiency, the proposed strategy is inferior to the reported technique [54].However, since the PLF approximation is carried out offline, the higher t c of the proposed method may be traded-off in favor of achieving superior modeling accuracy compared to the S-K method [54].

Experimental Validation
CFOAs are popular analog signal processing integrated circuits that offer a high slew rate, gain-bandwidth decoupling, and smaller power consumption compared to operational amplifiers [59].Due to their versatility, CFOAs have been widely employed to implement FO filters [19,32,41,51,60].
Discrete components-based circuit realization of the proposed PLFs is carried out by employing the CFOA as an active element in a follow-the-leader feedback topology [19].The circuit diagram for the generalized PLF and IPLF implementation is presented in Figure 8.The total numbers of CFOAs, capacitors, and resistors required to realize the proposed Nth-order approximant are N + 2, N, and 3N + 4, respectively.The proposed circuit resulted in a reduced component count compared to [54] where the PLF implementation requires N + 4 operational amplifiers, N capacitors, and 3N + 8 resistors.Thus, irrespective of the value of N, the circuit reported in [54] requires two extra op-amps and four additional resistors as compared to the proposed one.The general transfer function of the proposed circuit is given by ( 8): As a representative, the circuit implementation steps for the proposed PLFs and IPLFs of LP, HP, BP, and BS types, whose transfer functions for α = 0.5 are given by ( 9)-( 16), are presented below.
The design steps are described as follows: Step 1: Set N = 4 in (8).Therefore, six CFOAs, four capacitors, and 16 resistors are required to construct the circuit.The transfer function for the CFOA-based circuit is given by ( 17): Step 2: Select the desired center frequency for the filter.For instance, a center frequency of 1 kHz is used here.
Step 3: Set the values of R, R F , R out , and R in .Note that the resistor ratio R out /R in helps in gain adjustment.
Step 4: Compare the coefficients of (17) with the corresponding coefficients of the denormalized transfer functions from ( 9)-( 16) and determine the values of the remaining passive components for the filter.
Step 5: Choose the nearest values of the passive components from the E24 industrial series for the resistors and the E12 series for the capacitors.The passive components required to implement the proposed PLFs and IPLFs are presented in Tables 5 and 6, respectively.For better accuracy, R out was selected from the E48 series for the IHPPLF.
Table 5.Values of the passive components for the realization of the proposed PLFs (α = 0.5) (Note: All resistances are in kΩ, and all capacitances are in nF).Practical circuit implementations were carried out using the Analog Devices AD844ANtype CFOAs.Supply voltage of 12 volts for the chips was provided from the Agilent E3630A power supply.The magnitude-frequency and phase-frequency measurements were conducted using the OMICRON Lab Bode 100 network analyzer and displayed using the Bode Analyzer Suite software.In this regard, 801 logarithmically spaced frequency points in the range 10 Hz to 100 kHz were considered.The level of the testing harmonic signal was set to V PP = 1 V.

Filter
The receiver bandwidth of the analyzer was fixed at 100 Hz.The time-domain behaviors of the practical filters were observed on Agilent InfiniiVision DSO-X 2002A digital storage oscilloscope.The voltage 1 V PP (default value) was applied to the filter circuit from the Agilent 33521A function/arbitrary waveform generator during measurements of time-domain responses.Figure 9 presents the photograph of the hardware set-up used to experimentally validate the performance of the proposed filters with the display demonstrating the frequency responses for the BPPLF as a test case.

Measurement Results for the PLFs
The experimentally obtained magnitude and phase characteristics for the proposed LPPLF, HPPLF, BPPLF, and BSPLF are compared with the theoretical ones in Figure 10a-d, respectively.All the cases achieve excellent agreement with the ideal responses.The MARE values (determined using L = 801) attained for the practical PLFs exhibiting the LP, HP, BP, and BS behaviors are 0.1006, 0.9668, 0.1307, and 1.0031, respectively.
The time-domain responses of the proposed LPPLF and HPPLF measured at the half-power frequency ( f H ) of 1.11 kHz and 738 Hz, respectively, are shown in Figure 11a,b, respectively.The peak-to-peak output voltage (V OUT(PP) ) obtained for the LPPLF and HPPLF at f H are 700 mV and 680 mV, respectively.Figure 12a-c show the three timedomain input-output waveforms for the BPPLF (namely, BPPLF-a, BPPLF-b, and BPPLF-c) when the input signal is applied at the center frequency ( f 0 = 1.066 kHz), low half-power frequency ( f H,low = 353.9Hz), and high half-power frequency ( f H,high = 3.214 kHz), respectively.
The corresponding values of V OUT(PP) are attained as 1.01 V, 710 mV, and 720 mV, respectively.In Figure 13a,b, the two transient responses for the BSPLF (namely, BSPLF-a, and BSPLF-b) at the input signal frequencies of f H,low (614 Hz) and f H,high (1.47 kHz), respectively, are presented.Both cases yield 710 mV for V OUT(PP) .The values for f 0 , f H , f H,low and f H,high are considered with respect to the experimental measurements.

Measurement Results for the IPLFs
Comparisons between the theoretical and experimental magnitude and phase responses achieved for the proposed ILPPLF, IHPPLF, IBPPLF, and IBSPLF are presented in Figure 15a-d, respectively.The magnitude characteristics for all the cases demonstrate good agreement with the theoretical ones whereas the phase behavior deviates as the operating frequency approaches towards 100 kHz.The practical IPLFs with LP, HP, BP, and BS responses attain the MARE values of 0.2441, 6.8457, 0.2221, and 3.6490, respectively.

Conclusions
In this paper, we presented the optimal and stable rational approximation of the PLFs and their inverse functions.Design examples for the LP, HP, BP, and BS-type PLFs with α = {0.3,0.5, 0.7} and three different objective functions were presented.Statistical performance comparisons and Wilcoxon rank-sum test results revealed that no single objective function attained the best accuracy and solution quality consistency for all types of PLFs.Comparisons with published literature highlighted the improved accuracy for all the proposed PLFs.The drawback of obtaining unstable ILPPLF and IHPPLF for several design orders based on the reported S-K method was circumvented here through the incorporation of design constraints.
Since the PLF/IPLF design is an offline procedure, the inferior computational time of the proposed technique may be waived in favor of attaining the previously-mentioned benefits.CFOAs employed as active elements were used to practically realize the proposed PLFs and IPLFs for all four response types with α = 0.5.The experimental results exemplified an excellent agreement in the magnitude and phase responses with the ideal PLFs.The magnitude behavior for the IPLFs also attained proximity with the theoretical anticipations.For all measurements, the THD remained lower than 0.2%, and the SFDR exceeded 57.23 dBc for the PLFs.In the case of IPLFs, the THD and SFDR values were equal or smaller than 0.21% and larger than 54.72 dBc, respectively.
the Hurwitz determinants[57], and b N = 1.The locations of zeros are also restricted to avoid the right-half s-plane by incorporating the constraints in (

Figure 2 .
Figure 2. MARE comparison using boxplots for the (a) LPPLF, (b) HPPLF, (c) BPPLF, and (d) BSPLF.Note: The most accurate model (lowest MARE min ) for each case is indicated in blue.

Figure 3 .
Figure 3. (a) Magnitude and (b) phase-frequency responses of the proposed LPPLFs (dashed blue) and ILPPLFs (dashed green) as compared to the theoretical ones (solid black).

Figure 4 .
Figure 4. (a) Magnitude and (b) phase-frequency responses of the proposed HPPLFs (dashed blue) and IHPPLFs (dashed green) as compared to the theoretical ones (solid black).

Figure 5 .
Figure 5. (a) Magnitude and (b) phase-frequency responses of the proposed BPPLFs (dashed blue) and IBPPLFs (dashed green) as compared to the theoretical ones (solid black).

Figure 6 .
Figure 6.(a) Magnitude and (b) phase-frequency responses of the proposed BSPLFs (dashed blue) and IBSPLFs (dashed green) as compared to the theoretical ones (solid black).

Figure 7 .
Figure 7. Percentage improvement in absolute relative error of the proposed PLFs with respect to [54].

Figure 9 .
Figure 9. Photograph of the experimental set-up.

Figure
Figure 16a,b show the time-domain input-output waveforms of the practical ILPPLF and IHPPLF measured at the input frequency of f H = 1.212 kHz and 887.8 Hz, respectively.Measurements reveal that V OUT(PP) of 1.259 V and 1.421 V are, respectively, obtained for the ILPPLF and the IHPPLF.The time-domain responses for the IBPPLF (i.e., IBPPLF-a, IBPPLF-b, and IBPPLF-c) at the input frequencies of f 0 = 967 Hz, f H,low = 330.3Hz, and f H,high = 2.919 kHz are presented in Figure 17a-c.The values of V OUT(PP) achieved for these three cases are 0.966 V, 1.374 V, and 1.355 V, respectively.The input-output plots for the IBSPLF (i.e., IBSPLF-a and IBSPLF-b) at the input frequencies of f H,low = 650.3Hz and f H,high = 1.482 kHz are shown in Figure 18a,b, respectively.The measured values of V OUT(PP) at these two frequencies are obtained as 1.438V and 1.426 V, respectively.

Table 1 .
Transfer functions and frequency response expressions of the theoretical PLFs (ω 0 : pole frequency, and Q: quality factor).
b 0 ]; |H D (jω)| and |H P (jω)| denote the magnitude of the theoretical PLF and the proposed approximant, respectively; and the phase angles for the theoretical (∠H D (jω)) and proposed (∠H P (jω)) PLFs are expressed in radians.

Table 2 .
Statistical comparisons of the MARE for designed PLFs.

Table 3 .
Wilcoxon rank-sum test results for the designed PLFs.

Table 4 .
[54]arison of the most accurate proposed PLFs with the published literature[54]in terms of the MARE metric.

Table 6 .
Values of the passive components for the realization of the proposed IPLFs (α = 0.5) (Note: All resistances are in kΩ, and all capacitances are in nF).