# Dynamic Analysis for a Fractional-Order Autonomous Chaotic System

Authors: Jiangang Zhang; Juan Nan; Wenju Du; Yandong Chu; Hongwei Luo
Journal: Discrete Dynamics in Nature and Society
 (2016)

---

## Abstract

We introduce a discretization process to discretize a modified fractional-order optically injected semiconductor lasers model and investigate its dynamical behaviors. More precisely, a sufficient condition for the existence and uniqueness of the solution is obtained, and the necessary and sufficient conditions of stability of the discrete system are investigated. The results show that the system’s fractional parameter has an effect on the stability of the discrete system, and the system has rich dynamic characteristics such as Hopf bifurcation, attractor crisis, and chaotic attractors.


---

## Body

## 1. Introduction

The idea of fractional-order calculus (FOC) has been well known since the development of the regular calculus. The significant progress on FOC has been witnessed because the FOC has a wide range of applications in diffusion. In the past decades, chaotic systems have become a focal point of renewed interest for many researchers. And we can find these nonlinear systems in various natural and man-made systems, which are known to have great sensitivity to initial conditions. And we can find chaos phenomenon in various natural and man-made systems, and these chaotic systems have great sensitivity to initial conditions. Because differential equations with fractional order can be applied in many areas of science and engineering, they attracted many researchers’ attention and their complex behaviors have been widely studied in recent years. Nowadays, there is increasing interest in the subject of a fractional model which can give a more realistic interpretation of practical phenomena. Furthermore, many systems in interdisciplinary fields can be described by the fractional differential equations, such as turbulence, electromagnetism, signal processing, and quantum evolution of complex systems. It is also demonstrated that some fractional-order differential systems behave chaotically or hyperchaotically, such as fractional-order Chen system [1–3] and fractional-order Lorenz system [4, 5]. The fractional-order equations are more suitable than integer-order ones to describe the biological, economic, and social systems where memory effects are important [6]. More recently, it has been found that some differential systems with fractional order possess chaotic behavior [7–16]. And many Lorenz-like or Lorenz-based chaotic systems were proposed and investigated. Some classical 3D autonomous chaotic systems have three particular fixed points: one saddle and two unstable saddle foci [17]. And the other 3D chaotic systems have two unstable saddle foci [18, 19]. Yang and Chen [20] found another 3D chaotic system with three fixed points: one saddle and two stable equilibriums. However, an increasing number of three-dimensional chaotic systems have been found these years in many physical and engineering fields. Nonlinear dynamics of a semiconductor laser subject to optical injection is currently a hot research field due to its rich physics and complexity as well as its potential applications in communications [21, 22]. Momani et al. [23] applied the multistep generalized differential transform method to solve the fractional-order multiple chaotic FitzHugh-Nagumo (FHN) neurons model. They illustrated the algorithm by studying the dynamics of three coupled chaotic FHN neurons equations with different gap junctions under external electrical stimulation, and the fractional derivatives are described in the Caputo sense. The numerical simulation results show that only few terms are required to deduce the approximate solutions which are found to be accurate and efficient.Optically injected semiconductor lasers revealed amazingly rich behaviors like stable locking, coexistence of attractors, quasiperiodicity, instabilities, pulsations, and many routes to chaos such as period-doubling cascades, intermittency, breakup of tori, and homoclinic and heteroclinic tangencies [24–30], which deserve further systematic investigation. Wieczorek et al. [31] studied the dynamics and bifurcations of a single-mode semiconductor laser with optical injection and investigated the dependence of the dynamics on the injected field strength K and its detuning ω from the unperturbed laser resonant frequency and the linewidth enhancement factor α which seem to have the most significant influence on the dynamics. Chlouverakis and Adams [32] discussed the dynamics of a semiconductor laser subject to optical injection using the method of the Largest Lyapunov Exponent (LLE) and used this model as one of three examples to investigate the relation between the correlation (D
2) and the Kaplan-Yorke dimension (D
K
Y) of three-dimensional chaotic flows. Chu et al. [33] studied the dynamical properties of the semiconductor laser subject to optical injection and a new almost-Hamilton chaotic system with very high Lyapunov dimensions is constructed and investigated. The aim of this paper is to investigate the dynamical behavior of a discretization fractional order of a modified optically injected semiconductor lasers model.The paper is organized as follows. In Section2, we give some related preliminaries. We present the fractional-order semiconductor lasers model and introduce a discretization process to discretize a modified fractional-order optically injected semiconductor lasers model in Section 3. In Section 4, a sufficient condition for existence and uniqueness of the solution of the fractional-order system is investigated and the equilibrium points and their asymptotic stability in the fractional-order model and its discrete counterpart are also studied. In Section 5, we give some numerical simulations, which not only illustrate our results with the theoretical analysis, but also exhibit the complex dynamical behaviors of the discretization fractional order of a modified optically injected semiconductor lasers model. Finally, conclusions are given in Section 6.


## 2. Preliminaries

The Caputo definition of fractional derivative [34] is given as follows:(1)
D
α
f
t
=
I
l
-
α
f
l
t
,
α
>
0
,where l is the least integer which is not less than a and I
θ is the Riemann-Liouville integral operator of order h which is given by(2)
I
θ
u
t
=
1
Γ
θ
∫
0
t
t
-
τ
θ
-
1
u
τ
d
τ
,
θ
>
0
,where Γ
(
θ
) is Euler’s Gamma function. The operator D
α is termed “Caputo differential operator of order α.” In [35, 36], the geometric and physical interpretations of the fractional derivatives were investigated.The stability conditions and their applications to fractional-order differential equations were reported in [37–40]. The local stability of the equilibrium points of a linearized fractional-order system can be obtained from the following Matignon result [37]:(3)
arg
⁡
λ
i
>
α
π
2
,
i
=
1,2
,
3
,where λ
1, λ
2, and λ
3 are the eigenvalues of the system. Then, we consider the following nonlinear autonomous fractional-order system:(4)
D
α
X
t
=
f
X
t
,
X
0
=
X
0
,where X
(
t
)
=
(
x
1
,
x
2
,
x
3
)
T
∈
R
3 and f
:
R
3
→
R
3 is a nonlinear vector function in terms of X. The Jacobian matrix evaluated at the equilibrium point X
∗
=
(
x
1
∗
,
x
2
∗
,
x
3
∗
) is(5)
J
X
∗
=
∂
f
i
∂
x
j
i
j
X
=
X
∗
.Hence, we have the following lemma [41].Lemma 1.
If all the eigenvaluesλ
1, λ
2, and λ
3 of the equilibrium point X
∗ of system (4) satisfy Matignon’s condition (3), then X
∗ is locally asymptotically stable.


## 3. Fractional-Order Semiconductor Lasers Model and Its Discretization

Chlouverakis and Adams [32] present a modified optically injected semiconductor lasers model as follows:(6)
∂
E
x
∂
t
=
K
+
1
2
E
x
+
ω
-
1
2
α
n
E
y
,
∂
E
y
∂
t
=
-
ω
-
1
2
α
n
E
x
+
1
2
E
y
n
,
∂
n
∂
t
=
-
2
Γ
n
-
1
+
2
B
n
E
x
2
+
E
y
2
-
1
,where the complex electric field E
=
E
x
+
i
E
y and n is the normalized population inversion, K is the injection strength, ω is the detuning, and α is the linewidth enhancement factor. And the parameters B
=
ω
R
/
2
Γ
0 and Γ
=
Γ
N
/
2
ω
R
+
B, Γ
0 is the inverse photon lifetime, Γ
N is the inverse electron lifetime, and ω
R is the angular relaxation oscillation frequency. For simplification, in this paper, we established a dimensionless modified optically injected semiconductor lasers model by dimensionless method:(7)
x
˙
=
a
+
z
x
-
b
y
,
y
˙
=
z
b
x
-
c
y
,
z
˙
=
1
-
x
2
-
y
2
,where (
x
,
y
,
z
)
∈
R
3 are the state variables and (
a
,
b
,
c
)
∈
R
3 are real parameters. The bifurcation diagram of y
m
a
x versus control parameter c, the Lyapunov-exponent spectrum, and the Kaplan-Yorke dimension for specific values set (
a
=
0.4
,
b
=
3
) versus the bifurcation parameter c on the open interval (
-
8,1
) are shown in Figures 1(a)–1(c), respectively.Figure 1
(a) Bifurcation diagram for specific values set(
a
=
0.4
,
b
=
3
) versus c, (b) Lyapunov-exponent spectrum against c corresponding to the bifurcation diagram of (a), and (c) Kaplan-Yorke dimension versus the control parameter c corresponding to the bifurcation diagram of (a).
(a)
(b)
(c)System (7) is dissipative and exhibits a chaotic attractor when the parameter values a
=
0.4, b
=
3, and c
=
0.73. Figure 2 shows the chaotic attractor and the Poincaré map in x
-
z
{
y
=
0
} plane of chaotic system (7).Figure 2
(a) Time history and (b) Poincaré map inx
-
z
{
y
=
0
} plane of chaotic system (7).
(a)
(b)Next, we will investigate the following fractional order of a modified optically injected semiconductor lasers model:(8)
d
α
x
d
t
α
=
a
+
z
x
-
b
y
,
d
α
y
d
t
α
=
z
b
x
-
c
y
,
d
α
z
d
t
α
=
1
-
x
2
-
y
2
,where α is the fractional order satisfying 0
<
α
≤
1 and d
α
/
d
t
α is in the sense of the Caputo fractional derivative.And the discretization process of the fractional-order modified optically injected semiconductor lasers model is given as follows.Assume thatx
(
0
)
=
x
0, y
(
0
)
=
y
0, and z
(
0
)
=
z
0 are the initial conditions of system (8). So, the discretization of system (8) with piecewise constant arguments is given as(9)
D
α
x
t
=
a
+
z
t
s
s
x
t
s
s
-
b
y
t
s
s
,
D
α
y
t
=
z
t
s
s
b
x
t
s
s
-
c
y
t
s
s
,
D
α
z
t
=
1
-
x
t
s
s
2
-
y
t
s
s
2
.First, we lett
∈
[
0
,
s
), so t
/
s
∈
[
0,1
). Thus, we obtain(10)
D
α
x
t
=
a
+
z
0
x
0
-
b
y
0
,
D
α
y
t
=
z
0
b
x
0
-
c
y
0
,
D
α
z
t
=
1
-
x
0
2
-
y
0
2
,and the solution of (10) is reduced to(11)
x
1
t
=
x
0
+
I
α
a
+
z
0
x
0
-
b
y
0
=
x
0
+
t
α
α
Γ
α
a
+
z
0
x
0
-
b
y
0
,
y
1
t
=
y
0
+
I
α
z
0
b
x
0
-
c
y
0
=
y
0
+
t
α
α
Γ
α
z
0
b
x
0
-
c
y
0
,
z
1
t
=
z
0
+
I
α
1
-
x
0
2
-
y
0
2
=
z
0
+
t
α
α
Γ
α
1
-
x
0
2
-
y
0
2
.Second, lett
∈
[
s
,
2
s
), so t
/
s
∈
[
1,2
). Hence, we get(12)
D
α
x
t
=
a
+
z
1
x
1
-
b
y
1
,
D
α
y
t
=
z
1
b
x
1
-
c
y
1
,
D
α
z
t
=
1
-
x
1
2
-
y
1
2
,which have the following solutions:(13)
x
2
t
=
x
1
s
+
I
α
a
+
z
1
s
x
1
s
-
b
y
1
s
=
x
1
s
+
t
-
s
α
α
Γ
α
a
+
z
1
s
x
1
s
-
b
y
1
s
,
y
2
t
=
y
1
s
+
I
α
z
1
s
b
x
1
s
-
c
y
1
s
=
y
1
s
+
t
-
s
α
α
Γ
α
z
1
s
b
x
1
s
-
c
y
1
s
,
z
2
t
=
z
1
s
+
I
α
1
-
x
1
2
s
-
y
1
2
s
=
z
1
s
+
t
-
s
α
α
Γ
α
1
-
x
1
2
s
-
y
1
2
s
.Thus, after repeating the discretization processn times, we obtain(14)
x
n
+
1
t
=
x
n
n
s
+
I
α
a
+
z
n
n
s
x
n
n
s
-
b
y
n
n
s
=
x
n
n
s
+
t
-
n
s
α
α
Γ
α
a
+
z
n
n
s
x
n
n
s
-
b
y
n
n
s
,
y
n
+
1
t
=
y
n
n
s
+
I
α
z
n
n
s
b
x
n
n
s
-
c
y
n
n
s
=
y
n
n
s
+
t
-
n
s
α
α
Γ
α
z
n
n
s
b
x
n
n
s
-
c
y
n
n
s
,
z
n
+
1
t
=
z
n
n
s
+
I
α
1
-
x
n
2
n
s
-
y
n
2
n
s
=
z
n
n
s
+
t
-
n
s
α
α
Γ
α
1
-
x
n
2
n
s
-
y
n
2
n
s
,where t
∈
[
n
s
,
(
n
+
1
)
s
). For t
→
(
n
+
1
)
s, system (14) is reduced to(15)
x
n
+
1
=
x
n
+
I
α
a
+
z
n
x
n
-
b
y
n
=
x
n
+
s
α
α
Γ
α
a
+
z
n
x
n
-
b
y
n
,
y
n
+
1
=
y
n
+
I
α
z
n
b
x
n
-
c
y
n
=
y
n
+
s
α
α
Γ
α
z
n
b
x
n
-
c
y
n
,
z
n
+
1
=
z
n
+
I
α
1
-
x
n
2
-
y
n
2
=
z
n
+
s
α
α
Γ
α
1
-
x
n
2
-
y
n
2
.


## 4. Dynamical Behaviors of the Fractional-Order Semiconductor Lasers Model

### 4.1. Existence and Uniqueness of the Solution

The fractional order of the modified optically injected semiconductor lasers system can be written as(16)
D
α
X
t
=
F
X
t
,
t
∈
0
,
T
,
X
0
=
X
0
,where(17)
X
=
x
y
z
,
X
0
=
x
0
y
0
z
0
,
F
X
=
a
+
z
x
-
b
y
z
b
x
-
c
y
1
-
x
2
-
y
2
.Define the supremum norm as(18)
N
=
sup
t
∈
0
,
T
⁡
N
t
;then, the norm of the matrix M
=
[
m
i
j
[
t
]
] is defined by(19)
M
=
∑
i
,
j
sup
t
∈
0
,
T
⁡
m
i
j
t
.We investigate the existence and uniqueness of the solution in the regionΩ
×
(
0
,
T
], where(20)
Ω
=
x
,
y
,
z
:
max
⁡
x
,
y
,
z
≤
ξ
.Hence, we can get the solution of system (8) as follows:(21)
X
=
X
0
+
1
Γ
α
∫
0
t
t
-
τ
α
-
1
F
X
τ
d
τ
,so(22)
H
X
1
-
H
X
2
=
1
Γ
α
∫
0
t
t
-
τ
α
-
1
F
X
1
τ
-
F
X
2
τ
d
τ
.Thus, we get the following inequality:(23)
H
X
1
-
H
X
2
≤
1
Γ
α
∫
0
t
t
-
τ
α
-
1
F
X
1
τ
-
F
X
2
τ
d
τ
,which yields(24)
H
X
1
-
H
X
2
≤
T
α
Γ
α
+
1
max
⁡
-
ξ
e
-
ξ
,
c
ξ
e
-
c
ξ
,
0
X
1
-
X
2
,where(25)
L
=
T
α
Γ
α
+
1
max
⁡
-
ξ
e
-
ξ
,
c
ξ
e
-
c
ξ
,
0
;the mapping X
=
H
(
X
) is a contraction mapping if L
<
1, and the following theorem gives the sufficient condition for existence and uniqueness of the solution of system (8).Theorem 2.
The sufficient condition for existence and uniqueness of the solution of system (8) in the region Ω
×
(
0
,
T
] with initial conditions X
(
0
)
=
X
0 and t
∈
0
,
T is(26)
L
=
T
α
Γ
α
+
1
max
⁡
-
ξ
e
-
ξ
,
c
ξ
e
-
c
ξ
,
0
<
1
.


### 4.2. Stability of the Fixed Points of the Discrete System

In order to find the equilibrium points of system (8), we let(27)
d
α
x
d
t
α
=
0
,
d
α
y
d
t
α
=
0
,
d
α
z
d
t
α
=
0
,and from direct calculation we know that the system has at most five nonnegative equilibrium points:(28)
E
0
=
c
c
2
+
b
2
,
b
c
2
+
b
2
,
a
c
2
+
b
2
b
2
-
c
,
E
1
=
-
c
c
2
+
b
2
,
-
b
c
2
+
b
2
,
-
a
c
2
+
b
2
b
2
-
c
.By considering the Jacobian matrices for these fixed points and calculating their eigenvalues, we can investigate the local stability of these points based on the roots of the model’s characteristic equation [42]. The Jacobian matrix of system (15) is given by (29)
J
x
n
,
y
n
,
z
n
=
1
+
s
α
z
n
α
Γ
α
-
s
α
b
z
n
α
Γ
α
s
α
x
n
-
b
y
n
α
Γ
α
s
α
b
z
n
α
Γ
α
1
-
s
α
c
z
n
α
Γ
α
s
α
b
x
n
-
c
y
n
α
Γ
α
-
2
s
α
x
n
α
Γ
α
-
2
s
α
y
n
α
Γ
α
1
.Theorem 3.
For system (15), the following statements hold true:(i)
E
0 is unstable if a
c
>
0.
(ii)
E
0 is locally asymptotically stable if and only if 0
<
a
<
α
Γ
(
α
)
/
s
α
-
2
s
α
/
α
Γ
(
α
).Proof.
The Jacobian matrix evaluated at the fixed pointE
0 is(30)
J
E
0
=
1
-
s
α
a
α
Γ
α
0
s
α
α
Γ
α
0
1
+
s
α
c
a
α
Γ
α
0
2
s
α
α
Γ
α
0
1
.The eigenvalues of J
(
E
0
) are λ
2
=
1
+
s
α
c
a
/
α
Γ
(
α
) and λ
1,3
=
1
-
s
α
a
/
2
α
Γ
α
±
(
s
α
/
2
α
Γ
α
)
a
2
+
8. Since 0
<
α
≤
1 and s
>
0, s
α
/
α
Γ
(
α
), a
c
>
0, then λ
2
>
1 which implies that the fixed point E
0 is unstable. If 0
<
a
<
α
Γ
(
α
)
/
s
α
-
2
s
α
/
α
Γ
(
α
), the eigenvalues of J
(
E
0
) are λ
1,3
=
1
-
s
α
a
/
2
α
Γ
α
±
(
s
α
/
2
α
Γ
α
)
a
2
+
8
<
1 and λ
2
=
1
+
s
α
c
a
/
α
Γ
(
α
)
<
1, so the stability conditions λ
2
<
1 and λ
1,3
<
1 are satisfied if 0
<
a
<
α
Γ
(
α
)
/
s
α
-
2
s
α
/
α
Γ
(
α
).


#### 4.2.1. The Caseb
=
0

Theorem 4.
For system (15), the following statements hold true:(i)
E
1 is unstable if a
c
<
0.
(ii)
E
1 is locally asymptotically stable if and only if a
<
-
2
2.Proof.
The Jacobian matrix evaluated at the fixed pointE
1 is(31)
J
E
1
=
1
+
s
α
a
α
Γ
α
0
-
s
α
α
Γ
α
0
1
-
s
α
c
a
α
Γ
α
0
2
s
α
α
Γ
α
0
1
.The eigenvalues of J
E
1 are λ
2
=
1
-
s
α
c
a
/
α
Γ
(
α
) and λ
1,3
=
1
+
s
α
a
/
2
α
Γ
(
α
)
±
(
s
α
/
2
α
Γ
α
)
a
2
-
8. Since 0
<
α
≤
1 and s
>
0, s
α
/
α
Γ
(
α
), a
c
<
0, then λ
2
>
1 which implies that the fixed point E
0 is unstable. If a
<
-
2
2, the eigenvalues of J
(
E
1
) are λ
2
=
1
+
s
α
c
a
/
α
Γ
(
α
)
<
1 and λ
1,3
=
1
+
s
α
a
/
2
α
Γ
(
α
)
±
(
s
α
/
2
α
Γ
α
)
a
2
-
8
<
1, so the stability conditions λ
2
<
1 and λ
1,3
<
1 are satisfied if a
<
-
2
2.


#### 4.2.2. The Caseb
≠
0

Next, we will discuss the stability of the fixed pointE
1. The Jacobian matrix of system (15) evaluated at E
1 is(32)
J
E
1
=
γ
11
γ
12
γ
13
γ
21
γ
22
γ
23
γ
31
γ
32
γ
33
,where(33)
γ
11
=
1
+
s
α
a
c
2
+
b
2
α
b
2
-
c
Γ
α
,
γ
12
=
-
s
α
b
a
c
2
+
b
2
α
b
2
-
c
Γ
α
,
γ
13
=
-
s
α
b
2
-
c
α
Γ
α
c
2
+
b
2
,
γ
21
=
s
α
b
c
2
+
b
2
α
b
2
-
c
Γ
α
,
γ
22
=
1
-
s
α
c
a
c
2
+
b
2
α
b
2
-
c
Γ
α
,
γ
23
=
0
,
γ
31
=
-
2
s
α
c
α
c
2
+
b
2
Γ
α
,
γ
32
=
-
2
s
α
b
α
c
2
+
b
2
Γ
α
,
γ
33
=
1
.The characteristic equation ofE
1 is given by(34)
λ
3
+
A
1
λ
2
+
A
2
λ
+
A
3
=
0
,where(35)
A
1
=
2
-
V
1
,
A
2
=
-
3
+
V
2
-
V
3
-
V
4
+
V
5
,
A
3
=
-
1
-
V
6
+
V
2
-
V
7
-
V
8
,
V
1
=
s
α
c
a
c
2
+
b
2
α
b
2
-
c
Γ
α
,
V
2
=
2
s
2
α
c
b
2
-
c
α
2
Γ
2
α
c
2
+
b
2
,
V
3
=
s
α
a
c
2
+
b
2
α
b
2
-
c
Γ
α
,
V
4
=
s
2
α
b
2
a
c
2
+
b
2
α
2
b
2
-
c
2
Γ
2
α
,
V
5
=
s
2
α
c
a
2
c
2
+
b
2
α
2
b
2
-
c
2
Γ
2
α
,
V
6
=
-
2
s
3
α
b
2
+
2
a
c
2
α
3
c
2
+
b
2
Γ
3
α
,
V
7
=
s
2
α
b
2
-
a
c
a
c
2
+
b
2
α
2
b
2
-
c
2
Γ
2
α
,
V
8
=
s
α
1
-
c
a
c
2
+
b
2
α
b
2
-
c
Γ
α
.According to Jury’s criterion [42], the fixed point E
1 is locally asymptotically stable if(36)
1
+
A
1
+
A
2
+
A
3
>
0
,
1
-
A
1
+
A
2
-
A
3
>
0
,
1
-
A
2
+
A
1
A
3
-
A
3
2
>
0
,
1
+
A
2
-
A
1
A
3
-
A
3
2
>
0
.Thus, conditions (36) imply that the fixed point E
0 is locally asymptotically stable if the system’s parameters belong to the set(37)
a
,
b
,
c
,
α
,
s
:
-
1
+
V
6
-
V
2
+
V
7
+
V
8
∈
Ω
i
∀
i
=
1,2
,
3
,where(38)
Ω
1
=
V
1
-
V
2
+
V
3
+
V
4
-
V
5
,
-
4
+
V
1
+
V
2
-
V
3
-
V
4
+
V
5
,
Ω
2
=
2
-
V
1
-
2
-
V
1
2
-
4
-
4
+
V
2
-
V
3
-
V
4
+
V
5
2
,
2
-
V
1
-
2
-
V
1
2
-
4
-
4
+
V
2
-
V
3
-
V
4
+
V
5
2
,
Ω
3
=
-
2
+
V
1
-
2
-
V
1
2
-
4
-
4
+
V
2
-
V
3
-
V
4
+
V
5
2
,
-
2
+
V
1
+
2
-
V
1
2
-
4
-
4
+
V
2
-
V
3
-
V
4
+
V
5
2
.Lemma 5 (an explicit criterion of Hopf bifurcation in maps [43]).
LetX
∗ be a fixed point of the nth-order discrete-time dynamical system; the characteristic equation of the Jacobian matrix A
=
a
i
j
n
×
n is given as(39)
P
v
λ
=
λ
n
+
s
1
λ
n
-
1
+
⋯
+
s
n
-
1
λ
+
s
n
,where s
j
=
s
j
v
,
k, j
=
1
,
…
,
n, v is the bifurcation parameter, and k is the control parameter or another parameter to be determined. Consider the sequence of determinants Δ
0
±
v
,
k
=
1, Δ
1
±
v
,
k
,
…
,
Δ
n
±
v
,
k, where(40)
Δ
j
±
v
,
k
=
1
s
1
s
2
⋯
s
j
-
1
0
1
s
1
⋯
s
j
-
2
0
0
1
⋯
s
j
-
3
⋮
⋮
⋮
⋯
⋮
0
0
0
⋯
1
±
s
n
-
j
+
1
s
n
-
j
+
2
⋯
s
n
-
1
s
n
s
n
-
j
+
1
s
n
-
j
+
3
⋯
s
n
0
⋮
⋮
⋯
⋮
⋮
s
n
-
1
s
n
⋯
0
0
s
n
0
⋯
0
0
,
j
=
1
,
…
,
n
.However, the following conditions hold:(H1)
Eigenvalue assignmentΔ
n
-
1
-
v
0
,
k
=
0, P
v
0
1
>
0, -
1
n
P
v
0
-
1
>
0, Δ
n
-
1
-
v
0
,
k
>
0, Δ
n
-
1
+
v
0
,
k
=
0, Δ
j
±
v
0
,
k
>
0, j
=
n
-
3
,
n
-
5
,
…
,
1
o
r
2 when n is even (or odd, resp.).
(H2)
Transversality conditiond
Δ
n
-
1
-
v
0
,
k
/
d
v
≠
0.
(H3)
Nonresonance conditioncos
⁡
2
π
/
m
≠
φ or resonance condition cos
⁡
2
π
/
m
=
φ, where m
=
3,4
,
5
,
… and φ
=
1
-
0.5
P
v
0
1
Δ
n
-
3
-
v
0
,
k
/
Δ
n
-
2
+
v
0
,
k; then Hopf bifurcation occurs at v
0.


## 4.1. Existence and Uniqueness of the Solution

The fractional order of the modified optically injected semiconductor lasers system can be written as(16)
D
α
X
t
=
F
X
t
,
t
∈
0
,
T
,
X
0
=
X
0
,where(17)
X
=
x
y
z
,
X
0
=
x
0
y
0
z
0
,
F
X
=
a
+
z
x
-
b
y
z
b
x
-
c
y
1
-
x
2
-
y
2
.Define the supremum norm as(18)
N
=
sup
t
∈
0
,
T
⁡
N
t
;then, the norm of the matrix M
=
[
m
i
j
[
t
]
] is defined by(19)
M
=
∑
i
,
j
sup
t
∈
0
,
T
⁡
m
i
j
t
.We investigate the existence and uniqueness of the solution in the regionΩ
×
(
0
,
T
], where(20)
Ω
=
x
,
y
,
z
:
max
⁡
x
,
y
,
z
≤
ξ
.Hence, we can get the solution of system (8) as follows:(21)
X
=
X
0
+
1
Γ
α
∫
0
t
t
-
τ
α
-
1
F
X
τ
d
τ
,so(22)
H
X
1
-
H
X
2
=
1
Γ
α
∫
0
t
t
-
τ
α
-
1
F
X
1
τ
-
F
X
2
τ
d
τ
.Thus, we get the following inequality:(23)
H
X
1
-
H
X
2
≤
1
Γ
α
∫
0
t
t
-
τ
α
-
1
F
X
1
τ
-
F
X
2
τ
d
τ
,which yields(24)
H
X
1
-
H
X
2
≤
T
α
Γ
α
+
1
max
⁡
-
ξ
e
-
ξ
,
c
ξ
e
-
c
ξ
,
0
X
1
-
X
2
,where(25)
L
=
T
α
Γ
α
+
1
max
⁡
-
ξ
e
-
ξ
,
c
ξ
e
-
c
ξ
,
0
;the mapping X
=
H
(
X
) is a contraction mapping if L
<
1, and the following theorem gives the sufficient condition for existence and uniqueness of the solution of system (8).Theorem 2.
The sufficient condition for existence and uniqueness of the solution of system (8) in the region Ω
×
(
0
,
T
] with initial conditions X
(
0
)
=
X
0 and t
∈
0
,
T is(26)
L
=
T
α
Γ
α
+
1
max
⁡
-
ξ
e
-
ξ
,
c
ξ
e
-
c
ξ
,
0
<
1
.


## 4.2. Stability of the Fixed Points of the Discrete System

In order to find the equilibrium points of system (8), we let(27)
d
α
x
d
t
α
=
0
,
d
α
y
d
t
α
=
0
,
d
α
z
d
t
α
=
0
,and from direct calculation we know that the system has at most five nonnegative equilibrium points:(28)
E
0
=
c
c
2
+
b
2
,
b
c
2
+
b
2
,
a
c
2
+
b
2
b
2
-
c
,
E
1
=
-
c
c
2
+
b
2
,
-
b
c
2
+
b
2
,
-
a
c
2
+
b
2
b
2
-
c
.By considering the Jacobian matrices for these fixed points and calculating their eigenvalues, we can investigate the local stability of these points based on the roots of the model’s characteristic equation [42]. The Jacobian matrix of system (15) is given by (29)
J
x
n
,
y
n
,
z
n
=
1
+
s
α
z
n
α
Γ
α
-
s
α
b
z
n
α
Γ
α
s
α
x
n
-
b
y
n
α
Γ
α
s
α
b
z
n
α
Γ
α
1
-
s
α
c
z
n
α
Γ
α
s
α
b
x
n
-
c
y
n
α
Γ
α
-
2
s
α
x
n
α
Γ
α
-
2
s
α
y
n
α
Γ
α
1
.Theorem 3.
For system (15), the following statements hold true:(i)
E
0 is unstable if a
c
>
0.
(ii)
E
0 is locally asymptotically stable if and only if 0
<
a
<
α
Γ
(
α
)
/
s
α
-
2
s
α
/
α
Γ
(
α
).Proof.
The Jacobian matrix evaluated at the fixed pointE
0 is(30)
J
E
0
=
1
-
s
α
a
α
Γ
α
0
s
α
α
Γ
α
0
1
+
s
α
c
a
α
Γ
α
0
2
s
α
α
Γ
α
0
1
.The eigenvalues of J
(
E
0
) are λ
2
=
1
+
s
α
c
a
/
α
Γ
(
α
) and λ
1,3
=
1
-
s
α
a
/
2
α
Γ
α
±
(
s
α
/
2
α
Γ
α
)
a
2
+
8. Since 0
<
α
≤
1 and s
>
0, s
α
/
α
Γ
(
α
), a
c
>
0, then λ
2
>
1 which implies that the fixed point E
0 is unstable. If 0
<
a
<
α
Γ
(
α
)
/
s
α
-
2
s
α
/
α
Γ
(
α
), the eigenvalues of J
(
E
0
) are λ
1,3
=
1
-
s
α
a
/
2
α
Γ
α
±
(
s
α
/
2
α
Γ
α
)
a
2
+
8
<
1 and λ
2
=
1
+
s
α
c
a
/
α
Γ
(
α
)
<
1, so the stability conditions λ
2
<
1 and λ
1,3
<
1 are satisfied if 0
<
a
<
α
Γ
(
α
)
/
s
α
-
2
s
α
/
α
Γ
(
α
).


### 4.2.1. The Caseb
=
0

Theorem 4.
For system (15), the following statements hold true:(i)
E
1 is unstable if a
c
<
0.
(ii)
E
1 is locally asymptotically stable if and only if a
<
-
2
2.Proof.
The Jacobian matrix evaluated at the fixed pointE
1 is(31)
J
E
1
=
1
+
s
α
a
α
Γ
α
0
-
s
α
α
Γ
α
0
1
-
s
α
c
a
α
Γ
α
0
2
s
α
α
Γ
α
0
1
.The eigenvalues of J
E
1 are λ
2
=
1
-
s
α
c
a
/
α
Γ
(
α
) and λ
1,3
=
1
+
s
α
a
/
2
α
Γ
(
α
)
±
(
s
α
/
2
α
Γ
α
)
a
2
-
8. Since 0
<
α
≤
1 and s
>
0, s
α
/
α
Γ
(
α
), a
c
<
0, then λ
2
>
1 which implies that the fixed point E
0 is unstable. If a
<
-
2
2, the eigenvalues of J
(
E
1
) are λ
2
=
1
+
s
α
c
a
/
α
Γ
(
α
)
<
1 and λ
1,3
=
1
+
s
α
a
/
2
α
Γ
(
α
)
±
(
s
α
/
2
α
Γ
α
)
a
2
-
8
<
1, so the stability conditions λ
2
<
1 and λ
1,3
<
1 are satisfied if a
<
-
2
2.


### 4.2.2. The Caseb
≠
0

Next, we will discuss the stability of the fixed pointE
1. The Jacobian matrix of system (15) evaluated at E
1 is(32)
J
E
1
=
γ
11
γ
12
γ
13
γ
21
γ
22
γ
23
γ
31
γ
32
γ
33
,where(33)
γ
11
=
1
+
s
α
a
c
2
+
b
2
α
b
2
-
c
Γ
α
,
γ
12
=
-
s
α
b
a
c
2
+
b
2
α
b
2
-
c
Γ
α
,
γ
13
=
-
s
α
b
2
-
c
α
Γ
α
c
2
+
b
2
,
γ
21
=
s
α
b
c
2
+
b
2
α
b
2
-
c
Γ
α
,
γ
22
=
1
-
s
α
c
a
c
2
+
b
2
α
b
2
-
c
Γ
α
,
γ
23
=
0
,
γ
31
=
-
2
s
α
c
α
c
2
+
b
2
Γ
α
,
γ
32
=
-
2
s
α
b
α
c
2
+
b
2
Γ
α
,
γ
33
=
1
.The characteristic equation ofE
1 is given by(34)
λ
3
+
A
1
λ
2
+
A
2
λ
+
A
3
=
0
,where(35)
A
1
=
2
-
V
1
,
A
2
=
-
3
+
V
2
-
V
3
-
V
4
+
V
5
,
A
3
=
-
1
-
V
6
+
V
2
-
V
7
-
V
8
,
V
1
=
s
α
c
a
c
2
+
b
2
α
b
2
-
c
Γ
α
,
V
2
=
2
s
2
α
c
b
2
-
c
α
2
Γ
2
α
c
2
+
b
2
,
V
3
=
s
α
a
c
2
+
b
2
α
b
2
-
c
Γ
α
,
V
4
=
s
2
α
b
2
a
c
2
+
b
2
α
2
b
2
-
c
2
Γ
2
α
,
V
5
=
s
2
α
c
a
2
c
2
+
b
2
α
2
b
2
-
c
2
Γ
2
α
,
V
6
=
-
2
s
3
α
b
2
+
2
a
c
2
α
3
c
2
+
b
2
Γ
3
α
,
V
7
=
s
2
α
b
2
-
a
c
a
c
2
+
b
2
α
2
b
2
-
c
2
Γ
2
α
,
V
8
=
s
α
1
-
c
a
c
2
+
b
2
α
b
2
-
c
Γ
α
.According to Jury’s criterion [42], the fixed point E
1 is locally asymptotically stable if(36)
1
+
A
1
+
A
2
+
A
3
>
0
,
1
-
A
1
+
A
2
-
A
3
>
0
,
1
-
A
2
+
A
1
A
3
-
A
3
2
>
0
,
1
+
A
2
-
A
1
A
3
-
A
3
2
>
0
.Thus, conditions (36) imply that the fixed point E
0 is locally asymptotically stable if the system’s parameters belong to the set(37)
a
,
b
,
c
,
α
,
s
:
-
1
+
V
6
-
V
2
+
V
7
+
V
8
∈
Ω
i
∀
i
=
1,2
,
3
,where(38)
Ω
1
=
V
1
-
V
2
+
V
3
+
V
4
-
V
5
,
-
4
+
V
1
+
V
2
-
V
3
-
V
4
+
V
5
,
Ω
2
=
2
-
V
1
-
2
-
V
1
2
-
4
-
4
+
V
2
-
V
3
-
V
4
+
V
5
2
,
2
-
V
1
-
2
-
V
1
2
-
4
-
4
+
V
2
-
V
3
-
V
4
+
V
5
2
,
Ω
3
=
-
2
+
V
1
-
2
-
V
1
2
-
4
-
4
+
V
2
-
V
3
-
V
4
+
V
5
2
,
-
2
+
V
1
+
2
-
V
1
2
-
4
-
4
+
V
2
-
V
3
-
V
4
+
V
5
2
.Lemma 5 (an explicit criterion of Hopf bifurcation in maps [43]).
LetX
∗ be a fixed point of the nth-order discrete-time dynamical system; the characteristic equation of the Jacobian matrix A
=
a
i
j
n
×
n is given as(39)
P
v
λ
=
λ
n
+
s
1
λ
n
-
1
+
⋯
+
s
n
-
1
λ
+
s
n
,where s
j
=
s
j
v
,
k, j
=
1
,
…
,
n, v is the bifurcation parameter, and k is the control parameter or another parameter to be determined. Consider the sequence of determinants Δ
0
±
v
,
k
=
1, Δ
1
±
v
,
k
,
…
,
Δ
n
±
v
,
k, where(40)
Δ
j
±
v
,
k
=
1
s
1
s
2
⋯
s
j
-
1
0
1
s
1
⋯
s
j
-
2
0
0
1
⋯
s
j
-
3
⋮
⋮
⋮
⋯
⋮
0
0
0
⋯
1
±
s
n
-
j
+
1
s
n
-
j
+
2
⋯
s
n
-
1
s
n
s
n
-
j
+
1
s
n
-
j
+
3
⋯
s
n
0
⋮
⋮
⋯
⋮
⋮
s
n
-
1
s
n
⋯
0
0
s
n
0
⋯
0
0
,
j
=
1
,
…
,
n
.However, the following conditions hold:(H1)
Eigenvalue assignmentΔ
n
-
1
-
v
0
,
k
=
0, P
v
0
1
>
0, -
1
n
P
v
0
-
1
>
0, Δ
n
-
1
-
v
0
,
k
>
0, Δ
n
-
1
+
v
0
,
k
=
0, Δ
j
±
v
0
,
k
>
0, j
=
n
-
3
,
n
-
5
,
…
,
1
o
r
2 when n is even (or odd, resp.).
(H2)
Transversality conditiond
Δ
n
-
1
-
v
0
,
k
/
d
v
≠
0.
(H3)
Nonresonance conditioncos
⁡
2
π
/
m
≠
φ or resonance condition cos
⁡
2
π
/
m
=
φ, where m
=
3,4
,
5
,
… and φ
=
1
-
0.5
P
v
0
1
Δ
n
-
3
-
v
0
,
k
/
Δ
n
-
2
+
v
0
,
k; then Hopf bifurcation occurs at v
0.


## 4.2.1. The Caseb
=
0

Theorem 4.
For system (15), the following statements hold true:(i)
E
1 is unstable if a
c
<
0.
(ii)
E
1 is locally asymptotically stable if and only if a
<
-
2
2.Proof.
The Jacobian matrix evaluated at the fixed pointE
1 is(31)
J
E
1
=
1
+
s
α
a
α
Γ
α
0
-
s
α
α
Γ
α
0
1
-
s
α
c
a
α
Γ
α
0
2
s
α
α
Γ
α
0
1
.The eigenvalues of J
E
1 are λ
2
=
1
-
s
α
c
a
/
α
Γ
(
α
) and λ
1,3
=
1
+
s
α
a
/
2
α
Γ
(
α
)
±
(
s
α
/
2
α
Γ
α
)
a
2
-
8. Since 0
<
α
≤
1 and s
>
0, s
α
/
α
Γ
(
α
), a
c
<
0, then λ
2
>
1 which implies that the fixed point E
0 is unstable. If a
<
-
2
2, the eigenvalues of J
(
E
1
) are λ
2
=
1
+
s
α
c
a
/
α
Γ
(
α
)
<
1 and λ
1,3
=
1
+
s
α
a
/
2
α
Γ
(
α
)
±
(
s
α
/
2
α
Γ
α
)
a
2
-
8
<
1, so the stability conditions λ
2
<
1 and λ
1,3
<
1 are satisfied if a
<
-
2
2.


## 4.2.2. The Caseb
≠
0

Next, we will discuss the stability of the fixed pointE
1. The Jacobian matrix of system (15) evaluated at E
1 is(32)
J
E
1
=
γ
11
γ
12
γ
13
γ
21
γ
22
γ
23
γ
31
γ
32
γ
33
,where(33)
γ
11
=
1
+
s
α
a
c
2
+
b
2
α
b
2
-
c
Γ
α
,
γ
12
=
-
s
α
b
a
c
2
+
b
2
α
b
2
-
c
Γ
α
,
γ
13
=
-
s
α
b
2
-
c
α
Γ
α
c
2
+
b
2
,
γ
21
=
s
α
b
c
2
+
b
2
α
b
2
-
c
Γ
α
,
γ
22
=
1
-
s
α
c
a
c
2
+
b
2
α
b
2
-
c
Γ
α
,
γ
23
=
0
,
γ
31
=
-
2
s
α
c
α
c
2
+
b
2
Γ
α
,
γ
32
=
-
2
s
α
b
α
c
2
+
b
2
Γ
α
,
γ
33
=
1
.The characteristic equation ofE
1 is given by(34)
λ
3
+
A
1
λ
2
+
A
2
λ
+
A
3
=
0
,where(35)
A
1
=
2
-
V
1
,
A
2
=
-
3
+
V
2
-
V
3
-
V
4
+
V
5
,
A
3
=
-
1
-
V
6
+
V
2
-
V
7
-
V
8
,
V
1
=
s
α
c
a
c
2
+
b
2
α
b
2
-
c
Γ
α
,
V
2
=
2
s
2
α
c
b
2
-
c
α
2
Γ
2
α
c
2
+
b
2
,
V
3
=
s
α
a
c
2
+
b
2
α
b
2
-
c
Γ
α
,
V
4
=
s
2
α
b
2
a
c
2
+
b
2
α
2
b
2
-
c
2
Γ
2
α
,
V
5
=
s
2
α
c
a
2
c
2
+
b
2
α
2
b
2
-
c
2
Γ
2
α
,
V
6
=
-
2
s
3
α
b
2
+
2
a
c
2
α
3
c
2
+
b
2
Γ
3
α
,
V
7
=
s
2
α
b
2
-
a
c
a
c
2
+
b
2
α
2
b
2
-
c
2
Γ
2
α
,
V
8
=
s
α
1
-
c
a
c
2
+
b
2
α
b
2
-
c
Γ
α
.According to Jury’s criterion [42], the fixed point E
1 is locally asymptotically stable if(36)
1
+
A
1
+
A
2
+
A
3
>
0
,
1
-
A
1
+
A
2
-
A
3
>
0
,
1
-
A
2
+
A
1
A
3
-
A
3
2
>
0
,
1
+
A
2
-
A
1
A
3
-
A
3
2
>
0
.Thus, conditions (36) imply that the fixed point E
0 is locally asymptotically stable if the system’s parameters belong to the set(37)
a
,
b
,
c
,
α
,
s
:
-
1
+
V
6
-
V
2
+
V
7
+
V
8
∈
Ω
i
∀
i
=
1,2
,
3
,where(38)
Ω
1
=
V
1
-
V
2
+
V
3
+
V
4
-
V
5
,
-
4
+
V
1
+
V
2
-
V
3
-
V
4
+
V
5
,
Ω
2
=
2
-
V
1
-
2
-
V
1
2
-
4
-
4
+
V
2
-
V
3
-
V
4
+
V
5
2
,
2
-
V
1
-
2
-
V
1
2
-
4
-
4
+
V
2
-
V
3
-
V
4
+
V
5
2
,
Ω
3
=
-
2
+
V
1
-
2
-
V
1
2
-
4
-
4
+
V
2
-
V
3
-
V
4
+
V
5
2
,
-
2
+
V
1
+
2
-
V
1
2
-
4
-
4
+
V
2
-
V
3
-
V
4
+
V
5
2
.Lemma 5 (an explicit criterion of Hopf bifurcation in maps [43]).
LetX
∗ be a fixed point of the nth-order discrete-time dynamical system; the characteristic equation of the Jacobian matrix A
=
a
i
j
n
×
n is given as(39)
P
v
λ
=
λ
n
+
s
1
λ
n
-
1
+
⋯
+
s
n
-
1
λ
+
s
n
,where s
j
=
s
j
v
,
k, j
=
1
,
…
,
n, v is the bifurcation parameter, and k is the control parameter or another parameter to be determined. Consider the sequence of determinants Δ
0
±
v
,
k
=
1, Δ
1
±
v
,
k
,
…
,
Δ
n
±
v
,
k, where(40)
Δ
j
±
v
,
k
=
1
s
1
s
2
⋯
s
j
-
1
0
1
s
1
⋯
s
j
-
2
0
0
1
⋯
s
j
-
3
⋮
⋮
⋮
⋯
⋮
0
0
0
⋯
1
±
s
n
-
j
+
1
s
n
-
j
+
2
⋯
s
n
-
1
s
n
s
n
-
j
+
1
s
n
-
j
+
3
⋯
s
n
0
⋮
⋮
⋯
⋮
⋮
s
n
-
1
s
n
⋯
0
0
s
n
0
⋯
0
0
,
j
=
1
,
…
,
n
.However, the following conditions hold:(H1)
Eigenvalue assignmentΔ
n
-
1
-
v
0
,
k
=
0, P
v
0
1
>
0, -
1
n
P
v
0
-
1
>
0, Δ
n
-
1
-
v
0
,
k
>
0, Δ
n
-
1
+
v
0
,
k
=
0, Δ
j
±
v
0
,
k
>
0, j
=
n
-
3
,
n
-
5
,
…
,
1
o
r
2 when n is even (or odd, resp.).
(H2)
Transversality conditiond
Δ
n
-
1
-
v
0
,
k
/
d
v
≠
0.
(H3)
Nonresonance conditioncos
⁡
2
π
/
m
≠
φ or resonance condition cos
⁡
2
π
/
m
=
φ, where m
=
3,4
,
5
,
… and φ
=
1
-
0.5
P
v
0
1
Δ
n
-
3
-
v
0
,
k
/
Δ
n
-
2
+
v
0
,
k; then Hopf bifurcation occurs at v
0.


## 5. Numerical Simulations

Even if it is possible to obtain the analytical conditions, the processing is very difficult. So, in order to analyze the stability of the interior fixed points, we investigate the global dynamical behavior of system (15) through numerical simulations. And we show some attractors of system (15), such as asymptotic behaviors near the equilibrium points, periodic orbits, and strange attractors. In the following, we will use the numerical simulations to show the influence of different fractional-order parameters q and s on the stability of system (15).We fixeda
=
0.4 and b
=
3 and vary the parameters q
,
s, and c. Let s
=
0.01. Figures 3(a)–3(c) show the portraits of system (15) with q
=
0.87, q
=
0.93, and q
=
0.98, respectively. The numerical results show that with the increase of fractional-order parameter q the discrete system (15) is consistent with the fractional-order system (8).Figure 3
Trajectory of system (15) in phase space when a
=
0.4 and b
=
3.
(a)
c
=
-
0.1, q
=
0.87, and s
=
0.01
(b)
c
=
-
1.7, q
=
0.93, and s
=
0.01
(c)
c
=
-
0.8, q
=
0.98, and s
=
0.01
(d)
c
=
-
1.9, q
=
0.7, and s
=
0.01
(e)
c
=
-
1, q
=
0.8, and s
=
0.01
(f)
c
=
-
0.1, q
=
0.87, and s
=
0.01
(g)
c
=
-
0.8, q
=
0.9, and s
=
0.001
(h)
c
=
-
0.87, q
=
0.93, and s
=
0.001
(i)
c
=
-
0.1, q
=
0.97, and s
=
0.001
(j)
c
=
-
1.9, q
=
0.98, and s
=
0.001When we fixeds
=
0.01, Figures 3(d)–3(f) depict the attractors of system (15) with q
=
0.7, q
=
0.8, and q
=
0.87, respectively. And when s
=
0.001, we present the phase portraits of system (15) with q
=
0.9, q
=
0.93, and q
=
0.97, as shown in Figures 3(g)–3(i), respectively. From Figures 3(d)–3(i), we notice that the chaotic behavior of system (15) gradually becomes stable with the reduction of the parameter s and the gradual increase in the fractional-order parameter q.Figure3(j) is the strange attractor of system (15) with q
=
0.98 and s
=
0.001. Obviously, it confirms the above result; this is because Figure 3(j) shows that the reversibility of the above result, namely, the stable behavior of system (15), is destabilizing when s is small and decreasing the fractional-order parameter q.We let the time series length be equal to 3000. Figure4(a) shows the curve of log
⁡
C
(
r
,
m
)
-
log
⁡
r when the parameter m varies in [
1,18
], and the error bars are purely statistical. The curves attain a constant slope (the same for all curves) as r
→
0 and the slope should be equal to the information dimension. The slope curve of the time serials is shown in Figure 4(b). In Figure 4(c), with the various embedding dimension, correlation exponent is on the increase. When m
≥
8, the sharp increase to very big values indicates that D
2
=
1.1466. When embedding dimension increases to 7, the curve above it appears to converge and we compute the Kolmogorov entropy K
=
0.3096. From Figure 4(d), we found that the curves appear to converge to this curve at least up to the accuracy that we could hope to achieve in such numerical experiments. In general, results were found to be more sensitive to the number of iterations than to the averaging number.Figure 4
AtN
=
3000. (a) Correlation integral curve. Points pertaining to the same value of m are connected by lines. The values of m are m
=
1 (top curve), 4,6
,
…
,
18 (bottom curve); (b) slope curve of the time serials. An ordinate value is the slope of log
⁡
C
(
r
,
m
) and log
⁡
r; (c) correlation exponent curve D
2
=
1.1466; (d) the Kolmogorov entropy of system (15) attractor is equal to D
2
=
0.3096.
(a)
(b)
(c)
(d)From the above numerical experiments, we get the following results:(1)
When embedding dimension varies from 6 to 8, the correlation integral begins to gradually saturate. A correlation exponent of the different system verified by the number of the time series needs specific analysis.
(2)
If embedding dimension does not exceed the saturated values of the minimum embedding dimension, then the estimation of the correlation exponents is not up to the precision with which we could hope to extract the metric entropy.
(3)
The constancy of the slopes will break down forr smaller than some value r
0. With the incensement of the embedding dimension, the distances between points increase, and r
0 shifts to higher values of r. When the embedding dimension is to achieve a certain value, the Kolmogorov entropy reaches saturation and begins to converge.However, when the bifurcation diagram loses its continuity, this means that the state of the system is either quasiperiodic or chaotic, as shown in Figures5(a) and 5(b). And, according to Lemma 5, we select a
=
0.4, b
=
3, q
=
0.98, and s
=
0.01 and draw the phase diagram of system (15) with c
=
-
0.8, c
=
-
1, and c
=
-
1.9, as shown in Figures 6–8, respectively. From Figures 6–8, we know that system (15) loses its stability with the gradual decrease of bifurcation parameter c. That is to say, by choosing an appropriate bifurcation parameter, we can prove that the Hopf bifurcation occurs when the bifurcation parameter passes through the critical value.Figure 5
(a) Phase trajectory in 3D space witha
=
0.4, b
=
3, c
=
-
1.9, q
=
0.96, and s
=
0.001; (b) interior phase trajectory of (a).
(a)
(b)Figure 6
Phase diagram of system (15) with a
=
0.4, b
=
3, and c
=
-
0.8.
(a)
(b)Figure 7
Phase diagram of system (15) with a
=
0.4, b
=
3, and c
=
-
1.
(a)
(b)Figure 8
Phase diagram of system (15) with a
=
0.4, b
=
3, and c
=
-
1.9.
(a)
(b)


## 6. Conclusion

In this paper, we introduced a fractional order of modified optically injected semiconductor lasers model and discretized this system by using a new discretization technique. More precisely, a sufficient condition for existence and uniqueness of the solution of the proposed fractional-order system is investigated, and we also studied the local stability of the equilibrium of the discrete fractional-order semiconductor lasers system. Moreover, the results showed that the fractional parameterq has an effect on the stability of the discrete system. And the chaotic behavior of system (15) will be stabilized when reducing the value of parameter s and increasing the fractional-order parameter q. Meanwhile, the stable behavior of system (15) is destabilizing when decreasing the fractional-order parameter q. We also found that the discrete system exhibits much richer dynamical behaviors than its corresponding fractional-order counterpart, and the existence of Hopf bifurcation of discrete system has been verified. Finally, the numerical simulations are given to show the chaotic attractors of the fractional-order system and the richer dynamics of its discrete counterpart. In addition, we proposed a method to estimate the Kolmogorov entropy and use it to verify the state of the system. Future work on the topic may be extended to apply the chaos control (OGY method, feedback control method, and adaptive control method) or bifurcation control (the linear and nonlinear feedback method, washout-filter method, and frequency domain analysis and approximation method) on the proposed systems and include examples to demonstrate actual solutions of the systems.



---

*Source: 8712496-2016-07-21.xml*