\documentclass[reqno]{amsart}
\usepackage{hyperref}
\usepackage{graphicx} 

\AtBeginDocument{{\noindent\small
\emph{Electronic Journal of Differential Equations},
Vol. 2017 (2017), No. 128, pp. 1--15.\newline
ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu}
\thanks{\copyright 2017 Texas State University.}
\vspace{8mm}}

\begin{document}
\title[\hfilneg EJDE-2017/128\hfil Coexistence of synchronization types]
{Coexistence of some chaos synchronization types in fractional-order
differential equations}

\author[A. Ouannas, S. Abdelmalek, S. Bendoukha \hfil EJDE-2017/128\hfilneg]
{Adel Ouannas, Salem Abdelmalek, Samir Bendoukha}

\address{Adel Ouannas \newline
Laboratory of Mathematics, Informatics and Systems (LAMIS),
University of Larbi Tebessi, Tebessa, 12002 Algeria}
\email{ouannas.a@yahoo.com}

\address{Salem Abdelmalek \newline
Department of mathematics,
University of Tebessa 12002 Algeria}
\email{sallllm@gmail.com}

\address{Samir Bendoukha \newline
Electrical Engineering Department,
College of Engineering at Yanbu,
Taibah University, Saudi Arabia}
\email{sbendoukha@taibahu.edu.sa}

\dedicatory{Communicated by Mokhtar Kirane}

\thanks{Submitted March 10, 2017. Published May 10, 2017.}
\subjclass[2010]{34A08, 34H10, 34D06}
\keywords{Chaos synchronization; fractional-order systems; coexistence;
\hfill\break\indent fractional Lyapunov approach}

\begin{abstract}
 Referring to incommensurate and commensurate fractional systems, this
 article presents a new approach to investigate the coexistence of some
 synchronization types between non-identical systems characterized by
 different dimensions and different orders. In particular, the paper shows
 that complete synchronization (CS), anti-synchronization (AS) and inverse
 full state hybrid function projective synchronization (IFSHFPS) coexist
 when synchronizing a three-dimensional master system with a four-dimensional
 slave system. The approach is based on two new results involving stability
 theory of linear fractional systems and the fractional Lyapunov method.
 A number of examples are provided to highlight the applicability of the
 method.
\end{abstract}

\maketitle
\numberwithin{equation}{section}
\newtheorem{theorem}{Theorem}[section]
\newtheorem{lemma}[theorem]{Lemma}
\newtheorem{proposition}[theorem]{Proposition}
\newtheorem{definition}[theorem]{Definition}
\allowdisplaybreaks


\section{Introduction} \label{sec1}

Over the last few years, substantial efforts have been devoted to the study of
chaos synchronization in dynamical systems described by integer-order
differential equations \cite{Chen,Ott,Pec}. Different types of synchronization
have been proposed in the literature for continuous--time systems
\cite{Ouannas1,Ouannas2,Ouannas3,Ouannas4} as well as discrete-time
\cite{Ouannas6,Ouannas7,Ouannas5,Ouannas8}. Recently, a lot of attention has
been paid to dynamical systems described by fractional-order differential
equations \cite{Bashir,Kaslik,Tisdell}. Research studies have shown that
fractional-order systems, as generalizations of the more well--known
integer-order systems, may also have complex dynamics such as chaos and
bifurcation \cite{Cafagna+Grassi,Sayed,Rajagopal}. Some recent studies such as
\cite{LiYan,YanLi} have also shown that chaotic fractional-order systems can
be synchronized. However, since the subject is still relatively new, fewer
synchronization types have been introduced for fractional-order systems
compared to integer-order ones. It is important to note that most of the
approaches available in the literature are related to the synchronization of
identical fractional-order systems \cite{Zhang+Chen+Li+Kurths}. Very few
methods for synchronizing non-identical fractional-order chaotic systems
have been established, see \cite{Ouannas9}.

When studying the synchronization of chaotic systems, an interesting
phenomenon that may occur is the coexistence of several synchronization types.
In fact, the coexistence of these types between different dimensional chaotic
(hyperchaotic) systems remains entirely unexplored. Perhaps the most relevant
studies dedicated to the coexistence of synchronization types between two
chaotic systems include:

\begin{itemize}
\item \cite{Ouannas11}: the approach developed in this study proves rigorously
the coexistence of some synchronization types between discrete--time chaotic
(hyperchaotic) systems.

\item \cite{Ouannas12}: this study proposes two synchronization schemes of
coexistence for integer-order chaotic systems.

\item \cite{Ouannas13}: a robust method is applied to study the coexistence of
two generalized types of synchronization in fractional chaotic systems with
different dimensions. The coexistence of synchronization types can be used to
enhance the security in communications and chaotic encryption schemes.
\end{itemize}

This article investigates the coexistence of various synchronization types
between fractional chaotic (hyperchaotic) systems with different dimensions.
In particular, the paper shows that complete synchronization (CS)
\cite{Maheri+Arifin}, anti-synchronization (AS) \cite{Al-Sawalha-Al-Sawalha}
and inverse full state hybrid function projective synchronization (IFSHFPS)
\cite{Ouannas14} coexist between a three-dimensional fractional-order master
system and a four-dimensional fractional-order slave system. By exploiting
the stability theory of fractional linear systems, the coexistence of CS, AS
and IFSHFPS between two incommensurate fractional-order systems with
different dimensions is proved. Additionally, by using a fractional Lyapunov
approach, the coexistence of CS, AS and IFSHFPS is illustrated when the slave
system is of the commensurate fractional-order type. Numerical examples are
used to confirm the capability of the proposed approach in successfully
achieving the coexistence of these synchronization types in the commensurate
and incommensurate cases.

The paper is organized as follows: Section 2 lists some preliminaries relating
to fractional calculus and the stability of fractional systems. In Section 3,
the coexistence of CS, AS and IFSHFPS in fractional-order systems is
formulated. The main results of the study are presented in Section 4, followed
by some numerical examples in Section 5 that confirm the formulated problem. A
summary of the conclusions is, then, given in the last section.

\section{Preliminaries} \label{sec2}

\begin{definition}[\cite{Sam}] \label{def1}
The Riemann-Liouville fractional integral operator of
order $p>0$ of the function $f(t)$ is defined as
\begin{equation}
J^{p}f(t)=\frac{1}{\Gamma(p)}\int_{0}^{t}(t-\tau)^{p-1}f(\tau)d\tau,\quad
t>0, \label{eq1}
\end{equation}
where $\Gamma$ denotes Gamma function.
\end{definition}

\begin{definition}[\cite{Caputo}]\label{def2}
The Caputo fractional derivative of $f(t)$ is defined as
\begin{equation}
D_t^{p}f(t)=J^{m-p}\Big(\frac{d^{m}}{dt^{m}}f(t)\Big)
=\frac{1}{\Gamma(m-p)}\int_{0}^{t}\frac{f^{(m)}(\tau)}{(t-\tau)^{p-m+1}}d\tau,
\label{eq2}
\end{equation}
for $m-1<p\leq m$, $m\in \mathbf{N}$, $t>0$.
\end{definition}

\begin{lemma}[\cite{Miller}]\label{lemma1}
The Laplace transform of the Caputo fractional
derivative rule reads
\begin{equation}
\mathbf{L}(D_t^{p}f(t)  )  =s^{p}\mathbf{F}(s)
-\sum_{k=0}^{n-1}s^{p-k-1}f^{(k)  }(0),\quad (p>0,\;  n-1<p\leq n). \label{eq3}
\end{equation}
Particularly, when $0<p\leq1$, we have
\begin{equation}
\mathbf{L}(D_t^{p}f(t)  )  =s^{p}\mathbf{F}(s)  -s^{p-1}f(0)  . \label{eq4}
\end{equation}
\end{lemma}

\begin{lemma}[\cite{Pod}]\label{lemma2}
The Laplace transform of the Riemann-Liouville fractional integral rule satisfies
\begin{equation}
\mathbf{L}(J^{q}f(t))=s^{-q}\mathbf{F}(s),\quad (q>0). \label{eq5}
\end{equation}
\end{lemma}

\begin{lemma}[\cite{Matignon}] \label{lemma3}
The fractional-order linear system
\begin{equation}
D_t^{p_{i}}x_{i}(t)  =\sum_{j=1}^{n}a_{ij}x_{j}(t) ,\quad
 i=1,2,\dots ,n, \label{eq6}
\end{equation}
is asymptotically stable if all roots $\lambda$ of the characteristic
equation
\begin{equation}
\det(\operatorname{diag}(\lambda^{Mp_1},\lambda^{Mp_2},\dots ,\lambda
^{Mp_{n}})  -A)  =0, \label{eq7}
\end{equation}
satisfy $| \arg (\lambda )  | >\frac{\pi}{2M}$,
where $A=(a_{ij})  $ and $M$ is the least common multiple of the
denominators of $p_{i}$'s.
\end{lemma}

\begin{lemma}[\cite{Diyi-Chen}]\label{lemma4}
The trivial solution of the fractional order system
\begin{equation}
D_t^{p}X(t)  =F(X(t) )  , \label{eq8}
\end{equation}
where $X(t)  =(x_{i}(t)  )_{1\leq i\leq n},$ $p$ is
a rational number between $0$ and $1$, and $F:\mathbf{R}^{n}\to \mathbf{R}^{n}$
is asymptotically stable if there exists a positive definite
Lyapunov function $V(X(t)  )  $ such that
$D_t^{p}V(X(t)  )  <0$ for all $t>0$.
\end{lemma}

\begin{lemma}[\cite{Aguila-Camacho}]\label{lemma5}
For all $X(t)\in \mathbf{R}^{n}$, all $p\in ]0,1]  $ and all $t>0$,
\begin{equation}
\frac{1}{2}D_t^{p}(X^{T}(t)X(t))  \leq X^{T}(t)D_t^{p}(X(t))  . \label{eq9}
\end{equation}
\end{lemma}

\section{Problem formulation} \label{sec3}

We consider the master system given by
\begin{equation}
\begin{gathered}
D_t^{p_1}x_1(t)  =f_1(X(t)  ),\\
D_t^{p_2}x_2(t)  =f_2(X(t)  ),\\
D_t^{p_3}x_3(t)  =f_3(X(t)  ),
\end{gathered}\label{eq10}
\end{equation}
where $X(t)  =(x_1(t)$, $x_2(t)$, $x_3(t)  )  ^{T}$ is the state vector of the
master system \eqref{eq10}, $f_{i}:\mathbf{R}^{n}\to \mathbf{R}$,
$0<p_{i}<1$, and $D_t^{p_{i}}$ is the Caputo fractional derivative of order
$p_{i}$ for $i=1,2,3$. Also, consider the slave system defined as
\begin{equation}
\begin{gathered}
D_t^{q_1}y_1(t)  =\sum_{j=1}^{4}b_{1j}y_{j}(t)  +g_1(Y(t)  )  +u_1,\\
D_t^{q_2}y_2(t)  =\sum_{j=1}^{4}b_{2j}y_{j}(t)  +g_2(Y(t)  )  +u_2,\\
D_t^{q_3}y_3(t)  =\sum_{j=1}^{4}b_{3j}y_{j}(t)  +g_3(Y(t)  )  +u_3,\\
D_t^{q_4}y_4(t)  =\sum_{j=1}^{4}b_{4j}y_{j}(t)  +g_4(Y(t)  )  +u_4,
\end{gathered}\label{eq11}
\end{equation}
where $Y(t)=(y_1(t)$, $y_2(t)$, $y_3(t)$, $y_4(t))  ^{T}$ is the slave
system's state vector, $(b_{ij})  \in \mathbf{R}^{4\times4}$,
$g_{i}:\mathbf{R}^{n}\to \mathbf{R}$, $i=1,2,3,4$ are nonlinear functions,
$0<q_{i}<1,$ $D_t^{q_{i}}$ is the Caputo fractional derivative of order
$q_{i}$, and $u_{i}$, $i=1,2,3,4$ are controllers to be designed. Based on the
master--slave synchronizing system described by \eqref{eq10} and \eqref{eq11},
the following definition for the coexistence of different synchronization
types can be stated.

\begin{definition}\label{def3} \rm
Complete synchronization (CS), anti--synchronization (AS) and
inverse full state hybrid function projective synchronization (IFSHFPS)
co--exist in the synchronization of the master system \eqref{eq10} and the
slave system \eqref{eq11} if there exist controllers $u_{i}$ $(1\leq i\leq4)$
and given differentiable functions $\alpha_{i}(t)  $ $(1\leq
i\leq4)$ such that the synchronization errors:
\begin{equation}
\begin{gathered}
e_1(t)  =y_1(t)  -x_1(t)  ,\\
e_2(t)  =y_2(t)  +x_2(t)  ,\\
e_3(t)  =\sum_{j=1}^{4}\alpha_{j}(t)  y_{j}(t)  -x_3(t)  ,
\end{gathered}\label{eq12}
\end{equation}
satisfy
\begin{equation}
\lim_{t\to \infty}e_{i}(t)  =0, \label{eq12.1}
\end{equation}
for $i=1,2,3$.
\end{definition}

Before presenting the main result of this study, let us start by
rewriting the synchronization error problem \eqref{eq12}.
The system can be differentiated to yield
\begin{equation}
\begin{gathered}
D_t^{q_1}e_1(t)  =D_t^{q_1}y_1(t)-D_t^{q_1}x_1(t)  ,\\
D_t^{q_2}e_2(t)  =D_t^{q_2}y_2(t)+D_t^{q_2}x_2(t)  ,\\
\dot{e}_3(t)  =\sum_{j=1}^{4}\dot{\alpha}_{j}(t)
y_{j}(t)  +\sum_{j=1}^{4}\alpha_{j}(t)  \dot{y}_{j}(t)  -\dot{x}_3(t)  .
\end{gathered}\label{eq13}
\end{equation}
This can, then, be divided in two subsystems as follows
\begin{equation}
(D_t^{q_1}e_1(t)  ,D_t^{q_2}e_2(t)  )  ^{T}
=(B-C)  (e_1(t),e_2(t)  )  ^{T}+(u_1,u_2)  ^{T}+(R_1,R_2)  ^{T}, \label{eq14}
\end{equation}
and
\begin{equation}
\dot{e}_3(t)  =\alpha_3(t)  \dot{y}_3(t)  +R_3, \label{eq15}
\end{equation}
where $B=(b_{ij})  _{1\leq i;j\leq2}$, $C=(c_{ij})_{1\leq i;j\leq2}$
is a control matrix to be selected and
\begin{equation}
\begin{gathered}
\begin{aligned}
R_1&=(c_{11}-b_{11})  e_1(t)  +(c_{12}-b_{12})  e_2(t)
 +\sum_{j=1}^{4}b_{1j}y_{j}(t) \\
&\quad  +g_1(Y(t))  -D_t^{q_1}x_1(t)  ,
\end{aligned}\\
\begin{aligned}
R_2&=(c_{21}-b_{21})  e_1(t)  +(c_{22}-b_{22})  e_2(t)  
 +\sum_{j=1}^{4}b_{2j}y_{j}(t) \\
&\quad  +g_2(Y(t)  )  -D_t^{q_2}x_2(t)  ,
\end{aligned} \\
R_3=\sum_{j=1}^{4}\dot{\alpha}_{j}(t)  y_{j}(t)
+\sum_{\substack{j=1\\j\neq3}}^{4}\alpha_{j}(t)  \dot{y}_{j}(t)  -\dot{x}_3(t)  .
\end{gathered} \label{eq16}
\end{equation}


\section{Coexistence of synchronization types} \label{sec4}

In this section, we show that three different synchronization types can
coexist between the proposed systems \eqref{eq10} and \eqref{eq11} subject to
some conditions. In order to achieve synchronization between the master and
slave systems, we assume that $\alpha_3(t)  \neq0$ for all
$t\geq0$. Hence, we may now formulate the following theorem.

\begin{theorem} \label{theo1}
CS, AS and IFSHFPS coexist between the master system \eqref{eq10}
and the slave system \eqref{eq11} under the following conditions:
\begin{itemize}
\item[(i)] $\begin{pmatrix}
u_1\\
u_2
\end{pmatrix}
=-\begin{pmatrix}
R_1\\
R_2
\end{pmatrix}$,
\[
 u_3=-\sum_{i=1}^{4}b_{3j}y_{j}(t)-g_3
(Y(t))+J^{q_3}\big[  \frac{1}{\alpha_3(t)  }((
b_{33}-c)  e_3(t)  -R_3)  \big],
\]
 and $u_4=0$.

\item[(ii)] All roots of
\[
\det (\operatorname*{diag}(\lambda^{Mq_1},\lambda^{Mq_2})  +C-B)  =0
\]
satisfy
\[
| \arg (\lambda )  | >\frac{\pi}{2M},
\]
where $M$ is the least common multiple of the denominators of $q_1$ and
$q_2$.

\item[(iii)] The control constant $c$ is chosen such that
$b_{33}-c<0$.
\end{itemize}
\end{theorem}

\begin{proof}
To prove Theorem \ref{theo1}, we need to show that equations
\eqref{eq12.1} are satisfied. First of all, using (i), the error subsystem
\eqref{eq14} can be written in the form
\begin{equation}
D_t^{q}\hat{e}(t)  =(B-C)  \hat{e}(t)  ,\label{eq17}
\end{equation}
where
\[
D_t^{q}\hat{e}(t)  =(D_t^{q_2}e_1(t),D_t^{q_1}e_2(t)  )  ^{T}.
\]
If the feedback gain matrix $C$ is chosen according to (ii), then based on
Lemma \ref{lemma3}, we conclude that
\begin{equation}
\lim_{t\to +\infty}e_1(t)  =\lim_{t\to+\infty}e_2(t)  =0.\label{eq18}
\end{equation}

As for the third error, we use the controller $u_3$ to obtain the following
description for state $y_3(t) $
\begin{equation}
D_t^{q_3}y_3(t)  =J^{q_3}\big[  \frac{1}{\alpha
_3(t)  }((b_{33}-c)  e_3(t)-R_3)  \big]  . \label{eq19}
\end{equation}


Applying the Laplace transform to \eqref{eq19} and letting
\begin{equation}
\mathbf{F}(s)=\mathbf{L}(y_3(t)  )  , \label{eq20}
\end{equation}
we obtain,
\begin{equation}
s^{q_3}\mathbf{F}(s)-s^{q_3-1}y_3(0)
=s^{q_3-1}\mathbf{L}(\frac{1}{\alpha_3(t)  }((b_{33}-c)  e_3(t)  -R_3)  )  .
\label{eq21}
\end{equation}
Multiplying both sides of \eqref{eq21} by $s^{1-q_3}$ and applying the
inverse Laplace transform yields the equation
\begin{equation}
\dot{y}_3(t)  =\frac{1}{\alpha_3(t)  }((b_{33}-c)  e_3(t)  -R_3)  . \label{eq22}
\end{equation}
From  \eqref{eq22} and \eqref{eq15}, the dynamics of $e_3(t)  $ can be given by
\begin{equation}
\dot{e}_3(t)=(b_{33}-c)  e_3(t)  . \label{eq23}
\end{equation}
If $c$ is selected according to (ii), we obtain
\begin{equation}
\lim_{t\to \infty}e_3(t)=0. \label{eq24}
\end{equation}
Finally, from  \eqref{eq18} and \eqref{eq24}, we conclude that the master
system \eqref{eq10} and the slave system \eqref{eq11} are globally synchronized.
\end{proof}

The conditions derived in Theorem \ref{theo1} can be simplified in the case
where $q_1=q_2=q$. The following proposition shows the simplification.

\begin{proposition}\label{prop1}
Subject to $q_1=q_2=q$, condition (ii) of Theorem \ref{theo1}
may be replaced by the following condition:
\begin{itemize}
\item[(ii)] The control matrix $C$
is selected such that $B-C$ is a negative definite matrix.
\end{itemize}
\end{proposition}

\begin{proof}
If a Lyapunov function candidate is chosen as
$V(\hat{e}(t)  )  =\frac{1}{2}\hat{e}^{T}(t)  \hat{e}(t)  $,
then the time Caputo fractional derivative of order $q$ of $V$
along the trajectory of the system \eqref{eq17} may be stated as
\begin{equation}
D_t^{q}V(\hat{e}(t)  )
=D_t^{q}\big(\frac {1}{2}\hat{e}^{T}(t)  \hat{e}(t)  \big)  .
\label{eq25}
\end{equation}
Using Lemma \ref{lemma5} along with \eqref{eq25}, we obtain
\[
D_t^{q}V(e(t)  )   \leq \hat{e}^{T}(t)  D_t^{q}\hat{e}(t)
 =\hat{e}^{T}(t)(B-C)  \hat{e}(t).
\]
If we select the matrix $C$ such that $B-C$ is negative definite, we obtain
\[
D_t^{q}V(\hat{e}(t)  )  <0.
\]
From Lemma \ref{lemma4}, the zero solution of system \eqref{eq17} is globally
asymptotically stable, i.e
\begin{equation}
\lim_{t\to +\infty}e_1(t)  =\lim_{t\to
+\infty}e_2(t)  =0. \label{eq26}
\end{equation}
\end{proof}

\section{Numerical examples} \label{sec5}

In this section, we use numerical simulations to validate the theoretical
synchronization results proposed in the previous section given some examples
of nonlinear chaotic fractional systems.

\subsection*{Case I: $q_1\neq q_2$}
Consider as master the fractional version of the modified coupled dynamos
system proposed in \cite{Wang+Wang} and given by
\begin{equation}
\begin{gathered}
D^{p_1}x_1=-\alpha x_1+(x_3+\beta )  x_2,\\
D^{p_2}x_2=-\alpha x_2+(x_3-\beta )  x_1,\\
D^{p_3}x_3=x_3-x_1x_2.
\end{gathered}\label{eq27}
\end{equation}
System \eqref{eq27} can  exhibit chaotic behaviors when
$(p_1,p_2,p_3)=(0.9,0.93,0.96)$ and $(\alpha,\beta)=(2,1)$. Attractors
of the master system \eqref{eq27} are shown in Figure \ref{Fig1}.

\begin{figure}[htb]
\begin{center}
 \includegraphics[width=0.9\textwidth]{fig1} 
\end{center}
\caption{Chaotic attractors in 3-D of the master system \eqref{eq27} when 
$(p_1,p_2,p_3)=(0.9,0.93,0.96)$ and $(\alpha,\beta)=(2,1)$.}
\label{Fig1}
\end{figure}

Let us also consider the slave system 
\begin{equation}
\begin{gathered}
D^{q_1}y_1=\alpha (y_2-y_1)  +y_4+u_1,\\
D^{q_2}y_2=\gamma y_1-y_2-y_1y_3+u_2,\\
D^{q_3}y_3=-\beta y_3+y_1y_2+u_3,\\
D^{q_4}y_4=\delta y_4+y_2y_3+u_4,
\end{gathered} \label{eq28}
\end{equation}
where the vector controller is
\[
U=(u_1,u_2,u_3,u_4)  ^{T}\,.
\]
We observe that for 
$(u_1,u_2,u_3 ,u_4)  =(0,0,0,0)$, 
$(q_1,q_2,q_3,q_4)  =(0.94,0.96,0.97,0.99)  $
and $(\alpha,\beta,\gamma,\delta )  =(10,\frac{8}{3},28,-1)  $,
system \eqref{eq28} exhibits a hyperchaotic behavior, see \cite{Li+Wang+Luo}.
Attractors of the uncontrolled system \eqref{eq28} are shown in Figure
\ref{Fig2}.

\begin{figure}[htb]
\begin{center}
 \includegraphics[width=0.9\textwidth]{fig2} 
\end{center}
\caption{Hyperchaotic attractors in 3-D of the slave system \eqref{eq28}
 when $(u_1,u_2,u_3,u_4)=(0,0,0,0)$, 
$(q_1,q_2,q_3 ,q_4)  =(0.94,0.96,0.97,0.99) $ and 
$(\alpha,\beta,\gamma,\delta )  =(10,\frac{8}{3},28,-1)$.}
\label{Fig2}
\end{figure}

According to our approach presented in the previous sections, the error system
between the master \eqref{eq27} and slave \eqref{eq28} is defined as
\begin{equation}
\begin{gathered}
e_1=y_1-x_1,\\
e_2=y_2+x_2,\\
e_3=\alpha_1(t)  y_1+\alpha_2(t)
y_2+\alpha_3(t)  y_3+\alpha_4(t)
y_4-x_3,
\end{gathered}\label{eq29}
\end{equation}
where $\alpha_1(t)  =1$, 
$\alpha_2(t)  =\frac{1}{t^{2}+1}$, $\alpha_3(t)  =\exp (-t)  $ and
$\alpha_4(t)  =\sin t$. Using the notations defined in Section
\ref{sec3} above, we can write
\[
B=\begin{pmatrix}
-10 & 10\\
28 & -1
\end{pmatrix},\quad 
C=\begin{pmatrix}
0 & 10\\
28 & 0
\end{pmatrix},
\]
$b_{33}=-8/3$  and $c=0$.
According to Theorem \ref{theo1}, the controllers set of controllers $(
u_1,u_2,u_3,u_4)  $ may be designed as
\begin{equation}
\begin{gathered}
u_1=-10e_1+\alpha (y_2-y_1)  -D_t^{0.94}x_1,\\
u_2=-e_2+\gamma y_1-y_2-y_1y_3-D_t^{0.96}x_2,\\
\begin{aligned}
u_3&=\beta y_3-y_1y_2+\frac{1}{3}J^{0.97}
\Big(-\frac{8}{3}
e_3(t)  -\frac{1}{(t+1)  ^{2}}y_2+\exp (
-t)  y_3 \\
&\quad -(\cos t)  y_4-y_1-(\frac{1}{t+1})  y_2-(\sin t)  y_4+\dot{x}_3\Big),\\
u_4=0.
\end{aligned}\label{eq30}
\end{gathered}
\end{equation}

The roots of equation
\[
\det (\operatorname*{diag}(\lambda^{0.94M},\lambda^{0.96M})  +C-B)  =0
\]
are
\begin{gather*}
\lambda_1=10^{\frac{1}{0.94M}}(\cos \frac{\pi}{0.94M}+\mathbf{i}
\sin \frac{\pi}{0.94M}),\\
\lambda_2=\cos \frac{\pi}{0.96M}+\mathbf{i}\sin \frac{\pi}{0.96M},
\end{gather*}
where $M$ is the least common multiple of the denominators of $0.94$ and
$0.96$. It is easy to show that $| \arg (\lambda_{i})
| >\frac{\pi}{2M}$ for $i=1,2$. Hence, the conditions of Theorem
\ref{theo1} are satisfied and, consequently, systems \eqref{eq27} and
\eqref{eq28} are globally synchronized. The error system can be summarized by
the two subsystems
\begin{equation}
\begin{gathered}
D^{0.94}e_1=-10e_1,\\
D^{0.96}e_2=-e_2,
\end{gathered} \label{eq31}
\end{equation}
and
\begin{equation}
\dot{e}_3=-\frac{8}{3}e_3. \label{eq32}
\end{equation}
Fractional Euler integration and fourth order Runge-Kutta integration methods
have been used to solve systems \eqref{eq31} and \eqref{eq32}, respectively.
Time evolution of the errors $e_1,e_2$ and $e_3$ are shown in Figures
\ref{Fig3} and \ref{Fig4}, respectively.

\begin{figure}[htb]
\begin{center}
 \includegraphics[width=0.7\textwidth]{fig3} 
\end{center}
\caption{Time series of the synchronized error signals $e_1$ and 
$e_2$ between the master system \eqref{eq27} and the slave system \eqref{eq28}.}
\label{Fig3}
\end{figure}

\begin{figure}[htb]
\begin{center}
\includegraphics[width=0.7\textwidth]{fig4} 
\end{center}
 \caption{Time evolution of the error $e_3$ between the master system 
\eqref{eq27} and the slave system \eqref{eq28}.}
\label{Fig4}
\end{figure}

\subsection*{Case II: $q_1=q_2=q$}

Now, let us consider as master the fractional-order Liu system with the
hyperchaotic fractional-order Lorenz system as its slave. The master system
is 
\begin{equation}
\begin{gathered}
D^{p_1}x_1=a(x_2-x_1)  ,\\
D^{p_2}x_2=bx_1-x_1x_3,\\
D^{p_3}x_3=-cx_3+4x_1^{2}.
\end{gathered}\label{eq33}
\end{equation}
This system exhibits chaotic behavior when 
$(p_1,p_2,p_3)=(0.93,0.94,0.95)  $ and 
$(a,b,c)  =(10,40,2.5)  $ \cite{Liu+Liu+Liu+Liu}. The attractors for \eqref{eq33}
are shown in Figure \ref{Fig5}.

\begin{figure}[htb]
\begin{center}
\includegraphics[width=0.9\textwidth]{fig5} 
\end{center}
\caption{Chaotic attractors in 3-D of the master system \eqref{eq33} when 
$(p_1,p_2,p_3)  =(0.93,0.94,0.95)  $ and $(a,b,c)  =(10,40,2.5)$.}
\label{Fig5}
\end{figure}

The slave system is
\begin{equation}
\begin{gathered}
D^{q_1}y_1=0.56y_1-y_2+u_1,\\
D^{q_2}y_2=y_1-0.1y_2y_3^{2}+u_2,\\
D^{q_3}y_3=4y_2-y_3-6y_4+u_3,\\
D^{q_4}y_4=0.5y_3+0.8y_4+u_4,
\end{gathered} \label{eq34}
\end{equation}
where $u_1,u_2,u_3,u_4$ are the synchronization controllers. This
system, as illustrated in \cite{Zhou+Wei+Cheng}, exhibits a hyperchaotic
behavior when $(u_1,u_2,u_3,u_4)=(0,0,0,0)  $ and
$(q_1,q_2,q_3,q_4)  =(0.98,0.98,0.95,0.95)$. 
The attractors of \eqref{eq34} are shown in Figure \ref{Fig6}.

\begin{figure}[htb]
\begin{center}
\includegraphics[width=0.9\textwidth]{fig6} 
\end{center}
\caption{Hyperchaotic attractors in 3-D of the slave system \eqref{eq34}
 when when $(u_1,u_2,u_3,u_4)=(0,0,0,0)  $ and 
$(q_1,q_2,q_3,q_4)  =(0.98,0.98,0.95,0.95)$.}
\label{Fig6}
\end{figure}

The error system between the master \eqref{eq33} and slave \eqref{eq34} is 
\begin{equation}
\begin{array}
[c]{l}
e_1=y_1-x_1,\\
e_2=y_2+x_2,\\
e_3=\sin (t)  y_1+\cos (t)  y_2+\frac{1}
{t^{2}+1}y_3+y_4-x_3.
\end{array}
\label{eq35}
\end{equation}
In this case, based on the notation presented in the Section \ref{sec3}, we
 write
\[
B=\begin{pmatrix}
0.56 & -1\\
1 & 0
\end{pmatrix}, \quad
C=\begin{pmatrix}
1 & -1\\
1 & 1
\end{pmatrix},
\]
$b_{33}=-1$  and $c=0$,
and using Theorem \ref{theo1} and Proposition \ref{prop1}, the controllers can
constructed as
\begin{equation}
\begin{gathered}
u_1=-0.44e_1-0.56y_1+y_2-D_t^{0.98}x_1,\\
u_2=-e_2+-y_1+0.1y_2y_3^{2}-D_t^{0.98}x_2,\\
\begin{aligned}
u_3&=-4y_2+y_3+6y_4+J^{0.97}(t^{2}+1)  (-e_3
-\cos (t)  y_1+\sin (t)  y_2) \\
&\quad +J^{0.97}(t^{2}+1)  \Big(-\frac{2t}{(
t^{2}+1)  ^{2}}y_3-\sin (t)  y_1-\cos (t)
y_2-y_4+\dot{x}_3\Big)  ,
\end{aligned}\\
u_4=0.
\end{gathered}\label{eq36}
\end{equation}
We can show that $B-C$ is a negative definite matrix, which fulfills the
condition of Proposition \ref{prop1}. Therefore, systems \eqref{eq36} and
\eqref{eq37} are globally synchronized. The error systems is
\begin{equation}
\begin{gathered}
D^{0.98}e_1=-0.44e_1,\\
D^{0.98}e_2=-e_2,
\end{gathered} \label{eq37}
\end{equation}
and
\begin{equation}
\dot{e}_3=-e_3. \label{eq38}
\end{equation}
Again, similar to the previous example, the Fractional Euler and fourth order
Runge-Kutta integration methods have been used to solve systems \eqref{eq37}
and \eqref{eq38}, respectively. The time evolutions of $e_1,e_2$ and
$e_3$ are shown in Figures \ref{Fig7} and \ref{Fig8}.

\begin{figure}[htb]
\begin{center}
\includegraphics[width=0.7\textwidth]{fig7}
\end{center}
\caption{Time series of the synchronized error signals $e_1$ and $e_2$ 
between the master system \eqref{eq33} and the slave system \eqref{eq34}.}
\label{Fig7}
\end{figure}

\begin{figure}[htb]
\begin{center}
\includegraphics[width=0.7\textwidth]{fig8}
\end{center}
\caption{Time evolution of the error $e_3$ between the master system
 \eqref{eq33} and the slave system \eqref{eq34}.}
\label{Fig8}
\end{figure}

\subsection*{Conclusions}

This paper has proposed a new method to analyze the coexistence problem of
some fractional synchronization types. In particular, the approach developed
in this paper has proven the coexistence of complete synchronization (CS),
anti--synchronization(AS) and inverse full state hybrid function projective
synchronization (IFSHFPS) between a three-dimensional fractional-order
master system and a four-dimensional fractional-order slave system. It has
been shown that the approach presents the remarkable feature of being both
rigorous and applicable to a wide class of commensurate and incommensurate
systems with different dimensions and orders. The numerical examples reported
in the paper have clearly highlighted the capability of the proposed approach
in successfully achieving the coexistence of these synchronization types
between chaotic and hyperchaotic systems of different dimensions for both
commensurate and incommensurate fractional-order cases.

\begin{thebibliography}{99}


\bibitem{Aguila-Camacho} N. Aguila-Camacho, M. Duarte-Mermoud, J. Gallegos;
\emph{Lyapunov functions for fractional order systems}, Commun. Nonlin. Sci. Num.
Sim., Vol. 19 (9) (2014), pp. 2951--2957.

\bibitem{Bashir} A. Bashir, N. Sotirisk;
\emph{Existence and uniqueness of solutions
for caputo-hadamard sequential fractional order neutral functional
differential equations}, Elec. J. Diff. Eqs., Vol. 2017 (36)  (2017), pp. 1--11.

\bibitem{Cafagna+Grassi} D. Cafagna, G. Grassi;
\emph{Fractional-order systems without equilibria: The first example of 
hyperchaos and its application to
synchronization}, Chinese Physics B, Vol. 24 (8)  (2015),  080502-6.

\bibitem{Caputo} M. Caputo;
\emph{Linear models of dissipation whose $Q$ is almost
frequency independent--II}, 
Geophys J. R. Astron. Soc., Vol. 13 (1967), pp. 529--539.

\bibitem{Chen} G. Chen, X. Yu;
\emph{Chaos Control: Theory and Applications},
Springer-Verlag, Berlin, Germany, 2003.

\bibitem{Diyi-Chen} D. Chen, R. Zhang, X. Liu, X. Ma;
\emph{Fractional order Lyapunov stability theorem and its applications 
in synchronization of complex dynamical networks}, 
Commun. Nonlin. Sci. Num. Sim., Vol. 19 (12) (2014), pp. 4105--4121.

\bibitem{Sayed} A. M. A. El-Sayed, H. M. Nour, A. Elsaid, A. E. Matouk, A. Elsonbaty;
\emph{Dynamical behaviors, circuit realization, chaos control, and
synchronization of a new fractional order hyperchaotic system}, Applied Math.
Modeling, Vol. 40 (2016), pp. 3516--3534.

\bibitem{Kaslik} E. Kaslik, S. Sivasundaram;
\emph{Non-existence of periodic solutions in fractional-order dynamical 
systems and a remarkable difference between integer and fractional-order 
derivatives of periodic functions},
Nonlinear Analysis: Real World Applications, Vol. 13 (3), pp. 1489--1497.

\bibitem{LiYan} C. Li, J. Yan;
\emph{The synchronization of three fractional differential systems}, 
Chaos Solitons Fractals 22 (3) (2007), pp. 751-757.

\bibitem{Liu+Liu+Liu+Liu} C. Liu, T. Liu, L. Liu, K. Liu;
\emph{A new chaotic attractor}, Chaos Solitons Fractals, 
Vol. 22 (5) (2004), pp. 1031--1038.

\bibitem{Maheri+Arifin} M. Maheri, N. Md. Arifin;
\emph{Synchronization of two different fractional-order chaotic systems 
with unknown parameters using a robust adaptive nonlinear controller}, 
Nonlin. Dyn., Vol. 85 (2) (2016), pp. 825--838.

\bibitem{Matignon} D. Matignon;
\emph{Stability results of fractional differential
equations with applications to control processing}, 
IMACS, IEEE-SMC, France, 1996.

\bibitem{Miller} K. S. Miller, B. Ross;
\emph{An Introduction to The Fractional Calculus and Fractional 
Differential Equations}, John Wiley and Sons, New York, USA, 1993.

\bibitem{Al-Sawalha-Al-Sawalha} M. Mossa Al-Sawalha, Ayman Al-Sawalha;
\emph{Anti-synchronization of fractional order chaotic and hyperchaotic systems
with fully unknown parameters using modified adaptive control}. Open Phys.,
Vol. 14 (2016), pp. 304--313.

\bibitem{Ott} E. Ott, C. Grebogi, J. A. Yorke;
\emph{Controlling chaos}, Phys. Rev. Lett. 64 (1990), pp. 1196--1199.


\bibitem{Ouannas6} A. Ouannas;
\emph{A new generalized-type of synchronization for
discrete chaotic dynamical systems}, Journal of Computational and Nonlinear
Dynamics, Vol. 10 (6) (2015),  061019-5.

\bibitem{Ouannas1} A. Ouannas;
\emph{On inverse generalized synchronization of
continuous chaotic dynamical systems}, International Journal of Applied and
Computational Mathematics, Vol. 2 (1) (2016), pp. 1--11.

\bibitem{Ouannas2} A. Ouannas, M. M. Al-sawalha;
\emph{Synchronization between different dimensional chaotic systems 
using two scaling matrices}, Optik Int. J. for Light and Electron Optics, 
Vol. 127 (2016), pp. 959--963.

\bibitem{Ouannas3} A. Ouannas, M. M. Al-sawalha;
\emph{On $\Lambda-\phi$ generalized synchronization of chaotic dynamical systems 
in continuous-time}, Euro. Phys. J. Special Topics, Vol. 225 (1) (2016), 
pp. 187-196.

\bibitem{Ouannas9} A. Ouannas, M. M. Al-Sawalha, T. Ziar;
\emph{Fractional chaos synchronization schemes for different dimensional 
systems with non-identical fractional-orders via two scaling matrices}, 
Optik, Vol. 127 (20) (2016), pp. 8410--8418.

\bibitem{Ouannas12} A. Ouannas, A. T. Azar, V. Sundarapandian;
\emph{New hybrid synchronization schemes based on coexistence of various
types of synchronization between master-slave hyperchaotic systems},
Int. J. Comp. Apps. in Tech., Vol. 55 (2) (2017), pp. 112--120.

\bibitem{Ouannas13} A. Ouannas, A. T. Azar, S. Vaidyanathan;
\emph{A robust method for new fractional hybrid chaos synchronization},
Mathe. Meth. App. Sci., Vol. 40 (5) (2017), pp. 1804--1812.

\bibitem{Ouannas4} A. Ouannas, A. T. Azar, T. Ziar;
\emph{On inverse full state hybrid function projective synchronization
for continuous-time chaotic dynamical systems with arbitrary dimensions},
Diff. Eqs. and Dyn. Sys., pp. 1--14.

\bibitem{Ouannas7} A. Ouannas, G. Grassi;
\emph{Inverse Full State Hybrid Projective Synchronization for Chaotic
Maps with Different Dimensions}, Chinese Physics B,
Vol. 25 (9) (2016),  090503-6.

\bibitem{Ouannas11} A. Ouannas, G. Grassi;
\emph{A new approach to study coexistence of some synchronization types 
between chaotic maps with different dimensions},
Nonlin. Dyn. 86 (2) (2016), pp. 1319--1328.

\bibitem{Ouannas14} A. Ouannas, G. Grassi, T. Ziar, Z. Odibat;
\emph{On a Function Projective Synchronization Scheme between non-identical
Fractional-order chaotic (hyperchaotic) systems with different dimensions and orders},
Optik Int. J. Light and Electron Optics, Vol. 136 (2017), pp. 513--523.

\bibitem{Ouannas5} A. Ouannas, Z. Odibat;
\emph{Generalized synchronization of different dimensional chaotic
dynamical systems in discrete-time}, Nonlinear Dynamics, Vol. 81 (1) (2015),
 pp. 765--771.
\bibitem{Ouannas8} A. Ouannas, Z. Odibat, N. Shawagfeh;
\emph{Universal chaos synchronization control laws for general quadratic 
discrete systems}, Applied Mathematical Modeling Vol. 45 (2017), pp. 636--641.

\bibitem{Pec} L. M. Pecora, T. L. Carrol;
\emph{Synchronization in chaotic systems},
Phys. Rev. Lett. 64 (1990), pp. 821--824.

\bibitem{Zhou+Wei+Cheng} Z. Ping, W. Li-Jia, C. Xue-Feng;
\emph{A novel fractional-order hyperchaotic system and its synchronization}, 
Chin. Phys. B, Vol. 18 (7) (2009), pp. 2674--2679.

\bibitem{Pod} I. Podlubny;
\emph{Fractional Differential Equations}, Academic Press,
San Diego, CA, 1999.

\bibitem{Rajagopal} K. Rajagopal, L. Guessas, A. Karthikeyan, A. Sirinvasan, G. Adam;
\emph{Fractional order memristor no equilibrium chaotic system with its
adaptive sliding mode synchronization and genetically optimized fractional
order PID synchronization}, Complexity, Vol. 2017 (2017).

\bibitem{Sam} S. G. Samko, A. A. Klibas, O. I. Marichev;
\emph{Fractional Integrals and Derivatives: Theory and Applications}, 
Gordan and Breach, 1993.

\bibitem{Li+Wang+Luo} L. Tian-Zenga, W. Yua, L. Mao-Kanga;
\emph{Control of fractional chaotic and hyperchaotic systems based on a 
fractional order controller}, Chin. Phys. B, Vol. 23 (8) (2014),  080501-11.

\bibitem{Tisdell} C. C. Tisdell;
\emph{Basic existence and a priori bound results for
solutions to systems of boundary value problems for fractional differential
equations}, Elec. J. Diff. Eqs., Vol. 2016 (84) (2016), pp. 1--9.

\bibitem{Wang+Wang} X. Y. Wang, Y. J. He, M. J. Wang;
\emph{Chaos control of a fractional order modified coupled dynamos system},
 Nonlinear Analysis: Theory, Methods and Applications, Vol. 71 (12) (2009), 
pp. 6126--6134.

\bibitem{YanLi} J. Yan, C. Li;
\emph{On chaos synchronization of fractional differential equations},
 Chaos Solitons Fractals 32 (2) (2007), pp. 725-735.

\bibitem{Zhang+Chen+Li+Kurths} F. Zhang, G. Chen, C. Li, J. Kurths;
\emph{Chaos synchronization in fractional differential systems},
 Phil. Trans. R. Soc. 371, 2013.

\end{thebibliography}


\end{document}
