Abstract

A discreet fractional-order Cournot–Bertrand competition duopoly game is introduced based on the fractional-order difference calculus of the Caputo operator. The model is designed when players can make long memory decisions. The local stability of equilibrium points is discussed for the proposed model. Some numerical simulations explore the model’s bifurcation and chaos by employing bifurcation diagrams, phase portraits, maximal Lyapunov exponents, and time series. According to our findings, the fractional-order parameter has an effect on the game’s stability and dynamics.

1. Introduction

Game theory is one of the most interesting and complex topics that many researchers are interested in understanding. Game theory is concerned with predicting results for strategic games in which participants, for example, two or more firms competing on the market, have incomplete information on the intentions of others. It is known that the game theory is relevant to the study of corporate behavior in oligopolistic markets, for example, the decisions that companies must make in terms of production and pricing levels, as well as the amount of money invested in research and development. The decision-making mechanism has an important role to play in the process of adjusting the production and profits of firms. Firms typically use one of the following to change their market growth: naïve learning expectation, adaptive learning expectation, limited learning rationality, and local learning approximation. Discrete oligopoly dynamics based on company profit maximizations have been considered [17]. Furthermore, these models have been utilized to examine the dynamic characteristics of competitive markets, which has been classified as steady state, periodic, and chaotic [814].

Fractional calculus, particularly discrete fractional calculus, has attracted substantial interest in recent decades due to its extensive significance in a wide range of scientific disciplines, including complex systems with memory and heredity. Researchers turned their attention to a discreet fractional calculus and tried to develop a complete theoretical framework for this subject. This is due to the importance of this field in many real issues, such as discreet adaptive systems, biological growth systems, and digital engineering systems, all of which contain memory [1520]. The discrete difference analogues of classical Caputo and Riemann derivatives have been introduced in [21]. In addition, advances have been made in the study of fractional finite difference equations and the inclusion of fractional difference equations [2226]. Recently, the stability of fractional time systems in a variety of real-world problems has been investigated. These studies have shown that discrete fractional systems are more realistic to process real systems and have a rich dynamic compared to discrete systems with integer-order. Many studies have studied continuous fractional differentiation representing the effects of economic memory that have been presented [2730] and the references that exist in them as well. The fractional-order difference equation, which is a natural extension of the integer-order difference equation, has long-term memory effects that have been explored in a few studies [25, 3136]. Recently, the fractional difference calculus was used in the Cournot duopoly game [37] and the Bertrand duopoly game [38]. This reflects the long-term memory effects of Cournot–Bertrand dynamic games in the fractional-order form for the game. Xin et al. [37, 38] investigated the dynamics of the Cournot and Bertrand games, which indicated the market evolution of the two enterprises.

This work is especially interested in the novel discrete fractional-order Cournot–Bertrand duopoly game, which is a modification of the Cournot–Bertrand duopoly game with integer-order [39]. We aim to extend this game to a fractional case and to study the dynamics of the game. As we know, fractional-order calculus is a general form of integer-order calculus, so it has a higher representation for phenomena with a long memory. It can be shown that the fractional differentiation parameter is a vital indicator of the bifurcation path and the chaos that is created and disappeared. We will investigate the dynamics of the discreet fractional-order Cournot–Bertrand duopoly game such as the stability, bifurcation, and chaos of the proposed game. To analyze complexity of the game, explicit stability criteria [40, 41], asymptotic stability criteria [42] and the local stability regions of the boundary and Nash equilibrium points are provided through numerical simulation. The dynamic behavior of the proposed game is illustrated through an exploratory investigation of equilibrium point stability and numerical simulation. In this work, we are analyzing a Cournot–Bertrand duopoly game similar to Wang and Ma [39], but using a discrete fractional calculus. The equilibrium point structure dynamic reflects economic explanations for the proposed game’s market of two enterprises.

The work is organized as seen follows. Section 2 describes the market dynamics of two enterprises using a discrete fractional-order Cournot–Bertrand duopoly game. The Nash equilibrium local stability conditions are established in Section 3. Using numerical simulations, we investigate local bifurcations, maximal Lyapunov exponents, and phase portraits of complex dynamics in Section 4. Section 5 contains a summary of the findings as well as a few observations.

2. Preliminaries

In this section, some preliminaries about fractional-order difference calculus are introduced. On an arbitrary time scale, dynamic behaviors and applications of fractional difference models were investigated in the last decade where delta difference equation was used [4043].

Assume that a sequence is given, and the isolated time scale is represented in terms of real valued constant , i.e., , such that . Also, the difference operator is denoted by , where . Then, we summarize some of the basic definitions related to discrete fractional calculus as follows.

Definition 1. For , the fractional sum of order is given by [21]

Definition 2. The Caputo-like delta difference of order is defined by [21, 42]

Theorem 1. (See [21, 4043]). In order to solve the delta fractional difference equation,

As a result, the corresponding discrete integral equation iswhere

Remark 1. If , we rewrite discrete integral equation in the following numerical form:

Theorem 2 (See [4042]). The linear discrete-time fractional-order system is where , , , and . The zero equilibrium of system (10) is asymptotically stable iffor all the eigenvalues of matrix .

3. Discrete Fractional-Order Cournot–Bertrand Duopoly Game

According to traditional oligopoly models [1, 6], firms compete in the same strategic variable, such as output (Cournot) or price (Bertrand). A hybrid model, commonly known as the Cournot–Bertrand model [6, 39], permits certain enterprises to compete in output, while others compete in pricing. Real-world market observations that match Cournot–Bertrand behavior have bolstered the model’s validity and rapidly growing literature on advancements and applications. Long-term memory effects in dynamic oligopoly games are economically significant [34, 37, 38]. Therefore, we introduce the novel discrete fractional-order Cournot–Bertrand duopoly game, which is a modification of the Cournot–Bertrand duopoly game with integer order. As a consequence, it has a better representation of phenomena with a long memory in oligopoly game. The main goal is to establish out how the fractal parameter affects game dynamics including stability, bifurcation, and chaos.

We suggest a simple Cournot–Bertrand duopoly common oligopoly game [39]. Two enterprises, denoted by the letters , produced differentiated goods with perfect replacements and set their product pricing based on the same market rule. Assume that and denote the goods price and quantity output of firm for the period . The inverse demand functions for a variety of products 1 and 2 originate from the representative consumer maximization of the following utility function [39]:subjected to restrictions on the budget . Then, the inverse demand functions is defined as follows:where and represent the pricing of firm 1’s and firm 2’s items, respectively, and and are the quantities of products of company 1 and company 2. The parameter is the product differentiation between two firms. Products are homogeneous goods when , and each firm has a monopoly case when . The demand system can be written in two strategic variables and .

Consider that the two companies have the same unit cost and that the cost function has the same linear form:

Thus, the profit functions for firms are given by

In the classical dynamical Cournot–Bertrand duopoly game, to decide the corresponding profit-maximizing, every player erroneously believes that its rival’s pricing in period is the same as in period . Therefore, this type of game does not have a long memory effect. The traditional game will be introduced using discrete fractional-order calculus, and the two players will make decisions using a new dynamic adjustment mechanism with long memory and local marginal profit expectation. Thus, the marginal profit of two players is as follows [39]:

Assume that the two businesses have limited information about the market demand function and also their price decision is based on a dynamic adjustment process with limited rationality and a long-term memory effect of marginal profit. In the next step, the firm decides to raise (reduce) the price of its product based on if the long-term marginal profit is greater (less) than zero. As a result, we utilize the dynamic adjustment process for the Cournot–Bertrand game as follows:where is the speed of adjustment of firm and denotes a fractional-order number, indicating the long-term memory effect. Thus, the discrete fractional-order Cournot–Bertrand duopoly game is as follows:

Remark 2. When , the model (16) devolves to the Wang-Ma model [39]:

We will show that the model game parameters, especially the long-term memory effect, have an effect on the long-term dynamic response of our novel game when compared to the Wang-Ma game [39].

In the next sections, several theoretical features corresponding to game (16) are investigated.

4. The Equilibrium Points and Their Local Stability

We solve the following equation to find the equilibrium points of game system (16):

Their four equilibria are , , , and . In economics, its equilibria mean(i)The equilibrium is trivial equilibrium point(ii)The equilibrium point implies that the best quantity of the first player is if the second player sets its optimal product price . Likewise, the second player’s best price is if the player uses a zero-price approach. Clearly, is a border equilibrium point that corresponds to the pure monopoly.(iii)The equilibria indicates that the best quantity of the first player is if the second company determines the best good price . Likewise, the first company’s best pricing is if the company uses a zero-price approach. Clearly, is a border equilibrium point.(iv)The equilibrium indicates two enterprises will preserve their equilibrium quantity and pricing jointly since no enterprise can gain an advantage by deviating unilaterally from its own equilibrium. Clearly, the point is a Nash equilibrium. The complexity of system (16) will be explored. First, the Jacobian matrix of system (16) is computed at a given equilibrium point , and it can be expressed as

Then, the following theorems are presented to clarify linear stability for each equilibrium point in the model.

The trivial equilibrium point is unstable.

Proof. The Jacobian matrix’s eigenvalues at can be demonstrated to be , which indicate that one of them is always positive, and thus, the conditions of asymptotic stability in Theorem 2 are not satisfied.

The equilibrium point is asymptotically stable if

Proof. The eigenvalues of Jacobian matrix at are , which means that the conditions of asymptotic stability in Theorem 2 are satisfied if the second eigenvalue has a negative sign, and also, the modulus of the two eigenvalues is bounded by and , respectively. Figure 1 shows stability regions in some three and two-dimensional spaces of model’s parameters, whereas Figure 2 shows the resulting time series when the values of parameters are selected in stable regions of .
The equilibrium is asymptotically stable if and only if

Proof. The eigenvalues of Jacobian matrix at are , which means that the conditions of asymptotic stability in Theorem 2 are satisfied if the two eigenvalues have negative signs, and also, the modulus of the two eigenvalues is bounded by and , respectively.

However, detailed numerical examinations in space of promoters show that the aforementioned conditions cannot be simultaneously achieved at feasible values of parameters, and therefore, due to the impossibility of numerically satisfying the above conditions, the equilibrium point is unstable.

Finally, the Nash equilibrium point has long complicated expressions for its associated eigenvalues which renders numerical investigations of regions of stability inevitable. Figure 3 shows stability regions in some three and two-dimensional spaces of model’s parameters, whereas Figure 4 illustrates the resulting time series when the values of parameters are selected in stable regions of .

5. Numerical Simulations

In this section, the complex dynamic features of the discrete fractional Cournot–Bertrand model (16) are investigated using various methods such as bifurcation diagrams, phase portraits, and MLE. The effects of major model parameters are investigated. For the present fractional discrete model (16) using theorem (4), the system (16) can be numerically rewritten as follows:

The initial point is used in the following simulations. In particular, the complicated dynamics exhibited by the model are examined via using the bifurcation diagram, phase portraits, and maximal Lyapunov exponent (MLE). The Lyapunov exponent for a one-dimensional map can be calculated by calculating the average value for perturbations from the trajectory over a time interval. The Lyapunov exponents for an n-dimensional mapping can be obtained using the eigenvalues of the product of Jacobian matrices for integer-order systems. In order to approximate the values of Lyapunov exponents of the discrete fractional model (16), the Jacobian matrix method which has been introduced by Wu and Baleanu [44] can be employed [45, 46]. In the following part, the numerical analysis will look at the effects of the model’s main parameters, as well as the effects of long-term memory and adjustment speeds.

First, the influences of parameter in integer-order and fractional-order cases are explored. Figure 5 shows bifurcation diagrams and MLE plots versus parameter along with examples of phase portraits at some selected values of parameters.

Second, the influences of parameter in integer-order and fractional-order cases are explored. Figure 6 shows bifurcation diagrams and MLE plots versus parameter , along with examples of phase portraits at some selected values of parameters.

Third, the influences of parameter in integer-order and fractional-order cases are explored. Figure 7 shows bifurcation diagrams and MLE plots versus parameter along with examples of phase portraits at some selected values of parameters.

Finally, the effects of fractional-order are explored. Figure 8 shows bifurcation diagrams and MLE plots versus parameter along with examples of phase portraits at some selected values of parameters.

6. Conclusion

The new fractional-order Cournot–Bertrand game based on a long memory effect is proposed. The stability of the game’s equilibrium points, including the Nash equilibria, has been explored both qualitatively and numerically. Bifurcation, phase portraiture, time series, and maximal Lyapunov exponents’ diagrams have been used to analyze the complex dynamic characteristics of the proposed game. When we compared our new model to the Wang-Ma model [39], we found that the game parameters, especially the long-term memory influence, had an effect on the long-term dynamic response of our novel model. This is important for understanding the performance of the duopoly game with the long-term memory effect. According to our findings, the Cournot–Bertrand duopoly game with the long-term memory effect is more efficient than the duopoly game without long-term memory impact from economic viewpoint. As a consequence, we advise researchers to investigate the competitive games of long-term memory impact further.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors would like to express sincere gratitude to all referees and editors for their valuable comments and constructive suggestion to the earlier version of this study. They would also like to thank the Deanship of Scientific Research at Prince Sattam bin Abdulaziz University, Al-Kharj, Saudi Arabia.