1. Introduction
Oligopoly is a competition game where the market situation is controlled by a few competing firms producing similar or homogeneous products. The concept of static competition of duopoly game was first presented by Cournot [1]. The duopoly model is a competitive game where two companies control and have a dominant influence on the market. A typical feature that many articles share is the assumption of oligopoly model activities in order to maximize profit function as a single goal [2]. The game strategy of each company is related to each other and depends not only on market demand, but also on the planning expectations of the competitor. The pioneering contribution of Puu’s work [3] has shown the potential for the complex dynamics of the duopoly competition of Cournot game. In literature [4,5,6,7,8,9,10], the dynamic behaviors and chaos of Cournot models have been discussed on the basis of naive expectations. Subsequently, an analysis of the global dynamics of the Cournot Competition game with bounded rationality was presented by Bischi and Naimzada [11]. All of these models have shown that equilibrium points could lose their stability, leading to periodic or chaotic behavior. Recently, the stability and complexity of Cournot competition games with limited rationality attracted more attention of researchers.
In this paper, we analyze the Bischi–Naimzada duopoly game based on firms considering a profit maximizer under a sales constraint. Taking into account that the sales constraint is a more realistic economic situation than maximizing the only firm’s profit. Such an added objective motivates our contribution in the current paper and can improve and extend the region of stability for the Nash equilibrium point. Furthermore it generalizes other existing models in literature. We study the stability of the dynamic equilibrium point leading to an in-depth understanding of the dynamics of the game and the corresponding economic interpretations. We use stability criteria as well as the numerical simulation to determine the local and global stability of the equilibrium points as well as to discuss the internal complexity of the game. Indeed studying the model according to the parameters of the speed of adjustment mechanism and the parameters of the sales constraints that we propose is capable of leading the model to different cycles or forms of complex dynamics involving potential economic instability. On the other hand, small perturbations of related standards may lead to sudden jumps from one attraction basin to another, based on the complexity of the basin structures and the phenomena of coexistence even in the absence of chaos. We obtain the critical curves of the non-invertible map, and it turns out that a sequential change in game’s parameters leads to the collapse of game. Moreover, the behaviors of destabilization also involve a subcritical period doubling bifurcation, some periodic, aperiodic and chaos has been merged with an increase in the speed of adjustment parameters. In particular, we show that if the sales constraint parameters are taken to be high, this makes Nash equilibrium point goes back to local stability and has a reflective effect of increasing the model’s adjustment speeds.
This paper is divided as follows. In Section 2, we report related works in this direction of research. In Section 3, the Bischi–Naimzada duopoly game based on maximizing profits and sales constraints is introduced. In Section 4, the local stability conditions for the equilibrium points are determined. In Section 5, numerical simulation, local bifurcations, basin of attraction, critical curves and path to complex dynamics are investigated. Finally, the conclusion and discussion are set out in Section 6.
2. Literature Review
Research on the complex dynamics of competition games may be classified according to the type of model or the number of players or the type of players expected to be homogeneous and heterogeneous. Agiza et al. [12,13] analyzed the dynamics of two oligopolistic games based on a limited rationality with a linear and nonlinear demand function. Bischi et al. [14,15] investigated chaos and bifurcations in two dynamic oligopoly games. The oligopoly was presented in the economic models of heterogeneous expectations explored in [13,14]. A multi-market game based on Puu’s approach [16] has been introduced in [17]. Fanti et al. [18] examined the local and global dynamics of a homogeneous Cournot duopoly game with isoelastic demand function. Askar and Al-Khedhairi [19] analyzed the nonlinear duopoly games with product differentiation. Ma et al. [20] analyzed the complex behaviors and applied useful control of chaos to master-slave duopoly game. Agliari et al. [21] studied the dynamics of a differentiated Cournot duopoly game. Askar [22] explored the dynamics of Cournot duopoly game based on uncertainty cost. Andaluz et al. [23] studied the dynamic of triopoly Cournot and Bertrand games with differentiated products and heterogeneous expectations. Baiardi and Naimzada [24] investigated the dynamics of evolutionary Cournot oligopoly game with imitators. The dynamic behaviors of the two-stage Cournot duopoly model with an RD spillover effect have been studied in [25]. Al-Khedhairi [26] introduced the Cournot duopoly game with a generalized bounded rationality. The influences of asymmetric market information on the dynamics of duopoly game have been explored in [27]. Other interesting papers have been published in the literature, such as the works by Naimzada and Sbragia [28] which has studied oligopoly games with nonlinear demand and cost functions. Tramontana [29] has examined heterogeneous duopoly based on isoelastic demand. Askar [30] has analyzed Cournot–Bertrand game with differentiated goods. The dynamic of oligopoly games with goods differentiation discussed in [31]. The dynamics of differentiated Bertrand duopoly game investigated in [32]. The complex dynamics of delayed Bertrand duopoly game explored in [33]. A remanufacturing duopoly model with nonlinear cost studied in [34]. Effect of information asymmetry in the bounded Cournot duopoly game discussed by Ueda [35].
All competed firms in the works mentioned above were trying to maximize their profit functions without constraints as unique objectives. Baumol [36] presented an alternative model of increased sales as the company increases its sales on the basis of a minimum profit. Fischer [37,38] concluded that the idea of maximizing profits subject to minimum sales constraints was much more rational for business. More recently, Ahmed et al. [39,40,41] introduced dynamic multi-objective games as an alternative to firms which think with a single objective optimization. Mert [42] studied the duopoly model with two objectives to increase sales and maximize profits, giving Nash equilibrium a more stable when it comes to maximizing sales and profits rather than maximizing profit. Ibrahim [43] achieved the stability conditions for the Nash equilibrium point in the duopoly model with minimum sales constraints. Tian et al. [44] investigated coordination and control of multi-channel supply chain based on consumers’ channel preferences and sales effort. Pansera et al. [45] studied bifurcation analysis of a delayed duopoly game with R&D spillover and price competition. More recently, the book includes new directions for nonlinear dynamics in games, see [46]. In [47], the price and quantity that are traded off between buyer and seller have been determined in bilateral monopoly. Investigation of excessive competition and a lack of trust and their influences that may affect rational choice in coalitional game have been analyzed in [48].
3. Bischi–Naimzada Duopoly Game under Minimum Sales Constraint
Assuming a market dominated by two firms producing homogeneous or same products. Letqibe the output of firmi=1,2andQ=q1+q2 denotes the total output of the two firms. Suppose that the market demand function is linear in form [11]:
p=a−b(q1+q2)
where p indicates the price of the products. We assume that both firms adopt the following linear cost function,
Ci(qi)=ci qi
Since the total profit for the firm is the total revenue minus the total cost, the profit for each firm is given by
πiq1,q2=qi(a−bQ)−ci qi
Assuming that firms have a constant constraint on minimum sales constraint [43]. The firm i believes that it maximizes its profit under the condition on sales constraint, then we have the following optimization problem
maximizesqi(a−bQ)−ci qisubjecttoqi(a−bQ)≥Si,
It is equivalent to maximize the payoff function as the objective function of firm i, which has the form
Li=qi(a−bQ)−ci qi−μiqi(a−bQ)−Si,
where the parameterμiis positive and is associated with the sales constraint. By differentiatingLiwith respect toqi, we obtain
∂Li∂qi=(1−μi)(a−2bqi−bq−i)−ci=0,i=1,2.
The objective of each firm is to maximize its profit on the basis of minimum sales constraints. The objective is therefore to maximize the function defined in (5) which combines the function of profit with the constraints of sales. The Dynamic Game
We consider the situation where each firm adopts the same mechanism in adjusting its output in next period of time. We are of the opinion that these two firms are rationally limited, that is, they have limited information, each from the other, so that they will decide their outputs depending on the local estimation behavior of the objective function during the period [11,12,13]. Thus, the dynamic system that describes this game is in the form:
qi(t+1)=qi(t)+νi qi(t)∂Li∂qi,i=1,2.
whereνi>0represents the adjustment speed of firm i.
Substituting (6) in (7), the dynamics of the two firms’ outputs through a discreet dynamical system will be as follows:
q1(t+1)=q1(t)+ν1 q1[(1−μ1)(a−2bq1−bq2)−c1]q2(t+1)=q2(t)+ν2 q2[(1−μ2)(a−2bq2−bq1)−c2]
The game mentioned above and described by the map (8) is a generalization of the Bischi–Naimzada game [11] when there is a constraint on sales. The system described in (8) determines the trajectories of competing firms, and it is a noninvertible two-dimensional map as shown later. We are investigating the dynamics of the model by discussing the effect of the parameters on this system. In order to address the dynamical behavior (8), we analyze the stability regions and the bifurcation of the equilibrium points by which the points may be destabilized.
4. Equilibrium Points and Their Local Stability
To obtain the fixed points for the model (8), the fixed point conditionsq1(t+1)=q1(t)=q1∗andq2(t+1)=q2(t)=q2∗, are used,
q1[(1−μ1)(a−2bq1−bq2)−c1]=0,q2[(1−μ2)(a−2bq2−bq1)−c2]=0
Solving the algebraic system (9), we get the following four fixed points of the system (8):
E0=(0,0),E1=(a(1−μ1)−c12b(1−μ1),0),E2=(0,a(1−μ2)−c22b(1−μ2)),E∗=q1∗,q2∗=a(1−μ1)(1−μ2)+c2(1−μ1)−2c1(1−μ2)3b(1−μ1)(1−μ2),a(1−μ1)(1−μ2)+c1(1−μ2)−2c2(1−μ1)3b(1−μ1)(1−μ2)
E0,E1,E2are called boundary equilibria andE∗is the Nash equilibrium. They have all non-negative real values which hold under the following inequalities:
2c1(1−μ2)−c2(1−μ1)<a(1−μ1)(1−μ2),2c2(1−μ1)−c1(1−μ2)<a(1−μ1)(1−μ2).
It should be noted that the Nash equilibrium pointE∗ of the game (8) is the same of the Bischi–Naimzada game [30] whenμ1=μ2=0.
To study the local stability of those points we should calculate the Jacobian matrix of system (8) that is given by,
J=1+ν1[(1−μ1)(a−4bq1−bq2)−c1]−ν1(1−μ1)bq1−ν2(1−μ2)bq21+ν2[(1−μ2)(a−4bq2−bq1)−c2],
whose characteristic polynomial is as follows,
f(λ)=λ2−tr(J(E))λ+det(J(E)),
wheretr(J)anddet(J) are the trace and determinant of the matrix (12), respectively. The conditions for achieving the local asymptotic stability of the Nash equilibrium point are that the eigenvalues of the corresponding Jacobian matrix are within the unit circle. According to the Jury criterion [49,50], the conditions of local asymptotic stability of the Nash equilibrium point can be described in detail as
f(−1)=1+tr(J(E))+det(J(E))>0,f(1)=1−tr(J(E))+det(J(E))>0,det(J)<1.
Breaking any of the inequities mentioned in (14) while realizing the other two at the same time, leads to:
-
Flip bifurcation atf(−1)=0,
-
Transcritical or fold bifurcation atf(1)=0,
-
Neimark–Sacker bifurcation at det(J)=1.
It is clear that the system (8) is a nonlinear system and therefore it is difficult to obtain an explicit time-dependent solution. From this point of view, we will study the specific behaviors of this model by applying the stability theory around the equilibrium points.
Theorem 1.
The equilibrium pointE0of the two firms outside the market is
-
a repelling node ifa>max{c1−1ν1(1−μ1),c2−1ν2(1−μ2)},
-
an attracting node ifa<min{c1−1ν1(1−μ1),c2−1ν2(1−μ2)},
-
a saddle node ifa>c1−1ν1(1−μ1)anda<c2−1ν2(1−μ2)ora<c1−1ν1(1−μ1)anda>c2−1ν2(1−μ2),
-
a non-hyperbolic node ifa=c1−1ν1(1−μ1)ora=c2−1ν2(1−μ2).
Proof.
At the equilibrium pointE0the Jacobian matrix (12) is given by
J(E0)=1+ν1a(1−μ1)−c1001+ν2a(1−μ2)−c2.
AsJ(E0)is a diagonal matrix. Then the eigenvalues ofJ(E0)areλ1=1+ν1a(1−μ1)−c1andλ2=1+ν2a(1−μ2)−c2. Ifc1<1+ν1a(1−μ1)andc2<ν2a(1−μ2), thenE0is a repelling node. From the conditions shown on the eigenvalues, the rest of the results can be demonstrated based on the following theorem. □
Theorem 2.
When the Nash equilibrium point has positive coordinates, then the equilibrium pointsE1andE2are saddle points.
Proof.
SubstitutingE1in (12), we get
J(E1)=1−ν1[a(1−μ1)−c1]−ν1(a(1−μ1)−c1)201+ν22[a(1−μ1)(1−μ2)+c1(1−μ2)−2c2(1−μ1)].
The two eigenvalues of Jacobian matrix above areλ1=1−ν1[a(1−μ1)−c1]andλ2=1+ν22[a(1−μ1)(1−μ2)+c1(1−μ2)−2c2(1−μ1)]. If the Nash equilibrium point has positive coordinates, henceλ1<1andλ2>1.Therefore, the equilibrium pointE1is a saddle point. In a certain manner, we can demonstrate that theE2is a saddle point too. This completes the proof. □
In order to study the local stability of the Nash equilibriumE∗, we calculate the Jacobian matrix at this point in the form:
J(E∗)=1−2ν1b(1−μ1)q1∗−ν1b(1−μ1)q1∗−ν2b(1−μ2)q2∗1−2ν2b(1−μ2)q2∗.
Then the eigenvaluesλ1,λ2of this matrix are the roots of the following quadratic equation,
λ2−tr(J(E∗))λ+Det(J(E∗))=0,
wheretr(J(E∗))=2−2ν1b(1−μ1)q1∗−2ν2b(1−μ2)q2∗andDet(J(E∗))=(1−2ν1b(1−μ1)q1∗)(1−2ν2b(1−μ2)q2∗)−ν1 ν2 b2(1−μ1)(1−μ2)q1∗ q2∗.We findtr2(J(E∗))−4Det(J(E∗))=4b2[(ν1(1−μ1)q1∗−ν2(1−μ2)q2∗)2+ν1 ν2(1−μ1)(1−μ2)q1∗ q2∗]. Since all parameters are definitely positive then the following discriminant is positive.
tr2(J(E∗))−4Det(J(E∗))>0.
From above, it is easy to assume that the eigenvaluesλ1,λ2of this matrixJ(E∗)are real eigenvalues. By the way of calculation, we get that:
1−tr(J(E∗))+Det(J(E∗))=3ν1 ν2 b2(1−μ1)(1−μ2)q1∗ q2∗>0,
When the Nash equilibrium point has positive coordinates, it’s easy to see that the eigenvaluesλ1,λ2ofJ(E∗)are less than one. The local asymptotic stability ofE∗is the following necessary condition:
1+tr(J(E∗))+Det(J(E∗))=4−4ν1b(1−μ1)q1∗−4ν2b(1−μ2)q2∗+3ν1 ν2 b2(1−μ1)(1−μ2)q1∗ q2∗<0.
The condition (15) defines the stability region for the Nash equilibrium point of the game (10). We can conclude the above result in the following theorem:
Theorem 3.
The Nash equilibrium pointE∗of the game (8) is locally stable provided that:
4−4ν1b(1−μ1)q1∗−4ν2b(1−μ2)q2∗+3ν1 ν2 b2(1−μ1)(1−μ2)q1∗ q2∗<0.
It is clear that if the condition1+tr(J(E∗))+Det(J(E∗))=0satisfies the Nash equilibrium pointE∗,it will undergo a flip bifurcation starting on period doubling cascade. Inequality (16) governs the stability of the dynamic game (8) when it is drawn for any pair of the system’s parameters in the plane. The dynamical behaviors, phase portrait, bifurcation diagrams, critical curves and basin of attraction of model (8) will be shown in the next numerical simulation.
5. Numerical Simulation of Dynamic Game (8)
Since the system (8) has many parameters, we study in this section the influences of maximum price parameter a, the speed of adjustment parameterν1and the sale constraint parameterμ1on the stability of Nash equilibrium point keeping the other parameters values fixed. This numerical experiments include using some popular graphs such as largest Lyapunov exponent (LLE), 1D bifurcation diagram, phase plane of chaotic attractors, and attraction basins for some periodic cycles. These graphs are important when studying the complex dynamic characteristics of such systems. We begin our analysis by investigating the influences of the sale’s constraint parameterμ1on the dynamics of the map (8). This requires to consider the following parameters set,a=11.25,b=0.5,c1=0.1,c2=0.3,μ2=0.6,ν1=0.2andν2=0.7. Since the fixed pointE∗depends on the parametersa,b,c1,c2,μ1andμ2so it will be changed as we study the influence ofμ1and its stability will be affected too. Takingμ1 as the bifurcation parameter and at this parameters set, Figure 1a shows the 1D bifurcation diagram when varying this parameter. Experiment simulations show that the two parametersν1andμ1have an opposite impact on the stability of the fixed pointE∗. Taking the other parameters’ values as previously and changeμ1to0.000215the fixed point becomesE∗=7.866637992,6.566681002and the Jacobian gets,
−0.57299−0.78649−0.91934−0.83867
whose eigenvalues areλ+≈0.15481andλ−≈−1.56647withλ+<1andλ−>1. This means thatE∗is unstable saddle point. The dynamic behavior at this set of parameters’ values includingμ1 is displayed in Figure 1b. It shows a two-piece chaotic attractor around the saddle fixed point. The light grey denotes the basin ofE∗which includes holes of divergent and infeasible points. Asμ1increases further the dynamic is converted into a four-piece chaotic attractor atμ1=0.01075 . Figure 1c plots the basin of attraction of those four unconnected chaotic areas. Asμ1increases those four chaotic areas change into a period-16 cycle. It has a quite complicated attraction basin at the same parameters set andμ1=0.03 . Figure 1e,f display the attraction basin of period-4 and period-2 cycles. They are born at the same parameters’ set and forμ1=0.108andμ1=0.89respectively. Keeping the other parameters’ values fixed and changeμ1=0.828the fixed point becomesE∗=7.224806201,6.887596899and the Jacobian gets,
0.75147−0.12427−0.96426−0.92853
whose eigenvalues areλ+≈0.819996andλ−≈−0.997057withλ+<1andλ−<1. This means thatE∗is local stable point. From this discussion we can see that there is an inverse period-doubling bifurcation that is merged when taking the sale’s constraintμ1as the bifurcation parameter. Then the game (8) undergoes a chaos zone, quasi-period range, and comes back to a local stable state with an increasing inμ1. Noting that the system (8) will return to stability from a chaotic state after many iterations of this system. Therefore the two competed firms should control the degree of minimum sales constraint parameter so as to avoid a market chaotic state. Furthermore the game (8) becomes more stable if the minimum sales constraint is considered by the players. On the other hand, the speed of adjustmentν1has given interesting results on the game’s behavior. Numerical simulation examples show that when firms adjust the speed of a certain value, the Nash equilibrium point of the system (8) will lose its stability, the system will become unpredictable, and chaos may occur. It’s difficult for firms to deal with such a complicated scenario. Consequently, the speed of adjustment to a very appropriate range must be taken into account by firms in order to avoid decreasing market efficiency and profits. We can conclude that the game (8) tends to exhibit more stability when at least one firm operates beyond its minimum sales constraint. Furthermore, the numerical simulation shows that the parameter a at the parameters’ set,b=0.5,c1=0.2,c2=0.3,μ1=μ2=0.6,ν1=0.5andν2=0.65 makes the fixed point unstable due to flip bifurcation. Taking a as the bifurcation parameter, Figure 1g shows that the fixed point becomes stable for all a until a reaches the value on where period-2 cycle coexists. For further increase in a higher period cycles are born and routes to chaotic attractors are obtained. Ata=12.11 and the other set of parameters are fixed a two-piece chaotic attractor is coexisted as shown in Figure 1h. It is clear that its basin of attraction lies in a quadrilateral area that will be discussed later on in this section.
We should highlight here that the system (8) can generate unbounded or negative trajectories if the initial conditionq0,1,q0,2are assumed to be taken far from the origin. Ifqi0>1+νi(1−μi)a−cibνi(1−μi),i=1,2the first iteration of the system (8) gives negative values forqi(t+1). Therefore, successive iterations yield negative and decreasing values for the system (8) asqi(t+1)=qi(t)+νi qi[(1−μi)(a−2bqi−bqj)−ci]<qibeing(1−μi)a−ci>0if (11) holds. This means that attractors at finite distances are impossible to be globally attracting inR+2because their attraction basins can not be extended out the rectangle0,((1−μ1)a−c1)/bν1(1−μ1)×0,((1−μ2)a−c2)/bν2(1−μ2). The map (8) has another important distinctive feature that is each axisqi(t)=0givesqi(t+1)=0,i=1,2that meansqi(t)is mapped into itself and hence it is trapped. Settingqi(t)=0(the monopoly case) in (8) we get
qj(t+1)=1+νj(1−μj)a−cjqj(t)−2bνj(1−μj)qj2(t),j=1,2,i≠j
The Equation (17) conjugates the standard logistic map,y(t+1)=λy(t)(1−y(t))through the linear transformation
qj(t)=1+νj(1−μj)a−cj2bνj(1−μj)y(t)
where,λ=1+νj(1−μj)a−cj. This implies that the dynamic of (17) is governed by the well-known dynamics of standard logistic map. The map (17) is characterized as a unimodel map anddqj(t+1)dqj(t)=0gives the unique critical pointC−1as follows
qjC−1 =1+νj(1−μj)a−cj4bνj(1−μj),j=1,2
which conjugates the critical pointy=12for the logistic map. There are two fixed points for the map (17) given by
qjEo =0,qjEj =(1−μj)a−cj2b(1−μj)
which conjugate the fixed points,y=0andy=νj(1−μj)a−cj1+νj(1−μj)a−cj=1−1λof the logistic map. It is easy to see thatdqj(t+1)dqj(t)qj=qjEo >1and henceqjEo is unstable repelling point. The second fixed point of (20) is a stable attracting point provided that0<νj(1−μj)a−cj<2. Atνj(1−μj)a−cj>2the dynamic behavior of (17) may be periodic cycle or cyclic chaotic behavior whose attraction basin is bounded by the unstable repelling pointqjEo and its preimageO−1jwhose coordinate is
qjO−1 =1+νj(1−μj)a−cj2bνj(1−μj)
which conjugates the pointy=1in the standard logistic map. This implies that trajectories of the map (14) starting from an initial point selected out of the interval[0,qjO−1 ]are divergent to−∞. Numerically, we assume the set of parameters’ values,a=10,b=0.3,c=0.3andμ=0.5 . Figure 1a shows thatqjEj is locally asymptotically stable for all the value ofνtill it approachesνj((1−μj)a−cj)=2on where period-doubling bifurcation coexists. Takingν as the bifurcation parameter of (17), Figure 2a shows thatqjEj is locally stable for all the value ofνexcept at the pointν=2(1−μ)a−cwhere period-2 cycle is born due to flip bifurcation. Asν increases further chaotic attractors exist. Figure 2b shows a stable fixed point at the parameters’ values,a=10,b=0.3,c=0.3,μ=0.5andν=0.4. It is clear that the basin of the stable fixed point lies within the interval[0,qjO−1 ]. Increasingνtoν=0.635 gives rise to chaotic attractor with basin of attraction also lies within the same interval as shown in Figure 2c. On the other hand, when we takeμ as the bifurcation parameter, Figure 2d shows the bifurcation diagram for the map (17) at the parameters’ values,a=6,b=0.5,c=0.3andν=0.5. Atμ=1−νc−2νaa period-2 cycle is born. Asμincreases further toμ=0.25 the fixed point gets stable as shown in Figure 2e. Forμ=0.018a chaotic attractor exists and its basin of attraction is also within the interval[0,qjO−1 ] as shown in Figure 2f.
Now we study the dynamics of (8) at different sets of parameters’ values. Settinga=11,b=0.5,c1=0.2,c2=0.3,μ1=0.6,μ2=0.6,ν1=0.33andν2=0.71we have only the stableE∗ with its basin of attraction in Figure 3a. The light blue color consists of an open set of the pointsq1,q2in the phase plane whose trajectoriesT(q1,q2)belong to the basin of attraction ofE∗(which is denoted byB(E∗)). The grey color color denotes to the basin of infinity that contains all the points generating unbounded trajectories and is denoted byB(∞) . This basin includes the divergent and unfeasible points. Figure 3b shows the basin of attraction of a chaotic attractor that coexists at the parameters’ values,a=11,b=0.5,c1=0.2,c2=0.3,μ1=0.6,μ2=0.6,ν1=0.715andν2=0.6 while Figure 3c displays the basin of attraction of the period-8 cycle. As one can see from all those figures that their basins are bounded by line segments separateB(κ)fromB(∞)whereκrefers to a bounded attracting set such asE∗, periodic cycles or any complex chaotic attractors aroundE∗. At the parameters seta=10.99,b=0.5,c1=0.47,c2=0.3,μ1=0.6,μ2=0.6,ν1=0.72andν2=0.69 we give in Figure 3d another quite complicated attractor around the fixed point, as one can see that it appears atc1>c2. Now we assume another set of parameters’ values that isa=11,b=0.4,c1=c2=0.33,μ1=μ2=0.3,andν1=ν2=0.28. This set indicates equality of the marginal costs, speed of adjustments and sales constraints. Numerical simulations show that at this set a period-2 cycle is emerged with a basin of attraction that is bounded by the same previous segments. As the speed of adjustment parameters, the sales constraints have some influences on the dynamics of (8). For this reason we decreaseμ1to0.0086 to get four chaotic attractors that are represented in Figure 3e with the basin of attraction. Increasingμ1further to0.001 we get two unconnected chaotic areas represented in Figure 2f with their basins. In conclusion, we have seen thatB(∞) is an unconnected set which sometimes has unconnected regions represented by holes in Figure 3b,d,f.
Now we intend to identify the structure of attraction basin for the map (8) and show the regions on where the map’s phase plane can be divided. Since the map (8) is trapped in the point(0,0)so it is important to calculate the rank-1 preimages for this point. The preimages of point (0,0) can be derived by algebraically solving the following system,
q1+ν1 q1[(1−μ1)(a−2bq1−bq2)−c1]=0,q2+ν2 q2[(1−μ2)(a−2bq2−bq1)−c2]=0,
which gives four preimages,O−1(0)=(0,0),O−1(1)=(q1O−1 ,0),O−1(2)=(0,q2O−1 ),andO−1(3)=(q^1,q^2)where,
q^1=a(1−μ1)(1−μ2)+c2(1−μ1)−2c1(1−μ2)3b(1−μ1)(1−μ2)+2ν2(1−μ2)−ν1(1−μ1)3bν1 ν2(1−μ1)(1−μ2),q^2=a(1−μ1)(1−μ2)+c1(1−μ2)−2c2(1−μ1)3b(1−μ1)(1−μ2)+2ν1(1−μ1)−ν2(1−μ2)3bν1 ν2(1−μ1)(1−μ2).
Or,
q^1=q1∗+2ν2(1−μ2)−ν1(1−μ1)3bν1 ν2(1−μ1)(1−μ2),q^2=q2∗+2ν1(1−μ1)−ν2(1−μ2)3bν1 ν2(1−μ1)(1−μ2).
Now, the preimages of any pointp,0on the invariant axisq2=0can be obtained by algebraic solution of the system,
q1+ν1 q1[(1−μ1)(a−2bq1−bq2)−c1]=p,q2+ν2 q2[(1−μ2)(a−2bq2−bq1)−c2]=0.
The second Equation of (25) givesq2=0or,
1+ν2[(1−μ2)a−c2]−bν2(1−μ2)q1−2bν2(1−μ2)q2=0.
Withq2=0we get the following quadratic equation,
−2bν1(1−μ1)q12+1+ν1a(1−μ1)−c1q1−p=0,
which has the following discriminant,
Δ=1+ν1a(1−μ1)−c12−8bpv1(1−μ1).
This means the pointp,0has on the same axis no rank-1 preimages, two distinct rank-1 preimages or two coincided rank-1 preimages ifΔ<0,Δ>0orΔ=0respectively. It may have also four preimages two of which lie on the same axis and the other two on the line (26). The same discussion and observation are for any point in the form0,p. The above discussion gives the following proposition.
Proposition 1.
Let the line segmentsξ1andξ2beξ1=O−1(0) O−1(1)andξ2=O−1(0) O−1(2)then
ℑ=⋃n=0∞T−n(ξ1)∪⋃n=0∞T−n(ξ2),
be the set of all preimages of rank-n.
Figure 4 shows that the line segmentsξ1,ξ2and their rank-1 preimagesξ1−1,ξ2−1constitute a quadrilateral regionO−1(0) O−1(1) O−1(3) O−1(2)that is exactly the attraction basin of the chaotic attractor given in the figure. Those segments and their rank-1 preimages construct the whole boundary of ℑ.
Now we put the map (8) in the following form,
q´1=q1+ν1 q1[(1−μ1)(a−2bq1−bq2)−c1],q´2=q2+ν2 q2[(1−μ2)(a−2bq2−bq1)−c2]
where′refers to the time evolution. Solving algebraically (29) with respect toq1andq2one can get0,2and 4 real preimages. It means that the map (29) is noninvertible and the phase plane may be divided into three zonesZ0,Z2andZ4. In order to identify these zones we should calculate the critical curvesLCandLC−1.LC−1is obtained by solvingdet(J)=0that gives the following hyperbola,
q12+q22+q1 q2+Δ1 q2+Δ2 q1+Δ3=0,Δ1=ν1 ν25(1−μ1)(1−μ2)a−c2(1−μ1)−4c1(1−μ2)+ν1(1−μ1)+4ν2(1−μ2)−4bν1 ν2(1−μ1)(1−μ2),Δ2=ν1 ν25(1−μ1)(1−μ2)a−4c2(1−μ1)−c1(1−μ2)+4ν1(1−μ1)+ν2(1−μ2)−4bν1 ν2(1−μ1)(1−μ2),Δ3=aν1 ν2(1−μ1)(1−μ2)a−c2(1−μ1)−c1(1−μ2)+1−c1 ν1−c2 ν2+c1 c2 ν1 ν2−4bν1 ν2(1−μ1)(1−μ2)
NowLCis obtained fromLC=T(LC−1). The complicated form of the hyperbola makes us to plotLC−1at the parameters’ set,a=11,b=0.5,c1=0.2,c2=0.3,μ1=μ2=0.6,ν1=0.33andν2=0.71 . Figure 5. displays the zonesZ0,Z2,Z4and the critical curves. The fixed point at this set of parameters’ values belongs toZ2.
6. Conclusions
In this paper, we have developed the Bischi–Naimzada duopoly game where competitors can maximize local profits under minimal sales constraints. We have studied and calculated the stability conditions of equilibrium points of this game. On the basis of nonlinear dynamical systems theory, we have examined the effects of critical system parameters on the stability and complexity of the dynamic game model, both in terms of sales constraints and speeds of adjustment. Using numerical analysis, the qualitative behaviors of the model dynamics have been explored by bifurcation diagrams, phase diagrams and largest Lyapunov exponent. In addition, the phase level of the map which was divided into three zonesZ0,Z2andZ4as a result of the noninvertible game has been analyzed. Research and analytical simulations show that the sensitive dependence of the dynamics of the model on both the speeds of adjustment and the parameters of the sales constraints. It also showed that, in contrast to the sales constraint parameters with an increase which has a stabilizing effect, the adjustment speeds have an effect of instability on the dynamics of model. In our contribution, it has been found that the Bischi–Naimzada competition game model that takes into account sales constraints and maximizes the profit function is much more stable around the Nash equilibrium than the same model that does not take into account sales constraints. In oligopolistic games, limited rational expectations are important for the investigation of market dynamics subject to minimal sales constraints. These games are of great importance for the study of competition between airlines, banking, music, drinks and telecommunications companies. This is for simple reason that is firms in previous areas always want to maximize profits with a minimum of sales. The extension of this paper is to explore of heterogeneous behaviors in such competitive Cournot games on the basis of minimal sales. More realistic conditions need to be taken into account, such as team participation, random selection of expectations, reduced risk, increase revenue and chaos control. Such important issues will be investigated in this direction for future studies.
Author Contributions
Both authors have contributed equally to this work. All authors have read and agreed to the published version of the manuscript.
Funding
Deanship of Scientific Research at king Saud University, Research group NO (RG-1438-046).
Acknowledgments
The authors would like to extend their sincere appreciation to the Deanship of Scientific Research at king Saud University for its funding this Research group NO (RG-1438-046).
Conflicts of Interest
The authors declare no conflict of interest.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2020. This work is licensed under http://creativecommons.org/licenses/by/3.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
The Bischi–Naimzada game is a market competition between two firms with the objective of maximizing profits under limited information. In this article, we study a more generalized and realistic situation that takes into account the sales constraints. we generalize the economic model suggested by Bischi–Naimzada by introducing and studying the maximization of profits based on sales constraints. Our motivation in this paper is the studying of profit and sales constraints maximization and their influences on the game’s dynamics. The local stability of the equilibrium points of the proposed model is discussed. It examines how the dynamics of the proposed two-dimensional competition game model focusing on changes in both the speed of the adjustment and the sales constraint parameters. The map describing the game is proven to be noninvertible and yields many multi-stable, complex dynamics and the coexistence chaotic attractors may arise. The global behavior of the map is achieved by studying the critical curves. The numerical simulations demonstrate the coexistence of two attractors and complex structures of the attraction basins. Several examples are discussed in order to confirm all the analytical results obtained.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer