Using the fractional calculus approach, we present the Laplace analysis of an equivalent electrical circuit for a multilayered system, which includes distributed elements of the Cole model type. The Bode graphs are obtained from the numerical simulation of the corresponding transfer functions using arbitrary electrical parameters in order to illustrate the methodology. A numerical Laplace transform is used with respect to the simulation of the fractional differential equations. From the results shown in the analysis, we obtain the formula for the equivalent electrical circuit of a simple spectrum, such as that generated by a real sample of blood tissue, and the corresponding Nyquist diagrams. In addition to maintaining consistency in adjusted electrical parameters, the advantage of using fractional differential equations in the study of the impedance spectra is made clear in the analysis used to determine a compact formula for the equivalent electrical circuit, which includes the Cole model and a simple RC model as special cases.
- fractional calculus
- Nyquist and Bode diagrams
- numerical Laplace transform
According to Rigaud  who at the beginning of the 20th century began to study the structure of biological tissues based on their electrical properties, biological tissues are conductors and their resistance varies with frequency. The electrical property of any biological tissue depends on its intrinsic structure. In the case of human skin, the impedance can vary with the thickness and moisture content of the organ, the concentration and activity of sweat glands, injuries, age of subject and environmental factors such as temperature and humidity. Electrical impedance studies in biological systems, including human skin, generally, relate to direct measurements of impedance and phase angle as functions of frequency, voltage, or current applied [2, 3, 4, 5, 6]. In 1974, Burton  applied Bode analysis to measurements of impedance and phase angle of the skin. Through this method, a passive equivalent circuit can be used and considered as a "black box" to plot its impedance and phase angle versus frequency. The only necessary assumption is that the system consists entirely of linear passive elements [7, 8]. Although the resulting model is not necessarily unique , it describes the system with great precision in the range of frequencies studied.
Electrical impedance spectroscopy (EIS) has proven useful in the characterization of biomaterials by recording the behavior of their intrinsic properties while evaluating the bioelectrical response of the system as a sinusoidal excitation signal is applied . It has also been used to measure biological tissues in which the electrical impedance depends on water content and ionic conduction in the body. It should be noted that the terms “electrical impedance” and “resistance” are often used without distinction in literature . In terms of frequency, some researchers have reported that when the frequency is less than or equal to 10 kHz, the current does not cross the cell membrane, and thus the resistance obtained is relative only to the extracellular mass . This electrical conductivity is small in adipose tissue compared to fat-free tissue. The resistivity of components such as blood (150 Ωcm) or urine (40-200 Ωcm) is low, muscle (250-2000 Ωcm) and fat (2000-5000 Ωcm) are intermediate, and bone (10 000 Ωcm) is high. In the measurement of bioimpedance an electrical stimulus is applied and the response produced on a specific region of the body is analyzed. Usually, the stimulus is an alternating current signal of low amplitude intended to measure the electric field or potential difference generated between different parts of the tissue. The relationship between the data of the applied stimulus and the response obtained as a function of frequency provides the impedance spectrum of tissues studied .
Transfer function analysis is a mathematical approach to relate an input signal (or excitation) and the system response. The ratio formed by the pattern of the output and the input signal makes it possible to find the zeros and poles, respectively . A space state representation is a mathematical model of a physical system described by a set of inputs, outputs and state variables related by first order differential equations that are combined in a first order matrix differential equation. So as not to be affected by the number of inputs, outputs, and states, the variables are expressed as vectors, and algebraic equations are written in matrix form.
To analyze the behavior of the transfer function in the frequency domain, several graphical methods were used, such as Bode plots that provide a graphical representation of the magnitude and phase versus frequency of the transfer function . Nyquist plots are polar plots of impedance modulus and phase lag . One of the biggest drawbacks of working in polar coordinates is that the curve no longer retains its original shape when a modification to the system is made. For instance when the loop gain is altered, the magnitude-phase trace shifts vertically without distortion. However, when the properties of the phase are changed independently without affecting the gain, the magnitudephase trace is affected only in the horizontal direction.
The Cole impedance model was postulated in its final form by Kenneth Cole in 1940 . This impedance model is based on replacing the ideal capacitor in the Debye model  for a general element called constant phase element (CPE). The analysis of the Cole impedance model requires three things: (1) an equivalent circuit; (2) the development of the corresponding equations; and (3) a simulation which gives the complex impedance behavior.
The Cole model is represented by a series resistor (
Using the Cole model parameters, this model corresponds tissue impedance values with the values of the components of circuit models .
Fractional calculus (FC) is nearly as old as the conventional calculus, but is not very popular within many scientific and engineering communities. The idea of FC seems to be quite a strange topic and very hard to explain, due to the fact that it is not related to some important physical meaning, such as velocity and acceleration, unlike commonly used differential calculus. For this reason it has long been of interest only to mathematicians. On the other hand, many physical phenomena have "intrinsic" fractional order descriptions and so FC is necessary in order to explain them. In many applications FC provides more accurate models of the physical systems than ordinary calculus does. Because of its success in description of anomalous diffusion [16, 17, 18], non-integer order calculus both in one and multidimensional space, has become an important tool in many areas of physics, mechanics, engineering, and bioengineering [19, 20, 21, 22, 23, 24, 25, 26]. Fundamental physical considerations in favor of the use of models based on derivatives of non-integer order are given in Caputo & Mainardi and Westerlund [27-28]. Fractional derivatives provide an excellent instrument for the description of memory and hereditary properties of various materials and processes. This is the main advantage of FC in comparison with the classical integer-order models, in which such effects are in fact neglected. The other large field requiring the use of FC is the theory of fractals . The development of the theory of fractals has opened further perspective for the theory of fractional derivatives, especially in modeling dynamical processes in self-similar and porous structures.
To analyze the dynamical behavior of a fractional system it is necessary to use an appropriate definition of fractional derivative. In fact, the definitions of the fractional order derivative are not unique and there exist several definitions, including Grünwald-Letnikov, Riemann-Liouville, Weyl, Riesz, and the Caputo representation. In the Caputo case, the derivative of a constant is zero and we can properly define the initial conditions for the fractional differential equations so that they can be handled analogously to the classical integer case. Caputo derivative implies a memory effect by means of a convolution between the integer order derivative and a power of time [30, 31, , 32, 33]. For these reasons, in this paper we prefer to use the Caputo fractional derivative defined by:
From equation (3), using the initial conditions
This is consistent with the usual definition of the Laplace transform when
From (5) if
The Laplace transform is a useful tool for analyzing linear systems because it simplifies the problem of dealing with differential equations in the time domain by converting them into algebraic equations within the frequency domain. The numerical Laplace transform (NLT) is essentially a modified discrete Fourier transform (DFT) through a windowing function (Gibbs phenomenon) and a stability factor (aliasing) . Development of the NLT and its application to the analysis of systems has been well documented over the past 40 years [36-37]. However, its use has traditionally limited the analysis of problems where the result can be expressed in terms of simple functions that allow the use of tables of Laplace transforms . When using discrete techniques in the frequency domain computation time becomes an important factor since it requires a certain amount of time to transform the data from the frequency domain to the time domain or vice versa. However, by using the fast Fourier transform (FFT) the time necessary for computation is greatly reduced and as a result the techniques of analysis in the frequency domain become an attractive option.
Sheng investigated the validity of applying numerical inverse Laplace transform algorithms in fractional calculus . In a paper of Gómez an overview of a methodology based on the NLT and its application to analysis of a three physical system from the point of view of FC is discussed .
In this work for the multilayer biological system, we propose an analysis of an electrical circuit equivalent consisting of epidermis, dermis and the subcutaneous tissue. The response of the Bode diagram is interpreted and compared with the response to the modeling of integer order (using transfer function analysis) and fractional order (obtained through NLT techniques). The modeling of electrical impedance spectra applied to the study of experimental data of blood is made using the Cole model electrical impedance equations. The Nyquist diagrams (frequency response) are also compared for both integer and fractional order model responses.
In order to apply the general theory to the skin type system, we propose an equivalent electrical circuit for each component of the skin. Fig. 1 shows the schematic diagram to describe the proposed model of the skin, where
It is known that in biological systems, the module of the impedance is inversely dependent on frequency and does not show phenomena of conversion of electrical energy to magnetic energy [41-42]. Therefore, it is considered appropriate to model a biological system as only having capacitive and resistive behavior. Fig. 2 shows the circuit equivalent to the first layer of the system. Applying Kirchhoff's Laws the following representation of state equations is obtained:
Fig. 3 shows the corresponding Bode plot.
The circuit parameters are chosen by establishing that the array is not indeterminate. If any row of the matrix is zero, a branch of the circuit (Fig. 1) has been removed.
Following the above procedure, we next obtained the equations of state considering only the first two layers and finally for all three layers and their corresponding Bode plots, showing the effect of each layer. In equations (10) and (11) the matrices and vectors for
In each case, (12) and (13), the input passes through the circuit and reflects the voltage of the capacitor under study. Thus, the transfer function is related to a different output in each case for the same input. Fig. 4 shows the Bode diagram corresponding to the second and third layer.
Applying the formalism of fractional calculus to the circuit of Fig. 2 gives the fractional Laplace transform. For this purpose, the fractional time derivative operator of order γ can be written as follows
Expression (14) is not an ordinary (usual) time derivative operator because the dimensionality is
which is true if the parameter
Expression (15) is a time derivative in the usual sense, because its dimension is
The fractional differential equation for the circuit of Fig. 2 is
the parameter values appear as per unit values  and correspond to
The Laplace transform for (16) is
Applying NLT algorithm to the frequency response obtained is shown in Fig. 5. As a case study the values of γ=0.95 and γ=0.98 were used.
Fig. 6 shows the comparison for the Bode diagram for the integer order and the fractional order. The second and third layers had similar results.
Three samples were obtained from the blood transfusion center in Leon, Guanajuato.
All of the samples were clinically examined according to the blood transfusion protocol, to confirm they were free of contagious diseases such as HIV, hepatitis, and brucellosis. Cell populations were monitored through a counting study using a flow cytometer type Coulter. This was done so that cell mortality did not affect the results. To obtain the electrical impedance spectra a Solartron® 1260 was used. The fidelity of the frequency sweep for these tests was important since it shows the characteristic spectrum of the sample, which is necessary for a comparison with the electrical parameters of an equivalent circuit. The frequency range used was from 10Hz to 100kHz. The samples were placed in Bayer® strips. A frequency sweep with constant voltage amplitude of 25 mV was applied through two electrodes integrated in the Bayer® container. All test points were coated with silver and therefore had negligible contributions to polarization, and thus did not need to be compensated for. Once the measurements were represented in a Nyquist diagram, we obtained representations of equivalent electrical circuits via the software ZView®. The method of "instant adjustment" was used to fit the data to predefined circuit models. Fig. 7 shows the equivalent network for the measurement of erythrocytes, leukocytes and plasma.
Table 1 shows the values for the erythrocytes, leukocytes and plasma, for the samples 1, 2, and 3, respectively. The units of resistance are in ohms and the capacitors are farads.
Values for the erythrocytes, leukocytes and plasma.
To determine the equivalent equation, the impedance is found by the following formula in the complex frequency domain:
Applying Kirchhoff laws to the circuit of Fig. 7, we have:
Before applying the Laplace transform of (19) and (20) the following considerations must be made
Working with our definition of the fractional derivative , (21) and (22) become:
Substituting (23) and (24) into (20), we obtain:
Applying the Laplace transform to (19) and (25):
Finally from (26) and (27) we have:
Equation (29) is the result of the fractional temporal operator in the equation for the RC equivalent circuit; this general formula includes an arbitrary constant,
In this work we have analyzed in detail the transfer function of a multilayer system. Fig. 3 shows the Bode plot to the magnitude (top graph) and phase (bottom graph). In the first graph of Fig. 3, we see by increasing frequency, from the cutoff frequency (30 rad/s), the magnitude decreases at a rate of 20 dB/decade, while for frequencies below the cutoff magnitude it is almost constant.
This means that the current through the dermis is higher with decreasing frequency and is lower as the frequency increases. This means that the current through the dermis increases for higher frequencies. The attenuation in the first layer is solely due to the resistance of the electrodes. The current flowing through these resistors at low frequencies causes a Joule effect, which raises the temperature and therefore the kinetic energy of the molecules that make up the layer . As the frequency increases, other phenomena occur, such as displacement currents, ionization, polarization, and so on. Concerning phase, we have shown that by increasing frequency the displacement current and polarization are also increased (Fig. 3, bottom graph). This causes a decrease in phase 0° to -90° for very high frequencies.
There is not a significant decrease in the rate of change of the magnitude on the frequency, as can be seen from equations (9) and (10), where the difference in poles and zeros is two. By increasing the number of layers, the output voltage in the layer farthest from the source is attenuated at a faster rate due to the resistance of the previous layers and the increase in frequency. Thus the current in the deeper layers of the skin decreases more rapidly with increasing frequency. Fig. 4 shows that adding layers to the model changes the cutoff frequency – hence the shift in the phase – and the magnitude decreases from 20 dB/decade to 40 dB/decade. For the range of frequency 10 to 100 rad/s a shift in the magnitude and phase is shown, which implies that the proposed circuit better defines these frequencies.
In the first graph of Fig. 5, we see that by increasing frequency from the cutoff frequency (10 rad/s), the magnitude decreases at a rate of 6 dB/decade, while for frequencies below the cutoff magnitude it is almost constant. This means that the current through the dermis is very high with decreasing frequency and is very low when the frequency increases. The current flowing through these resistors at low frequencies causes a fractional Joule effect (behavior between a system conservative and dissipative), which raises the temperature and therefore the kinetic energy of the molecules that make up the layer. Concerning the phase Fig. 5 (bottom graph), we have that by increasing frequency the displacement current and polarization are very much increased.
With respect to the data of blood tissue, the families of curves of erythrocytes, leukocytes, and blood plasma are found to have similar parameters (
Exponent of the fractional differential equation that best fits the data
Fractional calculus has been used successfully to modify many existing models of physical processes. The representation of equivalent models in integer order derivatives provided a good approximation of the bioelectric response of the model. However, with the formal introduction of fractional calculus to the study of fractional derivative systems have a better approximation of this response. This is due in part to the nature of systems described by fractional calculus.
On the basis of Cole’s proposal to add a degree of extra freedom to solve the RC circuit for characterization purposes and to improve the correlation in the adjustment to experimental data, we have developed analytical arguments to derive this result based on integration in weighted individual relaxation processes. However, the distributions of relaxation times involve complex functions and are difficult to measure. This study has shown a pattern in which the Cole type behavior appears as a result of competition between a capacitive and resistive behavior within the sample, characterized by the fractional order derivative of the applied voltage. The new generalized model includes the Cole model and the simple RC circuit as particular cases.
The models for EIS are assumed to be linear in their first approximation. The electrical parameters only take a nonlinear behavior in the case of tissue damage due to the excessive power in the supply or in the presence of physical and chemical reactions in the sample induced by the input current (exothermic processes or release of electrons). The electrical conduction, even for alternating current, has impedance that depends on the temperature, which turns the frequency response into a function of temperature.
Exponent of the fractional differential equation that best fits the data
Values for the erythrocytes, leukocytes and plasma.