ORIGINAL RESEARCH article
A Continuous Time Bertrand Duopoly Game With Fractional Delay and Conformable Derivative: Modeling, Discretization Process, Hopf Bifurcation, and Chaos
- 1College of Economics and Management, Shandong University of Science and Technology, Qingdao, China
- 2Department of Management, Polytechnic University of Marche, Ancona, Italy
The purpose of this paper is three-fold. First, we present a discretization process to obtain numerical solutions of a conformable fractional-order system with delays. Second, we extend the classical Bertrand duopoly game with integer delays to that with fractional delays. Third, we extend the game based on ordinary differential derivative to that based on conformable fractional-order derivative. Finally, we analyze the local stability, Hopf bifurcation, and chaos of the proposed game model.
Bifurcation and chaos  are frequent phenomena in various scientific fields including economics. By analyzing bifurcation and chaos, we can reveal its evolutionary mechanism to control proposed systems [2–4]. Delay is another common phenomenon which usually occurs in various economic systems, such as the discrete duopoly game with integer delays [5–13], the monopoly with integer delay and bounded rationality [14, 15], the continuous time Cournot duopoly with integer delays [16, 17], the continuous Bertrand duopoly with integer delay [18, 19]. Naturally, we will consider what the Bertrand duopoly game with fractional delays should be.
As a natural extension of a classical ordinary differential equation, the fractional derivative is a derivative of arbitrary order, real or complex [20–28]. Similarly we also want to extend the integer delay to arbitrary delay, integer or fractional, real or complex. Thus the continuous time Bertrand duopoly game with fractional delay and conformable derivative would have stronger ability to represent complex problems than the corresponding game with integer delay and derivative.In other words, the fractional delay can also be conceived as a kind of conformable fractional delay.
Although researchers proposed various fractional operators, such as Riemann-Liouville, Caputo, and Grunwald-Letnikov [29–33], widely used in many fields including economics, physics, and engineering [34–38], these definitions of operators satisfy two characteristics: one is that they must be non-integer in form, and the other is that they must be non-local in essence. But non-locality is a double-edged sword. It has the advantage of long memory, but it also has the disadvantage of not satisfying the classical differential operator, such as the chain rule. The conformable operator is conceived as a natural extension of the classical differential operator [39–45], whose most important properties hold, such as the chain rule [21, 22, 46]. Certainly, the conformable derivative is only a kind of local definition of fractional derivative and is non-integer in form, but it has no non-locality. In fact, there are two opposing viewpoints about the definitions of fractional derivatives: one is two above-mentioned characteristics must hold simultaneously, but the other is the first characteristic is the key of the fractional derivative. The later opinion means that a derivative is fractional derivative only if the first characteristic holds. We will support the later opinion in this paper. The conformable operator has recently occurred in many scientific fields [23–27, 40, 47–50]. In this paper, we will introduce the conformable fractional-order derivative to the continuous time Bertrand duopoly game with integer delays. Thus we can obtain a conformable Bertrand duopoly game with fractional delay.
After we propose the continuous time Bertrand duopoly game with fractional delay and conformable derivative, we need suitable approaches to obtain its solutions. Though there are several schemes to solve a conformable fractional-order system without any delay [40, 41, 45, 50–72], they are not suitable for conformable system with fractional delays. Mohammadnezhad et al.  employed conformable Euler's method to obtain approximation solutions of conformable fractional differential equations. Inspired by the discretization process for other derivatives [74–78], we will present a simple discretization process for a conformable system with fractional delays. Our proposed method well coincides with the conformable Euler method . Using the proposed discretization scheme, we will detect the stability, Hopf bifurcations, and chaotic attractors of the continuous time Bertrand duopoly game with fractional delay and conformable derivative.
The remainder of this paper is organized as follows. In section 2, we show preliminaries of conformable calculus and stability conditions of a discrete system. In section 3, we propose a discretization process for conformable systems with fractional delays. In section 4, we present a continuous time Bertrand duopoly game with fractional delays and conformable derivative. In section 5, we analyze Nash equilibrium points, local stabilities, and Hopf bifurcation in the game model. In section 6, we employ 0-1 test algorithm to detect chaotic attractors. This paper concludes with a summary in section 7.
Definition 1. (See )
For a function f : [t0, ∞) → ℝ, its left conformable fractional derivative starting from t0 of order α ∈ (0, 1) is defined by
in which the function f is called as α-differentiable.
Definition 2. (See )
For a function f : [t0, ∞) → ℝ, its left conformable fractional integral starting from t0 of order α ∈ (0, 1) is defined by
where the integral is the usual Riemann improper integral.
Lemma 1. (See )
Suppose the derivative order α ∈ (0, 1) and the function f is α-differentiable at a point t0 > 0, then the left conformable fractional derivative satisfies:
Lemma 2. (Conformable fractional power series expansion)  Given an infinitely α−differentiable function f : [0, ∞) → ℝ, α ∈ (0, 1], f has the following conformable fractional power series expansion at a neighborhood of a point t0:
where denotes the application of the conformable fractional derivative n times.
Consider the following discrete dynamical system:
where X = (x1, x2, ··· , xn), and F = (f1, f2, ··· , fn) is C2 from ℝn to ℝn. Let Xeq be a fixed point of system (4), and λ1, λ2, ··· , λn be eigenvalues of Jacobian matrix . Then we can obtain the following characteristic polynomial of J(Xeq)
Lemma 3. (See ).
(i) If all the eigenvalues λi of J(Xeq) have |λi| < 1, then the periodic orbit is attracting.
(ii) If one eigenvalue λi0 of J(Xeq) has |λi0| > 1, then the periodic orbit is unstable.
(iii) If all the eigenvalues λi of J(Xeq) have |λi| > 1, then the periodic orbit is repelling.
3. Discretization Process of Conformable Fractional Systems With Delays
Theorem 1. (Conformable delay discretization by piecewise constant approximation)
Consider the following conformable fractional-order delay system
we obtain the following discretization of Equation (6)
where xn denotes xn(tn), , tn = nh, n = 0, 1, ··· , N.
Proof: By drawing on the discretization process proposed by Raheem and Salman , El-Sayed et al. , El-Sayed et al. , Agarwal et al. , and El-Sayed et al. , we discretize the fractional-order delay Equation (6) with piecewise constant approximation as follows
By using Lemma 1, we rewrite Equation (8) as follows.
which leads to
(i) Let n = 0, then t ∈ [0, h), we rewrite Equation (9) as follows
and the solution of Equation (10) is given by
(ii) Let n = 1, then t ∈ [h, 2h), we rewrite Equation (9) as follows
and the solution of Equation (11) is given by
(iii) By repeating the above process, we obtain the following solution of Equation (9)
Let t → (n + 1)h, we deduce the the following discretization
It is proved.
Remark 1. The presented conformable discretization process is well line with the conformable Euler's method proposed by Mohammadnezhad et al. .
Theorem 2. The conformable delay discretization process in Theorem 1 is convergent.
Proof: According to Lemma 2 and the conformable Euler's method , we obtain the following conformable fractional power series expansion:
where 0 < θn < 1 and .
So the conformable delay discretization process is convergent when step size h is small enough.
4.1. A Continuous Time Bertrand Duopoly Game With Fractional Delays and Conformable Derivative
Assume that each firm decides on a different product pricing strategy in a duopoly market, in which pi(t) indicates the i-th firm's price during time period t ∈ Z+, qi(t) represents the i-th firm's supply during time period t ∈ Z+, and qi and pi satisfies the common linear inverse demand function as follows.
where ai and b are positive constants. ai denotes the i-th firm's constant reservation price in the product market. b denotes supply margin effects on its price from itself and its rival. Suppose i-th firm's marginal costs is equal to ci, which is also a positive constant, as follows
Then the i-th firm's profit is
Thus we can obtain the i−th firm's marginal profit with respect to pi as follows
As we know, the higher a firm's marginal profit is, the greater its price adjustment range is. In other words, the a firm's price growth rate is proportional to its marginal profit . We get the gradient adjustment process of price as follows
where vi > 0 represents the i−th firm's price adjustment speed.
Inspired by Matsumoto et al.  who regarded the marginal profit as a delayed value, we rewrite Equation (17) with delay τ as follows
We can regard the conformable fractional-order derivative as a natural extension of the integer order form. Thus, we introduce the conformable derivative to Equation (18) and obtain the following Bertrand duopoly game with αi-order conformable derivative and fractional delay:
where α ∈ (0, 1], di = ai − ci, τ ≥ 0 denotes a fractional delay.
Remark 2. When α = 1, Equation (19) degenerates to Equation (18).
Remark 3. When α = 1 and τ = 0, Equation (19) degenerates to Equation (17).
4.2. Discretization Process
According to Theorem 1, we employ piecewise constant approximation to discretize the game model (18) as follows
Remark 4. If τ = mh, and m is integer, then we rewrite the game model (20) as follows
In the following, we only study the case m = 1. Then we rewrite the game model (21) as follows
5. Local Stability, Hopf Bifurcation, and Chaos
5.1. Local Stability
Fixed points of system (22) satisfy
By algebraic computation, we obtain the following fixed points: E1 = (0, 0, 0, 0), , , and E4 = (pe1, pe2, pe1, pe2), where , . There are some economic information as follows:
(i) At the first fixed point E1, no firm can gain anything by production if its opponent stops production.
(ii) At the second fixed point E2, stoping production is the best strategy for firm 1 if is the output strategy of firm 2, and vice versa.
(iii) At the third fixed point E3, stoping production is the best strategy for firm 2 if firm 1 adopts its output strategy , and vice versa.
(iv) At the fourth fixed point E4, firm 1 cannot obtain extra benefit from deviating unilaterally from its equilibrium output strategy if firm 2 adopts its equilibrium output strategy , and vice versa.
Points E1, E2, and E3 are all bounded equilibria . Thus, we only study the stability of non-bounded equilibrium point E4.
We obtain the following Jacobian matrix of system (23) evaluated at the fourth fixed point E4:
The characteristic equation of Jacobian matrix (24) is
According to Jury stability criterion, we directly obtain the following proposition for system (22).
Proposition 1. The fixed point E4 is locally asymptotically stable if the following stability constraints holds simultaneously: W(1) > 0, W(−1) > 0, |w4| < 1, , |k1| < k2, where
Set h = 0.01, d1 = 8, d2 = 6, b = 0.01, v1 = 0.05, v2 = 0.01, α = 0.2. We can obtain E1 = (0, 0, 0, 0), E2 = (0, 3, 0, 3), E3 = (4, 0, 4, 0), and E4 = (3.99, 2.98, 3.99, 2.98). Obviously, Proposition 1 holds, i.e., system (22) is stable, as shown in Figures 1, 2. As time goes on with α = 0.2, Figure 1 demonstrates outputs of firm 1 trend to an invariant value, and Figure 2 demonstrates system (22) convergences to a fixed point. Figures 1, 2 mutually validate the existence of fixed points E4 in system (22).
5.2. Hopf Bifurcation
According to the explicit criterion of Hopf bifurcation , we directly obtain the following theorem:
Theorem 3. In system (22), a Hopf bifurcation occurs at α = α* if the following conditions (H1)-(H3) hold:
(H1) Eigenvalue assignment W(1) > 0, W(−1) > 0,
(H3) Nonresonance condition cos(2π/m) ≠ ψ or resonance condition cos(2π/m) = ψ, where m = 3, 4, 5, ··· and .
If h = 0.01, d1 = 8, d2 = 6, b = 0.01, v1 = 0.05, v2 = 0.01, a Hopf bifurcation will occur at α* = 0.1767. When we adopt a small perturbation Δα = 0.0007, a sufficiently small positive real number, i.e. α = α* + Δα = 0.176, system (22) has a stable closed invariant curve around the fixed point E4, as shown in Figures 3, 4. As time goes on with α = 0.1767, Figure 3 demonstrates outputs of firm 1 trend to periodic fluctuation, and Figure 4 demonstrates system (22) convergences to an invariant closed curve. Figures 3, 4 mutually validate the existence of Hopf bifurcation in system (22).
6. 0-1 Test for Chaos
When h = 0.01, d1 = 8, d2 = 6, b = 0.01, v1 = 0.05, v2 = 0.01, α = 0.2, system (22) is chaotic, as shown in Figures 5, 6. Figure 5 shows that the time series of p1 irregularly evolutes with time. Figure 6 shows that a chaotic attractor exists in system (22).
We obtain the 0-1 test value K for chaos as follows:
where ξ = (1, 2, 3, …, round(N/10)), Δ = (Dc(1), Dc(2), …, Dc(round(N/10))), the covariance is defined with vectors ξ, Δ of length m by
For system (22) with α = 0.137, we obtain K = 0.9977 and plot a phase portrait in (μ, ν)-plane whose trajectories in the (μ, ν)−plane are Brownian-like, which indicate it is chaotic, as shown Figure 7.
As time goes on with α = 0.137, Figure 5 demonstrates outputs of firm 1 trend to irregular vibration, and Figure 6 demonstrates system (22) has chaotic attractors in (p1, p2, p4)-plane, and Figure 7 demonstrates the motion trajectories of system (22) are Brownian-like in (μ, ν)-plane. Figures 5–7 mutually validate the existence of chaos in system (22).
Introducing conformable derivative to a classical continuous Bertrand duopoly game with integer delays, we present a generalized game model with delays, named a continuous time Bertrand duopoly game with fractional delay and conformable derivative. We design a discretization process for conformable systems with fractional delays. Utilizing the proposed discretization process, we obtain a simplified model and get its solutions. At last, we examine the local stability, Hopf bifurcation, and chaos of the game model.
There are some future research directions that deserve our attention. We can apply the modeling approach and discretization process to other scientific researches, such as conformable partial differential equations with delays [90–94], conformable differential games with delays . In addition, we can also develop other discretization methods for conformable system with specific characteristics, such as stochastic, fuzzy, pulse factors.
All datasets generated for this study are included in the manuscript and/or the supplementary files.
BX and LG: conceptualization. BX, WP, and LG: methodology. BX and WP: software. BX and LG: validation. BX: formal analysis. BX and WP: writing-original draft preparation. BX and LG: writing-review and editing. BX and WP: visualization. BX: supervision, project administration, and funding acquisition.
This work is supported partly by Natural Science Foundation of Shandong Province (Grant No. ZR2016FM26), National Social Science Foundation of China (Grant No.16FJY008), and National Natural Science Foundation of China (Grant No.11801060).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The reviewer ZW declared a shared affiliation, with no collaboration, with the authors BX, WP to the handling editor at time of review.
The authors would like to express many thanks to referees for their time and efforts in reviewing this manuscript. Their helpful comments and constructive suggestions greatly improved this manuscript.
1. Wei Z, Zhu B, Yang J, Perc M, Slavinec M. Bifurcation analysis of two disc dynamos with viscous friction and multiple time delays. Appl Math Comput. (2019) 347:265–81. doi: 10.1016/j.amc.2018.10.090
4. El-Sayed A, Nour H, Elsaid A, Matouk A, Elsonbaty A. Dynamical behaviors, circuit realization, chaos control, and synchronization of a new fractional order hyperchaotic system. Appl Math Model. (2016) 40:3516–34. doi: 10.1016/j.apm.2015.10.010
10. Pecora N, Sodini M. A heterogenous Cournot duopoly with delay dynamics: Hopf bifurcations and stability switching curves. Commun Nonlinear Sci Numer Simul. (2018) 58:36–46. doi: 10.1016/j.cnsns.2017.06.015
11. Bignami F, Agliari A. Chaotic dynamics in a three-dimensional map with separate third iterate: the case of Cournot duopoly with delayed expectations. Chaos Soliton Fract. (2018) 110:216–225. doi: 10.1016/j.chaos.2018.03.023
12. Askar SS, Alnowibet K. Cooperation versus noncooperation: Cournot duopolistic game based on delay and time-dependent parameters. Chaos Soliton Fract. (2016) 91:580–4. doi: 10.1016/j.chaos.2016.08.004
20. Xin B, Peng W, Kwon Y. Modeling, discretization, and hyperchaos detection of conformable derivative approach to a financial system with market confidence and ethics risk. Adv Differ Equat. (2019) 2019:138. doi: 10.1186/s13662-019-2074-8
23. Acan O, Firat O, Keskin Y. Conformable variational iteration method, conformable fractional reduced differential transform method and conformable homotopy analysis method for non-linear fractional partial differential equations. Wave Rand Complex Media. (2018). doi: 10.1080/17455030.2018.1502485. [Epub ahead of print].
31. Wang F, Chen D, Zhang X, Wu Y. Finite-time stability of a class of nonlinear fractional-order system with the discrete time delay. Int J Syst Sci. (2017) 48:984–93. doi: 10.1080/00207721.2016.1226985
36. Kumar D, Agarwal RP, Singh J. A modified numerical scheme and convergence analysis for fractional model of Lienard's equation. J Comput Appl Math. (2018) 339:405–13. doi: 10.1016/j.cam.2017.03.011
38. Hajipour M, Jajarmi A, Baleanu D, Sun H. On an accurate discretization of a variable-order fractional reaction-diffusion equation. Commun Nonlinear Sci Numer Simul. (2019) 69:119–33. doi: 10.1016/j.cnsns.2018.09.004
39. Morales-Delgado V, Gómez-Aguilar J, Escobar-Jiménez R, Taneco-Hernández M. Fractional conformable derivatives of Liouville–Caputo type with low-fractionality. Phys A. (2018) 503:424–38. doi: 10.1016/j.physa.2018.03.018
40. Pérez J, Gómez-Aguilar J, Baleanu D, Tchier F. Chaotic attractors with fractional conformable derivatives in the Liouville–Caputo sense and its dynamical behaviors. Entropy. (2018) 20:384. doi: 10.3390/e20050384
41. Yépez-Martínez H, Gómez-Aguilar J, Atangana A. First integral method for non-linear differential equations with conformable derivative. Math Model Nat Phenom. (2018) 13:14. doi: 10.1051/mmnp/2018012
42. Eslami M, Khodadad FS, Nazari F, Rezazadeh H. The first integral method applied to the Bogoyavlenskii equations by means of conformable fractional derivative. Opt Quant Electron. (2017) 49:391. doi: 10.1007/s11082-017-1224-z
43. Rezazadeh H, Korkmaz A, Eslami M, Vahidi J, Asghari R. Traveling wave solution of conformable fractional generalized reaction Duffing model by generalized projective Riccati equation method. Opt Quant Electron. (2018) 50:150. doi: 10.1007/s11082-018-1416-1
45. Khodadad FS, Nazari F, Eslami M, Rezazadeh H. Soliton solutions of the conformable fractional Zakharov–Kuznetsov equation with dual-power law nonlinearity. Opt Quant Electron. (2017) 49:384. doi: 10.1007/s11082-017-1225-y
46. Abdeljawad T, Al-Mdallal QM, Jarad F. Fractional logistic models in the frame of fractional operators generated by conformable derivatives. Chaos Soliton Fract. (2019) 119:94–101. doi: 10.1016/j.chaos.2018.12.015
48. Rezazadeh H, Khodadad FS, Manafian J. New structure for exact solutions of nonlinear time fractional Sharma-Tasso-Olver equation via conformable fractional derivative. Appl Appl Math. (2017) 12:405–14.
53. Hosseini K, Bekir A, Ansari R. New exact solutions of the conformable time-fractional Cahn–Allen and Cahn–Hilliard equations using the modified Kudryashov method. Optik. (2017) 132:203–9. doi: 10.1016/j.ijleo.2016.12.032
55. Kumar D, Seadawy AR, Joardar AK. Modified Kudryashov method via new exact solutions for some conformable fractional differential equations arising in mathematical biology. Chin J Phys. (2018) 56:75–85. doi: 10.1016/j.cjph.2017.11.020
56. Srivastava HM, Günerhan H. Analytical and approximate solutions of fractional-order susceptible-infected-recovered epidemic model of childhood disease. Math Meth Appl Sci. (2019) 42:935–41. doi: 10.1002/mma.5396
59. Kartal S, Gurcan F. Discretization of conformable fractional differential equations by a piecewise constant approximation. Int J Comput Math. (2018) doi: 10.1080/00207160.2018.1536782. [Epub ahead of print].
60. Iyiola OS, Tasbozan O, Kurt A, Çenesiz Y. On the analytical solutions of the system of conformable time-fractional Robertson equations with 1-D diffusion. Chaos Soliton Fract. (2017) 94:1–7. doi: 10.1016/j.chaos.2016.11.003
62. He S, Sun K, Mei X, Yan B, Xu S. Numerical analysis of a fractional-order chaotic system based on conformable fractional-order derivative. Euro Phys J Plus. (2017) 132:36. doi: 10.1140/epjp/i2017-11306-3
63. Yokuş A. Comparison of Caputo and conformable derivatives for time-fractional Korteweg-de Vries equation via the finite difference method. Int J Modern Phys B. (2018) 32:1850365. doi: 10.1142/S0217979218503654
70. Liu S, Wang H, Li X, Li H. The extremal iteration solution to a coupled system of nonlinear conformable fractional differential equations. J Nonlinear Sci Appl. (2017) 10:5082–9. doi: 10.22436/jnsa.010.09.44
72. Morales-Delgado V, Gómez-Aguilar J, Taneco-Hernandez M. Analytical solutions of electrical circuits described by fractional conformable derivatives in Liouville-Caputo sense. AEU-Int J Elect Commun. (2018) 85:108–17. doi: 10.1016/j.aeue.2017.12.031
73. Mohammadnezhad V, Eslami M, Rezazadeh H. Stability analysis of linear conformable fractional differential equations system with time delays. Bol Soc Paran Math. (2020) 38:159–71. doi: 10.5269/bspm.v38i6.37010
88. Munoz-Pacheco J, Zambrano-Serrano E, Volos C, Jafari S, Kengne J, Rajagopal K. A new fractional-order chaotic system with different families of hidden and self-excited attractors. Entropy. (2018) 20:564. doi: 10.3390/e20080564
89. Bao BC, Wu PY, Bao H, Xu Q, Chen M. Numerical and experimental confirmations of quasi-periodic behavior and chaotic bursting in third-order autonomous memristive oscillator. Chaos Soliton Fract. (2018) 106:161–70. doi: 10.1016/j.chaos.2017.11.025
90. Yépez-Martínez H, Gómez-Aguilar J. Fractional sub-equation method for Hirota–Satsuma-coupled equation and coupled mKdV equation using the Atangana's conformable derivative. Wave Rand Complex Media. (2018). doi: 10.1080/17455030.2018.1464233. [Epub ahead of print].
91. Yépez-Martínez H, Gómez-Aguilar J, Baleanu D. Beta-derivative and sub-equation method applied to the optical solitons in medium with parabolic law nonlinearity and higher order dispersion. Optik. (2018) 155:357–65. doi: 10.1016/j.ijleo.2017.10.104
94. Eslami M, Rezazadeh H, Rezazadeh M, Mosavi SS. Exact solutions to the space–time fractional Schrödinger–Hirota equation and the space–time modified KdV–Zakharov–Kuznetsov equation. Opt Quant Electron. (2017) 49:279. doi: 10.1007/s11082-017-1112-6
Keywords: conformable calculus, fractional-order Bertrand game, fractional delay, Hopf bifurcation, 0-1 test for chaos
Citation: Xin B, Peng W and Guerrini L (2019) A Continuous Time Bertrand Duopoly Game With Fractional Delay and Conformable Derivative: Modeling, Discretization Process, Hopf Bifurcation, and Chaos. Front. Phys. 7:84. doi: 10.3389/fphy.2019.00084
Received: 11 April 2019; Accepted: 16 May 2019;
Published: 04 June 2019.
Edited by:Carla M. A. Pinto, Instituto Superior de Engenharia do Porto (ISEP), Portugal
Reviewed by:Devendra Kumar, University of Rajasthan, India
Zhen Wang, Shandong University of Science and Technology, China
Amin Jajarmi, University of Bojnord, Iran
Mostafa Eslami, University of Mazandaran, Iran
Francisco Gomez, Centro Nacional de Investigación y Desarrollo Tecnológico, Mexico
Copyright © 2019 Xin, Peng and Guerrini. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Baogui Xin, firstname.lastname@example.org