Advanced Nonlinear Dynamics of Population Biology and EpidemiologyView this Special Issue
Hopf Bifurcation Analysis in a Modified Price Differential Equation Model with Two Delays
The paper investigates the behavior of price differential equation model based on economic theory with two delays. The primary aim of this thesis is to provide a research method to explore the undeveloped areas of the price model with two delays. Firstly, we modify the traditional price model by considering demand function as a downward opening quadratic function, and supply and demand functions both depending on the price of the past and the present. Then the price model with two delays is established. Secondly, by considering the price model with one delay, we get the stable interval. Regarding another delay as a parameter, we studied the linear stability and local Hopf bifurcation. In addition, we pay attention to the direction and stability of the bifurcating periodic solutions which are derived by using the normal form theory and center manifold method. Afterwards, the study turns to simulate the results through numerical analysis, which shows that the provided method is valid.
Recent developments in mathematical economics and in problems of business administration have led to extensive use of differential equation model. Bifurcations and chaos always show in the contemporary literature of economics as basic concepts. The paper of Shuhe  was the pioneering work in studying price differential equation model which provides a dynamic system to investigate sectoral dynamics of an economy phenomenon. Since demand depends on the past price, further study on price model with delay can be found in . Reference  focused on the phenomenon of bifurcation which forms an integral part of qualitative approach to study dynamical systems. Similarly, [4–6], which investigated the local Hopf bifurcation and the existence of periodic solutions of price model, had important consequences for theoretical and empirical model building in economics. Reference  provided a brief survey of the literature of bifurcation detections in economic models. The existence of different bifurcate parameters leads to chaos and the causes of complicated phenomenon were argued in . In short, few people studied price model with delay; what is more, no results involving price model with two delays have occurred.
In order to illustrate the economic phenomena with price varying accurately, a reasonable mathematical model of price is needed. Thus, we introduce a traditional price differential equation model in : The meaning of parameters refers to . We modify the traditional price differential equation model by considering the following factors.
Firstly, we denote the correlation coefficient between demand and price rising rate by . We consider demand function as a downward opening quadratic function and supply function as a linear function:
Secondly, according to the cobweb theory, since manufacturers need a production cycle time from obtaining market information to adjust the production line, the role of price adjustment lags on the supply function of time. The purchase of the consumers also depends on the price change and decision whether to buy or not. Then, supply and demand functions both depend on the price of the past and the present. We introduce two delays , to denote supply and demand functions, which depend on the price of the past, respectively. A differential system with two delays for price differential equation model is transformed into the following form:
Let , , and notice that supply and demand functions both depend on the price of the past and the present; then the model is described by the following autonomous system: where , refers to . is the price at time , and is the amount of supply at time ; , , , , , .
Different from the previous work in [1, 2, 4], the purpose of this paper is to investigate the stability of the Hopf bifurcation and the direction of bifurcation periodic solution of a price differential equation model with two delays. The structure of the paper is as follows. In Section 2, linear stability and local Hopf bifurcations are studied by using qualitative methods. In Section 3, we regard as bifurcation parameter and consider (4) with in its stable interval; the direction of Hopf bifurcation and the stability of the bifurcation periodic solutions are derived by using the normal form theory and center manifold method. Afterwards, the presented numerical simulations in Section 4 have demonstrated the theoretical analysis.
2. Local Stability and Hopf Bifurcation
Obviously, system (4) always has an equilibrium . In the following, we will investigate the effect of the delay , on the dynamics of system (4).
Let , ; the linearization of system (4) at zero steady state is
For convenience, as in literature of , we denote , , . Then (5) become
The corresponding characteristic equation is
To study the stability of equilibrium of (4) and Hopf bifurcation, it is sufficient to analyze the distribution of the roots of (7). It is stable if all roots of (7) have negative real parts and unstable if one root has positive real part. In order to study the characteristic (7) with two delays, we first consider (7) with one delay. Without loss of generality, we choose as a parameter and employ Rouche's theorem in Cooke and Grossman ; we will find the stable interval for . Then we consider (7) with in its stable intervals. Using Rouche's theorem again, regard as a parameter; we will find the stable interval (depends on ) for . Then we obtain the stable interval for system (4).
Now we analyse the case when ; (7) becomes Using a procedure similar to , we make some hypotheses as follows: (H1),(H2),(H3),(H4), , and ,(H5)neither (H3) nor (H4).Obviously, when , (4) becomes a system of ODE. Under the hypothesis (H1), all roots of (7) have negative real parts if and only if ; under the hypothesis (H2), all roots of (8) have negative real parts if and only if and . Above all, either (H1) or (H2) holds for , and all roots of (7) have negative real parts.
Applying the lemma in  again, we obtain the following results.
Lemma 1. For (8), one has the following: (i)if (H3) holds and , then (8) has a pair of purely imaginary roots ;(ii)if (H4) holds and then (8) has a pair of imaginary roots of ;(iii)if (H5) holds and , then (8) has no purely imaginary root, where
Proof. Let be the root of (8); we thus have Separating the real and imaginary parts Eliminating the harmonic terms gives Obviously From (12) and by direct computation, we obtain , , . The result is straightforward.
Denote the minimum value of by ; that is, , and
To see if and are bifurcation values, we need to verify if the transversality conditions hold.
Lemma 2. The following transversality conditions hold.
Proof. Differentiating (8) with respect to yields then we have which satisfied and ; then for , , and for , . We complete the proof.
Thus, we get the distribution of the characteristic roots of (8).
Lemma 3. For (8), one has the following: (i)if (H3) and either (1) (H1) or (2) (H2) and hold, then when , all roots of (8) have negative real parts, and when , (8) has at least one root with positive real part;(ii)if (H4) and either (H1) or (H2) hold, then there are k switches from stability to instability; that is, when , , all roots of (8) have negative real parts, where , and , , and (8) has at least one root with positive real part.That is to say, under those conditions when , , system (4) is asymptotically stable, and system (4) undergoes a Hopf bifurcation when .
Then, we consider stable interval for in which all roots of (8) have negative real parts, regarding as a parameter.
Lemma 4. If all roots of (8) have negative real parts, then there exists a , such that when all roots of (7) have negative real parts.
Proof. All roots of (8) have negative real parts which means that system (4) is stable when . In what follows, we consider (4) with fixed in its stable interval, regarding as a parameter. Let be a root of (4); then we obtain
Suppose that . Since and , then (19) has at least one positive root. Without loss of generality, the roots of (19) are defined by . For every , there exists a sequence , such that (19) holds. The expression of and can be derived by (19) for fixed ; we will calculate them directly by the use of Mathematica software in Section 4; here we omit them. Let and let be the positive and simple root of (19) when . When , (19) has a pair of purely imaginary roots for . Then, as varies, the sum of the multiplicities of zeros in the open right half-plane can change only if a zero appears on or crosses the imaginary axis. In what follows, we assume that (H6) . Therefore, by the general Hopf bifurcation theorem for FDEs in Wei and Ruan , we get that when all roots of (7) have negative real parts.
Applying the discussion above and noticing that all roots of (8) have negative real parts, we know that there exist such that all roots of (7) with have negative real parts. The proof is complete.
Summarizing the above lemmas and literature of Hale , we obtain the following sufficient conditions for all characteristic roots of (8) to have negative real parts.
Theorem 5. Suppose that (H6) holds and either (H1) or (H2) is satisfied. (i)If (H3) holds, then for any there exists a such that when , all roots of (7) have negative real parts.(ii)If (H4) holds, then for any there exists a , such that when all roots of (7) have negative real parts, where and are defined by (10).(iii)If (H5) holds, then for any , there exists a , such that when all roots of (7) have negative real parts.That is to say, under the conditions that is stable interval, there exist a (depend on ) such that when , system (4) is asymptotically stable, and system (4) undergoes a Hopf bifurcation when , ; .
3. Direction and Stability of the Bifurcating Periodic Solutions
In the previous section, we obtain the conditions under which a family of periodic solutions bifurcate from the steady state and the equilibrium loses its stability when , ; for fixed , and the relationship between and can be derived by (19). Throughout this section, by using techniques of the normal form and center manifold theory due to Hale , we derive the algorithm for determining the direction of the Hopf bifurcations and the stability of the bifurcating periodic solutions at critical values on the center manifold.
Since the analysis is local, we regard , as bifurcation parameter. Choosing the space as and for , system (3) is transformed into FDE as with where , , and . Obviously, is continuous linear function mapping into . By the Riesz representation theorem, there exists a matrix whose elements are bounded variation functions in such that In fact, we choose where is a delta function.
For the operators and are defined as Let ; the adjoint operator corresponding to is defined as follows: Then system (20) can be written in the following form: , where for .
For and , define the adjoint bilinear: where .
Proposition 6. Let and be eigenvectors of and corresponding to and , respectively, satisfying and . Then where
Proof. We assume that is the eigenvector of corresponding to and is the eigenvector of corresponding to . It follows from the definition of , , (22), and (23) that we have and ; we have where is identity matrix; that is, By direct computation and considering , , we obtain and . Now, we calculate as follows: Since , we get Therefore, . This completes the proof.
Then, we construct the coordinates of the center manifold at . Let On the center manifold , we have where and and are local coordinates for the center manifold in the direction of and , respectively. Since , we have where We rewrite in abbreviated form as where
By (20) and (38), we obtain where Substituting (26) and (38) into on the center manifold and comparing the coefficients we get Comparing the coefficients with (41) gives that From (42), (43), and the definition of , we can derive the following equation: Solving for and , we get where , are two constant vectors, which can be determined by setting in .
By (33), , and noticing we have Then it is easy to obtain where Thus, from (39), it follows that Comparing the coefficient with (39), we have From (40), we get that ; that is, By the definition of and (42) we have Notice that Substituting (45) and (52) into (51), we obtain By direct computation, we obtain
Thus, we can compute the following values which determine the properties of bifurcating periodic solutions at the critical value : More specifically (see Hassard et al. ), determines the direction of the Hopf bifurcation: if (), then the Hopf bifurcation is forward (backward) and the bifurcating periodic solutions exist for (). determines the stability of the bifurcating periodic solutions: the bifurcating periodic solutions are stable (unstable) if (). determines the period of the bifurcating periodic solutions: the period increases (decreases) if ().
4. Numerical Simulation
In order to validate the theoretical analysis, we will present some numerical simulations. We take the following coefficients as an example: , , , , , , , , and , it is easy to obtain , , and . We choose ; then from (19) and by means of Mathematica software, we get and . Taking as the initial conditions. By Theorem 5 and the above results, we know equilibrium is locally asymptotically stable when as is illustrated in Figure 1. When , is unstable and periodic solutions occur from ; we take , the corresponding phase plots are shown in Figure 2.
Finally, the numerical simulation shows that it is a complex transformation process for the system changes from stable equilibrium to chaos.
Different from the previous work in [1, 2, 4], the main contribution of this paper lies in the following aspects. Firstly, we modify the traditional price differential equation model by considering demand function that is settled as a downward opening quadratic function and considering supply and demand functions that are both depending on the price of the past and the present. Then the price differential equation model with two delays is established. Secondly, to study the stability and Hopf bifurcation of system (4), we consider (4) with one delay and find the stable interval for . In the following, regarding as a parameter, we obtain the stable interval for . Then we get the stable interval for (4). In addition, we derive the algorithm for determining the direction of the Hopf bifurcation and the stability of the bifurcating periodic solutions by using techniques of normal form theory and center manifold method. Lastly, a numerical analysis confirms the effectiveness of our research results. The paper provided the preparative work for further discussion. For instance, we consider that supply cannot increase with the price without limit; the production capacity of enterprises and social resources are limited; then we can modify the supply function as a fractional linear function of price. The results in the paper enrich the toolbox for the qualitative analysis of mathematical economics and business administration.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
The authors are grateful to the referees for their helpful comments and constructive suggestions.
W. Shuhe, “Differential equation model and chaos,” Journal of China Science and Technology University, pp. 312–324, 1999.View at: Google Scholar
Z. Xi-fan, C. Xia, and C. Yun-qing, “A qualitative analysis of price model in differential equations of price,” Journal of Shenyang Institute of Aeronautical Engineering, vol. 21, no. 1, pp. 83–86, 2004.View at: Google Scholar
S. Banerjee and W. A. Barnett, “Bifurcation analysis of Zellner’s marshallain macroeconomic model,” Journal of Economic Dynamics and Control, vol. 35, no. 9, pp. 1577–1585.View at: Google Scholar
L. Tanghong and L. Zhenwen, “Hopf bifurcation of price Reyleigh equation with time delay,” Journal of Jilin University, vol. 47, no. 3, 2009.View at: Google Scholar
W. Yong and Z. Yanhui, “Stability and Hopf bifurcation of differential equation model of price with time delay,” Highlights of Sciencepaper Online, vol. 4, no. 1, 2011.View at: Google Scholar
Y. Zhai, H. Bai, Y. Xiong, and X. Ma, “Hopf bifurcation analysis for the modified Rayleigh price model with time delay,” Abstract and Applied Analysis, vol. 2013, Article ID 290497, 6 pages, 2013.View at: Publisher Site | Google Scholar | MathSciNet
O. I. Adeyemi and L. C. Hunt, “Modelling OECD industrial energy demand: asymmetric price responses and energy saving technical change,” Energy Economics, vol. 29, no. 4, pp. 693–709, 2007.View at: Publisher Site | Google Scholar
L. Tanghong and Z. Linhua, “Hopf and codimension two bifurcation for the price Rayleigh equation with two time delays,” Journal of Jilin University, vol. 50, no. 3, pp. 409–416, 2012.View at: Google Scholar | MathSciNet
K. L. Cooke and Z. Grossman, “Discrete delay, distributed delay and stability switches,” Journal of Mathematical Analysis and Applications, vol. 86, no. 2, pp. 592–627, 1982.View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
J. Wei and S. Ruan, “Stability and bifurcation in a neural network model with two delays,” Physica D, vol. 130, no. 3-4, pp. 255–272, 1999.View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
J. Hale, Theory of Functional Differential Equations, Springer, New York, NY, USA, 2nd edition, 1977.View at: MathSciNet
D. D. Hassard, N. D. Kazarinoff, and Y. H. Wan, Theory and Applications of Hopf Bifurcation, Cambridge University Press, Cambridge, UK, 1981.