# Amplitude Modulation and Synchronization of Fractional-Order Memristor-Based Chua's Circuit

Authors: A. G. Radwan; K. Moaddy; I. Hashim
Journal: Abstract and Applied Analysis
 (2013)

---

## Abstract

This paper presents a general synchronization technique and an amplitude modulation of chaotic generators. Conventional synchronization and antisynchronization are considered a very narrow subset from the proposed technique where the scale between the output response and the input response can be controlled via control functions and this scale may be either constant (positive, negative) or time dependent. The concept of the proposed technique is based on the nonlinear control theory and Lyapunov stability theory. The nonlinear controller is designed to ensure the stability and convergence of the proposed synchronization scheme. This technique is applied on the synchronization of two identical fractional-order Chua's circuit systems with memristor. Different examples are studied numerically with different system parameters, different orders, and with five alternative cases where the scaling functions are chosen to be positive/negative and constant/dynamic which covers all possible cases from conventional synchronization to the amplitude modulation cases to validate the proposed concept.


---

## Body

## 1. Introduction

Despite that the history of fractional calculus started in the same period of time as integer calculus, the major revolutions in this area have been discovered only during the last five decades where the realization, modeling, and numerical simulations were available [1–4]. Similarly, the chaotic systems have been studied heavily during the last four decades since they play an important role in industrial applications particularly in chemical reactions [5], biological systems [6], circuit theory [7–11], control [12], and security applications [13–16]. Recently much attention has been devoted to the search for better and more efficient methods for the control or determination of a solution, approximate or analytical, of chaotic systems.Antisynchronization is a phenomenon in which the state vectors of the synchronized systems have the same amplitude but opposite signs to those of the driving system. Therefore, the sum of two signals is expected to converge to zero when antisynchronization appears. Since the discovery of antisynchronization experimentally in the context of self-synchronization, it has been applied in many different fields, such as biological and physical systems, structural engineering, and ecological models [17]. Liu et al. [18] show that either synchronization or antisynchronization can appear depending on the initial conditions of the coupled pendula. Active control method is used to study the antisynchronization for two identical and nonidentical systems [19–21].The synchronization of fractional chaotic systems has started to attract much attention and has also raised up some problems [22–24]. Recently the consistency of the improvement of models based on fractional-order differential structure has increased in reputation in the research of dynamical systems [25–27]. Yu et al. studied the synchronization of three chaotic fractional-order Lorenz systems with bidirectional coupling [28], Odibat et al. [29] investigated the chaos synchronization of two identical systems via linear control, and Bhalekar and Daftardar-Gejji [30] demonstrated that two different fractional-order chaotic systems can be synchronized using active control.Recently, Radwan et al. [31] developed a framework to obtain approximate numerical solutions of the fractional-order Chua’s circuit with memristor using a nonstandard finite difference method. The most important advantage of using fractional order is to increase the chaotic range as proven from the stability analyses of this circuit in both the W-plane and s-plane. Moreover, another circuit from Chua’s family was studied and modified with the use of memristor with detailed responses by Petráš in [32]. In this paper, the concept of dynamic scaling between the master and slave systems in the fractional-order domain is introduced. The proposed concept covers the conventional synchronization up to the amplitude modulation technique. Different examples have been studied based on the nonlinear antisynchronization of two identical fractional-order chaotic Chua’s circuits with memristor. This generalized projective synchronization is based on the nonlinear control theory and Lyapunov stability theory.The paper is organized as follows. In Section2 we provide some mathematical models of the fractional calculus theory and describe the nonstandard finite difference scheme to solve fractional differential equations and in the last subsection we present the fractional Chua’s circuit with memristor. Nonlinear antisynchronization of two identical fractional-order chaotic Chua’s circuit systems is also introduced in Section 3, while the generalized projective synchronization is discussed in Section 4. Finally, conclusions are drawn in the last section.


## 2. Fractional Order System Memristor Based Chua’s Circuit

Recently, most of the dynamical systems based on the integer-order calculus have been modified into the fractional order domain due to the extra degrees of freedom which increase the flexibility of the design and can be used to precisely fit the experimental data much better than the integer-order modeling. The Caputo fractional derivative of orderα of a continuous function f:R+→R is defined as follows:
(1)Dαf(t)≡dαf(t)dtα={1Γ(m-α)∫0tf(m)(τ)(t-τ)α-m+1dτ,m-1<α<m,dmdtmf(t),α=m,
Where m is the first integer greater than α and Γ(·) is the gamma function and is defined by
(2)Γ(z)=∫0∞e-ttz-1dt,Γ(z+1)=zΓ(z).
The Grünwald-Letnikov method of approximation for the one-dimensional fractional derivative is as follows [33]:(3a)Dαx(t)=f(t,x),(3b)Dαx(t)=limh→0h-α∑j=0t/h(-1)j(αj)x(t-jh),where α>0 and D∝ denotes the fractional derivative. Also N=[t/h], and h is the step size. Therefore, (3a) and (3b) is discretized as follows:
(4)∑j=0n+1cjαx(t-jh)=f(tn,x(tn)),n=1,2,3,…,
Where tn=nh and cjα are the Grünwald-Letnikov coefficients defined as
(5)cjα=(1-1+αj)cj-1α,c0α=h-α,j=1,2,3,….
The nonstandard discretization technique [34] is a general scheme where we replace the step size h by a function φ(h). By applying this technique and using the Grünwald-Letnikov discretization method [2, 31], it yields the following relations:
(6)xn+1=-∑j=1n+1cjαxn+1-j+f1(tn+1,xn+1)c0α1,
Where c0α1=(φ1(h))-1 are functions of the step size h=Δt, with φ1(h)=h+O(h2), when h→0. The fractional-order Chua’s circuit [31] based on the new nanodevice element which is called memristor [35–38] is shown in Figure 1 with its q-ϕ relationship. The canonical memristor-based fractional-order Chua’s circuit can be written as follows:(7a)(DαxDβyDγzDu)=(-ak-101a0-b001c00000)(xyzu),(7b)k={m1|u|<1,m2|u|>1,where a,b,c,m1, and m2 are related to the circuit parameters [31]. The effect of the memristive property appears through the parameter k (the slope of the charge with respect to the flux in the memristor) which is a piecewise-linear function with two different slopes m1 and m2 depending on the value of the flux u (Figure 1(b)). It is clear that the system has four variables which are {x,y,z, and u} and the chaotic behavior exists in a wider range as proved in [31]. The chaotic behavior has been verified by the calculation of the maximum Lyapunov exponent (MLE). Figures 1(c) and 1(d) show the chaotic responses for two different cases where the α order changes from 1.0 to 1.1. It is clear from these subfigures that the output response is very sensitive to the fractional-order parameters. Therefore, these fractional-order parameters can be used for extra degrees of freedom and also to control the behavior of the strange attractors of the Chua’s circuit.(a) The canonical Chua’s circuit with memristor, (b) the piecewise-linear flux-controlled memristor, and (c) 3D shapes in theu-y-z plane when (α,β,γ) (a) (1.0, 0.9, 0.9) and (d) (1.1, 0.9, 0.9).
(a)
(b)
(c)
(d)


## 3. The Proposed Technique

From the previous section, the fractional-order memristor-based Chua’s circuit has ten parameters which are three fractional-orders{α,β,γ}, three system parameters {a,b,c}, and the initial conditions. Figure 2 shows the general block diagram of the proposed system assuming that the fractional-orders of the two systems are identical but the other seven parameters may be different. The idea is based on modifying the equations of the second system by adding control functions {u1,u2,u3,u4} on the four equations to minimize the functions x1+s1x2, y1+s2y2, z1+s3z2, and w1+s4w2 to approach zero. In this case, the output response of the second system is given by
(8)x2(t)=-x1(t)s1(t),y2(t)=-y1(t)s2(t),z2(t)=-z1(t)s3(t),w2(t)=-w1(t)s4(t).
Then, if s1=s2=s3=s4=-1, this will be the conventional synchronization, and if s1=s2=s3=s4=1, this will be the antisynchronization as will be discussed in the next subsection. Moreover, the scaling functions si(t), i=1,2,3,4, can have different values and time-dependent functions are discussed in the next subsections.Figure 2
The proposed block diagram.


### 3.1. Nonlinear Antisynchronization of Two Identical Fractional-Order Chaotic Chua’s Circuit Systems (si=1, i=1,2,3,4)

In this subsection we observe the antisynchronization behavior in two identical fractional-order chaotic Chua’s circuit systems via nonlinear control. The drive (master) and slave (response) systems are described, respectively, by the following equations. However, the initial condition on the drive system is different from that of the response system:(9a)Dαx1(t)=a1(y1-f(w1)x1),(9b)Dβy1(t)=z1-x1,(9c)Dγz1(t)=-b1y1+c1z1,(9d)Dw1(t)=x1,(10a)Dαx2(t)=a2(y2-f(w2)x2)+u1,(10b)Dβy2(t)=z2-x2+u2,(10c)Dγz2(t)=-b2y2+c2z2+u3,(10d)Dw2(t)=x2+u4,where the control function u=[u1;u2;u3;u4]T=0. Our goal is to determine the control function to make the first derivative of V(e), that is, V(e)<0. From (9a), (9b), (9c), and (9d) to (10a), (10b), (10c), and (10d), the error equations when the system parameters {a1,b1,c1}={a2,b2,c2}={a,b,c} are given as follows:(11a)Dαex=aey-a(f(w1)x1-f(w2)x2)+u1,(11b)Dβey=ez-ex+u2,(11c)Dγez=-bey+cez+u3,(11d)Dew=ex+u4,where ex=x2+x1, ey=y2+y1, ez=z2+z1, ew=w2+w1. In order to determine the controller, let
(12)u1=u1a+u1b,u1b=-a(f(w1)x1+f(w2)x2).
Then we rewrite system (11a), (11b), (11c), and (11d) in the following form:(13a)Dαex=aey+u1b,(13b)Dβey=ez-ex+u2,(13c)Dγez=-bey+cez+u3,(13d)Dew=ex+u4.Based on the Lyapunov stability theory, when the controller satisfies the assumption V(e)=0.5eTe, a positive definite function, the first derivative of V(e), V(e)<0, and the antisynchronization of two identical Chua’s circuit systems from different initial conditions is achieved, we choose u as follows:(14a)u1b=-ex-aey,(14b)u2=-ez+ex-ey,(14c)u3=-2cez+bey,(14d)u4=-ex-eu.Theorem 1.
The nonlinear controller achieves global antisynchronization between two identical Chua’s circuit systems master (9a), (9b), (9c), and (9d) and slave (10a), (10b), (10c), and (10d).Proof.
Take a Lyapunov function for (13a), (13b), (13c), and (13d) into consideration as
(15)V(e)=0.5eTe.
We then get the first derivative of V(e):(16a)V˙=ex(aey+u1b)+ey(ez-ex+u2)+ez(-bey+cez+u3)+ew(ex+u4),(16b)V˙=-ex2-ey2-bcez2-ew2,so V˙(e)<0 is satisfied. Since V˙(e) is a negative-definite function [39], the error states tend to zero at steady state. Therefore, the states of controlled response system and drive system are globally antisynchronized asymptotically.


### 3.2. Generalized Projective Synchronization of Two Identical Chua’s Circuits

In this subsection we study the synchronization between two identical fractional-order chaotic Chua’s circuits based on the generalized projective synchronization. The unknown termsu1, u2, u3, and u4 in (11a), (11b), (11c), and (11d) are active control functions to be determined. Define the error functions as
(17)(ex(t)ey(t)ez(t)ew(t))=(x1(t)+s1(t)x2(t)y1(t)+s2(t)y2(t)z1(t)+s3(t)z2(t)w1(t)+s4(t)w2(t)).
Then by similar procedure, the control functions are given by(18a)(u1(t)u2(t)u3(t)u4(t))=((a1f(w1(t))-kx)x1(t)-a1y1(t)s1(t)+(a2f(w2(t))-kx)x2(t)-a2y2(t)x1(t)-z1(t)-kyy1(t)s2(t)+x2(t)-z2(t)-kyy2(t)b1y1(t)-c1z1(t)-kzz1(t)s3(t)+b2y2(t)-c2z2(t)-kzz2(t)-x1(t)-kww1(t)s4(t)+-x2(t)-kww2(t)). Then
(18b)(DαexDβeyDγezDew)=(-kx0000-ky0000-kz0000-kw)(exeyezew)which are decaying functions as the values of kx, ky, kz, and kw are positive.


## 3.1. Nonlinear Antisynchronization of Two Identical Fractional-Order Chaotic Chua’s Circuit Systems (si=1, i=1,2,3,4)

In this subsection we observe the antisynchronization behavior in two identical fractional-order chaotic Chua’s circuit systems via nonlinear control. The drive (master) and slave (response) systems are described, respectively, by the following equations. However, the initial condition on the drive system is different from that of the response system:(9a)Dαx1(t)=a1(y1-f(w1)x1),(9b)Dβy1(t)=z1-x1,(9c)Dγz1(t)=-b1y1+c1z1,(9d)Dw1(t)=x1,(10a)Dαx2(t)=a2(y2-f(w2)x2)+u1,(10b)Dβy2(t)=z2-x2+u2,(10c)Dγz2(t)=-b2y2+c2z2+u3,(10d)Dw2(t)=x2+u4,where the control function u=[u1;u2;u3;u4]T=0. Our goal is to determine the control function to make the first derivative of V(e), that is, V(e)<0. From (9a), (9b), (9c), and (9d) to (10a), (10b), (10c), and (10d), the error equations when the system parameters {a1,b1,c1}={a2,b2,c2}={a,b,c} are given as follows:(11a)Dαex=aey-a(f(w1)x1-f(w2)x2)+u1,(11b)Dβey=ez-ex+u2,(11c)Dγez=-bey+cez+u3,(11d)Dew=ex+u4,where ex=x2+x1, ey=y2+y1, ez=z2+z1, ew=w2+w1. In order to determine the controller, let
(12)u1=u1a+u1b,u1b=-a(f(w1)x1+f(w2)x2).
Then we rewrite system (11a), (11b), (11c), and (11d) in the following form:(13a)Dαex=aey+u1b,(13b)Dβey=ez-ex+u2,(13c)Dγez=-bey+cez+u3,(13d)Dew=ex+u4.Based on the Lyapunov stability theory, when the controller satisfies the assumption V(e)=0.5eTe, a positive definite function, the first derivative of V(e), V(e)<0, and the antisynchronization of two identical Chua’s circuit systems from different initial conditions is achieved, we choose u as follows:(14a)u1b=-ex-aey,(14b)u2=-ez+ex-ey,(14c)u3=-2cez+bey,(14d)u4=-ex-eu.Theorem 1.
The nonlinear controller achieves global antisynchronization between two identical Chua’s circuit systems master (9a), (9b), (9c), and (9d) and slave (10a), (10b), (10c), and (10d).Proof.
Take a Lyapunov function for (13a), (13b), (13c), and (13d) into consideration as
(15)V(e)=0.5eTe.
We then get the first derivative of V(e):(16a)V˙=ex(aey+u1b)+ey(ez-ex+u2)+ez(-bey+cez+u3)+ew(ex+u4),(16b)V˙=-ex2-ey2-bcez2-ew2,so V˙(e)<0 is satisfied. Since V˙(e) is a negative-definite function [39], the error states tend to zero at steady state. Therefore, the states of controlled response system and drive system are globally antisynchronized asymptotically.


## 3.2. Generalized Projective Synchronization of Two Identical Chua’s Circuits

In this subsection we study the synchronization between two identical fractional-order chaotic Chua’s circuits based on the generalized projective synchronization. The unknown termsu1, u2, u3, and u4 in (11a), (11b), (11c), and (11d) are active control functions to be determined. Define the error functions as
(17)(ex(t)ey(t)ez(t)ew(t))=(x1(t)+s1(t)x2(t)y1(t)+s2(t)y2(t)z1(t)+s3(t)z2(t)w1(t)+s4(t)w2(t)).
Then by similar procedure, the control functions are given by(18a)(u1(t)u2(t)u3(t)u4(t))=((a1f(w1(t))-kx)x1(t)-a1y1(t)s1(t)+(a2f(w2(t))-kx)x2(t)-a2y2(t)x1(t)-z1(t)-kyy1(t)s2(t)+x2(t)-z2(t)-kyy2(t)b1y1(t)-c1z1(t)-kzz1(t)s3(t)+b2y2(t)-c2z2(t)-kzz2(t)-x1(t)-kww1(t)s4(t)+-x2(t)-kww2(t)). Then
(18b)(DαexDβeyDγezDew)=(-kx0000-ky0000-kz0000-kw)(exeyezew)which are decaying functions as the values of kx, ky, kz, and kw are positive.


## 4. Numerical Results

Nonstandard finite difference method withϕ(h)=1-e-h is used to solve the systems of differential equations. In addition, a time step size 0.01 is employed. We select the parameters of the master system of Chua’s circuit system as a1=4, b1=1, c1=0.65, m1=0.2, and m2=5, so these systems exhibit chaotic behavior with the following initial conditions:
(19)(x1(0)y1(0)z1(0)w1(0))=(0.010.020.010.05),(x2(0)y2(0)z2(0)w2(0))=(0.020.030.020.06).
To validate the proposed technique, five different cases of the scaling parameters s1(t), s2(t), s3(t), and s4(t) are discussed as shown in Figure 3 and we have the following.The five different cases of the parameterss1, s2, s3, and s4 as function of time where Case A; {s1,s2,s3,s4}={1,1,1,1} Case B; {s1,s2,s3,s4}={-2,-2,-2,-2} Case C; {s1,s2,s3,s4}={2,3,4,5} Case D; s1=s2=s3=s4=-0.5+ integer part of (t/20), Case E; s1=s2=s3=s4=1+(fractionpartof(t/50))/50.
(a)
(b)
(c)
(d)(i) Case A (
{
s
1
,
s
2
,
s
3
,
s
4
}
=
{
1,1
,
1,1
}
). In this case the scaling parameters are positive ones; then the output responses of system 2 are the antisynchronization of the system 1 output, that is, x2=-x1 and the same for other outputs (see Figures 4(a) and 4(b)). Moreover, when an additive random noise is added to both the X1 and Y1 with amplitude 0.1 under the same conditions, the output response tries to follow the changes as shown in Figure 4(c).The general synchronization of (a), (b) thexandyresponses when {α,β,γ,a1,b1,c1,a2,b2,c2}={0.96,0.96,0.96,4,1,0.65,4,1,0.65} for Case A, (c) the x response under the same condition with additive noise with amplitude 0.1, (d) the response when {α,β,γ,a1,b1,c1,a2,b2,c2}={0.96,0.96,0.96,4,1,0.65,4,1,0.65} for Case B, and (e) Case C.
(a)
(b)
(c)
(d)
(e)(ii) Case B (
{
s
1
,
s
2
,
s
3
,
s
4
}
=
{
-
2
,
-
2
,
-
2
,
-
2
}
). In this case the scaling parameters are equal and negative; then the output responses of system 2 are the half-synchronization of the system 1 output, that is, x2=0.5x1 and the same for other outputs (see Figure 4(d)).(iii) Case C (
{
s
1
,
s
2
,
s
3
,
s
4
}
=
{
2,3
,
4,5
}
). In this case the scaling parameters are different constants and positive; then the output responses of system 2 have different antisynchronization scale than system 1 output, that is, x2=-0.5x1, y2=-y1/3, z2=-z1/4, and w2=-w1/5 (see Figure 4(e)).(iv) Case D (
s
1
=
s
2
=
s
3
=
s
4
=
-
0.5
+
integer part of
(
t
/
20
)
). In this case the scaling parameters are time dependent (staircase) as steps (every 20 seconds); then the output responses of system 2 will start with double synchronization in the first 20 seconds and then double antisynchronization in the following 20 seconds, and the scale factor increases with time every 20 seconds (see Figure 5).Theyresponse when for Case D when (a) {α,β,γ,a1,b1,c1,a2,b2,c2}={0.96,0.96,0.96,4,1,0.65,4,1,0.65}, (b) {α,β,γ,a1,b1,c1,a2,b2,c2}={0.9,0.95,1,4,1,0.65,4,1,0.65}, (c) {α,β,γ,a1,b1,c1,a2,b2,c2}={0.9,0.95,1,4,1,0.65,4.5,1.2,0.55}, and (d) {α,β,γ,a1,b1,c1,a2,b2,c2}={0.95,0.95,0.95,4,1,0.65,4,1,0.65}.
(a)
(b)
(c)
(d)(v) Case E (
s
1
=
s
2
=
s
3
=
s
4
=
1
+
(
f
r
a
c
t
i
o
n
p
a
r
t
o
f
(
t
/
50
)
)
/
50
). In this case the scaling parameters are time dependent as they ramp up to 50 seconds and then reset and increase again. Since it is always positive, then the output response will be the antisynchronization of system 1 with a different scale (see Figure 6).((a), (b)) thex1-y1 and x2-y2 projections, and ((c), (d)) the x, y, and z responses when {α,β,γ,a1,b1,c1,a2,b2,c2}={0.95,0.95,0.95,4,1,0.65,4,1,0.65} for Case E.
(a)
(b)
(c)
(d)Figure5 shows the proposed technique under different orders and different system parameters where the output will follow the expected response; however the transient time until the synchronization happens may be increased. Figures 5(b) and 5(c) illustrate the response with different fractional orders where α=0.9, β=0.95, and γ=1.0, and Figure 5(c) shows the responses for different fractional orders and also system parameters {α,β,γ,a1,b1,c1,a2,b2,c2}={0.9,0.95,1,4,1,0.65,4.5,1.2,0.55}. Moreover, Figure 5(d) shows the system responses for equal fractional orders and different system parameters only. Therefore, many cases have been discussed when the synchronization parameters are changed by steps showing great matching with the expected results. Figure 6 illustrates the projections and system responses when the synchronization parameters change gradually and linearly resulting in a sawtooth waveform (Case E). It is clear from the previous responses for different system parameters, fractional orders, and time-independent and time-dependent synchronization parameters that the proposed method has many advantages and can be used for the synchronization of any two chaotic systems. There are many recent research articles which discuss the importance of communication techniques using chaotic signals [40, 41]. In the amplitude modulation, the amplitude of the output signal should be a function of the input signal. In Figure 7, we assumed that the information data s1(t)=s2(t)=s3(t)=s4(t)=1.5+sin(0.02πt). Therefore, the output will be different from system 1 output as shown from the strange attractors shown in Figure 7(a). The time waveform of the scaling functions, z1(t), and z2(t), is shown in Figure 7(b) where the modulation is clear. The demodulation can be done similarly by reversing the operation.The projection and responses for the amplitude modulation whens1(t)=s2(t)=s3(t)=s4(t)=1.5+sin(0.02πt).
(a)
(b)


## 5. Conclusions

This paper discussed the concept of a general time-dependent synchronization scheme based on a nonlinear controller and then applied this technique on the fractional-order Chua’s circuit with memristor. This nonlinear controller is based on the control theory and Lyapunov stability theory to achieve the required synchronization. Many examples including antisynchronization, synchronization, and both (as in Case D) between two identical fractional-order Chua’s circuits with same/different fractional orders and same/different system parameters are discussed using numerical simulations by the nonstandard finite method. This technique can be repeated for other chaotic systems in a similar way to achieve good results.



---

*Source: 758676-2013-11-28.xml*