Partial synchronization of different chaotic oscillators using robust PID feedback

Transcription

1 Chaos, Solitons and Fractals 33 (27) Partial synchronization of different chaotic oscillators using robust PID feedback Ricardo Aguilar-López a, *, Rafael Martinez-Guerra b a Departamento de Energía, Universidad Autónoma Metropolitana Azcapotzalco, San Pablo 18, Reynosa-Tamaulipas, Azcapotzalco, 22 México, D.F., Mexico b Departamento de Control Automático, CINVESTAV IPN, Apartado Postal 14-74, México, D.F. C.P. 736, Mexico Accepted 2 December 25 Abstract This work deals with the partial synchronization problem of two different chaotic oscillators considering model uncertainties in the slave system via control approach. The slave system is forced to follow the master signal via a linearizing controller based on model uncertainty reconstructor which leads to proportional integral derivative (PID) control structure. This reconstructor is related with a proportional derivative (PD) reduced-order observer, it would be considered as a sub-slave system for the original slave of the synchronization procedure. The asymptotic performance of the synchronization methodology is proven via the dynamic of the synchronization error. Numerical experiment illustrates the closed-loop behavior of the proposed methodology. Ó 26 Elsevier Ltd. All rights reserved. 1. Introduction As is well known, the study of the synchronization problem for chaotic oscillators has been very important from the non-linear science point of view, in particular, the applications to biology, medicine, cryptography, secure data transmission and so on. In general, the synchronization research has been focused on two areas: the first one, related with the employment of state observers, where the main applications lie on the synchronization of non-linear oscillators with the same model structure and order, but different initial conditions and/or parameters [1 5]; on the other hand, the use of control laws allows to achieve the synchronization between non-linear oscillators, with different structure and order, where the variable states of the slave system are forced to follow the trajectories of the master system. This approach can be seen as a tracking problem [6 1]; some authors design the controller from the dynamic of the synchronization error, because this approach allows transforming the tracking problem to a regulation problem with the origin (zero) as the corresponding set point [11]. Besides, several control approaches have used neural-network, fuzzy, adaptive, and sliding and other techniques [8]. Other traditional control methods [12] consider introducing an additive feedback controller, to force to the system to * Corresponding author. Tel.: ; fax: addresses: raguilar@correo.azc.uam.mx (R. Aguilar-López), rguerra@ctrl.cinvestav.mx (R. Martinez-Guerra) /$ - see front matter Ó 26 Elsevier Ltd. All rights reserved. doi:1.116/j.chaos

2 R. Aguilar-López, R. Martinez-Guerra / Chaos, Solitons and Fractals 33 (27) reach the desired reference (set point), i.e. kx(t) x sp (t)k!ast!1. The above mentioned methodologies are based on the cancellation of the non-linear terms of the chaotic systems in order to impose a desired behavior. Under the above philosophy, the non-linear geometric-differential control techniques have been successfully employed [13,14]. They correspond to systems that can be fully or partially linearized by a change of coordinates and/or state feedback following differential-geometric concepts [6]. Such class of non-linear systems can be linearized by a state feedback control, which cancels all the non-linearities assuming perfect knowledge of the mathematical model, producing global asymptotic stability [7]. A drawback of exact linearization techniques and the other model based controllers is that they rely on exact cancellation of non-linearities. In practice, exact knowledge of system dynamics is not possible. A more realistic situation is to know some nominal functions of the corresponding non-linearities, which are employed in the control design. However, the use of nominal model non-linearities can lead to performance degradation and even closed-loop instability. In fact, when the systems possess strong non-linearities, the standard linearizing, generic model, and active controllers, cannot cancel completely such non-linearities and instabilities can be induced. The worst case is when the knowledge of the non-linearities is very poor or null, such that, conventional linearizing techniques are inadequate. To avoid the above problems, the geometric approach for the design of non-linear controllers based on uncertainty observers has been employed; and these kinds of techniques show satisfactory capabilities for a wide range of systems [15,16]. The use of proportional observers coupled with linearizing controllers have been very successful, but the proportional observers have several problems, for example, they are very sensitive for noisy measurements, robustness issues are not completely saved. For these reasons a more sophisticated observers have been designed in order to generate more adequate open-loop and closed-loop performances. PI observers, sliding-mode, numeric, etc. have been developed [17,18]. In particular, in this paper is proposed an uncertainty proportional derivative reduced-order observer which is able to accelerate the estimation rate and then it is coupled to a linearizing controller to produce a robust proportional integral derivative (PID) structure. This methodology is applied to the synchronization of one signal of the slave system, the non-linear oscillators considered here have different order and structure and the examples proposed are de Duffing Chen systems and the Van der Pol Lorenz systems. 2. Problem description A generic representation of master system can be described by Eqs. (1a) and (1b) as follows: _X m ¼ f m ðx m ; S m Þ¼I m ðx m Þþ m ðx m ; S m Þ y m ¼ h m ðx m Þ¼CX m þ - m ð1aþ ð1bþ Here, X m 2 R n is the vector of states; S m 2 R is the message signal; I m ðþ : R nþ1! R n is a non-linear vector field; m ðþ : R nþ1! R n is a linear vector of its arguments; - m 2 R m is an additive bounded measurement noise and y 2 R m is the vector of transmitted states. Now, considering the following representation for the slave system: _X s ¼ f s ðx s ; UÞ ¼I s ðx s Þþ s ðx s ÞU y s ¼ h s ðx s Þ¼CX s þ - s ð2aþ ð2bþ where the vectors and the maps considered are defined in a similar way as mentioned above and U 2 R is the corresponding control command. Defining the m X s (synchronization error), the corresponding dynamic equation considering Eqs. (1) and (2) is as _ þ m ðx m ; S m Þ s ðx s ÞU ð3þ Now, let us consider the following hypothesis: H1. For the realized control input vector U(X(t)), ðkuðx ðtþþk 6 UÞ, the nominal closed-loop non-linear system ((1a) and (1b)) is quadratically stable; therefore there exists a Lyapunov function V P that satisfies ov o@ ½DIð@Þ þ mðx m ; S m Þ s ðx s ÞUŠ 6 b 1 k@k 2 ; Such that the system is a minimum phase. ov o@ 6 b 2k@k; b 1 ; b 2 >

3 574 R. Aguilar-López, R. Martinez-Guerra / Chaos, Solitons and Fractals 33 (27) H2. - is a vector function representing an external (maybe unknown), bounded perturbation k-k 2 D ¼ }<1; < D ¼ DT The normalized matrix D is introduced to ensure the possibility to deal with components of different physical structure and is assumed as given a priori; } represents the power of the corresponding perturbation. H3. The linear vector field () () is bounded, i.e. for n, k ()k 6 + < Synchronization methodology design Let us assume that there exists a coordinate transformation ðx; tþ¼wð@þ such that the system (3) can be locally transformed into canonical form given by _x i ¼ x iþ1 ; i ¼ 1; 2;...; r 1 _x r ¼ f ðx; tþþbðx; tþu ð4þ _t ¼ #ðx; tþ y ¼ x 1 Here, the following definitions were used: Bðx; tþ ¼BðxÞþDB, where f(x, t) and DB are considered model uncertainties related to the non-linear system, BðxÞ is the nominal value of the control input coefficient, u is the system input, x are the system states, y is the slave measured signal, r is the relative degree of the system. Now, it is possible to define a convenient change of variables gðx; uþ ¼f ðx; tþþdbu By substitution of Eq. (4 ) into Eq. (4), a new system is obtained _x i ¼ x iþ1 ; i ¼ 1; 2;...; r 1 _x r ¼ gðx; uþþbðxþu _t ¼ #ðx; tþ ð4 Þ ð5þ y ¼ x 1 In order to control this system, let us define the following nominal input output linearizing feedback control. u ¼ B 1 ðxþ½s g x r gðx; uþš with BðxÞ 6¼ ð6þ The controller defined by Eq. (6) guarantees exponential stability of non-linear systems with no uncertainties and perfect measurements [14], i.e., DB = and g(x,u) known, under the assumption that the named inner dynamics _t ¼ #ðx; tþ i.e. the state equations which are not controlled present a stable dynamic behavior (minimum phase system), such that in accordance with the control theory point of view the system is stabilizable. However, since the uncertainty term, g(x,u), is unknown and, moreover, is the function of the states, x and the control input, u this ideal control law is not causal and therefore is not realizable. Nonetheless, there is another way to develop an input output linearizing controller that is robust against uncertainties. The procedure shown below defines a method to estimate the uncertainty term, g(x, u). This approach is based on observer theory, where the uncertainty is only a function of the estimation error. Let us define the following dynamic sub-system under the assumption that system (5) has a relative degree r = 1. _x r ¼ g þ Bu _g ¼ Uðx; uþ ð7þ ð8þ The uncertain term, g, is considered as a new state and U(x,u) is a non-linear unknown function that describes g s dynamics. We can note that the uncertain term g is observable since g ¼ _s Bu ð9þ

4 R. Aguilar-López, R. Martinez-Guerra / Chaos, Solitons and Fractals 33 (27) From the system given by Eqs. (7) and (8), it can be seen that a standard observer structure design, i.e. a system copy plus output feedback, is not possible to construct since the term U is unknown and considering that variable state x r is the system output, let us propose the following uncertainty observer. _^g ¼ s 1 ðg ^gþþs 2 ð_g _^gþ ð1þ This uncertainty estimator a reduced-order observer, which infers the uncertain term g is obtained from the corresponding state equation. Note that the observer contains proportional and derivative actions; the aim of the derivative term is to improve the speed of the estimation algorithm, because it enhances the anticipatory and stabilizing effects of derivative actions. Substituting the estimate of the uncertain term in the ideal controller defined by Eq. (6), the following non-ideal controller is obtained: u ¼ B 1 ½s g x r ^gš ð11þ From the corresponding state equation g ¼ _x r Bðx r Þu; substituting it into Eq. (11) and taking the time derivative of the resulting equation the following expression is obtained: _g ¼ x r s g _x r þ _^g ð12þ Introducing this result in Eq. (1) yields: _^g ¼ s 2 x r þðs 1 s 2 s g Þ_x r s g s 1 x r ð13þ Since this controller uses an estimated value of the uncertainty, it cannot cancel the system non-linearities completely. Practical stability is achieved as long as the uncertainty estimation error is bounded [16]. Thus, the system trajectories remain inside a neighborhood close to the defined trajectory of reference. The final expression for the input output linearizing controller with uncertainty estimation can be obtained integrating the estimator (Eq. (13)) and substituting it into the non-ideal controller (Eq. (11)) to obtain Z t u ¼ B 1 ðxþ ðs g ðs 1 s 2 s g ÞÞx r þ s g s 1 x r ðhþdh s 2 _x r ð14þ Note that this controller (Eq. (14)) exhibits PID structure and is equivalent to the linearizing controller based on proportional derivative uncertainty observer (Eqs. (1) and (11)) Stability issues of the proposed observer In order to show the stability properties of the closed-loop system, firstly, a convergence analysis of the uncertainty observer has to be done. Proposition 1. Let us define g and ^g be the uncertainty term and its estimate. The dynamic system _^g ¼ s 1 ðg ^gþþ s 2 ð_g _^gþ is an asymptotic-type reduced-order observer for the systems ((7) and (8)). Proof 1. Let us define the uncertainty estimation error as e ¼ g ^g ¼ Dg ð15þ Now, the dynamic scalar equation of the estimation error is given by Eq. (16), according to Eqs. (8) and (1) as follows. _e ¼ s 1 Uðx; uþ e þ ð16þ 1 þ s 2 1 þ s 2 Integrating this last expression, it renders e ¼ e exp s Z t 1 s 1 Uðx; uþ t þ exp ðt sþ ds ð17þ 1 þ s 2 1 þ s 2 1 þ s 2 Now, consider the following assumption: A1. U(x,u) is bounded, kuk 6 W, with < W < 1. Considering the norms of both sides of Eq. (17) kek 6 ke k exp s Z t 1 s 1 kuðx; uþk t þ exp ðt sþ ds 1 þ s 2 1 þ s 2 1 þ s 2

5 576 R. Aguilar-López, R. Martinez-Guerra / Chaos, Solitons and Fractals 33 (27) Applying the assumption A1, the following expression is obtained: kek 6 exp s 1 t ke k W þ W 1 þ s 2 s 1 s 1 ð18þ In the limit, when t!1: kek 6 W s 1 ð19þ It is important to analyze the structure of Eq. (18) in order to obtain some characteristics of the proposed observer. As desired, in the limit when t!1, the estimation error remains around a closed-ball with radius proportional to W/ s 1, which can be made as small as desired by taking s 1 large enough. In order to improve the speed to the uncertainty observer convergence to the steady-state estimation error, two actions were taken. The first one is to consider the parameter s 1 large enough, which is necessary to obtain a small steady-state error, as it was mentioned above. The second one is faced with the influence of the parameter s 2 related to the derivative action of the uncertainty observer. If s 2! 1, the exponential term of the right side of Eq. (18) can be accelerated enough and consequently, the convergence of the uncertainty observer will exhibit better performance. Note that if the measurements of the system are corrupted by additive noise, i.e. s = x r + m, and this noise is considered bounded, kmk 6 X, a methodology similar to the one used to analyze the estimation error e can be applied in order to prove that the steady-state estimation error becomes kek 6 WþX s 1. This confirms robustness against noisy measurements. Now, it is possible to implement a non-ideal controller using the estimated uncertainty. In order to prove closed-loop stability, it is necessary to analyze the closed-loop equation of the state equation when the non-ideal control law is introduced Stability issues of the proposed controller Consider the state equation of the measured signal of the master system _x m ¼ f m ðx m Þþ m ðx m Þ ð2þ With the following assumption: _x m < M < 1 and the closed-loop of the state equation of the slave signal _x s ¼ Dg s g ðx s x m Þ Defining the synchronization error as s x m The corresponding dynamic is as _ m ; x s ; uþ ð22þ where h ¼ f m ðx m Þþ m ðx m Þ Dg Now, for the above considerations the following assumptions are satisfied. There exist s g 2 R and N 2 R þ such that N A1. khðx m ; x s ; uþk 6 N and lim t!t s g ¼. A2. lim t!t k expð R t s g drþk ¼, where t is large enough. Eq. (22) can be solved to obtain the error dynamics in the time domain (Eq. (23)). Z Z exp s g dt þ expð s g ðt rþþhðx m ; x s ; uþdr ð23þ Taking norms of both sides of Eq. (23) yields to an inequality (expression (24)), which is limited after assumptions A1 and A2.

6 R. Aguilar-López, R. Martinez-Guerra / Chaos, Solitons and Fractals 33 (27) R t Z lim sup exp R s g dt hðxm ; x s ; uþdr t!t 6 lim sup knk 6 kn k lim sup exp s g dt þ t!t t!t lim sup exp R s g dt t!t lim sup N R t exp R s g dr t!t 6 lim sup knðtþk 6 t!t lim sup exp R s g dt t!t Expression (25) is an example of the 1/1 case of uniform L hôpital s rule, which can be applied to solve the undefined quotient (Eq. (25)). Now, taking the limit when t! t it is possible to obtain a bound for the estimation error (Eq. (26)). R Nexp sg dt 6 lim sup k@ðtþk 6 lim sup t!t t!t exp R N s g dt ¼ lim sup ð26þ ksg k t!t ks g k ð24þ ð25þ and k@k! ð27þ 4. Results and discussion In this section, the results from the numerical experiments carried off to show the goodness of the proposed synchronization methodology are presented. Both of the two cases are considered as application examples; the first one is related with the synchronization of the Van der Pol (master system) Lorenz (slave system). The mathematical model of both non-linear oscillators is Van der Pol model: _x 1;m ¼ x 2;m _x 2;m ¼ lð1 x 2 1 Þx 2 x 2 x 1 ax 3 1 kx5 1 þ f cos xt ð28þ s m ¼ x 2;m where l ¼ :4; x ¼ :46; a ¼ 1; k ¼ :1; x ¼ :86; f ¼ 4:5 with the corresponding initial conditions x 1;m ðþ ¼ x 2;m ðþ ¼ Lorenz model: _x 1;s ¼ rðx 2;s x 1;s Þþu 1 _x 2;s ¼ qx 1;s x 2;s x 1;s x 3;s _x 3;s ¼ x 1;s x 2;s bx 3;s ð29þ s s ¼ x 1;s where r ¼ 1; q ¼ 28; b ¼ 8 3 and the corresponding initial conditions x 1;s ðþ ¼1 x 2;s ðþ ¼ x 3;s ðþ ¼ 5 The corresponding uncertain term g 1 = rx 2,s such that x 2,m is an unmeasured variable, the control gain was chosen as s g = 1, the observer gains were considered as s 1 = 1 and s 2 =.8. Fig. 1 is related with the synchronization proce-

7 578 R. Aguilar-López, R. Martinez-Guerra / Chaos, Solitons and Fractals 33 (27) X 2,s X 2,m, Time Fig. 1. Synchronization of Van der Pol and Lorenz oscillators. Solid line master signal (x 2,m ) and dotted line slave signal (x 2,s ). dure, it was turned on at 25 min, where the slave signal tracks the master signal with satisfactory agreement, in accordance with the theory developed. On the other hand, the systems Duffing Chen whose mathematical models are the following: Duffing model: _z 1;m ¼ z 2;m _z 2;m ¼ dx 2;m p 2 x 1;m cx 3 1 þ K 1 cosðx 1 t þ h 1 ÞþK 2 cosðx 2 t þ h 2 Þ ð3þ s m ¼ x 2;m where d ¼ :2; p ¼ :5; c ¼ 1; K 1 ¼ 11; K 2 ¼ ; h 1 ¼ h 2 ¼ ; x 1 ¼ 1; x 2 ¼ 4 with initial conditions z 1;m ðþ ¼ z 2;m ðþ ¼ Chen model: _z 1;s ¼ aðz 2;s z 1;s Þ _z 2;s ¼ðc aþz 1;s cz 2;s z 1;s z 3;s þ u 2 _z 3;s ¼ z 1;s z 2;s bz 3;s ð31þ s s ¼ z 2;s with parameters a ¼ 35; b ¼ 3; c ¼ 28 and initial conditions z 1;s ¼ 1 z 2;s ¼ 1 z 3;s ¼ 1

8 R. Aguilar-López, R. Martinez-Guerra / Chaos, Solitons and Fractals 33 (27) Here the uncertainty considered is the non-linear term g 2 = z 1,s z 3,s, same as in the above case, the terms z 1,s and z 3,s are unmeasured variables, for this case the observer gains were considered as s 1 = 5 and s 2 =.9 and the controller gain was chosen as s g = 1. The synchronization procedure was started after 1 min; can be observed in Fig. 2 that the slave signal tracks the master signal as expected without any problem. Figs. 3 and 4 are related with the respective control input efforts; note that an acceptable effort is required for the synchronization task Z 2,s Z 2,m, Time Fig. 2. Synchronization of Duffing and Chen oscillators. Solid line master signal (z 2,m ) and dotted line slave signal (z 2,s ) (u 1 ) input Control Time Fig. 3. Control input effort (u 1 ) for the synchronization of Van der Pol and Lorenz oscillators.

9 58 R. Aguilar-López, R. Martinez-Guerra / Chaos, Solitons and Fractals 33 (27) Control (u 2 ) input Time Fig. 4. Control input effort (u 2 ) for the synchronization of Duffing and Chen oscillators. 5. Conclusions A methodology to synchronize different non-linear oscillators is proposed. The methodology synchronizes one slave signal to one master signal and it is related with an uncertainty proportional derivative reduced-order observer coupled to linearizing controller, this procedure generates a robust PID feedback structure. Note that the employed observer can be interpreted as a sub-slave of the slave system, such that it provides valuable information (about uncertain terms) to the slave system. The proposed methodology is applied to two examples, the Van der Pol Lorenz oscillators and the Duffing Chen oscillators, where the synchronization shows an adequate performance, as can be seen in the numerical simulations. References [1] Solak E. A reduced-order observer for the synchronization of Lorenz systems. Phys Lett A 24;325: [2] Morgül O, Solak E. Observed based synchronization of chaotic systems. Phys Rev E 1996;54(5): [3] Liu Ling, Wu Xiaogang, Hu Hanping. Estimating system parameters of Chua s circuit from synchronizing signal. Phys Lett A 24;324: [4] Chen Shihua, Hu Ja, Wang Changping, Lü Jinhu. Adaptive synchronization of uncertain Rössler hyper chaotic system based on parameter identification. Phys Lett A 24;321:5 5. [5] Feki M. Observer-based exact synchronization of ideal and mismatched chaotic systems. Phys Lett A 23;39:53 6. [6] González J, Femat R, Alvarez-Ramírez J, Aguilar R, Barrón MA. A discrete approach to the control and synchronization of a class of chaotic oscillators. IEEE Trans Circ Syst, Part I 1999;46(9): [7] Femat R, Alvarez-Ramírez J, Castillo-Toledo B, González J. On robust chaos suppression of nonlinear oscillators: application to Chua s circuit. IEEE Trans Circ Syst, Part I 1999;46(9): [8] Bowong S. Stability analysis for the synchronization of chaotic systems with different order: application to secure communications. Phys Lett A 24;326: [9] Moukam Kakmeni FM, Bowong Samuel, Tchawoua Clement, Kaptouom Ernest. Chaos control and synchronization of a U 6 -Van der Pol oscillator. Phys Lett A 24;322: [1] Muraly K. Synchronization based signal transmission with heterogeneous chaotic systems. Int J Bifurcat Chaos 2;1(11): [11] Femat R, Solís-Perales G. On the chaos synchronization phenomena. Phys Lett A 1999;262:5 6. [12] Ott E, Grebogi C, York JA. Controlling chaos. Phys Rev Lett 199;64(11):

10 R. Aguilar-López, R. Martinez-Guerra / Chaos, Solitons and Fractals 33 (27) [13] Isidori A. Nonlinear control systems. Berlin (Germany): Springer-Verlag; [14] Nijmeijer H, van der Shaft A. Non-linear dynamical control systems. Springer-Verlag: New York; 199. [15] Aguilar R, Martínez-Guerra R, Poznyak A, Maya-Yescas R. Temperature control in catalytic cracking reactors via robust PID controller. J Process Control 22;12(6): [16] Aguilar-López R, Alvarez-Ramírez J. Sliding-mode control scheme for a class of continuous chemical reactors. IEE Proc Control Theory Appl 22;149(4): [17] Martínez-Guerra Rafael, Aguilar Ricardo, Poznyak Alexander. A new robust sliding-mode observer design for monitoring in chemical reactors. J Dyn Syst, Measure Control, ASME J 24;126(3): [18] Aguilar Ricardo, Martínez Rafael, Poznyak Alexander. PI observers for uncertainty estimation in continuous chemical reactors. In: Proceedings of the IEEE conference on control applications/isic 21, Ciudad de México, Mexico, September 21. p