Complexity of a Duopoly Game in the Electricity Market with Delayed Bounded Rationality

According to a triopoly game model in the electricity market with bounded rational players, a new Cournot duopoly game model with delayed bounded rationality is established. The model is closer to the reality of the electricity market and worth spreading in oligopoly. By using the theory of bifurcations of dynamical systems, local stable region of Nash equilibrium point is obtained. Its complex dynamics is demonstrated by means of the largest Lyapunov exponent, bifurcation diagrams, phase portraits, and fractal dimensions. Since the output adjustment speed parameters are varied, the stability of Nash equilibrium gives rise to complex dynamics such as cycles of higher order and chaos. Furthermore, by using the straight-line stabilizationmethod, the chaos can be eliminated. This paper has an important theoretical and practical significance to the electricity market under the background of developing new energy.


Introduction
In 1980s, chaos theory was first introduced into the economic research.Chaotic economists used the basic mathematic theory of chaos to improve the existing models of economic phenomena.The economic system is whether a chaotic system is a hot topic in the economic field.Bifurcation theory based on difference equation has been applied in all branches of chaos 1 .In recent years, a series of dynamic game models on the output decision Cournot model and price decision Bertrand model have been studied in related references.Agiza 2 and Kopel 3 have considered bounded rationality and established duopoly Cournot model with linear cost functions.From then on, the model has been extended to multioligopolistic market.Bischi et al. 4 suppose that firms determine their output based on the reaction functions, that is, all the players take adaptive expectation.Agiza and Elsadany 5 have improved the model that contains two-types of heterogeneous players: boundedly rational player and adaptive expectation player.Zhang et al. 6 have further improved the model with nonlinear cost functions.Matsumoto and Nonaka 7 have researched the complexity of Cournot model with linear cost functions.Ma and Ji 8 have constructed and considered a Cournot model in electric power triopoly with nonlinear inverse demand, and the model is further studied by Ji 9 based on heterogeneous players.Ma and Feng 10 have studied the chaotic behavior in retailer's demand model.Xin et al. 11 have researched the complex dynamics of an adnascent-type game model.Chen et al. 12 have used Bertrand model with linear demand functions to study the competition in Chinese telecommunications market.Sun and Ma 13 have introduced a Bertrand model with nonlinear demand functions in Chinese cold rolled steel market and researched the complexity and the control of the model.Yassen and Agiza 14 have considered a Cournot duopoly game and the model with delayed rationality.In these pieces of literatures, adjustment speed or other parameters are taken as bifurcation parameters, and complex results such as period doubling bifurcation, unstable period orbits, and chaos are found.
Economic dynamics seem to devote new interest to delay differential equations.This is because some economic phenomena cannot be described exhaustively with pure linear or nonlinear differential equations.Differential equations with time delay play an important role in economy, engineering, biology, and social sciences, because a great deal of problems may be described with their help.Based on the game model 8 , a new duopoly game model with delayed bounded rationality in the electricity market is obtained.The duopoly model with delayed bounded rationality is closer to the economic reality and is worth being used in oligopoly.Suppose the inverse demand function is nonlinear, and cost functions are one nonlinear and one linear.In this model, the bounded rational players regulate output speed according to marginal profit and decide the output.By theoretical analysis and numerical simulation, the stable region about the output adjustment speed parameters is derived.It is shown that the output adjustment speed lead to the chaos at a definite range.It has an important theoretical and applied significance to research the complexity of new style nonlinear dynamical system.
The paper is organized as follows.In Section 2, the dynamics of a game with delayed bounded rationality is presented.In Section 3, the existence, local stability, and bifurcation of the equilibrium points are also analyzed.Numerical simulations are used to show the complex characteristics of the system via computation of Lyapunov exponents, confirmation of the system sensitive dependence on initial conditions and calculation of the fractal dimension of the chaotic attractor.In Section 4, bifurcation and chaos control of the model is considered with the straight-line stabilization method.Finally, some conclusions are made.

Model
Suppose that there are two representative electricity enterprises the enterprise X represents the traditional coal electricity enterprise, the enterprise Y represents new energy enterprise such as water electricity, nuclear power, or solar energy, or wind energy, etc. in the electricity market, and they provide electric power to consumers through the electric utilities.The electricity enterprises X, Y make the optimal output decision and suppose the t-output is q i t , i 1, 2 , respectively.
At each period t, the price p is determined by the total output Q t q 1 t q 2 t .According to 8 , the inverse demand function is

2.1
We propose the cost function of the enterprise X is nonlinear.We propose the cost function of the enterprise Y is linear because the new energy enterprise has lower variable cost.So, the cost functions of the two enterprises are as follows:

2.2
It is assumed that the wheeling rate is γ.The profits of the companies are

2.3
As the game between the enterprises is a continuous and long-term repeated dynamic process in the electricity market, so the dynamic adjustment of this repeated Cournot duopoly game with bounded rational players is as follows: where α i , i 1, 2 is output adjustment speed parameters.Combining 2.3 , 2.4 , a dynamic duopoly game with bounded rationality has the following form:

2.5
According to 14 , there are two reasons for the occurrence of delayed structure in economic models: one is that decisions made by economic agents at time t is depended on the past observed variables by means of a prediction feedback, and the other is that the functional relationships describing the dynamics of the model both depend on the current state of the economy and in a nontrivial manner on past states.
Due to incomplete information and delayed decision, we propose that there is one step T 1 delay in the output of the mutual enterprises.Therefore, the dynamic game model 2.5 with delayed is as follows:

Complex Dynamics Character of System 2.6
The bounded rational player makes output decision on the basis of the marginal profit of the last period.The company decides to increase output if it has a positive marginal profit and decrease output if the marginal profit is negative.Thus, output adjustment speed parameter α i , i 1, 2 has an important effect on game results.In the following section, the effect of α i , i 1, 2 on dynamical behaviors of system 2.6 will be taken into account.

The Equilibrium Point and Stability Analysis
The bifurcation parameters are α 1 , α 2 , and the other parameters of system 2.6 are as follows: m 5.2, n 0.95, b 1 0.42, c 1 0.25, b 2 0.37, and γ 0.15.
By solving the following equations, the fixed points of system 2.6 can be obtained
To study the stability of system 2.6 , we rewrite it as a fourth-dimensional system in the form

3.2
The Jacobian matrix of 3.2 at the Nash equilibrium point p * is where 3.4 The characteristic polynomial of 3.2 is where A − j 33 j 44 2 , B 1 j 33 1 j 44 , C 0, D −j 32 j 41 .

3.6
Denote and ν 3 μ 2 μ 4 − μ 1 μ 3 ; the necessary and sufficient conditions for the local stability of Nash equilibrium can be gained by Jury test 15 :

3.7
Through computing the above equations, the local stable region of Nash equilibrium point can be obtained.The stable region of slash with positive α 1 , α 2 is shown in Figure 1.The Nash equilibrium is stable for the values α 1 , α 2 inside the stable region.The meaning of the stable region is that whatever initial outputs are chosen by the two electricity companies in the local stable region, they will eventually arrive at Nash equilibrium output after finite games.It is valuable to analyze the enterprises on accelerating the output adjustment speed to increase their profits.While output adjustment parameters do not change Nash equilibrium point, but once one party is adjusting output speed too fast and pushing α i , i 1, 2 out of the stable region, the system tends to become unstable and fall into chaos.Numerical simulation is used to analyze the characteristics of nonlinear dynamical system with the change of α i , i 1, 2 .Numerical results such as bifurcation diagrams, the largest Lyapunov exponent, strange attractors, sensitive dependence on initial conditions, and fractal structure will be researched.

The Output Adjustment Speed Effect on the System
Once company X accelerates output adjustment speed and pushes α 1 out of the stable region, the stability of Nash equilibrium point will change.Figure 2 demonstrates that the output evolution of the duopoly starts with equilibrium state, undergoes period doubling bifurcation to chaotic state, occurs period doubling bifurcation again and ends with chaotic state when the output adjustment speed α 1 increases and α 2 0.2.
When α 2 0.2, the diagrams of bifurcation and the largest Lyapunov exponent with α 1 increasing are shown in Figure 2.For 0 < α 1 < 0.2362, the largest Lyapunov exponent is negative and the system 2.6 is stable at Nash equilibrium point.For 0.2362 < α 1 < 0.2790, the largest Lyapunov exponent is negative and system 2.6 has an orbit of 2-cycle.For 0.2790 < α 1 < 0.3059 and 0.3329 < α 1 < 0.3427, the largest Lyapunov exponent is negative and system 2.6 has an orbit of 4-cycle.For 0.3427 < α 1 < 0.3476, the largest Lyapunov exponent is negative and system 2.6 has an orbit of 8-cycle.For 0.3059 < α 1 < 0.3280, the largest Lyapunov exponent just a little bigger than zero and system 2.6 has a double chaotic attractor as shown in Figures 3 a and 3 b .For 0.3280 < α 1 < 0.3329 and 0.3476 < α 1 < 0.35, the largest Lyapunov exponent is obviously positive and system 2.6 is in a chaotic state as shown in Figures 3 c and 3 d .
Similarly, Figure 4 shows a one-parameter bifurcation diagram and the largest Lyapunov exponent with respect to α 2 and when α 1 0.15.We can see that Nash equilibrium point is stable for 0 < α 2 < 0.2621, which implies that output of two firms is in equilibrium state.With α 2 increasing, the stability of equilibrium point changes, and output undergoes   Figure 5 illustrates chaos attractor with initial point q 0 1 0.6, q 0 2 0.5 when α 2 0.362, α 1 0.15.
The sensitive dependence on initial conditions is one of the important features of chaos.To verify whether system 2.6 depends on initial values sensitively, the relationships between output and time are shown in Figures 6 a , 6 b , 6 c , and 6 d when α 1 0.348, α 2 0.2 and α 1 0.15, α 2 0.362, respectively.At first, the difference is indistinguishable, but with the number of the game increasing, the difference between them is huge.This demonstrates that only a little difference between initial data will have a great impact on the results of the game.It further proves that system 2.6 falls into a chaotic state when α 1 0.348, α 2 0.2 and α 1 0.15, α 2 0.362.While the system is in a chaotic state, the market becomes volatile and it is difficult for electricity companies to plan long-term strategy.A slight adjustment of the initial output can have a great effect on the game results.
A fractal dimension is taken as a criterion to judge whether the system is chaotic.There are many specific definitions of fractal dimension and none of them should be treated as the universal one.According to 16 , we adopt the following definition of fractal dimension: where λ 1 > λ 2 >, . . ., λ n are the Lyapunov exponents and j is the maximum integer for which j 1 λ i > 0 and The Lyapunov exponents of system 2.6 are λ 1 0.067051, λ 2 −0.106206 when α 1 0.15, α 2 0.362.The largest Lyapunov exponent λ 1 is positive, which indicates system (t) q 2 0 = 0.5 q 2 0 = 0.51 d Figure 6: Sensitive dependence on initial conditions: the two orbits of a q 1 -coordinates for initial points are 0.62, 0.52 and 0.63, 0.52 ; b q 2 -coordinates for initial points are 0.62, 0.52 and 0.62, 0.53 ; c q 1coordinates for initial points are 0.6, 0.5 and 0.61, 0.5 ; d q 2 -coordinates for initial points are 0.6, 0.5 and 0.6, 0.51 .
2.6 is in a chaotic state.Fractal dimension illustrates that the chaotic motion has self-similar structure, that is to say, the chaotic motion follows a definite rule.This is an important difference between chaotic motion and stochastic motion.The fractal dimension of system 2.6 is d 2 − λ 1 /λ 2 ≈ 1.6313.The fractal dimension reflects the space density of the chaotic attractor 17 .The larger the dimension of the chaotic attractor is, the larger is the occupied space.Consequently, the structure is more compact, and the system is more complicated and vice versa.The fractal dimension of the 2D discrete system 2.6 is more than 1.5, so the occupied space is big and the structure is tight, which can be seen in Figure 5.

Chaos Control
The numerical simulation results show that the oligopolistic market becomes unstable and even falls into chaos when output adjustment speed parameter goes beyond the stable region.
All the players will be harmed and the market will become irregular when chaos occurs.So, nobody is able to make good strategies and decide reasonable output.To avert the risk, it is a good ideal for duopoly to maintain at Nash equilibrium output.Perturbation feedback and nonfeedback are two methods for the chaos control.Recently, Yang et al. 18 and Xu et al. 19 proposed a new control method, which is called the straight-line stabilization method.This method is adopted to control the chaos in this paper where |k| < 1 is the feedback control parameter and other parameters are the same as above.
Adding the external control signal 4.1 to system 2.6 , the controlled system is as follows:

4.2
It can be seen from Figure 7, at α 1 0.15, α 2 0.362 , that controlled system 4.2 stabilized at Nash equilibrium point when −0.9876 < k < −0.28.It demonstrates that the stable control of the specific goal can be realized even if the perturbation is very small.At k −0.5, α 1 0.15, the stable range of system 4.2 is 0.289 < α 2 < 0.436.The range of 2-cycle bifurcation is 0.436 < α 2 < 0.5171.The range of 4-cycle bifurcation is 0.5171 < α 2 < 0.5374.For 0.5374 < α 2 < 0.6084, system 4.2 is in the chaotic state.The details are shown in Figure 8, which indicates that once system 2.6 is under a chaotic state, system 4.2 can be efficient to eliminate chaos.

Conclusions
A new dynamics of nonlinear duopoly game in the electricity market with delayed bounded rationality is established in this paper.The stability of equilibria, bifurcation, and chaotic behavior of the duopoly game are investigated.It is found that bifurcation, chaos, and other complex phenomena occur when the output adjustment speed parameter changes.It is well known that the occurrence of chaos depends on the values of bifurcation parameters.The straight-line stabilization method is used to control the period-doubling bifurcation, unstable periodic orbits, and chaos.The system quickly arrived at the Nash equilibrium point when a small perturbation is applied in the chaos region.The research results have an important theoretical and practical significance to the electricity market under the view of developing new energy.This paper also shows guidance for electricity companies to formulate strategies of output and is helpful for the government to formulate relevant policies to macrocontrol economy.

Figure 1 :
The stable region of Nash equilibrium point about adjustment speed α 1 , α 2 .