# A Solution for the Generalized Synchronization of a Class of Chaotic Systems Based on Output Feedback

Authors: Carlos Aguilar-Ibanez; R. Martínez-Guerra; C. Pérez-Pinacho; E. García-Canseco; Miguel S. Suarez-Castanon
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

A solution to the output-feedbackgeneralized synchronization problem for two chaotic systems, namely, the master and the slave, is presented. The solution assumes that the slave is controlled by a single input, and the states of each system are partially known. To this end, both systems are expressed in their corresponding observable generalized canonical form, through their differential primitive element. The nonavailable state variables of both systems are recovered using a suitable Luenberger observer. The convergence analysis was carried out using the linear control approach in conjunction with the Lyapunov method. Convincing numerical simulations are presented to assess the effectiveness of the obtained solution.


---

## Body

## 1. Introduction

The synchronization of strictly different chaotic systems is a current challenge in Control Theory. In general, accomplishing and understanding this kind of synchronization, referred to as generalized synchronization (GS), are a difficult problem and are very important because its solution can be useful to actual applications. For instance, chaotic signals have been used in secure information transmission. In this regard, interesting secure communication methods based on GS can be found in the literature. In [1], a channel-independent chaotic secure communication scheme based on GS is proposed. A similar scheme, also based on GS, with a remarkable stability to noise and numerically validated using the Rössler system is developed in [2]. Two novel generalized chaos synchronization based secure communication schemes are introduced in [3, 4]. The first scheme is complemented with a transposition function and is tested using the Chen chaotic circuit; the second scheme is oriented to encrypt images and uses a theorem introduced by the authors, which is a generalization of GS to an array of differential equations. From a theoretical point of view, synchronization of strictly different chaotic systems helps to explain how some oscillatory systems, like the intrinsic neuron model, act cooperatively to develop tasks or solve problems [5]. This behavior also occurs in nature and can be responsible for the transition to low-dimensional behavior in systems with many degrees of freedom. A full review of GS is beyond the scope of this study; however, we mention two seminal works. In [6] the authors developed a method for decomposing chaotic systems such that linear GS can be achieved. Finally, we mention work [5], where the authors proposed a method used to detect and study generalized synchronization in drive-response systems. The method uses an identical response system to monitor the synchronized motions. We suggest to the interested reader the following list of remarkable works [7–16].Loosely speaking the GS phenomenon occurs when the trajectories of one system, through a functional mapping, are equal to trajectories of another. In other words, let us suppose that we have two systems,Σm and Σs,u, named, respectively, master and slave, with their corresponding trajectories, xm(t) and xs(t) [17]. That is, GS is accomplished [18] if there exists a differential primitive element that generates a mapping Hms from the trajectories xm(t) of Σm in the algebraic manifold M to the trajectories xs(t) of Σs,u in its space S, in which case it holds that Hms(xs(t))=xm(t). It is noteworthy that for identical systems the functional mapping Hms corresponds to the identity [19]. In our case, we assume that the behavior of the slave system Σs,u can be either chaotic or bounded and oscillatory. Evidently, to accomplish GS there must be a control action acting over the slave Σs,u (for simplicity, symbols Σm and Σs,u refer to the uncontrolled master system and the controlled slave system, resp. In the same way, xm and xs, are, resp., the states of systems Σm and Σs,u.)In this work, we solve the output-feedback GS problem for systemsΣm and Σs,u, when only a single output of each system is available or measurable, and only a single controller is used in the slave system. This restriction makes the GS problem more challenging and interesting than the works mentioned above. To solve it, the nonavailable state variables of systems Σm and Σs,u must be reconstructed from the measurable outputs. To overcome this issue, we exploit somead hoc differential algebraic properties found in some chaotic systems, by choosing the differentiable primitive element [20]. This primitive element allows us to express both systems in a generalized canonical form. Thence, we solve the GS problem as if it were an instance of the trajectory-tracking problem. To this end, it is necessary to split the slave controller into two parts: one devoted to feedback the output tracking error between the master and the slave systems and the other compensating for the underlying and nonavailable state variables of each system. We underscore that our control design basically consists of finding a suitable differentiable primitive element to express the master and slave systems in their corresponding canonical form, in conjunction with a kind of Luenberger observer to reconstruct the nonavailable variable states. Fortunately, there are in nature several chaotic systems that can be expressed in this form.The remainder of this paper is organized as follows. In Section2, we introduce some concepts necessary to establish the problem statement properly. In Section 3, we solve the output-feedback GS control problem. Numerical experiments to assess the effectiveness of our solution are presented in Section 4. The final remarks are given in Section 5.Notation. To simplify the development of this work, we adopt the following notation:(1)Ai=0100⋯00010⋯0⋮⋮⋮⋮⋮⋮00⋯⋯010000⋯0∈Ri×i;Bi=0⋮01∈Ri;Ci=1⋮00∈Ri.


## 2. Preliminaries

Before proceeding, we introduce some needful concepts [20–23]. From the theorem of differential primitive element [20] we posit a single element δ∈L, which is a differential primitive element, such that L=K(δ); that is, L is differentially generated by K and δ (K and L are differential fields).Consider the following nonlinear system:(2)x˙=Fx,u,where x∈Rn is the system state and u is the system input. Also, suppose that there exists an auxiliary variable y¯ (differential primitive element), defined as (3)y¯=∑iαixi,αi∈R,which allows us to write system (2) in the following form:(4)ξ˙=Anξ+Bn-£ξ+pu,where ξ=[ξ1,ξ2,…,ξn]T with ξi=y-(i-1), for i={1,…,n}, £(·) is a scalar nonlinear function of their arguments and may not be defined everywhere, and p is the gain constant of the single input of u. The above system representation is well known as the generalized controller canonical form, and variable y¯ is the measurable output of the system [21, 24].


### 2.1. Problem Statement

In this section we introduce formally the main problem of this work, which consists of solving the output-feedback GS problem for two chaotic systems, in a master-slave configuration, where the states of both systems are partially known. To this end, we propose a suitable definition for the GS problem, based on previous works [21]. It should be noted that our definition uses the notion of differential primitive element, already mentioned.Let us consider the following two nonlinear systems in a master-slave configuration, where the master system,Σm, is given by (5)x˙m=Fmxm;ym=hmxm,and the controlled slave system, Σs,u, is defined by(6)x˙s=Fsxs,uxs,xm;ys=hsxs,where xs=x1s,…,xns∈Rn, xm=x1m,…,xnm∈Rn, hs:Rn→R, hm:Rn→R, u:Rn×Rn→R, ym, and ys∈R. Fs, Fm, hs, and hm are assumed to be polynomial in their arguments. Before introducing the GS definition, we need to introduce the following two important assumptions, related to some structural properties of Σm and Σs,u and based on the result found in [21].Assumption A1. Suppose that the master system, Σm, given in (5), can be written as(7)z˙m=Anzm+Bn-£mzm,where zm=[zm1,…,zmn]T=Φm(xm) is the vector of function, with zmk=ym(k-1) for k={1,…,n}, ym is a differential primitive element for the master system, and £m(·) is a scalar nonlinear function. Additionally, the trajectories of free system (7) must be forward complete in some compact set in Rn.Assumption A2. Suppose that the slave system Σs,u, given in (6), can be written as(8)z˙s=Anzs+Bn-£szs+pu,where zs=[zs1,…,zsn]T=Φs(xs) with zsk=ys(k-1), for k={1,…,n}, ys is a differential primitive element for the slave system, p is a gain constant of the single input u, and £s(·) is a scalar nonlinear function.Notice that the functions vectors,Φm,Φs:Rn→Rn, are generated through the corresponding differentiable primitive element. This construction allows us to solve straightforwardly the GS problem. Finally, we mention that constructions (7) and (8) are in agreement with the canonical form previously introduced in [23, 25]. Now we can introduce the main definition of this study, as follows.Definition 1 (generalized synchronization).
Let us consider systems,Σm and Σs,u, under Assumptions A1 and A2, with Φs being invertible. The slave and master systems are said to be in a state of GS if there exists a differential primitive element that generates a transformation Hms:Rn→Rn with Hms=Φs-1∘Φm and there exist an algebraic manifold M=xs,xm∣xm=Hms(xs) and a compact set B⊂Rn×Rn with M⊂B such that their trajectories with initial conditions in B approach M as t→∞. This definition leads to the following criterion: (9)limt→∞⁡Hmsxs-xm=0.Remark 2.
As far as we know, the GS problem has been solved assuming that statesxs,xm are available for measurement. In this work, we moved a step forward, by considering the case when the states of both systems are partially known. To this end, the proposed controller is based on output-feedback stabilization, instead of using full-state control. It is important, because in several actual applications the whole state is not available or is unpractical or impossible to measure. For instance, in the classical Duffing system only the position is available. On the other hand, finding a general solution of this problem is rather difficult.As it can be seen original systems (5) and (6) are expressed in their corresponding observable generalized canonical form, through their corresponding differential primitive element. It is worthy to mention that, fortunately, several chaotic systems can be expressed in such observable generalized canonical form. Now, if Assumptions A1 and A2 hold, there exists a full-state feedback controller, u(zm,zs), such that limt→∞zm-zs=0 [21]. We underscore that this result assumes that the whole state of both systems, the master and the slave, is available for measurements.Now, we can introduce the main problem of this work.Problem Statement. Consider the uncontrolled master chaotic system (7) and the controlled slave chaotic system (8), satisfying, respectively, Assumptions A1 and A2, and the vector function Φs-1 existing at least locally. The main goal consists of proposing the slave system controller:(10)uys,y^s1,…,y^sn,ym,…,y^mn,which solves the GS problem, where y^(i) are estimations of the unknown y(i).As this shows, we are dealing with a more difficult problem configuration. Contrary to several previous works [5, 6, 26], we do not use any control action over the master chaotic system, while the slave controller only uses one measurable output from the slave and one from the master. This problem instance can be considered more appealing, because it resembles several actual implementations. On the other hand, it is less expensive to control a system using state estimations than using actual measurements.To solve this problem, it is necessary to reconstruct the nonavailable dynamics of both systems, the master and the slave. Therefore, we introduce a suitable version of the assumption introduced in [21].Assumption P1. Suppose that there exists a primitive element, ym, of the master system “Σm” (5) that transforms it into the new form “ΣmT”, given by(11)z˙m=Anzm+Bnϕmym+KmTzm+fmzm,where zm=Φs(xm), ϕm is a scalar nonlinear function that depends on the output ym, Km is a vector of constants, and fm(z) is a scalar Lipschitz function in z on the open set D⊂Rn. That is, (12)fmx-fmy≤γmx-y;∀x,y∈D.Assumption P2. Suppose that there exists a primitive element, ys, of the slave system “Σs,u” (6) that transforms it into the new form “ΣsT,u”, given by(13)z˙s=Anzs+Bnϕsys+KsTzs+fszs+pu,where zs=Φs(xs), ϕs is a scalar nonlinear function that depends on the output ys, Ks is a constant vector, and fs(zs) is a scalar Lipschitz function in z, with its corresponding constant γs. Finally, p≠0 is a fixed constant and u is the single input of the system.Comment 1. Assumptions A1 and A2 are needed to solve the GS problem, when the states of both systems, the master and the slave, are known. Assumptions P1 and P2 are needed to solve the GS problem when only a single output of each system can be measurable. As we can see, our goal is to move a step forward with respect to the previous works. On the other hand, there are many chaotic systems that admit the representation, given by (11) and (13).


## 2.1. Problem Statement

In this section we introduce formally the main problem of this work, which consists of solving the output-feedback GS problem for two chaotic systems, in a master-slave configuration, where the states of both systems are partially known. To this end, we propose a suitable definition for the GS problem, based on previous works [21]. It should be noted that our definition uses the notion of differential primitive element, already mentioned.Let us consider the following two nonlinear systems in a master-slave configuration, where the master system,Σm, is given by (5)x˙m=Fmxm;ym=hmxm,and the controlled slave system, Σs,u, is defined by(6)x˙s=Fsxs,uxs,xm;ys=hsxs,where xs=x1s,…,xns∈Rn, xm=x1m,…,xnm∈Rn, hs:Rn→R, hm:Rn→R, u:Rn×Rn→R, ym, and ys∈R. Fs, Fm, hs, and hm are assumed to be polynomial in their arguments. Before introducing the GS definition, we need to introduce the following two important assumptions, related to some structural properties of Σm and Σs,u and based on the result found in [21].Assumption A1. Suppose that the master system, Σm, given in (5), can be written as(7)z˙m=Anzm+Bn-£mzm,where zm=[zm1,…,zmn]T=Φm(xm) is the vector of function, with zmk=ym(k-1) for k={1,…,n}, ym is a differential primitive element for the master system, and £m(·) is a scalar nonlinear function. Additionally, the trajectories of free system (7) must be forward complete in some compact set in Rn.Assumption A2. Suppose that the slave system Σs,u, given in (6), can be written as(8)z˙s=Anzs+Bn-£szs+pu,where zs=[zs1,…,zsn]T=Φs(xs) with zsk=ys(k-1), for k={1,…,n}, ys is a differential primitive element for the slave system, p is a gain constant of the single input u, and £s(·) is a scalar nonlinear function.Notice that the functions vectors,Φm,Φs:Rn→Rn, are generated through the corresponding differentiable primitive element. This construction allows us to solve straightforwardly the GS problem. Finally, we mention that constructions (7) and (8) are in agreement with the canonical form previously introduced in [23, 25]. Now we can introduce the main definition of this study, as follows.Definition 1 (generalized synchronization).
Let us consider systems,Σm and Σs,u, under Assumptions A1 and A2, with Φs being invertible. The slave and master systems are said to be in a state of GS if there exists a differential primitive element that generates a transformation Hms:Rn→Rn with Hms=Φs-1∘Φm and there exist an algebraic manifold M=xs,xm∣xm=Hms(xs) and a compact set B⊂Rn×Rn with M⊂B such that their trajectories with initial conditions in B approach M as t→∞. This definition leads to the following criterion: (9)limt→∞⁡Hmsxs-xm=0.Remark 2.
As far as we know, the GS problem has been solved assuming that statesxs,xm are available for measurement. In this work, we moved a step forward, by considering the case when the states of both systems are partially known. To this end, the proposed controller is based on output-feedback stabilization, instead of using full-state control. It is important, because in several actual applications the whole state is not available or is unpractical or impossible to measure. For instance, in the classical Duffing system only the position is available. On the other hand, finding a general solution of this problem is rather difficult.As it can be seen original systems (5) and (6) are expressed in their corresponding observable generalized canonical form, through their corresponding differential primitive element. It is worthy to mention that, fortunately, several chaotic systems can be expressed in such observable generalized canonical form. Now, if Assumptions A1 and A2 hold, there exists a full-state feedback controller, u(zm,zs), such that limt→∞zm-zs=0 [21]. We underscore that this result assumes that the whole state of both systems, the master and the slave, is available for measurements.Now, we can introduce the main problem of this work.Problem Statement. Consider the uncontrolled master chaotic system (7) and the controlled slave chaotic system (8), satisfying, respectively, Assumptions A1 and A2, and the vector function Φs-1 existing at least locally. The main goal consists of proposing the slave system controller:(10)uys,y^s1,…,y^sn,ym,…,y^mn,which solves the GS problem, where y^(i) are estimations of the unknown y(i).As this shows, we are dealing with a more difficult problem configuration. Contrary to several previous works [5, 6, 26], we do not use any control action over the master chaotic system, while the slave controller only uses one measurable output from the slave and one from the master. This problem instance can be considered more appealing, because it resembles several actual implementations. On the other hand, it is less expensive to control a system using state estimations than using actual measurements.To solve this problem, it is necessary to reconstruct the nonavailable dynamics of both systems, the master and the slave. Therefore, we introduce a suitable version of the assumption introduced in [21].Assumption P1. Suppose that there exists a primitive element, ym, of the master system “Σm” (5) that transforms it into the new form “ΣmT”, given by(11)z˙m=Anzm+Bnϕmym+KmTzm+fmzm,where zm=Φs(xm), ϕm is a scalar nonlinear function that depends on the output ym, Km is a vector of constants, and fm(z) is a scalar Lipschitz function in z on the open set D⊂Rn. That is, (12)fmx-fmy≤γmx-y;∀x,y∈D.Assumption P2. Suppose that there exists a primitive element, ys, of the slave system “Σs,u” (6) that transforms it into the new form “ΣsT,u”, given by(13)z˙s=Anzs+Bnϕsys+KsTzs+fszs+pu,where zs=Φs(xs), ϕs is a scalar nonlinear function that depends on the output ys, Ks is a constant vector, and fs(zs) is a scalar Lipschitz function in z, with its corresponding constant γs. Finally, p≠0 is a fixed constant and u is the single input of the system.Comment 1. Assumptions A1 and A2 are needed to solve the GS problem, when the states of both systems, the master and the slave, are known. Assumptions P1 and P2 are needed to solve the GS problem when only a single output of each system can be measurable. As we can see, our goal is to move a step forward with respect to the previous works. On the other hand, there are many chaotic systems that admit the representation, given by (11) and (13).


## 3. Solving the GS Control Problem for a Class of Chaotic Nonlinear Systems

In this section we present a solution to the main problem of this study. To this end, we design a Lyapunov-based observer to recover the underlying dynamics of the master and slave systems,Σm and Σs,u. Afterwards, using the recovered dynamics by a static feedback u, we solve the GS problem.Reconstructing the Nonavailable Dynamic of
Σ
m
or
Σ
s
,
u
. Suppose that we have a nonlinear system, defined by(14)z˙=Anz+Bnϕy+fz+KTz+pu,y=CnTz,where y is the single measurable output, z∈Rn is the partially known state, ϕ(y) is any scalar function, and f(z) is a scalar Lipschitz function in z∈Rn, with γ as the Lipschitz constant, and K∈Rn constant. Now, consider the observer of this nonlinear system, defined as(15)z^˙=Anz^+Bnϕym+fz^+KTz^+pu-LnκCnTz^-y,where z^ is an estimation of z, Lnκ=[κ1,κ2,…,κn]T∈Rn is a vector of constants, and u is the system input. Notice that system (14) has a similar structure to that of the systems defined in Assumptions P1 and P2. The following proposition provides sufficient conditions to ensure that the estimation error z~=z-z^ converges to zero, as long as t→∞.Proposition 3.
Consider systems (14) and (15). If the vector of constants, Lnκ, is selected, according to (16)PnκMnκ+MnκTPnκ=-In,where Mnκ=(An+BnKT-LnκCnT) and Pnκ>0 satisfying(17)2γr¯n<1,where r¯n=PnκBn (r¯n is in fact the last row of matrix Pnκ), then, the error z~=z-z^ converges asymptotically and exponentially to zero.Proof.
From (11) and (15), it is easy to show that the dynamic of the error z~ is given by(18)z~˙=An+BnKT-LnκCnTz~+Bnfz-fz^.In order to analyze the convergence of the state z~, we propose the Lyapunov function V(z~)=z~TPnκz~, where Pnκ=PnκT>0, whose time derivative, along of (18), is given by(19)V˙z~=-z~2+2z~PnκBnfz-fz^.Now, it is quite easy to see that the following inequality holds: (20)z~PnκBnfz-fz^≤γPnκBnz~2.Therefore, from (19) and (20), we have that V˙ can be upper bounded, as(21)V˙z~≤-z~21-2γr¯n.Evidently, V˙<0, for all z~≠0, in the case that inequality (17) is fulfilled.The following remark allows us to propose the needed observers forΣmT and ΣsT,u.Remark 4.
Assumptions P1 and P2, in conjunction with Proposition3, suggest proposing the observer Σ^mT for the master system ΣmT (11) as(22)z^˙m=Anz^m+Bnϕmym+KmTz^m+fmz^m-Lnκmym-CnTz^m,where ym=zm1. Similarly, observer Σ^sT,u for the slave system can be proposed as(23)z^˙s=Anz^s+Bnϕsys+KsTz^s+fszs+pus-Lnκsys-CnTzswith ys=zs1, where Lnκm and Lnκs must be selected, such that the following matrices are obtained:(24)Mnκm=An+BnKmT-LnκmCnT;Mnκs=An+BnKsT-LnκsCnT,being Hurwitz and fulfilling the Lyapunov equation PnκiMnκi+MnκiTPnκi=-In, i={s,m}, with the following restrictions:(25)2γsPnκsBn<1;2γmPnκmBn<1.In the next section, we use the Lyapunov-based observer to solve the GS control problem.


### 3.1. Generalized Synchronization

Now we are ready to introduce the main proposition of this paper.Proposition 5.
Assuming that the observers, (22) and (23), satisfy the conditions in Remark 4, the following static feedback,(26)uz^s,z^m=1p-ϕ¯ym,ys-K¯z^m,z^s-f¯z^m,z^s-Kn,zz^m-z^s,where,(27)ϕ¯ym,ys≙ϕmym-ϕsys;K¯x,y≙KmTx-KsTy;f¯x,y≙fmx-fsy,ensures that the synchronization error ez=zs-zm exponentially converges to zero.Proof.
From Assumptions P1 and P2, we have that systemsΣm and Σs,u can be expressed according to (11) and (13), respectively. Computing the dynamic error, ez, we have(28)e˙z=Anez+Bnϕ¯ym,ys+K¯zm,zs+f¯zm,zs+pu.Therefore, introducing u(z^s,z^m), proposed in (26), into (28), we obtain(29)e˙z=Anez+Bn-K¯z^m,z^s+K¯zm,zs+Δz¯-Kn,zz^m-z^s,where z¯=(zmT,zsT,z^mT,z^sT) and Δ are defined by(30)Δz¯=f¯zm,zs-f¯z^m,z^s.Remembering that es=zs-z^s and em=zm-z^m, we have the following relations:(31)K¯z^m,z^s-K¯zm,zs=KmTz^m-KsTz^s-KmTzm+KsTzs;=-KmTem+KsTes.Similarly, (32)Kn,zz^m-z^s=Kn,zez+Kn,zes-em.Substituting these two relations into (29), we obtain(33)e˙z=An-BnKn,zez+BnKn,zes-em+KsTes-KmTem+BnΔz¯︷we,z¯.Using simple algebra, we have that Δ(z¯) can be upper bounded, as follows:(34)Δ=f¯zm,zs-f¯z^m,z^s=fmzm-fmz^m-fszs-fsz^s≤γmem+γses.Defining(35)A¯ns=An-BnKn,z,provided that A¯ns is Hurwitz, (33) can be written in the following simple way:(36)e˙z=A¯nsez+Bnwe,z¯,where w(e,z¯) is given by(37)we,z¯=BnKn,zes-em+KsTes-KmTem+BnΔz¯.However, according to (34), w(∗) can be upper bounded, as follows:(38)w∗≤Kn,zes-em+KsTes+KmTem+γmem+γses.By assumption, es and em converge exponentially and asymptotically to zero. Then, there are some α and β being strictly positives, such that w(∗)≤αe-βt. Consequently, the error ez of system (36) asymptotically and exponentially converges to zero.


## 3.1. Generalized Synchronization

Now we are ready to introduce the main proposition of this paper.Proposition 5.
Assuming that the observers, (22) and (23), satisfy the conditions in Remark 4, the following static feedback,(26)uz^s,z^m=1p-ϕ¯ym,ys-K¯z^m,z^s-f¯z^m,z^s-Kn,zz^m-z^s,where,(27)ϕ¯ym,ys≙ϕmym-ϕsys;K¯x,y≙KmTx-KsTy;f¯x,y≙fmx-fsy,ensures that the synchronization error ez=zs-zm exponentially converges to zero.Proof.
From Assumptions P1 and P2, we have that systemsΣm and Σs,u can be expressed according to (11) and (13), respectively. Computing the dynamic error, ez, we have(28)e˙z=Anez+Bnϕ¯ym,ys+K¯zm,zs+f¯zm,zs+pu.Therefore, introducing u(z^s,z^m), proposed in (26), into (28), we obtain(29)e˙z=Anez+Bn-K¯z^m,z^s+K¯zm,zs+Δz¯-Kn,zz^m-z^s,where z¯=(zmT,zsT,z^mT,z^sT) and Δ are defined by(30)Δz¯=f¯zm,zs-f¯z^m,z^s.Remembering that es=zs-z^s and em=zm-z^m, we have the following relations:(31)K¯z^m,z^s-K¯zm,zs=KmTz^m-KsTz^s-KmTzm+KsTzs;=-KmTem+KsTes.Similarly, (32)Kn,zz^m-z^s=Kn,zez+Kn,zes-em.Substituting these two relations into (29), we obtain(33)e˙z=An-BnKn,zez+BnKn,zes-em+KsTes-KmTem+BnΔz¯︷we,z¯.Using simple algebra, we have that Δ(z¯) can be upper bounded, as follows:(34)Δ=f¯zm,zs-f¯z^m,z^s=fmzm-fmz^m-fszs-fsz^s≤γmem+γses.Defining(35)A¯ns=An-BnKn,z,provided that A¯ns is Hurwitz, (33) can be written in the following simple way:(36)e˙z=A¯nsez+Bnwe,z¯,where w(e,z¯) is given by(37)we,z¯=BnKn,zes-em+KsTes-KmTem+BnΔz¯.However, according to (34), w(∗) can be upper bounded, as follows:(38)w∗≤Kn,zes-em+KsTes+KmTem+γmem+γses.By assumption, es and em converge exponentially and asymptotically to zero. Then, there are some α and β being strictly positives, such that w(∗)≤αe-βt. Consequently, the error ez of system (36) asymptotically and exponentially converges to zero.


## 4. Numerical Simulations

The goal of this example consists of solving the GS problem for the controlled slave Chua system and the free master Colpitts oscillator, provided that a single output of each system is available for measurements. On the other hand, we underscore that, in the present case, the trajectories of the Colpitts system are bounded in some compact set and exhibit chaotic or oscillatory behavior.Master System. Consider the state equations for the Colpitts oscillator Σm, given by(39)x˙1=-ame-x2+amx3+am,x˙2=bmx3,x˙3=-cmx1-cmx2-dmx3,where am=bmC2/C1, bm=I0/w0C2Vt, cm=Vt/w0LIo, and dm=R/Lw0 (to avoid confusion, the master parameters use the subscript m, and the slave parameters use the subscript s). It is well known that this oscillator exhibits, for certain values of its parameters, a chaotic behavior. We will show that, selecting ym=x2 as the differential primitive element, we can write it according to expression (11). We use xmT=(x1,x2,x3) to refer to the vector state of the Colpitts system. From ym=x2, we obtain the following representation:(40)zm=Φmxm=ymy˙my¨m=x2bmx3bm-cmx1-cmx2-dmx3,where Φm is invertible. Furthermore, the time derivative of zm can be written as(41)z˙m1z˙m2z˙m3=zm2zm3Ψnxm,where(42)Ψmxm=-bmcmx˙1-bmcmx˙2-bmdmx˙3.Notice that the term Ψm(xm) can be expressed in the coordinates of zm, as(43)Ψmzm=ambmcm-1+e-ym︷ϕmym+-amcm-bmcmz2m-dmz3m︷KmTzm.Notice that (41) and (43) are in agreement with the condition in Assumption P1. According to Proposition 3 and Remark 4, the observer can be proposed as(44)z^˙m=A3z^m+B3ϕmym+KmTz^m-L3κmC3Tz^m-zm,where L3κm has to be selected, such that M3κm=(A3+B3KmT-L3κmC3T) is Hurwitz. Using the following actual parameters reported in [27],(45)a=b=6.2723,c=0.0797,d=0.6898,it is easy to show that selecting L3κs=[6,11,6] allows us to obtain the following matrix:(46)M3κm=-610-1101-6-0.249-0.689,which is Hurwitz. Therefore, proposed observer (44) ensures that the estimation error em=zm-z^m converges exponentially and asymptotically to zero.For the numerical simulation, the experimental setup assumed the above defined parameter values, the master initial conditions asxm(0)=(1,-1,2), and the master observer as x^m(0)=0. The obtained results are shown in Figure 1, where we can see that the estimation error of each state variable settles to almost zero after nearly 4 seconds.Figure 1
Observation error for the Colpitts or master system in the transformed coordinates:em=zm-z^m.Slave System. Consider the Chua chaotic system as the controlled slave system, Σs,u. This system is formed by three linear energy-storage elements (an inductor and two capacitors), a linear resistor, and a single nonlinear resistor. It is well known that this circuit is described by the following differential equations:(47)x˙1=asx2-x1-vsx1+u,x˙2=x1-x2+x3,x˙3=-bsx2,with(48)vsx=m1x+m0-m12x+1-x-1,where parameters as, bs, m0, and m1 are chosen, such that system (47) exhibits a chaotic behavior when u=0. Notice that u is the input of the slave system that is necessary to achieve the GS problem. Particularly, in the case where u=0 and the parameters values are in the neighbourhood of [28],(49)bs=28;as=15.6;m0=-87;m1=-57,it is known that we have the so-called double scroll chaotic attractors. As before, we use symbol xsT=(x1,x2,x3) to indicate the states of slave system (47). In a similar fashion, we choose as the primitive element ys=x3. Therefore, to obtain the corresponding canonical form, we derive ys two times, obtaining(50)zs=Φsxs=ysy˙sy¨s=x3-bx2-bsx1-x2+x3,where Φs is an invertible map. From the above equation, the time derivative of zs is given by(51)z˙s1z˙s2z˙s3=zs2zs3Ψsxs-bsu,where(52)Ψsxs=-bsx˙1-x˙2+x˙3.Notice that Ψs(xs) can be rewritten in the coordinates of zs as(53)Ψsys,zs=KsTzs+asbsfsx1-bsu,where fs(x1)=vx(x1) and(54)x1=-zs1-zs2+zs3b;KsT=-ab-b-1-a.On the other hand, it is clear that vx is Lipchitz, with γs=m1+m0-m1 (see Appendix). Notice that (51) and (53) are in agreement with Assumption P2. Now, according to Proposition 3 and Remark 4, the observer can be proposed as(55)z^˙s=A3z^s+B3abvxx1-bu+KsTz^s-L3κsC3Tz^s-zs.Using L3κs=[24,26,9] and the list of data (49), we have the following matrix:(56)M3κs=-2410-2601-399-25-16,which is stable, and γs=1.1429. Solving the Lyapunov equation, (57)P3κsM3κs+M3κsTP3κs=-I3,it is easy to see that P3κsB3=0.308, ensuring that 2γs(0.308)<1. That is, the corresponding observer proposed in (55) ensures that the estimation error es=zs-z^s converges exponentially and asymptotically to zero.Similarly to the master system, the slave experimental setup uses the above defined parameter values, and the initial conditions of the slave and its observer were fixed asxs(0)=(1,0.5,0.4) and x^s(0)=0, respectively. The slave input was fixed as u=0. In Figure 2, it is easy to see that the errors are very close to zero after 3.5-second elapse.Figure 2
Observation error for the Chua or slave system in the transformed coordinates:es=zs-z^s.Colpitts Oscillator and Chua Circuit Generalized Synchronization. For this numerical simulation, we used the same experimental setup and the same observer as in the previous experiments. The setup for the slave controller was fixed according to Proposition 5. The used static feedback was(58)uz^s,z^m=1p-ϕ¯ym-K¯z^m,z^s-f¯zs-Kn,zz^m-z^s,where(59)ϕ¯ym=ambmcm-1+e-ym;K¯z^m,z^s=KmTz^m-KsTz^s;f¯zs=-asbsvs-zs1-zs2+zs3bs.The synchronization error of this experiment is shown in Figure3, where we can see that, after 9-second elapse, the difference between each pair of transformed coordinates is almost zero, for the case r=1.Figure 3
Generalized synchronization error between the free Colpitts oscillator and the driven Chua circuit in transformed coordinates.


## 5. Conclusions

In this paper, we proposed a solution to the output-feedback GS problem for a class of chaotic systems, namely, the free master and the controlled slave, assuming that the master system trajectories are bounded in a compact set. The problem was solved given that the slave was controlled using a single input, and a single output of each system was available or measurable. The solution consisted of expressing both systems in their corresponding observable generalized canonical form, through the differentiable primitive element. To this end, a suitable diffeomorphism transformation was applied, between the original state space and the corresponding canonical observability form of each system. Due to the fact that the controlled slave system uses its own whole state and the whole state of the master, a kind of a Luenberger observer for each system was proposed. Finally, the information obtained through the observers is used by the controlled slave system, solving the GS problem. The corresponding convergence analysis was carried out using the linear control approach and the Lyapunov method. The effectiveness of the obtained method was assessed by convincing numerical simulations, using the chaotic Colpitts oscillator as the master and the chaotic Chua circuit as the slave.



---

*Source: 848203-2015-10-26.xml*