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

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.

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 Section 2 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 𝑓∶𝑅⁺→𝑅 is defined as follows: 𝐷^𝛼𝑓(𝑡)≡𝑑^𝛼𝑓(𝑡)𝑑𝑡^𝛼=1Γ(𝑚−𝛼)^𝑡₀𝑓^(𝑚)(𝜏)(𝑡−𝜏)^{𝛼−𝑚+1}𝑑𝑑𝜏,𝑚−1<𝛼<𝑚,_𝑚𝑑𝑡^𝑚𝑓(𝑡),𝛼=𝑚, (1) Where 𝑚 is the first integer greater than 𝛼 and Γ(⋅) is the gamma function and is defined by Γ(𝑧)=^∞₀𝑒^−𝑡𝑡^𝑧−1𝑑𝑡,Γ(𝑧+1)=𝑧Γ(𝑧). (2) The Grünwald-Letnikov method of approximation for the one-dimensional fractional derivative is as follows [33]:𝐷^𝛼𝑥(𝑡)(),𝐷=𝑓𝑡,𝑥_𝛼𝑥(𝑡)=lim_ℎ→0ℎ^{−𝛼𝑡/ℎ}_𝑗=0(−1)^𝑗𝛼𝑗𝑥𝑡−𝑗ℎ, (3a)(3b)where 𝛼>0 and 𝐷^∝ denotes the fractional derivative. Also 𝑁=[𝑡/ℎ], and ℎ is the step size. Therefore, (3a) and (3b) is discretized as follows: ^𝑛+1_𝑗=0𝑐^𝛼_𝑗𝑥𝑡−𝑗ℎ=𝑓𝑡_𝑛,𝑥𝑡_𝑛,𝑛=1,2,3,…, (4) Where 𝑡_𝑛=𝑛ℎ and 𝑐^𝛼_𝑗 are the Grünwald-Letnikov coefficients defined as 𝑐^𝛼_𝑗=1−1+𝛼𝑗𝑐^𝛼_𝑗−1,𝑐^𝛼₀=ℎ^−𝛼,𝑗=1,2,3,…. (5) The nonstandard discretization technique [34] is a general scheme where we replace the step size ℎ by a function 𝜑(ℎ). By applying this technique and using the Grünwald-Letnikov discretization method [2, 31], it yields the following relations: 𝑥_𝑛+1=−∑^𝑛+1_𝑗=1𝑐^𝛼_𝑗𝑥_{𝑛+1−𝑗}+𝑓₁𝑡_𝑛+1,𝑥_𝑛+1𝑐^𝛼1₀, (6) Where 𝑐^𝛼1₀=(𝜑₁(ℎ))⁻¹ are functions of the step size ℎ=Δ𝑡, with 𝜑₁(ℎ)=ℎ+𝑂(ℎ²), when ℎ→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 𝑞−𝜙 relationship. The canonical memristor-based fractional-order Chua’s circuit can be written as follows:𝐷^𝛼𝑥𝐷^𝛽𝑦𝐷^𝛾𝑧=01𝑎0001𝑐00000𝑥𝑦𝑧𝑢,𝑚𝐷𝑢−𝑎𝑘−1−𝑏𝑘=¹|𝑢|𝑚<1,₂|𝑢|>1, (7a)(7b)where 𝑎,𝑏,𝑐,𝑚₁, and 𝑚₂ are related to the circuit parameters [31]. The effect of the memristive property appears through the parameter 𝑘 (the slope of the charge with respect to the flux in the memristor) which is a piecewise-linear function with two different slopes 𝑚₁ and 𝑚₂ depending on the value of the flux 𝑢 (Figure 1(b)). It is clear that the system has four variables which are {𝑥,𝑦,𝑧, and 𝑢} 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)

(b)

(c)

(d)

Figure 1 

(a) The canonical Chua’s circuit with memristor, (b) the piecewise-linear flux-controlled memristor, and (c) 3D shapes in the 𝑢-𝑦-𝑧 plane when (𝛼,𝛽,𝛾) (a) (1.0, 0.9, 0.9) and (d) (1.1, 0.9, 0.9).

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 {𝑎,𝑏,𝑐}, 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 {𝑢₁,𝑢₂,𝑢₃,𝑢₄} on the four equations to minimize the functions 𝑥₁+𝑠₁𝑥₂, 𝑦₁+𝑠₂𝑦₂, 𝑧₁+𝑠₃𝑧₂, and 𝑤₁+𝑠₄𝑤₂ to approach zero. In this case, the output response of the second system is given by 𝑥₂(𝑡)𝑥=−¹(𝑡)𝑠₁(𝑡),𝑦₂(𝑡)𝑦=−¹(𝑡)𝑠₂(𝑡),𝑧₂(𝑡)𝑧=−¹(𝑡)𝑠₃(𝑡),𝑤₂(𝑡)𝑤=−¹(𝑡)𝑠₄(𝑡). (8) Then, if 𝑠₁=𝑠₂=𝑠₃=𝑠₄=−1, this will be the conventional synchronization, and if 𝑠₁=𝑠₂=𝑠₃=𝑠₄=1, this will be the antisynchronization as will be discussed in the next subsection. Moreover, the scaling functions 𝑠_𝑖(𝑡), 𝑖=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 (𝑠_𝑖=1, 𝑖=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)(9b)(9c)(9d)𝐷^𝛼𝑥₂(𝑡)=𝑎₂𝑦₂−𝑓𝑤₂𝑥₂+𝑢₁,𝐷^𝛽𝑦₂(𝑡)=𝑧₂−𝑥₂+𝑢₂,𝐷^𝛾𝑧₂(𝑡)=−𝑏₂𝑦₂+𝑐₂𝑧₂+𝑢₃,𝐷𝑤₂(𝑡)=𝑥₂+𝑢₄, (10a)(10b)(10c)(10d)where the control function 𝑢=[𝑢₁;𝑢₂;𝑢₃;𝑢₄]^𝑇=0. Our goal is to determine the control function to make the first derivative of 𝑉(𝑒), that is, 𝑉(𝑒)<0. From (9a), (9b), (9c), and (9d) to (10a), (10b), (10c), and (10d), the error equations when the system parameters {𝑎₁,𝑏₁,𝑐₁}={𝑎₂,𝑏₂,𝑐₂}={𝑎,𝑏,𝑐} are given as follows:𝐷^𝛼𝑒_𝑥=𝑎𝑒_𝑦−𝑎𝑓𝑤₁𝑥₁−𝑓𝑤₂𝑥₂+𝑢₁,𝐷^𝛽𝑒_𝑦=𝑒_𝑧−𝑒_𝑥+𝑢₂,𝐷^𝛾𝑒_𝑧=−𝑏𝑒_𝑦+𝑐𝑒_𝑧+𝑢₃,𝐷𝑒_𝑤=𝑒_𝑥+𝑢₄, (11a)(11b)(11c)(11d)where 𝑒_𝑥=𝑥₂+𝑥₁, 𝑒_𝑦=𝑦₂+𝑦₁, 𝑒_𝑧=𝑧₂+𝑧₁, 𝑒_𝑤=𝑤₂+𝑤₁. In order to determine the controller, let 𝑢₁=𝑢_1𝑎+𝑢_1𝑏,𝑢_1𝑏=−𝑎𝑓𝑤₁𝑥₁+𝑓𝑤₂𝑥₂. (12) Then we rewrite system (11a), (11b), (11c), and (11d) in the following form:𝐷^𝛼𝑒_𝑥=𝑎𝑒_𝑦+𝑢_1𝑏,𝐷^𝛽𝑒_𝑦=𝑒_𝑧−𝑒_𝑥+𝑢₂,𝐷^𝛾𝑒_𝑧=−𝑏𝑒_𝑦+𝑐𝑒_𝑧+𝑢₃,𝐷𝑒_𝑤=𝑒_𝑥+𝑢₄. (13a)(13b)(13c)(13d)Based on the Lyapunov stability theory, when the controller satisfies the assumption 𝑉(𝑒)=0.5𝑒^𝑇𝑒, a positive definite function, the first derivative of 𝑉(𝑒), 𝑉(𝑒)<0, and the antisynchronization of two identical Chua’s circuit systems from different initial conditions is achieved, we choose 𝑢 as follows:𝑢_1𝑏=−𝑒_𝑥−𝑎𝑒_𝑦,𝑢₂=−𝑒_𝑧+𝑒_𝑥−𝑒_𝑦,𝑢₃=−2𝑐𝑒_𝑧+𝑏𝑒_𝑦,𝑢₄=−𝑒_𝑥−𝑒_𝑢. (14a)(14b)(14c)(14d)

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 𝑉(𝑒)=0.5𝑒^𝑇𝑒. (15) We then get the first derivative of 𝑉(𝑒):𝑉̇=𝑒_𝑥𝑎𝑒_𝑦+𝑢_1𝑏+𝑒_𝑦𝑒_𝑧−𝑒_𝑥+𝑢₂+𝑒_𝑧−𝑏𝑒_𝑦+𝑐𝑒_𝑧+𝑢₃+𝑒_𝑤𝑒_𝑥+𝑢₄,𝑉̇=−𝑒²_𝑥−𝑒²_𝑦−𝑏𝑐𝑒²_𝑧−𝑒²_𝑤, (16a)(16b)so 𝑉̇(𝑒)<0 is satisfied. Since 𝑉̇(𝑒) 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 terms 𝑢₁, 𝑢₂, 𝑢₃, and 𝑢₄ in (11a), (11b), (11c), and (11d) are active control functions to be determined. Define the error functions as 𝑒_𝑥(𝑡)𝑒_𝑦(𝑡)𝑒_𝑧(𝑡)𝑒_𝑤(𝑡)=𝑥₁(𝑡)+𝑠₁(𝑡)𝑥₂(𝑡)𝑦₁(𝑡)+𝑠₂(𝑡)𝑦₂(𝑡)𝑧₁(𝑡)+𝑠₃(𝑡)𝑧₂(𝑡)𝑤₁(𝑡)+𝑠₄(𝑡)𝑤₂(𝑡). (17) Then by similar procedure, the control functions are given by𝑢₁(𝑡)𝑢₂(𝑡)𝑢₃(𝑡)𝑢₄(𝑡)=𝑎₁𝑓𝑤₁(𝑡)−𝑘_𝑥𝑥₁(𝑡)−𝑎₁𝑦₁(𝑡)𝑠₁(𝑡)+𝑎₂𝑓𝑤₂(𝑡)−𝑘_𝑥𝑥₂(𝑡)−𝑎₂𝑦₂(𝑡)𝑥₁(𝑡)−𝑧₁(𝑡)−𝑘_𝑦𝑦₁(𝑡)𝑠₂(𝑡)+𝑥₂(𝑡)−𝑧₂(𝑡)−𝑘_𝑦𝑦₂(𝑡)𝑏₁𝑦₁(𝑡)−𝑐₁𝑧₁(𝑡)−𝑘_𝑧𝑧₁(𝑡)𝑠₃(𝑡)+𝑏₂𝑦₂(𝑡)−𝑐₂𝑧₂(𝑡)−𝑘_𝑧𝑧₂(𝑡)−𝑥₁(𝑡)−𝑘_𝑤𝑤₁(𝑡)𝑠₄(𝑡)+−𝑥₂(𝑡)−𝑘_𝑤𝑤₂(𝑡). (18a) Then 𝐷^𝛼𝑒_𝑥𝐷^𝛽𝑒_𝑦𝐷^𝛾𝑒_𝑧𝐷𝑒_𝑤=−𝑘_𝑥0000−𝑘_𝑦0000−𝑘_𝑧0000−𝑘_𝑤𝑒_𝑥𝑒_𝑦𝑒_𝑧𝑒_𝑤 (18b)which are decaying functions as the values of 𝑘_𝑥, 𝑘_𝑦, 𝑘_𝑧, and 𝑘_𝑤 are positive.

4. Numerical Results

Nonstandard finite difference method with 𝜙(ℎ)=1−𝑒^−ℎ 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 𝑎₁=4, 𝑏₁=1, 𝑐₁=0.65, 𝑚₁=0.2, and 𝑚₂=5, so these systems exhibit chaotic behavior with the following initial conditions: 𝑥₁(0)𝑦₁(0)𝑧₁(0)𝑤₁(0)=,𝑥0.010.020.010.05₂(0)𝑦₂(0)𝑧₂(0)𝑤₂(0)=.0.020.030.020.06 (19) To validate the proposed technique, five different cases of the scaling parameters 𝑠₁(𝑡), 𝑠₂(𝑡), 𝑠₃(𝑡), and 𝑠₄(𝑡) are discussed as shown in Figure 3 and we have the following.

(a)

(b)

(c)

(d)

Figure 3 

The five different cases of the parameters 𝑠1, 𝑠2, 𝑠3, and 𝑠4 as function of time where Case A; {𝑠1,𝑠2,𝑠3,𝑠4}={1,1,1,1} Case B; {𝑠1,𝑠2,𝑠3,𝑠4}={−2,−2,−2,−2} Case C; {𝑠1,𝑠2,𝑠3,𝑠4}={2,3,4,5} Case D; 𝑠1=𝑠2=𝑠3=𝑠4=−0.5+ integer part of (𝑡/20), Case E; 𝑠1=𝑠2=𝑠3=𝑠4=1+(fractionpartof(𝑡/50))/50.

(i) Case A ({s₁,s₂,s₃,s₄}={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, 𝑥₂=−𝑥₁ and the same for other outputs (see Figures 4(a) and 4(b)). Moreover, when an additive random noise is added to both the 𝑋₁ and 𝑌₁ with amplitude 0.1 under the same conditions, the output response tries to follow the changes as shown in Figure 4(c).

(a)

(b)

(c)

(d)

(e)

Figure 4 

The general synchronization of (a), (b) the𝑥and𝑦responses when {𝛼,𝛽,𝛾,𝑎1,𝑏1,𝑐1,𝑎2,𝑏2,𝑐2}={0.96,0.96,0.96,4,1,0.65,4,1,0.65} for Case A, (c) the 𝑥 response under the same condition with additive noise with amplitude 0.1, (d) the response when {𝛼,𝛽,𝛾,𝑎1,𝑏1,𝑐1,𝑎2,𝑏2,𝑐2}={0.96,0.96,0.96,4,1,0.65,4,1,0.65} for Case B, and (e) Case C.

(ii) Case B ({𝑠₁,𝑠₂,𝑠₃,𝑠₄}={−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, 𝑥₂=0.5𝑥₁ and the same for other outputs (see Figure 4(d)).

(iii) Case C ({𝑠₁,𝑠₂,𝑠₃,𝑠₄}={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, 𝑥₂=−0.5𝑥₁, 𝑦₂=−𝑦₁/3, 𝑧₂=−𝑧₁/4, and 𝑤₂=−𝑤₁/5 (see Figure 4(e)).

(iv) Case D (𝑠₁=𝑠₂=𝑠₃=𝑠₄=−0.5+integer part of  (𝑡/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).

(a)

(b)

(c)

(d)

Figure 5 

The𝑦response when for Case D when (a) {𝛼,𝛽,𝛾,𝑎1,𝑏1,𝑐1,𝑎2,𝑏2,𝑐2}={0.96,0.96,0.96,4,1,0.65,4,1,0.65}, (b) {𝛼,𝛽,𝛾,𝑎1,𝑏1,𝑐1,𝑎2,𝑏2,𝑐2}={0.9,0.95,1,4,1,0.65,4,1,0.65}, (c) {𝛼,𝛽,𝛾,𝑎1,𝑏1,𝑐1,𝑎2,𝑏2,𝑐2}={0.9,0.95,1,4,1,0.65,4.5,1.2,0.55}, and (d) {𝛼,𝛽,𝛾,𝑎1,𝑏1,𝑐1,𝑎2,𝑏2,𝑐2}={0.95,0.95,0.95,4,1,0.65,4,1,0.65}.

(v) Case E (𝑠₁=𝑠₂=𝑠₃=𝑠₄=1+(𝑓𝑟𝑎𝑐𝑡𝑖𝑜𝑛𝑝𝑎𝑟𝑡𝑜𝑓(𝑡/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)

(c)

(d)

Figure 6 

((a), (b)) the 𝑥1-𝑦1 and 𝑥2-𝑦2 projections, and ((c), (d)) the 𝑥, 𝑦, and 𝑧 responses when {𝛼,𝛽,𝛾,𝑎1,𝑏1,𝑐1,𝑎2,𝑏2,𝑐2}={0.95,0.95,0.95,4,1,0.65,4,1,0.65} for Case E.

Figure 5 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 {𝛼,𝛽,𝛾,𝑎₁,𝑏₁,𝑐₁,𝑎₂,𝑏₂,𝑐₂}={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 𝑠₁(𝑡)=𝑠₂(𝑡)=𝑠₃(𝑡)=𝑠₄(𝑡)=1.5+sin(0.02𝜋𝑡). 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, 𝑧₁(𝑡), and 𝑧₂(𝑡), is shown in Figure 7(b) where the modulation is clear. The demodulation can be done similarly by reversing the operation.

(a)

(b)

Figure 7 

The projection and responses for the amplitude modulation when 𝑠1(𝑡)=𝑠2(𝑡)=𝑠3(𝑡)=𝑠4(𝑡)=1.5+sin(0.02𝜋𝑡).

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.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.