# Chaos Control on a Duopoly Game with Homogeneous Strategy

Authors: Manying Bai; Yazhou Gao
Journal: Discrete Dynamics in Nature and Society
 (2016)

---

## Abstract

We study the dynamics of a nonlinear discrete-time duopoly game, where the players have homogenous knowledge on the market demand and decide their outputs based on adaptive expectation. The Nash equilibrium and its local stability are investigated. The numerical simulation results show that the model may exhibit chaotic phenomena. Quasiperiodicity is also found by setting the parameters at specific values. The system can be stabilized to a stable state by using delayed feedback control method. The discussion of control strategy shows that the effect of both firms taking control method is better than that of single firm taking control method.


---

## Body

## 1. Introduction

An oligopoly is a market form in which a market or industry is dominated by a small number of sellers. Thus, they have power to decide the market price and the yields. The behavior of one player inevitable influences the other players’ behavior. To maximize their profits, each oligopolist is likely to be aware of the actions of the others.Cournot, in 1838, introduced the first formal theory of oligopoly [1]. In Cournot model, each oligopolist assumes that other firms hold their outputs constant, and all oligopolists select a quantity based on others’ outputs to maximize profits. However, oligopoly game becomes complex when the players adopt dynamic strategy depending on their previous outputs and their rivals’ outputs.Expectations play an important role in modeling economic phenomena. A firm can choose its expectation rules to adjust its production. There exist three different firms’ expectations: naive, bounded rational, and adaptive [2]. The case of where both players adopt bounded rational expectation and that of one takes bounded rational expectation and the other takes adaptive expectation have been studied by many researchers [3–10]. However, there are few literatures that focus on the case of where both players adopt adaptive expectation.Recently, the dynamics of duopoly game has been studied by many researchers. Chaotic phenomena in duopoly game were first found by Puu [11]. Later, Kopel [12] studied the stability of the duopoly game with different demand functions and different cost functions. Chaos was also found in this paper. Agiza and Elsadany [2] studied the dynamics in the Cournot duopoly game with heterogeneous players; in which case, one player accepts naive expectation, and the other adopts bounded rational expectation. Ma and Ji [4] studied duopoly game with homogenous players in electric power industry, where both players adjust their outputs according to bounded rational rule. Sarafopoulos [5] assumed that one player accepts bounded rational expectation and the other uses adaptive expectation.Since chaotic phenomena were found in duopoly game, there are also many researchers focusing on the chaos control in duopoly game [6, 13–17]. However, to our best knowledge, there are few literatures that focus on the efficiency of chaos control. In this paper, we explore the dynamics of a homogeneous duopoly game where both players take adaptive expectation. According to the analysis of numerical simulations, the parameters’ effects on the stability of the system are obtained. To lead the system to stable state, we try to use DFC method. The efficiency of the case where only one company takes control measure and that of the case where both companies take control measure is compared.The paper is organized as follows. In Section2, the model and the parameters are introduced. In Section 3, we study the equilibriums and the stability of the equilibriums of the model. In Section 4, we give the numerical results about the reactions of the two players. In Section 5, we control the system to a stable situation. Section 6 gives the conclusions.


## 2. The Model

We assume that the duopoly players produce homogenous goods which are perfect substitutes and offer goods at discrete-time periods. Therefore, the duopoly players face the same market demands. They both choose adaptive expectation rule to decide the amounts of the goods in the next period as their response strategy. Assume that the total demand is reciprocal to pricep; therefore, the reverse demand function is (1)
p
=
1
q
1
+
q
2
,where q
i
(
i
=
1,2
) represents the quantity that firm i produced. The cost of firm i contains two parts: one is constant cost d
i and the other is variable cost. The per units cost corresponding to the variable cost is constant, and the value is c
i. Then, the total cost of firm i is(2)
C
i
q
i
=
c
i
q
i
+
d
i
.Now, we present a simple case which will be used to modify players’ expectation latter. In this case, both firms have the conception that the other firm will produce at time t
+
1 as it did at time t. Then, they will decide their outputs based on the amounts their rival player did at time t. Then, the expectation net profits of firm 1 and firm 2 at time t
+
1 can be expressed as(3)
π
1
=
p
∗
q
1
-
C
1
q
1
=
q
1
t
+
1
q
1
t
+
1
+
q
2
t
-
c
1
q
1
t
+
1
-
d
1
,
π
2
=
p
∗
q
2
-
C
2
q
2
=
q
2
t
+
2
q
1
t
+
q
2
t
+
1
-
c
2
q
2
t
+
1
-
d
2
.The marginal expectation profits of firm 1 and firm 2 at (q
1
,
q
2) are as follows:(4)
∂
π
1
∂
q
1
=
q
2
t
q
1
t
+
1
+
q
2
t
2
-
c
1
,
∂
π
2
∂
q
2
=
q
1
t
q
1
t
+
q
2
t
+
1
2
-
c
2
.The firms can make maximize profits when the marginal profits are zero. Then, the reaction outputs of the two firms with respect to their competitor’s last outputs are(5)
q
1
t
+
1
=
q
1
q
2
=
q
2
t
c
1
-
q
2
t
,
q
2
t
+
1
=
q
2
q
1
=
q
1
t
c
2
-
q
1
t
.In order to make sure that q
i
(
i
=
1,2
) is nonnegative, the reaction functions are modified as (6)
q
1
t
+
1
=
q
2
c
1
-
q
2
,
0
<
q
2
≤
1
c
1
,
0
,
q
2
>
1
c
1
,
(7)
q
2
t
+
1
=
q
1
c
2
-
q
1
,
0
<
q
1
≤
1
c
2
,
0
,
q
1
>
1
c
2
.On the condition that the two firms have adaptive expectation, the firms compute their outputs with weights between their own last outputs and their reaction outputs q
i
q
j
(
i
≠
j
). Then, the outputs they produce in the next period are as follows:(8)
q
1
t
+
1
=
1
-
v
1
∗
q
1
t
+
v
1
∗
q
1
q
2
,
q
2
t
+
1
=
1
-
v
2
∗
q
2
t
+
v
2
∗
q
2
q
1
,where q
1
q
2 and q
2
q
1 are as (6) and (7), respectively, and v
1 and v
2 are weights on firms’ reaction outputs, respectively. We focus on the dynamics of system (8) in the next section.


## 3. Equilibrium and Stability

By settingq
i
t
+
1
=
q
i
t
,
(
i
=
1,2
) in map (8), the equilibrium output points of the dynamic duopoly game can be obtained as the nonnegative solutions of the algebraic system:(9)
-
q
1
t
+
q
2
t
c
1
-
q
2
t
=
0
,
-
q
2
t
+
q
1
t
c
2
-
q
1
t
=
0
.The system has two equilibriums: one is E
1
∗
(
0,0
) and the other is E
2
∗
q
1
∗
,
q
2
∗
=
(
c
2
/
c
1
+
c
2
2
,
c
1
/
c
1
+
c
2
2
). E
1
∗ has no practical significance because both the outputs of two firms are zero. Hence, we just investigate Nash equilibrium E
2
∗.The Jacobian matrix of the two-dimensional map (8) at equilibrium E
2
∗ is(10)
J
E
=
1
-
v
1
1
2
v
1
1
c
1
q
2
∗
-
1
1
2
v
2
1
c
2
q
1
∗
-
1
1
-
v
2
.According to Agiza and Elsadany [3], equilibrium E
2
∗ is locally stable if the following conditions are held: (11)
1
-
T
+
D
>
0
,
1
+
T
+
D
>
0
,
1
-
D
>
0
,where(12)
T
=
2
-
v
1
-
v
2
,
D
=
v
1
v
2
-
v
1
-
v
2
-
1
4
v
1
v
2
1
c
1
q
2
∗
c
2
q
1
∗
+
1
2
v
1
1
c
1
q
1
∗
+
1
2
v
2
1
c
2
q
1
∗
=
v
1
v
2
-
v
1
-
v
2
-
1
4
v
1
v
2
c
1
+
c
2
2
c
1
c
2
+
1
2
v
1
c
1
+
c
2
c
1
+
1
2
v
2
c
1
+
c
2
c
2
.By substituting T and D into (11), the stable conditions of equilibrium E
2
∗ become(13)
2
c
1
+
c
2
v
1
c
2
+
v
2
c
1
-
v
1
v
2
c
1
-
c
2
2
-
4
c
1
c
2
>
0
,
2
c
1
+
c
2
v
1
c
2
+
v
2
c
1
-
v
1
v
2
c
1
-
c
2
2
-
8
c
1
c
2
v
1
+
v
2
+
12
c
1
c
2
>
0
,
-
2
c
1
+
c
2
v
1
c
2
+
v
2
c
1
+
v
1
v
2
c
1
-
c
2
2
+
4
c
1
c
2
v
1
+
v
2
+
1
>
0
.In view of (10), the eigenvalues associated with the equilibrium E
2
∗
(
q
1
∗
,
q
2
∗
) are (14)
λ
1
=
T
+
T
2
-
4
D
2
,
λ
2
=
T
-
T
2
-
4
D
2
.Now, suppose that (15)
T
2
-
4
D
<
0
,
λ
i
=
1
,
i
=
1,2
.Then, D
=
1; namely,(16)
v
1
v
2
-
v
1
-
v
2
-
1
4
v
1
v
2
c
1
+
c
2
2
c
1
c
2
+
1
2
v
1
c
1
+
c
2
c
1
+
1
2
v
2
c
1
+
c
2
c
2
=
1
.It follows from (16) that (17)
v
1
=
1
+
v
2
-
v
2
c
1
+
c
2
/
2
c
2
v
2
-
1
-
v
2
c
1
+
c
2
2
/
4
c
2
+
c
1
+
c
2
/
2
c
2
=
v
1
∗
.Then, q
1
∗
,
q
2
∗
,
v
1
∗ is a candidate for Neimark-Sacker bifurcation point [18].


## 4. Numerical Simulations

In this section, numerical simulation results corresponding to model (8) are presented. Bifurcation diagrams and phase portraits with respect to different parameters are used to show complex dynamical behaviors of the duopoly model.Figure1 presents the bifurcation diagram of system (8) with respect to parameter c
1 (per unit cost of firm 1) against variable q
1 for v
1
=
v
2
=
0.5
,
c
2
=
2. It is seen that the system is in periodic state for c
1
>
0.18. As c
1 decreases, periodic motion and chaotic motion occur alternatively. And the system is driven to chaos through quasiperiodic route. The main region where the system appears chaotic behavior is the range with c
1
∈
[
0.082
,
0.121
]. There exist many period orbits, such as period 7 orbit with c
1
∈
[
0.051
,
0.066
] and period 8 orbit with c
1
∈
[
0.014
,
0.017
], in model (8). For c
1
>
0.18, the stable output of firm 1 decreases with per unit cost c
1 increasing, which is consistent with the fact that the comparative superiority of firm 1 decreases as per unit cost c
1 increasing.Figure 1
Bifurcation diagram with respect to parameterc
1 against variable q
1 with 500 iterations of map (8) for v
1
=
v
2
=
0.5, c
2
=
2.Now, setv
2
=
0.7, c
1
=
0.14, and c
2
=
2. It follows from (17) that v
1
∗
≈
0.2958. Figure 2 shows the bifurcation diagram of map (8) with respect to parameter v
1. At v
1
∗
≈
0.2958, a Neimark-Sacker bifurcation occurs. Nash equilibrium E
2
∗ is locally stable for v
1
<
0.2985 and losses its stability for v
1
>
0.2985. The system evolves from Nash equilibrium E
2
∗ into chaos with parameter v
1 (weight on firm 1’s reaction function) increasing, through the mechanism of quasiperiodic route.Figure 2
Bifurcation diagram with respect to parameterv
1 against variable q
1 with 500 iterations of map (8) for v
2
=
0.7
,
c
1
=
0.14
,
c
2
=
2.Figure3 depicts the trajectories of outputs of firm 1 and firm 2 in the phase space (
q
1
,
q
2
) for v
1
=
0.5, v
2
=
0.5, c
2
=
2, and the initial point (0.4, 0.4). The red point is Nash equilibrium. A periodic attractor with c
1
=
0.13 is shown in Figure 3(b), while a quasiperiodic attractor with c
1
=
0.12 is shown in Figure 3(c). For c
1
=
0.19, the outputs of firm 1 and firm 2 converge to the Nash equilibrium, which is shown in Figure 3(a).Figure 3
(a) Phase portrait of map (8) for v
1
=
0.5
,
v
2
=
0.5
,
c
1
=
0.19
,
c
2
=
2. (b) Phase portrait of map (8) for v
1
=
0.5, v
2
=
0.5
,
c
1
=
0.13
,
c
2
=
2. (c) Phase portrait of map (8) for v
1
=
0.5
,
v
2
=
0.5
,
c
1
=
0.12
,
c
2
=
2.
(a)
(b)
(c)Figure4 depicts the trajectories of outputs of firm 1 and firm 2 in the phase space (
q
1
,
q
2
) for v
2
=
0.7, c
1
=
0.14, c
2
=
2, and the initial point (0.4, 0.4). The red point is Nash equilibrium. For v
1
=
0.3, the outputs of firm 1 and firm 2 converge to the Nash equilibrium, which is shown in Figure 4(a). Figure 4(b) shows a quasiperiodic attractor with v
1
=
0.37 and Figure 4(c) shows a periodic attractor with v
1
=
0.48.Figure 4
(a) Phase portrait of map (8) for v
1
=
0.3
,
v
2
=
0.7
,
c
1
=
0.14
,
c
2
=
2. (b) Phase portrait of map (8) for v
1
=
0.37
,
v
2
=
0.7
,
c
1
=
0.14, c
2
=
2. (c) Phase portrait of map (8) for v
1
=
0.48
,
v
2
=
0.7
,
c
1
=
0.14
,
c
2
=
2.
(a)
(b)
(c)


## 5. Chaos Control in the Dynamic Output System

Any firm does not want to adjust its outputs too frequently and largely. So it is important for firms to control their outputs to a stable process. Chaos control has been studied by many researchers since chaos has been found in economy [6, 13–17]. In DFC (delayed feedback control) method [19], the most common control function is (18)
u
t
=
K
q
t
-
q
t
-
1
.We consider the case that only firm 1 adopts control function to make the system stable, while firm 2 does not realize that it will generate chaos in the future on the condition that firm 1 does not take any measure to prevent this case. Then, system (8) can be rewritten as (19)
q
1
t
+
1
=
1
-
v
1
∗
q
1
t
+
v
1
∗
q
1
q
2
+
u
1
,
t
,
q
2
t
+
1
=
1
-
v
2
∗
q
2
t
+
v
2
∗
q
2
q
1
,where u
1
,
t
=
K
1
(
q
1
t
-
q
1
t
-
1
).To examine the effects of the control function, we compare the output process of firm 1 before and after adding the control function. Setv
1
=
0.85, v
2
=
0.7, c
1
=
0.2, and c
2
=
2 in model (8), which may exhibit chaos. The bifurcation diagram with respect to the control parameter K
1 is given in Figure 5. The system is stable for K
1
<
-
0.22. For K
1
=
-
0.3, the outputs processes of firm 1 before and after using control function are shown in Figure 6. It is seen that the control function effectively leads the outputs to a stable state.Figure 5
Bifurcation diagram with respect to parameterK
1 against variable q
1 with 500 iterations of map (19) for v
1
=
0.85
,
v
2
=
0.7
,
c
1
=
0.2
,
c
2
=
2.Figure 6
(a) Time series of output of firm 1 before using control function. (b) Time series of output of firm 1 after only firm 1 uses control function. (c) Time series of output of firm 1 after both firms use control function.
(a)
(b)
(c)Now, consider the case that both firms realize that they will generate chaos in the future if they do not take any measure to control it. Therefore, both firms use control function to stabilize their outputs. Then, system (8) can be rewritten as (20)
q
1
t
+
1
=
1
-
v
1
∗
q
1
t
+
v
1
∗
q
1
q
2
+
u
1
,
t
,
q
2
t
+
1
=
1
-
v
2
∗
q
2
t
+
v
2
∗
q
2
q
1
+
u
2
,
t
,where u
1
,
t
=
K
1
(
q
1
t
-
q
1
t
-
1
) and u
2
,
t
=
K
2
(
q
2
t
-
q
2
t
-
1
). The outputs of firm 1 are shown in Figure 6(c) for K
1
=
-
0.3 and K
2
=
-
0.3. By comparing Figure 6(b) with Figure 6(c), it is convenient to conclude that the effect of both firms taking control method is better than that of single firm taking control method.


## 6. Conclusion

In this paper, we have investigated a Cournot duopoly model where both the players decide their outputs weighting on their own previous outputs and the optimal outputs on the condition that their rival produces as their previous step. The Nash equilibrium and its local stability were analyzed. The numerical simulations show that the changes of marginal costc
1 and weight factor v
1 may lead the Nash equilibrium to be unstable and the system into chaotic state. The system can quickly arrive at the Nash equilibrium by taking DFC method with a suitable controlling parameter. The effect of both firms taking control method is better than that of single firm taking control method.In view of the fact that firms may not produce any goods if they suffer losses, it is needed to take firms’ profits into consideration in modeling the outputs of firms. This problem will be investigated in our future research.



---

*Source: 7418252-2016-07-12.xml*