Optimized Design of OTA-Based Gyrator Realizing Fractional-Order Inductance Simulator: A Comprehensive Analysis

: A detailed analysis of an operational transconductance ampliﬁer based gyrator implementing a fractional-order inductance simulator is presented. The inﬂuence of active element non-ideal properties on the gyrator operation is investigated and demonstrated by admittance characteristics and formulas for important values and cut-off frequencies in these characteristics. Recommendations to optimize the performance of the gyrator in terms of operation bandwidth, the range of obtainable admittance magnitude, and signal dynamic range are proposed. The theoretical observations are veriﬁed by PSpice simulations of the gyrator with LT1228 integrated circuit.


Introduction
Fractional-order (FO) systems, i.e., systems of non-integer order described by noninteger order differential and integral equations, are able to provide more arbitrary properties and characteristics in comparison with their integer-order counterparts [1,2]. Thus, FO systems are employed in many branches, such as signal processing, modeling various real-world phenomena, or system control [3]. Analog electrical circuits can also be designed with fractional order, and for their implementation, elements with FO impedance are frequently employed. These elements are also often referred to as fractional-order elements (FOEs) [4] or fractors [5]. The impedance of FOE can be expressed in the s-domain by the following non-integer power-law dependence where the real number α is fractional exponent or order of FOE and the positive constant F represents fractance having the unit Ω −1 ·sec α or F·sec α−1 ; whereas, here "F" is farad. After the substitution s = jω, the relation (1) can be written in the magnitude-phase format in the frequency domain: The roll-off of the impedance magnitude is −20α dB per decade of frequency, and the phase remains constant regardless of frequency. If α > 0, the FOE has a capacitive character and is labeled as FO capacitor or capacitive fractor. Its impedance magnitude decreases with increasing frequency, and the phase angle is negative. If α < 0, the FOE the paper [22] can also be mentioned. However, here a lossy FO inductor is implemented by a current feedback operational amplifier (CFOA) and FO capacitor. The analysis of parasitic properties is performed only at the level of symbolic relations, and the influence of individual CFOA parasitics on the properties of the FO inductor is not explained. This paper is organized as follows: In Section 2, the theory of OTA-based gyrator implementing FO inductance simulator is described. The influence of OTA parasitic properties on the gyrator function is analyzed, and dynamic properties are studied to prevent the amplifiers from overloading. Design guidelines and compensation possibilities are given to minimize the effect of the parasitic properties. Section 3 contains verification of the theoretical findings by PSpice simulations. The practical design of FO inductor for an electrical model of the human respiratory system is demonstrated as an example. Finally, our conclusions are given in Section 4.

Ideal OTA-Based Gyrator
Ideal OTA [23] (Figure 1) is a differential voltage-controlled current source with output current given by the following relation where g m is transconductance gain, which is usually electronically adjustable by an external DC current I SET or DC voltage V SET . Input and output terminals of an ideal OTA have infinite internal impedance. The gyrator circuit employing two OTAs is depicted in Figure 2 [23].
Appl. Sci. 2021, 11, x FOR PEER REVIEW 3 of 21 this work provides a FO inductor with electronically adjustable fractance, and it can be easily converted to a floating variant by using one of the OTAs with differential outputs. Among the works dealing with FO inductors, which contain an analysis of active element parasitic properties, the paper [22] can also be mentioned. However, here a lossy FO inductor is implemented by a current feedback operational amplifier (CFOA) and FO capacitor. The analysis of parasitic properties is performed only at the level of symbolic relations, and the influence of individual CFOA parasitics on the properties of the FO inductor is not explained. This paper is organized as follows: In Section 2, the theory of OTA-based gyrator implementing FO inductance simulator is described. The influence of OTA parasitic properties on the gyrator function is analyzed, and dynamic properties are studied to prevent the amplifiers from overloading. Design guidelines and compensation possibilities are given to minimize the effect of the parasitic properties. Section 3 contains verification of the theoretical findings by PSpice simulations. The practical design of FO inductor for an electrical model of the human respiratory system is demonstrated as an example. Finally, our conclusions are given in Section 4.

Ideal OTA-Based Gyrator
Ideal OTA [23] (Figure 1) is a differential voltage-controlled current source with output current given by the following relation where gm is transconductance gain, which is usually electronically adjustable by an external DC current ISET or DC voltage VSET. Input and output terminals of an ideal OTA have infinite internal impedance. The gyrator circuit employing two OTAs is depicted in Figure 2 [23]. this work provides a FO inductor with electronically adjustable fractance, and it can be easily converted to a floating variant by using one of the OTAs with differential outputs. Among the works dealing with FO inductors, which contain an analysis of active element parasitic properties, the paper [22] can also be mentioned. However, here a lossy FO inductor is implemented by a current feedback operational amplifier (CFOA) and FO capacitor. The analysis of parasitic properties is performed only at the level of symbolic relations, and the influence of individual CFOA parasitics on the properties of the FO inductor is not explained. This paper is organized as follows: In Section 2, the theory of OTA-based gyrator implementing FO inductance simulator is described. The influence of OTA parasitic properties on the gyrator function is analyzed, and dynamic properties are studied to prevent the amplifiers from overloading. Design guidelines and compensation possibilities are given to minimize the effect of the parasitic properties. Section 3 contains verification of the theoretical findings by PSpice simulations. The practical design of FO inductor for an electrical model of the human respiratory system is demonstrated as an example. Finally, our conclusions are given in Section 4.

Ideal OTA-Based Gyrator
Ideal OTA [23] (Figure 1) is a differential voltage-controlled current source with output current given by the following relation where gm is transconductance gain, which is usually electronically adjustable by an external DC current ISET or DC voltage VSET. Input and output terminals of an ideal OTA have infinite internal impedance. The gyrator circuit employing two OTAs is depicted in Figure 2 [23]. This gyrator provides inversion of the loading admittance Y as apparent from the following relation for its input admittance. A capacitor is most often connected as the loading element, which yields the inductive character of the input admittance. When connecting a FO capacitor with admittance Y(s) = s α F and 0 < α < 1 as the gyrator's load, the resulting input admittance is which represents FO inductor admittance with the magnitude and phase

Practical OTA-Based Gyrator
Parasitic properties of a real OTA are commonly modeled by shunt resistances and capacitances of the input and output terminals. Input resistance is usually very large, especially in CMOS technology; whereas, input capacitance is typically in the range of tens to the low hundreds of femtofarads when integrated on a substrate of standard CMOS processes. The output resistance of OTA generally reaches several hundreds of kiloohms and output capacitance tens to a few hundreds of femtofarads [24,25]. For some OTA implementations, the parasitic properties may be dependent on the OTA operating conditions, e.g., the terminal conductance may change depending on the set transconductance, which is also the case of the amplifier used in the computer simulation example in this article.

Influence of OTA Terminal Impedances
Considering the aforementioned OTA parasitic properties, the gyrator schematic from Figure 2 can be modified to the form in Figure 3.
This gyrator provides inversion of the loading admittance Y as apparent from the following relation for its input admittance.
A capacitor is most often connected as the loading element, which yields the inductive character of the input admittance. When connecting a FO capacitor with admittance Y(s) = s α F and 0 < α < 1 as the gyrator's load, the resulting input admittance is which represents FO inductor admittance with the magnitude and phase

Practical OTA-Based Gyrator
Parasitic properties of a real OTA are commonly modeled by shunt resistances and capacitances of the input and output terminals. Input resistance is usually very large, especially in CMOS technology; whereas, input capacitance is typically in the range of tens to the low hundreds of femtofarads when integrated on a substrate of standard CMOS processes. The output resistance of OTA generally reaches several hundreds of kiloohms and output capacitance tens to a few hundreds of femtofarads [24,25]. For some OTA implementations, the parasitic properties may be dependent on the OTA operating conditions, e.g., the terminal conductance may change depending on the set transconductance, which is also the case of the amplifier used in the computer simulation example in this article.

Influence of OTA Terminal Impedances
Considering the aforementioned OTA parasitic properties, the gyrator schematic from Figure 2 can be modified to the form in Figure 3. The parasitic properties of OTA1 output and OTA2 input respective to the ground are together modeled with GP1 and CP1 elements; whereas, the parasitics of OTA1 input and OTA2 output are modeled with GP2 and CP2. The relation for input admittance of the gyrator loaded with FO capacitor taking into account the OTA parasitic properties is The parasitic properties of OTA 1 output and OTA 2 input respective to the ground are together modeled with G P1 and C P1 elements; whereas, the parasitics of OTA 1 input and OTA 2 output are modeled with G P2 and C P2 . The relation for input admittance of the gyrator loaded with FO capacitor taking into account the OTA parasitic properties is This relation needs to be analyzed to determine the effect of the parasitic properties. To illustrate the meaning of (7), the passive circuit with the identical input admittance Y IN '(s) can be depicted in Figure 4.
This relation needs to be analyzed to determine the effect of the parasitic properties. To illustrate the meaning of (7), the passive circuit with the identical input admittance YIN'(s) can be depicted in Figure 4.  (7); the relations in the schematic express admittances of the particular elements. Analyzing (7) and Figure 4, the gyrator transforms the FO capacitor to FO inductor with admittance gm1gm2/(s α F), which is in agreement with the relation (5) valid for the ideal case. The parasitic properties CP1 and GP1 are also transformed by the gyrator and appear in the network in Figure 4 as an integer-order inductor with admittance gm1gm2/(sCP1) and resistor with conductance gm1gm2/GP1, respectively. The admittance of the parasitic elements GP2 and CP2 is according to (7) additional part of the input admittance YIN' (s) and therefore, these elements appear at the overall input of the structure in Figure  4.
Considering the above-mentioned typical values of the OTA parasitic properties, the format of the admittance magnitude and phase frequency characteristics resulting from (7) can be displayed as the asymptotic Bode plots in Figure 5. . Bode plots of (a) magnitude and (b) phase admittance characteristics according to (7). Blue parts denote ideal function without the influence of parasitics.
The lower cut-off frequency of the correct gyrator operation is approximately equal to Analyzing (7) and Figure 4, the gyrator transforms the FO capacitor to FO inductor with admittance g m1 g m2 /(s α F), which is in agreement with the relation (5) valid for the ideal case. The parasitic properties C P1 and G P1 are also transformed by the gyrator and appear in the network in Figure 4 as an integer-order inductor with admittance g m1 g m2 /(sC P1 ) and resistor with conductance g m1 g m2 /G P1 , respectively. The admittance of the parasitic elements G P2 and C P2 is according to (7) additional part of the input admittance Y IN ' (s) and therefore, these elements appear at the overall input of the structure in Figure 4.
Considering the above-mentioned typical values of the OTA parasitic properties, the format of the admittance magnitude and phase frequency characteristics resulting from (7) can be displayed as the asymptotic Bode plots in Figure 5.
This relation needs to be analyzed to determine the effect of the parasitic properties. To illustrate the meaning of (7), the passive circuit with the identical input admittance YIN'(s) can be depicted in Figure 4.  (7); the relations in the schematic express admittances of the particular elements. Analyzing (7) and Figure 4, the gyrator transforms the FO capacitor to FO inductor with admittance gm1gm2/(s α F), which is in agreement with the relation (5) valid for the ideal case. The parasitic properties CP1 and GP1 are also transformed by the gyrator and appear in the network in Figure 4 as an integer-order inductor with admittance gm1gm2/(sCP1) and resistor with conductance gm1gm2/GP1, respectively. The admittance of the parasitic elements GP2 and CP2 is according to (7) additional part of the input admittance YIN' (s) and therefore, these elements appear at the overall input of the structure in Figure  4.
Considering the above-mentioned typical values of the OTA parasitic properties, the format of the admittance magnitude and phase frequency characteristics resulting from (7) can be displayed as the asymptotic Bode plots in Figure 5. . Bode plots of (a) magnitude and (b) phase admittance characteristics according to (7). Blue parts denote ideal function without the influence of parasitics.
The lower cut-off frequency of the correct gyrator operation is approximately equal to Figure 5. Bode plots of (a) magnitude and (b) phase admittance characteristics according to (7). Blue parts denote ideal function without the influence of parasitics.
The lower cut-off frequency of the correct gyrator operation is approximately equal to At this frequency, the admittance magnitude of the gyrator load |ω α F| equals to G P1 . Below this frequency, the gyrator exhibits parasitic behavior with limited input admittance magnitude g m1 g m2 /G P1 and zero phase. This value represents the highest achievable admittance magnitude at the gyrator input. The value of G P2 has no effect below the frequency ω GP1 , as we assume G P2 g m1 g m2 /G P1 . Above the frequency ω GP1 , the operation of the gyrator can be considered correct, i.e., providing the admittance of FO inductor given by (5), which is marked by the blue color in Figure 5. The magnitude of the gyrator input admittance decreases with the roll-off rate −20α dB/dec, and the constant phase is −απ/2. Note that the admittance characteristics of an ideal FO inductor can be obtained by linearly extending the blue lines in Figure 5 over the entire frequency range.
Based on the relationships between the values of the elements in Figure 4, the upper cut-off frequency of the correct behavior of the gyrator may be equal to where the input admittance magnitude |Y IN '| decreases to the parasitic conductance G P2 . Above this frequency, the input admittance has the magnitude value G P2 and is purely real. Further, above the frequency the effect of the parasitic capacitance C P2 prevails, and the input admittance Y IN ' has purely capacitive character, which even more differs from the required FO inductor. As demonstrated in Figure 5 by the dashed part of the characteristic, when the conductance G P2 is lower than the higher value of the two admittances g m1 g m2 /(ω α F) and ωC P2 , the upper cut-off frequency of the correct gyrator behavior is In this case, the FO inductance appearing at the input port of the gyrator and the parasitic capacitance C P2 create a parallel resonant circuit, which can be a potential source of instability. Analyzing the relation (7), the real part of Y IN ' cannot be lower than G P2 , which means that the parasitic conductance G P2 ensures damping of the resonance. The value of α also has an impact on the gyrator stability. The more it is less than one, the more the resonance is damped.
The standard inductor in Figure 4 resulting from immittance inversion of the parasitic capacitance C P1 affects the input admittance of the gyrator negligibly. Its effect would be in decreasing the input admittance above a cut-off frequency that is higher than ω CP2 ; however, in this frequency band, the input admittance grows, due to the admittance of C P2 , which outweighs.
The following recommendations can be drawn from the previous analysis to improve the performance of the gyrator mainly as far as operation bandwidth, and range of obtainable admittance magnitude are concerned:

•
If there is a possibility to modify the fractance F, increase the values of both (g m1 g m2 ) and F such that the ratio (g m1 g m2 )/F remains constant. This decreases the lower cut-off frequency ω GP1 , and thus, expands the band of correct operation of the gyrator to lower frequencies without affecting the FO inductance. Moreover, the maximum obtainable magnitude of input admittance g m1 g m2 /G P1 increases. • From the OTA parasitic properties, the main effort should be to decrease G P1 , which also brings a decrease of ω GP1 and an increase of maximum input admittance magnitude. Adherence to the recommendation on decreasing G P1 will also efficiently increase ω GP2 , as due to the gyrator topology, we may assume G P1 ≈ G P2 . This can be achieved by choosing an OTA structure that offers low parasitic conductance at input and output terminals. An alternative way to reduce a parasitic conductance in the gyrator is connecting a negative conductance in parallel as described in Section 2.2.2.
To further optimize the gyrator operation at high frequencies in terms of increasing ω CP2 , it is recommended to choose an OTA structure with low C P2 . It can also be observed that for lower values of α, the bandwidth of the correct gyrator operation is increased, which is related to the lower slope of the admittance magnitude characteristic. Thus, a gyrator implementing a FO inductor provides higher bandwidth than a gyrator implementing a classic inductor.

Compensation of Parasitic Conductance
When a parasitic conductance of a gyrator node is to be decreased as described in the previous section, a simple compensation circuit with negative input conductance, shown in Figure 6, can be used. This circuit is connected with its COMP terminal to the node where the parasitic conductance is to be compensated.
the gyrator is connecting a negative conductance in parallel as described in Section 2.2.2. To further optimize the gyrator operation at high frequencies in terms of increasing ωCP2, it is recommended to choose an OTA structure with low CP2.
It can also be observed that for lower values of α, the bandwidth of the correct gyrator operation is increased, which is related to the lower slope of the admittance magnitude characteristic. Thus, a gyrator implementing a FO inductor provides higher bandwidth than a gyrator implementing a classic inductor.

Compensation of Parasitic Conductance
When a parasitic conductance of a gyrator node is to be decreased as described in the previous section, a simple compensation circuit with negative input conductance, shown in Figure 6, can be used. This circuit is connected with its COMP terminal to the node where the parasitic conductance is to be compensated. The conductance of the COMP terminal relative to ground equals to and can be set to an appropriate negative value by the electronic setting of gmC. Here GP is the total parasitic conductance of OTAC input and output terminals. If the implementation of OTAC is identical with the OTAs of the gyrator, all the parasitic conductances are equal (GP = GP1 = GP2). From (8), it is apparent that when the admittance GP1 is compensated in this way to almost zero value, the cut-off frequency ωGP1 can be theoretically shifted to very low values, and the maximum magnitude of the gyrator input admittance can be increased significantly. However, the demands on the accuracy of the compensation circuit increase, and it is important to prevent reaching a negative conductance of the node. The conductance GP2 is worth to be compensated when it affects the gyrator input admittance at high frequencies between ωGP2 and ωCP2, i.e., when the solid line (and not the dashed line) in Figure 5 holds. Remember that when α tends to unity, the conductance GP2 is an important damping factor of the above-mentioned LC resonance, and thus, it cannot be excessively reduced. It should also be noted that when connecting the circuit in Figure 6 to a gyrator node, the parasitic capacitance of this node increases by the own parasitic capacitance of the compensation circuit, which is equal to CP1 and CP2, assuming that all OTAs are the same. Thus, it is necessary to take this property into account in the relations containing the parasitic capacitance of the node being compensated. Note that the circuit in Figure 6 cannot be used for compensating the gyrator parasitic capacitance.

Dynamic Properties
Next to the non-ideal frequency analysis, it is also important to investigate the dynamic properties of the gyrator such that the amplifiers operate in linear mode and the The conductance of the COMP terminal relative to ground equals to and can be set to an appropriate negative value by the electronic setting of g mC . Here G P is the total parasitic conductance of OTA C input and output terminals. If the implementation of OTA C is identical with the OTAs of the gyrator, all the parasitic conductances are equal (G P = G P1 = G P2 ). From (8), it is apparent that when the admittance G P1 is compensated in this way to almost zero value, the cut-off frequency ω GP1 can be theoretically shifted to very low values, and the maximum magnitude of the gyrator input admittance can be increased significantly. However, the demands on the accuracy of the compensation circuit increase, and it is important to prevent reaching a negative conductance of the node. The conductance G P2 is worth to be compensated when it affects the gyrator input admittance at high frequencies between ω GP2 and ω CP2 , i.e., when the solid line (and not the dashed line) in Figure 5 holds. Remember that when α tends to unity, the conductance G P2 is an important damping factor of the above-mentioned LC resonance, and thus, it cannot be excessively reduced. It should also be noted that when connecting the circuit in Figure 6 to a gyrator node, the parasitic capacitance of this node increases by the own parasitic capacitance of the compensation circuit, which is equal to C P1 and C P2 , assuming that all OTAs are the same. Thus, it is necessary to take this property into account in the relations containing the parasitic capacitance of the node being compensated. Note that the circuit in Figure 6 cannot be used for compensating the gyrator parasitic capacitance.

Dynamic Properties
Next to the non-ideal frequency analysis, it is also important to investigate the dynamic properties of the gyrator such that the amplifiers operate in linear mode and the maximum voltage and current values of their terminals are not exceeded. As far as common OTAs are concerned, the most limiting is the non-linearity of their input transistor stage, which takes effect already at differential input voltages of several tens of millivolts. Let us denote the maximum allowed amplitude of OTA input voltage for linear operation V OTA,MAX . The output drive current of OTA is also limited to a maximum value I OTA,MAX ; whereas, it holds I OTA,MAX = g m V OTA,MAX . The value of V OTA,MAX , and also maximum allowable transconductance g m can be determined from the respective OTA documentation. The following two conditions must be met to avoid the non-linear operation of the gyrator:

•
The gyrator input voltage V IN must not be higher than the allowed input voltage range of OTA 1 Apparently, the allowed V IN according to (14) changes with frequency. It reaches the minimum value at the lowest operating frequency where |Y IN | is maximal. It can be concluded that the gyrator input voltage V IN must meet both the conditions (13) and (14) simultaneously when the linear operation of both OTAs is required. If the conditions for maximum gyrator input current I IN are required, they can be easily derived from (13)

Guidelines for Setting of Gyrator Parameters
Practical guidelines for gyrator parameter design were compiled based on the above observations. By following these rules, in addition to the correct input admittance, we also ensure its highest possible frequency and magnitude range and prevent the active elements from non-linear operation and overloading. The starting point of the design is the required value of input admittance Y IN specified, e.g., by F IN and α, or magnitude at specific frequency and phase, see (5) and (6). The first set of guidelines (labeled A) can be used if there is available any value of the fractance F of the FO capacitor loading the gyrator.
A.1. Choose the same g m1 and g m2 as the highest possible values for the correct functionality of the selected OTAs. The maximum transconductance values can be found in OTA documentation and are usually several hundreds of µS for CMOS and several tens of mS for BJT implementations. The choice of maximal transconductances extends the operating band to low frequencies. If this is not a priority, choose the transconductances equal to the required input admittance magnitude at the commonly used frequency, e.g., at the geometric mean of the operating band. A.2. Calculate the value of fractance F from the magnitude part of (6). If the value of F turns out too high for implementation, then it is possible to use the next set of rules labeled B, where a lower value of F can be set at the beginning. If F turns out too small, the only solution is to choose OTAs with higher maximum transconductances. A.3. Check the frequency range of the correct functionality of the gyrator, according to (8), (9), or (11). If the range is not sufficient, consider reducing parasitic conductance(s) by the compensation described in Section 2.2.2. A.4. To prevent the OTAs from the non-linear operation, ensure that the magnitude of the gyrator input voltage V IN meets both the conditions (13) and (14) simultaneously in the operating frequency band.
The next set of rules (labeled B) is applicable if a specific value of the fractance F should be set at the beginning. B.1. Set a specific value of the fractance F. If a cut-off frequency ω GP1 should be reached without compensation of parasitic conductance G P1 , the selected value F must be at least as determined from (8). B.2. Determine the required product g m1 g m2 for the selected F and given F IN from the magnitude part of (6). Evaluate whether the product is achievable by the selected OTAs. If not, F must be reduced for the design to continue (regardless of fulfilling (8) in the previous step). B.3. Choose g m2 as the highest possible value for the proper functionality of the selected OTA 2 . B.4. Calculate g m1 from the magnitude part of (6). If the ratio g m1 /g m2 turns out excessively low, then it is possible to decrease g m2 and an increase g m1 such that their product remains the same. However, keep in mind that decreasing g m2 may decrease the allowed gyrator input voltage, according to (14).
B.5. Check the frequency range of the correct functionality of the gyrator, according to (8), (9), or (11). If the range is not sufficient, consider reducing parasitic conductance(s) by the compensation described in Section 2.2.2. B.6. To prevent the OTAs from the non-linear operation, ensure that the magnitude of the gyrator input voltage V IN meets both the conditions (13) and (14) simultaneously in the operating frequency band.
To add more clarity, a flowchart of the above guidelines for the setting of gyrator parameters is enclosed in Figure 7.

Differential (Floating) Version of the Gyrator
If the OTA 2 with differential (balanced) output is available, it is possible to easily obtain a differential gyrator with floating input according to Figure 8. The describing relations (4) to (6) remain valid; however, the parasitics G P2 and C P2 must be considered with halved values G P2 /2 and C P2 /2 in the analyses, as the OTA terminal parasitic elements are connected in series here (through ground). These reduced parasitic values improve the functionality of the gyrator at high frequencies. Another advantage of the floating gyrator is that the loading element Y remains grounded.     Ensure that VIN meets both (13) and (14) End of the design

Simulation Results and Discussion
To verify the presented theoretical studies and also demonstrate the design optimization, the computer simulation of the gyrator behavior was carried out using OrCAD PSpice software. The commercially available integrated circuit LT1228 [26] was used as OTA to enable the readers the possibility of reproducing the results. This element provides transconductance g m controlled by an external DC current I SET ; whereas, it holds g m = 10I SET . An analysis of the LT1228 PSpice model parasitic properties revealed that the conductance of input pins G P,IN is variable and depends on g m with the relation G P,IN ≈ 0.005g m . It is, therefore, necessary to take into account that the functionality of the gyrator will be more affected by the parasitic conductance at higher values of g m . The input parasitic conductance is not included in the LT1228 PSpice model, i.e., it is considered zero. The parasitic conductance G P,OUT of output pin is approximately 0.23 nS and can be neglected in comparison with G P,IN for common values of g m . The parasitic output capacitance is 5.2 pF. Thus, the modeled parasitic conductances G P1 and G P2 of the gyrator correspond to 0.005g m2 and 0.005g m1 , respectively, and both the parasitic capacitances C P1 and C P2 are 5.2 pF.

Gyrator Simulation and Optimization
The input admittance magnitude of the gyrator chosen in this subsection is |Y IN | = 10 mS at 1 kHz, and the orders α are 0.25, 0.5, 0.75, 1, i.e., three of them representing FO and one classic inductance simulator. Two sets of element parameters will be used for demonstration. The first set (labeled set 1) will be designed only to provide the required input admittance, regardless of the previously proposed recommendations concerning the wide frequency and magnitude range. Then, the elements of set 2 will be designed regarding these recommendations to demonstrate the optimization of admittance characteristics. The schematic of the gyrator in the OrCAD Capture editor is depicted in Figure 9. Unless otherwise noted, the circuit is considered without the compensation part indicated by the dashed box.

Simulation Results and Discussion
To verify the presented theoretical studies and also demonstrate the design optimization, the computer simulation of the gyrator behavior was carried out using OrCAD PSpice software. The commercially available integrated circuit LT1228 [26] was used as OTA to enable the readers the possibility of reproducing the results. This element provides transconductance gm controlled by an external DC current ISET; whereas, it holds gm = 10ISET. An analysis of the LT1228 PSpice model parasitic properties revealed that the conductance of input pins GP,IN is variable and depends on gm with the relation GP,IN ≈ 0.005gm. It is, therefore, necessary to take into account that the functionality of the gyrator will be more affected by the parasitic conductance at higher values of gm. The input parasitic conductance is not included in the LT1228 PSpice model, i.e., it is considered zero. The parasitic conductance GP,OUT of output pin is approximately 0.23 nS and can be neglected in comparison with GP,IN for common values of gm. The parasitic output capacitance is 5.2 pF. Thus, the modeled parasitic conductances GP1 and GP2 of the gyrator correspond to 0.005gm2 and 0.005gm1, respectively, and both the parasitic capacitances CP1 and CP2 are 5.2 pF.

Gyrator Simulation and Optimization
The input admittance magnitude of the gyrator chosen in this subsection is |YIN| = 10 mS at 1 kHz, and the orders α are 0.25, 0.5, 0.75, 1, i.e., three of them representing FO and one classic inductance simulator. Two sets of element parameters will be used for demonstration. The first set (labeled set 1) will be designed only to provide the required input admittance, regardless of the previously proposed recommendations concerning the wide frequency and magnitude range. Then, the elements of set 2 will be designed regarding these recommendations to demonstrate the optimization of admittance characteristics. The schematic of the gyrator in the OrCAD Capture editor is depicted in Figure 9. Unless otherwise noted, the circuit is considered without the compensation part indicated by the dashed box. To show solely the influence of parasitic properties of the employed OTAs, ideal FO capacitors were used as gyrator load Y. These elements were modeled in OrCAD by GFREQ part representing voltage-controlled current source described by a table of frequency response in magnitude and phase domain [27]. Interpolation is performed between the entries; whereas, magnitude is interpolated logarithmically and phase linearly. By interconnecting the input and output of the GFREQ part as displayed in Figure 9, the resulting element admittance is described by a series of (input frequency in Hz, magnitude in dB, phase in degrees) triplets.
The OTA transconductances are chosen g m1 = g m2 = 1 mS in the first set of circuit parameters. To obtain the required |Y IN | = 10 mS at 1 kHz, the magnitude of loading admittance according to (4) is |Y| = 0.1 mS at 1 kHz. The orders of the loading admittance are α = 0.25, 0.5, 0.75, 1, i.e., three of them FO and one classic capacitor. The summary of this first set of element parameters is given in Table 1, along with the expected lower and upper cut-off frequencies determined by the relations (8), (9), and (11). Table 1. First set of element parameters for simulation of the circuit in Figure 9 and its expected cut-off frequencies. The second set of gyrator element parameters is designed following rules A from Section 2.3, which should provide optimization of both the frequency and dynamic range of the gyrator input admittance. The maximum transconductances of LT1228 OTA are chosen based on the datasheet g m1 = g m2 = 20 mS. To obtain the required |Y IN | = 10 mS at 1 kHz, the magnitude of loading admittance is |Y| = 40 mS at 1 kHz. The summary of this second set of element parameters with the expected cut-off frequency values is given in Table 2. Table 2. The second set of element parameters for simulation of the circuit in Figure 9 and its expected cut-off frequencies. The values of lower cut-off frequency f GP1 for the second set are lower than for the first set, as expected. As seen in the last column of Table 2 (values that are not in parentheses), the upper cut-off frequencies decreased for α = 0.75 and 1 compared with Table 1. This bandwidth degradation is caused by the parasitic conductance G P2 , which raised from 5 µS for the first element set to 100 µS for the second set, due to increasing the transconductance g m1 (it holds G P2 = 0.005g m1 ). The parasitic conductance G P2 can be decreased as described in Section 2.2.2, resulting in an increase of the upper cut-off frequency. The compensation circuit with negative admittance utilized for this purpose is marked with the dashed box in Figure 9. The transconductance g mC of the OTA in the compensation circuit was set to 99.5 µS, which leads to G P2 ≈ 1 µS and that is an even better value than for the first set of elements. It should be taken into account that the compensation circuit causes a doubling of the parasitic capacitance C P2 , therefore, the results after compensation have improved, but still do not reach the values as for the first set, see the frequencies in the last column of Table 2 in parentheses.
The simulated magnitude and phase-frequency characteristics of the gyrator input admittance are depicted in Figures 10-13 for the α values 0.25, 0.5, 0.75, 1, respectively. The solid black lines are valid for the first set of circuit parameters, the dotted green lines for the second set, and the dashed blue lines indicate ideal results. It is worth to mention again that the simulations in this subsection are performed with ideal FO capacitor loading the gyrator to solely illustrate the effect of OTA non-ideal properties. Common FO capacitors emulated by RC structures of reasonable complexity have limited accuracy in a limited frequency band [28], which would unnecessarily degrade the results, and it would not be clear whether the distortion is caused by non-idealities of OTA or FO capacitor.
values from Tables 1 and 2. Considering the maximum phase error ±2 degrees, the useful frequency ranges for the designed FO inductance simulators are listed in Table 3. Furthermore, the simulations confirm that the admittance magnitude at low frequencies is limited to the value gm1gm2/GP1, which is 200 mS for the first set of elements and 4 S for the second one. The constant magnitude section between ωGP2 and ωCP2 is not apparent in the simulations as the value GP2 is too low to take effect. This section would be visible for the set 2 magnitudes with α values 0.75 and 1 if no GP2 compensation were applied. The effect of the additional capacitance of the compensation circuit is demonstrated in Figures 12a and 13a by increasing the magnitude at high frequencies compared to set 1 results. The simulated results also confirm that the bandwidth of correct gyrator operation increases as α decreases. The simulated admittance characteristics in Figures 10-13 are in accordance with the Bode plots in Figure 5, and the important asymptotic magnitude and phase values are confirmed. The limited frequency band of the correct functionality of the gyrator can be observed, and the simulated magnitude cut-off frequencies correspond to the expected values from Tables 1 and 2. Considering the maximum phase error ±2 degrees, the useful frequency ranges for the designed FO inductance simulators are listed in Table 3.
Furthermore, the simulations confirm that the admittance magnitude at low frequencies is limited to the value g m1 g m2 /G P1 , which is 200 mS for the first set of elements and 4 S for the second one. The constant magnitude section between ω GP2 and ω CP2 is not apparent in the simulations as the value G P2 is too low to take effect. This section would be visible for the set 2 magnitudes with α values 0.75 and 1 if no G P2 compensation were applied. The effect of the additional capacitance of the compensation circuit is demonstrated in Figures 12a and 13a by increasing the magnitude at high frequencies compared to set 1 results. The simulated results also confirm that the bandwidth of correct gyrator operation increases as α decreases.  Figure 13. Simulated magnitude (a) and phase (b) characteristics of gyrator input admittance for α = 1.

Inductance Simulator for Fractional-Order Respiratory Model
Fractional-order electrical circuit models are useful in modeling the mechanical impedance of the respiratory system [11,12]. Mechanical impedance of the human respiratory system can be represented by an equivalent model given in Figure 14.

Inductance Simulator for Fractional-Order Respiratory Model
Fractional-order electrical circuit models are useful in modeling the mechanical impedance of the respiratory system [11,12]. Mechanical impedance of the human respiratory system can be represented by an equivalent model given in Figure 14.
The impedance of this system is expressed as The experimentally extracted values of (18) for human respiratory system are L = 0.0193 H·sec −0.5313 , α = 0.4687, R = 0.0147 Ω, C = 0.9921 F·sec −0.2729 , β = 0.7271 and the frequency range of the model validity is 4 Hz to 48 Hz [12]. The FO inductor with impedance s α L included in this model will be implemented by the discussed OTA-based gyrator R s α L 1 s C  Figure 14. FO impedance model of the human respiratory system [11].
The impedance of this system is expressed as The experimentally extracted values of (15) for human respiratory system are L = 0.0193 H·sec −0.5313 , α = 0.4687, R = 0.0147 Ω, C = 0.9921 F·sec −0.2729 , β = 0.7271 and the frequency range of the model validity is 4 Hz to 48 Hz [12]. The FO inductor with impedance s α L included in this model will be implemented by the discussed OTA-based gyrator in this section as an example of its practical utilization. The grounded version of the model will be demonstrated, as the floating version requires OTA with differential output (see Figure 8), which is not the case of the used integrated circuit LT1228. The value of F IN in accordance with (5) and (15) is equal to 1/L = 51.81 H −1 ·sec 0.5313 . This value represents a very high admittance magnitude of 11.43 S at the corner frequency of 4 Hz, and the implementation of the inductor is therefore challenging.
Rules A, in Section 2.3, will be used to design the gyrator element parameters. The transconductances are chosen g m1 = g m2 = 20 mS again as the maximum possible values for LT1228. To obtain the required |Y IN | = 11.43 S at 4 Hz, the magnitude of loading admittance is |Y| = 34.987 µS at 4 Hz, its order is α = 0.4687, and F = 7.72 µF·sec −0.5313 . Thanks to the narrow frequency band (4 Hz to 48 Hz) compared to the simulations in the previous subsection and to implement it by available components, the FO capacitor loading the gyrator was emulated by the Valsa structure [29] containing standard resistors and capacitors, as shown in Figure 15. Another reason to employ this RC emulator is that the GFREQ part cannot be used in transient analysis, which will demonstrate dynamic optimization of the gyrator at the end of this section. The element values of the Valsa RC circuit computed by the method described in Reference [29] for the mentioned values α, F, maximum phase error ±0.6 degrees, and the frequency range 0.8 Hz to 240 Hz (to have sufficient margin) are listed in Table 4. Figure 14. FO impedance model of the human respiratory system [11].
The impedance of this system is expressed as The experimentally extracted values of (18) for human respiratory system are L = 0.0193 H·sec −0.5313 , α = 0.4687, R = 0.0147 Ω, C = 0.9921 F·sec −0.2729 , β = 0.7271 and the frequency range of the model validity is 4 Hz to 48 Hz [12]. The FO inductor with impedance s α L included in this model will be implemented by the discussed OTA-based gyrator in this section as an example of its practical utilization. The grounded version of the model will be demonstrated, as the floating version requires OTA with differential output (see Figure 8), which is not the case of the used integrated circuit LT1228. The value of FIN in accordance with (5) and (18) is equal to 1/L = 51.81 H −1 ·sec 0.5313 . This value represents a very high admittance magnitude of 11.43 S at the corner frequency of 4 Hz, and the implementation of the inductor is therefore challenging.
Rules A, in Section 2.3, will be used to design the gyrator element parameters. The transconductances are chosen gm1 = gm2 = 20 mS again as the maximum possible values for LT1228. To obtain the required |YIN| = 11.43 S at 4 Hz, the magnitude of loading admittance is |Y| = 34.987 μS at 4 Hz, its order is α = 0.4687, and F = 7.72 μF·sec −0.5313 . Thanks to the narrow frequency band (4 Hz to 48 Hz) compared to the simulations in the previous subsection and to implement it by available components, the FO capacitor loading the gyrator was emulated by the Valsa structure [29] containing standard resistors and capacitors, as shown in Figure 15. Another reason to employ this RC emulator is that the GFREQ part cannot be used in transient analysis, which will demonstrate dynamic optimization of the gyrator at the end of this section. The element values of the Valsa RC circuit computed by the method described in Reference [29] for the mentioned values α, F, maximum phase error ±0.6 degrees, and the frequency range 0.8 Hz to 240 Hz (to have sufficient margin) are listed in Table 4.  The admittance magnitude and phase characteristics of the emulated FO capacitor (solid black lines) and the ideal values (dashed blue lines) are presented in Figure 16. It can be observed that the circuit operates according to the above-mentioned required parameters.  The admittance magnitude and phase characteristics of the emulated FO capacitor (solid black lines) and the ideal values (dashed blue lines) are presented in Figure 16. It can be observed that the circuit operates according to the above-mentioned required parameters.
The gyrator lower cut-off frequency computed by (8) is f GP1 = 37.6 Hz, which is not sufficiently low, as the FO inductor should operate from 4 Hz. This limitation is caused by the parasitic conductance G P1 = 100 µS, which will be decreased by the circuit with negative admittance. Thus, now the compensation circuit in the dashed box in Figure 9 is replaced and connected in parallel with the FO capacitor loading the gyrator. The transconductance g mC is 97 µS, and it leads to G P1 ≈ 3.5 µS. The expected cut-off frequency after compensation is f GP1 ≈ 0.03 Hz and maximum obtainable input admittance magnitude g m1 g m2 /G P1 ≈ 115 S, which is ten times higher than the required admittance magnitude of the implemented FO inductor at 4 Hz. Note that these values are valid for the gyrator loaded with an ideal FO capacitor. It can be expected that the results will be degraded after connecting the emulator of FO capacitor, which has limited accuracy of its admittance. The gyrator lower cut-off frequency computed by (8) is fGP1 = 37.6 Hz, which is not sufficiently low, as the FO inductor should operate from 4 Hz. This limitation is caused by the parasitic conductance GP1 = 100 μS, which will be decreased by the circuit with negative admittance. Thus, now the compensation circuit in the dashed box in Figure 9 is replaced and connected in parallel with the FO capacitor loading the gyrator. The transconductance gmC is 97 μS, and it leads to GP1 ≈ 3.5 μS. The expected cut-off frequency after compensation is fGP1 ≈ 0.03 Hz and maximum obtainable input admittance magnitude gm1gm2/GP1 ≈ 115 S, which is ten times higher than the required admittance magnitude of the implemented FO inductor at 4 Hz. Note that these values are valid for the gyrator loaded with an ideal FO capacitor. It can be expected that the results will be degraded after connecting the emulator of FO capacitor, which has limited accuracy of its admittance.
The simulated magnitude and phase frequency characteristics of the gyrator input admittance are shown in Figure 17. The solid black lines represent simulated results, and the blue lines indicate ideal values with the solid part denoting the frequency band of interest, i.e., 4 Hz to 48 Hz. The simulated magnitude and phase frequency characteristics of the gyrator input admittance are shown in Figure 17. The solid black lines represent simulated results, and the blue lines indicate ideal values with the solid part denoting the frequency band of interest, i.e., 4 Hz to 48 Hz.
In this frequency band, the maximum absolute magnitude error is 1.5 S, and phase error is 3.1 degrees when comparing the simulated results with the ideal ones. It is apparent that the error increases at the lower bound of the frequency band. The magnitude error can be easily corrected by fine increasing the product g m1 g m2 , which can be set electronically. The phase deviation can be reduced by fine-tuning the transconductance g mC ; however, the demands on the accuracy of setting g mC increase, and it is important to prevent overcompensating G P1 to a negative value.
The dynamic optimization was also performed based on the conditions introduced in Section 2.2.3. When following these conditions, only allowed voltages and currents will appear in the gyrator structure, and the amplifiers will operate in linear range without overloading. The maximum input voltage amplitude of both OTAs for their linear operation is considered V OTA,MAX = 50 mV. The maximum input voltage of the gyrator can be determined as the maximum value of V IN that meets both conditions (13) and (14) simultaneously. The maximum input current can then be computed using The following Table 5 shows the maximum allowed amplitudes of voltages and currents at the input (V IN , I IN ) of the gyrator, ensuring linear operation of both OTAs for three frequency points, in the beginning, inside (at 25 Hz), and at the end of the operating band. The computed amplitudes of voltages and currents at the gyrator load (V LOAD , I LOAD ) are also stated. In this frequency band, the maximum absolute magnitude error is 1.5 S, and phase error is 3.1 degrees when comparing the simulated results with the ideal ones. It is apparent that the error increases at the lower bound of the frequency band. The magnitude error can be easily corrected by fine increasing the product gm1gm2, which can be set electronically. The phase deviation can be reduced by fine-tuning the transconductance gmC; however, the demands on the accuracy of setting gmC increase, and it is important to prevent overcompensating GP1 to a negative value.
The dynamic optimization was also performed based on the conditions introduced in Section 2.2.3. When following these conditions, only allowed voltages and currents will appear in the gyrator structure, and the amplifiers will operate in linear range without overloading. The maximum input voltage amplitude of both OTAs for their linear operation is considered VOTA,MAX = 50 mV. The maximum input voltage of the gyrator can be determined as the maximum value of VIN that meets both conditions (13) and (14) simultaneously. The maximum input current can then be computed using IIN = VIN|YIN|. The following Table 5 shows the maximum allowed amplitudes of voltages and currents at the input (VIN, IIN) of the gyrator, ensuring linear operation of both OTAs for three frequency points, in the beginning, inside (at 25 Hz), and at the end of the operating band. The computed amplitudes of voltages and currents at the gyrator load (VLOAD, ILOAD) are also stated.  Applying the rules for dynamic optimization, the transient simulation with the maximum allowed input current amplitude I IN = 1 mA and frequency 25 Hz was carried out, and the resulting current and voltage waveforms are presented in Figure 18. The voltages are relative to ground, and currents are considered positive when flowing into the gyrator (I IN ) and into the loading FO capacitor (I LOAD ).
The delay between input current and voltage found by the simulation is 4.53 msec; the corresponding α is 0.453 which is close to the required value of 0.4687. The amplitude of the simulated input voltage 224 µV is slightly higher than the value 206 µV according to Table 5. This matches the fact that the gyrator input admittance magnitude is lower than the target value, which is also observed in Figure 17a. The current through the gyrator load reaches 4.42 µA, also slightly higher than the expected 4.13 µA. The simulated load voltage of 50.8 mV is almost equal to the value of 50 mV in Table 5. It is verified that the voltage and current amplitudes of the gyrator are optimized and do not cause the non-linear function of active elements.
It should be noted here that this particular FO inductor for the respiratory model is characterized by uncommonly high admittance magnitude at low frequency. Therefore, when it is excited by a current of 1 mA, a very small voltage of less than 300 µV appears at its terminals. Moreover, the gyrator load current is very small, of the order of microamperes. Thus, in the case of real implementation and measurement, these small electrical signals would probably be disturbed by noise. In this case, it would be more appropriate to use OTAs with higher transconductance and maximum output current or at least with a higher maximum voltage and current at the same transconductance. In general, it is suitable for the gyrator to select OTAs with transconductance values comparable to the required input admittance magnitude at the operating frequency. This will provide comparable voltage amplitudes at the input and load of the gyrator, as well as comparable currents at the same points.  Applying the rules for dynamic optimization, the transient simulation with the maximum allowed input current amplitude IIN = 1 mA and frequency 25 Hz was carried out, and the resulting current and voltage waveforms are presented in Figure 18. The voltages are relative to ground, and currents are considered positive when flowing into the gyrator (IIN) and into the loading FO capacitor (ILOAD). The delay between input current and voltage found by the simulation is 4.53 msec; the corresponding α is 0.453 which is close to the required value of 0.4687. The amplitude of the simulated input voltage 224 μV is slightly higher than the value 206 μV according to Table 5. This matches the fact that the gyrator input admittance magnitude is lower than the target value, which is also observed in Figure 17a. The current through the gyrator load reaches 4.42 μA, also slightly higher than the expected 4.13 μA. The simulated load voltage of 50.8 mV is almost equal to the value of 50 mV in Table 5. It is verified that

Conclusions
In this article, we presented a comprehensive analysis of the OTA-based gyrator implementing FO inductor to identify and evaluate the limitations caused by the active element real properties, particularly affecting the operational frequency bandwidth, range of the obtainable input admittance, and signal dynamic range. Based on this analysis, design recommendations were provided to overcome or minimize these limitations caused by the real properties of the assumed active elements. During the gyrator design, it is important to set sufficiently large OTA transconductances that are preferably comparable with the required gyrator input admittance magnitude. This also leads to sufficiently large fractance F, which reduces the circuit sensitivity to the parasitics G P1 and C P1 . To further improve the properties of the gyrator, it is possible to suppress the parasitic conductances by connecting a shunt active compensation circuit with negative conductance also designed with OTA. It is important to take into account that the FO capacitor connected to the gyrator also has limited accuracy and frequency bandwidth. Therefore, it is advisable to match the operating frequency bands of the FO capacitor and the gyrator to avoid unnecessary degradation of the FO inductor function. The proposed optimization steps were verified by OrCAD PSpice simulations while designing FO inductors with different orders; whereas, the FO inductor 0.0193 H·sec −0.5313 with order 0.4687 being part of fractional-order respiratory model operating in the frequency range 4 Hz to 48 Hz was also successfully emulated.