Preservation of Stability and Synchronization of a Class of Fractional-Order Systems

Abstract

We present sufficient conditions for the preservation of stability of fractional-order systems, and then we use this result to preserve the synchronization, in a master-slave scheme, of fractional-order systems. The systems treated herein are autonomous fractional differential linear and nonlinear systems with commensurate orders lying between 0 and 2, where the nonlinear ones can be described as a linear part plus a nonlinear part. These results are based on stability properties for equilibria of fractional-order autonomous systems and some similar properties for the preservation of stability in integer order systems. Some simulation examples are presented only to show the effectiveness of the analytic result.

(ProQuest: ... denotes non-US-ASCII text omitted.)

Recommended by Ricardo Femat

1, Departamento de Física y Matemáticas, Universidad Iberoamericana, Prol. Paseo de la Reforma 880, Lomas de Santa Fe, 01210 México, DF, Mexico

Received 18 April 2012; Revised 11 August 2012; Accepted 12 August 2012

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1. Introduction

The applications of fractional calculus to science and engineering have been growing in the last few years [ 1]; this is due in part to the properties of these operators. The applications, specifically, that involve fractional-order chaotic systems or their synchronization had been one of the principal subjects of investigation; some of these works are [ 2- 6]. There is also several works concerning chaotic systems or complex networks of integer order or their synchronization, for example, [ 7- 10]. There are many different works on the synchronization of fractional autonomous systems that can be described as a linear plus a nonlinear part [ 11- 14], in such works several schemes are proposed to ensure that the error dynamics satisfies the conditions from the celebrated theorem for autonomous commensurate differential systems with fractional order between 0 and 1 by [ 15]; this means that the error dynamics must hold a linear relation in order to achieve the synchronization. There is also a scheme proposed in [ 16], based on [ 17], where the dynamical system of the synchronization error can be nonlinear, which when viewed from the analytical point can be important because it does not restrict the error dynamics to be only linear.

There is some other interesting theme, the preservation of stability and synchronization, which is the main issue in this work. This problem can be stated as follows: if we have an original autonomous nonlinear system that can be described as a linear plus a nonlinear part whose origin is stable, we want to investigate some kinds of modifications that can occur to the fractional order, the linear part, and the nonlinear part in such a way that the origin of the modified system is also stable. This subject is important because such modifications can be interpreted as perturbations on the system. Note that the modification of the linear part of the vector field associated with the fractional differential equation modifies some local properties of the vector field at the point of equilibrium, in particular local stability. In [ 18], the authors developed two results for the preservation of stability of integer-order nonlinear systems; one of such results gives conditions for the preservation of stability between systems of different orders of the state vector but does not give direct insight on the transformations, and the other result gives more insight but in return is a little more restrictive because as part of the hypothesis it asks for diagonalizability of the linear part of the system. In [ 19], the authors have reached conditions for the preservation of stability for integer-order systems in the presence of nonlinear modifications to the Jacobian matrix; such modifications can be applied on the characteristic polynomial or in form of a nonlinear polynomial matrix evaluation.

The main objective of this work is to state under which conditions a certain family of transformations applied to the fractional order, the linear part, and the nonlinear part of an autonomous fractional differential system with commensurate order will preserve stability of the origin. It is important to point out that this analytical result is of relevance for its relation with robustness not for the use of an advanced controller in the stabilization or the synchronization. As far as the authors know, this problem has not been addressed for the case of fractional-order systems.

In Section 2, we present the definitions and some results on the stability of autonomous commensurate fractional-order systems. In Section 3, the main results are stated in form of propositions and corollaries. Based on these propositions, in Section 4, we present a methodology to illustrate how these results can be used and this is complemented by the application of this methodology in two examples of simulation presented in Section 5. Finally in Section 6, we present the obtained conclusions.

2. Preliminary Results

There are several definitions of a fractional derivative of order α ∈ ^{... +} [ 20- 22]. We will use the Caputo fractional operator because the meaning of the initial conditions for systems described using this operator is the same as for integer-order systems.

Definition 2.1 (Caputo fractional derivative).

The Caputo fractional derivative of order α ∈ ^{... +} of a function x is defined as (see [ 20]) [figure omitted; refer to PDF] where m -1 ...4; α <m , ^{d m} x ( τ ) /d ^{τ m} is the m th derivative of x in the usual sense, m ∈ ... , and Γ is the gamma function. (Throughout the paper, we use indistinctly ^{x ( α )} ...1; ^{x ( α )} ( t ) , x ...1;x ( t ) .)

We recall some previous results on the stability of autonomous commensurate fractional-order systems that are related to our study.

2.1. Autonomous Commensurate Fractional-Order Linear Systems Stability

Given an autonomous fractional-order system with state space representation [figure omitted; refer to PDF] where A ∈ ^{R n ×n} , B ∈ ^{R n ×m} , C ∈ ^{R p ×n} , the state vector x ∈ ^{R n} , the input vector u ∈ ^{R m} , and the output vector y ∈ ^{R p} .

Definition 2.2 (see [ 15]).

The fractional-order autonomous system ( 2.2) [figure omitted; refer to PDF] is said to be

(i) stable if and only if for all ... >0 ∃ δ = δ ( ... ) >0 , such that given || _{x 0} || < δ then ||x ( t ) || < ... for all t ...5;0 ;

(ii) asymptotically stable if and only if it is stable and _{lim t [arrow right] ∞} ||x ( t ) || =0 .

Firstly, we will introduce some results on fractional-order systems stability. First for 0 < α <1 , we have the celebrated Theorem [ 15] that gives us necessary and sufficient conditions for the asymptotic stability of the origin of a type of autonomous linear fractional-order systems; such conditions involve the argument of the eigenvalues of the system matrix.

Theorem 2.3.

The autonomous system [figure omitted; refer to PDF] is asymptotically stable if and only if |arg (spec (A ) ) | > α ( π /2 ) , where spec (A ) is the set of all the eigenvalues of A . Also, the state vector x decays towards 0 and meets the following condition: ||x || <N ^{t - α} , t >0 , α >0 .

And for 1 < α <2 , we have a similar result [ 23].

Theorem 2.4.

The autonomous fractional differential system [figure omitted; refer to PDF] with initial conditions ^{x (k )} ( _{t 0} ) = _{x k} (k =0 ,1 ) , with the Caputo derivative and where x ∈ ^{... n} , A ∈ ^{... n ×n} is asymptotically stable if and only if |arg (spec (A ) ) | > α ( π /2 ) . In this case, the components of the state decay towards 0 like ^{t - α -1} . Moreover, the system ( 2.5) is stable if and only if either it is asymptotically stable, or those eigenvalues which satisfy |arg (spec (A ) ) | = α ( π /2 ) have the same algebraic and geometric multiplicities.

2.2. Commensurate Fractional-Order Nonlinear Systems Stability

Given a commensurate fractional-order system with the Caputo fractional operator [figure omitted; refer to PDF] with initial condition x ( _{t 0} ) = _{x 0} , α ∈ (0,1 ) , f : [ _{t 0} , ∞ ) × Ω [arrow right] ^{... n} is piecewise continuous in t and locally Lipschitz in x on [ _{t 0} , ∞ ) × Ω , and Ω ⊂ ^{... n} is a domain that contains the origin x =0 .

The equilibrium point of ( 2.6) is defined as follows [ 24].

Definition 2.5.

The constant _{x e} is an equilibrium point of the fractional-order system ( 2.6) if and only if f (t , _{x e} ) =0 .

Without loss of generality, let the equilibrium point be x = _{x e} =0 . In this definition, we are considering that the result of the derivative of a constant is zero because we are using only the Caputo fractional operator.

Definition 2.6 (the Lyapunov stability).

The equilibrium point x =0 of the system ( 2.6) is said to be

(1) stable, if for all ... >0 ∃ δ >0 such that if || _{x 0} || < δ then ||x || < ... , for all t ...5;0 . Otherwise the equilibrium point is called unstable;

(2) asymptotically stable, if it is stable and in addition the following equality holds: [figure omitted; refer to PDF]

As a starting point for the construction of our own results, we can use the following result for the stability of the origin of commensurate fractional-order systems with 0 < α <1 [ 17].

Theorem 2.7.

Consider the n - dimensional nonlinear fractional-order dynamic system [figure omitted; refer to PDF] with a constant linear regular matrix A , a nonlinear function g (x ) of the states x , and 0 < α <1 . If

(1) the zero solution of ^{x ( α )} =Ax is asymptotically stable and α ρ (A ) >1 ;

(2) g (0 ) =0 and _{lim ||x || [arrow right]0} || ( g (x ) || / ||x || ) =0 , where ρ (A ) is the spectral radius of A ,

then x =0 , 0 ...4; _{t 0} ...4;t is a stable solution of the system ( 2.8).

The following result is valid for the asymptotic stability of systems with 1 < α <2 . Consider the n - dimensional nonlinear fractional-order dynamic system with the Caputo derivative [figure omitted; refer to PDF] under the initial conditions [figure omitted; refer to PDF] where x ∈ ... , matrix A ∈ ^{... n ×n} , and 1 < α <2 , g (t ,x ) : [ _{t 0} , ∞ ) × ^{... n} [arrow right] ^{... n} is a continuous function in which g (t ,0 ) =0 ; moreover, g (t ,x ) holds the Lipschitz condition with respect to x .

Theorem 2.8.

If the matrix A such that |arg (spec (A ) ) | ...0;0 , |arg (spec (A ) ) | > α ( π /2 ) , α +1 / ||A || <2 , and suppose that the function g (t ,x ) satisfies uniformly [figure omitted; refer to PDF] then the zero solution of ( 2.9) is asymptotically stable.

The proof of this theorem for the Caputo derivative follows from the proof of Theorem 3.3 in [ 23] and the application of Lemma 2.7 in [ 23] and Gronwall-Bellman inequality.

3. Preservation of Stability

So once given all these stability results, we need to give a definition for the preservation of stability in fractional-order systems in order to be in the possibility to state the conditions in form of a proposition.

Definition 3.1.

Given an asymptotically stable autonomous commensurate fractional-order linear system of the kind [figure omitted; refer to PDF] where A ∈ ^{... n ×n} , x ∈ ^{... n} , 0 < α <2 and A =P _{J A} ^{P -1} . If one has a transformation ψ : ^{R +} × ^{R n ×n} [arrow right] ^{R +} × ^{R n ×n} , namely, ψ ( α ,A ) = ( α β ,MA ) , such that the new system [figure omitted; refer to PDF] is also asymptotically stable, where MA =P _{J M J A} ^{P -1} , M ∈ ^{... n ×n} , for some matrix M =P _{J M} ^{P -1} , where _{J M} and _{J A} are Jordan matrices, then one says that ψ is an asymptotically stability preserving transformation for commensurate fractional-order autonomous linear systems.

We should notice that for the matrices M and A we are using the same matrix P . It is also worth to mention that, given the Jordan matrix _{J A} that corresponds to the matrix A , the Jordan matrix _{J M} , that represents the modifications, must have the same order and type of Jordan blocks that _{J A} . The reason behind this fact is that in several applications we will need that MA ∈ ^{... n ×n} . But it should also be noticed that another canonical form, instead of Jordan form, could result more convenient in the construction of state feedback controllers.

Now in a similar way we state the definition for this concept in commensurate fractional-order nonlinear systems.

Definition 3.2.

Given a commensurate fractional-order nonlinear system of the kind [figure omitted; refer to PDF] where A ∈ ^{... n ×n} , A =P _{J A} ^{P -1} , x ∈ ^{... n} , 0 < α <2 , g : [ _{t 0} , ∞ ) × Ω [arrow right] ^{... n} is piecewise continuous in t and locally Lipschitz in x on [ _{t 0} , ∞ ) × Ω , and Ω ⊂ ^{... n} is a domain that contains the origin and the origin itself is a stable solution of the system. If one has a transformation Ψ : ^{R +} × ^{R n ×n} × ^{C k} ( ^{R n} , ^{R n} ) [arrow right] ^{R +} × ^{R n ×n} × ^{C k} ( ^{R n} , ^{R n} ) , namely, Ψ ( α ,A ,g ( · ) ) = ( α β ,MA ,cg ( · ) ) , in such a way that in the new system [figure omitted; refer to PDF] the origin is also a stable solution, where c ∈R , MA =P _{J M J A} ^{P -1} , M ∈ ^{... n ×n} , for some matrix M =P _{J M} ^{P -1} , where _{J M} and _{J A} are Jordan matrices, then one calls to that transformation a stability preserving transformation for commensurate fractional-order nonlinear systems.

Remark 3.3.

In Definition 3.2, for the case where 0 < α <1 , the nonlinear part is considered as autonomous, that is, for the system ( 3.3), we have ^{x ( α )} =Ax +g (x ) , and for the modified system ( 3.4), we have ^{x ( α β )} =MAx +cg (x ) .

Now based on the Theorems 2.3, 2.4, 2.7, and 2.8, and the results from [ 18] for the preservation of stability for integer-order systems, the following criterion for the preservation of stability in autonomous commensurate fractional-order systems can be stated as follows.

Proposition 3.4.

Consider an autonomous commensurate fractional-order nonlinear system of the form [figure omitted; refer to PDF] with x ∈ ^{... n} , A ∈ ^{... n ×n} , g :D ⊂ ^{... n} [arrow right] ^{... n} is a continuous function, D is a neighborhood of the origin for 0 < α <1 . Let A ∈ ^{... n ×n} with the argument of its kth eigenvalue denoted by _{θ k} =arg ( _{λ k} (A ) ) . Given a transformation Ψ ( α ,A ,g ( · ) ) = ( α β ,MA ,cg ( · ) ) such that the new system is [figure omitted; refer to PDF] where c ∈R , M ∈ ^{... n ×n} , 0 < β ...4;1 , _{[varphi] k} =arg ( _{λ k} (M ) ) is the argument of the kth eigenvalue of M , A =P _{J A} ^{P -1} , M =P _{J M} ^{P -1} . Also let _{[varphi] a k} = - _{θ k} + α ( π /2 ) , _{[varphi] b k} = - _{θ k} - α π /2 , _{[varphi] max} = _{max k} { _{[varphi] a k} } , _{[varphi] min} = _{min k} { _{[varphi] b k} } , if [figure omitted; refer to PDF] for each k =1 ,2 , ... ,n , and if the system ^{x ( α )} =Ax is asymptotically stable, α ρ (A ) >1 , g (0 ) =0 , _{lim ||x || [arrow right]0} ||g (x ) || / ||x || =0 , and ρ (MA ) ...5; ρ (A ) , then one claims that such transformation is a stability preserving transformation for fractional-order systems of the kind of ( 3.5).

Proof.

Summarizing the initial hypothesis, the original system ( 3.5) holds the conditions from Theorem 2.7, so we have | _{θ k} | > α ( π /2 ) for k =1 ,2 , ... ,n .

By the hypothesis ρ (MA ) ...5; ρ (A ) , ρ (A ) >1 , and 0 < β ...4;1 , we have that α β ρ (MA ) >1 , and we have asked for g (x ) to hold g (0 ) =0 and _{lim ||x || [arrow right]0} ||g (x ) || / ||x || =0 . As a result we need the asymptotic stability of the system ^{x ( α β )} =MAx to hold the conditions of Theorem 2.7for the new system ( 3.6).

By the properties of the complex numbers, and based on the fact that _{J M} and _{J A} are Jordan matrices with the same structure and that MA =P _{J M J A} ^{P -1} , in order to assure that the system ^{x ( α β )} =MAx is asymptotically stable, we need for |arg (spec (MA ) ) | > α β ( π /2 ) to hold, so first we want for [figure omitted; refer to PDF] to hold. The last part of the hypothesis states that the inequality ( 3.7) holds. From the right part of ( 3.7), we know that given that each _{[varphi] k} is greater than _{[varphi] max} , we have that _{[varphi] k} > - _{θ k} + α ( π /2 ) . And similarly from the left part we have that - _{θ k} - α ( π /2 ) > _{[varphi] k} for any _{[varphi] k} . Then these two parts together give us precisely that ( 3.8) holds, and taking from the hypothesis that 0 < β ...4;1 , we have that α β ...4; α and therefore | _{[varphi] k} + _{θ k} | > α β ( π /2 ) , and thus the modified system holds all the conditions for the linear part from Theorem 3. From the demonstration of Theorem 3 given in [ 17], we can observe that cg (x ) also holds the corresponding conditions; therefore we claim that Ψ is a stability preserving transformation for the fractional-order autonomous systems of the form of ( 3.5).

Now we have a similar result for systems with fractional orders lying between 1 and 2.

Proposition 3.5.

Consider a partially autonomous commensurate fractional-order nonlinear system of the form [figure omitted; refer to PDF] with 1 < α <2 , x ∈ ^{... n} , A ∈ ^{... n ×n} , g : [ _{t 0} , ∞ ) ×D ⊂ ^{... n} [arrow right] ^{... n} is a continuous function, D is a neighborhood of the origin, and g (x ) holds the Lipschitz condition with respect to x .

Let A ∈ ^{... n ×n} with the argument of its kth eigenvalue denoted by _{θ k} =arg ( _{λ k} (A ) ) . Given a transformation Ψ ( α ,A ,g (t , · ) ) = ( α β ,MA ,cg (t , · ) ) such that the new system is [figure omitted; refer to PDF] where c ∈R , 1 / α < β <2 / α , M ∈ ^{... n ×n} , _{[varphi] k} =arg ( _{λ k} (M ) ) is the argument of the kth eigenvalue of M , A =P _{J A} ^{P -1} , M =P _{J M} ^{P -1} . Also let _{[varphi] a k} = - _{θ k} + α π /2 , _{[varphi] b k} = - _{θ k} - α π /2 , _{[varphi] max} = _{max k} { _{[varphi] a k} } , _{[varphi] min} = _{min k} { _{[varphi] b k} } , if [figure omitted; refer to PDF] for each k =1 ,2 , ... ,n , and if |arg (spec (A ) ) | ...0;0 , | arg(spec(A )) | > α ( π /2 ) , α +1 / ||A || <2 , β ...4;1 , β < α (2 -1 / ||MA || ) , _{lim x [arrow right]0} ||g (t ,x ) || / ||x || =0 is uniformly satisfied for t ∈ [ _{t 0} , ∞ ) , ||MA || ...5; ||A || , and |arg (spec (M ) ) | ...0;0 , then one claims that such transformation is a stability preserving transformation for fractional-order systems of the kind of ( 3.9).

Proof.

The first part of the hypothesis is that the original system ( 3.9) holds the conditions of Theorem 2.8; therefore, we have that | _{θ k} | > α ( π /2 ) , for k =1 ,2 , ... ,n , and that its origin is an asymptotically stable solution.

Another part of the hypothesis is that g (t ,x ) holds the conditions from Theorem 2.8, and from the proof of Theorem 2.8, we observe that such conditions also hold for cg (t ,x ) .

We have also asked for ||MA || ...5; ||A || , α +1 / ||A || <2 , and β < α (2 -1 / ||MA || ) to hold, and thus we have that α β +1 / ||MA || <2 . As a result, we only need |arg (spec (MA ) ) | > α β ( π /2 ) to hold the conditions from Theorem 2.8for the new system ( 3.10).

By very similar arguments as in the proof of Proposition 3.4, we have that ( 3.8) holds and we have asked for β ...4;1 to hold; therefore conditions from Theorem 2.8are satisfied and we claim that Ψ is a stability preserving transformation for the fractional-order autonomous systems of the form of ( 3.5).

The following corollaries are a direct consequence of Proposition 1 and Theorem 3.3 in [ 23].

Corollary 3.6.

Consider a linear autonomous commensurate fractional-order system of the form [figure omitted; refer to PDF] that holds the conditions from Theorem 2.4or Theorem 2.3accordingly, if one has a transformation Ψ ( α ,M ) = ( α β ,MA ) such that the new system is [figure omitted; refer to PDF] with all the variables and matrices as defined before, and if the inequality [figure omitted; refer to PDF] holds, with _{[varphi] k} , _{[varphi] min} , and _{[varphi] max} as defined before, then one claims that Ψ is an asymptotic stability preserving transformation for the autonomous commensurate fractional-order linear systems of the form of ( 3.12).

Remark 3.7.

If we take the case of α =1 for Proposition 3.4, or Proposition 3.5we will have similar conditions for the preservation of stability of the origin for an integer-order nonlinear system x =Ax +g (x ) . And if we take the case of α =1 for Corollary 3.6, we will have the conditions for the preservation of asymptotic stability of an integer-order linear system x =Ax , without the modification over the order, of course.

Remark 3.8.

Notice that the conditions from [ 18, 19] are different from the ones that we obtain choosing α =1 in Proposition 3.4and cannot be obtained as a particular case of Proposition 3.4.

4. Preservation of Synchronization

Given the conditions from Proposition 3.4, we want to illustrate how to use these criteria to ensure the preservation of stabilization and synchronization.

4.1. Preservation of Stabilization of Autonomous Commensurate Fractional-Order Nonlinear Systems

First for the stabilization of a system, we know from previous works that for an autonomous fractional commensurate order system of the form ^{x ( α )} =Ax +g (x ) +u , where g (x ) holds, the conditions from Theorem 2.7, and with A =P _{J A} ^{P -1} , we can choose a control u = - _{K 1} x , _{K 1} ∈ ^{R n ×n} , with _{K 1} ∈ ^{... n ×n} in such a way that for the system ^{x ( α )} = (A - _{K 1} )x +g (x ) , the origin is a stable solution. But in this particular case we want that A - _{K 1} =P ( _{J A} - _{J K 1} ) ^{P -1} ; therefore, we need to construct _{K 1} as _{K 1} =P _{J K 1} ^{P -1} ; the Jordan form _{J K 1} also has the restriction that its Jordan blocks must be of the same order and type as the ones of _{J A} .

Now for the modified system ^{x ( α β )} =MAx +cg (x ) +u , we will use a similar control (i.e., with the same _{K 1} ) defined as u = -M _{K 1} x , where M : =P _{J M} ^{P -1} holds the conditions from 1 and all the other matrices defined as before, in such a way that for the system ^{x ( α β )} = (MA -M _{K 1} )x +cg (x ) , the origin is also a stable solution.

4.2. Preservation of Complete Practical Synchronization of Autonomous Commensurate Fractional-Order Nonlinear Systems

First we need to describe the synchronization scheme. Let us consider two fractional-order systems as the master and the slave system, respectively, [figure omitted; refer to PDF] where _{x [physics M-matrix]} ∈ ^{R n} is the state vector of the master system, _{y [physics M-matrix]} : ∈ ^{R p} is the output of the master system, _{x ...AE;} ∈ ^{R n} is the state vector of the slave system, w ∈ ^{R n} is the control input, and _{y ...AE;} ∈ ^{R p} is the output of the slave system.

In this synchronization scheme, the master system represents the target dynamics, while the slave system represents the system to be controlled.

Let us consider that all the outputs are available, only to illustrate the effectiveness of the method by showing more states on the graphs of the synchronization error in the examples, of course the order of the output can be less than the order of the state vector. This consideration will lead us to the synchronization error: [figure omitted; refer to PDF] we must find a function w such that ||e (t ) || is bounded in a subset that contains the origin, because the result presented in [ 17] only gives conditions for the stability of the origin, not for asymptotic stability. Because of this fact, this approach is called complete practical synchronization.

We specifically choose the control [figure omitted; refer to PDF] where _{K 2} ∈ ^{R n ×n} , in such a way that the error dynamics is given by [figure omitted; refer to PDF]

As before, given that _{A ...AE;} = _{P A ...AE; J A ...AE;} ^{P _{A ...AE;} -1} , we will have _{A ...AE;} - _{K 2} = _{P A ...AE;} ( _{J A ...AE;} - _{J K 2} ) ^{P _{A ...AE;} -1} , and therefore we need to construct _{K 2} as _{K 2} = _{P A ...AE; J K 2} ^{P _{A ...AE;} -1} in such a way that the origin of the dynamic system of the error is stable. Again the Jordan form _{J K} has the restriction that its Jordan blocks must be of the same order and type as the ones of _{J A ...AE;} .

Then, we want to illustrate what kind of transformations can be applied to the master and slave systems in such a way that the same _{K 2} still stabilizes the origin of the modified synchronization error system. And to do this we define the modification matrices _{M [physics M-matrix]} : = _{P A [physics M-matrix] J M} ^{P _{A [physics M-matrix]} -1} and _{M ...AE;} : = _{P A ...AE; J M} ^{P _{A ...AE;} -1} that hold the conditions from Proposition 3.4. With these modifications applied to each of the systems in the following way: [figure omitted; refer to PDF] with the synchronization error defined as in ( 4.2), and the control defined as [figure omitted; refer to PDF] (the only change is that instead of the term _{K 2} e now is _{M ...AE; K 2} e ), one has that the autonomous commensurate fractional order dynamical system of the error is [figure omitted; refer to PDF] Given that one has constructed the matrix M in such a way that it holds the conditions from Proposition 3.4it is straightforward to prove that the origin of the new dynamic system of the error is also stable. Note that the modification of the linear part of the vector field associated to the fractional differential equation, modifies the manifold of synchronization, but not the stability.

5. Examples

5.1. Preservation of Stabilization for Fractional-Order Lorenz Systems

Let us take the Lorenz system with commensurate fractional-order α =0.97 that can be written as [ 5]: [figure omitted; refer to PDF] with x = [ _{x 1 x 2 x 3} ^{] T} , σ =10 , ρ =28 , β - =8 /3 , initial conditions x (0 ) = [ -9 -5 14 ^{] T} . This system is chaotic as it is claimed in [ 25].

The objective is to stabilize the system and then apply a modification that holds the conditions of Proposition 3.4, to illustrate the validity of the analytical results.

In order to do this the next, we choose u = - _{K 1} x . With [figure omitted; refer to PDF] the eigenvalues of the new matrix A - _{K 1} are _{λ 1} [approximate] -5.827 , _{λ 2} [approximate] -3.172 , and _{λ 3} [approximate] -6.666 , with this and the fact that g (x ) holds the conditions from Theorem 2.7(as it has been demonstrated in [ 17]), we know that the origin of the controlled system is a stable solution.

Now we are interested in verifying what will happen if we propose a modification such as α β =0.95 , β [approximate]0.97938 , c =0.8 , and [figure omitted; refer to PDF] it can be easily verified that M holds the conditions of Proposition 3.4. The eigenvalues of the modified system M (A - _{K 1} ) are _{λ 1} [approximate] -64.105 , _{λ 2} [approximate] -28.550 , and _{λ 3} [approximate] -69.333 . Given that this eigenvalues hold the conditions of Theorem 2.7and that 0.8g (x ) also holds the rest of the conditions, the origin of the modified controlled system is also a stable solution.

In Figures 1and 4, the simulation step and time were 0.004 and 50 s, respectively. In the first 25 s, the system was the original one ( u =0 ), and for the last 25 s the control u was activated, for the unmodified and modified systems (Figure 2).

Graphs of the states versus time of the controlled systems.

(a) The unmodified controlled system ^{x ( α )} = (A -K )x +g ( x )

[figure omitted; refer to PDF]

(b) The modified controlled system ^{x ( α )} =M (A -K )x +g (x )

[figure omitted; refer to PDF]

Phase plane of the unmodified and modified systems.

(a) The system ^{x ( α )} = (A -K )x +g (x )

[figure omitted; refer to PDF]

(b) The system ^{x ( α )} =M (A -K )x +g (x )

[figure omitted; refer to PDF]

5.2. Preservation of Complete Practical Synchronization for Fractional-Order Chen Systems

In this example specifically, we will make the synchronization of two fractional-order Chen systems with identical parameters and orders but different initial conditions, and because of these consideration, we will have that _{A [physics M-matrix]} = _{A ...AE;} , and therefore we drop the indices for the matrices A and M .

We use the structure of the Chen system with fractional-order α as presented in [ 3]. For the slave system, we have [figure omitted; refer to PDF] with w as defined in ( 4.6), _{x ...AE;} = [ _{x S 1 x S 2 x S 3} ^{] T} , a =35 , b =3 , d =28 , and α =0.975 , initial conditions _{x ...AE;} (0 ) = [3 0 10 ^{] T} .

Now, for the master system, we take [figure omitted; refer to PDF] with _{x [physics M-matrix]} = [ _{x M 1 x M 2 x M 3} ^{] T} , a =35 , b =3 , d =28 , α =0.975 , and initial conditions _{x [physics M-matrix]} (0 ) = [ -9 -5 14 ^{] T} .

Both master and slave systems are chaotic, as far as the conditions from [ 25]. For the control law ( 4.3), we have [figure omitted; refer to PDF] and for the modification we take c =0.9 , α β =0.96 , β [approximate]0.9846 , and [figure omitted; refer to PDF]

With _{K 2} as proposed, it is assured that the synchronization error system without modification holds the conditions of Theorem 2.7. It can be easily verified that M and 0.9g (x ) hold the conditions from Proposition 3.4, in such a way that we can assure that this is an example of preservation of synchronization.

Again the simulation step and time were 0.0028 and 50 s, respectively. In the first 25 s the system was the original one ( w =0 ), and for the last 25 s, the control w was activated, for the unmodified and modified systems.

In Figure 3, we can observe how the application of the control law, with the same _{K 2} , stabilizes the origin of the unmodified and the modified synchronization error systems.

Errors graph with the control law applied at t =25 s .

(a) Synchronization error of the systems without modification

[figure omitted; refer to PDF]

(b) Synchronization error of the modified systems

[figure omitted; refer to PDF]

Phase plane of the unmodified and modified synchronizations; the master systems are in blue, and the slave systems in black; the control law is applied at t =50 s .

(a) The unmodified synchronization

[figure omitted; refer to PDF]

(b) The modified synchronization

[figure omitted; refer to PDF]

In Figure 4, we can observe the way the slave system follows to the master system for the unmodified and the modified systems with the same _{K 2} . All the simulations were made using the algorithms presented in [ 26].

From several simulations, we have observed that under large variations in the parameters these transformations do not preserve chaos. This limits the possible variations in the transformations because, as it is well known, the chaos in dynamical systems is very sensitive to variations in the parameters.

6. Conclusions

As far as the authors know, this is the first time that the preservation of autonomous commensurate fractional order systems stability is made considering transformations that affect the fractional order, the linear part, and the nonlinear part of the vector field of the differential equation.

Furthermore, we have explained how these results can be used to ensure the preservation of stabilization and the preservation of synchronization of autonomous commensurate fractional-order systems and, through the presented examples, we have also showed the effectiveness of the results.

It is also worth to mention that there are some other results on the stability of fractional-order autonomous systems [ 25, 27, 28] that can be used in a similar way to Proposition 3.4to state the conditions for the preservation of asymptotic stability of the solutions, for 0 < α <1 .

References

[1] R. Hilfer Applications of Fractional Calculus in Physics , World Scientific Publishing, River Edge, NJ, USA, 2000.

[2] J. Lü, G. Chen, "Generating multiscroll chaotic attractors: theories, methods and applications," International Journal of Bifurcation and Chaos in Applied Sciences and Engineering , vol. 16, no. 4, pp. 775-858, 2006.

[3] C. Li, G. Chen, "Chaos in the fractional order Chen system and its control," Chaos, Solitons and Fractals , vol. 22, no. 3, pp. 549-554, 2004.

[4] C. Li, G. Chen, "Chaos and hyperchaos in the fractional-order Rössler equations," Physica A , vol. 341, no. 1-4, pp. 55-61, 2004.

[5] I. Grigorenko, E. Grigorenko, "Chaotic dynamics of the fractional Lorenz system," Physical Review Letters , vol. 91, no. 3, 2003.

[6] T. T. Hartley, C. F. Lorenzo, H. K. Qammer, "Chaos in a fractional order Chua's system," IEEE Transactions on Circuits and Systems I , vol. 42, no. 8, pp. 485-490, 1995.

[7] J. Lü, S. Yu, H. Leung, G. Chen, "Experimental verification of multidirectional multiscroll chaotic attractors," IEEE Transactions on Circuits and Systems I , vol. 53, no. 1, pp. 149-165, 2006.

[8] J. Lü, G. Chen, "A time-varying complex dynamical network model and its controlled synchronization criteria," IEEE Transactions on Automatic Control , vol. 50, no. 6, pp. 841-846, 2005.

[9] J. Zhou, J.-A. Lu, J. Lü, "Adaptive synchronization of an uncertain complex dynamical network," IEEE Transactions on Automatic Control , vol. 51, no. 4, pp. 652-656, 2006.

[10] J. Lü, X. Yu, G. Chen, D. Cheng, "Characterizing the synchronizability of small-world dynamical networks," IEEE Transactions on Circuits and Systems I , vol. 51, no. 4, pp. 787-796, 2004.

[11] Y. Tang, J.-A. Fang, "Synchronization of N -coupled fractional-order chaotic systems with ring connection," Communications in Nonlinear Science and Numerical Simulation , vol. 15, no. 2, pp. 401-412, 2010.

[12] G. Peng, Y. Jiang, F. Chen, "Generalized projective synchronization of fractional order chaotic systems," Physica A , vol. 387, no. 14, pp. 3738-3746, 2008.

[13] X. Y. Wang, J. M. Song, "Synchronization of the fractional order hyperchaos Lorenz systems with activation feedback control," Communications in Nonlinear Science and Numerical Simulation , vol. 14, no. 8, pp. 3351-3357, 2009.

[14] J. G. Lu, "Chaotic dynamics of the fractional-order Lü system and its synchronization," Physics Letters, Section A , vol. 354, no. 4, pp. 305-311, 2006.

[15] D. Matignon, "Stability results for fractional differential equations with applications to control processing," in Proceedings of the Computational Engineering in System Application, pp. 963-968, 1996.

[16] R. Martínez-Martínez, H. Jardón-Kojakhmetov, J. A. León, G. Fernández-Anaya Estabilización de Redes Complejas Fraccionarias de Sistemas de Lorenz y Sistemas de Chen , AMCA, 2009.

[17] X.-J. Wen, Z.-M. Wu, J.-G. Lu, "Stability analysis of a class of nonlinear fractional-order systems," IEEE Transsactions on Circuits and Systems-II , vol. 55, pp. 1178-1182, 2008.

[18] G. Fernández-Anaya, J. J. Flores-Godoy, R. Femat, J. J. Álvarez-Ramírez, "Preservation of stability and synchronization in nonlinear systems," Physics Letters, Section A , vol. 371, no. 3, pp. 205-212, 2007.

[19] D. Becker-Bessudo, G. Fernandez-Anaya, J. J. Flores-Godoy, "Preserving synchronization under matrix product modifications," Physica A , vol. 387, no. 26, pp. 6631-6645, 2008.

[20] I. Podlubny Fractional Differential Equations , vol. 198, Academic Press, San Diego, Calif, USA, 1999., 9th.

[21] K. S. Miller, B. Ross An Introduction to the Fractional Calculus and Fractional Differential Equations , John Wiley & Sons, New York, NY, USA, 1993.

[22] I. Podlubny, "Geometric and physical interpretation of fractional integration and fractional differentiation," Fractional Calculus & Applied Analysis , vol. 5, no. 4, pp. 367-386, 2002.

[23] F. Zhang, C. Li, "Stability analysis of fractional differential systems with order lying in (1 ,2 ) ," Advances in Difference Equations , vol. 2011, 2011.

[24] Y. Li, Y. Chen, I. Podlubny, "Mittag-Leffler stability of fractional order nonlinear dynamic systems," Automatica , vol. 45, no. 8, pp. 1965-1969, 2009.

[25] M. S. Tavazoei, M. Haeri, "Chaotic attractors in incommensurate fractional order systems," Physica D. Nonlinear Phenomena , vol. 237, no. 20, pp. 2628-2637, 2008.

[26] I. Petras Fractional Order Chaotic Systems, 2010, http://www.mathworks.com/matlabcentral/fileexchange/27336-fractional-order-chaotic-systems

[27] W. Deng, C. Li, J. Lü, "Stability analysis of linear fractional differential system with multiple time delays," Nonlinear Dynamics , vol. 48, no. 4, pp. 409-416, 2007.

[28] M. S. Tavazoei, M. Haeri, "A necessary condition for double scroll attractor existence in fractional-order systems," Physics Letters, Section A , vol. 367, no. 1-2, pp. 102-113, 2007.

[]