# Synchronization of Discrete-Time Chaotic Fuzzy Systems by means of Fuzzy Output Regulation Using Genetic Algorithm

Authors: Tonatiuh  Hernández Cortés; A. Verónica  Curtidor López; Jorge  Rodríguez Valdez; Jesús A.  Meda Campaña; Ricardo  Tapia Herrera; José de Jesús  Rubio
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

The synchronization of chaotic systems, described by discrete-time T-S fuzzy models, is treated by means of fuzzy output regulation theory. The conditions for designing a discrete-time output regulator are given in this paper. Besides, when the system does not fulfill the conditions for exact tracking, a new regulator based on genetic algorithms is considered. The genetic algorithms are used to approximate the adequate membership functions, which allow the adequate combination of local regulators. As a result, the tracking error is significantly reduced. Both the Complete Synchronization and the Generalized Synchronization problem are studied. Some numerical examples are used to illustrate the effectiveness of the proposed approach.


---

## Body

## 1. Introduction

A special nonlinear dynamical phenomenon, known as chaos, emerged in mid-1960s and reached applicable technology in the late 1990s and was considered as one of the three monumental discoveries of the twentieth century. On the other hand, fuzzy logic, a set theory and then an infinite-valued logic, gets a wide applicability in many industrial, commercial, and technical fields, ranging from control, automation, and artificial intelligence, just to name a few. Fuzzy logic and chaos had been considered by many researches and engineers as fundamental concepts and theories and their broad applicability in technology as well. The interaction between fuzzy logic and chaos has been developed for the last 20 years leading to research topics as fuzzy modeling of chaotic systems using Takagi-Sugeno models, linguistic descriptions of chaotic systems, fuzzy control of chaos, synchronization, and a combination of fuzzy chaos for engineering applications [1, 2].In the 1960s, Rechenberg [3] introduced “evolution strategies,” a method to optimize real-valued parameters for devices such as airfoils. This idea was further developed by Schwefel in [4]. Genetic algorithms (GAs) were initially developed by Bremermann [5] in 1958 but popularized and developed by Holland in the 1960s. In contrast with evolution strategies and evolutionary programming, Holland's idea was not to design algorithms to solve specific problems but rather to formally study the phenomenon of adaptation, as it occurs in nature, and develop ways in which the mechanisms of natural adaptation might be transferred into computer systems [6]. The genetic algorithm is presented as abstraction of biological evolution and theoretical framework for adaptation for moving from one population of “chromosomes” (e.g., strings of ones and zeros, or “bits”) to a new population by using a kind of “natural selection” together with the genetics-inspired operators of crossover, mutation, and inversion. Each chromosome consists of “genes” (e.g., bits); each gene is being an instance of a particular “allele” (e.g., 0 or 1). The operator selection chooses those chromosomes in the population that will be allowed to reproduce and those adjusted chromosomes produce more offspring than the less ones [7].According to Fogel and Aderson [8], Bremermann was the first to implement real-coded genetic algorithms as well as providing a mathematic model of GA known as the one-max function. In contrast to genetic algorithms, Evolutionary Strategies were initially developed for the purpose of parameter optimization. The idea was to imitate the principles of organic evolution in experimental parameter optimization for applications such as pipe bending or PID control for a nonlinear system [9].Synchronization of chaotic systems is one of the more exiting problems in control science and can be referred at least to Huygen’s observations [10]; it is understood as one of the trajectories of two autonomous chaotic systems, starting from nearly initial conditions and converging to the other, and remains as t→∞; in [11] it was reported that some kind of chaotic systems possesses a self-synchronization property. However, not all chaotic systems can be decomposed in two separate responses subsystems and be able to synchronize the drive system. The ideas of these works have led to improvement in many fields, such as communications [12], encrypted systems, the complex information processing within the human brain, coupled biochemical reactors, and earthquake engineering [13].Synchronization can be classified as follows: Complete Synchronization: it is when two identical chaos oscillators are mutually coupled and one drives to the other; Generalized Synchronization: it differs from the previous case by the fact that there are different chaos oscillators and the states of one are completely defined by the other; Phase Synchronization: it occurs when the coupled oscillators are not identical and have different amplitude that is still unsynchronized, while the phases of oscillators evolve to be synchronized [14]. It is worth mentioning that studies in synchronization of nonlinear systems have been reformulated based on the previous results from classical control theory such as [15–18].In this paper, the fuzzy output regulation theory and Takagi-Sugeno (T-S) fuzzy models are used to solve the Complete and Generalized Synchronization by using linear local regulators. Isidori and Byrnes [19] showed that the output regulation established by Francis could be extended for a nonlinear sector as a general case, resulting in a set of nonlinear partial differential equations called Francis-Isidori-Byrnes (FIB). Unfortunately these equations in many cases are too difficult to solve in a practical manner. For this reason in [20] the approach based on the weighted summation of local linear regulators is presented and in [21] the new membership functions in the regulator are approximated by soft computing techniques.So, the main contribution of the present work is to find a control law for synchronizing of chaotic systems described by discrete-time Takagi-Sugeno fuzzy models, first when the system fulfills the following: (1) the input matrix for all subsystems is the same and (2) the local regulators share the same zero error manifold π(wk). In this way, the results given in [20] are extended to the discrete-time domain. On the other hand, when the systemmaster-slave does not fulfill the aforementioned conditions, new membership functions are computed in order to enhance the performance of the fuzzy regulator. Such proposed membership functions are different from those given in the plant or exosystem and are tuned by using the GA. The tuning of the new membership functions, which is as generalized bell-shaped function, is given by optimization of the form parameter.The rest of the paper is organized as follows. In Section2 the discrete-time output regulation problem formulation is given with a brief review of the Takagi-Sugeno models and the discrete-time fuzzy regulation problem. In Section 3 the tuning of membership functions by means of GAs is thoroughly discussed. In Section 4 Complete and Generalized Synchronization with some examples are presented and finally, in Section 5, some conclusions are drawn.


## 2. The Discrete-Time Fuzzy Output Regulation Problem

Consider a nonlinear discrete-time system defined by(1)xk+1=fxk,ωk,uk,(2)yk=cxk,(3)ωk+1=sωk,(4)yref,k=qωk,(5)ek=hxk,ωk,where xk∈Rn is the state vector of the plant, wk∈W⊂Rs is the state vector of the exosystem, which generates the reference and/or the perturbation signals, and uk∈Rm is the input signal. Equation (5) refers to difference between output system of the plant (yk∈Rm) and the reference signal (yref,k∈Rm); that is, h(xk,ωk)=yk-yref,k=c(xk)-q(wk) and take into account that m≤n. Besides, it is assumed that f(xk,uk,wk), h(xk,wk), and s(wk) are analytical functions and also that f(0,0,0)=0, s(0)=0, and h(0,0)=0 [22].Clearly, by linearizing (1)–(5) around xk=0, one gets(6)xk+1=Axk+Buk+Pwk,yk=Cxk,wk+1=Swk,yref,k=Qwk,ek=Cxk-Qwk.Thus, the Nonlinear Regulator Problem [19, 23] consists of finding a controller uk=α(xk,wk), such that the closed-loop system xk+1=Axk+Bα(xk,0) has an asymptotically stable equilibrium point, and the solution of system (6) satisfies limk→∞ek=0.So, by definingπ(wk) as the steady-state zero error manifold and γ(wk) as the steady-state input, the following theorem gives the conditions for the solution of nonlinear regulation problem.Theorem 1.
Suppose thatwk+1=s(wk) is Poisson stable and there exists a gain K such that the matrix A+BK is stable and there exist mappings xss(t)=π(wk) and uss=γ(wk) with π(0)=0 and γ(0)=0 satisfying(7)πswk=fπwk,wk,γwk,0=hπwk,wk.Then the control signal for the nonlinear regulation is given by(8)uk=Kxk-πwk+γwk.Proof.
See [22–24].The equation set (7) is known as Discrete-Time Francis-Isidori-Byrnes (DTFIB) equations and linear counterpart is obtained when the mappings xss,k=π(wk) and uss,k=γ(wk) transform into xss,k=Πwk and uss,k=Γwk, respectively. Thus, the problem is reduced to solve linear matrix equations [25] given by(9)ΠS=AΠ+BΓ+P,0=CΠ-Q.


### 2.1. The Discrete-Time Output Fuzzy Regulation Problem

Takagi and Sugeno proposed a fuzzy model composed of a set of linear subsystems with IF-THEN rules capable of relating physical knowledge, linguistic characteristics, and properties of the system. Such a model successfully represents a nonlinear system at least in a predefined region of phase space [15]. The T-S model for the plant and exosystem can be described as follows [26]:Plant Model
Rulei:
IFz1,1,k is M1,1i and … and z1,p1,k is M1,p1i,(10)THENxk+1=Aixk+Biuk+Piwk,yk=Cixk,i=1,2,…,r1,where r1 is the number of rules in the model of the plant and the sets M1,ji are the fuzzy sets defined based on the previous dynamic knowledge of the system.Exosystem Model
Rulei:
IF
z
2,1
,
k is M2,1i and … and z2,p2,k is M2,p2i,(11)THENwk+1=Siwk,yref,k=Qiwk,i=1,2,…,r2,where r2 is the number of rules in the model of the exosystem and M2,ji are the fuzzy sets.Then, the regulation problem defined by (1)–(5) can be represented through the T-S discrete-time fuzzy model; that is, [20](12)xk+1=∑i=1r1h1,iz1,kAixk+Biuk+Piwk,wk+1=∑i=1r2h2,iz2,kSiwk,ek=∑i=1r1h1,iz1,kCixk-∑i=1r2h2,iz2,kQiwk,where xk∈Rn is the state vector of the plant, wk∈Rs is the state vector of the exosystem, uk∈Rm is the input signal, ek∈Rm, and h∗,i(z) is the normalized weight of each rule, 1 for the plant and 2 for the exosystem, which depends on the membership function for the premise variable z∗,k in M∗,ji; that is,(13)ω¯∗,iz∗,k=∏j=1p∗M∗,jiz∗,j,k,h∗,iz∗,k=ω¯∗,iz∗,k∑i=1r∗ω¯∗,iz∗,k,∑i=1r∗h∗,iz∗,k=1,h∗,iz∗,k≥0with z∗,k=[z∗,1,kz∗,2,k⋯z∗,p∗,k] as a function of xk and/or wk, i=1,…,r∗ and j=1,…,p.Thediscrete-time fuzzy output regulation problem consists of finding a controller uk=α(xk,wk), such that the closed-loop system with no external signal (14)xk+1=∑i=1r1h1,iz1,kAixk+Biαxk,0has an asymptotically stable equilibrium point.The solution of system (12) satisfies(15)limk→∞⁡ek=0.In order to achieve the synchronization of chaotic systems described by a T-S discrete-time fuzzy model it is necessary to fulfill (7) [27, 28]. Then(16)πswk=∑i=1r1h1,iz1,kAiπwk+Biγwk+Piwk,(17)0=∑i=1r1h1,iz1,kCiπwk-∑i=1r2h2,iz2,kQiwk,where π(wk) is the zero error steady-state manifold which becomes invariant by the effect of the steady-state input γ(wk).Assuming the mappingsπ(wk) and γ(wk) as(18)π~wk=∑i=1r1h1,iz1,k∑j=1r2h2,iz2,kΠi,jwk,(19)γ~wk=∑i=1r1h1,iz1,k∑j=1r2h2,iz2,kΓi,jwk,respectively, with Πi,j and Γi,j as a solution of r1·r2 lineal local problems,(20)ΠijSj=AiΠij+BiΓij+Pi,(21)0=CiΠi,j-Qi,for all i=1,…,r1 and j=1,…,r2, the following control law can be obtained [20, 22, 23]:(22)uk=∑h=1r1h1,iz1,kKixk-∑i=1r1h1,iz1,k∑j=1r2h2,jz2,kΠijwk+∑i=1r1h1,iz1,k∑j=1r2h2,jz2,kΓijwk.However, by substitution of π(wk) and γ(wk) in (16) and (17) and considering(1)
the steady-state zero error manifoldπ(wk)=Πwk, that is, Πij=Π,(2)
the input matricesBi=B and/or Γij=Γ,for all i=1,…,r1 and j=1,…,r2, the following control signal uk emerges:(23)uk=∑h=1r1h1,iz1,kKixk-Πwk+∑i=1r1h1,iz1,k∑j=1r2h2,jz2,kΓijwk.On the other hand, the existence of a fuzzy stabilizer of the formu=∑i=1r1h1,i(z1,k)Kixk, ensuring that the tracking error converges asymptotically to zero, can be obtained from the Parallel Distributed Compensator (PDC) [29, 30] or another stability analysis for T-S fuzzy models such as [31].Remark 2.
The control signal in (22) is given by the substitution of (18) and (19) in (16) and (17); the proposed controller provides the following advantages:(1)
All parameters included in the controller are known; this includes the membership functions of the plant and exosystem, which are well defined in the T-S fuzzy model.On the other hand,Πij and Γij come directly from solving of r1∗r2 local linear problems equivalent to solving the Francis equations; such problems can be easily solved by using programs like Matlab or Mathematica.(2)
In the case when ther1∗r2 local linear problems lead to Πij≠Π, then at least a bounded error is ensured.(3)
It is clear that whenΠij=Π, the term ∑i=1r1h1,i(z1,k)∑j=1r2h2,j(z2,k)Πijwk changes to Πwk leading to controller defined in (23).(4)
The following condition: the input matricesBi=B and/or Γij=Γ avoids the crossed terms in the solutions of (16) and (17), allowing the exact fuzzy output regulation.(5)
The proposed controller can be seen as a simple substitution of the aforementioned elements.
On the other hand, the following disadvantages can appear:(1)
If the condition is that the steady-state zero error manifoldπ(wk)≠Πwk, that is, Πij≠Π, then, it will be necessary to adjust the local regulator by means of new membership functions. Please refer to Section 3.(2)
As expected, the complexity of the controller increases according to the number of local subsystems.The following theorem provides the conditions for the existence of the exact fuzzy output regulator for a discrete-time T-S fuzzy models.Theorem 3.
The exact fuzzy output regulation with full information of systems defined as (12) is solvable if (a) there exists the same zero error steady-state manifold π(wk)=Π; (b) there exist uk=∑i=1r1hi,k(z1,k)Kixk for the fuzzy system; (c) the exosystem ωk+1=s(ωk) is Poisson stable, and the input matrices for all subsystems Bi are equal. Moreover, theExact Output Fuzzy Regulation Problem is solvable by the controller(24)uk=∑h=1r1h1,iz1,kKixk-Πwk+∑i=1r1h1,iz1,k∑j=1r2h2,jz2,kΓijwk.Proof.
From the previous analysis, the existence of mappingsπ(wk)=Πwk and γ(wk)=∑i=1r1h1,i(z1,k)∑j=1r2h2,j(z2,k)Γijwk is guaranteed when the input matrices for all subsystems Bi are equal, and the solution of r1·r2 lineal local problems is(25)ΠijSj=AiΠij+BiΓij+Pi,0=CiΠi,j-Qi,leading to Πij=Π for all i=1,…,r1 and j=1,…,r2.
On the other hand, the inclusion of condition (b) has been thoroughly discussed in [22, 23, 27, 28, 31–33], and it implies the existence of a fuzzy stabilizer.
Condition (c) ensures the nonexistence of crossed terms in the local Francis equations. Finally, condition (d) is introduced to avoid the fact that the reference signal converges to zero, which would turn the regulation problem into a simple stability problem. The rest of the proof follows directly from the previous analysis.


## 2.1. The Discrete-Time Output Fuzzy Regulation Problem

Takagi and Sugeno proposed a fuzzy model composed of a set of linear subsystems with IF-THEN rules capable of relating physical knowledge, linguistic characteristics, and properties of the system. Such a model successfully represents a nonlinear system at least in a predefined region of phase space [15]. The T-S model for the plant and exosystem can be described as follows [26]:Plant Model
Rulei:
IFz1,1,k is M1,1i and … and z1,p1,k is M1,p1i,(10)THENxk+1=Aixk+Biuk+Piwk,yk=Cixk,i=1,2,…,r1,where r1 is the number of rules in the model of the plant and the sets M1,ji are the fuzzy sets defined based on the previous dynamic knowledge of the system.Exosystem Model
Rulei:
IF
z
2,1
,
k is M2,1i and … and z2,p2,k is M2,p2i,(11)THENwk+1=Siwk,yref,k=Qiwk,i=1,2,…,r2,where r2 is the number of rules in the model of the exosystem and M2,ji are the fuzzy sets.Then, the regulation problem defined by (1)–(5) can be represented through the T-S discrete-time fuzzy model; that is, [20](12)xk+1=∑i=1r1h1,iz1,kAixk+Biuk+Piwk,wk+1=∑i=1r2h2,iz2,kSiwk,ek=∑i=1r1h1,iz1,kCixk-∑i=1r2h2,iz2,kQiwk,where xk∈Rn is the state vector of the plant, wk∈Rs is the state vector of the exosystem, uk∈Rm is the input signal, ek∈Rm, and h∗,i(z) is the normalized weight of each rule, 1 for the plant and 2 for the exosystem, which depends on the membership function for the premise variable z∗,k in M∗,ji; that is,(13)ω¯∗,iz∗,k=∏j=1p∗M∗,jiz∗,j,k,h∗,iz∗,k=ω¯∗,iz∗,k∑i=1r∗ω¯∗,iz∗,k,∑i=1r∗h∗,iz∗,k=1,h∗,iz∗,k≥0with z∗,k=[z∗,1,kz∗,2,k⋯z∗,p∗,k] as a function of xk and/or wk, i=1,…,r∗ and j=1,…,p.Thediscrete-time fuzzy output regulation problem consists of finding a controller uk=α(xk,wk), such that the closed-loop system with no external signal (14)xk+1=∑i=1r1h1,iz1,kAixk+Biαxk,0has an asymptotically stable equilibrium point.The solution of system (12) satisfies(15)limk→∞⁡ek=0.In order to achieve the synchronization of chaotic systems described by a T-S discrete-time fuzzy model it is necessary to fulfill (7) [27, 28]. Then(16)πswk=∑i=1r1h1,iz1,kAiπwk+Biγwk+Piwk,(17)0=∑i=1r1h1,iz1,kCiπwk-∑i=1r2h2,iz2,kQiwk,where π(wk) is the zero error steady-state manifold which becomes invariant by the effect of the steady-state input γ(wk).Assuming the mappingsπ(wk) and γ(wk) as(18)π~wk=∑i=1r1h1,iz1,k∑j=1r2h2,iz2,kΠi,jwk,(19)γ~wk=∑i=1r1h1,iz1,k∑j=1r2h2,iz2,kΓi,jwk,respectively, with Πi,j and Γi,j as a solution of r1·r2 lineal local problems,(20)ΠijSj=AiΠij+BiΓij+Pi,(21)0=CiΠi,j-Qi,for all i=1,…,r1 and j=1,…,r2, the following control law can be obtained [20, 22, 23]:(22)uk=∑h=1r1h1,iz1,kKixk-∑i=1r1h1,iz1,k∑j=1r2h2,jz2,kΠijwk+∑i=1r1h1,iz1,k∑j=1r2h2,jz2,kΓijwk.However, by substitution of π(wk) and γ(wk) in (16) and (17) and considering(1)
the steady-state zero error manifoldπ(wk)=Πwk, that is, Πij=Π,(2)
the input matricesBi=B and/or Γij=Γ,for all i=1,…,r1 and j=1,…,r2, the following control signal uk emerges:(23)uk=∑h=1r1h1,iz1,kKixk-Πwk+∑i=1r1h1,iz1,k∑j=1r2h2,jz2,kΓijwk.On the other hand, the existence of a fuzzy stabilizer of the formu=∑i=1r1h1,i(z1,k)Kixk, ensuring that the tracking error converges asymptotically to zero, can be obtained from the Parallel Distributed Compensator (PDC) [29, 30] or another stability analysis for T-S fuzzy models such as [31].Remark 2.
The control signal in (22) is given by the substitution of (18) and (19) in (16) and (17); the proposed controller provides the following advantages:(1)
All parameters included in the controller are known; this includes the membership functions of the plant and exosystem, which are well defined in the T-S fuzzy model.On the other hand,Πij and Γij come directly from solving of r1∗r2 local linear problems equivalent to solving the Francis equations; such problems can be easily solved by using programs like Matlab or Mathematica.(2)
In the case when ther1∗r2 local linear problems lead to Πij≠Π, then at least a bounded error is ensured.(3)
It is clear that whenΠij=Π, the term ∑i=1r1h1,i(z1,k)∑j=1r2h2,j(z2,k)Πijwk changes to Πwk leading to controller defined in (23).(4)
The following condition: the input matricesBi=B and/or Γij=Γ avoids the crossed terms in the solutions of (16) and (17), allowing the exact fuzzy output regulation.(5)
The proposed controller can be seen as a simple substitution of the aforementioned elements.
On the other hand, the following disadvantages can appear:(1)
If the condition is that the steady-state zero error manifoldπ(wk)≠Πwk, that is, Πij≠Π, then, it will be necessary to adjust the local regulator by means of new membership functions. Please refer to Section 3.(2)
As expected, the complexity of the controller increases according to the number of local subsystems.The following theorem provides the conditions for the existence of the exact fuzzy output regulator for a discrete-time T-S fuzzy models.Theorem 3.
The exact fuzzy output regulation with full information of systems defined as (12) is solvable if (a) there exists the same zero error steady-state manifold π(wk)=Π; (b) there exist uk=∑i=1r1hi,k(z1,k)Kixk for the fuzzy system; (c) the exosystem ωk+1=s(ωk) is Poisson stable, and the input matrices for all subsystems Bi are equal. Moreover, theExact Output Fuzzy Regulation Problem is solvable by the controller(24)uk=∑h=1r1h1,iz1,kKixk-Πwk+∑i=1r1h1,iz1,k∑j=1r2h2,jz2,kΓijwk.Proof.
From the previous analysis, the existence of mappingsπ(wk)=Πwk and γ(wk)=∑i=1r1h1,i(z1,k)∑j=1r2h2,j(z2,k)Γijwk is guaranteed when the input matrices for all subsystems Bi are equal, and the solution of r1·r2 lineal local problems is(25)ΠijSj=AiΠij+BiΓij+Pi,0=CiΠi,j-Qi,leading to Πij=Π for all i=1,…,r1 and j=1,…,r2.
On the other hand, the inclusion of condition (b) has been thoroughly discussed in [22, 23, 27, 28, 31–33], and it implies the existence of a fuzzy stabilizer.
Condition (c) ensures the nonexistence of crossed terms in the local Francis equations. Finally, condition (d) is introduced to avoid the fact that the reference signal converges to zero, which would turn the regulation problem into a simple stability problem. The rest of the proof follows directly from the previous analysis.


## 3. The Output Regulator by means of Local Regulators and Tuning of New Membership Functions

In this section, a discrete-time T-S fuzzy model is considered to solve the exact output regulation on the basis of linear local controllers. So, the main goal is to find a complete regulator based on the fuzzy summation of local regulators using adequate membership functions, such that the result given in Theorem3 can be relaxed [34]. These membership functions are not necessarily the same included in the fuzzy plant, as is described in (19). Thus, the steady-state input γ(wk) can be defined as(26)γwk=∑i=1r1μ1,iz1,k∑j=1r2μ2,jz2,kΓijwk,where μ1,i(z1,k) and μ2,j(z2,k) are new membership functions, such that the fuzzy output regulator obtained from local regulators provides the exact fuzzy output regulation. This approach requires the computation of the linear local controllers and the computation of the new membership functions. In this work, such functions are represented by the following expression:(27)μ1,ixk=11+xk-ci/ai2bi,∀i=1,…,r1,(28)μ2,jwk=11+wk-cj/aj2bj,∀j=1,…,r2.μ1,i and μ2,j are well known asgeneralized bell-shaped membership functions and the parameters ai, aj, bi, bj, ci, and cj determine the form, center, and amplitude, respectively. Therefore, from (22), the input can be defined by(29)uk=∑i=1r1h1,ixkkixk-πwk+∑i=1r1μ1,iwk∑j=1r2μ2,jwkΓijwk,because z∗,k is a function of xk and in steady-state xk=π(wk).Then, for tuning the membership functions (27) and (28), the parameters ai and aj will be optimized by means of genetic algorithms, ensuring the correct interpolation between the local linear regulators. The foregoing can be summarized in the control scheme depicted in Figure 1.Figure 1
Control scheme for fuzzy output regulation and genetic algorithms.To this end, the following algorithm should be carried out.Algorithm 4.
Main steps to take into account to solve the tuning membership functions problem are as follows [35]: (1)
Start with a randomly generated population of nl-bit chromosomes (candidate solutions to a problem). The traditional representation is binary as in Figure2.A binary string is called “chromosome.” Each position therein is called “gene” and the value in this position is named “allele.”(2)
Calculate the fitnessf(x) of each chromosome x in the population.(3)
Repeat the following steps until n offspring have been created:
(i)
Select a pair of parent chromosomes from the current population, the probability of selection being an increasing function of fitness. Selection is done “with replacement,” meaning that the same chromosome can be selected more than once to become a parent.(ii)
With probabilitypc (the “crossover probability” or “crossover rate”), crossover the pair at a randomly chosen point (chosen with uniform probability) to form two offspring. If no crossover takes place, form two offspring that are exact copies of their respective parents. (Note that here the crossover rate is defined to be the probability that two parents will cross over in a single point.)(iii)
Mutate the two offspring at each locus with probabilitypm (the mutation probability or mutation rate), and place the resulting chromosomes in the new population.(4)
Replace the current population with the new population.(5)
Go to step (2).Figure 2
Binary string commonly used in genetic algorithms.This process ends when the fitness function reaches a value less than or equal to a predefined bound, or when the maximum number of iterations is reached. Each iteration of this process is called a generation.In order to apply a genetic algorithm it requires the following five basic components:(i)
A representation of the potential solutions to the problem.(ii)
One way to create an initial population of possible solutions (typically a random process).(iii)
An evaluation function to play the role of the environment, classifying the solutions in terms of its “fitness.”(iv)
Genetic operators that alter the composition of the chromosomes that will be produced for generations.(v)
Values for the various parameters used by the genetic algorithm (population size, crossover probability, mutation probability, maximum number of generations, etc.).The aforementioned can be summarized by the flowchart depicted in Figure3.Figure 3
Flowchart of genetic algorithms.


## 4. Synchronization of T-S Discrete-Time Fuzzy Systems

The stabilization and synchronization of chaotic systems are two of the most challenging and stimulating problems due to their capabilities of describing a great variety of very interesting phenomena in physics, biology, chemistry, and engineering, to name a few. In this section, the regulation theory is used to synchronize chaotic systems described by T-S discrete-time fuzzy models. Both the drive system and the response system are modeled by the same attractor (Rössler’s equation) with the difference that response system can be influenced by an input signal. This type of synchronization is known asComplete Synchronization (CS) [36].Considering two Rössler chaotic oscillators asw˙=f(w) (drive system) and x˙=f(x,w,u) (response system), the ordinary differential equations of these systems are(30)w˙1=-w2+w3,w˙2=w1+aw2,w˙3=bw1-c-w1w3,Rössler  attractorDrive  system,x˙1=-x2+x3,x˙2=x1+ax2,x˙3=bx1-c-x1x3+u,Rössler  attractorResponse  system.According to [15], these systems can be exactly represented by means of the following two-rule T-S fuzzy models when a, b and c are constants and x1∈[c-d,c+d] with d>0. Thus(31)x˙t=∑i=12h1,ix1tA~ixt+B~iut,w˙t=∑i=12h2,iw1tS~iwt,et=∑i=12h1,ix1tC~ixt-∑i=12h2,iw1tQ~iwt,where(32)A~1=0-1-11a0b0-d,S~1=0-1-11aw0bw0-dw,A~2=0-1-11a0b0d,S~2=0-1-11aw0bw0dw,B~1=B~2=001T,C~i=Q~i=100.Besides, the membership functions for such systems are(33)h1,1x1=121+c-x1d,h1,2x1=121-c-x1d,h2,1w1=121+c-w1d,h2,2w1=121-c-w1d.Now, the continuous-time T-S fuzzy model can be converted to the following discrete counterpart by using [1](34)Gi=exp⁡AiTs=I+AiTs+Ai2Ts22!+⋯,Hi=∫0Tsexp⁡AiτBidτ=Gi-IAi-1Bi.Therefore, the discrete-time T-S fuzzy model is given by(35)xk+1=∑i=12h1,ix1,kAixk+Biuk+Piwk,wk+1=∑i=12h2,iw1,kSiwk,ek=∑i=12h1,ix1,kCixk-∑i=12h2,iw1,kQiwk,where(36)A1=1-Ts-TsTsaTs+10bTs01-dTs,S1=1-Ts-TsTsawTs+10bwTs01-dwTs,A2=1-Ts-TsTsaTs+10bTs0Tsd+1,S2=1-Ts-TsTsawTs+10bwTs0dwTs+1,B1=B2=00TsT,Ci=Qi=100.Besides, the membership functions for the T-S discrete-time fuzzy model are (37)h1,1x1,k=121+c-x1,kd,h1,2x1,k=121+c-x1,kd,h1,2w1,k=121+cw-w1,kdw,h2,2w1,k=121+cw-w1,kdw,with a=0.34, b=0.4, c=4.5, d=10,  aw=0.34, bw=0.4, cw=4.5, dw=10, and Ts as the sampling time.Figure4(b) shows the trajectory of the discrete-time version of the continuous-time T-S fuzzy Rössler model, with Ts=0.00357. It can be seen that the overall shape of the trajectory is similar to that in Figure 4(a).Figure 4
(a) Continuous T-S fuzzy model for a Rössler attractor and (b) discrete T-S fuzzy model for a Rössler attractor.
(a)
(b)Then, by using the approach derived in this work and from (20) and (17), the zero error steady-state manifold π(wk)=Π is(38)Π=100010001.Γij are(39)Γ1,1=000,Γ1,2=0020,Γ2,1=00-20,Γ2,2=000.On the other hand, the fuzzy stabilizer for this system is computed by means ofAckermann’s formula, and by locating the eigenvalues at (1)
subsystem[0.9980+0.0062i0.9980-0.0062i0.9982],(2)
subsystem[0.9980+0.0062i0.9980-0.0062i0.9982],the following gains are obtained:(40)K1=3.16290.92578.0354,K2=3.16290.9257-11.9646.Remark 5.
It is important to verify that the fuzzy feedback stabilizer is valid for the overall T-S fuzzy model, by checking that the eigenvalues of the interpolation regions are inside the unit circle also [31].Then, by setting the initial conditions asx1,k=5, x2,k=0, x3,k=6, w1,k=1, w2,k=0, and w1,k=0 and by applying the controller (23), the results depicted in Figures 5 and 6 are obtained.Figure 5
Control signal for the Complete Synchronization of two discrete-time T-S fuzzy Rössler systems.Figure 6
Tracking error signal for the Complete Synchronization of two discrete-time T-S fuzzy Rössler systems.The tracking for the drive states and response states is drawn in Figure7.Figure 7
States of drive and response tracking for a Complete Synchronization.
(a)
(b)
(c)


### 4.1. Generalized Synchronization

The discrete-time fuzzy synchronization problem is solvable when the conditions of Theorem3 are fulfilled. However, when two chaotic systems are different, fulfilling these conditions is not so common. This is because the local regulators have different zero error steady-state manifolds, in general [20]. From the regulation point of view, the problem can be seen as finding, if it is possible, a transformation π:wk→xk regarding mapping the trajectories of the drive system into the ones of the response systems; that is, xk=π(wk); this is known as Generalized Synchronization [37, 38] and satisfies (41)limk→∞⁡xkx0-wkw0=0,with x(0) and w(0) as initial conditions.Now, consider the chaotic drive system asw˙=f(w) and x˙=f(x,w,u) as the response system.The equations for the aforementioned systems are described as follows:(42)w˙1=αw2-w1,w˙2=rw1-w2-w1w3,w˙3=w1w2-βw3,Lorenz  attractorDrive  systemx˙1=ax2-x1,x˙2=c-ax1-x1x3+cx2+u,x˙3=x1x2-bx3,Chen's  attractorResponse  systemAs before these systems can be exactly represented by means of the following two-rule continuous-time T-S fuzzy models whenx1∈[Xmin⁡,Xmax⁡] and w1∈[M1,M2]. Therefore,(43)x˙t=∑i=12h1,ix1tA~ixt+B~iut,w˙t=∑i=12h2,iw1tS~iwt,et=∑i=12h1,ix1tC~ixt-∑i=12h2,iw1tQ~iwt,where(44)A~1=-aa0c-ac-Xmin⁡0Xmin⁡-b,S~1=-αα0r-1-M10M1-β,A~2=-aa0c-ac-Xmax⁡0Xmax⁡-b,S~2=-αα0r-1-M20M2-β,B~1=B~2=010T,C~i=Q~i=100.Then, by using (34), the discrete counterpart is obtained:(45)xk+1=∑i=12h1,ix1,kAixk+Biuk+Piwk,wk+1=∑i=12h2,iw1,kSiwk,ek=∑i=12h1,ix1,kCixk-∑i=12h2,iw1,kQiwk,where (46)A1=1-aTsaTs0-a-cTscTs+1-Xmin⁡Ts0TsXmin⁡1-bTs,S1=1-awTsawTs0cwTs1-Ts-M1Ts0M1Ts1-bwTs,A2=1-aTsaTs0-a-cTscTs+1-Xmax⁡Ts0TsXmax⁡1+bTs,S2=1-awTsawTs0cwTs1-Ts-M2Ts0M2Ts1-bwTs,B1=B2=0Ts0T,Ci=Qi=100,with a=35, b=3, c=28, aw=10, bw=8/3, and cw=28 and Ts as the sampling time. Notice that, in this case, the Generalized Synchronization problem consists of the tracking of w1 by x1. This can be inferred by the form of Ci and Qi. The membership functions for this system are defined as follows:Plant Membership Functions
(47)
h
1,1
x
1
,
k
=
-
x
1
,
k
+
X
m
a
x
X
m
a
x
-
X
min
⁡
,
h
1,2
x
1
,
k
=
x
1
,
k
-
X
min
⁡
X
m
a
x
-
X
min
⁡
.Exosystem Membership Functions
(48)
h
2,1
w
1
,
k
=
-
w
1
,
k
+
M
1
M
2
-
M
1
,
h
2,2
w
1
,
k
=
w
1
,
k
-
M
2
M
2
-
M
1
.They are depicted in Figure 8.Figure 8
Membership functions for Generalized Synchronization.Figures9(a) and 9(b) show the behavior of the two discrete-time T-S fuzzy models, with x0=[111]T and w0=[100]T. Then, by using the approach derived in this work and from (20) and (17), the zero error steady-state manifold for each subsystem is (49)Π1,1=1000.71430.285703.0041-0.01431.2861,Π1,2=1000.71430.285703.0041-0.0143-1.2861,Π2,1=1000.71430.28570-3.00410.0143-1.2861,Π2,2=1000.71430.28570-3.00410.01431.2861.Figure 9
(a) Discrete-time T-S fuzzy model for a Chen attractor and (b) discrete T-S fuzzy model for a Lorenz attractor.
(a)
(b)Γ
i
j are(50)Γ1,1=-102.2648-0.7142-30.0121,Γ1,2=-102.2648-0.714230.0121,Γ2,1=-102.2648-0.7142-30.0121,Γ2,2=-102.2648-0.714230.0121.On the other hand, the fuzzy stabilizer for this system is computed by means ofAckermann's formula, and by locating the eigenvalues at(1)
subsystem0.95960.99550.82,(2)
subsystem0.95960.99550.82,the following gains are obtained:(51)K1=-14.3228-214.9000-19.6566,K2=-14.3228-214.900019.6566.Notice that the stability is ensured in the fuzzy interpolation region by these gains, for all t≥0.As can be readily observed, conditions of Theorem3 are not fulfilled because the zero error steady-state manifold is not the same for the local regulators. However, the tracking error can be bounded by using the controller defined in (22). Figures 10–12 depicted the behavior of the controller (22) with x0=[111]T and w(0)=[100]T as the initial conditions.Figure 10
Control signal for the Generalized Synchronization of Chen-Lorenz discrete-time T-S fuzzy systems.Figure 11
Tracking errore1,k signal for the Generalized Synchronization of Chen-Lorenz discrete-time T-S fuzzy systems.Figure 12
States of drive and response system for a Generalized Synchronization.
(a)
(b)
(c)


### 4.2. Generalized Synchronization by Using Genetic Algorithms

The main objective of integrating genetic algorithms is to obtain values of the parameters (a, b, and c) of the new membership functions of the regulator and reducing, in this way, the tracking error. To this end, it is important to consider the following:(i)
Consider the representation of the potential solutions of the problem; the common representation is binary.(ii)
An initial population of possible solutions, in a random process, is selected (population size affects the efficiency and performance of GA). For this case, a population size equal to 20 chromosomes is chosen.(iii)
A fitness function, which indicates how good or bad a certain solution is, is defined. In this case the fitness function is the mean square tracking error, given by the following expression:(52)es=12x1,k-w1,k2.For the representation of possible solutions it is necessary to know the variables to optimize; for this case it hasai and aj parameters for i=1,2 and j=1,2 defined in (27) and (28); besides, bi, bj, ci, and cj will be fixed with a constant value; that is, b1=3, c1=30 for (27) considering μ1,2=1-μ1,1; b1=3, c1=20 for (28) also considering μ2,2=1-μ2,1. The abovementioned are represented in Figure 13.Figure 13
Initial membership functions for the proposed fuzzy regulator tuned by GAs.
(a)
(b)It is also important to know the interval value in which the variablea will be operating; for this case, such interval is x1,k∈[-3030]. Thus, the size of each variable, in bits, and the chromosome length can be computed by using the following expression [39]:(53)Svar=int↑log2⁡limupper-limlower10accurate,where Svar is the size of each variable in bits, the term (int⁡↑) is the decimal integer value, and limupper=30 and limlower=-30 with an accuracy in 3. Therefore, the variable size is equal to 16 bits, but there are two variables, a1 and a2; thus the chromosomes length is 32 bits. Therefore, the chromosomes are represented in Figure 14.Figure 14
Individual binary representation.The initial population is random binary. To assess each individual in the objective function is needed to decode, in this case, a real number. To this end, the following expression is used:(54)Var=Decvallimupper⁡-limlower2lengthbit-1+limlower.Var is a decimal value (phenotype) for each binary chromosome (genotype), Decval is the decimal value for each binary string, and length(bit) is equal to 16 bits.In order to compute the individual aptitude is necessary to introduce the value of each parameter in the bell-shaped membership functionsμ1 and μ2. The tracking error function (52) disregards the nonoptimal solution and allows the optimal performance by considering that x1,k tracks w1,k. In Table 1 the individual aptitude is depicted.Table 1
Individual aptitude.
Chromosome
Binary string
Decoded integer
Fitness function
1
01011010001000100010111010101101
f(−7.691493, −16.518837)
0.971109
2
11001011101111111110101000000111
f(15.386557, 21.537530)
6.840302
3
10100011000010001110111101010111
f(7.116228, 22.616648)
9.404465
4
01101110111100100011110001011000
f(−3.463890, −13.742489)
0.096779
5
01111110110010001000111010000000
f(−0.247166, 2.945754)
0.017387
6
01111110110010001000111010000000
f(−2.282414, −24.456703)
12.438373
7
01110100110000110000011110011001
f(−13.451286, 16.387121)
0.761323
8
00111101110001111101000010101100
f(4.320043, −25.338247)
13.686271
9
10010101010001000000001101000010
f(19.709377, 4.057404)
0.209690
10
10000010001011110001000101110111
f(0.443946, −22.452399)
9.159322
11
00111100101000111011110011100001
f(−13.682979, 12.366613)
0.000746
12
00001100011111011100101000101100
f(−23.463279, 15.066789)
0.000895
13
10010110000010001001010001111101
f(4.475563, 4.162142)
0.017576
14
10110011010110110100100101111010
f(10.432136, −11.074846)
0.001021
15
01010011101010101111101011001101
f(−9.005478, 24.944686)
13.030471
16
00100011011001100011000000001011
f(−18.809583, −16.241123)
0.399280
17
10010001100100001111110100111101
f(3.567834, 25.439811)
13.825474
18
10001011001010110100111011110001
f(2.268925, −9.964782)
0.003840
19
01100000011011101100000001001011
f(−6.412421, 13.060105)
0.039420
20
00001011110001010001010011000101
f(−23.609277, −21.781125)
5.130460Each fitness function value is converted into a set point or fitness value. To this end, the following expression is used:(55)fitnessval=F⋆-Fmin⁡,where F(⋆) is the fitness function value and Fmin⁡ is the minimum fitness function value.The two chromosomes with better fitness (elite chromosomes) can live and produce offspring in the next generation. There are many methods to select a pair of chromosomes; the most popular one is namedproportional selection method:(a)
Calculate the total fitness value(56)fitnessT=∑i=120fitnessi.(b)
Compute the probabilityPi for each chromosome(57)Pi=fitnessifitnessT.(c)
Calculate the cumulative probability for each chromosome(58)Qk=∑k=0iPk.The total fitness valuefitnessT is equal to 86.03420. The probabilities Pi and Qk are shown in Table 2.(d)
Generate a random numberr in the range [0,1].(e)
IfQi-1<r≤Qi, then select the chromosome to be the one of the parents.(f)
Repeat (d) and (e) to obtain the other parent.Table 2
Individual aptitude.
Chromosome
P
i
Q
k
1
0.011
0.011
2
0.080
0.091
3
0.109
0.200
4
0.001
0.201
5
0.000
0.201
6
0.145
0.346
7
0.009
0.355
8
0.159
0.514
9
0.002
0.516
10
0.106
0.623
11
0.000
0.623
12
0.000
0.623
13
0.000
0.623
14
0.000
0.623
15
0.151
0.775
16
0.005
0.779
17
0.161
0.940
18
0.000
0.940
19
0.000
0.940
20
0.060
1.000Apply crossover operation on the selected pair, if they have been chosen for crossover (based on probability of crossoverpc=1.0). The most applied crossover operation is single point crossover. Based on the probability of bit mutation pm=0.01, flip the correspondent bit if selected for mutation. At this point, the process of producing a pair of offspring from two selected parents is finished.The elite chromosomes of the previous population are not subject to mutation.In the following simulations the statesx1,k, control signal uk, and tracking error e1,k depicted in Figures 16–18, respectively, are obtained after applying (29) and by replacing the new membership functions adjusted by GAs. The final membership functions after tuning the form parameter by means of GA are depicted in Figure 15. It can be readily observed that these new membership functions are different from the original ones (see Figure 13). Notice that, in this example, the tracking error (Figure 18) is less than the error obtained by using the approach discussed in Section 4. This is due to the new membership functions which allow reducing considerably the tracking error. This suggests that the approach presented in this work may be improved by tuning the parameters of center and amplitude in the new membership functions; however this study is not completed yet.Figure 15
Final membership functions for the proposed fuzzy regulator tuned by GAs.
(a)
(b)Figure 16
Statesx1,k and w1,k for the Generalized Synchronization of Chen-Lorenz discrete-time T-S fuzzy systems.Figure 17
Control signal for the Generalized Synchronization of Chen-Lorenz discrete-time T-S fuzzy systems.Figure 18
Tracking errore1,k signal for the Generalized Synchronization of Chen-Lorenz discrete-time T-S fuzzy systems.


### 4.3. Pseudofuzzy Generalized Synchronization by Using GAs

In this section the Generalized Synchronization by using GAs and linear local regulator design will be addressed. However, the crossed terms within the local regulators are arbitrarily removed from the design process. Thus, the steady-state inputγ(wk) can be defined as(59)γwk=∑i=1r1μ1,iz1,kΓiiwk,where μ1,i(z1,k) are new membership functions, such that the fuzzy output regulator obtained from local regulators provides the exact fuzzy output regulation. This approach requires the computation of the linear local controllers and the computation of the new membership functions; such functions are represented by (27) and (28). Therefore the final control system is defined by(60)uk=∑i=1r1h1,ixkkixk-πwk+∑i=1r1μ1,iwkΓiiwk,because z∗,k is a function of xk and in steady-state xk=π(wk).For this case the initial membership functions proposed are shown in Figure19. As it can be seen, these membership functions are not properly fuzzy because they do not fulfill the convex sum; that is, ∑i=1r1μi(wk)≠0. Even so, these functions will be adjusted by GAs in order to ensure the output regulation by using input (60). The following simulation provides a better behavior in the Lorenz-Chen discrete-time fuzzy system by using different membership function and tuning by GAs. Such results are depicted in Figures 20–22.Figure 19
Initial membership functions for the pseudofuzzy regulator tuned by GAs.Figure 20
Statesx1,k and w1,k for the Generalized Synchronization of Chen-Lorenz discrete-time T-S fuzzy systems.Figure 21
Control signal for the Generalized Synchronization of Chen-Lorenz discrete-time T-S fuzzy systems.Figure 22
Tracking errore1,k for the Generalized Synchronization of Chen-Lorenz discrete-time T-S fuzzy systems.Notice that the final membership function depicted in Figure23 is different from the initial proposed in Figure 19 and also notice that the sum of these new interpolation functions is not equal to one. For that reason this approach is called Pseudofuzzy Generalized Synchronization.Figure 23
Final interpolation functions for the pseudofuzzy regulator tuned by GAs.


## 4.1. Generalized Synchronization

The discrete-time fuzzy synchronization problem is solvable when the conditions of Theorem3 are fulfilled. However, when two chaotic systems are different, fulfilling these conditions is not so common. This is because the local regulators have different zero error steady-state manifolds, in general [20]. From the regulation point of view, the problem can be seen as finding, if it is possible, a transformation π:wk→xk regarding mapping the trajectories of the drive system into the ones of the response systems; that is, xk=π(wk); this is known as Generalized Synchronization [37, 38] and satisfies (41)limk→∞⁡xkx0-wkw0=0,with x(0) and w(0) as initial conditions.Now, consider the chaotic drive system asw˙=f(w) and x˙=f(x,w,u) as the response system.The equations for the aforementioned systems are described as follows:(42)w˙1=αw2-w1,w˙2=rw1-w2-w1w3,w˙3=w1w2-βw3,Lorenz  attractorDrive  systemx˙1=ax2-x1,x˙2=c-ax1-x1x3+cx2+u,x˙3=x1x2-bx3,Chen's  attractorResponse  systemAs before these systems can be exactly represented by means of the following two-rule continuous-time T-S fuzzy models whenx1∈[Xmin⁡,Xmax⁡] and w1∈[M1,M2]. Therefore,(43)x˙t=∑i=12h1,ix1tA~ixt+B~iut,w˙t=∑i=12h2,iw1tS~iwt,et=∑i=12h1,ix1tC~ixt-∑i=12h2,iw1tQ~iwt,where(44)A~1=-aa0c-ac-Xmin⁡0Xmin⁡-b,S~1=-αα0r-1-M10M1-β,A~2=-aa0c-ac-Xmax⁡0Xmax⁡-b,S~2=-αα0r-1-M20M2-β,B~1=B~2=010T,C~i=Q~i=100.Then, by using (34), the discrete counterpart is obtained:(45)xk+1=∑i=12h1,ix1,kAixk+Biuk+Piwk,wk+1=∑i=12h2,iw1,kSiwk,ek=∑i=12h1,ix1,kCixk-∑i=12h2,iw1,kQiwk,where (46)A1=1-aTsaTs0-a-cTscTs+1-Xmin⁡Ts0TsXmin⁡1-bTs,S1=1-awTsawTs0cwTs1-Ts-M1Ts0M1Ts1-bwTs,A2=1-aTsaTs0-a-cTscTs+1-Xmax⁡Ts0TsXmax⁡1+bTs,S2=1-awTsawTs0cwTs1-Ts-M2Ts0M2Ts1-bwTs,B1=B2=0Ts0T,Ci=Qi=100,with a=35, b=3, c=28, aw=10, bw=8/3, and cw=28 and Ts as the sampling time. Notice that, in this case, the Generalized Synchronization problem consists of the tracking of w1 by x1. This can be inferred by the form of Ci and Qi. The membership functions for this system are defined as follows:Plant Membership Functions
(47)
h
1,1
x
1
,
k
=
-
x
1
,
k
+
X
m
a
x
X
m
a
x
-
X
min
⁡
,
h
1,2
x
1
,
k
=
x
1
,
k
-
X
min
⁡
X
m
a
x
-
X
min
⁡
.Exosystem Membership Functions
(48)
h
2,1
w
1
,
k
=
-
w
1
,
k
+
M
1
M
2
-
M
1
,
h
2,2
w
1
,
k
=
w
1
,
k
-
M
2
M
2
-
M
1
.They are depicted in Figure 8.Figure 8
Membership functions for Generalized Synchronization.Figures9(a) and 9(b) show the behavior of the two discrete-time T-S fuzzy models, with x0=[111]T and w0=[100]T. Then, by using the approach derived in this work and from (20) and (17), the zero error steady-state manifold for each subsystem is (49)Π1,1=1000.71430.285703.0041-0.01431.2861,Π1,2=1000.71430.285703.0041-0.0143-1.2861,Π2,1=1000.71430.28570-3.00410.0143-1.2861,Π2,2=1000.71430.28570-3.00410.01431.2861.Figure 9
(a) Discrete-time T-S fuzzy model for a Chen attractor and (b) discrete T-S fuzzy model for a Lorenz attractor.
(a)
(b)Γ
i
j are(50)Γ1,1=-102.2648-0.7142-30.0121,Γ1,2=-102.2648-0.714230.0121,Γ2,1=-102.2648-0.7142-30.0121,Γ2,2=-102.2648-0.714230.0121.On the other hand, the fuzzy stabilizer for this system is computed by means ofAckermann's formula, and by locating the eigenvalues at(1)
subsystem0.95960.99550.82,(2)
subsystem0.95960.99550.82,the following gains are obtained:(51)K1=-14.3228-214.9000-19.6566,K2=-14.3228-214.900019.6566.Notice that the stability is ensured in the fuzzy interpolation region by these gains, for all t≥0.As can be readily observed, conditions of Theorem3 are not fulfilled because the zero error steady-state manifold is not the same for the local regulators. However, the tracking error can be bounded by using the controller defined in (22). Figures 10–12 depicted the behavior of the controller (22) with x0=[111]T and w(0)=[100]T as the initial conditions.Figure 10
Control signal for the Generalized Synchronization of Chen-Lorenz discrete-time T-S fuzzy systems.Figure 11
Tracking errore1,k signal for the Generalized Synchronization of Chen-Lorenz discrete-time T-S fuzzy systems.Figure 12
States of drive and response system for a Generalized Synchronization.
(a)
(b)
(c)


## 4.2. Generalized Synchronization by Using Genetic Algorithms

The main objective of integrating genetic algorithms is to obtain values of the parameters (a, b, and c) of the new membership functions of the regulator and reducing, in this way, the tracking error. To this end, it is important to consider the following:(i)
Consider the representation of the potential solutions of the problem; the common representation is binary.(ii)
An initial population of possible solutions, in a random process, is selected (population size affects the efficiency and performance of GA). For this case, a population size equal to 20 chromosomes is chosen.(iii)
A fitness function, which indicates how good or bad a certain solution is, is defined. In this case the fitness function is the mean square tracking error, given by the following expression:(52)es=12x1,k-w1,k2.For the representation of possible solutions it is necessary to know the variables to optimize; for this case it hasai and aj parameters for i=1,2 and j=1,2 defined in (27) and (28); besides, bi, bj, ci, and cj will be fixed with a constant value; that is, b1=3, c1=30 for (27) considering μ1,2=1-μ1,1; b1=3, c1=20 for (28) also considering μ2,2=1-μ2,1. The abovementioned are represented in Figure 13.Figure 13
Initial membership functions for the proposed fuzzy regulator tuned by GAs.
(a)
(b)It is also important to know the interval value in which the variablea will be operating; for this case, such interval is x1,k∈[-3030]. Thus, the size of each variable, in bits, and the chromosome length can be computed by using the following expression [39]:(53)Svar=int↑log2⁡limupper-limlower10accurate,where Svar is the size of each variable in bits, the term (int⁡↑) is the decimal integer value, and limupper=30 and limlower=-30 with an accuracy in 3. Therefore, the variable size is equal to 16 bits, but there are two variables, a1 and a2; thus the chromosomes length is 32 bits. Therefore, the chromosomes are represented in Figure 14.Figure 14
Individual binary representation.The initial population is random binary. To assess each individual in the objective function is needed to decode, in this case, a real number. To this end, the following expression is used:(54)Var=Decvallimupper⁡-limlower2lengthbit-1+limlower.Var is a decimal value (phenotype) for each binary chromosome (genotype), Decval is the decimal value for each binary string, and length(bit) is equal to 16 bits.In order to compute the individual aptitude is necessary to introduce the value of each parameter in the bell-shaped membership functionsμ1 and μ2. The tracking error function (52) disregards the nonoptimal solution and allows the optimal performance by considering that x1,k tracks w1,k. In Table 1 the individual aptitude is depicted.Table 1
Individual aptitude.
Chromosome
Binary string
Decoded integer
Fitness function
1
01011010001000100010111010101101
f(−7.691493, −16.518837)
0.971109
2
11001011101111111110101000000111
f(15.386557, 21.537530)
6.840302
3
10100011000010001110111101010111
f(7.116228, 22.616648)
9.404465
4
01101110111100100011110001011000
f(−3.463890, −13.742489)
0.096779
5
01111110110010001000111010000000
f(−0.247166, 2.945754)
0.017387
6
01111110110010001000111010000000
f(−2.282414, −24.456703)
12.438373
7
01110100110000110000011110011001
f(−13.451286, 16.387121)
0.761323
8
00111101110001111101000010101100
f(4.320043, −25.338247)
13.686271
9
10010101010001000000001101000010
f(19.709377, 4.057404)
0.209690
10
10000010001011110001000101110111
f(0.443946, −22.452399)
9.159322
11
00111100101000111011110011100001
f(−13.682979, 12.366613)
0.000746
12
00001100011111011100101000101100
f(−23.463279, 15.066789)
0.000895
13
10010110000010001001010001111101
f(4.475563, 4.162142)
0.017576
14
10110011010110110100100101111010
f(10.432136, −11.074846)
0.001021
15
01010011101010101111101011001101
f(−9.005478, 24.944686)
13.030471
16
00100011011001100011000000001011
f(−18.809583, −16.241123)
0.399280
17
10010001100100001111110100111101
f(3.567834, 25.439811)
13.825474
18
10001011001010110100111011110001
f(2.268925, −9.964782)
0.003840
19
01100000011011101100000001001011
f(−6.412421, 13.060105)
0.039420
20
00001011110001010001010011000101
f(−23.609277, −21.781125)
5.130460Each fitness function value is converted into a set point or fitness value. To this end, the following expression is used:(55)fitnessval=F⋆-Fmin⁡,where F(⋆) is the fitness function value and Fmin⁡ is the minimum fitness function value.The two chromosomes with better fitness (elite chromosomes) can live and produce offspring in the next generation. There are many methods to select a pair of chromosomes; the most popular one is namedproportional selection method:(a)
Calculate the total fitness value(56)fitnessT=∑i=120fitnessi.(b)
Compute the probabilityPi for each chromosome(57)Pi=fitnessifitnessT.(c)
Calculate the cumulative probability for each chromosome(58)Qk=∑k=0iPk.The total fitness valuefitnessT is equal to 86.03420. The probabilities Pi and Qk are shown in Table 2.(d)
Generate a random numberr in the range [0,1].(e)
IfQi-1<r≤Qi, then select the chromosome to be the one of the parents.(f)
Repeat (d) and (e) to obtain the other parent.Table 2
Individual aptitude.
Chromosome
P
i
Q
k
1
0.011
0.011
2
0.080
0.091
3
0.109
0.200
4
0.001
0.201
5
0.000
0.201
6
0.145
0.346
7
0.009
0.355
8
0.159
0.514
9
0.002
0.516
10
0.106
0.623
11
0.000
0.623
12
0.000
0.623
13
0.000
0.623
14
0.000
0.623
15
0.151
0.775
16
0.005
0.779
17
0.161
0.940
18
0.000
0.940
19
0.000
0.940
20
0.060
1.000Apply crossover operation on the selected pair, if they have been chosen for crossover (based on probability of crossoverpc=1.0). The most applied crossover operation is single point crossover. Based on the probability of bit mutation pm=0.01, flip the correspondent bit if selected for mutation. At this point, the process of producing a pair of offspring from two selected parents is finished.The elite chromosomes of the previous population are not subject to mutation.In the following simulations the statesx1,k, control signal uk, and tracking error e1,k depicted in Figures 16–18, respectively, are obtained after applying (29) and by replacing the new membership functions adjusted by GAs. The final membership functions after tuning the form parameter by means of GA are depicted in Figure 15. It can be readily observed that these new membership functions are different from the original ones (see Figure 13). Notice that, in this example, the tracking error (Figure 18) is less than the error obtained by using the approach discussed in Section 4. This is due to the new membership functions which allow reducing considerably the tracking error. This suggests that the approach presented in this work may be improved by tuning the parameters of center and amplitude in the new membership functions; however this study is not completed yet.Figure 15
Final membership functions for the proposed fuzzy regulator tuned by GAs.
(a)
(b)Figure 16
Statesx1,k and w1,k for the Generalized Synchronization of Chen-Lorenz discrete-time T-S fuzzy systems.Figure 17
Control signal for the Generalized Synchronization of Chen-Lorenz discrete-time T-S fuzzy systems.Figure 18
Tracking errore1,k signal for the Generalized Synchronization of Chen-Lorenz discrete-time T-S fuzzy systems.


## 4.3. Pseudofuzzy Generalized Synchronization by Using GAs

In this section the Generalized Synchronization by using GAs and linear local regulator design will be addressed. However, the crossed terms within the local regulators are arbitrarily removed from the design process. Thus, the steady-state inputγ(wk) can be defined as(59)γwk=∑i=1r1μ1,iz1,kΓiiwk,where μ1,i(z1,k) are new membership functions, such that the fuzzy output regulator obtained from local regulators provides the exact fuzzy output regulation. This approach requires the computation of the linear local controllers and the computation of the new membership functions; such functions are represented by (27) and (28). Therefore the final control system is defined by(60)uk=∑i=1r1h1,ixkkixk-πwk+∑i=1r1μ1,iwkΓiiwk,because z∗,k is a function of xk and in steady-state xk=π(wk).For this case the initial membership functions proposed are shown in Figure19. As it can be seen, these membership functions are not properly fuzzy because they do not fulfill the convex sum; that is, ∑i=1r1μi(wk)≠0. Even so, these functions will be adjusted by GAs in order to ensure the output regulation by using input (60). The following simulation provides a better behavior in the Lorenz-Chen discrete-time fuzzy system by using different membership function and tuning by GAs. Such results are depicted in Figures 20–22.Figure 19
Initial membership functions for the pseudofuzzy regulator tuned by GAs.Figure 20
Statesx1,k and w1,k for the Generalized Synchronization of Chen-Lorenz discrete-time T-S fuzzy systems.Figure 21
Control signal for the Generalized Synchronization of Chen-Lorenz discrete-time T-S fuzzy systems.Figure 22
Tracking errore1,k for the Generalized Synchronization of Chen-Lorenz discrete-time T-S fuzzy systems.Notice that the final membership function depicted in Figure23 is different from the initial proposed in Figure 19 and also notice that the sum of these new interpolation functions is not equal to one. For that reason this approach is called Pseudofuzzy Generalized Synchronization.Figure 23
Final interpolation functions for the pseudofuzzy regulator tuned by GAs.


## 5. Conclusions

A fuzzy output regulator for discrete-time systems, based on the combination of linear regulators combined by different membership functions, has been presented. Synchronization of discrete-time chaos attractors can be possible; by means of fuzzy output regulation, sufficient conditions for the controller are given. However, when the conditions can not be fulfilled, new membership functions in the regulator are included; these ones are optimized by genetic algorithms. The main advantage is that membership functions, which allow the proper combination of the local regulators, can be easily obtained by means of GAs. As a consequence, the presented result allows a very precise synchronization of chaotic systems described by T-S discrete-time fuzzy models on the basis of local regulators.Complete and Generalized Synchronization are used to illustrate the applicability of the proposed approach. Besides, the method proposed in this work avoids the disadvantage of constructing an exact fuzzy regulator based on overall T-S fuzzy system which may result to be very large. Instead, the given approach offers a simple way to design the complete regulator based on local regulators but with membership functions optimized by soft computing.



---

*Source: 198371-2015-12-09.xml*