Dynamic Analysis and Chaos Control of Bertrand Triopoly Based on Differentiated Products and Heterogeneous Expectations

Price competition has become a universal commercial phenomenon nowadays. ,is paper considers a dynamic Bertrand price game model, in which enterprises have heterogeneous expectations. By the stability theory of the dynamic behavior of the Bertrand price game model, the instability of the boundary equilibrium point and the stability condition of the internal equilibrium point are obtained. Furthermore, bifurcation diagram, basin of attraction, and critical curve are introduced to investigate the dynamic behavior of this game. Numerical analysis shows that the change ofmodel parameters in a dynamic system has a significant impact on the stability of the system and can even lead to complex dynamic behaviors in the evolution of the entire economic system. ,is kind of complex dynamic behavior will cause certain damage to the stability of the whole economic system, causing the market to fall into a chaotic state, which is manifested as a kind of market disorder competition, which is very unfavorable to the stability of the economic system.,erefore, the chaotic behavior of the dynamical system is controlled by timedelay feedback control and the numerical analysis shows that the effective control of the dynamical system can be unstable behavior and the rapid recovery of the market can be stable and orderly.


Introduction
e oligopoly theory is based on two classic economic models developed by Cournot [1] and Bertrand [2]. Cournot competition is named after French economist Antoine Augustin Cournot. In 1838, Cournot argued that decisionmaking in the oligopoly market would hinge on qualities in Cournot competition, as a firm has the market power over quantities but is aware of rivals' actions. us, the firm can deploy its strategies based on the symmetry information it has and finally maximize its profit. In 1883, Joseph Bertrand developed a game theory model that concerns price competition in a duopoly by employing the concepts of Nash equilibrium with market prices as a strategic factor. erefore, there is a distinct difference between Bertrand's and Cournot's economic models in terms of a firm's pricing strategy under the perfect competition, namely, pricing of marginal cost. Overall, the two economic models have made great contributions to the oligopoly markets, which are the cornerstone of the subsequent improvements and developments.
Bounded rationality and Puu's imperfect information are the two main frameworks for studying the oligopoly market. Enterprises with bounded rationality would update their own strategies based on the discrete time series and marginal benefit. Under such adjustment mechanism, enterprises do not need to have a comprehensive understanding of demand and cost functions. Instead, all they need is to decide whether they need to react with output's tiny change based on the estimated marginal benefit. e adjustment mechanism related to discrete time series has been extensively studied by numerous researchers.
Meanwhile, some researchers believe that the decisionmaking progress of discrete time series is more down to earth than the tangible economic system, which is less likely to revise the output decision real time. Elsewhere, Puu has introduced the idea of the so-called "Puu's imperfect information" with major benefit as being close to reality. An enterprise does not need to comprehend its profit function. Instead, it only needs to understand its previous profit and prices. is provides an even more convenient method to analyze monopoly competition.
Both Naїve and Bertrand assume an enterprise uses the latest available information. Each enterprise believes its rival provides the same quantity or price and the deviation did not exist in the past period of time. Under the static background, we can conclude that the monopoly is a state of global stable balance. However, based on the rule of expectation, the marginal cost of linear demand function and constant as well as the increasing number of firms play important roles in the Bertrand duopoly. Particularly, if there are more than three firms in the market, the Bertrand-Nash balance would become instable. e paradox would be generalized y changing different hypothesis, via taking different demand function [3,4] and cost function [5,6] into account. us, the research will start from the perspective of bounded rationality and imperfect information in order to transform the Bertrand duopoly from the static to dynamic forms. Furthermore, we will provide an effective chaotic control regarding the complexity of the model in order to explain how a monopoly market transforms from chaos to stability.
In an oligopoly market, each firm needs to take its and rivals' strategies into account. us, we see dynamics in an oligopoly competition. Furthermore, the research on linear dynamic system in oligopoly has generated the expectation rules that are more down to earth than those in Naïve expectation. In other words, we have seen expectation rules based on marginal benefit [7][8][9] and self-adaption [10]. In most of the literature, researchers focus on the competitive heterogeneity. Firms adopt the same set of expectation rules. However, a firm's isomerism expectation rules, such as competitive duopoly anticipated by isomers, should be closer to reality [11][12][13]. Recently, Fanti and Gori also analyze the dynamic model of Cournot competition [14]. Other literature also concern the analysis of competition between two firms, such as the research on three firms in oligopoly [15][16][17]. Currently, there is limited research on multiple heterogeneity firms in oligopoly and dynamic competition model. e study of nonlinear difference equations by Elsadany et al. believes that four heterogeneity enterprises would produce perfect substitute products in the gaming in an oligopoly. ey state that the stability from Cournot-Nash equilibrium hinges on the adjustment speed from gradient firms [18][19][20][21][22]. Meanwhile, adaptive expectant could stabilize the process of gaming. Elsewhere, Tramontana's research [19,[23][24][25][26] indicates that inverse demand function of elasticity, constant marginal cost, and growing heterogeneity firms would boost the range of stable region's parameter. Andaluz and Jarne have conducted a research on different products from three heterogeneity firms. Among them, one engages in price and quantity competition based on a different linear dynamic system in oligopoly and a firm's Naїve expectation. e other two firms revise their own judgement based on a gradient criterion. e research indicates that Nash equilibrium would increase stability if the product is more independent given there are complements and substitutions for the product.
Since the beginning of the twentieth century, we have witnessed accelerated development of study on the complexity of chaos in natural science. Currently, the study has yielded abundant results. anks to the cross-disciplinary sturdy, the study on the complexity of chaos for social science and business management has been gradually understood. Many scholars have attempted to broaden their understanding about the complexity of chaos, control, and application of the simplest economic and management model. More recently, some researchers have begun to use dynamic game method to study realistic management questions, such as service value, channel conflict, social influence, and spillover effect [27][28][29][30][31].
In this paper, the research will extend the models based on the previous research made by Andaluz and Jarne [32][33][34][35][36]. We will attempt to construct a nonlinear dynamic Bertrand model based on the heterogeneity in expectations (namely, Naïve expectations and bounded rationality) and make use of the advanced tool of complex system theory to extend the analysis of coopetition relationship by addressing the following research questions: (i) We analyze the local stability of the Nash equilibrium under price competition but assuming a nonlinear dynamic system of a differentiated triopoly. (ii) We have constructed nonlinear dynamic Bertrand model based on differentiated products and heterogeneous expectations and analyze the Nash equilibrium's partial stability. is will be based on the hypothesis that inverse demand function and cost functions are both linear functions along with participants' bounded rationality. By employing the stability of discrete dynamical system, we analyze boundary equilibrium of three dimensions' dynamical system. In the meantime, we provide the stable parameter of interior equilibrium based on Schur-Cohn principles. (iii) e equilibrium stability of nonlinear dynamic Bertrand model would come from numerical simulation. e result of analyzing nonlinear dynamic Bertrand model's equilibrium stability via numerical simulation would discuss the parameters' impact on the dynamic system's stability and system stability. e time-delayed feedback control could effectively control the system's chaotic behavior [37,38].
In order to solve these problems, we consider the price of the company's product as a decision variable. is paper considers both the market price function and the production cost function with knowledge absorption capabilities being linear functions, and a dynamic Bertrand price with heterogeneous expectations for participating companies is constructed. A game model and and a dynamic repetitive game model are proposed under bounded rationality. To the best of our knowledge, this article is the first to use analytical and numerical tools to consider duopoly repetitive games in dynamic Bertrand prices under conditions of limited competition. In particular, we adopt a gradient adjustment 2 Discrete Dynamics in Nature and Society mechanism in economics literature to capture players' bounded rationality [16]. is paper is organized as follows. In Section 2, we construct the production of three enterprises with differentiated products and independent competitive pricing in an economic monopoly industry. We theoretically analyze the system through the theory of complexity in Section 3. In Section 4, we use numerical simulation to demonstrate the various dynamic behaviors of the three-dimensional system and by time-delay feedback control to control the chaotic behavior of dynamical systems and numerical analysis show that effective and efficient control of power system can be unstable behavior. Finally, the main conclusions of this study are given.

Model Specification
is paper considers the production of three enterprises with differentiated products and independent competitive pricing in an economic monopoly industry. q i denotes the output of the enterprise i,p i denotes the product price of the enterprise i, and the consumer utility function is of the form where u(q 1 , q 2 , q 3 ) represents the second strict concave function. During the research, we use the following form: To simplify calculations, let By equations (1) and (2), the inverse requirement function can be obtained as follows: From equation (3), the corresponding requirements can be obtained as follows: , , .
From equation (4), we can know that when d � 0, the products are the most independent of each other and each enterprise is a monopolist in the market relative to other products; when − (1/2) < d < 0, it indicates that the product is a complementary product; when 0 < d < 1, it indicates that the products are substitutes for each other; when d � 1, it means the exact same product.
Assume that the linear production cost function of enterprise i is as follows: where c i > 0 represents the parameters related to the production technology level of the enterprise. e higher the level of production technology is the smaller c i will be. r i ≥ 0 represents the parameter of the enterprise knowledge absorption capacity effect, which indicates the absorption of external knowledge and technology by the enterprise i and thus the positive cost externality of the enterprise i. erefore, the profit that can be obtained from the enterprise i can be expressed as follows: Discrete Dynamics in Nature and Society en the marginal profit of enterprise i in the t period is as follows: Now, suppose that the decision rule of enterprise 1 is based on marginal profit, the decision rule of enterprise 2 is to follow the expected value of Naїve, and the decision rule of enterprise 3 is to adopt the maximized expected value of Naїve under the adaptive expectation plan. So, how to adjust the price of enterprise 1 mainly depends on whether the current marginal profit is positive or negative in t + 1 the period. erefore, enterprise 1 adjusts the decision mechanism over time as follows: where α(p 1 (t)) � αp 1 (t)(α > 0) indicates the speed at which enterprise 1 adjusts the price based on marginal profit. e decision rules for the enterprise 2 to adjust product prices over time are as follows: Substituting equation (6) into equation (9) yields the maximum expected value of Naїve in the following form: Finally, suppose that enterprise 3 adopts the adaptive expectation decision rule to adjust the product price (i.e., when it determines the price of the product in the next period, it is based on its own price decision in the previous period and the Naїve expected value decision in the previous period); it can be expressed in the following form: where 0 < β < 1 represents the coefficient of adaptive adjustment. If β � 0, then the level of the price remains permanent; when β � 1 it indicates the price level is the maximum expected value of Naїve.
From equations (8), (10), and (11), we can obtain the following dynamic system: A nonlinear discrete dynamical system (12) is obtained based on the assumption that the market inverse demand function and the production cost function are both linear forms. System (12) describes a price competition game relationship between three oligarchic enterprises with product differentiation and heterogeneous expectations. Under the condition of bounded rationality, the game subject will adjust the investment strategy of each period according to the marginal profit of the enterprise.

Local Stability of Equilibrium Points
In this section, firstly, the equilibrium points of the nonlinear discrete system (12) are calculated. en, through qualitative en, the stability characteristic equilibrium point of system (12) is analyzed qualitatively. Let p i (t + 1) � p i (t), i � 1, 2, 3.
erefore, system (12) can be converted into the following form: By solving equation (13), the two equilibrium points of the complex dynamic system (12) are obtained as follows: where It is easy to know that E 1 is the boundary equilibrium point and E 2 is the inner equilibrium point. Consider the realistic economic significance of system (12) equilibrium point representation.
erefore, only the case where the equilibrium point is nonnegative is discussed in the paper.
Since a, c i , and r i are all positive parameters, both E 1 and E 2 are greater than zero.
Since the local stability of the equilibrium solution of the complex dynamic system (12) depends on the eigenvalues of the Jacobian matrix, first calculate the Jacobian matrix corresponding to system (12) as follows: Discrete Dynamics in Nature and Society Lemma 1 (cf. [39]): in n-dimensional discrete dynamical system (17), the system equilibrium point x * is stable under the condition that all eigenvalues of the Jacobian matrix J(x * ) of the right-end function of the dynamic system (17) satisfy |λ n | < 1.
Next, the stability of the two equilibrium points of the complex dynamical system (12) will be further analyzed based on Lemma 1.

Proposition 1.
e boundary equilibrium point E 1 is unstable (saddle point).
Proof: To prove this result, we calculate the Jacobian matrix (14) at the boundary equilibrium point E 1 as follows: e eigenvalues of the Jacobian matrix J(E 1 ) are obtained by calculation as follows: From the previous assumptions, we can know the parameters α > 0 and − (1/2) < d < 1, so we can derive λ 1 > 1. erefore, according to Lemma 1, the discrete complex dynamic system (12) Nash equilibrium stability determination condition can be known that E 1 is an unstable boundary equilibrium. In addition, due to 0 < β < 1, we can know |λ 2,3 | < 1.
e local stability of the Nash equilibrium solution is discussed below. e Jacobian matrix (16) at the inner point equilibrium solution E 2 is expressed as follows: So, the characteristic polynomial of J(O) is as follows.
Let the characteristic polynomial of the matrix J(E 2 ) be G(λ), and the specific form of G(λ) is as follows: G(λ) � λ 3 + g 1 λ 2 + g 2 λ 1 + g 3 . en, by calculation, the coefficients of the polynomial G can be obtained as follows: It can be known from the Schur-Cohn stability criterion that if the characteristic polynomial G(λ), that is, all the eigenvalues of the Jacobian matrix J(E 2 ), is located in the unit circle on the complex plane, the coefficients of the 6 Discrete Dynamics in Nature and Society polynomial need to satisfy the conditions in the Schur-Cohn stability decision. erefore, the system parameters are substituted into it and simplified: e Bertrand-Nash equilibrium solution E 2 of the complex dynamical system (12) is locally asymptotically stable if there is Proof. It is known from the positive and negative assumptions of the system parameters that G(1) > 0 and 3 − g 2 > 0 are clearly established in stable conditions.
In addition, we can also know The white area in Figure 1 is shown, ∀α > 0, The red area in Figure 1 is shown.
When the parameter (d, β) is in the white area, the thresholds α 1 and α 2 are compared, and α 1 < α 2 is known.
From the above conclusions, we can conclude that the Nash equilibrium is locally asymptotically stable by adjusting the threshold of the following speeds: e local stability of Nash equilibrium is inferred from this proposition which largely depends on the adjustment speed of enterprise product price decision. If the company is overresponsive (i.e., the value of α is high), the Nash Discrete Dynamics in Nature and Society 7 equilibrium loses stability and complex dynamic behavior may occur throughout the system. erefore, it is also necessary to analyze the influence of other parameters on the threshold α C to see if there are other more interesting dynamic behaviors. en, we will discuss the effect of β and d on the threshold α C .
For (22), the partial derivative is obtained for β. e result is as follows: ere is a hypothesis that β � 1 represents the Nash expected value of the product price level. en, the higher the value of the product price adjustment speed parameter β is, the more the decision-makers of enterprise 3 tend to change the product price of the previous period, which leads to the phenomenon of system instability. Now discuss the impact of product differentiation, and obtain the following equation by calculation:  ) .

(28)
In Figure 2, (b) is magnification of (a) in the 0.16 ≤ d ≤ 0.22 interval. As can be seen from Figure 2, the curve (22) divides the parameter plane into two regions: the blue region represents zα C /zd > 0 and the pink region represents zα C /zd < 0.
In the complementary commodity (d < 0), the complementarity of the product is reduced (i.e., d tends to 0), and the Nash equilibrium is increased. On the contrary, this monotonic relationship does not appear in the substitute product. In this case, the threshold for adjusting the speed is the largest and has a certain degree of substitutability (for the fixed β); if the degree of difference between the products is increased, the stability of the Nash equilibrium is difficult to be guaranteed (note: d is true when the value approaches 0 on the right side). In the mutual substitution, this result is related to the complementary intertemporal strategy of the product.

Numerical Simulation and Chaos Control
In order to better demonstrate the complex dynamic evolution behavior of the nonlinear dynamic system (12), the numerical simulation of the system will be carried out. We will give the evolution trajectory diagram, bifurcation diagram, attractor chaotic diagram, system sensitivity dependence of the nonlinear dynamic system, and so forth. Because the chaotic behavior is irregular and disordered, its trajectory cannot be judged. In addition, the generation of chaos will also have an unpredictable impact on the management system, so this Discrete Dynamics in Nature and Society section uses the feedback control method to control the nonlinear dynamic system and gives theoretical proof.

Numerical Simulation.
In this section, the dynamic evolution of the nonlinear dynamical system (12) is numerically simulated, and the dynamic evolution process is visualized by simulation results such as bifurcation diagram, stable region, singular attractor, and the speed at which the system converges to the equilibrium point. Convenient for comparative analysis, based on the previous parameter assumptions, some of the system parameter values in this section are determined as follows: d � 0.5 and β � 0.5. Figure 3 shows that as the adjustment speed parameter α of enterprise 1 increases, the oscillation of the Nash equilibrium of system (12) causes more and more complex dynamic behaviors. e attractor in Figure 3(a) is Nash equilibrium, and Figures 3(b) and 3(c) show the complex behavior exhibited by the system when the period-doubling bifurcation leads to chaos. Figure 4 shows the phase diagram of singular attractor corresponding to Figure 3(c). Figure 5 shows bifurcation diagram in which the Nash equilibrium of the nonlinear dynamic system (12) loses stability and the complexity of the singular attractors increases as the adjustment speed α of enterprise 1 increases. When the value of α gradually increases from 0, the Nash equilibrium point of the nonlinear dynamic system (12) begins to become unstable at α � 0.4831 and enters the chaotic state after bifurcation. When α > 0.65, the evolution of the nonlinear dynamic system (12) collapses due to negative price. From the economic point of view, participants can withdraw from the commodity market for explanation.
In Figure 6, the effect of the parameter d on the Nash equilibrium of the nonlinear dynamic system (12) is illustrated by simulation in an alternative commodity. e greater the difference between commodities is (i.e., d approaches 0), the small the competition between market products is, and the Nash equilibrium point of the nonlinear dynamic system (12) is more stable. On the contrary, when the product differentiation is reduced (i.e., d approaches 1), the Nash equilibrium point of the nonlinear dynamic system (12) becomes unstable, which also indicates that the product differentiation is small, resulting in more intense market competition and unstability. It can be explained that increasing the degree of product differentiation makes the enterprise to improve its own competitive advantage, and increasing product differentiation requires enterprises to have more investment in technology and innovation. Product differentiation is also an important factor in regulating market stability. In addition, when the enterprise responds faster to the market under the same conditions, it will increase the market instability. If individuals adjust their own decisions and get more profits in the enterprise group,   other individuals will compare and follow the learning to adjust their strategies, which will lead to the unstable phenomenon of the market. Figure 7 shows bifurcation diagram in which the Nash equilibrium of the nonlinear dynamic system (12) loses stability and the complexity of the singular attractors increases as the difference d of product changes. When the value of α gradually increases from 0.5, the Nash equilibrium point of the nonlinear dynamic system (12) begins to become unstable at d � 0.5755. With the appearance of the bifurcation phenomenon, after the period-doubling bifurcation, it finally evolves into a chaotic state. Figure 8 shows the iterative evolution of the state variable orbit over time for the dynamical system (12) under two slightly different initial conditions. It can be seen from Figure 8 that as the number of iterations increases, the orbits under the two initial values of the same state variable are gradually separated, and the nonlinear dynamic system (12) shows significant sensitivity dependence on the initial value.

Chaos Control.
From the numerical simulation, we can know that the adjustment speed has a great influence on the stability of the nonlinear dynamic system (12) [40][41][42]. e dynamic behavior of system (12) will be more complicated if the model parameters cannot be located to the desired stable region. In the real economic system, chaos is unpredictable and unpredictable, so it needs to be avoided or controlled to make the dynamic system run better and better. Next, the control method is introduced by the time-delay feedback control to control the chaotic behavior of system (12). First, system (12) is embedded in the controller k(q 1 (t) − q 1 (t + 1)), where K > 0 is the control parameter. en the controlled system (12) changes to the following form: e Jacobian matrix of the controlled system (29) is as follows: It can be seen from Figure 5 that when enterprise 1 adjusts the parameter α > 0.65, the nonlinear dynamic system (12) starts to become unstable and accompanied the occurrence of chaotic behavior, so in the control system, let α � 0.7. Substitute the parameters into the Jacobian matrix and calculate the following form: e eigenvalue of the entire matrix is less than 1, which indicates that the set of parameter values of the dynamic system (29) is stable, so the chaotic control of the control system (12) tends to stable orbit. As can be seen from the stability conditions, the eigenvalues of the Jacobian matrix need to satisfy |λ i | < 1, (i � 1, 2, 3), thus solving k > 0.4025. erefore, when k increases, the controlled system (29) will gradually reach stable state.
In Figure 9, with the increase of the control parameter k, system (29) gradually tends to be stable from chaos and quasi-period. When k > 0.4025, system (29) finally evolves to be stable. When k � 0.1 and k � 0.5, their respective characteristic values are |λ 1 | � 0.9665, |λ 2 | � 0.5576, and |λ 3 | � 0.0095 and |λ 1 | � 0.5615, |λ 2 | � 0.482, and |λ 3 | � 0.0275. Figure 10 shows the behavior of the track stability of the control system (29) with other parameters unchanged. ese two graphs can find stable regions. As the feedback strength value increases, the chaotic behavior can be quickly controlled to stable orbit. is means that enterprise participants use time-delay feedback control to adjust strategies to quickly restore stable, orderly market competition.

Conclusion
In this paper, the dynamic properties of the three oligarchic games with product differentiation and heterogeneity are studied. Under the price competition, by studying the local stability of the equilibrium points, the following results were found.
First, under price competition, the behavior of gradient firms may constitute a source of instability. is phenomenon occurs when participants are more sensitive to market information and make their adjustments too fast. is result confirms the conclusions presented in the literature (see [14,15,21,43], etc.).
Second, the adaptive expectations of enterprises have an impact on the stability or instability of the game model, because if enterprises are more inclined to change the number of previous periods, Nash equilibrium will lose local stability.
is result is consistent with the conclusions in [14,43] in the absence of product differentiation.
ird, we discuss the effect of product differentiation on the local stability of the Nash equilibrium. We found that the product differentiation is larger; the Nash equilibrium will become unstable in complementary products reproducing the previous results proposed in the literature (see [15] for the duopoly case and [21] for a differentiated triopoly with linear dynamic system). However, the monotonic relationship between product differentiation and stability is different in alternative products. e threshold for adjusting the speed is the largest and has a certain degree of substitutability (for the fixed adaptive adjustment factor β); if the degree of difference between products is increased, the stability of the Nash equilibrium is difficult to guarantee (note d, the case where the right side approaches 0 is established). In the mutual substitution, this result is related to the complementary intertemporal strategy of the product. As a result, reducing the degree of product differentiation may undermine the Nash equilibrium.
Finally, the chaotic behavior of the dynamic system is controlled by time-delay feedback control and the numerical analysis shows that the unstable behavior of the dynamic system can be effectively and quickly controlled, which can quickly restore the stability and order of the market and provide reference for enterprise decision-makers to adjust strategies (in line with the conclusions of [21]).

Data Availability
e numerical simulations data used to support the findings of this study were supplied by Liuwei Zhao under license and so cannot be made freely available. Requests for access to these data should be made to Liuwei Zhao (e-mail address: 136901672@qq. com).

Conflicts of Interest
e author declares that there are no conflicts of interest.