Next Article in Journal
Idleness-Aware Dynamic Power Mode Selection on the i.MX 7ULP IoT Edge Processor
Next Article in Special Issue
Surrogate Assisted Optimization for Low-Voltage Low-Power Circuit Design
Previous Article in Journal
Evidence of Limitations of the Transconductance-to-Drain-Current Method (gm/Id) for Transistor Sizing in 28 nm UTBB FD-SOI Transistors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Implementation of a Fractional-Order Electronically Reconfigurable Lung Impedance Emulator of the Human Respiratory Tree

1
Department of Physics, Electronics Laboratory, University of Patras, GR-26504 Rio Patras, Greece
2
Nanoelectronics Integrated Systems Center (NISC), Nile University, Giza 12588, Egypt
3
Department of Electrical & Computer Engineering, University of Calgary, Calgary, AB T2N 1N4, Canada
*
Author to whom correspondence should be addressed.
Current address: Department of Electrical & Computer Engineering, University of Sharjah, Sharjah 27272, UAE.
J. Low Power Electron. Appl. 2020, 10(2), 18; https://0-doi-org.brum.beds.ac.uk/10.3390/jlpea10020018
Submission received: 28 April 2020 / Revised: 11 May 2020 / Accepted: 13 May 2020 / Published: 16 May 2020
(This article belongs to the Special Issue Low-Power CMOS Analog and Digital Circuits and Filters)

Abstract

:
The fractional-order lung impedance model of the human respiratory tree is implemented in this paper, using Operational Transconductance Amplifiers. The employment of such active element offers electronic adjustment of the impedance characteristics in terms of both elements values and orders. As the MOS transistors in OTAs are biased in the weak inversion region, the power dissipation and the dc bias voltage of operation are also minimized. In addition, the partial fraction expansion tool has been utilized, in order to achieve reduction of the spread of the required time-constants and scaling factors. The performance of the proposed scheme has been evaluated, at post-layout level, using MOS transistors models provided by the 0.35 μ m Austria Mikro Systeme technology CMOS process, and the Cadence IC design suite.

1. Introduction

Fractional calculus is a branch of mathematical analysis which forms a more general description of the conventional integer-order calculus. As it is well known, integer-order calculus deals with differentiation that, given a function f x , is described by the expression d α f x d x α , where α is an integer number ( α Z ) . In fractional calculus, the same expression describes the differentiation, yet the order is now a real ( α R ) number, leading to a wider range of values. The idea of fractional derivative dates back to the end of the 17-th century; it is not until the past few decades, however, when fractional calculus started being used for a variety of scientific problems and applications. The cause of its increasing application is that, as mentioned above, it forms a generalization of the traditional calculus, and thus describes more real problems with higher accuracy [1,2,3,4]. One of the most popular methods, that makes the analysis of a real system easier, is the Laplace transform. This is an efficient way to describe and solve equations in frequency domain with simple calculations and without transformations to time domain during the process. The Laplace transform of the derivative of a function, f ( t ) , in the fractional domain is: L 0 D t α f ( t ) = s α F ( s ) for zero initial conditions [1]. Employing this definition, the impedance in the fractional-domain has the general form:
Z ( s ) = k s α ,
where the variable a represents the order of the element. For 0 < α < 1 the element is called fractional-order capacitor, or sometimes Constant Phase Element (CPE), with k C α being the pseudo-capacitance in Farad/sec 1− α . In the case that 1 < α < 0 the expression in (1) represents a fractional-order inductor with k L α being the pseudo-inductance in Henry/sec 1− α . When α = 0 , 1 , 1 , the element is a resistor, capacitor, and inductor, respectively.
Fractional-order modeling considers the effects of the long-term memory, and due to the offered extra degrees of freedom, it is suitable for modeling, analyzing, and synthesizing electrical, chemical, and biological systems [5,6,7,8,9,10,11,12]. In particular, circuit emulators of individual parts of the human body and even the complete human body play an important role for fast and accurate diagnosis of different diseases and morphological changes arising from environmental impacts such as long time exposure to pollution and contamination. For example, the complete human body impedance model studied in [13,14] offers circuits more accurate than the international IEC 60990 safety standard. It is impossible to test the safety of any designed instrument on real humans and therefore, tests are conducted on circuit models. In a similar manner, circuit models for specific human organs have been widely proposed and used to characterize changes that occur in these organs as a result of diseases, pollution or other types of damages. In this regards, the clear relationship between fractional-order models and fractal geometry has already been highlighted in different articles (see for example [15]). Fractal geometry is especially visible in the tissue structure of the lungs [16]. However, the large variations in the circuit model parameters from one human to the other require developing electronically-adjustable circuits that can be re-configured easily to reflect these changes and allow future integration in smart point of care bio-medical devices (see [17]).
Lungs are located in the chest on either side of the heart in the rib cage. Their purpose is to provide oxygen to the blood [18]. Their shape is conical with a narrow rounded apex at the top, and a broad concave base that rests on the convex surface of the diaphragm. Due to the presence of the heart, the left lung is slightly smaller than the right. Lungs are part of the lower respiratory track beginning from trachea and followed by bronchi and bronchioles. The bronchial tree is a series of passages that supplies air to the alveoli of the lungs. It begins with the trachea, which divides into a left and right bronchus. Broncholies are divided into alveolar ducts that create the alveolar sacs, which contains the alveoli, where gas exchanging happens [19]. Oxygen breathed in and carbon dioxide breathed out. Due to its airway’s branches self similarities (fractal geometry), lungs are a great candidate for the application of fractional-order calculus theory [20,21,22,23,24,25,26]. At a low frequency range inheritance is insignificant, so each airway can be modeled by a gamma cell resistance and capacitance. In biology, when a circuit shares the same impedance with a tissue in a specific frequency range, then it is called equivalent circuit. Although the model that is used to characterize viscoelastic behavior is composed of two resistors and one capacitor, in the case of lungs the electrical equivalent impedance of each airway is considered to be consisted of one resistor and one capacitor [27]. Theoretically, the air flow will be the electric current and the air pressure will be the voltage, respectively.
Regarding to the above description the electrical equivalent, shown in Figure 1, consists of a series connection of a passive resistor and a fractional-order capacitor. The equivalent impedance is given by (2)
Z e q ( s ) = R + 1 C α s α ,
with R being an ohmic resistance, C α being the pseudo-capacitance and 0 < α < 2 being the order of the element [16].
The (scaled) values of the resistances and pseudo-capacitances, as well as of the order of the fractional-order capacitors, taken from seven subjects in the range [0.9 Hz, 6 Hz] are summarized in Table 1.
Table 1 of the manuscript shows that the measured circuit parameters for 7 individuals have a variation in pseudo-capacitance of more than 14 times (from 17.5 nFarad/s0.05 to 238 nFarad/s0.42). Also, the variation in the dispersion coefficient from 0.6 to nearly double the value at 1.1 indicates the importance of re-configurability of any proposed circuit model in order to reflect the measured values of any particular individual, which is the focus of this work.
Owing to the fact that fractional-order capacitors are not available in the market [28], their behavior is approximated through the employment of RC networks [29], or the employment of the gyration techniques implemented using active elements [10]. The first solution is easily realizable, but it does not offer the capability of programming because the values of all passive elements must be changed, in order to change the characteristics of the element. On the other hand, the second method offers electronic adjustment of the characteristics of the element, but this is achieved at the expense of circuit complexity and active component count. The contribution, made in this work, is that a novel solution for implementing emulators of fractional-order capacitors with order in the range (0, 2), which offers significant reduction of the spread of time-constants and scaling factors with regards to the straightforward implementation [30], is presented. This also leads to a significant reduction of the total capacitance area and the total power dissipation. All the above have been achieved keeping: (a) the fully electronic adjustment of the characteristics of the capacitor, including its value as well as its order, and (b) the capability of emulating fractional-order capacitors with order greater than one by the same core.
The paper is organized as follows: the procedure for emulating the behavior of the fractional-order capacitors is presented in Section 2, while the circuit implementation of the proposed scheme is given in Section 3. The performance of the intermediate stages and of the whole system are evaluated at post-layout level in Section 4, using the Cadence IC design suite and MOS transistor models provided by the Austria Mikro Systeme (AMS) 0.35 μ m CMOS process.

2. Emulation of Fractional-Order Capacitors

According to [30], the emulation of a fractional-order capacitor is preformed through the utilization of the Functional Block Diagram (FBD) depicted in Figure 2. The equivalent impedance is given by the expression in (3)
Z C P E ( s ) = 1 g m , V I τ · s q ,
with g m , V I being the transconductance parameter of the voltage-to-current ( V / I ) converter, τ being a time-constant associated to the unity-gain frequency ( ω 0 ) of the differentiator according to the formula: ω 0 = 1 / τ , and 0 < q < 1 being the order of the differentiator. The value of the emulated pseudo-capacitance has the form of (4)
C α = g m , V I · τ q .
Due to the restriction on the values range of the variable q, only the cases #1 to #6 (i.e., α = q < 1 ) are implementable. The last case in Table 1, where the order of the fractional-order capacitor is greater than one (i.e., α = 1.134 ), cannot be emulated by the present form of the FBD in Figure 2. In order to overcome this obstacle, the modified FBD in Figure 3 will be utilized, using an extra integer-order differentiator stage between the fractional-order differentiator and the V / I converter and, thus, α = q + 1 [31].
As the equivalent in Figure 1 is valid over the frequency range f = [0.9 Hz, 6 Hz], an appropriate selection of unity-gain frequency of the differentiator ( f 0 ) could be 2.5 Hz. The values of the pseudo-capacitance and the order α of the fractional-order elements, along with the values of the order q of the fractional-order differentiator and the transconductance g m , V I of the FBD in Figure 3, for the seven cases of emulation are summarized in Table 2.
Taking also into account the results provided in [32], the 2nd-order Continued Fraction Expansion (CFE) approximation will be utilized for approximating the behavior of the fractional-order differentiator with an error less than 5% in the frequency range of interest. The resulted rational transfer function has the form
τ · s q A 2 s 2 + A 1 s + A 0 s 2 + B 1 s + B 0 ,
where in (5) the coefficients A i ( i = 0 , 1 , 2 ) and B j ( j = 1 , 2 ) are positive real numbers provided in [33].
The implementation of (5) can be performed using the multi-feedback Inverse-Follow-the-Leader- Feedback (IFLF) structure demonstrated in Figure 4, where the realized transfer function is
H ( s ) = K 2 s 2 + ( K 1 τ 1 ) s + K 0 τ 1 · τ 2 s 2 + ( 1 τ 1 ) s + 1 τ 1 · τ 2 .
The values of scaling factors K i ( i = 0 , 1 , 2 ) and time-constants τ j ( j = 1 , 2 ) are calculated by equating the coefficients of the numerator and denominator of (5) and (6) [10]. The resulted values are summarized in Table 3. Considering the spread (i.e., the ratio of the maximum and minimum value) of both scaling factors and time-constants, it is derived that the spread of scaling factors for all possible cases is equal to 12, 174, while for the time-constants the corresponding value is 388. Taking into account that in the case of employing OTAs these scaling factors will be controlled by dc currents, it is readily obtained that an extremely wide range of bias currents is required towards this goal.
Another possible solution for overcoming this practical problem is the utilization of the Partial Fraction Expansion (PFE) method introduced in [34]. According to this, the transfer function in (5) can be decomposed in the form of (7)
τ · s q = K 0 + r 1 s p 1 + r 2 s p 2 + + r n s p n ,
where r j ( i = 1 , 2 , , n ) are the residues of (5), p j ( i = 1 , 2 , , n ) are the poles of (5) and K 0 = A 2 .
The expression in (7) can be alternatively written as
H ( s ) = K 0 + j = 1 2 K j 1 1 + τ j s ,
where K j = r j p j and τ i = 1 p j . The implementation of (8) can be performed by the functional block diagram given in Figure 5. Using (5), (7) and (8), the resulted values of times constants are summarized in Table 4. The spread of scaling factors is now significantly reduced into the value 884 and the spread of time-constants is now equal to 362, more or less the same as in the previous case. Therefore, the PFE approximation tool is preferable than the CFE tool with regards to reduction of the spread of scaling factors.

3. Circuit Implementation

In order to approximate the fractional-order differentiator stage, which is required in the functional block diagram in Figure 3, the OTA-C topology depicted in Figure 6 will be utilized. This topology implements the transfer function in (8); the values of scaling factors, provided in Table 4, are performed through an appropriate scaling of the corresponding transconductances and, also, their signs are implemented through an appropriate configuration of the input terminals of the corresponding OTAs. The realized time-constants are given by the general expression in (9)
τ j = C j g m j j = 1 , 2 ,
with g m j and C j being the transconductance of the OTAs and the capacitance associated to the jth- integration stage.
The integer-order differentiator is implemented by the topology demonstrated in Figure 7 [35], and the realized transfer function is given by (10)
H d ( s ) = τ · s = C d g m , d · s .
The V / I conversion stage is implemented by the multiple-output OTA configuration depicted in Figure 8, with the conversion gain being equal to 1 / g m , V I .
In order to implement the whole model, which is provided in Figure 1, a floating resistor emulator is also required. This is realized by the topology, given in Figure 9, where the realized resistance is R = 1 / g m , R .
A possible OTA topology, in order to implement the topologies in Figure 6, Figure 7, Figure 8 and Figure 9, is that presented in Figure 10, which is capable of handling input voltages with maximum amplitude 50 mV, with a total harmonic distortion less than 1% [36]. Considering MOS transistors biased in the sub-threshold region, the realized transconductance is
g m = 5 9 · I B n V T ,
with I B being the bias current, 1 < n < 2 being the sub-threshold slope factor, and V T (≃26 mV at 27 C) being the thermal voltage. Using (9) and (11), the resulted expression of the realized time-constants is that given by (12)
τ j = 9 5 · n C j V T I B j j = 1 , 2 .
According to (11) and (12), the temperature variations affect both the implemented values of transconductance and time-constant due to the dependence of thermal voltage from the temperature. This can be compensated by employing proportional to the absolute temperature (PTAT) dc current sources [37].

4. Simulation Results

Using the expression in (12) and considering that C 1 = 5 pF, C 2 = 70 pF, then the calculated values of the dc bias currents, required for realizing the time-constants in Table 4, are provided in Table 5. The scaling factors { K 0 , K 1 , K 2 } were implemented by multiplying the dc bias current I B by these factors, as it is derived from (11). Using (10) and (12), and considering that C d = 25 pF, the dc bias current of the OTAs employed for implementing an integer-order differentiator with time-constant 63.69 ms (i.e., f 0 = 2.5 Hz) was equal to 22.1 pA. In addition, the values of the dc bias currents for implementing the V / I converter stage and the floating resistors, calculated using (4), (11) and Table 1, are given in Table 5. The evaluation of the behavior of the proposed lung impedance model will be performed using the Cadence IC design suite and MOS transistor models provided by the Austria Mikro Systeme (AMS) 0.35 μ m CMOS process. Considering that V D D = V S S = 0.75 V, the MOS transistors aspect ratios of the OTAs used in Figure 6, Figure 7 and Figure 8, selected in such way that all transistors will operate in the sub-threshold region, are provided in Table 6.
The layout design of the proposed lung impedance emulator is demonstrated in Figure 11, where the silicon area is 555.35 μ m × 485.75 μ m. As a first-step the correct operation of the building blocks of the model in Figure 1 was verified. The obtained impedance magnitude and phase responses of the CPE emulator are given in Figure 12. The characteristic values of the gain and phase of the fractional-order differentiator at the center frequency of the approximation (i.e., f 0 = 2.5 Hz) as well as the values of the implemented resistances of the employed V / I converter stage in Figure 8 and resistor emulator in Figure 9 are summarized in Table 7, with the corresponding theoretically predicted values (i.e., derived using ideal circuit elements) given in parentheses. The gain and phase of the integer-order differentiator in Figure 7 at frequency f 0 = 2.5 Hz, which is used in the case #7 for implementing the order α = 1 + q = 1.134 , were 1 and +88.3° with the theoretical values being 1 and +90°, respectively.
The second step is the consideration of the total circuit of the lung impedance emulator in Figure 11. The obtained impedance magnitude and phase responses are provided in Figure 13, where the corresponding theoretically predicted plots are given by dashes. The most important performance characteristics are given in Table 8, along with the theoretical values given between parentheses. Considering the case R = 6.44 M Ω , C a = 238 nF/s 1 α and α = 0.5808, the derived statistical plots about the magnitude and phase of the impedance at the center frequency of the approximation (i.e., f 0 = 2.5 Hz) are provided in Figure 14.
These results are derived, by employing the Monte-Carlo analysis offered by the Analog Design Environment of the Cadence IC design suite, for a number of N = 500 runs and considering variations of the channel width and length, gate oxide thickness, threshold voltage, carrier mobility, and substrate and channel doping concentration. The values of standard deviation were 0.09 M Ω and 0.25 , with the nominal values being 6.89 M Ω and 5.62 , confirming that the proposed structure has reasonable sensitivity characteristics.
It must be mentioned at this point that, in the simulations, the tuning operation is performed by altering the bias current of the transconductors. In the case of a practical implementation with large number of circuit elements, this can be performed by using an appropriate digital circuitry as in [38,39,40] or using an auxiliary analog system as that proposed in [41].

5. Conclusions

The proposed OTA-C emulator of the lung impedance in the human respiratory tree offers design flexibility and versatility, in the sense that the parameters are fully controllable by appropriate dc bias currents. Thus, the same topology is capable of implementing the results derived from subjects without any restriction about their condition, as well as about the experimental setup. The main contribution with regards to the literature is that the spread of dc bias currents is minimized, resulting in a significant reduction of the dc power consumption. Future research work includes the development of simpler structures where all the aforementioned benefits will be preserved, as well as the employment of the presented concept, which is general, in the implementation of various types of electrical equivalents in biomedical applications. It is important to note that developing sensors for different lung diseases has become an increasingly important topic due to its huge economic impact [42]. In this regards, bio-impedance-based sensors are a competitive solution and the proposed emulators could be an attractive candidate for the implementation of such systems.

Author Contributions

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

Funding

This research received no external funding.

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). This article is based upon work from COST Action CA15225, a network supported by COST (European Cooperation in Science and Technology).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
AMSAustria Mikro Systeme
CFEContinued Fraction Expansion
CMOSComplementary Metal-Oxide-Semiconductor
CPEConstant Phase Element
FBDFunctional Block Diagram
ICIntegrated Circuits
MOSMetal-Oxide-Semiconductor
OTAOperational Transconductance Amplifier
PFEPartial Fraction Expansion

References

  1. Podlubny, I. Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications; Elsevier: Amsterdam, The Netherlands, 1998. [Google Scholar]
  2. Monje, C.A.; Chen, Y.; Vinagre, B.M.; Xue, D.; Feliu-Batlle, V. Fractional-Order Systems and Controls: Fundamentals and Applications; Springer Science & Business Media: Berlin, Germany, 2010. [Google Scholar]
  3. Elwakil, A.S. Fractional-order circuits and systems: An emerging interdisciplinary research area. IEEE Circuits Syst. Mag. 2010, 10, 40–50. [Google Scholar] [CrossRef]
  4. Das, S.; Pan, I. Fractional oRder Signal Processing: Introductory Concepts and Applications; Springer Science & Business Media: Berlin, Germany, 2011. [Google Scholar]
  5. Freeborn, T.J. A survey of fractional-order circuit models for biology and biomedicine. IEEE J. Emerg. Sel. Top. Crcuits Syst. 2013, 3, 416–424. [Google Scholar] [CrossRef]
  6. Gómez-Aguilar, J.; Yépez-Martínez, H.; Escobar-Jiménez, R.; Astorga-Zaragoza, C.; Reyes-Reyes, J. Analytical and numerical solutions of electrical circuits described by fractional derivatives. Appl. Math. Model. 2016, 40, 9079–9094. [Google Scholar] [CrossRef]
  7. Sun, H.; Zhang, Y.; Baleanu, D.; Chen, W.; Chen, Y. A new collection of real world applications of fractional calculus in science and engineering. Commun. Nonlinear Sci. Numer. Simul. 2018, 64, 213–231. [Google Scholar] [CrossRef]
  8. Psychalinos, C.; Elwakil, A.S.; Radwan, A.G.; Biswas, K. Guest editorial: Fractional-order circuits and systems: Theory, design, and applications. Circuits Syst. Signal Process. 2016, 35, 1807–1813. [Google Scholar] [CrossRef] [Green Version]
  9. Psychalinos, C.; Elwakil, A.S.; Allagui, A.; Tepljakov, A. Special Issue on recent advances in the design and applications of fractional-order circuits and systems. AEU Int. J. Electron. Commun. 2017, 100, 132–135. [Google Scholar] [CrossRef]
  10. Tsirimokou, G.; Psychalinos, C.; Elwakil, A. Design of CMOS Analog Integrated Fractional-Order Circuits: Applications in Medicine and Biology; Springer: Berlin, Germany, 2017. [Google Scholar] [CrossRef]
  11. Ionescu, C.; Lopes, A.; Copot, D.; Machado, J.T.; Bates, J. The role of fractional calculus in modeling biological phenomena: A review. Commun. Nonlinear Sci. Numer. Simul. 2017, 51, 141–159. [Google Scholar] [CrossRef]
  12. Psychalinos, C.; Elwakil, A.S.; Maundy, B.; Radwan, A.G. Editorial note: Special Issue on the Design and implementation of fractional-order circuits and systems in real-world applications. Microelectron. J. 2019, 90, 154–155. [Google Scholar] [CrossRef]
  13. De Santis, V.; Beeckman, P.A.; Lampasi, D.A.; Feliziani, M. Assessment of human body impedance for safety requirements against contact currents for frequencies up to 110 MHz. IEEE Trans. Biomed. Eng. 2010, 58, 390–396. [Google Scholar] [CrossRef]
  14. Freeborn, T.J.; Elwakil, A.S.; Maundy, B. Compact wide frequency range fractional-order models of human body impedance against contact currents. Math. Probl. Eng. 2016, 2016, 4967937. [Google Scholar] [CrossRef]
  15. Butera, S.; Di Paola, M. A physically based connection between fractional calculus and fractal geometry. Ann. Phys. 2014, 350, 146–158. [Google Scholar] [CrossRef]
  16. Ionescu, C.M.; Machado, J.T.; De Keyser, R. Modeling of the lung impedance using a fractional-order ladder network with constant phase elements. IEEE Trans. Biomed. Circuits Syst. 2010, 5, 83–89. [Google Scholar] [CrossRef] [PubMed]
  17. Bigelow, M.E.G.; Jamieson, B.G.; Chui, C.O.; Mao, Y.; Shin, K.S.; Huang, T.J.; Huang, P.H.; Ren, L.; Adhikari, B.; Chen, J.; et al. Point-of-care technologies for the advancement of precision medicine in heart, lung, blood, and sleep disorders. IEEE J. Transl. Eng. Health Med. 2016, 4, 1–10. [Google Scholar] [CrossRef] [PubMed]
  18. Moore, K.L.; Dalley, A.F.; Agur, A.M. Clinically Oriented Anatomy; Lippincott Williams & Wilkins: Philadelphia, PA, USA, 2013. [Google Scholar]
  19. Drake, R.; Vogl, A.W.; Mitchell, A.W. Gray’s Anatomy for Students E-Book; Elsevier Health Sciences: Amsterdam, The Netherlands, 2009. [Google Scholar]
  20. Ionescu, C.M.; Segers, P.; De Keyser, R. Mechanical properties of the respiratory system derived from morphologic insight. IEEE Trans. Biomed. Eng. 2008, 56, 949–959. [Google Scholar] [CrossRef]
  21. Ionescu, C.M.; Muntean, I.; Tenreiro-Machado, J.; De Keyser, R.; Abrudean, M. A theoretical study on modeling the respiratory tract with ladder networks by means of intrinsic fractal geometry. IEEE Trans. Biomed. Eng. 2009, 57, 246–253. [Google Scholar] [CrossRef]
  22. Ionescu, C.; Desager, K.; De Keyser, R. Fractional order model parameters for the respiratory input impedance in healthy and in asthmatic children. Comput. Methods Programs Biomed. 2011, 101, 315–323. [Google Scholar] [CrossRef]
  23. Ionescu, C.; Kelly, J.F. Fractional calculus for respiratory mechanics: Power law impedance, viscoelasticity, and tissue heterogeneity. Chaos Solitons Fractals 2017, 102, 433–440. [Google Scholar] [CrossRef]
  24. Assadi, I.; Charef, A.; Copot, D.; De Keyser, R.; Bensouici, T.; Ionescu, C. Evaluation of respiratory properties by means of fractional order models. Biomed. Signal Process. Control 2017, 34, 206–213. [Google Scholar] [CrossRef]
  25. Copot, D.; De Keyser, R.; Derom, E.; Ionescu, C. Structural changes in the COPD lung and related heterogeneity. PLoS ONE 2017, 12, e0177969. [Google Scholar] [CrossRef]
  26. Copot, D.; De Keyser, R.; Derom, E.; Ortigueira, M.; Ionescu, C.M. Reducing bias in fractional order impedance estimation for lung function evaluation. Biomed. Signal Process. Control 2018, 39, 74–80. [Google Scholar] [CrossRef]
  27. Jossinet, J. Elementary electrodynamics. Technol. Health Care 2008, 16, 465–474. [Google Scholar] [CrossRef] [PubMed]
  28. Caponetto, R.; Di Pasquale, G.; Graziani, S.; Murgano, E.; Pollicino, A. Realization of Green Fractional Order Devices by using Bacterial Cellulose. AEU Int. J. Electron. Commun. 2019, 112, 152927. [Google Scholar] [CrossRef]
  29. Valsa, J.; Vlach, J. RC models of a constant phase element. Int. J. Circuit Theory Appl. 2013, 41, 59–67. [Google Scholar] [CrossRef]
  30. Tsirimokou, G.; Psychalinos, C.; Elwakil, A.S.; Salama, K.N. Electronically tunable fully integrated fractional- order resonator. IEEE Trans. Circuits Syst. II Express Briefs 2017, 65, 166–170. [Google Scholar] [CrossRef] [Green Version]
  31. Kapoulea, S.; Tsirimokou, G.; Psychalinos, C.; Elwakil, A.S. Generalized Fully Adjustable Structure for Emulating Fractional-Order Capacitors and Inductors of Orders less than Two. Circuits Syst. Signal Process. 2020, 39, 1797–1814. [Google Scholar] [CrossRef]
  32. Tsirimokou, G.; Kartci, A.; Koton, J.; Herencsar, N.; Psychalinos, C. Comparative study of discrete component realizations of fractional-order capacitor and inductor active emulators. J. Circuits Syst. Comput. 2018, 27, 1850170. [Google Scholar] [CrossRef]
  33. Krishna, B. Studies on fractional order differentiators and integrators: A survey. Signal Process. 2011, 91, 386–426. [Google Scholar] [CrossRef]
  34. Bertsias, P.; Psychalinos, C.; Maundy, B.J.; Elwakil, A.S.; Radwan, A.G. Partial fraction expansion-based realizations of fractional-order differentiators and integrators using active filters. Int. J. Circuit Theory Appl. 2019, 47, 513–531. [Google Scholar] [CrossRef]
  35. Mohan, P.A. VLSI Analog Filters: Active RC, OTA-C, and SC; Springer Science & Business Media: Berlin, Germany, 2012. [Google Scholar]
  36. Corbishley, P.; Rodriguez-Villegas, E. A nanopower bandpass filter for detection of an acoustic signal in a wearable breathing detector. IEEE Trans. Biomed. Circuits Syst. 2007, 1, 163–171. [Google Scholar] [CrossRef]
  37. Gray, P.R.; Hurst, P.J.; Lewis, S.H.; Meyer, R.G. Analysis and Design of Analog Integrated Circuits; John Wiley & Sons: Hoboken, NJ, USA, 2009. [Google Scholar]
  38. Tsirimokou, G.; Psychalinos, C.; Elwakil, A.S. Digitally programmed fractional-order Chebyshev filters realizations using current-mirrors. In Proceedings of the IEEE International Symposium on Circuits and Systems (ISCAS), Lisbon, Portugal, 24–27 May 2015; pp. 2337–2340. [Google Scholar]
  39. Bertsias, P.; Khateb, F.; Kubanek, D.; Khanday, F.A.; Psychalinos, C. Capacitorless digitally programmable fractional-order filters. AEU Int. J. Electron. Commun. 2017, 78, 228–237. [Google Scholar] [CrossRef]
  40. Tsirimokou, G.; Sotner, R.; Jerabek, J.; Koton, J.; Psychalinos, C. Programmable analog array of fractional-order filters with CFOAs. In Proceedings of the 40th International Conference on Telecommunications and Signal Processing (TSP), Barcelona, Spain, 5–7 July 2017; pp. 706–709. [Google Scholar]
  41. Kapoulea, S.; Psychalinos, C.; Elwakil, A.S.; Radwan, A.G. One-terminal electronically controlled fractional- order capacitor and inductor emulator. AEU Int. J. Electron. Commun. 2019, 103, 32–45. [Google Scholar] [CrossRef]
  42. Gholizadeh, A.; Voiry, D.; Weisel, C.; Gow, A.; Laumbach, R.; Kipen, H.; Chhowalla, M.; Javanmard, M. Toward point-of-care management of chronic respiratory conditions: Electrochemical sensing of nitrite content in exhaled breath condensate using reduced graphene oxide. Microsyst. Nanoeng. 2017, 3, 1–8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Fractional-order electrical equivalent of the impedance of the lung in the human respiratory tree [16].
Figure 1. Fractional-order electrical equivalent of the impedance of the lung in the human respiratory tree [16].
Jlpea 10 00018 g001
Figure 2. Functional block diagram for emulating a fractional-order capacitor of order α = ( 0 , 1 ) [30].
Figure 2. Functional block diagram for emulating a fractional-order capacitor of order α = ( 0 , 1 ) [30].
Jlpea 10 00018 g002
Figure 3. Functional block diagram for emulating a fractional-order capacitor of order α = ( 0 , 2 ) [31].
Figure 3. Functional block diagram for emulating a fractional-order capacitor of order α = ( 0 , 2 ) [31].
Jlpea 10 00018 g003
Figure 4. FBD of an Inverse-Follow-the-Leader-Feedback topology for approximating a fractional-order differentiator using the Continued Fraction Expansion tool [30].
Figure 4. FBD of an Inverse-Follow-the-Leader-Feedback topology for approximating a fractional-order differentiator using the Continued Fraction Expansion tool [30].
Jlpea 10 00018 g004
Figure 5. FBD for approximating a fractional-order differentiator using the Partial Fraction Expansion tool [34].
Figure 5. FBD for approximating a fractional-order differentiator using the Partial Fraction Expansion tool [34].
Jlpea 10 00018 g005
Figure 6. Implementation of the fractional-order differentiation stage, required for the implementation of the functional block diagram in Figure 3.
Figure 6. Implementation of the fractional-order differentiation stage, required for the implementation of the functional block diagram in Figure 3.
Jlpea 10 00018 g006
Figure 7. Implementation of the integer-order differentiation stage, required for the implementation of the functional block diagram in Figure 3 [35].
Figure 7. Implementation of the integer-order differentiation stage, required for the implementation of the functional block diagram in Figure 3 [35].
Jlpea 10 00018 g007
Figure 8. Implementation of the V / I converter stage using an appropriately configured multiple- output OTA.
Figure 8. Implementation of the V / I converter stage using an appropriately configured multiple- output OTA.
Jlpea 10 00018 g008
Figure 9. Implementation of a floating resistor using an appropriately configured multiple-output OTA.
Figure 9. Implementation of a floating resistor using an appropriately configured multiple-output OTA.
Jlpea 10 00018 g009
Figure 10. CMOS OTA circuitry used for implementing the topologies in Figure 6, Figure 7, Figure 8 and Figure 9 [36].
Figure 10. CMOS OTA circuitry used for implementing the topologies in Figure 6, Figure 7, Figure 8 and Figure 9 [36].
Jlpea 10 00018 g010
Figure 11. Layout design of the proposed emulator of the lung impedance model of the human respiratory tree in Figure 1 (fractional-order differentiator in red polygon, integer-order differentiator in green polygon, V / I converter in yellow polygon, and active resistor in blue polygon).
Figure 11. Layout design of the proposed emulator of the lung impedance model of the human respiratory tree in Figure 1 (fractional-order differentiator in red polygon, integer-order differentiator in green polygon, V / I converter in yellow polygon, and active resistor in blue polygon).
Jlpea 10 00018 g011
Figure 12. Post-layout impedance (a) magnitude, and (b) phase frequency responses of the CPE emulator implemented according to the FBDs in Figure 2 and Figure 3.
Figure 12. Post-layout impedance (a) magnitude, and (b) phase frequency responses of the CPE emulator implemented according to the FBDs in Figure 2 and Figure 3.
Jlpea 10 00018 g012
Figure 13. Post-layout frequency responses of the lung impedance (a) magnitude, and (b) phase.
Figure 13. Post-layout frequency responses of the lung impedance (a) magnitude, and (b) phase.
Jlpea 10 00018 g013
Figure 14. Statistical histograms about the sensitivity of the (a) magnitude, and (b) phase of the lung impedance, at the center frequency of the approximation ( f 0 = 2.5 Hz).
Figure 14. Statistical histograms about the sensitivity of the (a) magnitude, and (b) phase of the lung impedance, at the center frequency of the approximation ( f 0 = 2.5 Hz).
Jlpea 10 00018 g014
Table 1. Values of the elements of the fractional-order electrical equivalent of the impedance of the lung in the human respiratory tree [16].
Table 1. Values of the elements of the fractional-order electrical equivalent of the impedance of the lung in the human respiratory tree [16].
Case R ( M Ω ) C α (nF/s 1 α )Order of Element ( α )
#16.442380.5808
#24.29290.40.7978
#36.021350.802
#49.542060.81
#54.891780.9349
#65.9317.50.95
#76.06261.134
Table 2. Values of the elements of the fractional-order emulator in Figure 3.
Table 2. Values of the elements of the fractional-order emulator in Figure 3.
Case C α (nF/s 1 α )Order of Element ( α ) Order of Differentiator ( q ) g m , VI ( μ S)
#12380.58080.58081.18
#290.40.79780.79780.81
#31350.8020.8021.23
#42060.810.811.92
#51780.93490.93492.34
#617.50.950.950.24
#7261.1340.1340.59
Table 3. Values of scaling factors and time-constants for approximating fractional-order differentiators of the FBD in Figure 3, using the FBD in Figure 4.
Table 3. Values of scaling factors and time-constants for approximating fractional-order differentiators of the FBD in Figure 3, using the FBD in Figure 4.
CaseOrder of Differentiator ( q ) K 0 K 1 K 2 τ 1 (ms) τ 2 (ms)
#10.58080.14616.865.17114.3
#20.79780.048120.692.3085.14
#30.8020.047121.292.2584.65
#40.810.044122.52.1583.71
#50.93490.012181.90.7170.09
#60.950.0091109.60.5468.56
#70.1340.66811.512.9209.5
Table 4. Values of scaling factors and time-constants for approximating fractional-order differentiators of the FBD in Figure 3, using the FBD in Figure 5.
Table 4. Values of scaling factors and time-constants for approximating fractional-order differentiators of the FBD in Figure 3, using the FBD in Figure 5.
CaseOrder of Differentiator ( q ) K 0 K 1 K 2 τ 1 (ms) τ 2 (ms)
#10.58086.86−6.12−0.595.43108.9
#20.797820.69−20.24−0.42.3782.78
#30.80221.29−20.85−0.392.3182.34
#40.8122.5−22.07−0.382.2181.50
#50.934981.9−81.73−0.160.7169.37
#60.95109.6−109.4−0.120.5468.02
#70.1341.5−0.51−0.3213.8195.7
Table 5. Values of dc bias currents for approximating the fractional-order differentiator, with unity gain frequency f 0 = 2.5 Hz, in Figure 6 and implementing the V / I converter and resistor in Figure 8 and Figure 9.
Table 5. Values of dc bias currents for approximating the fractional-order differentiator, with unity gain frequency f 0 = 2.5 Hz, in Figure 6 and implementing the V / I converter and resistor in Figure 8 and Figure 9.
CaseDifferentiator V / I ConverterResistor
Order ( q ) I B 1 (pA) I B 2 (pA) I B (pA) I B , VI (mA) I B , R (mA)
#10.580851.736.1350668.7
#20.7978118.547.535045.613.1
#30.802121.647.735069.19.3
#40.81127.148.2300107.85.9
#50.9349395.556.7700131.311.5
#60.9552057.860013.59.5
#70.13420.320.160033.29.3
Table 6. Aspect ratio values of MOS transistors of the OTA in Figure 10 used for implementing the topologies in Figure 6, Figure 7, Figure 8 and Figure 9.
Table 6. Aspect ratio values of MOS transistors of the OTA in Figure 10 used for implementing the topologies in Figure 6, Figure 7, Figure 8 and Figure 9.
TransistorFigure 6Figure 7Figure 8Figure 9
( μ m/ μ m)( μ m/ μ m)( μ m/ μ m)( μ m/ μ m)
M1n-M2n8/151/51/12/1
M3n-M4n40/155/55/110/1
M5n-M7n25/1525/515/1510/5
M1p-M2p10/50.5/525/105/15
Table 7. Simulated and theoretical (in parentheses) values of the magnitude and phase (at f 0 = 2.5 Hz) of the differentiator in Figure 6, of the impedance of the V / I converter in Figure 8, and resistor in Figure 9, used for implementing the FBDs in Figure 2 and Figure 3.
Table 7. Simulated and theoretical (in parentheses) values of the magnitude and phase (at f 0 = 2.5 Hz) of the differentiator in Figure 6, of the impedance of the V / I converter in Figure 8, and resistor in Figure 9, used for implementing the FBDs in Figure 2 and Figure 3.
CaseOrder ( q ) Differentiator V / I ConverterResistor
Gain arg ( Z ) ( ) 1 / g m , VI ( M Ω ) R ( M Ω )
#10.58081 (1)53.7 (52.3)0.85 (0.84)6.45 (6.44)
#20.79781 (1)71.8 (71.8)1.22 (1.23)4.27 (4.29)
#30.8021 (1)71.7 (72.18)0.81 (0.81)6.01 (6.02)
#40.811 (1)72.7 (72.9)0.51 (0.52)9.54 (9.54)
#50.93491 (1)83.7 (84.14)0.43 (0.43)4.89 (4.89)
#60.951 (1)84.2 (85.5)0.41 (0.42)5.92 (5.93)
#70.1341 (1)12.5 (11.5)1.69 (1.69)6.04 (6.06)
Table 8. Simulated and theoretical (in parentheses) values of impedance magnitude and phase of the model in Figure 1, at the center frequency of the approximation ( f 0 = 2.5 Hz).
Table 8. Simulated and theoretical (in parentheses) values of impedance magnitude and phase of the model in Figure 1, at the center frequency of the approximation ( f 0 = 2.5 Hz).
CaseOrder of Element ( α ) Z ( M Ω ) arg ( Z ) ( )
#10.58086.9 (6.9)−5.64 (−6.13)
#20.79784.70 (4.71)−13.12 (−14.7)
#30.8026.25 (6.24)−6.26 (−7.3)
#40.819.68 (9.66)−2.68 (−3)
#50.93494.96 (4.94)−5.35 (−4.96)
#60.957.83 (7.42)−32.5 (−34)
#71.1345.92 (5.96)−17.1 (−16.2)

Share and Cite

MDPI and ACS Style

Kaskouta, E.; Kapoulea, S.; Psychalinos, C.; Elwakil, A.S. Implementation of a Fractional-Order Electronically Reconfigurable Lung Impedance Emulator of the Human Respiratory Tree. J. Low Power Electron. Appl. 2020, 10, 18. https://0-doi-org.brum.beds.ac.uk/10.3390/jlpea10020018

AMA Style

Kaskouta E, Kapoulea S, Psychalinos C, Elwakil AS. Implementation of a Fractional-Order Electronically Reconfigurable Lung Impedance Emulator of the Human Respiratory Tree. Journal of Low Power Electronics and Applications. 2020; 10(2):18. https://0-doi-org.brum.beds.ac.uk/10.3390/jlpea10020018

Chicago/Turabian Style

Kaskouta, Elpida, Stavroula Kapoulea, Costas Psychalinos, and Ahmed S. Elwakil. 2020. "Implementation of a Fractional-Order Electronically Reconfigurable Lung Impedance Emulator of the Human Respiratory Tree" Journal of Low Power Electronics and Applications 10, no. 2: 18. https://0-doi-org.brum.beds.ac.uk/10.3390/jlpea10020018

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

Article Metrics

Back to TopTop