Next Article in Journal
Entropy-Based Modeling of Velocity Lag in Sediment-Laden Open Channel Turbulent Flow
Next Article in Special Issue
Computational Complexity
Previous Article in Journal
Exergy Analysis of a Syngas-Fueled Combined Cycle with Chemical-Looping Combustion and CO2 Sequestration
Previous Article in Special Issue
Normalized Minimum Error Entropy Algorithm with Recursive Power Estimation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Entropy Complexity and Stability of a Nonlinear Dynamic Game Model with Two Delays

1
College of Management and Economics, Tianjin University, Tianjin 300072, China
2
School of Management Science and Engineering, Anhui University of Finance and Economics, Bengbu 233030, China
*
Authors to whom correspondence should be addressed.
Submission received: 13 June 2016 / Revised: 17 August 2016 / Accepted: 17 August 2016 / Published: 30 August 2016
(This article belongs to the Special Issue Computational Complexity)

Abstract

:
In this paper, a duopoly game model with double delays in hydropower market is established, and the research focus on the influence of time delay parameter on the complexity of the system. Firstly, we established a game model for the enterprises considering both the current and the historical output when making decisions. Secondly, the existence and stability of Hopf bifurcation are analyzed, and the conditions and main conclusions of Hopf bifurcation are given. Thirdly, numerical simulation and analysis are carried out to verify the conclusions of the theoretical analysis. The effect of delay parameter on the stability of the system is simulated by a bifurcation diagram, the Lyapunov exponent, and an entropic diagram; in addition, the stability region of the system is given by a 2D parameter bifurcation diagram and a 3D parameter bifurcation diagram. Finally, the method of delayed feedback control is used to control the chaotic system. The research results can provide a guideline for enterprise decision-making.

1. Introduction

As the demand for power increases, the supply is proportionately increased. However, the cost of thermal power, which is generated by non-renewable resources, is extremely high. Therefore, the development of hydroelectric power is a more sensible choice. In this paper, the game process of power output of two hydropower enterprises is considered. In practice, there is an obvious relationship between the amount of electricity generated and the amount of water reserves. As a result, the water reserves, the current output, and the output of the time delay are mainly considered to make an electricity output plan. Therefore, it is very meaningful to study the game competition of the electricity output of the two-oligopoly enterprises.
Electricity is an energy source of great significance, related to country development and societal stability, so has always been the focus of research. As for studies on energy in relation to economic strategies and the environment, Weidou [1] analyzed China’s energy policies and provided reasonable suggestions for the supply and consumption of China’s energy, which would ensure the sustainable development of China’s social energy economy. Nwaobi [2] proposed an economic model based on the emission reduction policy, and gave an empirical analysis of Nigeria as an example. Vera et al. [3] are in the process of analyzing the problem of 3E. They established a national energy index system to develop the energy policy and realize the sustainable development of society, the economy, the environment, and energy. Omri [4] analyzed the relationship between economic growth and a group of parameters: energy consumption, electricity consumption, nuclear consumption, and renewable energy consumption.
In recent years, the oligopoly game has always been the frontier research area of renewable energy, especially the nonlinear dynamic model, which is closely related to real life. Using the entropy theory and chaos theory to establish an oligopolistic competition game model, establishing a measuring method for uncertainty, and eliminating uncertainty all have broad applicability.
Dajka et al. [5] considered a two-player quantum game in the presence of a thermal decoherence modeled in terms of a rigorous Davies approach. This shows how the energy dissipation and pure decoherence affect the payoffs of the players of the game. Harré et al. [6] considered the issue about how changes in the underlying incentives can move us from an optimal economy to a sub-optimal economy, meanwhile making it impossible to collectively navigate our way to a better strategy without forcing us to pass through a socially undesirable “tipping point”.
The power market of the three-oligopoly game was analyzed by Ma [7], where the complex dynamic characteristics of the system are studied, and the dynamic behavior of the game is given, providing an excellent practical guideline that is of great significance. Batabyal [8] took the international trade of renewable resources as the background, did corresponding analysis using the Stackelberg differential model, and pointed out that the policy tool has an important role in promoting energy conservation. This study has important implications for the conservation of renewable energy. As for the Bertrand model, Sun et al. [9] proposed a three-oligopoly game model, which is based on the cold rolling market of China and studied the complex dynamic characteristics of the game process. Basing on the Markov information structure, Halkos et al. [10] put forward a renewable energy and non-renewable energy Nash game, and carried out the analysis of its strategy by using the utility function, which is of great significance in practice. Liao [11] analyzed the role of developing hydro-energy, wind energy, nuclear power, and so forth. Yoon and Ratti [12] examined the effect of energy price uncertainty on firm-level investment with data on U.S. manufacturing firms.
Chaos analysis and applications in dynamical systems are observed in many practical applications in engineering, biology, and economics [13,14,15]. Sun and Tian established an energy resource demand-supply system based on the background of the real energy resources demand and supply in the East and the West of China [16], and have obtained a series of findings [17,18].
In this paper, we study the problem of the oligopoly game in the market of hydroelectric power using nonlinear theory and complex dynamics theory. We analyze the complex dynamic characteristics of the system and study the influence of the time delay and weight on the system. Moreover, an effective method for controlling a chaotic system is carried out.
This paper is organized as follows. In Section 2, a continuous differential duopoly game mode with two delays is established. In Section 3, we focus on analyzing the existence and stability of Hopf bifurcation. In Section 4, we carry out numerical simulation and analysis by using the methods of the attractor, bifurcation diagram, Lyapunov exponent, entropic, etc. to study the influence of delay and weight on the stability of the system. In Section 5, we elaborate on delayed feedback control, which is an effective control method of chaos. Finally, the conclusions of this paper are given in Section 6.

2. The Model

In the market of hydroelectric power, we want to find the total economic variables: electrovalence, market demand, and water reserves. The conditions need to be met when the market reaches equilibrium. So, in this paper, we focus on the influence of the time delay parameter on the dynamic characteristics of the system.
Assumptions are listed as follows:
  • The electricity market is composed of two hydropower enterprises.
  • The output of electricity quantity is only related to the amount of water reserves and electrovalence.
  • Electricity price is determined by the output of electricity quantity of two enterprises and social demand.
  • In order to generate electricity and ensure water storage safety, the minimum water reserve is X min the maximum water reserve is X max . Therefore, the water reserves X must satisfy X ( X min , X max ) .
The inverse demand function for the electricity market is p ( q ) = p ¯ q , where p is electrovalence, q = i = 1 N q i is the demand of electrical energy, and p ¯ is the highest price of the electrical energy [9].
In the electricity market, the cost function C ( X ) is related to its own water reserves X [10,11]. X min is the water reserve when C ( X min ) = p ¯ . Obviously, x = X X min is the effective water reserve. When X X min , hydro power stations will suspend power generation. The unit cost of water reserves is c ( x ) = p ¯ c x . If water is renewable, and the growth rate of the water reserve is k , then x ˙ = k x . So, the water reserves should meet x ˙ = k x q .
In the process of the market game, every player in the game has to achieve its own profit maximization:
π = [ p ( q ) c ( x i ) ] q i = ( c x i q ) q i .
So, the problem of the game strategy for each oligopoly can be expressed as the optimization problem:
max q i 0 e θ t ( c x i q ) q i d t
s . t . x ˙ i = k x i q i ,
where θ is discount factor, θ > 0 , q i 0 , x i 0 and i = 1 , 2 .
Solving the optimization problem (1) with its Hamiltonian function:
H i = ( c x i q 2 ) q i + λ i ( k x i q i ) ,
where λ i is a co-state variable. Solving:
H i q i = 0 ,   λ ˙ i = H i x i + θ λ i .
Simultaneously we can obtain:
q ˙ i + q ˙ = ( θ k ) ( q ˙ i + q ˙ ) + c ( 2 k θ ) x i ,   i = 1 , 2 .
The sum of the two parts of Equation (4) is obtained, and, using q = i = 1 N q i and x = i = 1 N x i ( N = 2 ) , we can derive:
q ˙ = ( θ k ) q + c ( 2 k θ ) N + 1 x .
With Equations (4) and (5), we can calculate the game strategy for each oligopoly as:
q ˙ i = ( θ k ) q i + c ( 2 k θ ) x i c ( 2 k θ ) 3 x .
In the actual operation of the electricity market, in order to ensure the sustainable development of electric energy, as well as the stability of the electrovalence, the electricity supply of the market should be fixed in a certain period of time, to ensure that the supply of electricity is not too large or too small, which leads to fluctuations in market price. Suppose that in a certain period of time, the largest electricity supply is a , so q ˙ i positive correlation with v i ( a q ) , v i is the change rate of electricity supply.
Therefore, integrating Equations (5) and (6) and v i ( a q ) , we can obtain a double oligopoly game model for the hydroelectric power market:
{ x ˙ 1 = k x 1 q 1 x ˙ 2 = k x 2 q 2 q ˙ 1 = v 1 ( a q 2 q 1 ) [ ( θ k ) q 1 c ( 2 k θ ) x 2 + 2 3 c ( 2 k θ ) x 1 ] q ˙ 2 = v 2 ( a q 1 q 2 ) [ ( θ k ) q 2 c ( 2 k θ ) x 1 + 2 3 c ( 2 k θ ) x 2 ] .
In the actual operation process of the electric power enterprise, it is not only necessary to consider the current market demand; we also need to consider the market demand of τ time early, which makes the final decision results closer to the actual situation. In this paper, we assume that the duopoly considers time delay. The delay parameters are τ 1 and τ 2 respectively. That is:
q i ( t ) = μ i q i ( t ) + ( 1 μ i ) q i ( t τ i ) ,
and 0 < μ i < 1 , i = 1 , 2 is the weight of the current period price.
We know that the time delay does not affect the system’s solution. Therefore, the equilibrium point E * ( 0 , 0 , 0 , 0 ) of Equation (7) is also the equilibrium point with a time delay system. So, in this paper, we analyze the delay decision with the equilibrium point of E * ( 0 , 0 , 0 , 0 ) . Then, from Equations (7) and (8), the linear system with time delay is:
{ x ˙ 1 = k x 1 q 1 , x ˙ 2 = k x 2 q 2 , q ˙ 1 ( t ) = a v 1 [ ( θ k ) ( μ 1 q 1 ( t ) + ( 1 μ 1 ) q 1 ( t τ 1 ) ) c ( 2 k θ ) x 2 ( t ) + 2 3 c ( 2 k θ ) x 1 ( t ) ] , q ˙ 2 ( t ) = v 2 ( a μ 1 q 1 ( t ) ( 1 μ 1 ) q 1 ( t τ 1 ) μ 2 q 2 ( t ) ( 1 μ 2 ) q 2 ( t τ 2 ) )        [ ( θ k ) ( μ 2 q 2 ( t ) + ( 1 μ 2 ) q 2 ( t τ 2 ) ) c ( 2 k θ ) x 1 ( t ) + 2 3 c ( 2 k θ ) x 2 ( t ) ] .

3. Stability and Existence of Hopf Bifurcation

Equation (9) is linearized at the equilibrium point E * ( 0 , 0 , 0 , 0 ) by using s Jacobian matrix. So Equation (9) can be converted to:
{ x ˙ 1 = k x 1 q 1 , x ˙ 2 = k x 2 q 2 , q ˙ 1 = A 1 q 1 ( t ) + B 1 ( 1 μ 1 ) q 1 ( t τ 1 ) + C 1 x 1 ( t ) + D 1 x 2 ( t ) , q ˙ 2 = μ 2 q 2 ( t ) + A 2 ( 1 μ 2 ) q 2 ( t τ 2 ) ( a ν 2 M + M ) x 1 ( t ) + a ν 2 N x 2 ( t )      B 2 μ 2 q 1 ( t ) q 2 ( t ) B 2 ( 1 μ 2 ) q 1 ( t ) q 2 ( t τ 2 ) + ν 2 μ 1 q 1 ( t ) ν 2 μ 1 q 1 ( t ) N x 2 ( t )      C 2 μ 2 q 1 ( t τ 1 ) q 2 ( t ) C 2 ( 1 μ 2 ) q 1 ( t τ 1 ) q 2 ( t τ 2 ) + D 2 M q 1 ( t τ 1 ) x 1 ( t )      D 2 N q 1 ( t τ 1 ) x 2 ( t ) E 2 q 2 ( t ) 2 E 2 ( 1 μ 2 ) q 2 ( t ) q 2 ( t τ 2 ) + ν 2 μ 2 M q 2 ( t ) x 1 ( t )      ν 2 μ 2 N q 2 ( t ) x 2 ( t ) F 2 μ 2 q 2 ( t ) q 2 ( t τ 2 ) F 2 ( 1 μ 2 ) q 2 ( t τ 2 ) 2      + G 2 M q 2 ( t τ 2 ) x 1 ( t ) G 2 N q 2 ( t τ 2 ) x 2 ( t ) ,
where
A 1 = a ν 1 μ 1 ( θ k ) ,   B 1 = a ν 1 ( θ k ) ,   C 1 = 2 3 a c ν 1 ( 2 k θ ) ,
D 1 = a ν 1 c ( 2 k θ ) ,   M = c ( 2 k θ ) ,   N = 2 3 c ( 2 k θ ) ,   A 2 = a ν 2 ( θ k ) ,
B 2 = ν 2 μ 1 ( θ k ) ,   C 2 = ν 2 ( 1 μ 1 ) ( θ k ) ,   D 2 = ν 2 ( 1 μ 1 ) ,
E 2 = ν 2 μ 2 ( θ k ) ,   F 2 = ν 2 ( 1 μ 2 ) ( θ k ) ,   G 2 = ν 2 ( 1 μ 2 ) .
The characteristic equation of Equation (10) is
λ 4 + f 13 λ 3 + f 12 λ 2 + f 11 λ + f 10 + ( f 23 λ 3 + f 22 λ 2 + f 21 λ + f 20 ) e λ τ 1 + ( f 33 λ 3 + f 32 λ 2 + f 31 λ + f 30 ) e λ τ 2 + ( f 42 λ 2 + f 41 λ + f 40 ) e λ ( τ 1 + τ 2 ) = 0 ,
where
f 13 = 2 k μ 2 A 1 ,
f 12 = C 1 + k 2 + 2 A 1 k + A 1 μ 2 + 2 k μ 2 + N a v 2 ,
f 11 = D 1 μ 1 v 2 A 1 k 2 k 2 μ 2 C 1 k C 1 μ 2 2 A 1 k μ 2 N a A 1 v 2 N a k v 2 ,
f 10 = C 1 k μ 2 M D 1 + A 1 k 2 μ 2 C 1 k μ 1 v 2 M a D 1 v 2 + N a C 1 v 2 + N a A 1 k v 2 ,
f 23 = B 1 μ 1 B 1 ,
f 22 = 2 B 1 k 2 B 1 k μ 1 + B 1 μ 2 B 1 μ 1 μ 2 ,
f 21 = B 1 k 2 μ 1 B 1 k 2 N a B 1 ν 2 2 B 1 k μ 2 + N a B 1 μ 1 ν 2 + 2 B 1 k μ 1 μ 2 ,
f 20 = B 1 k 2 μ 2 B 1 k 2 μ 1 μ 2 + N a B 1 k v 2 N a B 1 k μ 1 v 2 ,
f 33 = A 2 μ 2 A 2 ,
f 32 = A 1 A 2 2 A 2 k μ 2 + 2 A 2 k A 1 A 2 μ 2 ,
f 31 = A 2 k 2 + A 2 k 2 μ 2 A 2 C 1 2 A 1 A 2 k + 2 A 1 A 2 k μ 2 + A 2 C 1 μ 2 ,
f 30 = A 1 A 2 k 2 + A 2 C 1 k A 2 C 1 k μ 2 A 1 A 2 k 2 μ 2 ,
f 42 = A 2 B 1 + A 2 B 1 μ 1 μ 2 A 2 B 1 μ 1 A 2 B 1 μ 2 ,
f 41 = 2 A 2 B 1 k μ 2 2 A 2 B 1 k + 2 A 2 B 1 k μ 1 2 A 2 B 1 k μ 1 μ 2 ,
f 40 = A 2 B 1 k 2 A 2 B 1 k 2 μ 1 A 2 B 1 k 2 μ 2 + A 2 B 1 k 2 μ 1 μ 2 .

3.1. Case 1 τ 1 = 0 , τ 2 > 0

For τ 1 = 0 , Equation (11) reduces to:
λ 4 + A 3 λ 3 + B 3 λ 2 + C 3 λ + D 3 + ( E 3 λ 3 + F 3 λ 2 + G 3 λ + H 3 ) e λ τ 2 = 0
with
A 3 = f 13 + f 23 ,   B 3 = f 12 + f 22 ,   C 3 = f 11 + f 21 ,   D 3 = f 10 + f 20 ,
E 3 = f 33 ,   F 3 = f 32 + f 42 ,   G 3 = f 31 + f 41 ,   H 3 = f 30 + f 40 .
Let λ = i ω 1 ,   ( ω 1 > 0 ) be the root of Equation (12). Then, we can get:
{ ( F 3 ω 1 2 H 3 ) sin ω 1 τ 2 + ( G 3 ω 1 E 3 ω 1 3 ) cos ω 1 τ 2 = A 3 ω 1 3 C 3 ω 1 ( G 3 ω 1 E 3 ω 1 3 ) sin ω 1 τ 2 ( F 3 ω 1 2 H 3 ) cos ω 1 τ 2 = ω 1 4 D 3 .
From Equation (13), we can obtain:
cos ω 1 τ 2 = ( F 3 A 3 E 3 ) ω 1 6 + ( A 3 G 3 H 3 + C 3 E 3 ) ω 1 4 + ( D 3 F 3 C 3 G 3 ) ω 1 2 D 3 H 3 ( G 3 ω 1 E 3 ω 1 3 ) 2 + ( F 3 ω 1 2 H 3 ) 2 .
Squaring both sides, adding both equations, and regrouping by powers of ω 1 , we can get:
ω 1 8 + ( A 3 2 E 3 2 ) ω 1 6 + ( 2 D 3 2 A 3 C 3 + 2 G 3 E 3 F 3 2 ) ω 1 4 + ( C 3 2 G 3 2 + 2 H 3 F 3 ) ω 1 2 + D 3 2 H 3 2 = 0 .
Let r 1 = ω 1 2 , then Equation (15) becomes:
r 1 4 + ( A 3 2 E 3 2 ) r 1 3 + ( 2 D 3 2 A 3 C 3 + 2 G 3 E 3 F 3 2 ) r 1 2 + ( C 3 2 G 3 2 + 2 H 3 F 3 ) r 1 + D 3 2 H 3 2 = 0 .
If the values of the coefficient are given, it is easy to get the root of Equation (16). In order to give the main results in this paper, we make the following assumption: (H1) Equation (16) has at least one positive root.
If condition (H1) holds, we know that Equation (15) has at least one positive root ω 10 , such that Equation (12) has a pair of purely imaginary roots ± i ω 10 . The corresponding critical value of the delay is
τ 20 = 1 ω 10 arccos ( F 3 A 3 E 3 ) ω 10 6 + ( A 3 G 3 H 3 + C 3 E 3 ) ω 10 4 + ( D 3 F 3 C 3 G 3 ) ω 10 2 D 3 H 3 ( G 3 ω 10 E 3 ω 10 3 ) 2 + ( F 3 ω 10 2 H 3 ) 2 .
Next, take the derivative with respect to τ 2 in Equation (12), we can obtain:
[ d λ d τ 2 ] 1 = ( 4 λ 3 + 3 A 3 λ 2 + 2 B 3 λ + C 3 ) e λ τ 2 + ( 3 E 3 λ 2 + 2 F 3 λ + G 3 ) λ ( E 3 λ 3 + F 3 λ 2 + G 3 λ + H 3 ) τ 2 λ .
Thus
Re [ d λ ( τ 20 ) d τ 2 ] λ = i ω 10 1 = Q 1 Q 3 + Q 2 Q 4 Q 1 2 + Q 2 1 ,
where
Q 1 = E 3 ω 10 4 G 3 ω 10 2 ,   Q 2 = H 3 ω 10 F 3 ω 10 3 ,
Q 3 = 4 ω 10 3 sin ω 10 τ 20 3 A 3 ω 10 2 cos ω 10 τ 20 2 B 3 ω 10 sin ω 10 τ 20 3 E 3 ω 10 2 + G 3 + C 3 cos ω 10 τ 20 ,
Q 4 = 2 B 3 ω 10 cos ω 10 τ 20 3 A 3 ω 10 2 sin ω 10 τ 20 4 ω 10 3 cos ω 10 τ 20 + 2 F 3 ω 10 + C 3 sin ω 10 τ 20 .
If condition (H2): Q 1 Q 3 + Q 2 Q 4 0 , then Re [ d λ ( τ 20 ) d τ 2 ] λ = i ω 10 1 0 . According to the Hopf bifurcation theorem in [19], we obtain the following results.
Theorem 1. 
If the conditions (H1)–(H2) hold, the equilibrium point E * of Equation (9) is asymptotically stable for τ 2 [ 0 , τ 20 ) and unstable for τ 2 > τ 20 ; Equation (9) undergoes a Hopf bifurcation when τ 2 = τ 20 .

3.2. Case 2 τ 1 > 0 , τ 2 > 0

We consider τ 1 as a parameter, and τ 2 in its stability region τ 2 [ 0 , τ 20 ) . At this point, the characteristic equation is Equation (11).
Let λ = i ω 2 ,   ( ω 2 > 0 ) is the root of Equation (11). Then, we can get:
{ A 4 cos ω 2 τ 1 + B 4 sin ω 2 τ 1 = C 4 A 4 sin ω 2 τ 1 B 4 cos ω 2 τ 1 = D 4 ,
where
A 4 = f 21 ω 2 f 23 ω 2 3 + f 42 ω 2 2 sin ω 2 τ 2 + f 41 ω 2 cos ω 2 τ 2 f 40 sin ω 2 τ 2
B 4 = f 22 ω 2 2 f 20 + f 42 ω 2 2 cos ω 2 τ 2 f 41 ω 2 sin ω 2 τ 2 f 40 cos ω 2 τ 2
C 4 = f 13 ω 2 3 f 11 ω 2 + f 33 ω 2 3 cos ω 2 τ 2 f 32 ω 2 2 sin ω 2 τ 2 f 31 ω 2 cos ω 2 τ 2 + f 30 sin ω 2 τ 2
D 4 = f 12 ω 2 2 ω 2 4 f 10 + f 33 ω 2 3 sin ω 2 τ 2 + f 32 ω 2 3 cos ω 2 τ 2 f 31 ω 2 sin ω 2 τ 2 f 30 cos ω 2 τ 2 .
From Equation (19), we can obtain:
cos ω 2 τ 1 = A 4 C 4 B 4 D 4 A 4 2 + B 4 2 .
Similar to Section 3.1, we can get:
ω 2 8 + A 5 ω 2 7 + B 5 ω 2 6 + C 5 ω 2 5 + D 5 ω 2 4 + E 5 ω 2 3 + F 5 ω 2 2 + G 5 ω 2 + H 5 = 0 ,
where
A 5 = 2 f 33 sin ω 2 τ 2 2 f 32 cos ω 2 τ 2 ,
B 5 = f 13 2 + 2 f 33 f 13 cos ω 2 τ 2 2 f 12 + f 33 2 f 23 2 ,
C 5 = 2 f 32 f 13 sin ω 2 τ 2 f 32 f 33 sin 2 ω 2 τ 2 + 2 f 33 f 12 sin ω 2 τ 2     + 2 f 32 f 12 cos ω 2 τ 2 + 2 f 31 sin ω 2 τ 2 + 2 f 42 f 23 sin ω 2 τ 2 ,
D 5 = 2 f 13 f 11 2 f 33 f 11 cos ω 2 τ 2 2 f 31 f 13 cos ω 2 τ 2 2 f 31 f 33 cos 2 ω 2 τ 2 + f 12 2       + 2 f 10 + f 32 2 2 f 31 sin 2 ω 2 τ 2 f 31 f 32 sin 2 ω 2 τ 2 + 2 f 30 cos ω 2 τ 2 + 2 f 21 f 23       + 2 f 41 f 23 cos ω 2 τ 2 f 22 2 f 42 2 2 f 42 f 22 cos ω 2 τ 2 ,
E 5 = f 31 f 32 sin 2 ω 2 τ 2 + 2 f 32 f 11 sin ω 2 τ 2 + 2 f 30 f 13 sin ω 2 τ 2 2 f 33 f 10 sin ω 2 τ 2       2 f 32 f 10 cos ω 2 τ 2 2 f 31 f 12 sin ω 2 τ 2 + 2 f 30 f 32 2 f 42 f 21 sin ω 2 τ 2 2 f 40 f 23 sin ω 2 τ 2 ,
F 5 = f 11 2 + 2 f 31 f 11 cos ω 2 τ 2 2 f 10 f 12 + f 31 2 2 f 30 f 12 cos ω 2 τ 2 f 21 2 + 2 f 20 f 22      + 2 f 42 f 20 cos ω 2 τ 2 f 41 2 + 2 f 41 f 22 sin ω 2 τ 2 + 2 f 40 f 22 cos ω 2 τ 2 + 2 f 42 f 40 ,
G 5 = f 30 f 31 sin 2 ω 2 τ 2 2 f 30 f 11 sin ω 2 τ 2 + 2 f 31 f 10 sin ω 2 τ 2 + f 30 f 31 sin 2 ω 2 τ 2       2 f 41 f 21 cos ω 2 τ 2 + 2 f 40 f 21 sin ω 2 τ 2 2 f 41 f 20 sin ω 2 τ 2 ,
H 5 = f 10 2 + f 30 2 + 2 f 30 f 10 cos ω 2 τ 2 f 20 2 f 40 2 2 f 40 f 20 cos ω 2 τ 2 .
We assume that (H3): Equation (21) has k positive roots, they are ω 2 i , i = 1 , 2 , ... , k . The corresponding delay parameter values are:
τ 1 i = 1 ω 2 i arccos A 4 C 4 B 4 D 4 A 4 2 + B 4 2 + 2 j π ω 2 i   ( i = 1 , 2 , ... , k ; j = 0 , 1 , 2 , ... , ) .
We let
τ 10 = min { τ 1 i ( j ) | i = 1 , 2 , ... , k ; j = 0 , 1 , ... , } .
When τ 1 = τ 10 , Equation (11) has a pair of purely imaginary roots ± i ω 20 .
Next, taking the derivatives for τ 1 in Equation (11), we can get:
[ d λ d τ 1 ] 1 = A 6 e λ ( τ 1 + τ 2 ) + B 6 e λ τ 2 + C 6 e λ τ 1 + D 6 E 6 e λ τ 2 + F 6 τ 1 λ ,
where
A 6 = 4 λ 3 + 3 f 13 λ 2 + 2 f 12 λ + f 11 ,   B 6 = 3 f 23 λ 2 + 2 f 22 λ + f 21 ,
C 6 = τ 2 f 33 λ 3 + ( 3 f 33 τ 2 f 32 ) λ 2 + ( 2 f 32 τ 2 f 31 ) λ + f 31 τ 2 f 30 ,
D 6 = τ 2 f 42 λ 2 τ 2 f 41 λ τ 2 f 40 + 2 f 42 + f 41 ,
E 6 = f 23 λ 4 + f 22 λ 3 + f 21 λ 2 + f 20 λ ,   F 6 = f 42 λ 3 + f 41 λ 2 + f 40 λ .
Thus
Re [ d λ ( τ 10 ) d τ 1 ] λ = i ω 20 1 = Q 5 Q 7 + Q 6 Q 8 Q 5 2 + Q 6 2 ,
Q 5 = f 23 ω 20 4 cos ω 20 τ 2 + f 22 ω 20 3 sin ω 20 τ 2 f 21 ω 20 2 cos ω 20 τ 2 f 20 ω 20 sin ω 20 τ 2 f 41 ω 20 2 ,
Q 6 = f 23 ω 20 4 sin ω 20 τ 2 f 22 ω 20 3 cos ω 20 τ 2 f 21 ω 20 2 sin ω 20 τ 2 + f 20 ω 20 cos ω 20 τ 2 f 42 ω 20 3 + f 40 ω 20 ,
Q 7 = 4 ω 20 3 sin ω 20 ( τ 10 + τ 2 ) 3 f 13 ω 20 2 cos ω 20 ( τ 10 + τ 2 ) 2 f 12 ω 20 sin ω 20 ( τ 10 + τ 2 )      + f 11 cos ω 20 ( τ 10 + τ 2 ) 3 f 23 ω 20 2 cos ω 20 τ 2 2 f 22 ω 20 sin ω 20 τ 2 + f 21 cos ω 20 τ 2      τ 2 f 33 ω 20 3 sin ω 20 τ 10 ( 3 f 33 τ 2 f 32 ) ω 20 2 cos ω 20 τ 10 ( 2 f 32 τ 2 f 31 ) ω 20 sin ω 20 τ 10      + ( f 31 τ 2 f 30 ) cos ω 20 τ 10 + τ 2 f 42 ω 20 2 τ 2 f 40 + 2 f 42 + f 41 ,
Q 8 = 4 ω 20 3 cos ω 20 ( τ 10 + τ 2 ) 3 f 13 ω 20 2 sin ω 20 ( τ 10 + τ 2 ) + 2 f 12 ω 20 cos ω 20 ( τ 10 + τ 2 )       + f 11 sin ω 20 ( τ 10 + τ 2 ) 3 f 23 ω 20 2 sin ω 20 τ 2 + 2 f 22 ω 20 cos ω 20 τ 2 + f 21 sin ω 20 τ 2       + τ 2 f 33 ω 20 3 cos ω 20 τ 10 ( 3 f 33 τ 2 f 32 ) ω 20 2 sin ω 20 τ 10 + ( 2 f 32 τ 2 f 31 ) ω 20 cos ω 20 τ 10       + ( f 31 τ 2 f 30 ) sin ω 20 τ 10 τ 2 f 41 ω 20 .
If condition (H4): Q 5 Q 7 + Q 6 Q 8 0 , then Re [ d λ ( τ 10 ) d τ 1 ] λ = i ω 20 1 0 . According to the Hopf bifurcation theorem in [19], we have the following results.
Theorem 2. 
For τ 2 [ 0 , τ 20 ) , if the conditions (H3)–(H4) hold, then the equilibrium E * of Equation (9) is asymptotically stable for τ 1 [ 0 , τ 10 ) and unstable when τ 1 > τ 10 . Equation (9) has a Hopf bifurcation at τ 1 = τ 10 .

4. Numerical Simulation and Analysis

In order to examine the correctness of the theoretical analysis, the numerical simulation and analysis are carried out in this part. We set a = 2 , θ = 0.036 , c = 0.1 , ν 1 = 25 , ν 2 = 1.3 , μ 1 = 0.4 , μ 2 = 0.6 , k = 0.05 .The initial values of x 1 , x 2 , q 1 , q 2 are 0.1, 0.2, 0.3, and 0.4, respectively. Then Equation (9) becomes:
{ x ˙ 1 = 0.05 x 1 q 1 , x ˙ 2 = 0.05 x 2 q 2 , q ˙ 1 = ( 0.28 q 1 0.42 q 1 ( t τ 1 ) 0.32 x 2 + 0.21 x 1 , q ˙ 2 = ( 2.6 0.52 q 1 0.78 q 1 ( t τ 1 ) 0.78 q 2 0.52 q 2 ( t τ 2 ) )       ( 0.0084 q 2 0.0056 q 2 ( t τ 2 ) 0.0064 x 1 + 0.00427 x 2 ) ; .

4.1. The Influence of τ 2 on the Stability of Equation (25)

For τ 1 = 0 , from Equations (15)–(18) we can get ω 10 = 21.462 , τ 20 = 2.829 , Q 1 Q 3 + Q 2 Q 4 = 132.746 0 , so (H1)–(H2) hold. By Theorem 1 we can know that when τ 2 [ 0 , τ 20 ) , Equation (25) is asymptotically local stable; when τ 2 > τ 20 , Equation (25) is unstable. Figure 1, Figure 2 and Figure 3 show the above properties.
Figure 1a is a bifurcation diagram about τ 2 of Equation (25). It can be seen that as τ 2 increases, the system undergoes bifurcation at τ 2 = 2.829 . That is to say, the system will lose its steady state, even descending into chaos. Figure 1b represents the largest Lyapunov exponent with respect to τ 2 . A Lyapunov index value less than 0 indicates that the system is stable; a value greater than 0 indicates that the system is instable, even chaotic; when the value equals 0, the system undergoes bifurcation. From Figure 1b we can see that the system has bifurcation at τ 2 = 2.829 . This is consistent with the conclusions of the theoretical analysis.
Figure 2 and Figure 3 further show the influence of τ 2 on the stability of the system. When τ 2 is on both sides of the value of τ 20 , the stability of the system changes significantly. Figure 2 is stable for τ 2 = 2 < τ 20 = 2.829 ; Figure 3 is unstable as τ 2 = 3.5 > τ 20 = 2.829 .
Entropy is a tool to measures the degree of chaos in a system. Greater entropy means that the system is more chaotic, while smaller entropy means the system is more stable. It can be seen from Figure 4 that entropy is stable when τ 2 < 2.829 . The starting point of the rapid increase of entropy is at τ 2 = 2.829 . As τ 2 increases, the entropy increases gradually. This shows that with a greater τ 2 , the system is more chaotic and it will take longer for it to return to stability.

4.2. The Influence of τ 1 on the Stability of Equation (25)

When we let τ 2 = 1.5 [ 0 , τ 20 ) , from Equations (21)–(24) we can obtain ω 20 = 18.492 , τ 10 = 2.648 , Q 1 Q 3 + Q 2 Q 4 = 81.735 0 , so that (H3)–(H4) hold. Using Theorem 2 we find that when τ 1 [ 0 , τ 10 ) , Equation (25) is asymptotically local and stable; as τ 1 > τ 10 , Equation (25) is unstable. Figure 5, Figure 6 and Figure 7 show the above properties.
The analysis process is similar to case 1. Figure 5 shows the evolution of the system in detail with the increase of τ 1 . When τ 1 = 2.648 , the system undergoes bifurcation. The bigger τ 1 is, the more chaotic the system becomes. Therefore, the enterprise must take reasonable time delay parameters to ensure the system is in a stable state.
Figure 6 and Figure 7 show that the system has a bifurcation when τ 10 = 2.648 . It is also the demarcation point between the stable and unstable conditions of the system.
Basin of attraction is an effective tool to measure the stability of the system. When the basin area is larger, the system is more stable. From Figure 8, it is known that with an increase of τ 1 , the area of the basin decreases gradually, which means the stable elements in the system are reduced and the unstable factors are increased.

4.3. The Influence of τ 1 and τ 2 on the Stability of Equation (25)

In this part, we study the effects of τ 1 and τ 2 on the stability of the system. Let τ 1 [ 1 , 3 ] and τ 2 [ 1 , 3 ] . This ensures that a combination of τ 1 and τ 2 appears bifurcation in the region of [ 1 , 3 ] × [ 1 , 3 ] plane. From Figure 9 we see that with an increase of τ 1 and τ 2 , the system gradually loses its stability and enters a chaotic state. When the system is in a state of chaos (the blue region in Figure 10), the maximum value of q 1 is 82.15 for τ 1 = 3 , τ 2 = 3 ; the minimum value of q 1 is −34.93 for τ 1 = 2.75 , τ 2 = 3 . When the system is in a stable state (the green region in Figure 10), the value of q 1 is basically stable at 1.284. Therefore, in order to maintain the stability of the system, it is necessary to keep τ 1 and τ 2 in the green region of Figure 10.
From Section 4.2, we know that when τ 1 = 2.648 and τ 2 = 1.5 , the system undergoes bifurcation. It can be seen from Figure 10 that the point (2.648, 1.5) is on the boundary between stability and instability. This is consistent with the theoretical analysis in Section 4.2.
The analysis is similar to Figure 4. When τ 1 and τ 2 are in the green region of Figure 10, the entropy is essentially the same, and the system is in a stable state. As τ 1 and τ 2 are in the blue region of Figure 10, the entropy increases sharply, which means that the system is instable and the degree of chaos is increasing. The above properties are shown in Figure 11.

4.4. The Influence of τ 1 , τ 2 and μ 1 on the Stability of Equation (25)

This part focuses on the influence of τ 1 , τ 2 and μ 1 on the stability of q 1 . It can be seen from Figure 12 that the changes in τ 1 and μ 1 have a greater impact on the stability of q 1 , while τ 2 has little effect on q 1 . With the increase of τ 1 , the system becomes unstable, and the system will be changed from an instable state to a stable state with the increase of μ 1 . That is, when ( τ 1 , τ 2 , μ 1 ) is in the green region of Figure 13, the system is stable; when ( τ 1 , τ 2 , μ 1 ) is in the blue region of Figure 13, the system is not stable, and even in a state of chaos. Therefore, when the hydropower enterprises are making power output decisions, they must ensure the combination of parameter values is in the green area in Figure 13.

5. Chaos Control

Through the previous analysis, we know that if the parameter values are not reasonable, the system will be instable, even chaotic. A system in a state of chaos will lead to market fluctuations. Therefore, some measures must be taken to maintain the stability of the system. The instable system can be controlled to restore it to a stable state.
From the analysis in Section 4.2, we know that when τ 1 = 3 , τ 2 = 1.5 , the system is chaotic; the effect is shown in Figure 7. From Figure 10 of Section 4.3, τ 1 = 3 , τ 2 = 1.5 is in the blue region; it can be clearly seen that the system is chaotic. We consider the following control system using the delayed feedback control method:
{ x ˙ 1 = 0.05 x 1 q 1 , x ˙ 2 = 0.05 x 2 q 2 , q ˙ 1 = ( 0.28 q 1 0.42 q 1 ( t τ 1 ) 0.32 x 2 + 0.21 x 1 , q ˙ 2 = ( 2.6 0.52 q 1 0.78 q 1 ( t τ 1 ) 0.78 q 2 0.52 q 2 ( t τ 2 ) )       ( 0.0084 q 2 0.0056 q 2 ( t τ 2 ) 0.0064 x 1 + 0.00427 x 2 ) μ ( q 2 q 2 ( t τ 2 ) ) . .
Only adjusting the control parameter μ while the other parameters are kept constant to realize the control of chaos. Similar to the derivation process in Section 3, the system can obtain bifurcation point μ = 0.07909 for τ 1 = 3 , τ 2 = 1.5 . From Figure 14, we can see that Equation (26) changes from a chaotic state to a steady state as μ increases. When μ = 0.07909 , the system undergoes bifurcation. That is to say, when μ ( 0 , 0.07909 ) , the system is chaotic; when μ ( 0.07909 , 0.2 ] , the system is stable.
In order to verify the correctness of the above analysis, set μ = 0.03 < 0.07909 ; according to the analysis, the system should be chaotic at this time. Figure 15 proves the correctness of the conclusion. When we let μ = 0.15 > 0.07909 , it can be seen that the system is stable at this time in Figure 16; that is, the chaotic system has been effectively controlled.

6. Conclusions

A game model of the duopoly with two delays in the hydropower market is established. The influence of time delay parameters on the stability of the game model is studied. Firstly, the existence and stability of Hopf bifurcation are studied, and the conditions of bifurcation are given. The analysis is carried out in two aspects: τ 1 = 0 , τ 2 > 0 and τ 1 > 0 , τ 2 > 0 . Secondly, the numerical simulation and analysis are carried out on the theoretical derivation. The specific analysis is developed from four aspects: the influence of τ 2 on the stability of the system for τ 1 = 0 simulated by using time series, attractor, bifurcation diagram, Lyapunov exponent, and entropic; τ 2 fixed, the influence of τ 1 on the stability of the system discussed by using time series and basin of attraction; the influence of τ 1 and τ 2 on the stability of the system displayed by 3D surface chart, 2D parameter bifurcation diagram, and 3D entropic diagram; the influence of τ 1 , τ 2 , and μ 1 on the stability of the system simulated through a 4D cubic chart and 3D parameter bifurcation diagram. The stability regions of the system are given for each aspect of the specific analysis. The conclusions of the above analysis can provide decision-making guidelines for enterprises to maintain market stability. Finally, the chaotic system is controlled effectively by the method of delayed feedback control, which can efficiently return a chaotic system to a stable state. For future study, we can consider the influence of more factors on the stability of the system, so that it is closer to reality.

Acknowledgments

We thank the reviewers for their careful reading and helpful comments on the revision of the paper. The research was supported by the National Natural Science Foundation of China (Nos. 71571131 and 61273231) and the Doctoral Fund of the Ministry of Education of China (Grant No. 20130032110073).

Author Contributions

Zhihui Han built the economic model; Junhai Ma provided economic interpretation of the conclusions; Fengshan Si performed mathematical derivation; Wenbo Ren carried out numerical simulation; they wrote this research articles together. All authors have read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Weidou, N.; Johansson, T.B. Energy for Sustainable Development in China. Energy Policy 2004, 32, 1225–1229. [Google Scholar] [CrossRef]
  2. Nwaobi, G.C. Emission Policies and the Nigerian Economy: Simulations from a Dynamic Applied General Equilibrium Model. Energy Econ. 2004, 26, 921–936. [Google Scholar] [CrossRef] [Green Version]
  3. Vera, I.; Langlois, L. Energy Indicators for Sustainable Development. Energy 2007, 32, 875–882. [Google Scholar] [CrossRef]
  4. Omri, A. An international literature survey on energy-economic growth nexus: Evidence from country-specific studies. Renew. Sustain. Energy Rev. 2014, 38, 951–959. [Google Scholar] [CrossRef]
  5. Dajka, J.; Łobejko, M.; Sładkowski, J. Payoffs and Coherence of a Quantum Two-Player Game in a Thermal Environment. Entropy 2015, 17, 7736–7751. [Google Scholar] [CrossRef]
  6. Harré, M.S.; Bossomaier, T. Strategic islands in economic games: Isolating economies from better outcomes. Entropy 2014, 16, 5102–5121. [Google Scholar] [CrossRef]
  7. Ma, J.H.; Ji, W.Z. Complexity of repeated game model in electric power triopoly. Chaos Solitons Fractals 2009, 40, 1735–1740. [Google Scholar] [CrossRef]
  8. Batabyal, A.A.; Beladi, H. A Stackelberg Game Model of Trade in Renewable Resources with Competitive Sellers. Rev. Int. Econ. 2006, 14, 136–147. [Google Scholar] [CrossRef]
  9. Sun, Z.H.; Ma, J.H. Complexity of triopoly price game in Chinese cold rolled steel market. Nonlinear Dyn. 2012, 67, 2001–2008. [Google Scholar] [CrossRef]
  10. Halkos, G.; Papageorgiou, G. Extraction of non-renewable resources: A differential game approach. Arch. Econ. Hist. 2008, 1, 5–22. [Google Scholar]
  11. Liao, G.C. A novel evolutionary algorithm for dynamic economic dispatch with energy saving and emission reduction in power system integrated wind power. Energy 2011, 36, 1018–1029. [Google Scholar] [CrossRef]
  12. Yoon, K.H.; Ratti, R.A. Energy price uncertainty, energy intensity and firm investment. Energy Econ. 2011, 33, 67–78. [Google Scholar] [CrossRef]
  13. Chang-Jian, C.W.; Chang, S.M. Bifurcation and chaos analysis of spur gear pair with and without nonlinear suspension. Nonlinear Anal. Real World Appl. 2011, 12, 979–989. [Google Scholar] [CrossRef]
  14. Lü, J.H.; Chen, G.; Celikovský, S. Bridge the gap between the Lorenz system and the Chen system. Int. J. Bifurc. Chaos 2002, 12, 2917–2926. [Google Scholar] [CrossRef]
  15. Celikovský, S.; Chen, G. On a generalized Lorenz canonical form of chaotic systems. Int. J. Bifurc. Chaos 2002, 12, 1789–1812. [Google Scholar] [CrossRef]
  16. Sun, M.; Tian, L.X.; Fu, Y.; Qian, W. Dynamics and adaptive synchronization of the energy resource system. Chaos Solitons Fractals 2007, 31, 879–888. [Google Scholar] [CrossRef]
  17. Sun, M.; Tian, L.X.; Zeng, C.Y. The energy resources system with parametric perturbations and its hyperchaos control. Nonlinear Anal. Real World Appl. 2009, 10, 2620–2626. [Google Scholar] [CrossRef]
  18. Sun, M.; Wang, X.F.; Chen, Y.; Tian, L.X. Energy resources demand-supply system analysis and empirical research based on non-linear approach. Energy 2011, 36, 5460–5465. [Google Scholar] [CrossRef]
  19. Hassard, B.D.; Kazarinoff, N.D.; Wan, Y.H. Theory and Applications of Hopf Bifurcation; Cambridge University Press: Cambridge, UK, 1981. [Google Scholar]
Figure 1. The influence of τ 2 on the stability of Equation (25) when τ 1 = 0 . (a) Bifurcation diagram; (b) the Largest Lyapunov exponent.
Figure 1. The influence of τ 2 on the stability of Equation (25) when τ 1 = 0 . (a) Bifurcation diagram; (b) the Largest Lyapunov exponent.
Entropy 18 00317 g001
Figure 2. Equation (25) is stable when τ 2 = 2 < τ 20 = 2.829 for τ 1 = 0 . (a) Time series plot; (b) attractor.
Figure 2. Equation (25) is stable when τ 2 = 2 < τ 20 = 2.829 for τ 1 = 0 . (a) Time series plot; (b) attractor.
Entropy 18 00317 g002
Figure 3. Equation (25) is instable when τ 2 = 3.5 > τ 20 = 2.829 for τ 1 = 0 . (a) Time series plot; (b) attractor.
Figure 3. Equation (25) is instable when τ 2 = 3.5 > τ 20 = 2.829 for τ 1 = 0 . (a) Time series plot; (b) attractor.
Entropy 18 00317 g003
Figure 4. Entropy diagram with an increase of τ 2 for τ 1 = 0 .
Figure 4. Entropy diagram with an increase of τ 2 for τ 1 = 0 .
Entropy 18 00317 g004
Figure 5. The influence of τ 1 on the stability of Equation (25) when τ 2 = 1.5 . (a) Bifurcation diagram; (b) the Largest Lyapunov exponent.
Figure 5. The influence of τ 1 on the stability of Equation (25) when τ 2 = 1.5 . (a) Bifurcation diagram; (b) the Largest Lyapunov exponent.
Entropy 18 00317 g005
Figure 6. Equation (25) is stable when τ 1 = 2 < τ 10 = 2.648 for τ 2 = 1.5 . (a) Time series plot; (b) attractor.
Figure 6. Equation (25) is stable when τ 1 = 2 < τ 10 = 2.648 for τ 2 = 1.5 . (a) Time series plot; (b) attractor.
Entropy 18 00317 g006
Figure 7. Equation (25) is instable when τ 1 = 3 > τ 10 = 2.648 for τ 2 = 1.5 . (a) Time series plot; (b) attractor.
Figure 7. Equation (25) is instable when τ 1 = 3 > τ 10 = 2.648 for τ 2 = 1.5 . (a) Time series plot; (b) attractor.
Entropy 18 00317 g007
Figure 8. The change of basin of attraction with an increase of τ 1 for τ 2 = 1.5 . (a) τ 1 = 2.3 ; (b) τ 1 = 2.5 ; (c) τ 1 = 2.7 .
Figure 8. The change of basin of attraction with an increase of τ 1 for τ 2 = 1.5 . (a) τ 1 = 2.3 ; (b) τ 1 = 2.5 ; (c) τ 1 = 2.7 .
Entropy 18 00317 g008aEntropy 18 00317 g008b
Figure 9. The influences of τ 1 and τ 2 on q 1 .
Figure 9. The influences of τ 1 and τ 2 on q 1 .
Entropy 18 00317 g009
Figure 10. 2D parameter bifurcation in the ( τ 1 , τ 2 ) plane, where different colors represent different price regions: stability region (green) and chaotic region (blue). For interpretation of the references to color in this figure, the reader is referred to the web version of this article.
Figure 10. 2D parameter bifurcation in the ( τ 1 , τ 2 ) plane, where different colors represent different price regions: stability region (green) and chaotic region (blue). For interpretation of the references to color in this figure, the reader is referred to the web version of this article.
Entropy 18 00317 g010
Figure 11. Entropy diagram with increase of τ 1 and τ 2 .
Figure 11. Entropy diagram with increase of τ 1 and τ 2 .
Entropy 18 00317 g011
Figure 12. The influence of τ 1 , τ 2 and μ 1 on q 1 .
Figure 12. The influence of τ 1 , τ 2 and μ 1 on q 1 .
Entropy 18 00317 g012
Figure 13. 3D parameter bifurcation in the ( τ 1 , τ 2 , μ 1 ) plane, where different colors represent different regions of q 1 : stability region (green) and chaotic region (blue). For interpretation of the references to color in this figure, the reader is referred to the web version of this article.
Figure 13. 3D parameter bifurcation in the ( τ 1 , τ 2 , μ 1 ) plane, where different colors represent different regions of q 1 : stability region (green) and chaotic region (blue). For interpretation of the references to color in this figure, the reader is referred to the web version of this article.
Entropy 18 00317 g013
Figure 14. The influence of μ on the stability of Equation (26). (a) Bifurcation diagram; (b) the Largest Lyapunov exponent.
Figure 14. The influence of μ on the stability of Equation (26). (a) Bifurcation diagram; (b) the Largest Lyapunov exponent.
Entropy 18 00317 g014
Figure 15. Equation (26) is instable when μ = 0.03 < 0.07909 for τ 1 = 3 , τ 2 = 1.5 . (a) Time series plot; (b) attractor.
Figure 15. Equation (26) is instable when μ = 0.03 < 0.07909 for τ 1 = 3 , τ 2 = 1.5 . (a) Time series plot; (b) attractor.
Entropy 18 00317 g015
Figure 16. Equation (26) is stable when μ = 0.15 > 0.07909 for τ 1 = 3 , τ 2 = 1.5 . (a) Time series plot; (b) attractor.
Figure 16. Equation (26) is stable when μ = 0.15 > 0.07909 for τ 1 = 3 , τ 2 = 1.5 . (a) Time series plot; (b) attractor.
Entropy 18 00317 g016

Share and Cite

MDPI and ACS Style

Han, Z.; Ma, J.; Si, F.; Ren, W. Entropy Complexity and Stability of a Nonlinear Dynamic Game Model with Two Delays. Entropy 2016, 18, 317. https://0-doi-org.brum.beds.ac.uk/10.3390/e18090317

AMA Style

Han Z, Ma J, Si F, Ren W. Entropy Complexity and Stability of a Nonlinear Dynamic Game Model with Two Delays. Entropy. 2016; 18(9):317. https://0-doi-org.brum.beds.ac.uk/10.3390/e18090317

Chicago/Turabian Style

Han, Zhihui, Junhai Ma, Fengshan Si, and Wenbo Ren. 2016. "Entropy Complexity and Stability of a Nonlinear Dynamic Game Model with Two Delays" Entropy 18, no. 9: 317. https://0-doi-org.brum.beds.ac.uk/10.3390/e18090317

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop