# Nonlinear Complex Dynamics of Carbon Emission Reduction Cournot Game with Bounded Rationality

Authors: LiuWei Zhao
Journal: Complexity
 (2017)

---

## Abstract

Based on the hypothesis of participant’s bounded rationality, our study formulated a novel Cournot duopoly game model of carbon emission reduction and, subsequently, analyzed the dynamic adjustment mechanism of emission reduction for enterprises. The existence and stability of the equilibrium solution of game are further discussed by the nonlinear dynamics theory. Our findings revealed that the parameters have key significance on the dynamic properties of the system. However, when the adjustment speed gets too large, the system loses the original stability and vividly demonstrates complex chaos phenomenon. Higher market prices in carbon trading have an outstanding impact on the stability of the system, which easily leads to system instability. Our study further controlled the chaos behavior of the power system by the delay feedback control. The results of the numerical analysis depict that the unstable behavior of the dynamic system can be controlled efficiently and quickly, in the quest to restore back a stable and orderly market. Our novel method is proved to have provided decision makers with effective solution to market instability.


---

## Body

## 1. Introduction

Massive economic growth has resulted in severe environmental issues such as environmental pollution and global climate warming. There is, therefore, huge concerns arising from the United Nations Framework Convention on Climate Change (UNFCCC) and various countries, to develop a series of related policies and regulations to curb these environmental menaces. Carbon emission trading (carbon trading for short) is one of the key market mechanism approaches used to reduce global greenhouse, thus carbon dioxide emission reduction. On 9 May 1992, an environmental treaty on the issue of climate change was negotiated and passed by international governments and the United Nations Framework Convention. In December 1997 in Tokyo, Japan, the “Kyoto Protocol” was issued. The Kyoto Protocol launched a market mechanism as a new solution to the problems of greenhouse gas emission reduction with a focus on the reduction of carbon dioxide since carbon emission has been the principal greenhouse gas. This entails trading carbon as a commodity: emitting carbon dioxide beyond the limits generates a potential cost; and conversely, those who are within a defined scope hold a potential value.Despite the fact that China is a developing country and not subject to stricter obligation of reducing emissions under the Kyoto Protocol, this global climate changes and the carbon trading in developed countries around the world has brought tremendous changes in China. Consequently, China is playing an active role in the establishment of carbon trading pilots in larger cities, such as Beijing, Shanghai, and other provinces and cities. China has already initiated such pilot projects and far advanced in its implementation. The main participants in the carbon trading market are the industrial enterprises. This carbon trading impacts on enterprises’ production processes, technology, profits, and other critical aspects. Therefore, it is of great significance to study and analyze the influence of carbon trading on the production, decision-making, and profits of the enterprise, which would provide sufficient theoretical basis for the perfection and implementation of the right carbon trading policy.The main concerns about carbon emissions trading can be categorized into three aspects: firstly, the initial allocation of carbon emissions rights [1, 2]; secondly, the price of carbon emissions trading [3]; and thirdly the economic benefits arising from carbon emissions trading [4]. With regard to the development of carbon emissions trading, the industry structure of the enterprises has a great impact on the carbon trading. For instance, in the oligopolistic industry market, the carbon trading price is more susceptible to the impact of enterprises’ transactions, in consequence, making the carbon emissions trading prices deviate from the market’s equilibrium price. Therefore, the research on the game behavior of carbon emissions trading in oligopoly market has become a vital problem that needs to be solved urgently.The two classical models in the theory of oligopoly are those of Cournot [5, 6] and Bertrand [7, 8]. In Cournot’s model, the firms choose to compete on output quantity, and in Bertrand’s model, they choose to compete on price. Both models can be interpreted as static games where decisions are made simultaneously and where each firm maximizes its own profit, in a context of perfect and complete information [9, 10]. This paper introduces carbon emission trading in the classic Cournot output game framework and formulates a novel Cournot duopoly game model of carbon emission reduction based on the hypothesis of participant's bounded rationality. The dynamic adjustment mechanism of emission reduction is also analyzed in this paper. According to the hypothesis that the market inverse demand function and the cost function are linear functions and the participants have bounded rationality, our study constructs the nonlinear dynamic Cournot game model to analyze the local stability of the Nash equilibrium. The stability of the boundary equilibrium for the four-dimensional dynamical system is further analyzed based on the stability of the discrete system. Consequently, the parameters conditions of interior equilibrium and stability are obtained according to Schur-Cohn criterion. Based on the analysis of all equilibrium’s stability of the nonlinear investment game system, our study adopts the numerical simulation results to analyze the influence of the model parameters on the stability and complexity of the dynamic system. The chaos behavior of the system is effectively controlled by the delayed feedback control method.This article is organized as follows: in Section2, the model is established based on the hypothesis of bounded rationality. In Section 3, we discuss the existence and local stability of the equilibrium points in the system. In Section 4, we demonstrate the dynamic features of this system with numerical simulations, including bifurcation diagram, phase portrait, and sensitive dependence on initial conditions. In Section 5, delayed feedback control is utilized to stabilize the chaotic behaviors of the system.


## 2. Model Hypotheses and Establishment

Our study assumes that the two production enterprises develop and produce homogeneous products. Each enterprise's action strategy is to choose the amount of investment in emission reduction during various periods and meet the requirements: two participating companies develop the corresponding emission reduction investment strategy in the discrete timeline. TheKi(t-1) reflects the company’s capital stock in the period of the t-1, xi(t) for its single-stage carbon emission reduction in the period of t, and emission reduction cost coefficient is βi(β>0), so the enterprise’s single-stage emission reduction investment is βixi(t). Due to the existence of a certain amount of capital depreciation, investment capital stock Ki(t-1) flows into the next economic period, and it remains θKi(t-1), among which θ (0<θ<1) stands for the residual ratio. Thus, we can reach that the relational expression in the adjacent two periods of investment capital stock of enterprise i is(1)Kit=θKit-1+βixit,i=1,2.For each particular enterprise’s products amount in a certain period of time, production potential is determined by the accumulation amount of corporate investment capital in the periodKi(t). Our study utilizes a linear form qi(t)=BiKi(t) to represent [11], among which positive constant Bi denotes the technological level of enterprises for innovative production on emerging products. Thus, we can reach that qi(t)=BiθKi(t-1)+βixi(t),i=1,2.Our study assumes that the sales price of this kind of products on the market is a linear inverse demand functionp(t)=a-bQ(t), a,b>0, and Qt=q1(t)+q2(t) represents the sum of the two firms' supply in the market. At the same time, the firm's linear production cost function Ci(qi(t))=ciqi(t)(i=1,2), marginal costs c1 and c2 are positive numbers.The relationship between the amount of carbon dioxide produced by the enterprise during the production process and its productionqi(t) is expressed as a linear relationship in the form of kiqi(t). ki(ki>0) is the coefficient of carbon dioxide produced by the enterprise i. Therefore, the enterprise’s total carbon emission is kiqi(t)-xi(t) in the single period of t. Each enterprise in the t period has a certain quota of carbon emissions rights-yi(t), and the enterprise’s carbon emissions trading volume is yi(t)-(kiqi(t)-xi(t)), and the market price of carbon trading is P. If yi(t)-(kiqi(t)-xi(t))>0, it means that the company sells the remaining carbon emissions; on the other hand, if the company needs to buy the lacking carbon emissions, yi(t)-(kiqi(t)-xi(t))=-(yj(t)-(kjqj(t)-xj(t))).Based on the above assumptions about the function relationship among the variables, the profit of the enterprisei in the period t is calculated as follows:(2)πix1t,x2t=qitpt-ci-βixit+Pyit-kiqit-xit,i=1,2.The specific expression of each variable is substituted into (2) and reaches first partial derivative about xi(t) on πi(x1(t),x2(t)), and the marginal profits of the two enterprises are(3a)φ1t=∂π1x1t,x2t∂x1t=P+β1-bB1B2θK2t-1+β2x2t+aB1-c1B1-Pk1B1-1-2bB12θK1t-1+β1x1t,(3b)φ2t=∂π2x1t,x2t∂x2t=P+β2-bB1B2θK1t-1+β1x1t+aB2-c2B2-Pk2B2-1-2bB22θK2t-1+β2x2t.In reality, due to the asymmetry of market information and the constraints of their own conditions, enterprises can generally reach only limited degree of rationality in the decision-making and cannot make a complete control on the future market demand; therefore, the two production enterprises will rely on their own local knowledge about the marginal profits to adjust the emission reduction strategies during the next period. In other words, if the business observes the positive marginal profit in the current period, then it will increase emission reduction in the period oft+1 and vice versa. As a result, the enterprise's dynamic adjustment mechanism can be expressed as(4)xit+1=xit+αixitφit,i=1,2,where αi(xi(t)) is a positive adjustment function and reflects the company’s adjustment range of emission reduction based on the marginal profits. Our study considers the form of the adjustment function in linear form as αi(xi(t))=αixi(t), the coefficient αi>0 when the enterprise adjusts the rate of emission reduction strategy according to the marginal profit signal. Thus, (4) is(5)xit+1=xit+αixitφit,i=1,2.By connecting (1), (3a), (3b), and (5), in order to obtain the same expression by replacing Ki(t-1) with Ii(t), as shown from (6), we obtained a four-dimensional discrete dynamical model by replacing Kit-1 with Ii(t) as follows:(6)x1t+1=x1t+α1x1tP+β1-bB1B2θI2t+β2x2t+aB1-c1B1-Pk1B1-1-2bB12θI1t+β1x1tx2t+1=x2t+α2x2tP+β2-bB1B2θI1t+β1x1t+aB2-c2B2-Pk2B2-1-2bB22θI2t+β2x2tI1t+1=θI1t+β1x1tI2t+1=θI2t+β2x2t.The discrete dynamical system (6) is based on the assumption that the market inverse demand function and the production cost function are all linear forms, describing the competition model of the emission reduction strategy between the two production enterprises. The game participants with bounded rational expectation will adjust the emission reduction strategy in different periods according to the marginal profits.


## 3. Analysis of Equilibrium Point’s Stability

In the dynamical system (6), let xi(t+1)=xi(t) and Ii(t+1)=Ii(t)(i=1,2); we then get(7)x1tP+β1-bB1B2θI2t+β2x2t+aB1-c1B1-Pk1B1-1-2bB12θI1t+β1x1t=0x2tP+β2-bB1B2θI1t+β1x1t+aB2-c2B2-Pk2B2-1-2bB22θI2t+β2x2t=01-θβ1I1t=x1t1-θβ2I2t=x2t.In solving (7), we obtained the four equilibrium points of the dynamical system (6) as follows:(8)E0=0,0,0,0,E1=1-θP-β1+B1β1a-c1-Pk12bB12β12,0,P-β1+B1β1a-c1-Pk12bB12β1,0,E2=0,1-θP-β2+B2β2a-c2-Pk22bB22β22,0,P-β2+B2β2a-c2-Pk22bB22β1E∗=1-θW13bB12β12B2β2,1-θW23bB1B22β1β22,W13bB12B2β1β2,W23bB1B22β1β2,where W1=2B2P-β1β2+B1β1β21+B2a-2c1+c2-2Pk1+Pk2-P, W2=B1β12P+B2a+c1-2c2+Pk1-2Pk2-2β2-B2P-β1β2.Easy to know,E0, E1, and E2 are the boundary equilibrium points; E∗ is the only interior point equilibrium. Since the real economic significance of the system equilibrium point is taken into account, only the case when the equilibrium point is nonnegative is discussed. According to the assumption that the parameters b,β1,β2,B1,B2, and θ are positive parameters, so parameters need to meet the following conditions when E1,E2, andE∗are greater than zero:(9)P-β1+B1β1a-c1-Pk1>0,P-β2+B2β2a-c2-Pk2>0,W1>0,W2>0.In the following analysis, the nonnegative condition (9) is assumed to be true.


### 3.1. Boundary Equilibrium Point Stability

To discuss the stability of the equilibrium points(x1,x2,I1,I2) of the dynamical system (6), we first calculated the corresponding Jacobian matrix:(10)Jx1,x2,I1,I2=1-α1M1-bB1B2x1α1β1β2-2bθB12x1α1β1-bθB1B2x1α1β1-bB1B2x2α2β1β21-α2M2-bθB1B2x2α2β2-2bθB22x2α2β2β10θ00β20θ,where M1=β1+2bB12β1(θI1+2x1β1)-P-B1β1a-c1-Pk1-bB2(θI2+x2β2), M2=β2+2bB22β2(θI2+2x2β2)-B2β2a-c2-Pk2-bB1(θI1+x1β1)-P.According to the Schur-Cohn stability criterion [12], when all the eigenvalues of the characteristic polynomial corresponding to the Jacobian matrix (10) are in the unit circle on the complex plane, that is, the modulus of any eigenvalue is less than 1, the equilibrium point (x1,x2,I1,I2) is asymptotically stable.Proposition 1.
Boundary equilibrium point is unstable.Proof.
SubstitutingE0=(0,0,0,0) into the general formula (10) of the system Jacobian matrix, the Jacobi matrix of system (6) at the equilibrium point is expressed as follows:(11)JE0=1+α1B1a-c1β1-β1+P1-B1k1β100001+α2+B2a-c2β2-β2+P1-B2k2β200β10θ00β20θ.
We can reach four characteristic roots by calculation, namely,(12)λ1=λ2=θ,λ3=1+α1B1a-c1β1-β1+P1-B1k1β1,λ4=1+α2B2a-c2β2-β2+P1-B2k2β2.
From the nonnegative condition (9) of the equilibrium point and the adjustment coefficient αi>0, the latter two characteristic roots both satisfy λ3,4>1. Therefore, the equilibrium point is unstable.Proposition 2.
E
1 and E2 are unstable equilibrium points.Proof.
Jacobian matrix (10) in the specific form at the boundary equilibrium point is(13)JE1=1+α11-θH11-θB2α1H1β22B1β11-θθα1H1β11-θθB2α1H12B1β101+α2W22B1β100β10θ00β20θ,where H1=1-B1a-c1-Pk1β1-P.
By calculating the eigenvalue of the Jacobian matrixJ(E1),λ1=θ,λ2=1+α2W2/2B1β1,λ3,4=1/2H2±H22-4θ, where H2=1+θ+1-θα11-B1a-c1-Pk1β1-P.
From the previous parameters thatB1,α2>0, λ2>1 can be pushed by the inequality equation (9). Thus, according to discrete dynamical system’s equilibrium point stability determination conditions in the second chapter, we can know that E1 is unstable boundary equilibrium. Similarly, the border equilibrium E2 is also unstable.


### 3.2. Interior Point Equilibrium Stability

This section focuses on the stability of interior point equilibriumE∗. The Jacobian matrix (10) of interior point equilibrium E∗ is expressed as follows:(14)JE∗=1-21-θα1W13B2β2-1-θα1W13B1β121-θθα1W13B2β1β2-1-θθα1W13B1β1β2-1-θα2W23B2β21-21-θα2W23B1β1-1-θθα2W23B2β1β221-θθα2W23B1β1β2β10θ00β20θ.The characteristic polynomial of the matrixJ(E∗) is given as P(λ), and P(λ)=λ4+p1λ3+p2λ2+p3λ+p4, and the polynomial coefficients by the numerical calculation are as follows:(15)p1=-21+θ-1-θ3W2α23B1β1+W1α1B2β2,p2=-1+θW1α121+3θB1β1-1-θW2α2+B231+4θ+θ2B1β1-21+2θ-3θ2W2α2β23B1B2β1β2,p3=-2θ3B2β22θ2-3+θW2α2+31+θB1β1-1-θW1α133+2θB1β1-81-θW2α29B1B2β1β2,p4=θ2.From the Schur-Cohn stability criterion [12], if there is an associated characteristic polynomial P(λ), that is, all eigenvalues of the Jacobian matrix J(E∗) are located in the unit circle on the complex plane, then the coefficients of the polynomial need to satisfy the conditions in the Schur-Cohn stability criteria at the same time. The system parameters are added to it and simplified:(16)P1=1+p1+p2+p3+p4=α1α2θ-12W1W23B1B2β1β2>0,P-1=1-p1+p2-p3+p4=B1B2β1β241+θ31+θ-N1-1-θα1P-β181+θ-N2+21-θB22α221+θ-1-θα1P-β1P-β1β22+B12β12N341+θ-N13B1B2β1β2>0,1±p4=1±θ2>0,where N1=1-θα22P+B2a+c1-2c2+Pk1-2Pk2-2β2,N2=1-θα25P+B2a+4c1-5c2+4Pk1-5Pk2-5β2,N3=1-θα1P-1+B2a-2c1+c2-2Pk1+Pk2β2.From the positive and negative assumptions of the system parameters and inequality (9), we can learn that the stability conditions P(1)>0 and 1±p4>0 are tenable for sure. Therefore, the sufficient conditions for the asymptotic stability of the internal point equalization E∗ of the dynamical system (6) are classified as follows.WhenP(-1)=1-p1+p2-p3+p4>0 and M3±>0, (17)p1+p3<1+p2+p4p2-p2p4+p4-p43-p1p3+p12p4<1-p42+p2p4-p2p42+p1p3p4-p32.All the eigenvalues of J(E∗) satisfy the λi<1(i=1,2,3,4), so that the equilibrium point E∗ is stable when the condition (17) is satisfied. From the above equilibrium point analysis, we can see that the stability of the known equilibrium points is dependent on the value of the system parameters.


## 3.1. Boundary Equilibrium Point Stability

To discuss the stability of the equilibrium points(x1,x2,I1,I2) of the dynamical system (6), we first calculated the corresponding Jacobian matrix:(10)Jx1,x2,I1,I2=1-α1M1-bB1B2x1α1β1β2-2bθB12x1α1β1-bθB1B2x1α1β1-bB1B2x2α2β1β21-α2M2-bθB1B2x2α2β2-2bθB22x2α2β2β10θ00β20θ,where M1=β1+2bB12β1(θI1+2x1β1)-P-B1β1a-c1-Pk1-bB2(θI2+x2β2), M2=β2+2bB22β2(θI2+2x2β2)-B2β2a-c2-Pk2-bB1(θI1+x1β1)-P.According to the Schur-Cohn stability criterion [12], when all the eigenvalues of the characteristic polynomial corresponding to the Jacobian matrix (10) are in the unit circle on the complex plane, that is, the modulus of any eigenvalue is less than 1, the equilibrium point (x1,x2,I1,I2) is asymptotically stable.Proposition 1.
Boundary equilibrium point is unstable.Proof.
SubstitutingE0=(0,0,0,0) into the general formula (10) of the system Jacobian matrix, the Jacobi matrix of system (6) at the equilibrium point is expressed as follows:(11)JE0=1+α1B1a-c1β1-β1+P1-B1k1β100001+α2+B2a-c2β2-β2+P1-B2k2β200β10θ00β20θ.
We can reach four characteristic roots by calculation, namely,(12)λ1=λ2=θ,λ3=1+α1B1a-c1β1-β1+P1-B1k1β1,λ4=1+α2B2a-c2β2-β2+P1-B2k2β2.
From the nonnegative condition (9) of the equilibrium point and the adjustment coefficient αi>0, the latter two characteristic roots both satisfy λ3,4>1. Therefore, the equilibrium point is unstable.Proposition 2.
E
1 and E2 are unstable equilibrium points.Proof.
Jacobian matrix (10) in the specific form at the boundary equilibrium point is(13)JE1=1+α11-θH11-θB2α1H1β22B1β11-θθα1H1β11-θθB2α1H12B1β101+α2W22B1β100β10θ00β20θ,where H1=1-B1a-c1-Pk1β1-P.
By calculating the eigenvalue of the Jacobian matrixJ(E1),λ1=θ,λ2=1+α2W2/2B1β1,λ3,4=1/2H2±H22-4θ, where H2=1+θ+1-θα11-B1a-c1-Pk1β1-P.
From the previous parameters thatB1,α2>0, λ2>1 can be pushed by the inequality equation (9). Thus, according to discrete dynamical system’s equilibrium point stability determination conditions in the second chapter, we can know that E1 is unstable boundary equilibrium. Similarly, the border equilibrium E2 is also unstable.


## 3.2. Interior Point Equilibrium Stability

This section focuses on the stability of interior point equilibriumE∗. The Jacobian matrix (10) of interior point equilibrium E∗ is expressed as follows:(14)JE∗=1-21-θα1W13B2β2-1-θα1W13B1β121-θθα1W13B2β1β2-1-θθα1W13B1β1β2-1-θα2W23B2β21-21-θα2W23B1β1-1-θθα2W23B2β1β221-θθα2W23B1β1β2β10θ00β20θ.The characteristic polynomial of the matrixJ(E∗) is given as P(λ), and P(λ)=λ4+p1λ3+p2λ2+p3λ+p4, and the polynomial coefficients by the numerical calculation are as follows:(15)p1=-21+θ-1-θ3W2α23B1β1+W1α1B2β2,p2=-1+θW1α121+3θB1β1-1-θW2α2+B231+4θ+θ2B1β1-21+2θ-3θ2W2α2β23B1B2β1β2,p3=-2θ3B2β22θ2-3+θW2α2+31+θB1β1-1-θW1α133+2θB1β1-81-θW2α29B1B2β1β2,p4=θ2.From the Schur-Cohn stability criterion [12], if there is an associated characteristic polynomial P(λ), that is, all eigenvalues of the Jacobian matrix J(E∗) are located in the unit circle on the complex plane, then the coefficients of the polynomial need to satisfy the conditions in the Schur-Cohn stability criteria at the same time. The system parameters are added to it and simplified:(16)P1=1+p1+p2+p3+p4=α1α2θ-12W1W23B1B2β1β2>0,P-1=1-p1+p2-p3+p4=B1B2β1β241+θ31+θ-N1-1-θα1P-β181+θ-N2+21-θB22α221+θ-1-θα1P-β1P-β1β22+B12β12N341+θ-N13B1B2β1β2>0,1±p4=1±θ2>0,where N1=1-θα22P+B2a+c1-2c2+Pk1-2Pk2-2β2,N2=1-θα25P+B2a+4c1-5c2+4Pk1-5Pk2-5β2,N3=1-θα1P-1+B2a-2c1+c2-2Pk1+Pk2β2.From the positive and negative assumptions of the system parameters and inequality (9), we can learn that the stability conditions P(1)>0 and 1±p4>0 are tenable for sure. Therefore, the sufficient conditions for the asymptotic stability of the internal point equalization E∗ of the dynamical system (6) are classified as follows.WhenP(-1)=1-p1+p2-p3+p4>0 and M3±>0, (17)p1+p3<1+p2+p4p2-p2p4+p4-p43-p1p3+p12p4<1-p42+p2p4-p2p42+p1p3p4-p32.All the eigenvalues of J(E∗) satisfy the λi<1(i=1,2,3,4), so that the equilibrium point E∗ is stable when the condition (17) is satisfied. From the above equilibrium point analysis, we can see that the stability of the known equilibrium points is dependent on the value of the system parameters.


## 4. Numerical Simulation

In order to intuitively analyze the dynamic behavior of model (6), we used the numerical simulation to visually describe the dynamic evolution process of the discrete dynamical system (6) with the change of the model parameters, focusing on the market price of carbon trading Pand adjusting speed α’s influence on the dynamic nature of the system. Our study assumes that the parameters are taken as follows: a=5,b=1,β1=0.4,β2=0.5,k1=0.4,k2=0.6,B1=0.6,B2=0.8,c1=0.3,c2=0.5,θ=0.35.Fixed enterprise 2’s strategy adjustment speed isα2=0.3. Figure 1 shows the bifurcation diagram of system (6) with enterprise 1 adjusting the speed α1, and the carbon trading market price is P, and P is 1.5 and 2, respectively. It is easy to see that, in Figures 1(a) and 1(b), with the value α1 increase, the system from the internal equalization changes period doubling bifurcation and eventually enters the chaotic state. Figure 1 not only shows the different paths of chaos in the system, but also shows that the larger the market price of carbon trading P is, the earlier the state variable orbit enters into the chaotic state, and the system has the stronger instability. Figure 2 gives a more detailed description of the system orbit’s change, which is a two-dimensional phase diagram corresponding to the values of the orbit in the different states in Figure 1.Figure 1
Bifurcation diagrams for system (6) with respect to the adjustment rate α1.
(a)
P
=
1.5
(b)
P
=
2Figure 2
Phase portraits for Figure1 with various values of P and α1.
(a)
P
=
1.5
(b)
P
=
2Figure3 clearly validates this conclusion by calculating the equilibrium point stability condition (17) through the program and plotting the stable region of the system on the plane (α1,α2). Comparing the two subgraphs of Figure 3, it can be seen that the increase of P leads to the decrease in the stability of system (6), and the dynamic evolution of the system becomes more unstable.Figure 3
The stability region of dynamic system (6) in (α1,α2)-plane.
(a)
P
=
1.5
(b)
P
=
2We selected the adjustment coefficientα1=1.23 of the system in the chaotic state of Figure 1(b). The initial values of the system are (x1(0),x2(0),I1(0),I2(0))=(0.2,0.2,0.2,0.2) and (0.20001,0.2,0.2,0.2), and the evolution diagram of the system state variables with the time is plotted. The partial image is shown in Figure 4. As the number of iterations increases, the same state variable is gradually separated under the influence of two initial conditions and subsequently conducts motion according to the respective orbit; thus, system (6) has sensitive dependence on the initial condition.Figure 4
Sensitivity of the dynamic system (6) to initial conditions when losing stability.
(a)
(b)In order to highlight the effect of the market price of carbon trading on the stability of the orbital motion of the dynamical system (6), the bifurcation results of the system with different adjustment coefficients α1 are depicted in Figure 5. Figure 5 shows the doubling cycle bifurcation phenomenon. From Figure 5, the system becomes more unstable with the increase of P. This suggests that increases in carbon trading’s market prices lead to the increase of carbon stocks on the market, hence, facilitating environmental protectionism mechanisms which are environmentally beneficial. Companies are also actively working for energy-saving and emission reduction. When the price continues to rise, it results in the increase of inventories in the market. However, beyond a certain extent, the carbon stocks create market’s oversupply, hence influencing the stability of the market. Although, from an economic point of view, this phenomenon is in line with the laws of the market, this phenomenon depicts clearly an undesirable infinite loop under the harmonious development and environmental protection as it fails to achieve the desired environmental protection.Figure 5
Bifurcation diagram for dynamic system (6) with respect to the residual rate P.
(a)
α
1
=
1.15
(b)
α
1
=
1.55


## 5. Chaos Control

From the numerical simulation results of the evolution process in system (6), it can be seen that the adjustment speed and the marketing price of carbon trading have a great influence on the stability of the system. The chaotic behavior of the dynamical system (6) is controlled by the delayed feedback control method, so that the system can be restored to the stable motion state.Adding the feedback control itemKx1(t)-x1(t+1) to system (6) and simplifying the system, we obtained the following form of controlled dynamical system:(18)x1t+1=x1t+α11+Kx1tP+β1-bB1B2θI2t+β2x2t+aB1-c1B1-Pk1B1-1-2bB12θI1t+β1x1tx2t+1=x2t+α2x2tP+β2-bB1B2θI1t+β1x1t+aB2-c2B2-Pk2B2-1-2bB22θI2t+β2x2tI1t+1=θI1t+β1x1tI2t+1=θI2t+β2x2t.The stability of the controlled system (18) is analyzed. It is easy to know that the controlled system has the same equilibrium point as the discrete dynamical system (6). The Jacobian matrix corresponding to the controlled system (18) is(19)Jx1,x2,I1,I2=1-α1M11+K-bB1B2x1α1β1β21+K-2bθB12x1α1β11+K-bθB1B2x1α1β1K+1-bB1B2x2α2β1β21-α2M2-bθB1B2x2α2β2-2bθB22x2α2β2β10θ00β20θ.As shown in Figure1(b), when the adjustment coefficient is about 1.21, system (6) enters the chaotic state. The expression of the Jacobian matrix (19) in the interior point equalization E∗ is given under the original parameter value:(20)J=1-1.374531+K-0.687281+K-1.202721+K-0.481091+K-0.12610.7478-0.11034-0.176540.400.35000.500.35.Our study sets the characteristic polynomial of the matrix (20) as P(λ)=λ4+u1λ3+u2λ2+u3λ+u4, and when the coefficients of the matrix satisfy the Schur-Cohn criterion, then(21)1+u2+u4>u1+u31±u4>01-u42+u2u4-u2u42+u1u3u4-u32>u2-u2u4+u4-u43-u1u3+u12u4.The characteristic root of matrix (18) is located within the unit circle on the complex plane, and it can be seen that the internal equilibrium point of the controlled system (18) is stable under the settled parameter values, and the chaotic motion of system (6) can also be adjusted to the expected stable orbit, and it can consequently reach the value range of feedback gain intensity.The effectiveness of control on the chaos of the discrete dynamical system (6) is reached by numerical simulation. Figure 6 shows that if we add the feedback control term Kx1(t)-x1(t+1) into the dynamic iterative equation on the emission reduction of the controlled system (6), the unstable orbit can be controlled as long as the feedback gain intensity K is greater than 0.231. The process that chaotic behaviors of the controlled system under the different gain intensities evolve from the initial value to the stable orbit is expressed in Figure 7. By comparison, the time the state variable reaches stable state in Figure 7(b) is shorter than that in Figure 7(a). The carbon trading market also gets stabilized more quickly.Figure 6
Bifurcation diagram of system (18) with respect to the controlling factor k.Figure 7
The evolution path of system (18) at different gain intensities.
(a)
K
=
0.21
(b)
K
=
0.45


## 6. Conclusion

Based on the hypothesis of participant's bounded rationality, our study formulated a novel Cournot duopoly game model of carbon emission reduction and subsequently analyzed the dynamic adjustment mechanism of emission reduction of the participating enterprises. The main idea in our model is that each firm's decision is to choose its carbon emission in each period according to the marginal profit observed from the previous period. We have established a corresponding dynamics of players' carbon emission adjustment and done a detailed dynamic analysis for it. There are three boundary equilibriums and a unique interior equilibrium in this system. By analyzing the stability of each equilibrium point of the discrete dynamical model system, the conditions of the asymptotic equilibrium of the interior point were obtained according to the Schur-Cohn stability criterion. Furthermore, the numerical simulation model of this paper intuitively shows the dynamic evolution process and the complexity of the system. The simulation results revealed that the parameters influence the dynamic properties of the system. It is observed that the equilibrium of the system may lose stability via different bifurcations, either flip bifurcation or Neimark-Sacker bifurcation. Moreover, when the adjustment speed gets too high, the system is likely to lose the original stability and even show complex chaos easily. Higher carbon trading market prices affect the stability of the system and may lead to system instability. Finally, the delay control method has been proved to effectively control the system in the quest to restore back a stable and an equilibrium market.



---

*Source: 8301630-2017-10-12.xml*