
Citation: Darrell B. Roundy, April Hulet, Bruce A. Roundy, Ryan R. Jensen, Jordan B. Hinkle, Leann Crook, Steven L. Petersen. Estimating pinyon and juniper cover across Utah using NAIP imagery[J]. AIMS Environmental Science, 2016, 3(4): 765-777. doi: 10.3934/environsci.2016.4.765
[1] | Mahtab Mehrabbeik, Fatemeh Parastesh, Janarthanan Ramadoss, Karthikeyan Rajagopal, Hamidreza Namazi, Sajad Jafari . Synchronization and chimera states in the network of electrochemically coupled memristive Rulkov neuron maps. Mathematical Biosciences and Engineering, 2021, 18(6): 9394-9409. doi: 10.3934/mbe.2021462 |
[2] | Prasina Alexander, Fatemeh Parastesh, Ibrahim Ismael Hamarash, Anitha Karthikeyan, Sajad Jafari, Shaobo He . Effect of the electromagnetic induction on a modified memristive neural map model. Mathematical Biosciences and Engineering, 2023, 20(10): 17849-17865. doi: 10.3934/mbe.2023793 |
[3] | Andrey Olypher, Jean Vaillant . On the properties of input-to-output transformations in neuronal networks. Mathematical Biosciences and Engineering, 2016, 13(3): 579-596. doi: 10.3934/mbe.2016009 |
[4] | Pshtiwan Othman Mohammed, Hari Mohan Srivastava, Sarkhel Akbar Mahmood, Kamsing Nonlaopon, Khadijah M. Abualnaja, Y. S. Hamed . Positivity and monotonicity results for discrete fractional operators involving the exponential kernel. Mathematical Biosciences and Engineering, 2022, 19(5): 5120-5133. doi: 10.3934/mbe.2022239 |
[5] | Jian Huang, Zhongdi Cen, Aimin Xu . An efficient numerical method for a time-fractional telegraph equation. Mathematical Biosciences and Engineering, 2022, 19(5): 4672-4689. doi: 10.3934/mbe.2022217 |
[6] | Guowei Wang, Yan Fu . Spatiotemporal patterns and collective dynamics of bi-layer coupled Izhikevich neural networks with multi-area channels. Mathematical Biosciences and Engineering, 2023, 20(2): 3944-3969. doi: 10.3934/mbe.2023184 |
[7] | Bingrui Zhang, Jin-E Zhang . Fixed-deviation stabilization and synchronization for delayed fractional-order complex-valued neural networks. Mathematical Biosciences and Engineering, 2023, 20(6): 10244-10263. doi: 10.3934/mbe.2023449 |
[8] | Amit Kumar, Jehad Alzabut, Sudesh Kumari, Mamta Rani, Renu Chugh . Dynamical properties of a novel one dimensional chaotic map. Mathematical Biosciences and Engineering, 2022, 19(3): 2489-2505. doi: 10.3934/mbe.2022115 |
[9] | Hwayeon Ryu, Sue Ann Campbell . Stability, bifurcation and phase-locking of time-delayed excitatory-inhibitory neural networks. Mathematical Biosciences and Engineering, 2020, 17(6): 7931-7957. doi: 10.3934/mbe.2020403 |
[10] | Chenxi Huang, Qianqian Zhang, Sanyi Tang . Non-smooth dynamics of a SIR model with nonlinear state-dependent impulsive control. Mathematical Biosciences and Engineering, 2023, 20(10): 18861-18887. doi: 10.3934/mbe.2023835 |
Computational neuroscience studies the function of the brain from different molecular, cellular and behavioral levels [1]. The mathematical neuronal models are used to evaluate the qualitative function of neurons [2,3]. A wide variety of neural activities, such as different firing patterns, can be represented by neuronal models [4,5,6,7,8]. Neuronal models are generally divided into continuous time and discrete time. Continuous time neuronal models are described by the ordinary differential equations (ODE) such as Hodgkin-Huxley [9], Fitzhugh-Nagumo [10,11], and Hindmarsh-Rose [12,13]. For example, Njitack and colleagues examined the effect of a memristive autapse on the dynamics of a two-dimensional Hindmarsh–Rose neuron [14]. In another work, the effect of a two-dimensional FitzHugh–Nagumo, along with a 3D Hindmarsh–Rose neuronal model, was evaluated by coupling a multistability memristive synapse [15]. Models of continuous time neurons have been extensively studied dynamically, while less attention has been paid to the map-based discrete time models [16]. The distinctive features of discrete-time neuron maps are their high potential for complex behavior modeling, memory savings, and simplifying large neural network calculations [17,18]. Therefore, many scientists have tried to introduce discrete-time neuron maps [19]. For example, with the Euler discretization of the Hodgkin-Huxley neuronal model, the Izhikevich neuronal model was introduced [20]. The Nagumo-Sato discrete neuronal model was proposed in 1972 [21], and in 1990 a modified version of this model was introduced[22]. A one-dimensional map-based neuronal model of self-interaction was proposed in 1997 by Pasemann [23]. In 2001, Rulkov proposed a simple two-dimensional turbulent neuron model, one of the first phenomenological models to show a variety of basic neural behaviors (bursts of spikes, spikes, silence) [24]. Mathematical modeling of bursts and some spikes is impossible in the context of two-dimensional continuous-time systems, but Rulkov's two-dimensional model can demonstrate these behaviors [25].
The dynamic behaviors of the coupled Rulkov neuron map have been studied in recent research [26,27,28]. Sun et al. investigated the complete synchronization of Rulkov chaotic neural networks using master stability function analysis [29]. Franović investigated the effects of synaptic time delay on Rulkov map neurons with chemical coupling [30]. Rakshit et al. examined the stabilization, neural synchronization, and dynamics of two coupled Rulkov neurons with electrical and chemical synaptic interactions [31]. Wagemakers et al. introduced an electronic implementation of the Rulkov neuron map [32]. Li et al. proposed a discrete memristive neuron model based on the Rulkov model, which can better describe the actual firing activity of biological neurons due to biophysical memory [33]. Cheng et al. studied the synchronization dynamics of two heterogeneous chaotic Rulkov neurons with electrical synapses [34]. It is notable that the target for synchronization stability means energy balance between neurons and the neurons reach energy balance after complete synchronization [35,36].
Discrete fractional calculus has been used in the dynamic analysis of map-based systems. The fractional-order derivatives are relatively more accurate than actual order derivatives because they are an effective tool for describing the effect of memory on processing [37]. Memory effects in fractional order indicate that all previous models are involved in determining the system modes [38]. Over the past several years, researchers have focused on incorporating the fractions in various applications of neural networks, physics, biology, etc. [39,40,41,42,43]. The discrete fraction order equation is also a powerful tool in modeling chaotic systems [44,45,46]. For example, Guo-Cheng et al. reported the delayed logistic equation using the discrete fractional calculus approach and the corresponding discrete perturbation [47]. Elsonbaty et al. proposed a susceptible-infected-treatment-recovered-susceptible (SITRS) fractional discrete model to simulate the coronavirus (COVID-19) epidemic [48]. Kassim et al. developed a secure image transfer scheme based on synchronizing turbulent time-discrete fractional systems [49]. Fractional derivative properties can be expressed in the biological sense as follows [50]:
1. The fractional-order derivative represents the inheritance characteristics more accurately compared to integer-order models.
2. FO neuron models correctly represent biological properties in the presence of noise, which full-order models fail to do.
3. Dynamic firing activities of multiple time scales for a single neuron can be treated with FO derivatives.
4. Firing frequency responses are better in fractional designs than in integer models.
5. Dynamic excitatory properties reveal several neurocomputational properties that can be achieved in FO systems, which enrich functional neural mechanisms.
In general, human evolution is done by two genetic systems based on DNA and another based on the transmission of information through the functions of the nervous system. Mechanisms underlying neuron-based inheritance include hippocampal neurogenesis and memory and learning processes that generate new neural complexes and alter brain structure and function [51]. From the concept of neuroscience, the discrete fractional order models the neuron's hereditary behaviors as memory in the dynamic model. Given the above, this paper aims to assist researchers in computational neuroscience by introducing a discrete fraction order derivative to the Rulkov neuronal model. The dynamics of the fractional order map is analyzed completely and, in the end, its synchronization is examined. Synchronization of neurons has become an interesting topic in recent years [26,52,53] due to its critical applications in the neural processes [54,55]. This paper includes the following sections: Section 2 introduces the Rulkov neuron map, demonstrating the system's dynamic and biological behaviors. Section 3 presents a discrete fractional order Rulkov neuron map and dynamically evaluates the effect of fractional order and model parameters. In section 4, the fixed points are calculated, and their stability is examined. Section 5 investigates the synchronization of two fractional-order Rulkov neurons in different fractional orders, and finally, Section 6 gives conclusions.
In this section, Rulkov neuron map is introduced as a discrete system. Rulkov neuron map is a discrete biological neuronal model proposed by Nikolai F. Rulkov in 2001 [24,25]. This model, in which n indicates the discrete time, is defined as follows:
xn+1=α1+x2n+yn, | (2.1) |
yn+1=yn−μ(xn−σ), | (2.2) |
where the variable xn is the membrane potential of neurons, and the variable yn in the model is a system-slowing variable due to the minimal value of μ (0<μ<<1). Unlike the variable xn, the variable yn has no clear biological meaning [30]. The parameter σ is an external dc current to neurons, and the parameter α is a nonlinear map parameter. The value of μ is always considered 0.001. By changing the parameters σ and α of the system, different behaviors such as silence (σ=−2,α=6.22), bursting (σ=−1,α=4.1) and chaotic firing (σ=−1,α=5.7) can be seen which are shown in Figure 1 for random initial conditions. Phase diagram for burst mode and chaotic firing mode are shown in parts (d) and (e) of Figure 1.
In continuation, the dynamic analysis of the Rulkov neuron map is performed through bifurcation diagrams and the maximum Lyapunov exponents (MLE) calculations. The positive MLE proves the existence of chaotic behavior. The bifurcation diagram is plotted by changing the parameter α in the range of 4 to 20 with random initial conditions in the range of [−1,1] for the variable xn in part (a) and for the variable yn in part (b) of Figure 2. Also, the MLE of the model for parameter α is shown in part (c) of Figure 2 (σ=−0.1,μ=0.001). It can be seen that with the change of α in the interval [4,4.725], the neuron's firing is period one (MLE is zero) and then enters chaos in α=4.725 until α=8.576. At α=8.576, the firing changes to periodic and again by period doubling enters another chaotic region in α=14.688, and finally exits from chaos in α=19.280. There is also a periodic window in the interval [16.213,16.327]. The bifurcation diagram for variables xn and yn with respect to variation of parameter σ in the range of [−2.5,0] are shown in parts (d) and (e) of Figure 2 (α=4.1, and μ=0.001). The corresponding MLE is drawn in part (f) of Figure 2. By changing σ, the dynamic behavior is periodic for the interval [−2.5,−1.632], [−1.632,−0.145], and then switches to chaotic. The bifurcation diagrams by changing parameter μ in the interval [0,1] for α=4.1, and σ=−0.1 are shown in parts (g) and (h) of Figure 2 for variables xn and yn. By changing the μ parameter, the chaotic behavior is seen in the intervals [0,0.087] and [0.377,0.805].
In this section, Rulkov neuron map is reformulated as a fractional-order discrete system using the Caputo left difference operator. Discrete fraction order maps have recently been recognized as a new powerful tool in modeling and describing the dynamics of complex nonlinear discrete systems [56,57]. Different numerical methods have been developed to solve discrete fraction order maps [58,59]. In this paper, we examine the following Rulkov neuron map with discrete fraction order:
Δqax(t)=α1+x(t+q−1)2+y(t+q−1), | (3.1) |
Δqay(t)=y(t+q−1)−μ(x(t+q−1)−σ), | (3.2) |
t∈Na+1−q,x(0)=c1,y(0)=c2,0<q, |
where Δqa is a function of the Caputo delta difference on the left, Na=[a,a+1,a+2,…] represents the time scale isolated (a∈R and constant). To obtain the fractional order form of the Rulkov neuron map, we first define the necessary definitions from the discrete fractional calculus and generalize the results to the Rulkov neuron map.
Definition 1. When u:Na→R, q>0 and q∉N, the sum of the fraction of order q is defined as follows:
Δ−qau(t)≔1Γ(q)∑t−qs=a(t−s−1)(q−1)u(s),t∈Na+q, | (3.3) |
where a is the starting point, and t(q) is the falling function, which is defined as follows:
t(q)=Γ(t+1)Γ(t−q+1)=t(t−1)…(t−q+1), | (3.4) |
where Γ(t) is Gamma function, defined as follows [60]:
Γ(t)=∫∞0xq−1e−xdx | (3.5) |
Definition 2. For q>0, q∈Na and u(t) defined on Na, the difference in the Caputo delta fraction can be expressed as follows [58]:
Δqau(t):=Δ−(m−q)aΔmu(t)=1Γ(m−q)∑t−(m−q)s=a(t−s−1)(m−q−1)Δmu(s), | (3.6) |
where m=[q]+1 and m=[q]+1.
Theorem 1. For the delta fraction difference equation
Δqau(t)=f(t+q−1,u(t+q−1)),Δku(a)=uk,m=[q]+1,k=0,…,m−1, | (3.7) |
the equivalent discrete integral equation can be obtained as follows:
u(t)=u0(t)+1Γ(q)∑t−qs=a+m−q(t−s−1)(q−1)f(s+q−1,u(s+q−1)),t∈Na+m, | (3.8) |
where the initial iteration u0(t) is expressed as follows:
u0(t)=∑m−1k=0(t−a)(k)k!Δku(a). | (3.9) |
Detailed descriptions of the above results are given in reference [61].
Rulkov's neural map of discrete fraction order (Eq. 3.1, and Eq. 3.2) can be written as the following classical map:
Δx(n)=α1+x(n)2+y(n),x(0)=c1, | (3.10) |
Δy(n)=y(n)−μ(x(n)−σ),y(0)=c2. | (3.11) |
According to Theorem 1, we can write the discrete fractional for Rulkov's model in the following form:
x(t)=x(0)+1Γ(q)∑t−qs=1−q(t−s−1)(q−1)(α1+(x(s+q−1))2+y(x(s+q−1))), | (3.12) |
y(t)=y(0)+1Γ(q)∑t−qs=1−q(t−s−1)(q−1)(y(s+q−1)−μ(x(s+q−1)−σ)). | (3.13) |
In conclusion, the numerical formula of the n-dependent discrete Rulkov fractional neuronal model is defined as follows:
x(n)=x(0)+1Γ(q)∑ni=1Γ(n−i+q)Γ(n−i+1)(α1+x2i−1+yi−1), | (3.14) |
y(n)=y(0)+1Γ(q)∑ni=1Γ(n−i+q)Γ(n−i+1)(yi−1−μ(xi−1−σ)). | (3.15) |
The difference between the maps of the integers of Eq. 2.1, and Eq. 2.2 and the fractional order of Eq. 3.14, and Eq. 3.15 is a discrete core function and x(n) and y(n) refer to the past information x(0),…,x(n−1) and y(0),…,y(n−1). As a result, the current state depends on all past forms, which express the memory effects. We define the proposed model in Eq. 3.14, and Eq. 3.15 as the fractional-order of discrete Rulkov neuron map (FORNM). Also, in Eq. 3.14, and Eq. 3.15, the following equation is used:
Γ(n−i+q)Γ(n−i+1)=elnlnΓ(n−i+q)−Γ(n−i+1). | (3.16) |
By changing the parameters α, σ, and μ of the FORNM system, different behaviors such as silence (α=3.15, σ=−2, μ=0.2, q=0.001), bursts of spikes (α=3.5, σ=−2, μ=0.2, q=0.001), and chaotic firing (α=6, σ=−1, μ=0.3, q=0.001) can be seen which are shown in Figure 3. When the membrane potential reaches a threshold, the neuron fires and generates a signal that travels to other neurons, which in turn increases or decreases their potential in response to the signal. The neuron model that fires when the threshold is crossed is called a spiking neuron. Neuron spiking occurs when up-down oscillations accompany this spike train [62,63]. According to the stated definition, bursts of spikes occurs in part (b) of Figure 3, and this spiking behavior has been previously identical in [24,31]. The phase diagram for the bursts of spikes mode is shown in part (d) and for the chaotic firing in part Ⅰ.
In the following, the dynamical analysis of the FORNM is done by plotting the bifurcation diagrams versus the fractional order (q) and the system parameters at different intervals.
The dynamic behavior of the FORNM is firstly investigated by changing the parameters of the Rulkov neuron map. It was shown that there are chaotic, periodic, and fixed-point dynamics in the system for different values of parameters α, σ, and μ. Figure 4 shows the bifurcation diagrams wherein the left, middle, and right columns, the fractional order is set to q=0.001, q=0.01, and q=0.02, respectively. The bifurcation parameter in parts (a-c) is α in the interval [0,12] for σ=−1 and μ=0.3, and in parts (d-f) is σ in the interval [−5,0] for α=4 and μ=0.3, and in parts (g-i) is μ in the Interval [0.15,0.5] for α=4 and σ=−1. The initial conditions are selected randomly. It can be observed that the dynamics are strongly dependent on q. Generally, it can be concluded that by increasing q, the chaotic behavior is observed for higher bifurcation parameters, and the region of chaotic behavior is shrunken.
Another important factor in FORNM is the fractional order parameter, which can affect the overall dynamics. Figure 5 presents the bifurcation diagrams of the model according to the fractional order for different parameter values. The bifurcation diagram for α=1.9, σ=−1, μ=0.3 in part (a) of Figure 5 shows chaotic behavior for the range [0,0.0102] and periodic behavior for the range [0.0102,0.1]. In part (b), in which α=4.9, σ=−1, μ=0.3 the chaos is seen in the range of [0,0.0180]. The bifurcation diagram in terms of q for σ=−2, α=4, and μ=0.3 is shown in part (c) that has chaotic behavior for the interval [0,0.00389], while for σ=−0.833 in part (d), the chaotic behavior happens in [0,0.01928]. The bifurcation diagram for μ=0.33, α=4, and σ=−1 is shown in part (e) of Figure 5. In this case, the chaotic behavior is in the range [0,0.01809]. In part (f) where μ=0.45, two chaotic regions in [0,0.2928] and [0.1860,0.2] are observed. As shown in Figure 5, increasing the parameters' values increases the chaotic region.
In this section, we calculate the equilibrium points of the proposed FORNM of Eq. 3.14, and Eq. 3.15 and evaluate them in terms of stability. The equilibrium points are calculated by:
α1+x(n)2+y(n)=0, | (4.1) |
y(n)−μ(x(n)−σ)=0, | (4.2) |
i.e.
y(n)=−α1+x(n)2, | (4.3) |
y(n)=μ(x(n)−σ). | (4.4) |
Putting two equations together:
−α1+x(n)2=μ(x(n)−σ), | (4.5) |
i.e.
μ(x(n)−σ)(1+x(n)2)+α=0, | (4.6) |
i.e.
x(n)3−σx(n)2+x(n)−σ+αμ=0. | (4.7) |
According to the method of solving the equation of degree 3 [64], by solving Eq. 4.7, we obtain a=−σ, b=1, c=−σ+αμ, so p=1−σ23, and q=αμ−2σ3−2σ327. The eigenvalues of Eq. 4.7 are equal to Δ=(2σ327+2σ3−αμ)24−(σ23−1)327, so the equilibrium points of the model in terms of Δ states are equal to:
● If Δ>0,
x=σ3+(σ3−α2μ+√(2σ327+2σ3−αμ)24−(σ23−1)327+σ327)13+(σ3−α2μ−√(2σ327+2σ3−αμ)24−(σ23−1)327+σ327)13, | (4.8) |
i.e.
y=μ((σ3−α2μ+√(2σ327+2σ3−αμ)24−(σ23−1)327+σ327)13−2σ3+(σ3−α2μ−√(2σ327+2σ3−αμ)24−(σ23−1)327+σ327)13), | (4.9) |
● If Δ=0,
x1=σ3−2(α2μ−σ3−σ327)13, | (4.10) |
x2=x3=σ3+(α2μ−σ3−σ327)13, | (4.11) |
i.e.
y1=−μ(2σ3+(α2μ−σ3−σ327)13), | (4.12) |
y2=x3=−μ(2σ3−(α2μ−σ3−σ327)13), | (4.13) |
● If Δ<0,
x1=σ3−2√3sinsin((3√3(2σ327+2σ3−αμ)2(σ23−1)32)3)√(σ23−1)3, | (4.14) |
x2=σ3−2√3sinsin(π3−(3√3(2σ327+2σ3−αμ)2(σ23−1)32)3)√(σ23−1)3, | (4.15) |
x3=σ3−2√3sinsin(π6−(3√3(2σ327+2σ3−αμ)2(σ23−1)32)3)√(σ23−1)3, | (4.16) |
i.e.
y1=−μ(2σ3+2√3sinsin((3√3(2σ327+2σ3−αμ)2(σ23−1)32)3)√(σ23−1)3), | (4.17) |
y2=−μ(2σ3+2√3sinsin(π3−(3√3(2σ327+2σ3−αμ)2(σ23−1)32)3)√(σ23−1)3), | (4.18) |
y3=−μ(2σ3+2√3sinsin(π6−(3√3(2σ327+2σ3−αμ)2(σ23−1)32)3)√(σ23−1)3). | (4.19) |
Thus, the equilibrium points of the system are calculated, and we use the following theorem to examine the stability conditions of the equilibrium points of the fractional order system.
Theorem 2. The zero equilibrium point of the system:
ΔqaF(t)=MF(t+q−1), | (4.20) |
is asymptotically stable when the eigenvalues λi,i=1,2,…,n of the matrix M have the following two conditions [65]:
|Arg(λi)|>qπ2, | (4.21) |
|λi|<(2coscos[|Arg(λi)|−π2−q])q. | (4.22) |
Therefore, we calculate the Jacobian matrix of Eq. 3.1, and Eq. 3.2, which is as follows:
J=[2αx(n)(x(n)2+1)21−μ1]. | (4.23) |
The eigenvalues of the Jacobin matrix are calculated for equilibrium points at different parameters, and the conditions of Theorem 2 are obtained numerically for the eigenvalues. The parameters that satisfy the Theorem 2 requirements make Eq. 3.1, and Eq. 3.2 always stable.
The stability regions of the equilibrium points of the FORNM model have been investigated by varying the parameters in the model. The fractional order in the FORNM model is examined in the first step. Stability regions of the FORNM model are shown in Figure 6 in the three-dimensional space of α−σ−μ parameters and the two-dimensional space of parameters α−σ, α−μ, and σ−μ. The fractional order in parts (a-c) are q=0.001, q=0.01, q=0.025 respectively. It is shown that the stable regions decrease by increasing the value of fractional order.
In the second step, the effects of q−σ−μ parameters with α=15 constant in part (a), the effects of α−q−μ parameters with σ=−5 constant in part (b), and the effects of α−σ−q parameters with constant μ=0.3 in part (c) on the stability of the FORNM model are investigated and shown in Figure 7. In part (a) of Figure 7, the stability of the system is lost as the value of σ increases toward zero, and the stability of the system is lost as μ increases toward one. In part (b) of Figure 7, the stability of the system increases with the increase of α value, and the stability of the system decreases with the increase of μ. In part (c) of Figure 7, by increasing the value of α and σ, stability is created in the system.
Another important aspect of the dynamical analysis is the synchronization in a chaotic state. Two coupled fractional-order Rulkov neuron maps can be formulated as follows:
x1(n)=x1(0)+1Γ(q)∑ni=1Γ(n−i+q)Γ(n−i+1)(α1+x1(i−1)2+y1(i−1))+dx(x2(n)−x1(n)), | (4.24) |
y1(n)=y1(0)+1Γ(q)∑ni=1Γ(n−i+q)Γ(n−i+1)(y1(i−1)−μ(x1(i−1)−σ))+dy(y2(n)−y1(n)), | (4.25) |
x2(n)=x2(0)+1Γ(q)∑ni=1Γ(n−i+q)Γ(n−i+1)(α1+x2(i−1)2+y2(i−1))+dx(x1(n)−x2(n)), | (4.26) |
y2(n)=y2(0)+1Γ(q)∑ni=1Γ(n−i+q)Γ(n−i+1)(y2(i−1)−μ(x2(i−1)−σ))+dy(y1(n)−y2(n)), | (4.27) |
where subscript 1 corresponds to neuron one and subscript two corresponds to neuron 2, dx is the coupling strength of variable x(n), and dy is the coupling strength of variable y(n). For both neuron models, the initial conditions are randomly chosen between -1 and 1. To determine the synchronization level of two neurons, the average error is calculated in the form of the following equation:
Error=√1N∑Nj=1(x1(n)−x2(n))2+(y1(n)−y2(n))2, | (4.28) |
where N represents the number of time series data samples. We investigate the effects of fractional order changes on the synchronization of neurons. Figure 8 shows the synchronization error by changing the coupling strengths dx and dy in different fractional orders. It has been proven that two integer-order Rulkov maps cannot reach synchronization when connected by electrical synapses [66]. Here, it can be seen that the synchronization cannot occur in the fractional-order neurons too. It should be noted that the systems become unstable for coupling coefficients higher than 1. Unsynchronized neuron maps based on previous papers show that increased activity in certain areas is accompanied by decreased activity in other areas [67], and this pattern exists during rest and sleep [68].
This paper proposed the discrete fractional order Rulkov neuron map in Caputo's concept. The dynamic analysis of the fractional model was performed, and its synchronization was evaluated. Adding discrete fraction calculus to the Rulkov neuron map can consider the effect of memory on neural model dynamics. Memory effects refer to the fact that system states are determined in fractional order by all previous states. At first, the Rulkov neuron map was examined in terms of phase plane, bifurcation diagram, and Lyapunov exponent. It was shown that the biological behaviors of the Rulkov neuron map, such as silence, bursting, and chaotic firing, also exist in the discrete fractional order Rulkov neuron map. Hence, the Rulkov neuron map of discrete fractional order shows the same behaviors as the Rulkov map of integer order. However, changing the fractional order can result in bifurcation. In the next step, the stability regions of the system from a theoretical and numerical point of view were examined and it was shown that increasing the order of the fraction decreases the stability in the FORNM system. In the last stage, the synchronization behavior of two FORNM systems was investigated. The results showed that the fractional-order systems do not achieve synchronization with an increase in the coupling coefficient, which is similar to the integer-order systems. Table 1 shows the reason for the superiority of the result of Rulkov fractional neuron map to the previous works.
Ref | Year | Model | Description |
[69] | 2021 | Nagumo–Sato discrete neuron map | A fractional-order version of the one-dimensional neuron map is proposed. The dynamic behavior of the model has been analyzed by drawing bifurcation diagrams and the Lyapunov power diagram. The appearance of the spiral wave in the two-dimensional network for the fractional order version has been investigated in terms of the effect of various parameters such as the amplitude and frequency of stimuli, the coupling strength, and the fractional order parameter. |
[70] | 2021 | Rulkov Neuron Map | A theoretical investigation is done on a fractional order version of the Rulkov neural model, which is only the time series analysis of the proposed model, and its dynamic behavior is not evaluated. |
[71] | 2022 | Rulkov Neuron Map | In this research, the model has been implemented by applying discrete memristor on two-dimensional Rulkov neuron map. The dynamic behaviors of the discrete memristor-based neuron have been analyzed by experiments including phase diagram, bifurcation, and spectral entropy complexity algorithm. In this system, the intended memory, which is the inheritance of information transmission, is not considered. |
[72] | 2020 | Rulkov Neuron Map | A theoretical investigation for the asymptotic stability and instability of two-dimensional independent linearly incommensurate systems of fractional order Caputo difference equations has been carried out. In this article, only the time series analysis of the proposed model is done and its dynamic behavior is not evaluated. |
- | - | Proposed Model | In this research, the fractional order discrete Rulkov neuron map is analyzed in terms of dynamic behavior and synchronization. Biological behaviors of the Rulkov neuron map, such as silence, bursting, and chaotic firing, are also present in its discrete fractional order version. The bifurcation diagrams of the proposed model are investigated under the influence of neuronal model parameters and fractional order. The stability regions of the system are obtained theoretically and numerically, and it is shown that increasing the order of the fractional order causes the reduction of the stable regions. In general, the hereditary behavior of the neural model is proposed using discrete fractional order. |
For future work, discrete fractional order derivatives can be implemented on DNA-based genetic systems to investigate hereditary factors, so that human evolution can be described from these systems. For the next research topic, different types of discrete fractional order derivatives can be evaluated on the Rulkov neuron map. Due to the fact that the Rulkov discrete neuron map has no synchronization power, different controllers can be used to check synchronization.
This work is partially funded by Centre for Nonlinear Systems, Chennai Institute of Technology, India, vide funding number CIT/CNS/2023/RP/005.
The authors declare that they have no conflict of interest.
[1] | Miller RF, Tausch RJ (2001) The role of fire in pinyon and juniper woodlands: a descriptive analysis. In: K. E. M. Galley and T. P. Wilson [eds]. Proceedings of the Invasive Species Workshop: the Role of Fire in the Control and Spread of Invasive Species. Fire Conference 2000: the First National Congress on Fire Ecology, Prevention, and Management. Miscellaneous Publication No. 11, Tall Timber Research Station, Tallahassee, FL, USA. p. 15-30. |
[2] | Miller RF, Bates JD, Svejcar TJ, et al. (2005) Biology, ecology, and management of western juniper. Corvallis, OR, USA: Oregon State University, Technical Bulletin 152. 77 p. |
[3] | Miller RF, Svejcar TJ, Rose JA (2000) Impacts of western juniper on plant community composition and structure. J Range Manage 53: 574-585. |
[4] | Bates JD, Miller RF, Svejcar TJ (2005) Long-term successional trends following western juniper cutting. Rangeland Ecol Manag 58: 533-541. |
[5] | Pierson FB, Williams CJ, Kormos PR, et al. (2010) Hydrologic vulnerability of sagebrush steppe following pinyon and juniper encroachment. Rangeland Ecol Manag 63: 614-629. |
[6] |
Roundy BA, Miller RF, Tausch RJ, et al. (2014) Understory cover responses to piñon-juniper treatments across tree cover gradients in the Great Basin. Rangeland Ecol Manag 67: 482-494. doi: 10.2111/REM-D-13-00018.1
![]() |
[7] |
Knick ST, Hanser SE, Leu M (2014) Ecological Scale of Bird Community Response to Pinyon-Juniper Removal. Rangeland Ecol Manag 67: 553-562. doi: 10.2111/REM-D-13-00023.1
![]() |
[8] | Gruell GE (1999) Historical and modern roles of fire in pinyon-juniper. In: Steve B. Monsen and Richard Stevens [EDS]. Proceedings: Ecology and management of pinyon-juniper communities within the Interior West; 15-18 September 1997; Brigham Young University, Provo, Utah, USA. Ogden, UT, USA: U.S. Department of Agriculture, Forest Service, Rocky Mountain Research Station, RMRS-P-9. p. 24-28. |
[9] |
Miller RF, Ratchford J, Roundy BA, et al. (2014) Response of conifer-encroached shrublands in the Great Basin to prescribed fire and mechanical treatments. Rangeland Ecol Manag 67: 468-481. doi: 10.2111/REM-D-13-00003.1
![]() |
[10] | Young KR, Roundy BA, Eggett DL (2014) Mechanical mastication of Utah juniper encroaching sagebrush steppe increases inorganic soil N. Appl Environ Soil Sci, 2014. |
[11] | Young KR, Roundy BA, Bunting SC, et al. (2015) Utah juniper and two- needle piñon reduction alters fuel loads. Int J Wildland Fire 24: 236-248. |
[12] | Hulet A, Roundy BA, Petersen SL, et al. (2013) Assessing the relationship between ground measurements and object-based image analysis of land cover classes in pinyon and juniper woodlands. Photogramm Eng Rem S 79: 799-808. |
[13] | Hulet A, Roundy BA, Petersen SL, et al. (2013) An object-based image analysis of pinyon and juniper woodlands treated to reduce fuels. Environ Manag 53: 660-671. |
[14] | Hulet A, Roundy BA, Petersen SL, et al. (2014) Cover estimations using object-based image analysis rule sets developed across multiple scales in pinyon and juniper woodlands. Rangeland Ecol Manag 67: 318-327. |
[15] | Madsen MD, Zvirzdin DL, Davis BD, et al. (2011) Feature extraction techniques for measuring pińon and juniper tree cover and density, and comparison with field-based management surveys. Environ Manag 47: 766-776. |
[16] | Hulet A, Roundy BA, Petersen SL, et al. (2014) Utilizing NAIP imagery to estimate tree cover and biomass of pinyon and juniper woodlands. Rangeland Ecol Manag 67: 563-572. |
[17] | Weisberg PJ, Lingua E, Pillai RB (2007) Spatial patterns on pinyon-juniper woodland expansion in Central Nevada. Rangeland Ecol Manag 60: 115-124. |
[18] | Bybee J, Roundy BA, Young KR, et al. (2016) Vegetation response to pinon and juniper tree shredding. Rangeland Ecol Manag 69: 224-234. |
[19] | Booth DT, Cox SE, Fifield C, et al. (2005) Image analysis compared with other methods for measuring ground cover. Arid Land Res Manag 19: 91-100. |
[20] | Ko D, Bristow N, Greenwood D, et al. (2009) Canopy cover estimation in semiarid woodlands: comparison of field-based and remote sensing methods. Forest Sci 55: 132-141. |
[21] |
Laliberte AS, Fredrickson EL, Rango A (2007) Combining decision trees with hierarchical object-oriented image analysis for mapping arid rangelands. Photogramm Eng Rem S 73: 197-207. doi: 10.14358/PERS.73.2.197
![]() |
[22] |
Laliberte AS, Rango A, Herrick JE, et al. (2007) An object-based image analysis approach for determining fractional cover of senescent and green vegetation with digital plot photography. J Arid Environ 69: 1-14. doi: 10.1016/j.jaridenv.2006.08.016
![]() |
[23] |
Laliberte AS, Rango A (2009) Texture and scale in object-based analysis of subdecimeter resolution unmanned aerial vehicle (UAV) imagery. IEEE T Geosci Remote 47: 761-770. doi: 10.1109/TGRS.2008.2009355
![]() |
[24] | Davies KW, Petersen SL, Johnson DD, et al. (2010) Estimating juniper cover from national agriculture imagery program (NAIP) imagery and evaluating relationships between potential cover and environmental variables. Rangeland Ecol Manag 63: 630-637. |
[25] | Visual Information Solutions (2008) ENVI Feature Extraction Module User’s Guide. Available from: http://www.exelisvis.com/portals/0/pdfs/envi/feature_extraction_module.pdf. |
[26] | Blundell S, Opitz D, Morris M, et al. (2008) Feature Analyst V5.0. Paper presented at The ASPRS Annual Conference, Portland, Oregon, USA. |
[27] | Trimble (2011) eCognition Developer 8.64.1 reference book, version 8.64.1, Trimble Germany GmbH, München, Germany. |
[28] | McIver JD, Brunson M (2014) Multidisciplinary, multisite evaluation of alternative sagebrush steppe restoration treatments: the SageSTEP project. Rangeland Ecol Manag 67: 435-439. |
[29] | Bybee J, Roundy BA, Young KR, Hulet A, Roundy DB, Crook L, Aanderud Z, Eggett D, Cline NL (2016) Vegetation response to pinon and juniper tree shredding. Rangeland Ecol Manag 69: 224-234. |
[30] | Baatz M, Schäpe A (2000) Multiresolution segmentation—an optimization approach for high quality multi-scale image segmentation. In: J. Strobl, et al. [eds]. Angewandte Geographische Informationsverarbeitung XII. Wichmann, Heidelberg, Germany. p. 12-23. |
[31] | Congalton RG (2001) Accuracy assessment and validation of remotely sensed and other spatial information. Int J Wildland Fire 10: 321-328. |
[32] | Landis J, Koch G (1977) The measurement of observer agreement for categorical data. Biometrics 33: 159-174. |
Ref | Year | Model | Description |
[69] | 2021 | Nagumo–Sato discrete neuron map | A fractional-order version of the one-dimensional neuron map is proposed. The dynamic behavior of the model has been analyzed by drawing bifurcation diagrams and the Lyapunov power diagram. The appearance of the spiral wave in the two-dimensional network for the fractional order version has been investigated in terms of the effect of various parameters such as the amplitude and frequency of stimuli, the coupling strength, and the fractional order parameter. |
[70] | 2021 | Rulkov Neuron Map | A theoretical investigation is done on a fractional order version of the Rulkov neural model, which is only the time series analysis of the proposed model, and its dynamic behavior is not evaluated. |
[71] | 2022 | Rulkov Neuron Map | In this research, the model has been implemented by applying discrete memristor on two-dimensional Rulkov neuron map. The dynamic behaviors of the discrete memristor-based neuron have been analyzed by experiments including phase diagram, bifurcation, and spectral entropy complexity algorithm. In this system, the intended memory, which is the inheritance of information transmission, is not considered. |
[72] | 2020 | Rulkov Neuron Map | A theoretical investigation for the asymptotic stability and instability of two-dimensional independent linearly incommensurate systems of fractional order Caputo difference equations has been carried out. In this article, only the time series analysis of the proposed model is done and its dynamic behavior is not evaluated. |
- | - | Proposed Model | In this research, the fractional order discrete Rulkov neuron map is analyzed in terms of dynamic behavior and synchronization. Biological behaviors of the Rulkov neuron map, such as silence, bursting, and chaotic firing, are also present in its discrete fractional order version. The bifurcation diagrams of the proposed model are investigated under the influence of neuronal model parameters and fractional order. The stability regions of the system are obtained theoretically and numerically, and it is shown that increasing the order of the fractional order causes the reduction of the stable regions. In general, the hereditary behavior of the neural model is proposed using discrete fractional order. |
Ref | Year | Model | Description |
[69] | 2021 | Nagumo–Sato discrete neuron map | A fractional-order version of the one-dimensional neuron map is proposed. The dynamic behavior of the model has been analyzed by drawing bifurcation diagrams and the Lyapunov power diagram. The appearance of the spiral wave in the two-dimensional network for the fractional order version has been investigated in terms of the effect of various parameters such as the amplitude and frequency of stimuli, the coupling strength, and the fractional order parameter. |
[70] | 2021 | Rulkov Neuron Map | A theoretical investigation is done on a fractional order version of the Rulkov neural model, which is only the time series analysis of the proposed model, and its dynamic behavior is not evaluated. |
[71] | 2022 | Rulkov Neuron Map | In this research, the model has been implemented by applying discrete memristor on two-dimensional Rulkov neuron map. The dynamic behaviors of the discrete memristor-based neuron have been analyzed by experiments including phase diagram, bifurcation, and spectral entropy complexity algorithm. In this system, the intended memory, which is the inheritance of information transmission, is not considered. |
[72] | 2020 | Rulkov Neuron Map | A theoretical investigation for the asymptotic stability and instability of two-dimensional independent linearly incommensurate systems of fractional order Caputo difference equations has been carried out. In this article, only the time series analysis of the proposed model is done and its dynamic behavior is not evaluated. |
- | - | Proposed Model | In this research, the fractional order discrete Rulkov neuron map is analyzed in terms of dynamic behavior and synchronization. Biological behaviors of the Rulkov neuron map, such as silence, bursting, and chaotic firing, are also present in its discrete fractional order version. The bifurcation diagrams of the proposed model are investigated under the influence of neuronal model parameters and fractional order. The stability regions of the system are obtained theoretically and numerically, and it is shown that increasing the order of the fractional order causes the reduction of the stable regions. In general, the hereditary behavior of the neural model is proposed using discrete fractional order. |