Abstract

Two different time delay structures for the dynamical Cournot game with two heterogeneous players are considered in this paper, in which a player is assumed to make decision via his marginal profit with time delay and another is assumed to adjust strategy according to the delayed price. The dynamics of both players output adjustments are analyzed and simulated. The time delay for the marginal profit has more influence on the dynamical behaviors of the system while the market price delay has less effect, and an intermediate level of the delay weight for the marginal profit can expand the stability region and thus promote the system stability. It is also shown that the system may lose stability due to either a period-doubling bifurcation or a Neimark-Sacker bifurcation. Numerical simulations show that the chaotic behaviors can be stabilized by the time-delayed feedback control, and the two different delays play different roles on the system controllability: the delay of the marginal profit has more influence on the system control than the delay of the market price.

1. Introduction

A monopoly market is the case where a trade is completely controlled by a small number of firms. The few firms produce the same or homogeneous products and they must take into account all the information in the market and the actions of the competitors. Cournot [1] first introduced the game model which gives a mathematical description of the dynamic competition in a duopoly market. In a classic Cournot model, each firm is assumed to have naïve expectation and thus guesses that the opponent’s output remains at the same level as in the previous period and then chooses a optimal production strategy in the current period. The classic research has paid much attention to the stability and the complex phenomena in the dynamics of a Cournot game played by players with naïve expectation (e.g., [28]).

Rather than the naïve expectation, a so-called boundedly rationality based on players’ marginal profits has received great attention in recent years. Bischi and Naimzada [9] first gave the general formula of bounded rationality in a duopoly game, where each firm uses the local knowledge of its marginal profit to adjust strategy in a new period: increasing its output if it perceives a positive marginal profit and decreasing its production if the perceived marginal profit is negative. A great deal of research work has been done for the dynamical Cournot game played by homogeneous players with such a kind of rationality (e.g., [912]).

Besides the work on the models with homogeneous players, there is another branch of literature that is interested in the games with heterogeneous players. The dynamic duopoly model with one boundedly rational player and one naïve player has been studied by Agiza and Elsadany [13]. By modifying the cost function in [13], Zhang et al. [14] studied the case with nonlinear cost function. A duopoly game with heterogeneous players and isoelastic demand function has been studied by Tramontana [15] and Angelini et al. [16]. Ma and Xie [17] considered the changing demand to develop the dynamic game models for the two scenarios and analyzed the model’s dynamic behavior. In the model with two heterogeneous players, Fan et al. [18] supposed that the naïve player does not know the rival’s output and adjusts his production according to the previous market price. Ding et al. [19] studied a two-team game played by one team consisting of two boundedly rational players and one team consisting of a naïve player.

In the work on the dynamic Cournot model, time delay has also attracted the attention of many scholars. Many researches examine the effect of delay on dynamics. One-time delay, two-time delay, continuously distributed time delay, and geometric delay are systematically reviewed and studied in [2022], which offer results on existence, stability, and local bifurcations of the equilibrium points. Besides, the work in [2328] studied the dynamical Cournot game with delayed bounded rationality and found that time delay can increase the system stability and delay the occurrence of complex behaviors. Ma and Si [29] also confirmed that the stability of the system depends on the delay and weight. The result showed that the stability of price was closely related to the parameters and the reasonable region of price benefited the firm profit. Further, Ma and Wu [30] studied the influence of parameter change on the market stability and chaos on the sensitivity level. By comparing parameter basin plots, the influences of decision parameters on the dynamics behaviors of the two models are further analyzed [31]. In [2325, 27, 28] time delay is done for the output variable so that decision makers’ expectations are based on the delayed data of their outputs, which implicitly means that any player must know each opponent’s output in the previous periods. These models may not be suitable for the market with incomplete information where a decision maker has too limited knowledge to get the opponents’ output data [26]. It is also pointed out [26] that any producer is able to know his own marginal profit, which is the profit of the last unit production and can be observed in its accounting. Then Ding et al. [26] considered a case that each player makes adjustment strategy according to the delayed marginal profit, the one that averages his own previous marginal profits with different weights. Besides the information of marginal profit, in this work we will consider another fact that the market price history is usually a common knowledge to every player. In the market, there may be players who simply make decision according to an expected price, for example, an average from several previous periods.

In this work, we reconsider the duopoly model in Ding et al. [26]. But different from [26] that only considers homogeneous players and only discusses the time delay structure built for the marginal profit, this paper revises the model and assumes that in the duopoly market the two players are heterogeneous and they consider different time delay structures. One player is still assumed to set delay on the marginal profit and thus adjusts strategy by weighting his marginal profits in the previous periods, as done in [26]. But another player is assumed to make decision according to an expected price that averages the price history with proper weights. For this kind of duopoly game with different delay structures, this paper establishes the dynamical system and discusses its dynamic properties form theoretical and numerical approaches.

2. The Model

2.1. Two Different Delay Structures

We assume that there are two firms producing homogeneous goods for sale and the total supply determines the market price through an inverse demand function . If denotes the cost function of firm , then the profit of firm will be given by

Suppose that firm 1 adjusts its output with the same time delay structure as discussed in [26]. That is, firm 1 sets delay on the marginal profit , and by weighting the marginal profit history firm 1 updates its strategy in period according to the following equation (see [26]):where is a weight coefficient to the delayed period , , and is a positive constant representing the relative adjustment rate. As done in [26], we also consider a delay structure with one step: , then (3) can be written aswhere is the weight coefficient assigned to the nondelayed period and hence is the one to the delayed period .

Next, we consider firm 2 that is supposed to have a time delay structure built for the price history. In a real market, the market price is usually a common knowledge to all players so that the information of the price history may be perceived by all of them. In our model, we suppose that in period firm 2 makes an expectation for the market price, which averages the above price history with different weight coefficients and hence takes its form aswhere is the weight coefficient and . For player 2, we also consider the case of one step delay; that is,where is the weight to the nondelayed period and is associated with the delayed period . From such a time delay structure, firm 2 may expect a profit for its output choice :

It is supposed that in period firm 2 will make its decision in accordance with a best reply to the expected profit . That is, to maximize its expected profit firm 2 will choose its output strategy according to the following formula:

2.2. The Dynamical System

For a duopoly game with two homogeneous players, Section 2.1 has built two different time delay structures and obtained (4) and (9), which describe two-dimensional and discrete dynamics:

To give a specific form for dynamics (10), in our model we also assume linear inverse demand function [5, 9, 10, 13, 14, 18, 23, 24, 26, 3234]:and a quadratic form of cost function [10, 14, 24, 26]:where , , and are all positive. Then the profit of firm is given byfrom which we obtain the marginal profit of firm 1 in period :And from firm ’s expected price in (7), the market price expression in (11), and the cost function in (12), we obtain firm ’s expected profit as follows:which gives a maximization solution

Taking (14) and (16) into (10), we obtain a dynamical system with one step delay:

To study dynamics (17) conveniently, we write for and for ; then we rewrite (17) as a four-dimensional system as follows:

3. Stability of the Equilibriums

Letting () in system (18), we get two nonnegative equilibrium points: where , and . is a boundary point and is an interior equilibrium.

At an equilibrium point , the Jacobian matrix of (18) takes its form aswhere

We know that the local stability of an equilibrium point is determined by the eigenvalues of the above matrix. That is, an equilibrium will be stable if the inequality holds for every eigenvalue of the Jacobian matrix ; it will be unstable if there exists an eigenvalue such that . And in a discrete system, there maybe exist a critical case that holds for each eigenvalue and holds for at least one eigenvalue . In a critical case, the system may be stable or unstable, and it may be needed to check the high order items to discuss the system stability.

Taking the expression of the boundary point into (20), we have that the Jacobian matrix at is simply given by which has four eigenvalues: Since holds evidently, we conclude that the boundary equilibrium is unstable.

Next we study the local stability of the interior equilibrium . Notice that, at the interior point , the item must be zero. Then from (20), we have of which the characteristic polynomial , where

The local stability of is determined by the roots of : if all the roots lie inside the unit disk (i.e., holds for ever root ), then will be asymptotically locally stable. About a polynomial , the following three statements give the equivalent conditions for all the roots lie inside the unit disk (Schur-Cohn Criterion, see, e.g., [35]):(i)(ii)(iii)The determinants of the matrices and the matrices are all positive, where

In our model, and ; then is obviously positive and by calculating we get which tells that holds also. So we conclude that in our model the equilibrium will be asymptotically stable if and . These conditions can be reduced to the following inequality systems:(S1)(S2)

If the inequalities in (S1) and (S2) are all satisfied, then the equilibrium will be locally asymptotically stable in the evolution of (18). If these conditions fail to hold, then the system may lose stability and even exhibit much complicated dynamical behaviors, which will be shown in the following section of numerical simulation.

4. Numerical Simulation

In this section, the constants , , , and are fixed and the different influence of the two delay weights and on the system stability are simulated. The complicated dynamical behaviors are shown when the system loses its stability. Simply, in the next numerical simulations, the constants , , , and are fixed as , , , and .

4.1. The Influence of the Delay Weight

In Figure 1, for different and the same with respect to the adjustment speed , the delay weight (by player 1 for the marginal profit) is fixed as and the bifurcation diagrams are plotted for different (the delay weight by player 2 for the market price): . We find that the system will be stable when the value of the adjustment speed stays at a low level, while when increases, the system will lose stability through bifurcations. Comparing the three diagrams in Figure 1, there is little difference in the values of for the system to lose stability.

To study the influence of the delay weight further, the stability regions of the system are shown in Figure 2 for four different values of . From Figure 2, we see that the stability regions for different have little distinction for different .

4.2. The Influence of the Delay Weight

Figure 3 fixes the weight and plots the bifurcation diagrams for different : From Figure 3, we see that the system will be stable when the value of the adjustment speed stays at a low level, while when increases, the system will lose stability through bifurcations. Comparing the three diagrams in Figure 1 done for different , we find that the point causing the stability loss is much more distinguishable.

The stability regions of the system are shown in Figure 4 for four different values of ). By numerical simulations, Figure 4 shows clearly that the delay weight has great effect on the system stability. It is also observed that the stability region for an intermediate level of is larger than the one for the nondelayed case and the full delayed case . That is, an intermediate marginal profit delay weight will promote the system stability (which is a result similar to the one obtained in the homogeneous-players model [26]).

Comparing Figures 1 and 2, we say that the marginal profit delay weight has more influence on the dynamical behavior of the system but the market price delay weight has less effect. And comparing Figures 3 and 4, it is clear that the delay weight plays little role in the system stability while the delay weight has great effect.

4.3. Complexity of the System

As is shown in Figures 1 and 3, the system may lose stability through a process of bifurcation when the stability conditions are broken. From the two figures, we also find that the system stability loss occurs in two different ways. Indeed, in Figures 1(a), 1(b), 1(c), and 3(c), the system loses stability through a period-doubling bifurcation, whereas in Figures 3(a) and 3(b) the stability loss is evidently due to a Neimark-Sacker bifurcation. For this observation, we can give visually displays by two-dimensional phase portraits. For instance, four phase portraits of the trajectories and associated with Figure 1(b) are plotted in the first line of Figure 5 and the ones associated with Figure 3(a) are plotted in the second line. The first line of phase portraits evidently shows that in the stability loss process there are periods , which is a doubling bifurcation process and at last drives the system to chaos. The second line shows a typical Neimark-Sacker bifurcation process, where closed invariant curves take place when the system loses stability.

5. Chaos Control

We know that chaos is not desirable in a real economic system and it is often hoped to be controlled so that the dynamic system could run in a stable status. In this section, we show that the usually used method of time-delayed feedback control (e.g., [32, 3638]) can also be used to control the chaos in our model.

Using the time-delayed feedback control method, we adjust (18) by inserting a control into the right hands of the first two equations, where is the controlling coefficient. Then we get a controlled system as follows:which can be rewritten as

It is obvious that controlled system (29) and (30) has the same equilibrium points as original system (18). At the interior equilibrium , the Jacobian matrix of system (30) takes its form as

By a similar approach in Section 3 to get the stability conditions (S1)-(S2) for the original system (18), we can also give the stability conditions of controlled system (30) as follows:(S3)(S4),where

To show the availability of the control method, we consider the case . Figure 1(c) shows that for this parameter group original system (18) exhibits chaotic behaviors in its evolution process. For controlled system (30) with this group of parameters unchanged and , Figure 6 plots a bifurcation diagram with respect to the controlling coefficient . It is observed that the controlled system will get out of chaos and be asymptotically locally stable once the controlling coefficient is increased to about 0.34. Theoretically, if we take these parameter values into the stability conditions (S3) and (S4), we can obtain that the condition for the system which is locally asymptotically stable will be .

Below we numerically check the regions determined by stability conditions (S3) and (S4) for more chaos control situations. Figure 1 shows that when , , and , the original system will be chaotic; Figure 2 shows that when , , and , system (18) will also exhibit complicated behaviors. But stability conditions (S3) and (S4) tell that controlled system (30) will be stable, once the controlling coefficients and are set to meet the inequalities in (S3) and (S4). As numerical solutions of the inequality system (S3) and (S4), Figures 7 and 8 plot the regions for all the cases mentioned. Figures 7 and 8 mean that in each case the system chaos will be controlled if the controlling coefficients and together lie within the plotted regions. Figure 7 is done for different levels of and Figure 8 is done for different levels of . By comparing the three cases in Figure 7, we observe that there is little difference among the stability regions. That is, there is little change in the system controllability when is set at varied levels. However, from Figure 8 we see that has great influence on the system controllability: the stability region for a more intermediate is much larger so that the system can be more easily controlled. By these observations we draw one more conclusion that in our model the two different delays play also different roles on the system controllability: the delay of the marginal profit has more influence on the system control than the delay of the market price, and an intermediate level of delay for the marginal profit has a positive effect on the chaos control.

6. Conclusion

In this paper, we reconsider the duopoly game with time delay and formulate a model with two different delay structures. The main consideration is that the price history may be observed by all players in the market and there may be players making decision according to the historic information of the market price. In this work, we modify the model in [26] (where the two players are assumed to be homogeneous so that both players make decision according to the profit average with time delay). We assume that there are two different kinds of producers: one sets time delay on its marginal profit so that it adjusts its output by the approach discussed in [26]; but another sets time delay on the price variable and makes decision in accordance with the averaged data from the price history.

The dynamics of the players’ output adjustments has been established, the equilibrium solutions are obtained, and the stability of the equilibriums is mathematically analyzed. The boundary equilibrium is proved to be unstable and the stability conditions for the interior equilibrium are obtained. Numerical simulations are done to show the influence of the main model parameters on the system stability and on the complicated behaviors of the system when losing stability. It is observed that there is difference in the roles played by the two delay structures. The time delay for the marginal profit has more influence on the dynamical behaviors of the system but the market price delay has less effect. An intermediate level of the delay weight for the marginal profit can promote the system stability and thus reduce the occurrence of complex behavior such as bifurcation and chaos. It is shown that the system stability losing may be due to either period-doubling bifurcations or Neimark-Sacker bifurcations. Numerical simulations show that the chaotic behaviors can be stabilized by the time-delayed feedback control, and the two different delays play also different roles in the system controllability: the delay for the marginal profit has more influence on the system control than the delay for the market price, and an intermediate level of delay for the marginal profit has a positive effect on the chaos control.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

Financial support by the National Natural Science Foundation of China (nos. 71171098, 51306072, 71690242, and 71273120) is gratefully acknowledged.