Next Article in Journal
Process Modeling, Optimization and Cost Analysis of a Sulfur Recovery Unit by Applying Pinch Analysis on the Claus Process in a Gas Processing Plant
Next Article in Special Issue
A Hybrid Localized Meshless Method for the Solution of Transient Groundwater Flow in Two Dimensions
Previous Article in Journal
Multi-UAV Coverage Path Planning Based on Hexagonal Grid Decomposition in Maritime Search and Rescue
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Vibrational Resonance and Electrical Activity Behavior of a Fractional-Order FitzHugh–Nagumo Neuron System

College of Mechanics and Materials, Hohai University, Nanjing 211100, China
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(1), 87; https://doi.org/10.3390/math10010087
Submission received: 13 December 2021 / Revised: 22 December 2021 / Accepted: 23 December 2021 / Published: 27 December 2021

Abstract

:
Making use of the numerical simulation method, the phenomenon of vibrational resonance and electrical activity behavior of a fractional-order FitzHugh–Nagumo neuron system excited by two-frequency periodic signals are investigated. Based on the definition and properties of the Caputo fractional derivative, the fractional L1 algorithm is applied to numerically simulate the phenomenon of vibrational resonance in the neuron system. Compared with the integer-order neuron model, the fractional-order neuron model can relax the requirement for the amplitude of the high-frequency signal and induce the phenomenon of vibrational resonance by selecting the appropriate fractional exponent. By introducing the time-delay feedback, it can be found that the vibrational resonance will occur with periods in the fractional-order neuron system, i.e., the amplitude of the low-frequency response periodically changes with the time-delay feedback. The weak low-frequency signal in the system can be significantly enhanced by selecting the appropriate time-delay parameter and the fractional exponent. In addition, the original integer-order model is extended to the fractional-order model, and the neuron system will exhibit rich dynamical behaviors, which provide a broader understanding of the neuron system.

1. Introduction

Vibrational resonance (VR) has attracted considerable attention in the field of nonlinear sciences in the last twenty years. Based on the study of stochastic resonance (SR) [1], VR is firstly proposed by Landa and McClintock [2]. When replacing the noise in SR with an appropriate high-frequency signal, the weak low-frequency signal can be greatly amplified, which is similar to the typical “Inverted Bell” resonance phenomenon in SR, and is named as VR. Compared with the noise in SR, the high-frequency signal in VR is more controllable. Biharmonic signals are common in various fields, such as acoustics [3], optics [4], engineering [5], neuroscience [6], etc. Recently, the research hotspots of VR have changed from classic bistable systems [7,8,9] to fractional-order systems [10], delay systems [11], and network dynamical systems [12].
Due to the complex definition and lack of corresponding application background, the research on fractional calculus has been limited to the field of mathematics for a long time. However, compared with integer calculus, the power-law characteristic of complex social and physical phenomena can be accurately approximated through fractional calculus. Hence, this theory is gradually used in viscoelastic materials [13], electrification process [14], control theory [15], and neuron models [16], etc. To the best of our knowledge, the research on VR for fractional-order systems is quite few, and most of which are limited to bistable and multi-stable systems [17,18,19]. However, it is found that fractional calculus has its own unique advantages in describing certain neuronal characteristics. For example, fractional-order differentiation can be used to account for the firing rate of neocortical pyramidal neurons when stimulated by sinusoidal current [20]. Anastasio et al. [21] thought that the net output of the motor neurons in the visual system is consistent with fractional-order differentiation relative to eye position. Therefore, it is of great significance to study VR in the fractional-order neuron system.
In this study, we consider the FitzHugh–Nagumo (FHN) neuron model. As one of the simplest mathematical models for disclosing the dynamical behavior of neurons, the FHN model is widely used in integer-order systems for studying VR [22,23]. In brain activities [24], the neurons may exhibit two quite different time scales, and then, it is reasonable to reveal the mechanism of weak signal detection of neuron by VR. Since time delay is inherent in the neuron system, great progress has been made in the research on the effect of time delay on VR. Adding the time-delay feedback to the recovery variable of the FHN neuron model, the phenomenon of multiple VR can be induced in the neuron [25]. It is found that the bifurcation point and equilibrium point change periodically with the increase of time delay in a fractional order quantic oscillator system, and the output can be enhanced by selecting appropriate time delay [26].
Inspired by the above-mentioned ideas, the effects of fractional order and time delay on VR and dynamical behavior are studied in the fractional-order FHN neuron model. The VR can be induced in the fractional-order neuron model without strict requirement for the amplitude of high-frequency signal. It can be found that multiple VR occurs in the neuron with the increase of delay. Compared with integer-order model, the fractional-order FHN neuron exhibits rich electrical activities. The remainder of this paper is organized as follows: In Section 2, the fractional-order FHN neuron model excited by two periodic signals is briefly introduced. In Section 3, the VR in the fractional-order FHN model without and with time delay is studied. In Section 4, the effects of fractional exponent and time delay on the dynamical behavior in the fractional-order FHN model are discussed. Several conclusions are given in Section 5.

2. Fractional-Order FitzHugh–Nagumo Neuron Model

There are three common definitions of fractional calculus, namely Riemann–Liouville (R-L) definition, Grünwald–Letnikov (G-L) definition, and the Caputo definition. With the property of supersingularity, R-L definition is mainly used for the analysis of mathematical theory and is not convenient for engineering and physical modeling. G-L definition can be regarded as the extension of the limit form of integral calculus difference definition, which is widely applied to early numerical calculation. As its initial condition is the form of integer calculus, and it has clear physical meaning, the Caputo definition is used in this paper.
The fractional derivative of univariate function f ( t ) is defined as
d α f ( t ) d t α = D t α t 0 f ( t ) = 1 Γ ( n α ) t 0 t f ( n ) ( τ ) d τ ( t τ ) α n + 1 ,   ( n 1 < α < n ) ,
where Γ ( n α ) is the gamma function, t 0 and t are the lower limit and upper limit of the definite integral, n is the minimum positive integer greater than α , and f ( n ) ( τ ) is the n t h derivative of the function f ( τ ) . Under the joint excitation of harmonic signals, the fractional-order FHN neuron model with time delay (the model comes from the model involved in Ref. [27]) is given by the following form:
ε d α x d t α = x ( t ) x ( t ) 3 3 y ( t ) , d α y d t α = x ( t ) + a + f cos ( ω t ) + F cos ( Ω t ) + K ( y ( t τ ) y ( t ) ) ,
where x ( t ) and y ( t ) represent the fast-varying membrane potential and the slow-varying recovery variable of neuronal cells, respectively; α ( 0 < α 1 ) is the fractional exponent; K is the strength of time-delay feedback. τ 0 is the delay parameter, and Equation (2) degenerates into the fractional-order FHN neuron model without delay when τ = 0 . ε = 0.05 is the time scale ratio, which is chosen to ensure that the membrane potential x ( t ) evolves faster than the recovery variable y ( t ) . f cos ( ω t ) and F cos ( Ω t ) represent the low-frequency signal and the high-frequency signal, respectively, which satisfy f < < 1 , ω < < Ω . The value of a determines the behavior of the system under the conditions that f = F = τ = 0 and α = 1 . If | a | > 1.0 , the system is excitable and has only a stable fixed point; if | a | < 1.0 , a limit cycle in the system arises. Here, the parameter a = 1.05 is chosen to make the system in the excitable state [28]. The variables in the model are dimensionless.
We consider the fractional derivative of x ( t ) defined with the Caputo fractional derivative,
d α x ( t ) d t α = f ( x , t ) .
The discrete format of fractional-order L1 algorithm is [29]:
d α x ( t ) d t α ( d t ) α Γ ( 2 α ) [ k = 0 N 1 [ x ( t k + 1 ) x ( t k ) ] [ ( N k ) 1 α ( N 1 k ) 1 α ] ] ,
where t k = k Δ t . Combining the right sides of Equations (3) and (4) and solving for x at time t N , it can be concluded that the discrete format of Equation (3) is
x ( t N ) ( d t ) α Γ ( 2 α ) f ( x , t ) + x ( t N 1 ) [ k = 0 N 2 [ x ( t k + 1 ) x ( t k ) ] [ ( N k ) 1 α ( N 1 k ) 1 α ] ] .
where the Markov term weighted by the gamma function is given by
( d t ) α Γ ( 2 α ) f ( x , t ) + x ( t N 1 )
and the memory trace is given by
k = 0 N 2 [ x ( t k + 1 ) x ( t k ) ] [ ( N k ) 1 α ( N 1 k ) 1 α ]
The memory trace integrates information of all the previous activities and has a memory effect, which is the typical property of fractional-order system. When α = 1 , the memory trace has no effect, and Equation (5) degenerates into classical Euler algorithm.
According to Equation (5), the discrete format of Equation (2) can be obtained as follows:
x ( t N ) 1 ε ( d t ) α Γ ( 2 α ) ( x ( t N 1 ) ( x ( t N 1 ) ) 3 3 y ( t N 1 ) ) + x ( t N 1 ) [ k = 0 N 2 [ x ( t k + 1 ) x ( t k ) ] [ ( N k ) 1 α ( N 1 k ) 1 α ] ] , y ( t N ) ( d t ) α Γ ( 2 α ) ( x ( t N 1 ) + a + f cos ( ω t N 1 ) + F cos ( Ω t N 1 ) + K ( y ( t N m ) y ( t N 1 ) ) ) + y ( t N 1 ) [ k = 0 N 2 [ y ( t k + 1 ) y ( t k ) ] [ ( N k ) 1 α ( N 1 k ) 1 α ] ] ,
where m is the number of discrete points caused by time delay.

3. Vibrational Resonance in the Fractional-Order FHN Neuron Model

3.1. VR in the Fractional-Order FHN Neuron Model without Time Delay

The response amplitude Q at the low-frequency signal is usually used as the index to measure the VR, which is defined by:
Q = Q s 2 + Q c 2 ,
with
Q s = 2 n T 0 n T x ( t ) sin ( ω t ) d t , Q c = 2 n T 0 n T x ( t ) cos ( ω t ) d t ,
where T = 2 π / ω and n is a positive integer. For the numerical simulation, the parameters are selected as ε = 0.05 ,   K = 0.2 ,   f = 0.01 ,   n = 50 .
Under fixed parameters ω = 0.5 ,   Ω = 5 , the curves of VR in the fractional-order FHN neuron without time delay are plotted in Figure 1a. It can be seen that different fractional exponents correspond to different phenomena of the VR. In Figure 1b, with the decrease of α , the region of the VR in the neuron changes, and the amplitude of low-frequency Q corresponding to the optimal F gradually decreases. Compared with the integer-order FHN neuron model, it is found that the low-frequency signal can also be amplified without strict requirement for the amplitude of high-frequency F in the fractional-order neuron. For example, in Figure 1b, the integer-order FHN neuron model fails to reach the state of the VR for F = 0.15 , while for α = 0.96 , the VR occurs for the same value of F . Hence, it is shown that the VR in the FHN neuron can be induced by adjusting appropriate fractional exponent and the high-frequency force.
In order to further discuss the effect of α on the VR, the response amplitude Q max versus α for F [ 0.1 , 0.25 ] is depicted in Figure 2a. For a certain range of F , although the response amplitude Q max is not a strictly monotonic increasing function of α , in general, with the increase of α , the response amplitude Q can be optimized under the appropriate high-frequency force. However, when the value of F is fixed, the function of Q versus α shows different monotonic characteristics, as shown in Figure 2b. For F = 0.11 , the curve of response amplitude is approximately a straight line, which indicates that VR does not occur with the change of α , while for F = 0.15 or F = 0.2 , the response amplitude Q is a nonlinear function of α , and the low-frequency signal can be significantly enhanced by selecting an appropriate fractional-order α compared with the integer-order system.

3.2. Multiple VR in Fractional-Order FHN Neuron Model with Delay

For the fractional-order FitzHugh–Nagumo neuron with time delay, the numerical result of functional curves of the response amplitude Q versus τ can be obtained by combining Equation (8) with Equation (9). The response amplitude Q versus time delay τ in the fractional-order FHN neuron excited by two frequency signals is given in Figure 3. Figure 3e,f are partial, enlarged views of Figure 3a,d, respectively. From Figure 3a–d, it is seen that multiple resonance occurs in the neuron with the increase of the time-delay parameter. Therefore, the response amplitude Q can reach the maximum by selecting appropriate time-delay parameters, and the response amplitude Q can be greatly amplified compared with the neuron without time delay. Another notable phenomenon is that the response amplitude Q is periodic with the change of τ . From Figure 3a–d, it is clear that that the period of response amplitude Q is 2 π / ω , while Figure 3e,f shows that the response amplitude Q varies with another period of 2 π / Ω . Therefore, the response amplitude Q in the fractional-order FHN neuron with time-delay feedback presents two different periods, namely fast period 2 π / Ω and slow period 2 π / ω , which exactly correspond to the period of high-frequency signal and low-frequency signal. Utilizing the periodicity of response amplitude Q and selecting appropriate time-delay parameters are helpful to realizing the effective control of the fractional-order FHN neuron. The only regret is that the parameters involved in the system are difficult to be optimized quickly to satisfy the requirement of resonance, which is consistent with the conclusion in Ref. [30].

4. Dynamical Behavior of Fractional-Order FHN Neuron Model

4.1. Effect of the Fractional-Order

In this subsection, the dynamical behavior of the fractional-order FHN neuron without time-delay feedback is studied. Figure 4 shows the bifurcation diagram of interspike interval (ISI) of the neuron versus the fractional-order α , where α [ 0.6 , 1 ] . From Figure 4, when α is small, it can be seen that the fractional-order FHN neuron does not fire, and when the α is relatively large, various periodic and irregular firing patterns appear in the neuron. It can be found that the fractional-order FHN neuron displays complex dynamics by adjusting a single parameter α , which is also observed in other neurons [31,32].
In order to intuitively depict the effect of the fractional-order α on the dynamical behavior of the FHN neuron, the evolution of the membrane potential x ( t ) is given in Figure 5. From Figure 5a, when α = 0.6 , the neuron is in the quiescent state. As can be seen from Figure 5b,c, when α = 0.792 and α = 0.84 , the firing pattern of the neuron is period 1 bursting. Obviously, their interspike intervals are different, which are four times and two times 2 π / ω , respectively. In Figure 5a, when α = 0.934 , the discharge rhythm of the FHN neuron alternates between period 1 and period 2. From Figure 5a–d, it is found that the fractional-order FHN neuron is more active with the increase of α , which is consistent with that revealed in Figure 4.

4.2. Effect of the Time Delay

In order to study the effect of time-delay feedback on the dynamical behavior of the fractional-order FHN neuron, α = 0.84 in Figure 5c is selected. Figure 6 depicts the bifurcation diagram of interspike interval (ISI) of the neuron versus the time-delay parameter τ , from which it can be seen that the fractional-order FHN neuron with time-delay feedback exhibits rich dynamical behaviors. The time series of the membrane potential x ( t ) for different values of τ are given in Figure 7. When τ = 0 , the neuron fires with period 1, which corresponds to the case of the neuron without time delay in Figure 5c. When τ = 5.6 , it is seen that the firing pattern of the neuron is period 2 bursting. It can be seen from Figure 7c that a new firing pattern appears in the neuron, which regularly alters between period 2 and period 1. The multiple spiking of the firing pattern also appears in the neuron, as shown in Figure 7d. Thus, bursting patterns and interspike intervals of the fractional-order FHN neuron can be controlled by adjusting α or τ .

5. Conclusions

The vibrational resonance and electrical activity behaviors in the fractional-order FHN neuron are studied in this paper. When the original integer-order model is extended to the fractional-order model, the region of the VR will be wider, so the requirement for the amplitude of high-frequency signal can be relaxed, and then, the VR phenomenon can be induced by choosing the appropriate fractional-order α . Introducing time-delay feedback, it is found that the phenomenon of vibrational resonance appears with two different periods in the fractional FHN neuron, which are exactly equal to the periods of two frequency signals. The fractional-order FHN neuron exhibits rich electrical activities, and its bursting patterns and interspike interval can be controlled by adjusting α or τ .

Author Contributions

Conceptualization, J.-W.M. and D.-L.H.; methodology, J.-W.M. and D.-L.H.; validation, J.-W.M. and D.-L.H.; formal analysis, J.-W.M.; investigation, J.-W.M.; resources, D.-L.H.; data curation, J.-W.M.; writing—original draft preparation, J.-W.M.; writing—review and editing, D.-L.H.; visualization, D.-L.H.; supervision, D.-L.H.; project administration, D.-L.H.; funding acquisition, D.-L.H. All authors have read and agreed to the published version of the manuscript.

Funding

This work in the present paper was supported by the National Natural Science Foundation of China (Grant No. 11502067), and Natural Science Foundation of Jiangsu Province of China (Grant No. BK20191295).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Landa, P.S.; Mcclintock, P. Vibrational resonance. J. Phys. A 2000, 33, L433–L438. [Google Scholar] [CrossRef]
  2. Mcnamara, B.; Wiesenfeld, K. Theory of stochastic resonance. Phys. Rev. A 1989, 39, 4854–4869. [Google Scholar] [CrossRef] [PubMed]
  3. Maksimov, A.O. On the subharmonic emission of gas bubbles under two-frequency excitation. Ultrasonics 1997, 35, 79–86. [Google Scholar] [CrossRef]
  4. Volkov, E.I.; Ullner, E.; Zaikin, A.A.; Kurths, J. Oscillatory amplification of stochastic resonance in excitable systems. Phys. Rev. E 2003, 68, 026214. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Gherm, V.E.; Zernov, N.N.; Lundborg, B.; Västberg, A. The two-frequency coherence function for the fluctuating ionosphere: Narrowband pulse propagation. J. Atmos. Sol.-Terr. Phys. 1997, 59, 1831–1841. [Google Scholar] [CrossRef]
  6. Victor, J.D.; Conte, M.M. Two-frequency analysis of interactions elicited by Vernier stimuli. Vis. Neurosci. 2000, 17, 959–973. [Google Scholar] [CrossRef] [Green Version]
  7. Yao, C.; Liu, Y.; Zhan, M. Frequency-resonance-enhanced vibrational resonance in bistable systems. Phys. Rev. E 2011, 83, 061122. [Google Scholar] [CrossRef]
  8. Abusoua, A.; Daqaq, M.F. Experimental evidence of vibrational resonance in a mechanical bistable twin-well oscillator. J. Comput. Nonlinear Dyn. 2018, 3, 061002. [Google Scholar] [CrossRef]
  9. Chizhevsky, V.N. Vibrational higher-order resonances in an overdamped bistable system with biharmonic excitation. Phys. Rev. E 2014, 90, 042924. [Google Scholar] [CrossRef]
  10. Yang, J.; Sanjuan, M.; Liu, H. Bifurcation and resonance in a fractional Mathieu-Duffing oscillator. Eur. Phys. J. B 2015, 88, 310. [Google Scholar] [CrossRef]
  11. Yang, Z.; Ning, L. Vibrational resonance in a harmonically trapped potential system with time delay. Pramana 2019, 92, 1–12. [Google Scholar] [CrossRef]
  12. Uzuntarla, M.; Yilmaz, E.; Wagemakers, A.; Ozer, M. Vibrational resonance in a heterogeneous scale free network of neurons. Commun. Nonlinear Sci. Numer. Simul. 2015, 22, 367–374. [Google Scholar] [CrossRef]
  13. Meral, F.C.; Royston, T.J.; Magin, R. Fractional calculus in viscoelasticity: An experimental study. Commun. Nonlinear Sci. Numer. Simul. 2010, 15, 939–945. [Google Scholar] [CrossRef]
  14. Matouk, A.E. Chaos, feedback control and synchronization of a fractional-order modified Autonomous Van der Pol–Duffing circuit. Commun. Nonlinear Sci. Numer. Simul. 2011, 16, 975–986. [Google Scholar] [CrossRef]
  15. Chaudhary, N.I.; Raja, M.; Khan, A.R. Design of modified fractional adaptive strategies for Hammerstein nonlinear control autoregressive systems. Nonlinear Dyn. 2015, 82, 1811–1830. [Google Scholar] [CrossRef]
  16. Khanday, F.A.; Kant, N.A.; Dar, M.R.; Zulkifli, T.Z.A.; Psychalinos, C. Low-Voltage low-power integrable CMOS circuit implementation of integer and fractional-order FitzHugh-Nagumo Neuron model. IEEE Trans. Neural Netw. Learn. Syst. 2019, 30, 2108–2122. [Google Scholar] [CrossRef]
  17. Yang, J.H.; Zhu, H. Vibrational resonance in Duffing systems with fractional-order damping. Chaos 2012, 22, 013112. [Google Scholar] [CrossRef]
  18. Mbong, T.L.M.D.; Siewe, M.S.; Tchawoua, C. The effect of the fractional derivative order on vibrational resonance in a special fractional quintic oscillator. Mech. Res. Commun. 2016, 78, 13–19. [Google Scholar] [CrossRef]
  19. Qin, T.; Xie, T.; Luo, M.; Deng, K. Vibrational resonance in fractional-order overdamped multistable systems. Chin. J. Phys. 2017, 55, 546–555. [Google Scholar] [CrossRef]
  20. Lundstrom, B.N.; Higgs, M.H.; Spain, W.J.; Fairhall, A.L. Fractional differentiation by neocortical pyramidal neurons. Nat. Neurosci. 2008, 11, 1335–1342. [Google Scholar] [CrossRef]
  21. Anastasio, T.J. The fractional-order dynamics of brainstem vestibulo-oculomotor neurons. Biol. Cybern. 1994, 72, 69–79. [Google Scholar] [CrossRef]
  22. Ullner, E.; Zaikin, A.; García-Ojalvo, J.; Bascones, R.; Kurths, J. Vibrational resonance and vibrational propagation in excitable systems. Phys. Lett. A 2003, 312, 348–354. [Google Scholar] [CrossRef]
  23. Deng, B.; Wang, J.; Wei, X.; Tsang, K.M.; Chan, W.L. Vibrational resonance in neuron populations. Chaos 2010, 20, 013113. [Google Scholar] [CrossRef]
  24. Izhikevich, E.M. Simple model of spiking neurons. IEEE Trans. Neural Netw. 2003, 14, 1569–1572. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Hu, D.L.; Yang, J.H.; Liu, X.B. Delay-induced vibrational multiresonance in FitzHugh-Nagumo system. Commun. Nonlinear Sci. Numer. Simul. 2012, 17, 1031–1035. [Google Scholar] [CrossRef]
  26. Guo, W.; Ning, L. Vibrational resonance in a fractional order quintic oscillator system with time delay feedback. Int. J. Bifurc. Chaos 2020, 30, 2050025. [Google Scholar] [CrossRef]
  27. Abdelouahab, M.S.; Lozi, R.P.; Chen, G. Complex Canard Explosion in a Fractional-Order FitzHugh–Nagumo Model. Int. J. Bifurc. Chaos 2019, 29, 1950111. [Google Scholar] [CrossRef]
  28. Wu, D.; Zhu, S. Stochastic resonance in FitzHugh–Nagumo system with time-delayed feedback. Phys. Lett. A 2008, 372, 5299–5304. [Google Scholar] [CrossRef]
  29. Teka, W.; Marinov, T.M.; Santamaria, F. Neuronal spike timing adaptation described with a fractional leaky integrate-and-fire model. PLoS Comput. Biol. 2014, 10, e1003526. [Google Scholar] [CrossRef] [PubMed]
  30. Wu, C.; Yang, J.; Huang, D.; Liu, H.; Hu, E. Weak signal enhancement by the fractional-order system resonance and its application in bearing fault diagnosis. Meas. Sci. Technol. 2019, 30, 035004. [Google Scholar] [CrossRef]
  31. Teka, W.W.; Upadhyay, R.K.; Mondal, A. Spiking and bursting patterns of fractional-order Izhikevich model. Commun. Nonlinear Sci. Numer. Simul. 2018, 56, 161–176. [Google Scholar] [CrossRef]
  32. Mondal, A.; Sharma, S.K.; Upadhyay, R.K.; Mondal, A. Firing activities of a fractional-order FitzHugh-Rinzel bursting neuron model and its coupled dynamics. Sci. Rep. 2019, 9, 15721. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Response amplitude Q as a function of α and F . (a) The three-dimensional surface of response amplitude Q versus α and F ; (b) response amplitude Q versus F for different values of α .
Figure 1. Response amplitude Q as a function of α and F . (a) The three-dimensional surface of response amplitude Q versus α and F ; (b) response amplitude Q versus F for different values of α .
Mathematics 10 00087 g001
Figure 2. Relationship curve between response amplitude Q and fractional-order α . (a) Response amplitude Q max versus α ; (b) response amplitude Q versus α for three different values of F .
Figure 2. Relationship curve between response amplitude Q and fractional-order α . (a) Response amplitude Q max versus α ; (b) response amplitude Q versus α for three different values of F .
Mathematics 10 00087 g002
Figure 3. Relationship curve between response amplitude Q and τ : (a) ω = 0.5 ,   Ω = 5 ,   α = 0.96 ,   F = 0.123 ; (b) ω = 0.5 ,   Ω = 5 ,   α = 0.85 ,   F = 0.138 ; (c) ω = π / 4 ,   Ω = 12 ,   α = 0.7 ,   F = 0.395 ; (d) ω = π / 4 ,   Ω = 12 ,   α = 0.6 ,   F = 0.298 ; (e) ω = 0.5 ,   Ω = 5 ,   α = 0.96 ,   F = 0.123 ; (f) ω = π / 4 ,   Ω = 12 ,   α = 0.6 ,   F = 0.298 .
Figure 3. Relationship curve between response amplitude Q and τ : (a) ω = 0.5 ,   Ω = 5 ,   α = 0.96 ,   F = 0.123 ; (b) ω = 0.5 ,   Ω = 5 ,   α = 0.85 ,   F = 0.138 ; (c) ω = π / 4 ,   Ω = 12 ,   α = 0.7 ,   F = 0.395 ; (d) ω = π / 4 ,   Ω = 12 ,   α = 0.6 ,   F = 0.298 ; (e) ω = 0.5 ,   Ω = 5 ,   α = 0.96 ,   F = 0.123 ; (f) ω = π / 4 ,   Ω = 12 ,   α = 0.6 ,   F = 0.298 .
Mathematics 10 00087 g003
Figure 4. Bifurcation diagram of interspike interval (ISI) versus the fractional-order α ; other parameters are ω = 0.5 ,   Ω = 5 ,   F = 0.15 .
Figure 4. Bifurcation diagram of interspike interval (ISI) versus the fractional-order α ; other parameters are ω = 0.5 ,   Ω = 5 ,   F = 0.15 .
Mathematics 10 00087 g004
Figure 5. Time evolution of the membrane potential x ( t ) for different values of α ; other parameters are ω = 0.5 ,   Ω = 5 ,   F = 0.15 ; (a) α = 0.6 ; (b) α = 0.792 ; (c) α = 0.84 ; (d) α = 0.934 .
Figure 5. Time evolution of the membrane potential x ( t ) for different values of α ; other parameters are ω = 0.5 ,   Ω = 5 ,   F = 0.15 ; (a) α = 0.6 ; (b) α = 0.792 ; (c) α = 0.84 ; (d) α = 0.934 .
Mathematics 10 00087 g005aMathematics 10 00087 g005b
Figure 6. Bifurcation diagram of interspike interval (ISI) versus the time-delay τ ; other parameters are α = 0.84 ,   ω = 0.5 ,   Ω = 5 ,   F = 0.15 .
Figure 6. Bifurcation diagram of interspike interval (ISI) versus the time-delay τ ; other parameters are α = 0.84 ,   ω = 0.5 ,   Ω = 5 ,   F = 0.15 .
Mathematics 10 00087 g006
Figure 7. Time evolution of the membrane potential x ( t ) for different values of τ ; other parameters are α = 0.84 ,   ω = 0.5 ,   Ω = 5 ,   F = 0.15 ; (a) τ = 0 ; (b) τ = 5.6 ; (c) τ = 25 ; (d) τ = 31.8 .
Figure 7. Time evolution of the membrane potential x ( t ) for different values of τ ; other parameters are α = 0.84 ,   ω = 0.5 ,   Ω = 5 ,   F = 0.15 ; (a) τ = 0 ; (b) τ = 5.6 ; (c) τ = 25 ; (d) τ = 31.8 .
Mathematics 10 00087 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mao, J.-W.; Hu, D.-L. Vibrational Resonance and Electrical Activity Behavior of a Fractional-Order FitzHugh–Nagumo Neuron System. Mathematics 2022, 10, 87. https://doi.org/10.3390/math10010087

AMA Style

Mao J-W, Hu D-L. Vibrational Resonance and Electrical Activity Behavior of a Fractional-Order FitzHugh–Nagumo Neuron System. Mathematics. 2022; 10(1):87. https://doi.org/10.3390/math10010087

Chicago/Turabian Style

Mao, Jia-Wei, and Dong-Liang Hu. 2022. "Vibrational Resonance and Electrical Activity Behavior of a Fractional-Order FitzHugh–Nagumo Neuron System" Mathematics 10, no. 1: 87. https://doi.org/10.3390/math10010087

APA Style

Mao, J. -W., & Hu, D. -L. (2022). Vibrational Resonance and Electrical Activity Behavior of a Fractional-Order FitzHugh–Nagumo Neuron System. Mathematics, 10(1), 87. https://doi.org/10.3390/math10010087

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