A new approach on the modelling, chaos control and synchronization of a fractional biological oscillator

This research aims to discuss and control the chaotic behaviour of an autonomous fractional biological oscillator. Indeed, the concept of fractional calculus is used to include memory in the modelling formulation. In addition, we take into account a new auxiliary parameter in order to keep away from dimensional mismatching. Further, we explore the chaotic attractors of the considered model through its corresponding phase-portraits. Additionally, the stability and equilibrium point of the system are studied and investigated. Next, we design a feedback control scheme for the purpose of chaos control and stabilization. Afterwards, we introduce an efficient active control method to achieve synchronization between two chaotic fractional biological oscillators. The efficiency of the proposed stabilizing and synchronizing controllers is verified via theoretical analysis as well as simulations and numerical experiments.


Introduction
Chaos is one of the most prominent features of complex dynamical systems whose state variables are highly dependent on their initial conditions. This property may lead to the divergence of nonlinear chaotic systems; therefore, the analysis of such phenomena is of great importance from both theoretical and practical points of view [1].
Because of the wide appearance of chaotic systems in different fields of study such as secure communications, physics, biology and economy, a growing interest has been devoted to the control and synchronization of chaotic processes [2]. In [3], a new systematic approach within the use of scalar transmitted signal was designed to synchronize a class of hyperchaotic systems. In [4], a synchronization approach was presented for chaotic/hyperchaotic systems by formulating a linear feedback control problem. In [5], two identical hyperchaotic systems were synchronized by using a nonlinear control technique. In [6], a parameter observer was designed in order to identify unknown parameters in hyperchaotic systems, which was needed for designing a state-feedback controller. The author in [7] controlled and synchronized chaotic systems using state-dependent Riccati equations. In [8], a newly introduced fractional-order hyperchaotic system was controlled and synchronized by choosing an appropriate Lyapunov function and using a state-feedback control computation. In [9], Lorenz and Chen systems were considered in the presence of parameter uncertainty and controlled by using an adaptive dual synchronization controller. In [10], an adaptive terminal sliding mode control was applied to the chaotic Chen and hyperchaotic Lorenz systems. In [11], considering the stability theory of Lyapunov, the authors synchronized three nonidentical systems of different dimensions.
Recently, the investigation of chaos in biology has attracted the attention of many scientists due to the unpredictable behaviour of biological processes. In [12], the authors controlled the chaotic behaviour of tumor cells using a non-feedback loop. In [13], the author designed an optimal drug delivery schedule for the chaotic behaviour of tumor cells by minimizing Hamiltonian function; then he controlled the same system considering uncertainty in the chaotic model. In [14], the author analysed the chaotic behaviour of Lotka-Voltera biological systems. In [15], the authors synchronized two fractional chaotic maps by using a feedback control method. In [16], a fuzzy adaptive controller was designed to achieve a projective synchronization for a class of fractional chaotic systems with input nonlinearities.
Nowadays, a noticeable number of sensational researches have considered fractionalorder systems due to their memory-oriented features, a fact which makes such systems more realistic to describe real-world dynamics [17][18][19][20][21][22][23]. Thus, it is very practical to investigate and analyse chaos in fractional-order systems in order to simulate the complex behaviour of realistic phenomena. In addition, efficient stabilization and synchronization techniques are needed to be developed in order to overcome the chaotic behaviour of such fractional processes. Thus, the problem of chaos control and synchronization for complex fractional-order systems has been of great importance to a wide audience. More to the point, efficient numerical methods should be extended to deal with controversial problems in fractional calculus like modelling and control [24][25][26]. Inspired by the aforementioned argument, this paper introduces a new fractional model for a chaotic biological oscillator. We study the stability and equilibrium of the new system and analyse its chaotic behaviour by using phase-portrait responses. Then we propose a state-feedback control to stabilize the fractional model and apply an active control strategy to synchronize two identical chaotic oscillators. The applicability of the proposed controllers is also verified through theoretical analysis as well as some simulation results.
The reminder of this paper is organized as follows. In Sect. 2, some preliminary results are presented. In Sect. 3, the new mathematical fractional model of a chaotic biological oscillator is introduced, and its stability and equilibrium point are investigated. Afterwards, a state-feedback control is designed to stabilize the new model, and an active controller is introduced in order to achieve synchronization. Finally, the paper is closed by some concluding remarks in the last section.

Preliminaries
In this section, we recall the basic definition of Caputo fractional derivative and introduce its corresponding integral operator. The Laplace transform and the antiderivative property of this definition are also presented in the sequel.

Definition 2.1 ([27]) Let q ∈ (0, 1) and f be an integrable function. Then the qth-order
Caputo fractional derivative of f is defined by The associated qth-order Riemann-Liouville fractional integral of f is also defined as whenever the integral exists.
For the qth-order Caputo derivative when q ∈ (0, 1), the Laplace transform is introduced by where . In addition, the fractional derivative (1), together with the integral operator (2), satisfies the antiderivative property such that

The new fractional model
In the recent study [28], the authors employed an autonomous model to describe the interaction between enzyme and substrate mechanism. In [28], the aforesaid model was introduced in the form of the following system of ordinary differential equations: where the coefficients α, β, γ , ε, η, κ are real constants. As mentioned in [28], the model (5) exhibits chaotic attractors when the parameters (α, β, γ , ε, η, κ) take the values (5, 24, -0.05, 2.001, 2.55, 1.7), respectively. However, the integer-order model (5) suffers from the lack of memory effect, while hereditary property is the intrinsic feature of many complex biological systems. To overcome this drawback, the concept of fractional calculus is used to include memory in the model formulation. To do so, we extend the model (5) by substituting the integer-order derivatives with fractional Caputo ones defined in (1). Moreover, we take into account the new auxiliary parameter θ in order to keep away from dimensional mismatching [29]. Thus, we formulate the new fractional biological oscillator by Meanwhile, please note that the model (6) is the extension of a generalized Van der Pol oscillator; hence, this model can describe various biological and bio-physical phenomena.
Now, we list two well-known lemmas discussing the stability of fractional-order systems.

Lemma 3.1 ([30])
The equilibrium point f * of the fractional-order system (6) is locally asymptotically stable if the following inequality holds: where q is the fractional order, and λ i is the ith eigenvalue of J * at the point f * .
Then the equilibrium f * is locally asymptotically stable for all q ∈ (0, 1) if the corresponding Routh-Hurwitz table has elements with the same sign in its first column.

Simulation results
In the following, we simulate the complex behaviour of the new modified model (6)   Remark 3.1 Concerning the numerical simulation, here we applied the predictorcorrector method [32] in which an FFT algorithm is employed to evaluate the discrete convolutions [33]. The accuracy and the convergence of this scheme were studied in [34].
Also, contrary to the classical implementation in which the computational cost is proportional to N 2 , here we have N log(N 2 ) for the computational cost when the solution is evaluated at N time-points.

Figure 2
The phase planes of the fractional chaotic system (6) for q = 0.92

Chaos control and synchronization
In this section, designing a state-feedback controller, we compensate the undesirable chaotic behaviour of the fractional-order biological system (6). Then we introduce an active control approach to synchronize two identical fractional chaotic oscillators.

Chaos control
In order to control the chaotic behaviour of the fractional-order model (6), we consider the following modified controlled system: where u i ≥ 0 is the state-feedback control gain. The objective is to design u i in such a way that the state variables f 1 , f 2 , f 3 , f 4 of the controlled system (14) converge to the steadystate f * = (0, 0, 0, 0) asymptotically. Considering the system parameters (α, β, γ , ε, η, κ, θ ) = (5, 24, -0.05, 2.001, 2.55, 1.7, 0.99), we compute the Jacobian matrix J * at the point Then we derive

Synchronization
In this section, two identical fractional-order chaotic systems are synchronized by an active control scheme. To do so, the master and the slave systems are, respectively, consid-

Figure 6
The phase planes of the fractional chaotic system (6) for q = 1  Figure 7 The time-domain responses of the controlled system (14) for q = 0.9

Figure 8
The time-domain responses of the controlled system (14) for q = 0.92

Figure 9
The time-domain responses of the controlled system (14) for q = 0.94

Figure 10
The time-domain responses of the controlled system (14) for q = 0.96

Figure 11
The time-domain responses of the controlled system (14) for q = 0.98

Figure 12
The time-domain responses of the controlled system (14) for q = 1 where f i and g i are the states of the master and the slave systems, respectively, and u i is the control input. The objective is to design u i in such a way that the synchronization error goes to zero asymptotically. To reach this aim, we define the error of synchronization by e i (t) = g i (t)f i (t). Then the error dynamics will be Considering the control variables as one attains where v i represents the active control input. Now, we select where A is a constant matrix and selected in order to satisfy the stability condition in Lemma 3.1 for the error dynamics (23). Note that different choices exist for A; an appropriate one can be

Figure 13
The state trajectories of the synchronized systems (19) and (20) for q = 0.9

Figure 14
The state trajectories of the synchronized systems (19) and (20)  The state trajectories of the synchronized systems (19) and (20) for q = 0.94

Figure 16
The state trajectories of the synchronized systems (19) and (20) for q = 0.96

Figure 17
The state trajectories of the synchronized systems (19) and (20) for q = 0.98

Figure 18
The state trajectories of the synchronized systems (19) and (20) for q = 1 In addition, synchronization speed can be controlled by changing the main diagonal of matrix A.
Remark 4.1 Synchronization control can be regarded as a special tracking problem. When the master system is fixed in advance, the future dynamics of the master system are known, which can be considered as preview information. In this case, the speed of synchronization can be improved by using the preview control technique discussed in the new study [35].
To verify the theoretical synchronization approach mentioned above, we consider the initial states as

Conclusions
The objective of this paper was to investigate the chaotic behaviour of a fractional biological system and compensate this undesirable behaviour by an efficient control strategy.
In the first step, we discussed and explored the stability and the equilibrium of the considered model. Then we sketched the associated chaotic attractors in Figs. 1-6. Next, we designed a stabilizing state-feedback controller to overcome the undesirable chaotic behaviour. Afterwards, we introduced an active control design to synchronize two identical fractional biological oscillators. The validity of the chaos control and the synchronization approach was verified theoretically, and this analysis was also confirmed by some simulations in Figs. 7-18. Finally, as the future direction of this manuscript, we will consider more practical models, which will also be more complicated from the mathematical point of view.