Research on a Cournot–Bertrand Game Model with Relative Profit Maximization

1is paper considers a Cournot–Bertrand game model based on the relative profit maximization with bounded rational players. 1e existence and stability of the Nash equilibrium of the dynamic model are investigated.1e influence of product differentiation degree and the adjustment speed on the stability of the dynamic system is discussed. Furthermore, some complex properties and global stability of the dynamic system are explored. 1e results find that the higher degree of product differentiation enlarges the stable range of the dynamic system, while the higher unit product cost decreases the stable range of price adjustment and increases the one of output adjustment; period cycles and aperiodic oscillation (quasi-period and chaos) occur via period-doubling or Neimark–Sacker bifurcation, and the attraction domain shrinks with the increase of adjustment speed values. By selecting appropriate control parameters, the chaotic system can return to the stable state. 1e research of this paper is of great significance to the decision-makers’ price decision and quantity decision.


Introduction
Oligarchic market is one type of market structures controlled by several players who produce homogeneous products and take actions on the basis of their rival reactions for seeking profit maximization. Oligopoly game is extensively studied by many scholars; especially, in recent years, the dynamics of the oligopoly game has been extensively studied as a branch of mathematical economics.
A large significant number of previous articles have analyzed the stability of discrete dynamic oligarchic models as well as the dynamic behaviors with players having the naive rule, adaptive rule, and bounded rational expectation rule, respectively. e dynamic Cournot game was studied by many researchers, such as Kopel [1], Agiza [2], Bischi et al. [3], Bischi and Lamantia [4], Agiza and Elsadany [5], Fanti and Gori [6], Fanti et al. [7], Zhu et al. [8], Gori et al. [9], and Elsadany [10]. Most of the aforementioned works analyzed the adjustment process of quantity of the players with bounded rationality and found the periodic cycles and chaotic behaviors occurred in the Cournot model when the adjustment speed is faster than a certain value.
Many scholars researched on the Bertrand game and explored its complex features using nonlinear dynamic theory. Ma and Guo [11] analyzed the impacts of information on the dynamical price game, including the changes of the stability region and attraction basin of the Nash equilibrium point. Yi and Zeng [12] developed a duopoly Bertrand model with quadratic function and analyzed the stability region, chaotic behaviors, and chaos attractors of the Bertrand model. Elsadany and Awad [13] investigated a duopolistic Bertrand competition market with environmental taxes when the public firm is privatized and not privatized. e dynamical behaviors of the models are studied by numerical simulation and period-doubling bifurcation and chaos appeared. In addition, the research about the dynamic Stackelberg game was also studied by many scholars, such as Shi et al. [14], Peng and Lu [15], Peng et al. [16], and Askar [17].
In recent years, the Cournot-Bertrand game received more attention by some scholars. Tremblay and Tremblay [18] established a Cournot-Bertrand duopoly model with differentiated products and discussed the Nash equilibrium of the model. Askar [19] presented a Cournot-Bertrand model based on a nonlinear price function and analyzed the stability of the Nash equilibrium point. Naimzada and Tramontana [20] considered a dynamic Cournot-Bertrand duopoly model with product differentiation and studied the role of the best dynamic response for the stability of the equilibrium point. Wang and Ma [21] constructed a mixed duopoly game model with limited information, and the existence and stability of the Nash equilibrium point are investigated. Ma and Pu [22] researched the chaotic behaviors the Cournot-Bertrand duopoly model using nonlinear dynamics theory. Different from the above works, the Cournot-Bertrand game model developed in this paper is based on relative profit maximization of policymakers. e relative performance or relative profits as a competition mechanism have paid attention by some literature studies. ere are cases in point, on the one hand, even if a person earns big money, he is not happy enough and may be disappointed if his brother/sister or close friend earns bigger money. On the other hand, even if he is very poor and his neighbor is poorer, he may be consoled by that fact [23]. e firm's relative profit is the difference between its absolute profit and that of its competitors. Pursuing the relative profit or utility is considered to be rational based on human nature. e firms consider their relative profits rather than absolute profits as the business objective is closer to reality.
Miller and Pazgal [24] argued that business people are likely to be inherently competitive, caring more about the relative position and status than about absolute profits, and researched the equilibria of profit maximizing owners selecting different types of managers with various types of competition (a variety of attitudes toward relative performance). e equilibrium with the quantity setting behaviors and price setting behaviors of the duopoly game under relative profit maximization in the static state was explored in depth by Satoh and Tanaka [23], Hattori and Tanaka [25], Tanaka [26,27], and so on. Some scholars have studied the game model with relative profit maximization under dynamic decision-making. Fanti et al. [7] developed the dynamics of the nonlinear Cournot duopoly game with managerial delegation and bounded rational players. Furthermore, the on-off intermittency and blow-out bifurcations as well as coexistence of attractors are observed. Lu [28] studied the influence of the relative profit maximization on endogenous timing in the Stackelberg game with public and private oligopoly. Elsadany [10] developed a dynamic Cournot duopoly model with relative profit maximization and externality cost functions. Askar [17] investigated the dynamic characteristic of a duopoly game of two firms based on their relative profit maximization.
To the best of our knowledge, little literature has studied the Cournot-Bertrand game with relative profit maximization using nonlinear dynamic theory. It is a very interesting item to study the dynamic characteristic of the Cournot-Bertrand game model based on the relative profit maximization.
Our theoretical contribution is as follows. e first contribution is to construct a dynamic Cournot-Bertrand game model, in which one firm adopts price as his decision variable and the other adopts quantity as his decision variable based on the relative profit maximization. e second contribution is to study the influence of parameters changing on the stability and profitability of the dynamic Cournot-Bertrand game model and to obtain some interesting phenomenon, such as existing period-doubling bifurcation and Neimark-Sacker bifurcation.
is paper is organized as follows. In Section 2, the dynamic Cournot-Bertrand model with the relative profit maximization is given. In Section 3, the stability of Nash equilibrium points of the model is analyzed. Numerical simulations are carried out to present the complex features of the dynamic system in Section 4. Section 5 gives the global stability analysis of the dynamic system using basins of attraction. In Section 6, the variable feedback control method is adopted to remove the chaos of the model. Finally, some conclusions are concluded in Section 7.

Model Construction
is paper considers two firms (firm 1 and firm 2) in the same market with differentiation products (x 1 and x 2 ) with linear inverse demand function and linear cost function. Firm 1 regards the quantity of the products as the strategy variable, while firm 2 considers the price of the products as the strategy variable.
at is to say, firm 1 competes in product quantity (q 1 ) and firm 2 competes in product price (p 2 ). Both firms will take actions based on the information from their rival and market for seeking relative profit maximization. is paper mainly focuses on contributing to construct a dynamic Cournot-Bertrand game model based on relative profit maximization and investigate the dynamic behaviors of players and the complex properties of the system. e inverse demand functions of products are as follows (for more details, refer to [29,30]): where p 1 , p 2 denote the retail prices of products from firm 1 and firm 2, respectively, and q 1 and q 2 are the product quantities of firm 1 and firm 2. Moreover, d ∈ [0, 1] denotes the degree of product differentiation between two firms. When d approaches one, the products of two firms become less differentiated; when d � 1, the two products are completely homogeneous; and when d � 0, the structure of the market transfers to the monopolist from the duopoly. e cost functions are in linear form: where c > 0 is the unit product cost.
Considering that the strategic variables are different between two firms, the inverse demand function of firm 1 and the demand function of firm 2 are given, respectively: Complexity en, the absolute profit functions of two firms are as follows: Based on previous assumptions, the relative profits of the players are defined as the difference between their absolute profits and those of other players. e relative profits of firm 1 and firm 2 are represented by Φ 1 , Φ 2 , respectively, given by e first partial derivative of relative profit maximization of firm 1 with respect to q 1 is Similarly, the first condition of the relative profit maximization of firm 2 with respect to p 2 is given by e second condition of the relative profit maximization of firm 1 and firm 2 is easily calculated as follows: en, the relative profit functions of two firms are concave, and unique maximum value exists.
According to the actual market, the players in the market cannot get all the information they need in collecting and processing. Hence, both firms express bounded rational behaviors when making strategy decisions of their own. e optimal actions of two firms based on their rival's reaction are derived, respectively: For achieving the relative profit maximization, the firms will choose the optimal actions based on the marginal relative profit. If the marginal relative profit is positive (negative) in the current time period, the firms will increase (decrease) their quantity (price) in the next time period. erefore, the discrete dynamical adjustment mechanism of the Cournot-Bertrand game model is described by where α, β > 0 represents the adjustment speed of firm 1 and firm 2, respectively. Substituting equations (8) and (9) into equation (11), dynamical system (11) can be obtained as Consequently, the time evolution of the Cournot-Bertrand game with bounded rational expectation is built. e next section analyzes the stability and complex features of dynamic system (12).

Stability Analysis of the Dynamic System (12)
In this section, we emphasis on analyzing the stability of the Nash equilibrium points of dynamic system (12). Letting q 1 (t + 1) � q 1 (t) and p 2 (t + 1) � p 1 (t), four equilibrium solutions are computed as follows: In order to ensure that the equilibrium solutions are nonnegative, conditions should be satisfied as follows: Lemma 1. E 1 , E 2 , and E 3 are the unstable equilibrium points, and E 4 is a unique Nash equilibrium solution of dynamic system (12).
At E 2 , the Jacobian matrix is given as follows: Two eigenvalues of matrix J(E 2 ) are obtained: erefore, E 2 is a saddle point. By the same way, we can prove that E 3 is a saddle point, and E 4 is the only positive Nash equilibrium solution. is completes the proof. en, the Jacobian matrix at the positive Nash equilibrium point E 4 is where According to Jury's conditions, the dynamic system (12) is asymptotically stable if and only if it satisfies the following conditions: 1 − N > 0, More precisely, for inequality (19), when |N| > 1, the Hopf (in the continuous system) or Neimark-Sacker bifurcation occurs; when P(1) < 0, it means one of the real eigenvalues is bigger than 1, and then the tangent, pitch, or transcritical bifurcation will take place; and when P(−1) < 0, it denotes one of the real eigenvalues is smaller than −1, and then the flip (period-doubling) bifurcation will be raised. In addition, if the aforementioned three conditions are not fulfilled at the same time, then the dynamic system (12) will lose its stability and become unstable, so the conditions above defined the surfaces in the parameter space on which kinds of bifurcations will appear.
However, solving the inequality equation (19) is very complicated. Next, we give the stable region of the dynamic system (12) through the parameter basins. According to the current situation, we take the parameter values as follows: e parameter basin map is also called 2D bifurcation map, which can clearly show the path of the system from the stable state to chaotic state. e parameter basins are drawn in Figure 1 with d � 0.2, in which different colors represent different period cycles, for example, red (stable), blue (2-period), green (3-period), yellow (4-period), orange (5-period), light blue (7-period), brown (8-period), magenta (chaos), and white (divergence). From Figure 1, when the values of the adjustment speed of two firms are in the stable region (red region), the dynamic system (12) finally stabilizes at the Nash equilibrium point after many iterations. When the values of the adjustment speeds of two firms escape from the stable region (red region), the Nash equilibrium point will become unstable, perioddoubling cycles or Neimark-Sacker bifurcations occur and even evolve into chaos finally. Furthermore, it is important that the player will drop out the market when the adjustment speed is greater than a certain value (white region). Figure 2 gives the parameter basins of dynamic system (12) when d and c take different values. Comparing the sizes of stability regions (red region) in Figures 1 and 2, we find that higher product differentiation enlarges the stable region of the dynamic system (12), and the higher unit product cost decreases the stable region of price adjustment and increases the stable region of output adjustment. So, two firms should adjust the parameters according to the actual market conditions so that the dynamic system (12) is in a stable state.

Numerical Simulations
In this section, the numerical simulation is carried out in depth to analyze the influence of the price speed adjustment and the degree of differentiated products on the stability of the Nash equilibrium point of the dynamic system (12). Moreover, the dynamic complex behaviors of the dynamic system (12) are analyzed by simulations, such as bifurcations, chaos, the largest Lyapunov exponent (LLE), strange attractors, and sensitive dependence on initial conditions.

e Effect of the Adjustment Speed on the Stability of the Dynamic System.
e dynamic characteristics of the system can be well shown by the one-dimensional bifurcation diagram and the LLE with one-parameter changing when keeping other parameters fixed. Figure 3 4 Complexity shows the bifurcation diagram and the corresponding LLE with α varying from 0 to 3.2 when β � 1. Along with the increase of α from 0 to 2.355, the dynamic system (12) is in the stable state; when α > 2.355, the dynamic system (12) losses its stability and periodic-doubling occurs, including period-2 and period-4. Finally, system (5)     system (12) becomes chaotic finally through period-doubling bifurcations, and the positive LLEs indicate the appearing of chaos. As can be seen from Figure 5(a), the bifurcation diagram with respect to α is quite different from Figure 3 when β � 2.
e dynamic system (12) remains in the stable state when α < 2.290, with α increasing, and when exceeds 2.290, the dynamic system (12) becomes unstable, Neimark-Sacker bifurcation appears and then enters into the chaotic state at last. As shown in Figure 5(b), the dynamic behaviors of dynamic system (12) become more complicated when β is fixed at 2.3, and the Nash equilibrium point is in a 2-period state as α varies from 0 to 1.342; when 1.342 < α < 1.981, the dynamic system (12) returns to the stable state from the period-2 state, and then the system undergoes Neimark-Sacker bifurcation and enters into chaos when α > 1.981. In addition, as one can find that there is an intermittent odd cycle (period-7) in the chaos in Figures 5(a) and 5(b), it implies that the so-called topological chaos is created.
From the perspective of economics, when the adjustment speed of firms exceeds a certain value, the Nash equilibrium of the system will lose its stability, and the system becomes unpredictable and disorder; it is hard for firms to deal with   Complexity the complex scenario like that. Hence, firms should keep the adjustment speed in a proper range in order to avoid reducing efficiency of the market and profits of themselves.
Some strange attractors of the system are plotted in Figure 6 with respect to (q 1 , p 2 ) in the two-dimension map under different parameter combination values of the adjustment speed of firms. e limit cycle and coexistence of several attractors are displayed in Figures 6(a)-6(c) when α � 2.5 and β has different values. e strange attractors with the fractal structure are exhibited in Figure 6(d) which implies existence of chaotic behaviors.
An important feature of the chaotic system is the sensitivity to initial values. Here, keeping p 2 unchanged and α � 2.5, β � 2.615, Figure 7 shows the differences of the evolution process of the dynamic system (12) when the initial value of q 1 only changes, 0.001. We can see that the values of q 1 and p 2 have no difference in the first 71 time iterations, but after that, the values of q 1 and p 2 show great difference. at is to say, the small difference of the initial value will cause great deviation after much iteration, which provides us with the enlightenment that decision-makers should be more cautious in choosing the initial values of decision variables.

e Effect of Degree of Product Differentiation and Unit Product Cost on System Stability.
e bifurcation diagram and the corresponding LLE with respect to d are drawn in Figure 8 when α and β are fixed at 1.5; the dynamic system (12) is in the stable state at first, and as d increases, we can see that the product quantity of firm 1 is unchanged until d � 0.935, while the price of firm 2 is sharply decreased at the same time; then, a Neimark-Sacker bifurcation occurs when d > 0.935, and the dynamic behaviors of players fall into the  Complexity quasi-periodic state; it is hard for the bounded rational players to predict the trajectory of its quantity or price and make a long competition strategy. In addition, the bifurcation diagram with respect to d is shown in Figure 9 when α � 1.5, β � 2.4. e dynamic system (12) enters the stable state from the chaotic state, and then enters the chaotic state when d ∈ (0.25, 0.85), so it is an efficient way to remove chaos by controlling d in a proper range of value. As d > 0.877, the dynamic system (12) exhibits Neimark-Sacker bifurcation and stays in the quasi-period state. Compared with Figure 8, the dynamic behaviors of the market become much more complicated in the chaotic state. e LLE plots in Figures 8 and 9 are used to verify existence of stable cycles, periodic cycles, and chaotic behaviors. One can obtain a better understanding of the effects of degree of product differentiation on stability of the Nash equilibrium point and complex characteristics of the dynamic system. e bifurcation diagrams with respect to c are drawn in Figure 10 when the dynamic system (12) is in the stable state and unstable state. From Figure 10(a), when 0 < c ≤ 0.88, the output of firm 1 decreases first and then goes into the unstable state; when c > 0.88, the output of firm 1 equals to zero; when 0 < c ≤ 1, the price of firm 2 increases first and then goes into the chaotic state through flip bifurcation. From Figure 10(b), we can find that when the dynamic system (12) is in the unstable state, the dynamic system (12) will go into chaos early with c increasing.
From the above analysis, we can get the following conclusions: (1) when the dynamic system (12) is in the stable state, a higher substitution between products may destabilize the stability of the Nash equilibrium point and lead to appearance of aperiodic cycles. In addition, the product quantity of firm 1 remains in a certain value and unchanged, while the price of firm 2 is decreased with the increase in degree of product differentiation. It denotes that the player who takes price as the strategy variable is more sensitive to changes of degree of product differentiation, which is different from the classic Cournot-Bertrand game discussed in previous articles. When the dynamic system (12) is in the chaotic state and the degree of product differentiation tending to 0 or 1 , periodic cycle and chaos will appear in system (12). (2) When the dynamic system (12) is in the stable state, as the unit product cost exceeds a certain value, firm 1 will withdraw from the market; when the dynamic system (12) is in the chaotic state, the increase of unit product cost of two firms will make the dynamic system (12) go into chaos early. So, two firms should control the degree of product differentiation and unit product cost so as to avoiding the market to go into a chaotic state.

Global Stability of the Dynamic System (12)
An effective method for analyzing the influences of parameter changing on the global stability is the basins of attraction in which there include attraction domain and escaping area. It is the set of initial conditions; if the initial variable is taken from the attraction domain, the system will emerge the same attractor after a series of iterations. If the initial variable is outside the basins of attraction, the system will fall into divergence at last. By fixing the parameter values of the dynamic system (12)as mentioned above, the basins of attraction about q 1 and p 2 of the dynamic system (12) are shown in Figure 11 under different values of adjustment speed, in which the red region denotes the stable attraction domain and the white region denotes the escape area. We found that the attraction domain shrank with the increase of adjustment speed values.
Parameter's value has a great influence on the global stability of the dynamic system. Figure 12 shows the basins of attraction of q 1 and p 2 with c and d having different values when α � β � 1.5. Comparing Figures 11 and 12, we find that the basins of attraction decrease with the increase of unit product cost and degree of product differentiation.
From an economic point, the initial values of variables of two firms should be in basins of attraction in order to maintain market stability.

Chaos Control
Generally, the chaos of the economic system is regarded as an obstacle for players to make a long business strategy and achieve their own business objective. In recent years, the methods in controlling bifurcation and chaos are proposed by many scholars such as the OGY method, the variable feedback control method, and the delay feedback control (DFC) method. Consequently, this paper uses the variable feedback control, which is widely applied in supply chain management and insurance market to delay and eliminate the chaos of dynamic system (12).
We assume that the discrete dynamic model of system (6) takes the following form: en, the controlled system can be described as where k is the control parameter, which can be regarded as government regulation or the learning ability of the players. Figure 13 shows the bifurcation diagram and the LLE of the controlled system (21) when α � 2.5 and β � 2.615. One can see that chaos, period-doubling bifurcation, and 7- period cycle occur twice; when 0 < k < 0.129, then controlled system (21) undergoes chaos range, quasi-period range, and returns to a stable state with k increasing. e corresponding LLE becomes negative when k > 0.378, which verifies that the controlled system (21) returns to the stable state. In Figure 14, the trajectory of q 1 and p 2 is plotted when α � 2.5, β � 2.615 and k � 0.38; after many iterations, the control system (21) gradually returns to the stable state from the chaotic state.

Conclusion
is paper devotes to establish a dynamic Cournot-Bertrand game model based on relative profit maximization. e stability of the dynamic system is analyzed, and the influences of variable adjustment speeds, unit product cost, and product differentiated degree on complex characteristics of the dynamic system are investigated. e results showed that the higher degree of product differentiation enlarges the stable range of the dynamic system, while higher unit product cost decreases the stable range of price adjustment and increases the one of output adjustment. Some complex properties and global stability of the dynamic system are explored, and it is found that period cycles and aperiodic oscillation (quasi-period and chaos) occur via period-doubling or Neimark-Sacker bifurcation. In addition, the limit cycle and coexistence of several attractors will appear in some cases, and the attraction domain shrank with the increase of adjustment speed values.
e study of complexity of the Cournot-Bertrand game model based on relative profit objectives has theoretical and practical significance. Nonetheless, there are some factors that are not taken into account, for example, if the risk attitude of decision-makers is taken into account, the Cournot-Bertrand game model will be more realistic. Firms can consider dual-channel to expand market share. We believe that the viewpoint of this paper will lay a dynamic foundation for future research in these directions.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.

Authors' Contributions
Huang Yi-min and Guo Yan-yan revised the paper, Li Qiuxiang provided research methods, and Zhang Yu-hao wrote the original draft. Complexity 11