Abstract

This article investigates the dynamics of a mixed triopoly game in which a state-owned public firm competes against two private firms. In this game, the public firm and private firms are considered to be boundedly rational and naive, respectively. Based on both quantity and price competition, the game’s equilibrium points are calculated, and then the local stability of boundary points and the Nash equilibrium points is analyzed. Numerical simulations are presented to display the dynamic behaviors including bifurcation diagrams, maximal Lyapunov exponent, and sensitive dependence on initial conditions. The chaotic behavior of the two models has been stabilized on the Nash equilibrium point by using the delay feedback control method. The thresholds under price and quantity competition are also compared.

1. Introduction

The research on chaos theory and bifurcation theory based on different dynamical systems is the earliest chaotic dynamics branch with the most widespread application [15]. In recent years, mathematical models which describe the oligopoly dynamics of firms’ competition have got much attention due to their important characteristics in analyzing competition patterns and market control. Literature includes many intensive studies on duopoly games (for example, see [69] for Cournot competition, [10, 11] for Bertrand competition, [12, 13] for Stackelberg competition, and [1416] for Cournot–Bertrand case).

Compared with the duopoly game models, the analysis of oligopolies with more than two firms has been less addressed in the literature. In practice, such a game is close to the economic reality and it is widely deployed in oligopolies. However, analyzing the dynamics of such a game is a complex task. Therefore, various investigations were performed previously by considering Cournot triopoly game with three homogeneous players. For example, Puu [17] considered three naive players in the game model, in which an iso-elastic demand function and constant marginal cost were assumed. Agiza et al. [18] proposed Kopel triopoly Cournot game model; they assumed that the competitors did not immediately offer the optimal quantity computed on the basis of the (bounded rationality) profit maximization problem, but that they adjusted the previous quantity in the direction of the computed one. Tu and Wang [19] presented a dynamic master-slave Cournot triopoly game model with homogeneous bounded normal players. They analyzed the effects of changes in the adjustment speed parameters on the system dynamics. Alnowibet et al. [20] modeled a dynamic Cournot triopoly game that was composed of three homogeneous bounded rational players, and the demand function was characterized by log-concavity.

Additionally, Cournot triopoly games considering heterogeneous players were also investigated. Ma and Ji [21] considered a Cournot triopoly game with a square inverse demand. Ji [22] further investigated the game model based on heterogeneous expectations in electric power triopoly. Elabbasy et al. [23] studied the dynamical behavior of the triopoly game with heterogeneous players with linear cost function. The three players were considered to be boundedly rational, adaptive, and naive. Then, Elabbasy et al. [24] further generalized the triopoly game with heterogeneous players with nonlinear cost function. Based on a quantity competition, Askar and Alshamrani [25] studied four different models, namely, cooperative Cournot triopoly, rational Bertrand triopoly, rational Cournot triopoly, and Puu triopoly.

On the other hand, there are few studies on Bertrand triopoly game. For instance, in Sun and Ma [26], a Bertrand triopoly model was developed. They considered the bounded rational expectations, and local stability of the Nash equilibrium was studied. Furthermore, Sun and Ma [27] introduced a Bertrand triopoly model using nonlinear demand functions. Ma and Wu [28] studied the effect of delayed decision on the stability of a Bertrand triopoly model. Zhao [29] constructed a nonlinear dynamic Bertrand triopoly model based on the heterogeneity in expectations, namely, naive expectation and bounded rationality. The instability of the boundary equilibrium point and the stability condition of the internal equilibrium point were obtained.

As exceptions, in a differentiated triopoly model with heterogeneous firms, the local stability of the Nash equilibrium under both quantity and price competition was analyzed in [30, 31]. We can also cite Elsadany et al. [32], who considered an oligopoly game with four heterogeneous firms producing perfect substitute goods.

However, these results mentioned above depend on the assumption that all firms are private and profit-maximizers. Therefore, they may not apply to the increasingly important and popular mixed oligopolies, in which state-owned public firms compete with private firms. Analyses of mixed oligopolies date to Merrill and Schneider [33]. In most countries, there exist state-owned public firms that have substantial influence on their market competitors, such as the airline, steel, insurance, hospital, and banking industries. In the studies of mixed oligopolies, it is generally assumed that a public firm maximizes social surplus, while private firms maximize profits. For examples of mixed oligopolies and recent developments in research in this field, see [34, 35] and the references therein for more details.

This paper proposes a mixed triopoly game with heterogeneous players. Firm 0 is a state-owned public firm adopting bounded rationality, and its payoff is the social surplus. Firm is a private firm, and its payoff is its own profit. We assume the private firms are both naive players. We analyze the local stability of the Nash equilibrium under both quantity and price competition. Numerical simulations are used to provide experimental evidence for the complex behavior of the evolution of the system, including bifurcation diagrams, maximal Lyapunov exponent, and sensitive dependence on initial conditions. A feedback control scheme is adopted to overcome the uncontrollable behavior of the dynamical system occurred due to chaos. To the best of our knowledge, this article is the first to use analytical and numerical tools to consider the dynamics of mixed triopoly game.

The remainder of the paper is organized as follows. The description of the mixed triopoly game with heterogeneous players is described in Section 2. In Section 3, we shall study the existence and local stability of fixed points of the three-dimensional dynamics system under Cournot competition. Dynamical behaviors under some change of control parameters of the game are investigated by numerical simulations, and chaos control of the system is also given. The dynamics under price competition are analyzed in Section 4. Finally, a conclusion is drawn in Section 5.

2. Model Specification

We adopt a standard differentiated oligopoly with a linear demand. The quasi-linear utility function of the representative consumer iswhere is the quantity produced by the public firm, is the quantity produced by the private firm, and is the consumption of outside goods provided competitively (with a unitary price). Parameters and are positive constants and represents the degree of product differentiation; a smaller indicates a larger degree of product differentiation. For , the inverse demand function is given bywhere and denote firm ’s price and quantity, respectively. From (2), the corresponding direct demand function is given by

The marginal production costs are constant. Let denote the state-owned public firm 0’s marginal cost. To simplify the calculations, we assume that two private firms have the same marginal cost with . Furthermore, we assume that the equilibrium quantities of public and private firms are strictly positive under both Bertrand and Cournot competition. Let ; this assumption is satisfied if and only if and .

Since firm 0 is a public firm, its payoff is the social surplus, given by

Firm is a private firm, and its payoff is its own profit: .

3. Cournot Competition

3.1. Local Stability of Equilibrium Points

The first-order conditions for public and private firms are, respectively,

From (5) and (6), we obtain the following reaction functions for public and private firms:

We assume that the public firm 0 uses bounded rationality, hence does not have a complete knowledge of the demand function of the market, and builds its output decision on the basis of the expected marginal payoff . If the marginal payoff is positive (negative), it increases (decreases) its production at the next period output. Then, the dynamical equation of player 0 has the formwhere is a positive parameter which represents the relative speed adjustment. By using equation (5), the dynamical equation of the boundedly rational player is

The private firm is a naive player; he computes his outputs using the reaction function in equation (7); the dynamical equation of player has the form

Then, the dynamical mixed triopoly game with partial heterogeneous players is given by

To study the qualitative behavior of the solutions of the discrete dynamical system (11), we define the equilibrium points of the mixed triopoly game as a nonnegative fixed point of the dynamical system (11). By setting in (11), it is concluded that the dynamical system has two equilibrium points:where

Note that is the boundary equilibrium point and is the Nash equilibrium point. The local stability of the equilibrium points depends on the eigenvalues of the Jacobian matrix given by

Theorem 1. The boundary equilibrium point is a saddle point.

Proof. At the boundary equilibrium point , the Jacobian matrix takes the formIt is easy to calculate that the three eigenvalues are , , and . By the condition , we conclude that . On the other hand, implies that . Then, is a saddle point of the dynamical system.
Next, we discuss the local stability of the Nash equilibrium point. The Jacobian matrix evaluated at takes the formwhere the constant is determined byIt is known that is asymptotically stable if and only if all the roots of the characteristic equationhave magnitudes of eigenvalues less than one, in whichFrom Schur–Cohn stability criterion, the coefficients of the characteristic polynomial need to satisfy the following conditions [29, 30]:

Theorem 2. The Nash equilibrium is locally asymptotically stable provided thatwhere is defined as in equation (17).

Proof. For , it can be calculated thatFurthermore, the inequality is obvious by noting that . On the other hand, after some modifications, we obtainIt concludes that the second inequality in (20) holds ifFinally, some tedious rearrangements lead toIt is easy to see that is a quadratic equation in . When , for all . In the case of , the discriminant of is positive; then, has two real solutions and . We can easily obtain that one root and the other root . From these results, the threshold of ensuring the local asymptotic stability of is given by .
From Theorem 2, we can see that the local stability of Nash equilibrium point largely depends on the adjustment speed of the public firm 0’s quantity decision. If the value of is high (company is overresponsive), the Nash equilibrium will lose stability and dynamical behavior may occur.

3.2. Numerical Simulations

In this section, the dynamical evolution of the nonlinear dynamical system (11) is numerically simulated, and the dynamical evolution process is visualized by simulation results such as stable region, bifurcation diagram, maximum Lyapunov exponent, strange attractor, and sensitive dependence on initial conditions.

We fixed the parameter values as , , and . The stability region of Nash equilibrium in the -plane is shown in Figure 1(a). It shows that the stability region increases (decreases) as increases when is smaller (larger) than 0.35. When the parameter takes values out of the stability region, the equilibrium point loses its stability due to flip bifurcation, which is presented in Figure 1(b).

Figure 2(a) depicts the characteristic behavior of a period-doubling bifurcation, and one can get that the Nash equilibrium point is locally stable provided that . Period-2 cycle occurs at the market when as indicated by the two branches. After that, both branches split simultaneously, yielding a period-4 cycle. A cascade of further period doubling occurs when the first firm increases , yielding period-8, period-16, and so on. The equilibrium undergoes a flip bifurcation at , and the system exhibits chaotic behavior while approaches 1.737.

Figure 2(b) shows the corresponding maximum Lyapunov exponents of system (11), which is a good indicator for bifurcation and chaos. It is observed that the period-doubling bifurcation arises as reaches the value of , and then the system turns into chaos as increases.

Three-dimensional phase portraits for different values of are shown in Figure 3, which can give a more detailed description for the orbits of system (11). It shows a flip bifurcation process, where chaos occurs in the end of periods , and the strange attractors are shown in the fourth figure.

Figure 4 demonstrates the sensitivity of system (11) to initial conditions, which is one of the main characteristics of chaotic behavior. Figure 4 plots two orbits initially from the slightly deviated points and , respectively. The blue curves represent outputs of the public firm and the red ones represent outputs of the private firm. It is seen clearly that even if the initial production of the public firm alters a little, great impact will emerge in all firms after a series of iterations.

3.3. Chaos Control

The occurrence of chaotic dynamics in the triopoly game is unacceptable and the firms are hoping to find some methods to control chaos to avoid the complexity. Many methods have been used to control chaos in oligopoly games. In the papers [24, 28, 29], it has been presented how the delay feedback control (DFC) method can be applied to control chaos in different economic models. We apply this technique to control the chaotic behavior for the present triopoly game.

We modify the first equation of system (11) by adding the control action , where is the control parameter. Then, the controlled system has the following form:

The Jacobian matrix of controlled system (26) is as follows:

At the parameter values (k, δ, a0, a1, β) = (2, 0.4, 2, 2.5, 0.5), the Jacobian matrix (27) has the form

Applying Schur–Cohn conditions, matrix (28) has eigenvalues with an absolute less than one when . It means that the system is stable around the Nash equilibrium point under the condition of . From Figure 5, one can see that with the increase of the control parameter , system (26) gradually tends to be stable from chaos and quasi-period. When , system (26) finally evolves to be stable.

Figure 6 depicts the behavior of controlled system (26) for and , respectively, with initial values . From the two graphs, we can find that as the feedback strength value increases, the chaotic behavior can be quickly controlled to stable orbit. Therefore, the above control method is able to control chaos in the game when is large enough, and the market game can switch from chaotic trajectories to an equilibrium state.

4. Bertrand Competition

4.1. Local Stability of Equilibrium Points

In this section, we discuss the Bertrand model in which all firms choose prices. From equation (3), the first-order conditions for public and private firms are, respectively,

From the first-order conditions, we obtain the following reaction functions:

Now, suppose that the decision rule of the public firm 0 is based on marginal profit, and the decision rule of the private firms follows the expected value of naive. Therefore, how to adjust the price of firm 0 mainly depends on whether the current marginal profit is positive or negative in the -th period. Suppose that the public firm adjusts the decision mechanism over time as follows:where indicates the speed at which firm 0 adjusts the price based on marginal profit. Suppose that both private firms adopt the naive expected value decision rule to adjust the product price; it can be expressed in the following form:

From equations (31) and (32), we can obtain the following three-dimensional dynamical system based on price competition:

By letting in (18) for , we have

By solving (34), we get two equilibrium points for dynamical system (33) as follows:where , .

Now, we shall discuss the dynamical behavior of system (33) around and . As discussed in Section 3.1, we return to the three-dimensional map (33). Let be the Jacobian matrix of system (33) corresponding to the state variables ; then, we have

Next, the stability of the boundary equilibrium point of the complex dynamical system (33) will be given in Theorem 3.

Theorem 3. The boundary equilibrium point is a saddle point.

Proof. At the boundary equilibrium point , the Jacobian matrix takes the formThe three eigenvalues can be calculated as , , and . By the condition , we conclude that . On the other hand, implies that .
Now we consider the stability properties of . The Jacobian matrix at takes the formwhere , . The necessary and sufficient conditions for to be asymptotically stable are all eigenvalues of the Jacobian matrix located in the unit circle of the complex plane. The characteristic equation of is as follows:where the coefficients are , , .
Again, according to Schur–Cohn stability criterion, we need to verify the four inequalities discussed as in (20). First, for and , we obtainThe inequality is obvious. After some careful calculations, it is deduced that the inequalityholds if public firm 0’s adjustment speed satisfiesFinally, we consider a quadratic equation in defined byFor , it is not hard to see that the discriminant of is positive and has two roots with andComparing the thresholds and implies that . Therefore, we eventually obtain the following result.

Theorem 4. The Nash equilibrium is locally asymptotically stable provided thatwhere the constant is determined by

As expected, the local stability of the Nash equilibrium in Bertrand game is determined by the adjustment speed of gradient player. If the value of is high (company is overresponsive), the complex dynamical behavior may occur. Moreover, different from the quantity competition case, the threshold of the adjustment speed depends on an additional parameter . The greater the value of is, the higher the stability of the Nash equilibrium is (see Figure 7 for details).

4.2. Numerical Simulations

In this section, we provide some numerical evidence for the chaotic behavior of the dynamical game (33). To study the local stability properties of the Nash equilibrium point, it is convenient to consider the following set of parameters: , , and .

According to equation (42), we know that the stability region in the plane changes as varies. For fixed , the stability region becomes small with increasing. Moreover, it is deduced that the greater the value of is, the larger the stability of the Nash equilibrium is, regardless of the value of the parameter .

Figures 8(a) and 8(b) display the bifurcation diagram with respect to and the corresponding maximum Lyapunov exponents of system (33), respectively. For values of lower than 0.182, the Nash equilibrium point is locally stable. As increases, system (33) becomes unstable and complex dynamical behavior occurs, including period bifurcations and chaos. This means that for large value of speed of adjustment of boundedly rational player, dynamical game (33) converges always to complex dynamics.

A more specific representation is shown in Figure 9. Figure 9 plots the price evolution with time when the system is in the stable region, 2,4-period bifurcations, and chaotic state, respectively. When , the price of all firms will be stable in equilibrium state gradually after a series of fluctuations. When , the complex behavior can be observed and the system is in chaos.

The three-dimensional chaotic attractors are given in Figure 10, which exhibit fractal structure. When all parameters are kept fixed and only the adjustment speed of public firm varies, one can see that the structure of the triopoly game becomes complicated through period-doubling bifurcations. In addition, more complex bounded attractors are created around Nash equilibrium point, which are periodic cycles of high order or chaotic attractors.

Figures 11(a) and 11(b) show the two orbits of public firm’s price and private firm’s price, respectively, where the blue ones start from the initial point and the red ones start from the initial point . From the two figures, it is clear that the time series are indistinguishable at the beginning. However, after a number of iterations, the difference between them builds up rapidly.

The comparisons between the thresholds and are depicted in Figure 12. It is shown that the value of the adjustment speed where the Nash equilibrium loses its stability is lower under Bertrand competition than under Cournot competition for small . More precisely, for all if (see Figure 12(a)). As stated in [30], it can be understood as the consequences of the stability of a fiercer competition under Bertrand competition than under Cournot competition regardless of the degree of product differentiation for smaller . It is more complicated in the case of ; Figure 12(b) shows us the relationship between the two critical values for .

4.3. Chaos Control

As discussed in Section 3.3, by adding the control action to system (33), the controlled system has the following form:

The Jacobian matrix of the controlled system (47) is as follows:

At the parameter values (ρ, δ, α, c0, c1, β) = (0.27, 0.4, 5.5, 3.5, 3, 0.5), the Jacobian matrix (48) takes the form

Again, using Schur–Cohn conditions, the matrix (49) has eigenvalues with an absolute less than one when . From Figure 13(a), one can see that the system is controlled from a chaotic state to a stable state when . Figure 13(b) shows the behavior of the controlled system for , starting from initial values . One can conclude that a controlled behavior converges to the fixed point when is large enough.

5. Conclusion

In this paper, we study the dynamics in a triopoly game with product differentiation, in which a state-owned public firm competes with two private firms. We investigate the local stability of the equilibrium points under both quantity and price competition. Numerical simulations of the two kinds of games have been analyzed by means of period-doubling bifurcations, strange attractors, maximal Lyapunov exponent, and sensitive dependence on initial conditions. We have stabilized the chaotic behavior of the models to stable fixed points by the delay feedback control method. We find that the boundedly rational player accelerates the adjustment speed of the output quantity (price), and it leads to the instability of the system and makes the system go to a chaotic region. The comparison between Cournot and Bertrand competition implies that the Nash equilibrium is more stable in a quantity setting for smaller .

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This study was supported by the Humanities and Social Science Project of Ministry of Education (grant no. 20YJA910001) and the Foundation of Educational Department of Liaoning Province (W201783664).