Open access peer-reviewed chapter

Complex Dynamical Behavior of a Bounded Rational Duopoly Game with Consumer Surplus

By Wei Zhou and Tong Chu

Submitted: March 25th 2019Reviewed: June 3rd 2019Published: July 19th 2019

DOI: 10.5772/intechopen.87200

Downloaded: 447


In this chapter, we assume that two bounded rational firms not only pursue profit maximization but also take consumer surplus into account, so the objections of all the firms are combination of their profits and the consumer surplus. And then a dynamical duopoly Cournot model with bounded rationality is established. The existence and stability of the boundary equilibrium points and the Nash equilibrium of the model are discussed, respectively. And then the stability condition of the Nash equilibrium is given. The complex dynamical behavior of the system varies with parameters in the parameter space is studied by using the so-called 2D bifurcation diagram. The coexistence of multiple attractors is discussed through analyzing basins of attraction. It is found that not only two attractors can coexist, but also three or even four attractors may coexist in the established model. Then, the topological structure of basins of attraction and the global dynamics of the system are discussed through invertible map, critical curve and transverse Lyapunov exponent. At last, the synchronization phenomenon of the built model is studied.


  • bifurcation
  • chaos
  • duopoly
  • consumer surplus
  • synchronization

1. Introduction

Oligopoly is a market between perfect monopoly and perfect competition [1]. With the application of chaos theory and nonlinear dynamic system into oligopoly models, the static game evolves into a dynamic game. Especially in recent years, with the rapid development of computer technology, a powerful tool has been provided for dealing with the complex nonlinear problems. And hence, the economists and the mathematicians can simulate the complex dynamical behavior of oligopoly market by using computer technology. Recently, a large number of scholars have improved the oligopoly models, and introduced bounded rationality (see [2, 3]), incomplete information [4], time delay [5], market entering and entering barriers [6], differentiated products [7] and other factors [8, 9] into the classical oligopoly models, and the bifurcation and chaos phenomenon were founded in the process of repeated game.

However, all of the above discussions are mainly based on private enterprises, which pursuit the maximization of their own profits. In fact, the public ownership enterprises, which always aim at maximizing social welfare, and mixed ownership enterprises, which always aim at maximizing the weighted average of the social welfare and their own profits, are also widespread in the real economic environment. De Fraja and Delbono [10] found that the social welfare might be higher when a public ownership enterprise is a profit-maximizer rather than a social-welfare-maximizer. Matsumura [11] proposed that the social welfare could be improved through partial privatization of public enterprises. The research of Fujiwara [12] suggested that partial privatized public enterprises are more efficient than private enterprises. Elsadany and Awad [13] explored the complex dynamical behavior of competition between two partial public enterprises under the assumption of bounded rationality. However, the global dynamical behavior and synchronization behavior of semi-public enterprises, which corporate social responsibility into their objectives, are rarely studied. In this chapter, the occurrence of synchronization, the coexistence of attractors and the global dynamic of a duopoly game corporate consumer surplus are mainly discussed.


2. The model

Considering a duopolistic market where two firms produce homogeneous goods. In order to study the long-term behaviors of the duopoly market with quantity competition, we briefly present the economic setup leading to the final model in this chapter. The price and quantity of product of firm iare given by piand qirespectively, with i=1,2. We also assume the existence of a continuum of identical consumers which have preferences toward q1and q2.

Following Dixit [14] and Singh and Vives [15], we suppose that the utility function used in this chapter is quadratic and can be given by,


where q1, q2are the quantity of goods produced by firm 1 and firm 2, respectively. a>0represents the maximum price of a unit’s commodity, b>0represents the amount of its price decreases when the price of the product increases by one unit.

Suppose that the budget constraint of consumer is,


where p1and p2denote the prices of goods produced by firm 1 and firm 2, respectively. And Mdenotes the budget of the consumers on the product. The utility function of consumers is maximized under the budget constraint, and then the inverse demand function of the two firms can be obtained as,


This chapter discusses homogenous products, so here it is assumed that all these two players have the identical marginal cost. Therefore, the cost function of firm 1 and firm 2 are same and can be given by, Cq=cq, where c>0denotes the marginal cost of the goods and a>calways holds. Then, the profits of firm i, i=1,2can be obtained as follows,


In the real market, there are a lot of firms, who not only pursue their own profits but also take corporate social responsibility into account. A large number of empirical studies have shown how the introduced corporate social responsibility affects firm’s performance, where we interpret corporate social responsibility as either consumer surplus (for short CS) or social welfare (for short SW). In this chapter we take CS into account to analyze which firms have an incentive to exhibit corporate social responsibility as a means of maximizing their profits in a Cournot competition. Based on the above assumptions and the definition of consumer surplus, CS can be written as,


where p¯pais the price variable.

According to the above assumptions, the objective function of the firm ican be given as,


where αirepresents the weight of the consumer surplus in the objective function of firm i, and 0αi1always holds. By substituting (4) and (5) into (6), the objection function of firm ican be given by,


And the first-order condition of the objection function (7) is given as,


It is now significant to specify the information set of both players regarding the objection functions, to determine the behaviors of the players with the change of time. We assume a discrete time tZ+dynamic setting, where two firms with bounded rationality make decisions at the same time. That is, all firms do not have complete knowledge of their competitors’ decisions and the market demands. So it can only use the local estimation of the steepest slope of the objection function at period tto determine the output at period t+1. By following Bischi et al. [16] and Fanti et al. [17], the adjustment mechanism of quantities with the change of time of firm ican be obtained as,


where vi>0,i=1,2is an adjustment parameter of firm i. The firm iwill increase its output at period t+1, if Oitqit>0. But the firm iwill reduce its output at period t+1, if Oitqit>0. By substituting (8) into (9), we can get a two-dimensional map as,


Since the output of a firm cannot be negative, the initial conditions of map Tbelong to


By setting qit+1=qit,i=1,2in system (10), the fixed points of the system are obtained. Besides the trivial equilibrium E0=00, system (10) admits the following non-trivial fixed points (boundary equilibrium points),


and the only Nash equilibrium is


where A1=1α1, A2=1α2, B1=3α12, B2=3α22, C1=2α11and C2=2α21. The positivity of Eii=12and Eis ensured by requiring S=S1S2, where


3. Stability properties

The local stability analyses of system (10) near the fixed points are too difficult to carry on. For the sake of analyzing the local stability of the system, we firstly let α1=α2=αin system (10). And the Jacobian matrix of map Tat any fixed point q1q2can be given as,


Then all the equilibrium points are substituted into the Jacobian matrix (12). According to the characteristic values of the Jacobian matrix evaluated at each equilibrium, the type and stability of the equilibrium can be analyzed and the following results can be obtained.

Proposition 1. The equilibrium point E0is always an unstable node.

Proof. It is clear that the Jacobian matrix of map T, evaluated at the boundary equilibrium point E0can be written as,


The eigenvalues of JE0are given by λi=1+vi1αac,i=1,2. Since (11) holds and vi>0,i=12then λi>1,i=12which implies that E0is an unstable node.

Proposition 2.E1is a saddle point, when v2<21αac. And E1is an unstable node, when v2>21αac.

Proof. By substituting the equilibrium E1into (11), the Jacobian matrix of map Tevaluated at the boundary equilibrium point E1can be written as,


The eigenvalues of JE1are given by λ1=1+v11αacα13α2and λ2=1v21αac. Since (11) holds and v1>0, then λ1>1always meets. According to (11) and v2>0, we can deduce that λ2<1if and only if v2<21αac. That is, when v2<21αacmeets, the equilibrium E1is a saddle point. Similarly, the equilibrium E1is an unstable node, when v2>21αac.

Similar to the case of the equilibrium E1, it can be also proved that E2is a saddle point when v1<21αac, and E2is an unstable node when v1>21αac.

For the purpose of research of the local stability near the Nash equilibrium, we should compute the Jacobian matrix evaluated at the Nash equilibrium Eas,


It can be seen that the form of the Jacobian matrix is so complex. In order to simplify the calculation, let


Then the trace and the determinant of the Jacobian matrix evaluated at the Nash equilibrium Ecan be given as,


According to Jury condition, if we substitute the specific mathematical expressions of q1, q2into the above two equations, then the following set of inequalities can be gotten through a complex calculation,


Since all the equilibrium points are non-negative when the parameters meet 0α<35, a>cand vi>0,i=12. So we can get A>0, B<0, B+D<0and BD<0, then the set of inequalities (14) are equivalent to


Then the stability region of the Nash equilibrium can be obtained by substituting A, Band Dinto inequalities (15), which are given as,


The stability condition of the Nash equilibrium gives a parameters region, in which the Nash equilibrium is always stable. For the sake of better analysis of the stability of the Nash equilibrium under different set of parameters, a useful tool called “two-dimensional bifurcation diagram” (also called 2-D bifurcation diagram) is employed. From (16), we can find that the stability region of the Nash equilibrium is related to the difference of parameters aand c, that is ac. So we only discuss the values of ac, rather than the values of aand cin the rest of this chapter.

Figure 1 is a two-dimensional bifurcation diagram of system (10) with a set of fixed parameter ac, band α. Figure 1a is a two-dimensional bifurcation diagram when ac=2.1, α=0.50608821and b=2.00232504. We can observe that there are two different routes to chaos, when the parameters are chosen as this set of parameters. The system enters chaos through flip bifurcation when the parameter v1v2passes through green, yellow and light green from the brown region to the black region. It means that when the firms change their speed of adjustment according to the path, a periodic fluctuation of system (10) will happen. That is, the period motion will increase exponentially until it enters chaos. And the chaotic behavior of this system can be understood as the confusion of the market competition, and one of the two firms may be even out of the market with an increasing speed of adjustment. But if the parameter passes through the green region from the brown region to the black region, the system will first undergo a flip bifurcation, and then enters quasi-period motion through a Neimark-Sacker bifurcation. The system enters quasi-period from period-2 when firms determine their speed of adjustment along this path. Figure 1b is a partial amplification of Figure 1a, and we can observe that there are many scattered points of different colors, which is caused by the coexistence of multi-attractors with different period.

Figure 1.

2D bifurcation diagram in thev1v2orα1α2parameters plane. (a)ac=2.1,α=0.50608821,b=2.00232504, (b) enlarged the square region inFigure 1a, (c)ac=0.88017028,α=0.27445462,b=0.52714274 and (d) enlarged the square region in theFigure 1c.

Figure 2a shows the coexistence of attractors with the parameters chosen as v1=2.44, and v2=2.45, where the scatter points are shown in Figure 1b. we can observe a period-6 cycle coexisting with a period-4 cycle. Figure 1c is the two-dimensional bifurcation diagram, where the fixed parameters are given as ac=0.88017028, α=0.27445462and b=0.52714274. At this set of parameters, the system enters chaos through a flip bifurcation. Figure 1d is a partial enlargement of Figure 1c. Similarly, the parameter space of Figure 1d is also chosen according to the area with scattered points of Figure 1c. Figure 2b shows the coexisted attractors and their basins of attraction at this set of parameters. We can observe that there are three attractors coexisting.

Figure 2.

(a) A period-8 cycle coexists with a period-8 cycle withac=2.1,α=0.50608821,b=2.00232504,v1=2.44, and v2=2.45 and (b) a attractor on the diagonal and a period-4 cycle coexists with a period-8 cycle withac=0.88017028,α=0.27445462, v1=4.16,b=0.52714274, and v2=4.13.

Figure 3 shows a series of two-dimensional bifurcation diagram under different parameters. It shows a very beautiful gallery, from which we can enjoy the system (10) with full complex dynamics phenomenon. We can observe from Figure 3 that the difference between the maximum price of the unit product aand the marginal cost caffects the size of the stable region. The weight of the consumer surplus αaffects the shape of the two-dimensional bifurcation diagram, while the parameter balmost has hardly effect on the two-dimensional bifurcation of the system. Therefore, the game can be balanced more quickly by reducing the difference ac. In Figure 3a, it can be observed that the chaotic area surrounded by the period-4 area is like a “hand” and the 8-period area is like a small “bottle” raised by the beautiful hand. It is observed that the shape of Figure 3b is similar to Figure 3a due to a tiny adjustment of parameter α. Since the difference between the parameters aand cis reduced, the stability area of the Nash equilibrium becomes larger. Figure 3c is like a “volcanic eruption.” It can be observed that there is an inward cave in the diagonal that is like a “crater.” As the parameters vary, the hole in Figure 3c continues to sink inward and become larger and larger. From Figure 3d we can see that the period-4 arrives at a quasi-periodic motion directly. Therefore, we should change the weight of consumer surplus αslightly in order to maintain the market fluctuations not fierceness.

Figure 3.

2D bifurcation diagram in the (v1,v2) parameters plane, (a) with parametersac=1.83,α=0.278,b=0.0012, (b) with parametersac=1.52,α=0.277,b=0.348, (c) with parametersac=1.66,α=0.27,b=0.19 and (d) with parametersac=1.99,α=0.39,b=0.22.

4. Global dynamics and synchronization

The type and stability of the equilibrium points have been analyzed as above. And the boundary equilibrium E1and the boundary equilibrium E2are in symmetrical positions with respect to the main diagonal line Δ=q1q2R+2:q1=q2. It is also clear that the unique Nash equilibrium Eof system (10) is located on the main diagonal Δ. So we mainly study the dynamical behavior of the system on the diagonal. We choose the initial conditions near the diagonal, and the phenomenon via finite iteration back to the diagonal is called synchronization. The synchronization of chaotic systems was quite interesting and unexpected. In fact, due to the nonlinear system usually has sensitive dependence on initial conditions, a property which implies that the slightly change of initial conditions will lead to an exponential difference between the trajectories of two identical systems, making it impossible for two separated and even identical systems to synchronize. Therefore, the small coupling between two chaotic oscillators makes the system asymptotically converge to same trajectory, which is worth studying.

Subsequently, we assume that both firms have the same speed of adjustment. It means that the latter discussion is based on v1=v2=vand α1=α2=α. In this case the two players are identical; the system Tcan then be rewritten as follows,


It can be proved that the map Thas symmetry property, i.e., there exists a map S:q1q2q2q1, which makes TS=ST. The symmetry property of the map Timplies that the diagonal Δis a one-dimensional sub-manifold of system (17), i.e., TΔΔ. However, the phenomenon of synchronization occurs when the diagonal Δis invariant one-dimensional submanifold of the system (17). Therefore, the phenomenon of synchronization of the system can be analyzed by studying the invariant set. We can also use critical curve and noninvertible map to describe the global dynamical behaviors of a 2-dimensional map.

4.1 Critical curve and noninvertible map

We divide the discrete dynamical system into invertible and noninvertible. The invertible discrete dynamical system refers that an image q1q2of the map Tis correspond to the only preimage q1q2. The noninvertible discrete dynamical system implies that the map Tis multi-valued, i.e., the image of Thas more than one preimages. In a noninvertible discrete dynamical system, the curve that divides the phase space into regions with a different number of rank-1 preimages is called critical curve, denoted by LC. And the regions can be represented by Zi,iN. For example, a point belonging to area Z0has no preimage and a point belonging to area Z2has two preimages. Let us denote the rank-1 preimages of critical curve LCunder map Tas LC1. The set LCis the 2-dimensional generalization of the critical value or local extremes of 1-dimensional noninvertible map. Its preimages LC1are corresponding to local extreme point (critical point) in the one-dimensional noninvertible map. Since the map (17) is a continuously differentiable map, LC1belongs to the locus of points where the Jacobian determinant of Tvanishes, i.e., LC1=q1q2R2detDT=0. In this case, curve LC1can be determined by the following equation,


LCis the rank-1 image of LC1under map T. That is, LC=TLC1. The noninvertible properties play a significant role in analyzing the global behavior of a nonlinear discrete dynamical model. So the critical curve is a powerful tool for us to study these complex structures. Using the segment of critical curve as well as their preimages of any rank, and we will get the boundary of the basins of attraction as shown in Figure 10f.

4.2 Invariant sets

The dynamics of the system on the diagonal is studied by analyzing the invariant sets. Firstly, we can prove that the coordinates are invariant sets of map T. Let q2t=0, then we can obtain q2t+1=0, and the first equation of (17) can be rewritten as,


It is easy to verify that the dynamics on the axis q2is also controlled by the map (19). It means that the system (17) can be regarded as a 1-dimensional map at the coordinate axes. The map (19) is topologically conjugate to the standard logistic map xt+1=ωxt1xtthrough a linear transformation, which is given as,


and the parameter ωcan be presented as ω=1+v1αac. Thus the nonlinear dynamics of system (17) on the invariant axes can be analyzed through the standard logistic map.

It can also be proved that the diagonal Δis an invariant set of system (17), i.e., the trajectory starting from the diagonal Δwill stay forever on it. Therefore, the dynamical behavior of system (17) can be analyzed through the map Twhich is restricted to the diagonal. If we let q1=q2=q, then the dynamics generated by Ton the diagonal Δcan be analyzed through the following map,


Similarly, through the following linear transformation


we can also prove that the map (21) is topologically conjugate to the standard logistic map yt+1=μyt1yt, where


Through the standard logistic map, we can easily analyze the dynamical behavior of the two-dimensional map Ton the diagonal Δ. Under this situation, the Nash equilibrium Eof the system (10) is identical with the fixed point of map TΔ. Since μ=1+v1αac, we take different values of the bifurcation parameter μof the logistic map, and Figure 4a gives different bifurcation curves of the system on the parameter plane αv. The flip bifurcation occurs when the system parameter vequals v=21αac, and the Nash equilibrium Eloses its stability and forms a period-2 cycle around E. At v=61αac, the period-2 cycle generates a period-4 cycle after a flip bifurcation. When μ3.5699, the standard period doubling cascade ends and the system enters chaos. When v>31αac, the general trajectory of the map TΔis divergent.

Figure 4.

(a) Bifurcation curves on the parameter plane (α,v) related to the bifurcation values of mapTΔforac=2 and (b) 2D bifurcation diagram in the parameter plane (α,v) forac=2 andb=0.4.

As shown in Figure 4b, which the parameters is the same as Figure 4a, a two-dimensional bifurcation diagram of the system with vand αis obtained. Since it has been proved that the map TΔis topologically conjugate to the logistic map, we can find that the bifurcation curves of the two graphs are the same. In Figure 4a, the curve C1is correspond to the equation v=21αac, and the region below it represents the set of points of vand αat 1<μ<3. In this region, the fix point is stable. That is, the synchronization trajectory converges to the Nash equilibrium point. In Figure 4b, it corresponds to the period-1 region below the green region. When the point above the curve C2passes through the curve Cin Figure 4a, the system goes into chaos through a period doubling cascade. In Figure 4a, the curve Cis correspond to the equation v=31αac, where the parameters vand αof the upper area lead to the divergence of the trajectory, corresponding to the black area in Figure 4b. Through the above analysis, the following proposition can be derived,

Proposition 3. If we let v1=v2=v, the parameters a,band care fixed for system (17). Then, a threshold α=13vacof the weight of consumer surplus αor a threshold v=31αacof the speed of adjustment vdoes exist such that synchronized trajectories of the system (17) are divergent when α0αor vv+.

In order to analyze the effect of any slight perturbation of one parameter on the system, we study the transverse stability of an attractor Aof map T. And the Jacobian matrix of map Ton the diagonal can be obtained as follow,


Then, the characteristic values of the Jacobian matrix Jqqevaluated at any point on the diagonal are given by,


where the corresponding eigenvectors are 11and 11, respectively. And the eigenvalue λis related to the invariant manifolds on the diagonal.

It is assumed that a period-k cycle (q1q1)(q2q2)(qkqk)embedded into the invariant set Δof the map Tis correspond to the cycle q1q2qkof the map TΔwhen the synchronized phenomenon occurs, the two multipliers are given as,


Since the stability conditions of the period-k cycle on the diagonal Δof system (17) is same with the one-dimensional map TΔ, here we only study the transverse stability of the one-dimensional map TΔ. Under this situation, the transverse eigenvalue evaluated at the Nash equilibrium point Eis given by


Through Eq. (27), we can draw the following conclusions directly. That is, when all the parameters satisfy 0<vac1α2+5α<3, the Nash equilibrium Eis transversely attractive.

As we know that an attractor Aof Tis asymptotically stable if and only if all the trajectories that belong to attractor Aare transversely attractive. To study the stability of the attractor, we can calculate its transverse Lyapunov exponent as,


where q0Aand qiis a generic trajectory generated by the map TΔ. If the initial condition q0belongs to a period-k cycle, then Λ=lnλk. In this case, if Λ<0, then the period-k cycle is transversely stable. When the initial condition q0belongs to a generic aperiodic trajectory embedded in the chaotic attractors, then the transverse Lyapunov exponent Λis the natural transverse Lyapunov exponent Λnat. Since many unstable cycles along the diagonal are embedded in the chaotic attractor A, a spectrum of transverse Lyapunov exponents can be determined by the inequality


If all cycles embedded in Aare transversely stable (Λmax<0) then Ais asymptotically stable in the Lyapunov sense. If some cycles embedded in the chaotic attractor Aare transversely unstable (Λmax>0and Λnat<0) then Ais not stable in the Lyapunov sense, but it is a stable Milnor attractor. So we can look for the Milnor attractors by transverse Lyapunov exponents.

Figure 5 gives the natural transverse Lyapunov exponent and the bifurcation diagram with the fixed parameter vwhen ac=5.15, b=0.1911895and α=0.4. Later, we will exhibit the attractors and their basins of attraction corresponding to different values vunder this set of parameters, and analyze the changes of attractors and their basins of attraction when the parameter vvaries.

Figure 5.

Bifurcation diagram for the restriction of the mapTΔto the invariant diagonal and the corresponding transverse Lyapunov exponent forv[0.78, 0.91] givenac=5.15,b=0.1911895 andα=0.4.

4.3 Global bifurcation and basins of attraction

A closed invariant set Ais a attractor which means that it is asymptotically stable, i.e., a neighborhood Uof Adoes exist such that TUUand Tmq1q2Awhen m,q1q2A. We also define a asymptotically stable invariant set as attractor. A basin of attraction may contain one or more attractors that may coexist with a set of repel points that produce either intermittent chaos or a blurry boundary. The basin of attraction of attractor Ais the set of those initial conditions that cause the trajectory to converge to A, i.e.,


For the sake of analyzing the topological structure of the basin of attraction BA, we study the boundary of BAfirstly. Suppose that the map Thas a unique attractor Aat finite distance, let BAbe the boundary of the basin BA, then it is also the boundary of the basin of infinity Bgenerated by unbounded trajectories. Firstly, we take the dynamics of system (17) into account and restrict it to the invariant axis. When v1αac<3, if the initial conditions belong to the interval εi=001i,i=12, according to map (17), we can obtain the bounded trajectories along the invariant axes, where 01i,i=12is the rank-1 preimage of the origin. It has been obtained previously that the dynamical behavior of system (17) on the coordinate axis is governed by the map (19), so that 01i,i=12can be computed by the following algebraic system


The result is given as,


Since ε1and ε2are the segments of the boundary BA, and BAis also the boundary of the basin of infinity B, their rank-kpreimages Tkεi,i=12also belong to BA. We can compute the rank-1preimages of a point P=p0ε1or P=0pε2that belongs to εi=001i,i=1,2, according to the algebraic system as follows


We can easily obtain the rank-1 preimages of the origin, which are O11=0110and O12=0012, there are itself and O13besides, i.e.,


Through the discussion above, we can get the following propositions,

Proposition 4. Let 1<v1αac<3and εi=001i,i=1,2be the segments of the coordinate axes qi,i=1,2, then we can obtain the boundary of BAas follow,


Basins of attraction may be connected or not. The connected basins of attraction are divided into simple connected and complex one, and the complex connected basins of attraction means the existence of holes. If Ais a connected attractor, the direct basin of attraction D0of Ais the largest connected area of the entire attractor domain Dcontaining A. The system (17) has the coexistence of attractors in a set of given parameters, the basin of attraction Drefers to the union of the domain of attraction of all attractors in such a situation.

Figure 6 shows the coexistence of attractors and their basins of attraction for given parameters ac=5.3, b=0.234and v=0.85. In Figure 6a, the parameter αis chosen as α=0.4, there are two attractors coexisting, one is a Milnor attractor Alocated on the diagonal and the other consisting of 4-piece chaos attractor is in symmetrical positions with respect to the diagonal, i.e., F=i=14Fi. The basin of attraction is composed of the union of the attractive domain of two attractors. The attractive domain of the Milnor attractor Ais the complex connected set, and the attraction domain of the attractor Fis non-connected set. And the boundary of the 4-piece chaos attractor is just contact with the critical curve, and it is because of this contact that the system undergoes a global bifurcation. Figure 6b is the attractor and the attractive basin at α=0.387275after the global bifurcation occurred in the Figure 6a. We can find 4 attractors coexisted in this figure the attractor Fin Figure 6a undergoes a global bifurcation and turns into 3 period-4 cycles, and the attractor Ais also a Milnor attractor. The attraction domain of the period-4 cycle is composed of some complex connected sets being in symmetrical positions with respect to the diagonal. As is shown in Figure 6b, there are many holes in the attracting domain, but it is a non-connected set. The basin of attraction of the Milnor attractor lying on the diagonal is still a complex connected set.

Figure 6.

Basin of attractions for parametersac=5.3,b=0.234 andv=0.85. (a)α=0.4, a four-cyclic chaotic attractor coexists with a Milnor attractor on the diagonal and (b)α=0.387275, a four-cyclic chaotic attractor has undergone a global bifurcation and three period-4 attractors coexist with a Milnor attractor.

We have analyzed the global bifurcations that occur when the attractor’s boundary contact to the critical curve, and we discuss the global bifurcation when the attractor contacts to the boundary of its basin of attraction. We also denote global bifurcation as “boundary crisis,” the attractor is destroyed when it contacts to its basin of attraction. Figure 7 shows the coexistence of attractors and their basins of attraction corresponding to the parameter ac, when the parameters are chosen as α=0.5023335, b=0.4and v1=v2=0.85. We can see that as the difference between the maximum price aof a unit commodity and the marginal cost c increases, a period-4 cycle turns into 4-piece chaos attractor being in symmetrical positions with respect to the diagonal and finally merges into 2-piece chaos attractor. However, the 2-piece attractor being in symmetrical positions with respect to the diagonal grows larger as the parameter acincreasing, until it contacts to its basin’s boundary, and eventually occurs a global bifurcation, causing itself and its basin are destroyed until it disappears. We can also see its “ghost” in Figure 7d. This means that trajectories of the initial conditions that belong to the basin of attraction spend a long number of steps in the region occupied by the former attractor before converging to the other attractor. Figure 8 is the bifurcation diagram of the system at this set of parameters, and the bifurcation parameter is chosen as ac.

Figure 7.

Basin of attractions for parametersb=0.4,α=0.502335 andv=0.85 (a)ac=6, a two-piece chaotic attractor coexists with a period-4 cycle, (b)ac=6.1, the 4-piece chaotic attractor coexist with a two-piece chaotic attractor, (c)ac=6.1896, the 2-piece chaotic attractor formed by the a 4-piece chaotic attractor coexist with a two-piece chaotic attractor, and the many holes created by the global bifurcation and (d)ac=6.26, a two-piece chaotic has a contact with its basin’s boundary, and it is destroyed.

Figure 8.

One-dimensional bifurcation diagram with respect toacfor the set of parameters inFigure 7.

Figures 6 and 7 give two different global bifurcations, such bifurcations which can be restored clearly by numerical simulation method only. With the set of parameters in Figure 9 being identical to Figure 4, we select different speed of adjustment to analyze the change of attractors and their basins of attraction. We can observe that as the speed of adjustment changes from 0.79 to 0.9, the period-4 cycle being in symmetrical positions with respect to the diagonal of the system generates smooth limit cycle via a Neimark-Sacker bifurcation, as shown in Figure 9b, and the limit cycle becomes non-smooth gradually, and finally forms four-piece chaotic attractor, as shown in Figure 9d. The basin of attraction shrinks as the speed of adjustment vincreasing. It is implied that when both firms choose a lower speed of adjustment, they can reach the balance easily in the game. However, the period-2 cycle embedded in the diagonal becomes a period-4 cycle, period-8 cycle, etc. That is, a flip-bifurcation happens. And finally a Milnor attractor forms with the increasing speed of adjustment. Its basin of attraction increases with the increasing speed of adjustment gradually. In the bifurcation diagram of Figure 5b, we can observe the process of entire bifurcation process.

Figure 9.

Basin of attractions for parameterac=5.15,b=0.1911895 andα=0.4. (a)v=0.79, an attracting four-period cycle coexists with the two-period cycle on the diagonal, (b)v=0.83, a period-4 cycle has undergone a Neimark-Sacker bifurcation and an attractor formed by four smooth curves coexist with a period-8 cycle, (c)v=0.872, a four-cyclic chaotic attractor formed by the four smooth curves coexists with the Milnor attractor and (d)v=0.9, a four-piece chaotic attractor exists outside the diagonal.

4.4 Synchronization

In this section we study the formation mechanism of the synchronization trajectories. The trajectories starting from different initial conditions return to the diagonal eventually, i.e., q10q20. A tdoes exist such that q1t=q2twhen t>t, and we define such trajectories as synchronization. However, when the diagonal Δis an invariant sub-manifold, synchronized dynamics occur. We have proved that the map Tcan be obtained by two identical one-dimensional coupling maps, and the synchronization trajectory can be controlled by a map TΔwhich is topologically conjugate to the standard logistic map. When we choose the parameters as v=1, ac=5and b=0.234, the weight αvaries, and we can observe that the dynamic behavior of system is controlled by the attractor on the diagonal in Figure 10. When α=0.4837015, we can observe a Milnor attractor in Figure 10a. This means that cycle embedded in the diagonal are transversely unstable and blowout phenomenon occurs when the trajectory is near diagonal. The trajectory converges to the unique Milnor attractor embedded in the diagonal after experiencing a long transient. Figure 10b shows that the evolution of q1q2versus time and synchronization is observed after a long transient. This is a typical on-off intermittency phenomenon. We can observe 16-piece chaos attractor being in symmetrical positions with respect to the diagonal in Figure 10c, when αincreases to 0.485092. Figure 10d shows a chaotic attractor when αdecreases to 0.437609, then synchronization occurs. As shown in Figure 10f, we adopt the trial-and-error method, with suitable part of LC1taken as the starting part of Γ=ALC1to obtain the boundary of the chaotic attractor Aand the entire basin of attraction in Figure 10e, i.e., the boundary of the chaotic attractor Ais A=k=16TkΓ.

Figure 10.

Parameter values are chosen asv=1,ac=5 andb=0.234. (a) Four-piece Milnor attractor of system T belonging to the diagonal forα=0.48365, (b) the displacement q1 − q2 versus time for the same parameters as in (a), (c)α=0.485092, a 16-cyclic chaotic attractor is in symmetrical positions with respect to the diagonal, (d)α=0.437609, a trajectory in the phase space (q1, q2) whose transient part is out of diagonal that synchronizes along the Milnor attractor in the long run, (e)α=0.476955, a two-cyclic chaotic attractor coexists with a period-2 cycle and (f) boundary of the chaotic area obtained byA=k=16TkΓ.

5. Conclusion

In this chapter, the nonlinear dynamics of a Cournot duopoly game with bounded rationality is investigated. Unlike the existing literature, we suppose that the two firms not only pursue profit maximization but also take consumer surplus into account. Meanwhile, the objection of firms is supposed as the weighted sum of profit and consumer surplus. Based on the theory of gradient adjustment, all the firms adjust the output of next period according to the estimation of “marginal goal.” The existence and stability of fixed points are analyzed. It is found that the boundary equilibrium point is always unstable, no matter what the parameters of the system are satisfied. At the same time, with the two-dimensional bifurcation diagram as the tool, the stability of the Nash equilibrium is analyzed. We found that the Nash equilibrium will lose its stability when the speed of adjustment of firms is too large, which maybe lead the market into chaos. The stability region of the Nash equilibrium will be only affected by the weight of consumer surplus. And the parameters acand bhave hardly effect on the stability region of the Nash equilibrium. Meanwhile, we found that the two-dimensional bifurcation diagram have a beautiful fractal structure, but there are also many scattered points which is due to the coexistence of multiple attractors of the system through numerical simulation. By selecting corresponding parameters in the two-dimensional bifurcation diagram with scattered points, we draw the corresponding basin of attraction, and found the model not only has two attractors coexistence phenomenon, but also has 3, or even 4 attractors coexistence phenomenon.

Moreover, with the theory of invertible mapping and the critical curves of the system, the topological structure of basin of attraction is analyzed. By calculating the transverse Lyapunov exponent, the weak chaotic attractor of the system in the sense of Milnor is found, and the synchronization of the system is further studied. If we fix the other parameters of the system, and only change the weight of the firm to the consumer surplus, we can find on-off intermittency phenomenon and synchronization phenomenon. With the increasing of α, the synchronization phenomenon is vanished and a 16-piece chaotic attractor being in symmetrical position with respect to the main diagonal is produced. Under another set of parameters, and the parameter αis chosen as the bifurcation parameter. Through numerical simulation, it can be found that when the critical curve contact with the boundary of the basin, a global bifurcation is obtained. The global bifurcation makes the basin of attraction of the attractor non connected. In addition, if we fixed parameters of the system, and change the values of the parameters aconly, we find another global bifurcation called “boundary crisis,” i.e., when the attractor contact with its boundary of the basin of attraction, one of the attractors and its basin of attraction will be destroyed.

© 2019 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Wei Zhou and Tong Chu (July 19th 2019). Complex Dynamical Behavior of a Bounded Rational Duopoly Game with Consumer Surplus, Research Advances in Chaos Theory, Paul Bracken, IntechOpen, DOI: 10.5772/intechopen.87200. Available from:

chapter statistics

447total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Nonlinear Dynamical Regimes and Control of Turbulence through the Complex Ginzburg-Landau Equation

By Joël Bruno Gonpe Tafo, Laurent Nana, Conrad Bertrand Tabi and Timoléon Crépin Kofané

Related Book

First chapter

Classical and Quantum Conjugate Dynamics – The Interplay Between Conjugate Variables

By Gabino Torres-Vega

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us