# Generalized Synchronization of Nonlinear Chaotic Systems through Natural Bioinspired Controlling Strategy

Authors: Shih-Yu Li; Shi-An Chen; Chin-Teng Lin; Li-Wei Ko; Cheng-Hsiung Yang; Heng-Hui Chen
Journal: Abstract and Applied Analysis
 (2015)

---

## Abstract

A novel bioinspired control strategy design is proposed for generalized synchronization of nonlinear chaotic systems, combining the bioinspired stability theory, fuzzy modeling, and a novel, simple-form Lyapunov control function design of derived high efficient, heuristic and bioinspired controllers. Three main contributions are concluded: (1) apply the bioinspired stability theory to further analyze the stability of fuzzy error systems; the high performance of controllers has been shown in previous study by Li and Ge 2009, (2) a new Lyapunov control function based on bioinspired stability theory is designed to achieve synchronization without using traditional LMI method, which is a simple linear homogeneous function of states and the process of designing controller to synchronize two fuzzy chaotic systems becomes much simpler, and (3) three different situations of synchronization are proposed; classical master and slave Lorenz systems, slave Chen’s system, and Rossler’s system as functional system are illustrated to further show the effectiveness and feasibility of our novel strategy. The simulation results show that our novel control strategy can be applied to different and complicated control situations with high effectiveness.


---

## Body

## 1. Introduction

Since Pecora and Carroll first obtained the synchronization of two identical chaotic systems under different initial conditions [1], synchronization of chaotic systems has received significant attention. Nowadays, more and more techniques of chaos synchronization were proposed, such as active control [2–4], backstepping control method [5–8], linear error feedback control [9–11], adaptive control [12–16], and sliding mode control [17–20].In our real world, various kinds of nonlinear chaotic biosystems describing natural phenomenon are found to have some states always positive. It means these states are always in the first quadrant, such as the three species prey-predator systems [15, 16, 21], double Mackey-Glass systems [17, 22, 23], energy communication system in biological research [19, 24], and virus-immune system [20]. For the three species prey-predator systems, which consist of two competing preys and one predator can be described by the following set of nonlinear differential equations:(1)dxdT=r11-k1-1x-k1-1c2y-Φ1x,yz,dydT=r21-k2-1c1x-k2-1y-Φ2x,yz,dzdT=e1Φ1x,yz+e2Φ2x,yz-αz,where α, ri, ki, ei, and ci, i=1,2 are the model parameters assuming only positive values, and the functions Φi(x,y), i=1,2 represent the densities of the two prey species and z represents the density of the predator species. The predator z consumes the preys x,y according to the response functions [25]. The prescribed model characterized by nonlinear response since amount of food consumed by predator per unit time depends upon the available food sources from the two preys x and y. And for the virus-immune system, a mathematical model is considered as follows:(2)dTdt=s-μ1T+rT1-T+IΓ-kVT,dIdt=kVT-μ2I,dVdt=μ3NI-kVT-μ4V,where T, I, and V represent the population concentrations of uninfected, infected target cells and virus, respectively. We denote by the s constant supply of target cells from its precursor, where s is the rate at which new T cells are created. These cells have finite life time and μ1 represent the average death rates of these cells, μ2 and μ3 represent the death rate of the infected cells before viral production commences, and μ4 is the viral clearance rate constant. These target cells are assumed to grow logistically with specific growth rate r and carrying capacity Γ. Above all, their states are all positive.As a result, a natural bioinspired stability theory is proposed for generalized synchronization in this paper. It means that there exists a given functional relationship between the states of the master and that of the slave. Via using the natural bioinspired stability theory, the new Lyapunov function is a simple linear homogeneous function of states and introduces less simulation error.In this paper, a new Lyapunov control function based on the bioinspired stability theory is proposed to synchronize fuzzy chaotic systems. This new Lyapunov control function is designed as a simple linear homogeneous function of states and the process of designing controller to synchronize two fuzzy chaotic systems becomes much simpler. Furthermore, we design three different synchronization situations to show the effectiveness of our strategy—(1)CASE I, synchronization of identical systems: synchronization of master fuzzy Lorenz system and slave fuzzy Lorenz system, (2)CASE II, synchronization of two different systems: synchronization of master fuzzy Lorenz system and slave fuzzy Chen system, and (3)CASE III, synchronization of complicated different systems: synchronization of master fuzzy Lorenz system and complicated fuzzy chaotic system (slave fuzzy Chen system and functional fuzzy Rossler’s system).The layout of the rest of the paper is as follows. In Section2, the natural bioinspired stability theory is introduced. In Section 3, modeling of chaotic systems is proposed. In Section 4, three simulation examples are given. In Section 5, conclusions are given.


## 2. Natural Bioinspired Controlling Strategy

### 2.1. Definition of the Stability on Partial Region

Consider the differential equations of disturbed motion of a nonautonomous system in the normal form(3)dxsdt=Xst,x1,…,xn,s=1,…,n,where the function Xs is defined on the intersection of the partial region Ω (shown in Figure 1) and(4)∑sxs2≤Hand t>t0, where t0 and H are certain positive constants. Xs, which vanishes when the variables xs are all zero, is a real valued function of t, x1,…,xn. It is assumed that Xs is smooth enough to ensure the existence and uniqueness of the solution of the initial value problem. When Xs does not contain t explicitly, the system is autonomous.Figure 1
Bioinspired partial regionsΩ and Ω1.Obviously,xs=0 (s=1,…,n) is a solution of (3). We are interested in the asymptotical stability of this zero solution on partial region Ω (including the boundary) of the neighborhood of the origin which in general may consist of several subregions (Figure 1).Definition 1.
For any given numberε>0, if there exists a δ>0, such that on the closed given partial region Ω when(5)∑sxs02≤δ,s=1,…,nfor all t≥t0, the inequality(6)∑sxs2<ε,s=1,…,nis satisfied for the solutions of (3) on Ω, then the disturbed motion xs=0 (s=1,…,n) is stable on the partial region Ω.Definition 2.
If the undisturbed motion is stable on the partial regionΩ, and there exists a δ′>0, so that on the given partial region Ω when(7)∑sxs02≤δ′,s=1,…,n,the equality(8)limt→∞⁡∑sxs2=0is satisfied for the solutions of (3) on Ω, then the undisturbed motion xs=0 (s=1,…,n) is asymptotically stable on the partial region Ω.The intersection ofΩ and region defined by (7) is called the region of attraction.Definition of Functions
V
(
t
,
x
1
,
…
,
x
n
). Let us consider the functions V(t,x1,…,xn) given on the intersection Ω1 of the partial region Ω and the region(9)∑sxs2≤h,s=1,…,nfor t≥t0>0, where t0 and h are positive constants. We suppose that the functions are single-valued and have continuous partial derivatives and become zero when x1=⋯=xn=0.Definition 3.
If there existst0>0 and a sufficiently small h>0, so that on partial region Ω1 and t≥t0, V≥0 (or ≤0), then V is a positive (or negative) semidefinite, in general semidefinite, function on the Ω1 and t≥t0.Definition 4.
If there exists a positive (negative) definitive functionW(x1⋯xn) on Ω1, so that on the partial region Ω1 and t≥t0(10)V-W≥0or-V-W≥0,then V(t,x1,…,xn) is a positive definite function on the partial region Ω1 and t≥t0.Definition 5.
IfV(t,x1,…,xn) is neither definite nor semidefinite on Ω1 and t≥t0, then V(t,x1,…,xn) is an indefinite function on partial region Ω1 and t≥t0. That is, for any small h>0 and any large t0>0, V(t,x1,…,xn) can take either positive or negative value on the partial region Ω1 and t≥t0.Definition 6 (bounded functionV).
If there existst0>0, h>0, so that on the partial region Ω1, we have(11)Vt,x1,…,xn<L,where L is a positive constant, then V is said to be bounded on Ω1.Definition 7 (function with infinitesimal upper bound).
IfV is bounded, and for any λ>0, there exists μ>0, so that on Ω1 when ∑sxs2≤μ, and t≥t0, we have(12)Vt,x1,…,xn≤λthen V admits an infinitesimal upper bound on Ω1.


### 2.2. Bioinspired GYC Theorem of Stability and of Asymptotical Stability on Partial Region

Theorem 8.
If there can be found a definite functionV(t,x1,…,xn) on the partial region for (3) and the derivative with respect to time based on these equations is(13)dVdt=∂V∂t+∑s=1n∂V∂xsXs,then, it is a semidefinite function on the partial region whose sense is opposite to that of V, or if it becomes zero identically, then the undisturbed motion is stable on the partial region.Proof.
Let us assume for the sake of definiteness thatV is a positive definite function. Consequently, there exists a sufficiently large number t0 and a sufficiently small number h<H, such that on the intersection Ω1 of partial region Ω and(14)∑sxs2≤h,s=1,…,nand t≥t0, the following inequality is satisfied(15)Vt,x1,…,xn≥Wx1,…,xn,where W is a certain positive definite function which does not depend on t. Besides that, (13) may assume only negative or zero value in this region.
Letε be an arbitrarily small positive number. We will suppose that in any case ε<h. Let us consider the aggregation of all possible values of the quantities x1,…,xn, which are on the intersection ω2 of Ω1 and(16)∑sxs2=ε,and let us designate by l>0 the precise lower limit of the function W under this condition. By virtue of (10), we will have(17)Vt,x1,…,xn≥lforx1,…,xnon⁡ω2.
We will now consider the quantitiesxs as functions of time which satisfy the differential equations of disturbed motion. We will assume that the initial values xs0 of these functions for t=t0 lie on the intersection Ω2 of Ω1 and the region(18)∑sxs2≤δ,where δ is so small that(19)Vt0,x10,…,xn0<l.
By virtue of the fact thatV(t0,0,…,0)=0, such a selection of the number δ is obviously possible. We will suppose that in any case the number δ is smaller than ε. Then the inequality(20)∑sxs2<εbeing satisfied at the initial instant will be satisfied, in the very least, for a sufficiently small t-t0, since the functions xs(t) very continuously with time. We will show that these inequalities will be satisfied for all values t>t0. Indeed, if these inequalities were not satisfied at some time, there would have to exist such an instant t=T for which this inequality would become an equality. In other words, we would have(21)∑sxs2T=ε,and consequently, on the basis of (19)(22)VT,x1T,…,xnT≥l.
On the other hand, sinceε<h, the inequality (9) is satisfied in the entire interval of time [t0,T], and consequently, in this entire time interval dV/dt≤0. This yields(23)VT,x1T,…,xnT≤Vt0,x10,…,xn0,which contradicts (20) on the basis of (19). Thus, the inequality (23) must be satisfied for all values of t>t0.
Finally, we must point out that from the view-point of mathematics, the stability on partial region in general does not be related logically to the stability on whole region. If an undisturbed solution is stable on a partial region, it may be either stable or unstable on the whole region and vice versa. In specific practical problems, we do not study the solution starting withinΩ2 and running out of Ω.Theorem 9.
If in satisfying the conditions of Theorem8 the derivative dV/dt is a definite function on the partial region with opposite sign to that of V and the function V itself permits an infinitesimal upper limit, then the undisturbed motion is asymptotically stable on the partial region.Proof.
Let us suppose thatV is a positive definite function on the partial region and that, consequently, dV/dt is negative definite. Thus on the intersection Ω1 of Ω and the region defined by (9) and t≥t0 there will be satisfied not only the inequality (10), but also the following inequality: (24)dVdt≤-W1x1,…,xn,where W1 is a positive definite function on the partial region independent of t.
Let us consider the quantitiesxs as functions of time which satisfy the differential equations of disturbed motion assuming that the initial values xs0=xs(t0) of these quantities satisfy the inequalities (18). Since the undisturbed motion is stable in any case, the magnitude δ may be selected so small that for all values of t≥t0 the quantities xs remain within Ω1. Then, on the basis of (24) the derivative of function V(t,x1(t),…,xn(t)) will be negative at all times and, consequently, this function will approach a certain limit, as t increases without limit, remaining larger than this limit at all times. We will show that this limit is equal to some positive quantity different from zero. Then for all values of t≥t0 the following inequality will be satisfied:(25)Vt,x1t,…,xnt>α,where α>0.
SinceV permits an infinitesimal upper limit, it follows from this inequality that(26)∑sxs2t≥λ,s=1,…,n,where λ is a certain sufficiently small positive number. Indeed, if such a number λ did not exist, that is, if the quantity ∑sxs(t) was smaller than any preassigned number no matter how small, then the magnitude V(t,x1(t),…,xn(t)), as follows from the definition of an infinitesimal upper limit, would also be arbitrarily small, which contradicts (25).
If for all values oft≥t0 the inequality (26) is satisfied, then (24) shows that the following inequality will be satisfied at all times:(27)dVdt≤-l1,where l1 is positive number different from zero which constitutes the precise lower limit of the function W1(t,x1(t),…,xn(t)) under condition (26). Consequently, for all values of t≥t0 we will have(28)Vt,x1t,…,xnt=Vt0,x10,…,xn0+∫t0tdVdtdt≤Vt0,x10,…,xn0-l1t-t0,which is, obviously, in contradiction with (25). The contradiction thus obtained shows that the function V(t,x1(t),…,xn(t)) approached zero as t increase without limit. Consequently, the same will be true for the function W(x1(t),…,xn(t)) as well, from which it follows directly that(29)limt→∞⁡xst=0,s=1,…,n,which proves the theorem.Discussion and Statement. Theorems 8 and 9 in Bioinspired GYC Theorem of Stability of Asymptotical Stability on Partial Region provides a new way to investigate the stability problem—to solve the stability problem in the first quadrant, calledPartial Region, which is inspired via the biological behavior in nature and is used to control the states of a system which exist in the first quadrant to achieve the goal system with a set of simple controllers through designing a simpler and more convenient Lyapunov function.


## 2.1. Definition of the Stability on Partial Region

Consider the differential equations of disturbed motion of a nonautonomous system in the normal form(3)dxsdt=Xst,x1,…,xn,s=1,…,n,where the function Xs is defined on the intersection of the partial region Ω (shown in Figure 1) and(4)∑sxs2≤Hand t>t0, where t0 and H are certain positive constants. Xs, which vanishes when the variables xs are all zero, is a real valued function of t, x1,…,xn. It is assumed that Xs is smooth enough to ensure the existence and uniqueness of the solution of the initial value problem. When Xs does not contain t explicitly, the system is autonomous.Figure 1
Bioinspired partial regionsΩ and Ω1.Obviously,xs=0 (s=1,…,n) is a solution of (3). We are interested in the asymptotical stability of this zero solution on partial region Ω (including the boundary) of the neighborhood of the origin which in general may consist of several subregions (Figure 1).Definition 1.
For any given numberε>0, if there exists a δ>0, such that on the closed given partial region Ω when(5)∑sxs02≤δ,s=1,…,nfor all t≥t0, the inequality(6)∑sxs2<ε,s=1,…,nis satisfied for the solutions of (3) on Ω, then the disturbed motion xs=0 (s=1,…,n) is stable on the partial region Ω.Definition 2.
If the undisturbed motion is stable on the partial regionΩ, and there exists a δ′>0, so that on the given partial region Ω when(7)∑sxs02≤δ′,s=1,…,n,the equality(8)limt→∞⁡∑sxs2=0is satisfied for the solutions of (3) on Ω, then the undisturbed motion xs=0 (s=1,…,n) is asymptotically stable on the partial region Ω.The intersection ofΩ and region defined by (7) is called the region of attraction.Definition of Functions
V
(
t
,
x
1
,
…
,
x
n
). Let us consider the functions V(t,x1,…,xn) given on the intersection Ω1 of the partial region Ω and the region(9)∑sxs2≤h,s=1,…,nfor t≥t0>0, where t0 and h are positive constants. We suppose that the functions are single-valued and have continuous partial derivatives and become zero when x1=⋯=xn=0.Definition 3.
If there existst0>0 and a sufficiently small h>0, so that on partial region Ω1 and t≥t0, V≥0 (or ≤0), then V is a positive (or negative) semidefinite, in general semidefinite, function on the Ω1 and t≥t0.Definition 4.
If there exists a positive (negative) definitive functionW(x1⋯xn) on Ω1, so that on the partial region Ω1 and t≥t0(10)V-W≥0or-V-W≥0,then V(t,x1,…,xn) is a positive definite function on the partial region Ω1 and t≥t0.Definition 5.
IfV(t,x1,…,xn) is neither definite nor semidefinite on Ω1 and t≥t0, then V(t,x1,…,xn) is an indefinite function on partial region Ω1 and t≥t0. That is, for any small h>0 and any large t0>0, V(t,x1,…,xn) can take either positive or negative value on the partial region Ω1 and t≥t0.Definition 6 (bounded functionV).
If there existst0>0, h>0, so that on the partial region Ω1, we have(11)Vt,x1,…,xn<L,where L is a positive constant, then V is said to be bounded on Ω1.Definition 7 (function with infinitesimal upper bound).
IfV is bounded, and for any λ>0, there exists μ>0, so that on Ω1 when ∑sxs2≤μ, and t≥t0, we have(12)Vt,x1,…,xn≤λthen V admits an infinitesimal upper bound on Ω1.


## 2.2. Bioinspired GYC Theorem of Stability and of Asymptotical Stability on Partial Region

Theorem 8.
If there can be found a definite functionV(t,x1,…,xn) on the partial region for (3) and the derivative with respect to time based on these equations is(13)dVdt=∂V∂t+∑s=1n∂V∂xsXs,then, it is a semidefinite function on the partial region whose sense is opposite to that of V, or if it becomes zero identically, then the undisturbed motion is stable on the partial region.Proof.
Let us assume for the sake of definiteness thatV is a positive definite function. Consequently, there exists a sufficiently large number t0 and a sufficiently small number h<H, such that on the intersection Ω1 of partial region Ω and(14)∑sxs2≤h,s=1,…,nand t≥t0, the following inequality is satisfied(15)Vt,x1,…,xn≥Wx1,…,xn,where W is a certain positive definite function which does not depend on t. Besides that, (13) may assume only negative or zero value in this region.
Letε be an arbitrarily small positive number. We will suppose that in any case ε<h. Let us consider the aggregation of all possible values of the quantities x1,…,xn, which are on the intersection ω2 of Ω1 and(16)∑sxs2=ε,and let us designate by l>0 the precise lower limit of the function W under this condition. By virtue of (10), we will have(17)Vt,x1,…,xn≥lforx1,…,xnon⁡ω2.
We will now consider the quantitiesxs as functions of time which satisfy the differential equations of disturbed motion. We will assume that the initial values xs0 of these functions for t=t0 lie on the intersection Ω2 of Ω1 and the region(18)∑sxs2≤δ,where δ is so small that(19)Vt0,x10,…,xn0<l.
By virtue of the fact thatV(t0,0,…,0)=0, such a selection of the number δ is obviously possible. We will suppose that in any case the number δ is smaller than ε. Then the inequality(20)∑sxs2<εbeing satisfied at the initial instant will be satisfied, in the very least, for a sufficiently small t-t0, since the functions xs(t) very continuously with time. We will show that these inequalities will be satisfied for all values t>t0. Indeed, if these inequalities were not satisfied at some time, there would have to exist such an instant t=T for which this inequality would become an equality. In other words, we would have(21)∑sxs2T=ε,and consequently, on the basis of (19)(22)VT,x1T,…,xnT≥l.
On the other hand, sinceε<h, the inequality (9) is satisfied in the entire interval of time [t0,T], and consequently, in this entire time interval dV/dt≤0. This yields(23)VT,x1T,…,xnT≤Vt0,x10,…,xn0,which contradicts (20) on the basis of (19). Thus, the inequality (23) must be satisfied for all values of t>t0.
Finally, we must point out that from the view-point of mathematics, the stability on partial region in general does not be related logically to the stability on whole region. If an undisturbed solution is stable on a partial region, it may be either stable or unstable on the whole region and vice versa. In specific practical problems, we do not study the solution starting withinΩ2 and running out of Ω.Theorem 9.
If in satisfying the conditions of Theorem8 the derivative dV/dt is a definite function on the partial region with opposite sign to that of V and the function V itself permits an infinitesimal upper limit, then the undisturbed motion is asymptotically stable on the partial region.Proof.
Let us suppose thatV is a positive definite function on the partial region and that, consequently, dV/dt is negative definite. Thus on the intersection Ω1 of Ω and the region defined by (9) and t≥t0 there will be satisfied not only the inequality (10), but also the following inequality: (24)dVdt≤-W1x1,…,xn,where W1 is a positive definite function on the partial region independent of t.
Let us consider the quantitiesxs as functions of time which satisfy the differential equations of disturbed motion assuming that the initial values xs0=xs(t0) of these quantities satisfy the inequalities (18). Since the undisturbed motion is stable in any case, the magnitude δ may be selected so small that for all values of t≥t0 the quantities xs remain within Ω1. Then, on the basis of (24) the derivative of function V(t,x1(t),…,xn(t)) will be negative at all times and, consequently, this function will approach a certain limit, as t increases without limit, remaining larger than this limit at all times. We will show that this limit is equal to some positive quantity different from zero. Then for all values of t≥t0 the following inequality will be satisfied:(25)Vt,x1t,…,xnt>α,where α>0.
SinceV permits an infinitesimal upper limit, it follows from this inequality that(26)∑sxs2t≥λ,s=1,…,n,where λ is a certain sufficiently small positive number. Indeed, if such a number λ did not exist, that is, if the quantity ∑sxs(t) was smaller than any preassigned number no matter how small, then the magnitude V(t,x1(t),…,xn(t)), as follows from the definition of an infinitesimal upper limit, would also be arbitrarily small, which contradicts (25).
If for all values oft≥t0 the inequality (26) is satisfied, then (24) shows that the following inequality will be satisfied at all times:(27)dVdt≤-l1,where l1 is positive number different from zero which constitutes the precise lower limit of the function W1(t,x1(t),…,xn(t)) under condition (26). Consequently, for all values of t≥t0 we will have(28)Vt,x1t,…,xnt=Vt0,x10,…,xn0+∫t0tdVdtdt≤Vt0,x10,…,xn0-l1t-t0,which is, obviously, in contradiction with (25). The contradiction thus obtained shows that the function V(t,x1(t),…,xn(t)) approached zero as t increase without limit. Consequently, the same will be true for the function W(x1(t),…,xn(t)) as well, from which it follows directly that(29)limt→∞⁡xst=0,s=1,…,n,which proves the theorem.Discussion and Statement. Theorems 8 and 9 in Bioinspired GYC Theorem of Stability of Asymptotical Stability on Partial Region provides a new way to investigate the stability problem—to solve the stability problem in the first quadrant, calledPartial Region, which is inspired via the biological behavior in nature and is used to control the states of a system which exist in the first quadrant to achieve the goal system with a set of simple controllers through designing a simpler and more convenient Lyapunov function.


## 3. Modeling of Chaotic Systems

In this section, the well-famous Takagi-Sugeno fuzzy model is applied to model the classical Lorenz system, Chen system, and Rossler’s system for further designing the natural bioinspired controller.


### 3.1. Fuzzy Modeling of Master Lorenz System

For master Lorenz system(30)x˙1=ax2-x1,x˙2=cx1-x1x3-x2,x˙3=x1x2-bx3,where a, b, c are the parameters. When a=10, b=8/3, c=28, and initial states are (0.5,1,5), the dynamic behavior is chaotic. Assume that x1∈[-d1,d1] and d1>0; then Lorenz system can be exactly represented by T-S fuzzy model as follows:(31)Rule1:IFxisM1,THENX˙t=A1Xt,(32)Rule2:IFxisM2,THENX˙t=A2Xt,where(33)X=x1,x2,x3T,A1=-aa0c-1-d10d1-b,A2=-aa0c-1d10-d1-b,M1x=121+x1d1,M2x=121-x1d1,and d1=30. M1 and M2 are fuzzy set of Lorenz system. Here, we call (31) the first liner subsystem under the fuzzy rule and (32) the second liner subsystem under the fuzzy rule. The final output of the fuzzy Lorenz system is inferred as follows and the chaotic behavior is shown in Figure 2(34)X˙t=∑i=12HiAiXt,where(35)H1=M1M1+M2,H2=M2M1+M2.We further summarize all the parameters in Section 3.1 into Table 1.Table 1
System
Parameters
Initial states
Fuzzy Lorenz system (master)
a
=
10, b=8/3, c=28, d1=30
(x1, x2, x3) = (0.5, 1, 5)Figure 2
Chaotic behavior of master fuzzy Lorenz system.


### 3.2. Fuzzy Modeling of Slave Lorenz System

For slave Lorenz system(36)y˙1=ay2-y1+u1,y˙2=cy1-y1y3-y2+u2,y˙3=y1y2-by3+u3,where a, b, c are the parameters. When a=10, b=8/3, c=28, and initial states are (0.6,3,10), the dynamic behavior is chaotic. Assume that y1∈[-d2,d2] and d2>0; then Lorenz system can be exactly represented by T-S fuzzy model as follows:(37)Rule1:IFyisM~1,THENY˙t=A~1Yt+U1,(38)Rule2:IFyisM~2,THENY˙t=A~2Yt+U2,where(39)Y=y1,y2,y3T,A~1=-aa0c-1-d20d2-b,A~2=-aa0c-1d20-d2-b,M~1x=121+y1d2,M~2x=121-y1d2,U1=u11,u12,u13T,U2=u21,u22,u23Tand d2=20. M~1 and M~2 are fuzzy set of Lorenz system. Here, we call (37) the first linear subsystem under the fuzzy rule and (38) the second linear subsystem under the fuzzy rule. As a result, U1 and U2 are controllers of the first and second subsystems. The final output of the fuzzy Lorenz system is inferred as follows and the chaotic behavior without using controllers is shown in Figure 3(40)Y˙t=∑i=12H~iA~iYt+Ui,where(41)H~1=M~1M~1+M~2,H~2=M~2M~1+M~2.We further summarize all the parameters in Section 3.2 into Table 2.Table 2
System
Parameters
Initial states
Fuzzy Lorenz system (slave)
a
=
10, b=8/3, c=28, d2=20
(y1, y2, y3) = (0.6, 3, 10)Figure 3
Chaotic behavior of slave fuzzy Lorenz system.


### 3.3. Fuzzy Modeling of Chen System as Slave System

For Chen system as slave system(42)z˙1=ez2-z1+u1,z˙2=g-ez1-z1z3+gz2+u2,z˙3=z1z2-fz3+u3,where e, f, g are the parameters. When e=35, f=3, g=28, and initial states are (0.5,1,5), the dynamic behavior is chaotic. Assume that z1∈[-d3,d3] and d3>0; then Chen system can be exactly represented by T-S fuzzy model as follows:(43)Rule1:IFzisN1,THENZ˙t=B1Zt,(44)Rule2:IFzisN2,THENZ˙t=B2Zt.Here(45)Z=z1,z2,z3T,B1=-ee0g-eg-d30d3-f,B2=-ee0g-egd30-d3-f,N1z=121+z1d3,N2z=121-z1d3,U1=u11,u12,u13T,U2=u21,u22,u23Tand d3=30. N1 and N2 are fuzzy set of Lorenz system. Here, we call (43) the first liner subsystem under the fuzzy rule and (44) the second liner subsystem under the fuzzy rule. As a result, U1 and U2 are controllers of the first and second subsystems. The final output of the fuzzy Chen system is inferred as follows and the chaotic behavior without using controllers is shown in Figure 4(46)Z˙t=∑i=12IiBiZt+Ui,where(47)I1=N1N1+N2,I2=N2N1+N2.We further summarize all the parameters in Section 3.3 into Table 3.Table 3
System
Parameters
Initial states
Fuzzy Chen system (slave)
e
=
35, f=3, g=28, d3=30
(z1, z2, z3) = (0.5, 1, 5)Figure 4
Chaotic behavior of slave fuzzy Chen system.


### 3.4. Fuzzy Modeling of Rossler’s System as Functional System

For Rossler’s system as functional system(48)w˙1=-w2-w3,w˙2=w1+hw2,w˙3=i+w1w3-jw3,where e, f, g are the parameters. When h=0.2, i=0.2, j=5.7, and initial states are (10,15,11), the dynamic behavior is chaotic. Assume that w1∈[-d4,d4] and d4>0; then Chen system can be exactly represented by T-S fuzzy model as follows:(49)Rule1:IFwisO1,THENW˙t=C1Wt+c1,(50)Rule2:IFwisO2,THENW˙t=C2Wt+c2.Here(51)W=w1,w2,w3T,C1=0-1-11h00d4-j,C2=0-1-11h00-d4-j,c1=00i,c2=00iO1z=121+w1d4,O2z=121-w1d4,and d4=10. O1 and O2 are fuzzy set of Lorenz system. Here, we call (49) the first liner subsystem under the fuzzy rule and (50) the second liner subsystem under the fuzzy rule. The final output of the fuzzy Rossler’s system is inferred as follows and the chaotic behavior is shown in Figure 5(52)Z˙t=∑i=12JiCiZt+ci,where(53)J1=O1O1+O2,J2=O2O1+O2.We further summarize all the parameters in Section 3.4 into Table 4.Table 4
System
Parameters
Initial states
Fuzzy Rossler’s system (functional system)
h
=
0.2, i=0.2, j=5.7, d4=10
(w1, w2, w3) = (10, 15, 11)Figure 5
Chaotic behavior of fuzzy Rossler’s system as functional system.


## 3.1. Fuzzy Modeling of Master Lorenz System

For master Lorenz system(30)x˙1=ax2-x1,x˙2=cx1-x1x3-x2,x˙3=x1x2-bx3,where a, b, c are the parameters. When a=10, b=8/3, c=28, and initial states are (0.5,1,5), the dynamic behavior is chaotic. Assume that x1∈[-d1,d1] and d1>0; then Lorenz system can be exactly represented by T-S fuzzy model as follows:(31)Rule1:IFxisM1,THENX˙t=A1Xt,(32)Rule2:IFxisM2,THENX˙t=A2Xt,where(33)X=x1,x2,x3T,A1=-aa0c-1-d10d1-b,A2=-aa0c-1d10-d1-b,M1x=121+x1d1,M2x=121-x1d1,and d1=30. M1 and M2 are fuzzy set of Lorenz system. Here, we call (31) the first liner subsystem under the fuzzy rule and (32) the second liner subsystem under the fuzzy rule. The final output of the fuzzy Lorenz system is inferred as follows and the chaotic behavior is shown in Figure 2(34)X˙t=∑i=12HiAiXt,where(35)H1=M1M1+M2,H2=M2M1+M2.We further summarize all the parameters in Section 3.1 into Table 1.Table 1
System
Parameters
Initial states
Fuzzy Lorenz system (master)
a
=
10, b=8/3, c=28, d1=30
(x1, x2, x3) = (0.5, 1, 5)Figure 2
Chaotic behavior of master fuzzy Lorenz system.


## 3.2. Fuzzy Modeling of Slave Lorenz System

For slave Lorenz system(36)y˙1=ay2-y1+u1,y˙2=cy1-y1y3-y2+u2,y˙3=y1y2-by3+u3,where a, b, c are the parameters. When a=10, b=8/3, c=28, and initial states are (0.6,3,10), the dynamic behavior is chaotic. Assume that y1∈[-d2,d2] and d2>0; then Lorenz system can be exactly represented by T-S fuzzy model as follows:(37)Rule1:IFyisM~1,THENY˙t=A~1Yt+U1,(38)Rule2:IFyisM~2,THENY˙t=A~2Yt+U2,where(39)Y=y1,y2,y3T,A~1=-aa0c-1-d20d2-b,A~2=-aa0c-1d20-d2-b,M~1x=121+y1d2,M~2x=121-y1d2,U1=u11,u12,u13T,U2=u21,u22,u23Tand d2=20. M~1 and M~2 are fuzzy set of Lorenz system. Here, we call (37) the first linear subsystem under the fuzzy rule and (38) the second linear subsystem under the fuzzy rule. As a result, U1 and U2 are controllers of the first and second subsystems. The final output of the fuzzy Lorenz system is inferred as follows and the chaotic behavior without using controllers is shown in Figure 3(40)Y˙t=∑i=12H~iA~iYt+Ui,where(41)H~1=M~1M~1+M~2,H~2=M~2M~1+M~2.We further summarize all the parameters in Section 3.2 into Table 2.Table 2
System
Parameters
Initial states
Fuzzy Lorenz system (slave)
a
=
10, b=8/3, c=28, d2=20
(y1, y2, y3) = (0.6, 3, 10)Figure 3
Chaotic behavior of slave fuzzy Lorenz system.


## 3.3. Fuzzy Modeling of Chen System as Slave System

For Chen system as slave system(42)z˙1=ez2-z1+u1,z˙2=g-ez1-z1z3+gz2+u2,z˙3=z1z2-fz3+u3,where e, f, g are the parameters. When e=35, f=3, g=28, and initial states are (0.5,1,5), the dynamic behavior is chaotic. Assume that z1∈[-d3,d3] and d3>0; then Chen system can be exactly represented by T-S fuzzy model as follows:(43)Rule1:IFzisN1,THENZ˙t=B1Zt,(44)Rule2:IFzisN2,THENZ˙t=B2Zt.Here(45)Z=z1,z2,z3T,B1=-ee0g-eg-d30d3-f,B2=-ee0g-egd30-d3-f,N1z=121+z1d3,N2z=121-z1d3,U1=u11,u12,u13T,U2=u21,u22,u23Tand d3=30. N1 and N2 are fuzzy set of Lorenz system. Here, we call (43) the first liner subsystem under the fuzzy rule and (44) the second liner subsystem under the fuzzy rule. As a result, U1 and U2 are controllers of the first and second subsystems. The final output of the fuzzy Chen system is inferred as follows and the chaotic behavior without using controllers is shown in Figure 4(46)Z˙t=∑i=12IiBiZt+Ui,where(47)I1=N1N1+N2,I2=N2N1+N2.We further summarize all the parameters in Section 3.3 into Table 3.Table 3
System
Parameters
Initial states
Fuzzy Chen system (slave)
e
=
35, f=3, g=28, d3=30
(z1, z2, z3) = (0.5, 1, 5)Figure 4
Chaotic behavior of slave fuzzy Chen system.


## 3.4. Fuzzy Modeling of Rossler’s System as Functional System

For Rossler’s system as functional system(48)w˙1=-w2-w3,w˙2=w1+hw2,w˙3=i+w1w3-jw3,where e, f, g are the parameters. When h=0.2, i=0.2, j=5.7, and initial states are (10,15,11), the dynamic behavior is chaotic. Assume that w1∈[-d4,d4] and d4>0; then Chen system can be exactly represented by T-S fuzzy model as follows:(49)Rule1:IFwisO1,THENW˙t=C1Wt+c1,(50)Rule2:IFwisO2,THENW˙t=C2Wt+c2.Here(51)W=w1,w2,w3T,C1=0-1-11h00d4-j,C2=0-1-11h00-d4-j,c1=00i,c2=00iO1z=121+w1d4,O2z=121-w1d4,and d4=10. O1 and O2 are fuzzy set of Lorenz system. Here, we call (49) the first liner subsystem under the fuzzy rule and (50) the second liner subsystem under the fuzzy rule. The final output of the fuzzy Rossler’s system is inferred as follows and the chaotic behavior is shown in Figure 5(52)Z˙t=∑i=12JiCiZt+ci,where(53)J1=O1O1+O2,J2=O2O1+O2.We further summarize all the parameters in Section 3.4 into Table 4.Table 4
System
Parameters
Initial states
Fuzzy Rossler’s system (functional system)
h
=
0.2, i=0.2, j=5.7, d4=10
(w1, w2, w3) = (10, 15, 11)Figure 5
Chaotic behavior of fuzzy Rossler’s system as functional system.


## 4. Numerical Simulation Results of the Bioinspired Strategy

In this section, there are three mainCASEs for illustrations, Case 1: the generalized synchronization error function is ei=(xi-yi+100), (i=1,2,3), Case 2: the generalized synchronization error function is ei=(xi-zi+100), (i=1,2,3), and Case 3: the generalized synchronization error function is ei=xi-zi+wi+100, (i=1,2,3). A set of simple-form and efficient Lyapunov function based on the natural Bioinspired stability theory is directly used to achieve generalized synchronization of two fuzzy chaotic systems.Case 1.
The generalized synchronization error function isei=(xi-yi+100), (i=1,2,3).
The addition of 100 makes the error dynamics always happen in first quadrant. Our goal isyi=xi+100; that is, (54)limt→∞⁡ei=limt→∞⁡xi-yi+100=0i=1,2,3.
The error and error dynamics are(55)e=e1e2e3=x1-y1+100x2-y2+100x3-y3+100,e˙=e˙1e˙2e˙3=x˙1-y˙1x˙2-y˙2x˙3-y˙3=∑i=12HiAiXt-∑i=12H~iA~iYt+Ui=H1ax2-x1-H~1ay2-y1+u11+H2ax2-x1-H~2ay2-y1+u21H1cx1-d1x3-x2-H~1cy1-d2y3-y2+u12+H2cx1+d1x3-x2-H~2cy1+d2y3-y2+u22H1d1x2-bx3-H~1d2y2-by3+u13+H2-d1x2-bx3-H~2-d2y2-by3+u23.
Before control action, the error dynamics always happens in first quadrant as shown in Figure6. By the natural bioinspired stability theory, one can choose a Lyapunov function in the form of a positive definite function in first quadrant:(56)V=e1+e2+e3>0.
Its time derivative through error dynamics (55) is(57)V˙=e˙1+e˙2+e˙3=H1ax2-x1-H~1ay2-y1+u11+H2ax2-x1-H~2ay2-y1+u21+H1cx1-d1x3-x2-H~1cy1-d2y3-y2+u12+H2cx1+d1x3-x2-H~2cy1+d2y3-y2+u22+H1d1x2-bx3-H~1d2y2-by3+u13+H2-d1x2-bx3-H~2-d2y2-by3+u23.Choose (58)u11=H1H~1ax2-x1-ay2-y1+e1,u12=H1H~1-d1x3+cx2--d2y3+cy2+e2,u13=H1H~1d1x2-bx3-d2y2-by3+e3,u21=H2H~2ax2-x1-ay2-y1+e1,u22=H2H~2d1x3+cx2-d2y3+cy2+e2,u23=H2H~2-d1x2-bx3--d2y2-by3+e3.We obtain(59)V˙=-e1-e2-e3<0,which is negative definite function in the first quadrant. The numerical results are shown in Figures 7 and 8. In Figures 7 and 8, the master and slave systems are two identical Lorenz systems with different initial conditions, where the parameters a=10, b=8/3, c=28, d1=30, d2=20, initial states are (x1,x2,x3)=(0.5,1,5) for master system, and (y1,y2,y3)=(0.6,3,10) for slave system. The control procedure is arranged that, in the initial 30 sec, the system is in a nature situation without any control. The designed controllers operate after 30 sec to show the high efficiency of our bioinspired controllers. In Figure 7, it is clear that the error states e1, e2, e3 almost immediately smoothly approach zero after the controllers put in at 30 sec. Also, in Figure 8, the states of the slave system are all synchronized to the states of the target system with high performance.Figure 6
Phase portraits of error dynamics for Case1.Figure 7
Time histories of error dynamics for Case1—the error states e1, e2, e3 almost immediately smoothly approach zero after the controllers put in at 30 sec.Figure 8
Time histories ofx1, x2, x3 and y1, y2, y3 for Case 1—the states of the slave system are all synchronized to the states of the target system with high performance.Case 2.
The generalized synchronization error function isei=(xi-zi+100), (i=1,2,3).
The addition of 100 makes the error dynamics always happen in first quadrant. Our goal iszi=xi+100; that is, (60)limt→∞⁡ei=limt→∞⁡xi-zi+100=0i=1,2,3.
The error and error dynamics are(61)e=e1e2e3=x1-z1+100x2-z2+100x3-z3+100,e˙=e˙1e˙2e˙3=x˙1-z˙1x˙2-z˙2x˙3-z˙3=∑i=12HiAiXt-∑i=12IiBiZt+Ui=H1ax2-x1-I1ez2-z1+u11+H2ax2-x1-I2ez2-z1+u21H1cx1-d1x3-x2-I1g-ez1-d3z3+gz2+u12+H2cx1+d1x3-x2-I2g-ez1+d3z3+gz2+u22H1d1x2-bx3-I1d3z2-fz3+u13+H2-d1x2-bx3-I2-d3z2-fz3+u23.
Before control action, the error dynamics always happens in first quadrant as shown in Figure9. By the natural bioinspired stability theory, one can choose a Lyapunov function in the form of a positive definite function in first quadrant:(62)V=e1+e2+e3>0.Its time derivative through error dynamics (61) is(63)V˙=e˙1+e˙2+e˙3=H1ax2-x1-I1ez2-z1+u11+H2ax2-x1-I2ez2-z1+u21+H1cx1-d1x3-x2-I1g-ez1-d3z3+gz2+u12+H2cx1+d1x3-x2-I2g-ez1+d3z3+gz2+u22+H1d1x2-bx3-I1d3z2-fz3+u13+H2-d1x2-bx3-I2-d3z2-fz3+u23.Choose (64)u11=H1I1ax2-x1-ez2-z1+e1,u12=H1I1cx1-d1x3-x2-g-ez1-d3z3+gz2+e2,u13=H1I1d1x2-bx3-d3z2-fz3+e3,u21=H2I2ax2-x1-ez2-z1+e1,u22=H2I2cx1+d1x3-x2-g-ez1+d3z3+gz2+e2,u23=H2I2-d1x2-bx3--d3z2-fz3+e3.We obtain(65)V˙=-e12-e22-e32<0,which is negative definite function in the first quadrant. The numerical results are shown in Figures 10 and 11. In Figures 10 and 11, Lorenz system with initial state (x1,x2,x3)=(0.5,1,5) in (34) is considered as master system and Chen system with initial state (z1,z2,z3)=(0.5,1,5) in (46) is referring to the slave system, where the parameters a=10, b=8/3, c=28, e=35, f=3, g=28, d1=30, d3=30. The control procedure is arranged that, in the initial 30 sec, the system is in a nature situation without any control. The designed controllers operate after 30 sec to show the high efficiency of our bioinspired controllers. In Figure 10, it is clear that the error states e1, e2, e3 almost immediately smoothly approach zero after the controllers put in at 30 sec. Also, in Figure 11, the states of the slave system are all synchronized to the states of the target system with high performance.Figure 9
Phase portraits of error dynamics for Case2.Figure 10
Time histories of error dynamics for Case2—the error states e1, e2, e3 almost immediately smoothly approach zero after the controllers put in at 30 sec.Figure 11
Time histories ofx1+100, x2+100, x3+100, and z1, z2, z3 for Case 2—the states of the slave system are all synchronized to the states of the target system with high performance.Case 3.
The generalized synchronization error function isei=xi-zi+wi+100, (i=1,2,3), where the states (z1,z2,z3) of the slave Chen system in (46) are going to achieve generalized synchronization of a complicated designed goal system—the states (x1,x2,x3) of Lorenz system in (32) with the states (w1,w2,w3) of the functional system in (52).
The addition of 100 makes the error dynamics always happen in first quadrant. Our goal iszi=xi+wi+100; that is, (66)limt→∞⁡ei=limt→∞⁡xi-zi+wi+100=0i=1,2,3.
The error and error dynamics are(67)e=e1e2e3=x1-z1+w1+100x2-z2+w2+100x3-z3+w3+100,e˙=e˙1e˙2e˙3=x˙1-w˙1x˙2-w˙2x˙3-w˙3=∑i=12HiAiXt+∑i=12JiCiZt+ci-∑i=12IiBiZt+Ui=H1ax2-x1+J1-w2-w3-I1ez2-z1+u11+H2ax2-x1+J2-w2-w3-I2ez2-z1+u21H1cx1-d1x3-x2+J1w1+hw2-I1g-ez1-d3z3+gz2+u12+H2cx1+d1x3-x2+J2w1+hw2-I2g-ez1+d3z3+gz2+u22H1d1x2-bx3+J1i+d4w3-jw3-I1d3z2-fz3+u13+H2-d1x2-bx3+J2i-d4w3-jw3-I2-d3z2-fz3+u23.
Before control action, the error dynamics always happens in first quadrant as shown in Figure12. By the natural bioinspired stability theory, one can choose a Lyapunov function in the form of a positive definite function in first quadrant:(68)V=e1+e2+e3>0.Its time derivative through error dynamics (67) is(69)V˙=e˙1+e˙2+e˙3=H1ax2-x1+J1-w2-w3-I1ez2-z1+u11+H2ax2-x1+J2-w2-w3-I2ez2-z1+u21+H1cx1-d1x3-x2+J1w1+hw2-I1g-ez1-d3z3+gz2+u12+H2cx1+d1x3-x2+J2w1+hw2-I2g-ez1+d3z3+gz2+u22+H1d1x2-bx3+J1i+d4w3-jw3-I1d3z2-fz3+u13+H2-d1x2-bx3+J2i-d4w3-jw3-I2-d3z2-fz3+u23.Choose (70)u11=H1I1ax2-x1+J1I1-w2-w3-ez2-z1+e1,u12=H1I1cx1-d1x3-x2+J1I1w1+hw2-g-ez1-d3z3+gz2+e2,u13=H1I1d1x2-bx3+J1I1i+d4w3-jw3-d3z2-fz3+e3,u21=H2I2ax2-x1+J2I2-w2-w3-ez2-z1+e1,u22=H2I2cx1+d1x3-x2+J2I2w1+hw2-g-ez1+d3z3+gz2+e2,u23=H2I2-d1x2-bx3+J2I2i-d4w3-jw3--d3z2-fz3+e3.We obtain(71)V˙1=-e1-e2-e3<0,which is negative definite function in the first quadrant. The numerical results are shown in Figures 13 and 14. In Figures 13 and 14, a complicated goal system—Lorenz system with initial state (x1,x2,x3)=(0.5,1,5) in (34) combined with the functional system with initial state (w1,w2,w3)=(10,15,11) in (52), is considered as the master system and Chen system with initial state (z1,z2,z3)=(0.5,1,5) in (46) is referring to the slave system, where the parameters a=10, b=8/3, c=28, e=35, f=3, g=28, h=0.2, i=0.2, j=5.7, d1=30, d3=30, d4=10. The control procedure is arranged that, in the initial 30 sec, the system is in a nature situation without any control. The designed controllers operate after 30 sec to show the high efficiency of our bioinspired controllers. In Figure 12, it is clear that the error states e1, e2, e3 almost immediately smoothly approach zero after the controllers put in at 30 sec. Also, in Figure 13, the states of the slave system are all synchronized to the states of the target system with high performance.Figure 12
Phase portraits of error dynamics for Case3.Figure 13
Time histories of error dynamics for Case3—the error states e1, e2, e3 almost immediately smoothly approach zero after the controllers put in at 30 sec.Figure 14
Time histories ofx1+w1+100, x2+w3+100, x3+w3+100, and z1, z2, z3 for Case 3—the states of the slave system are all synchronized to the states of the target system with high performance.


## 5. Conclusions

In our daily life, there are lots of natural systems whose certain states are always positive, that is, exist in the first quadrant. Consequently, a novel control strategy by using bioinspired stability theory is proposed to achieve generalized fuzzy chaos synchronization. In this paper, a new Lyapunov function based on bioinspired control strategy is proposed to directly achieve chaos synchronization of chaotic systems. Via using this strategy, the new Lyapunov function used is a simple linear homogeneous function of states and the process of synchronization of two fuzzy chaotic systems becomes more simple and the controllers operate after 30 sec; the error statese1, e2, e3 almost immediately smoothly approach zero. Classical Lorenz, Chen, and Rossler’s systems are used in numerical examples and are given to demonstrate the effectiveness of the proposed new strategy. Also, this paper provides a new sight of view for researchers to reference and we are going to further apply the novel bioinspired controlling strategy to other real-world applications in the following stage, such as energy communication system in biological research, virus-immune system, and HIV virus-immune system.



---

*Source: 725674-2015-04-30.xml*