Chaos in a Fractional-Order Single-Machine Infinite-Bus Power System and Its Adaptive Backstepping Control

This paper has numerically studied the dynamical behaviors of a fractional-order single-machine infinite-bus (FOSMIB) power system. Periodic motions, perioddoubling bifurcations and chaotic attractors are observed in the FOSMIB power system. The existence of chaotic behavior is affirmed by the positive largest Lyapunov exponent (LLE). Based on the fractional-order backstepping method, an adaptive controller is proposed to suppress chaos in the FOSMIB power system. Numerical simulation results demonstrate the validity of the proposed controller.


Introduction
As a mathematical branch with a history of over 300 years, fractional calculus and its applications to physics and engineering have attracted increasing attentions in recent years [1] [2].Fractional calculus provides a good instrument to describe the memory, hereditary, non-locality and self-similarity properties of various materials and processes.Many chaotic systems, such as Lorenz system [3], Chua's system [4], Duffing system [5], Rössler system [6], Chen system [7] and so on, still remain chaotic when their equations become fractional.
Chaotic phenomena have been observed in power systems during the past few decades [8]- [13].Chaos causes electromechanical oscillations to behave randomly, which are harmful to the secure and stable operation of power systems, and even produce undesired negative consequences, such as angle divergence, voltage collapse and system splitting [14].So far, almost all the studies of dynamics of power systems are concerned Z. H. Liang, J. F. Gao 123 with the integer-order models, and there are little research results on fractional modeling and control design of power systems.Tan et al. studied the dynamics of a fractional-order interconnected power system and found that the system became chaotic when the fractional order is no less than 0.88 [15].Sun and Li investigated the chaotic and bifurcation phenomena in a fractional-order three-bus power system and the existence of chaos was demonstrated for different orders [16].
In this paper, we numerically investigate the chaotic dynamics of a fractional-order single-machine infinite-bus (FOSMIB) power system.Period-doubling bifurcation and chaos are observed in FOSMIB power system and the existence of chaos is confirmed by evaluating the largest Lyapunov exponent (LLE).Based on the fractional-order backstepping method, an adaptive controller is presented to suppress chaos in the FOSMIB power system, and the effectiveness of the proposed controller is proved by the numerical simulation results.
The rest of the paper is organized as follows.Some definitions and lemmas about fractional calculus are introduced in Section 2. The dynamics of the FOSMIB power system are analyzed in Section 3.An adaptive controller is designed using the fractional-order backstepping method to suppress chaos in the FOSMIB power system in Section 4. Finally, conclusions are addressed in Section 5.

Preliminaries
There are several different definitions of fractional derivatives.The most appropriate one for practical problems is the Caputo definition.The Caputo fractional derivative is given by where m is integer and ( ) Γ ⋅ is the Gamma function.
The Caputo fractional derivative satisfies the following properties: where C, 1 k and 2 k are real constants.
[17]- [19] Consider the fractional-order system where ( ) ∈ .The equilibrium point * x of system (3) is locally asymptotically stable if all the eigenvalues λ of the Jacobian matrix Lemma 2. [20] Let ( ) x t ∈  be a continuous differentiable function.Then, at any instant the following inequality holds A continuous function Lemma 3. (Fractional-order extension of Lyapunov direct method [22]) Let be an equilibrium point of the nonautonomous fractional-order system ( ) ( ) ( ) with initial condition ( ) are class-K functions.Then , , where ( ) and ⋅ denotes an arbitrary norm.

The FOSMIB Power System
In [12] Chen et al. analyzed the angle dynamics of the classical single-machine infinite-bus (SMIB) power system, which is governed by the so-called swing equation where M is the moment of inertia, D is the damping constant, max P is the maximum power of generator and sin , the SMIB power system is chaotic.Here, we consider the fractional-order single-machine infinite-bus (FOSMIB) power system , sin sin , where 0 1 q < ≤ is the fractional order.When 1 q = , system (11) is the original integer-order SMIB power system.
The autonomous system (11) (as 0 f = ) has two equilibrium points: ( ) . For the equilibrium point O, the Jacobian matrix is and its eigenvalues are In both cases, ( ) . According to Lemma 1, O is asymptotically stable.For the equilibrium point E, the Jacobian matrix is and its eigenvalues are It can be seen that 1 0 λ > and ( ) . In accordance with Lemma 1, E is unstable.

Dynamic Analysis of the FOSMIB Power System
In this section, we use the Adams-Bashforth-Moulton predictor-corrector algorithm proposed by Diethelm et al. in [22]- [24] to solve the FOSMIB power system (11).The , and vary f from 2.4 to 3.5.The corresponding bifurcation diagram is plotted in Figure 1(a), from which a period-doubling route to chaos can be found.To confirm chaos, the largest Lyapunov exponent (LLE) is calculated using Wolf algorithm [25] and plotted in Figure 1 , respectively.After a cascade of period-doubling bifurcations, the system loses its stability and enters chaos at 2.65 f = .As f increases further, the system becomes stable again via inverse period-doubling bifurcations.Now, let and vary q from 0.87 to 1.The resulting bifurcation diagram is plotted in Figure 3(a), which indicates period-doubling bifurcations and chaos.The fractional-order SMIB power system is chaotic over most of the range , where the LLEs are positive as shown in Figure 3(b).The phase portraits for different values of q are plotted in Figure 4.With the increase of q from 0.87, period-1, period-2 and period-4 orbits are obtained at 0.88 q = , 0.893 q = and 0.913 q = , respectively.As q increases further, after a cascade of period-doubling bifurcations, a chaotic attractor is obtained at 0.92 q = .

Adaptive Backstepping Control of Chaos
In this section, an active controller is designed using fractional-order backstepping method to suppress chaos in the FOSMIB power system and stabilize it to the unstable equilibrium point ( ) π,0 E .

Controller Design
Consider the controlled FOSMIB power system   where ( ) and the parameter f is unknown.The backstepping design procedure consists of two steps.
Step 1. Define 1 π e x = − .Its derivative is given by where 2 1 e y α = − , 1 α is the virtual control to be defined later.
Select the candidate Lyapunov function as Now, applying Lemma 2, it can be found that ( ) Define the virtual control 1 α as where 1 c is a positive constant, which leads to 2 1 1 e y c e = + .Substituting Equation (20) into Equation (17) and inequality (19), we have Step 2. The derivative of 2 e is expressed as ( ) where f is the estimate of f .Choose the candidate Lyapunov function as ( ) where k is a positive constant, which can adjust the speed of the adaptive law.Using Lemma 2, it can be found that ( ) sin sin 1 ˆŝin .
Choose the control input and the adaptive law as ( ) ( ) According to Lemma 3, the closed-loop error system is asymptotically stable at the origin ( ) 0, 0 .It means that, with the proposed controller and adaptive law, the FOSMIB power system is asymptotically stable at the equilibrium point ( ) 2 π, 0 E .

Simulation Results
In the simulation, the fractional order q is equal to 0.95.The parameters of system (16) are taken as The closed-loop system consisted of Equations (( 16), ( 26) and ( 27)) is solved by using the predictor-corrector algorithm.The simulation results are shown in Figure 5.
The time-domain waveforms the states of the controlled system ( 16) are shown in Figure 5(a) and Figure 5(b).The FOSMIB power system has experienced chaotic behavior before the controller is put into effect.By activating the controller u at 20 s t = , the chaotic behavior is suppressed and the controlled system converges to the equilibrium point ( ) 2 π, 0 E quickly.The parameter estimate f is converged to f as shown in Figure 5(c) and the controller u is bounded as shown in Figure 5(d).From Figure 5, it can be seen that the proposed controller is feasible for suppressing chaos in the FOSMIB power system.

Conclusion
In this paper, we have numerically investigated the FOSMIB power system.The parameter f and the fractional order q are selected as bifurcation parameters respectively.
Complex dynamical behaviors, such as periodic orbits, period-doubling bifurcations and chaotic attractors, are observed in the FOSMIB power system.The LLE is calculated using Wolf algorithm to confirm the existence of chaos.Furthermore, by exploiting the fractional-order backstepping method, we propose an adaptive controller to suppress chaos in the FOSMIB power system.The effectiveness of the presented controller is verified by numerical simulation results.
of the machine.Let x θ = and y θ =  , then Equation (9) can be rewritten as dynamics are numerically analyzed by means of bifurcation diagrams, phase portraits and Lyapunov exponents.In the following simulations, parameter f is chosen as bifurcation parameter and the other parameters are fixed at 0
f = .The parameters of the controller (26) and the adaptive law (27) are chosen as

Figure 5 .
Figure 5.The time-domain waveforms of the controlled system (16).