#### Abstract

Now dealing with power system became the most important arts. In this paper we report the dynamic analysis of a fractional-order power system with parameter . To the best know of our knowledge, that was the first report about bifurcation analysis of the fractional order power system. So first we discuss the dynamic analysis with different fractional order and different parameters. Furthermore we will establish its numerical simulations which are provided to demonstrate the feasibility and efficacy of our analysis.

#### 1. Introduction

For a long time, many researchers have been growing interest in investigating the potential use of dynamics in applications, taking power system as an example. Under normal and emergency operating conditions, power systems always can show rich nonlinear dynamic phenomena. From a mathematical point of view, differential equations can be used to describe the dynamic behavior of the power system. The power system is usually run under equilibrium conditions, but in engineering practice, the stability of the equilibrium point form sometimes changes with the range of small perturbations. After 1980s, with the development of the power system, even if a small perturbation can also cause the system to lose stability, due to the changing of some control parameters. Researchers focus on the study of the critical state of the system.

Along with the 300-year-old history of the fractional calculus, its applications to electrical engineering have received more attention. Typically, chaotic systems are still chaotic even when their differential equations become fractional-order differential systems [1–9]. We all know there are many results on fractional-order system, but the bifurcation of fractional-order nonlinear systems has not been well studied. Obviously there are two kinds of methods which have been used in the previous paper to solve fractional-order differential systems: the frequency-domain methods [10] and the time-domain methods [11]. Some researcher found it is sometimes invalid to research some chaotic systems who used the frequently method to investigate chaos [12, 13]. However, we give up the first method and choice the time domain method to make the numerical simulations in this paper. The main idea was introduced by Diethelm et al. [14] and has been used by the following paper [15–17].

In this paper, we report the first investigation of a fractional-order power system using the time domain method. Recently, a simplified 4D power system was reported by Chiang et al. [18]. Its chaotic and bifurcation analysis have been investigated and got well results [19, 20]. But to the best know of our knowledge, that was the first report about bifurcation analysis of the fractional-order system. According to the different values we will show the bifurcation diagrams of the system which chaos exists.

The paper is organized as follows. In Section 2, the numerical algorithm for the fractional-order power system is presented. In Section 3, the chaotic behavior and bifurcations of the system are studied. Finally, we summarize the results.

#### 2. Numerical Algorithm for Fractional-Order Power System

##### 2.1. System Description [21]

In this paper we consider the following power system in [21] (see Figure 1). Here we use as the varying parameter to derive the corresponding attractor.

**(a)**

**(b)**

Now we consider the fraction order power system by here are variables and (), its order is denoted by . When , , system (1) becomes the integer order power system in [21].

##### 2.2. Fractional Derivative [22]

There are several definitions of a fractional-order differential system. In the following, we introduce the most common one of them [22]: with ; that is, is the first integer which is not less than , is the -order derivative in the usual sense, and () is the -order Reimann-Liouville integral operator with expression Here stands for Gamma function, and the operator is generally called “-order Caputo differential operator” [22].

Next, we consider the fractional-order differential equation with related initial conditions

Combined with the methods [11, 14] and the above definition of fractional-order differential equation. Set , the fractional-order power system can be rewritten in which

#### 3. Bifurcations Analysis of the Fractional-Order Power System

##### 3.1. Integer-Order Power System

In the following discussion, we mainly use the simulation method to analyze the fractional-order power system by calculating the largest Lyapunov exponent using the Wolf algorithm [6]. Because the nonperiodicity solutions and sensitive dependence on initial conditions are important reason of chaos, here, we let the initial condition be defined as [21] and vary (). According to our numerical simulations, we have visually founded the difference between the integer order and the fractional order. Here we let (), that is, integer order, the 3D trajectory in the power system is shown in Figure 1. The related bifurcation diagram of the parameter is shown in Figure 1, and the largest Lyapunov exponent diagram is shown in Figure 2.

##### 3.2. Dynamic Analysis of the Fractional-Order System

Next we correspond to discuss the fractional-order power system. The fractional orders () are equal and fixed with at some value which is from 0 to 1. On the basis of Section 2, first we find that the range of parameter is changing along with the changing of fractional order. When we let () be an example while the parameter is changed from 0 to 11.4 [21]. For a step size in is 0.01, we clearly got the range of in which the solutions of the fractional-order power system are bounded. We take as an example, the 2D trajectory is shown in Figure 3 which clearly finds the difference between two orders 0.99 and 1. The difference is that the former generates chaotic phenomena and the later happens that after a period of offloading shock of system voltage, ultimately stabilize at a certain value, corresponding to the equilibrium point is asymptotically stable focus. Further, we let (); we clearly got the range of . The 2D trajectory is shown in Figure 4 in which and .

**(a)**

**(b)**

##### 3.3. Dynamic Analysis with Different Order

(1)Fix , , and let vary. We calculate numerically the fractional-power system for with a varied step size equal to 0.01. The related Largest Lyapunov exponent is shown in Figure 5.(2)Fix , , and let vary. We calculate numerically the fractional-power system for with a varied step size equal to 0.01. The related Largest Lyapunov exponent is shown in Figure 5.(3)Fix , , and let vary. We calculate numerically the fractional-power system for with a varied step size equal to 0.01. The related Largest Lyapunov exponent is shown in Figure 6.(4)Fix , , and let vary. We calculate numerically the fractional-power system for with a varied step size equal to 0.01. The related Largest Lyapunov exponent is shown in Figure 6.

**(a)**

**(b)**

**(a)**

**(b)**

#### 4. Conclusions

Through the above discussion, we successfully have numerically studied the existence of the bifurcations type of the fractional order power system with the difference orders. According to the different initial conditions and fractional order, the largest Lyapunov exponent values and the bifurcation maps are given to verify the rationality of the analysis. Of course, the next work on this topic should include the depth study in chaos control and synchronization.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgment

The present work was partially supported by the National Natural Science Foundation of China under Grant no. 61273128.