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 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.

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, and are real constants.

Lemma 1. [17] - [19] Consider the fractional-order system

where and. The equilibrium point of system (3) is locally asymptotically stable if all the eigenvalues of the Jacobian matrix satisfy

Lemma 2. [20] Let be a continuous differentiable function. Then, at any instant the following inequality holds

A continuous function is referred as class-K if it is strictly increasing and [21] .

Lemma 3. (Fractional-order extension of Lyapunov direct method [22] ) Let be an equilibrium point of the nonautonomous fractional-order system

with initial condition. Assume that is a Lyapunov candidate and are class-K functions. Then is asymptotically stable if the following conditions hold

where and denotes an arbitrary norm.

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, is the maximum power of generator and is the power of the machine.

Let and, then Equation (9) can be rewritten as

where, and are positive constant parameters. When, , and, the SMIB power system is chaotic.

Here, we consider the fractional-order single-machine infinite-bus (FOSMIB) power system

where is the fractional order. When, system (11) is the original integer-order SMIB power system.

The autonomous system (11) (as) has two equilibrium points: and. 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 and. In accordance with Lemma 1, E is unstable.

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 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, ,. The initial conditions are selected as,.

First, let, 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(b). The FOSMIB power system is chaotic over most of the range, where the LLEs are positive. The phase portraits for different values of f are plotted in Figure 2. With the increase of f from 2.4, period-1, period-2 and period-4 orbits are obtained at, and, respectively. After a cascade of period-doubling bifurcations, the system loses its stability and enters chaos at. 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, and, respectively. As q increases further, after a cascade of period-doubling bifurcations, a chaotic attractor is obtained at.

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.

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.

The work was supported by the Natural Science Foundation of Henan Province, China (Grant No. 14A120005) and Excellent Young Scientist Development Foundation of Zhengzhou University, China (Grant No. 1421319086).

Liang, Z.H. and Gao, J.F. (2016) Chaos in a Fractional-Or- der Single-Machine Infinite-Bus Power System and Its Adaptive Backstepping Control. International Journal of Modern Nonlinear Theory and Application, 5, 122-131. http://dx.doi.org/10.4236/ijmnta.2016.53013