Abstract and Applied Analysis

Volume 2013 (2013), Article ID 962738, 12 pages

http://dx.doi.org/10.1155/2013/962738

## Dynamics Analysis of a Class of Delayed Economic Model

^{1}College of Mathematics and Computing Science, Changsha University of Science and Technology, Changsha, Hunan 410114, China^{2}School of Business, Central South University, Changsha, Hunan 410083, China

Received 14 January 2013; Accepted 20 February 2013

Academic Editor: Zhichun Yang

Copyright © 2013 Chuangxia Huang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

This investigation aims at developing a methodology to establish stability and bifurcation dynamics generated by a class of delayed economic model, whose state variable is described by the scalar delay differential equation of the form . At appropriate parameter values, linear stability and Hopf bifurcation including its direction and stability of the economic model are obtained. The main tools to obtain our results are the normal form method and the center manifold theory introduced by Hassard. Simulations show that the theoretically predicted values are in excellent agreement with the numerically observed behavior. Our results extend and complement some earlier publications.

#### 1. Introduction

Trade cycles, business cycles, and fluctuations in the price and supply of various commodities have attracted the attention of economists for well over one hundred years and possible more than thousands of years [1]. In the case of the most models discussed earlier in the literature, it is assumed that each economic agent has instantaneous information about its own as well as its rivals' behavior. Many authors often attributed these fluctuations to instantaneous information factors. This assumption is mathematically convenient but does not fully describe real economic situations in which there are always time delays between the times when information is obtained and the times when the decisions are implemented. In recent years it has been recognized in continuous-time economic dynamics that a delay differential equation is useful to describe the periodic and aperiodic behavior of economic variables [2–7]. Time delays usually cause the models to generate not only periodic cycles but also chaotic behavior for certain values of the shape parameter of the production function. With the infinite dimensionality created by a fixed-time delay, even a single first-order equation is transformed into an equation with a sufficient number of degrees of freedom to permit the occurrence of complex dynamics involving chaotic phenomena. This finding indicates that fixed-time delay models of a dynamic economy may explain various complex dynamic behaviors of the economic variables. For decades, a lot of effort has been devoted to deepen the understanding of economic complexity including chaotic behavior, stability, and basins of attractions. In [2], Matsumoto and Szidarovszky and in [3], Akio and Ferenc investigate a continuous-time neoclassical growth model with time delay and study the dynamics of the delayed model. In [4], Bélair and Mackey develop a model of price adjustment with production delays. In [5], Zhang et al. consider a differential-algebraic biological economic system with time delay and harvesting where the dynamics is logistic with carrying capacity proportional to prey population. Howroyd and Russel detect the stability conditions of delay output adjustment processes in a general N-firm oligopoly with fixed time delays [6]. Matsumoto and Szidarovszky introduce a fixed delay in production and a mound-shaped production function into the neoclassical one-sector growth model and show the birth of complex dynamics [2]. Considering about a delay in the production process, Li establish the single commodity price-inventory control model [7].

It is well known that persistent oscillations are one of the most ubiquitous forms by which economic phenomena may be observed [8]. Limit cycles are the simplest nonlinear phenomena, for example, they are the simplest example of how the interaction between economic forces may compel a system to abandon its steady state and start to steadily oscillate. Just as reported by Manfredi and Fanti [9], the detection of stable oscillations, for example, stable limit cycles, in continuous-time systems, is intimately related with the notion of Hopf bifurcation. Hopf bifurcation is important in economics. There are at least three reasons. First, it is always the outcome of a fully endogenous interaction between (nonlinear) economic forces. Second, it is a local bifurcation, thus much in spirit with the common belief of our science by which economic systems are generally close to their equilibrium state. Third, because it implieslocal oscillations, which are the normal route through which disequilibrium manifests itself when the equilibrating forces operating in the economy are relaxed (e.g., the adjustment process of a Walrasian market). For instance, when oscillations persist in a market normally in equilibrium (in the absence of stochastic and seasonal perturbations), it is very likely that these oscillations are the outcome of a Hopf bifurcation. So it is necessary for us to research the Hopf bifurcation in the economic system.

In the case of the most models discussed in the mentioned literature, only one delay appears in the models. Considering that the consumer memory plays an important role in the process of economic activities, just as pointed out by Li [7], bringing another delay might be better candidates for some purposes and would be of great interest. In this paper, we will generalize an economic model with the help of bringing two delays. In specifying how consumer behavior affects commodity demand, we will assume that the behavior is influenced not only by the instantaneous price, but also by the information regarding past prices. As the fact that the quantity supplied may not increase infinitely during the price increasing, similar to [7], we also assume the supply function to be fractional linear function. This investigation aims at developing a methodology to establish stability and bifurcation dynamics generated by the new delayed economic model.

The organization of the rest of this paper is as follows. In Section 2, we establish an economic model with two delays combing with the model considered in [7]. In Section 3, we take the delays as the parameters and use the distribution theory of the transcendental equation root [10] and the theory of Hopf bifurcation about functional differential equation [11] to discuss the stability of the equilibrium point for economic system and the existence of Hopf bifurcation. In Section 4, we apply the normal form theory and center manifold theory to investigate the bifurcation direction and the stability of periodic solution. In Section 5, an example with numerical simulations is arranged to illustrate the obtained results.

#### 2. Characterization of a Generalized Delayed Economic Model

Considering a single commodity market, the quantity of supplied and demanded can be regarded as the function of time, namely, and . The inventory and the level of inventory are recorded, respectively, as and .

Let denote the price at time , so that the rate of price increase is in proportion to the difference between and , namely, where is a positive real number depending on the speed of price adjustment, recording as

In the traditional cobweb model, demand function is a function of price. If we consider price as the only factor that influences the quantity demanded, there will be certain limitations to reflect the regularity of price change. We should consider other factors influencing the demand such as the rate of the price increase. In [7], Li assume the demand function as where , , represents the sensitive degree of consumers to the increase of commodity price. is the level of price relying on the rate of increase. In terms of consumer behavior, we consider a class of consumers who base buying decisions on the past prices and recent price. Therefore, in this paper, we introduce a delay in demand function and set up the demand function as shown later

In generally, supply function is monotone increasing about price, but consider that as price goes up, the supply could not unlimitedly increase, one can assume supply function as a fractional linear function as the following: where , , .

Noticing the delay in the production process, supply function should be a function of past price, therefore, we can introduce another delay in supply function and record it as follows: where , , , , are constants and , , , .

Substituting (2) to (1), calculating the derivation of both sides about time in (1), one can get
Substituting (4) and (6) to (7), we can establish a single commodity price inventory control model with two delays take the following form:
Let denote the Banach space of continuous and differentiable mapping from into equipped with the Supremum Norm for . The initial condition of (8) is
We will consider the following basic assumptions to further investigate the stability and bifurcation dynamics of model (8).(H_{1}) The inequality hold: .(H_{2}) The inequality hold: , where is a positive equilibrium point.(H_{3}) The inequality hold: , where is a positive equilibrium point.

#### 3. Stability Analysis and the Existence of Hopf Bifurcation

At first, we will show that system (8) has only one positive equilibrium point under some assumption. We state the following theorem.

Theorem 1. *If the inequality (H _{1}) holds: , then system (8) has only one positive equilibrium point.*

*Proof. *Without loss of generality, we may assume that
then model (8) can be rewritten as the following:

Assume to be the equilibrium point of system (8), one can show that is the equilibrium point of system (11). Therefore, one can obtain
Obviously, if , we have
Then (8) has only one positive equilibrium point. This completes the proof.

In real life, is the equilibrium price. With the help of coordinate translation system (11) can be further rewritten as the following form: Then the linearized system at (0,0) is The characteristic equation of the linearized system (16) at (0,0) takes the following form:

It is well known that the equilibrium (0,0) is asymptotically stable if all roots of the characteristic equation (17) have negative real parts. Now we reach the position to study the distribution of the roots of (17). We will consider three cases as follows: (Case 1): , ; (Case 2): , ; and (Case 3): .

*Case 1. *,
.

Proposition 2. *If , and the inequality (H _{2}) hold: , then the equilibrium point (0,0) of system (15) is asymptotically stable.*

*Proof. *As the inequality (H_{2}) holds, then the characteristic equation (17) turns to be
It is obvious that ; from Hurwitz criterion, all roots of this equation have negative real parts; therefore, the equilibrium point (0,0) of system (15) is asymptotically stable.

*Case 2 (). *If , , then the characteristic equation (17) takes the following form:
Let ; then we obtain the following results.

Lemma 3. *If the inequality (H _{3}) holds: , and , then has the only pair of purely imaginary roots , where*

*Proof. *If is a root of (19), then
Separating the real and imaginary parts, we have
Then we can obtain
The root of (23) can be expressed as follows:
Then we obtain
Let ; if is not the only pair of purely imaginary roots, then we have
From (17), one can get
That is to say, , which contradicts (25). This completes the proof.

Denote ; then the root of (19) satisfies: , .

Lemma 4. *If one chooses
**
and has the only pair of purely imaginary roots , then one has
*

*Proof. *Taking the derivative of with respect to , we get
Therefore
This completes the proof.

Based on the lemmas presented previous and the classical Hopf-Bifurcation-Theorem (see, [11, pages 245–249]), we have the following result.

Theorem 5. *Assume that the inequalities (H _{1})–(H_{3}) hold, then one has the following result.*(1)

*If , all roots of (19) have negative real parts. Namely, the equilibrium (0,0) of system (15) is locally asymptotically stable.*(2)

*If , (19) have a pair of purely imaginary roots , all the other roots have negative real parts. That is to say, system (15) undergoes a Hopf bifurcation at .*(3)

*If , (19) has roots of positive real parts. Namely the equilibrium (0,0) of system (15) is unstable.*

*Case 3 (). *We now discuss the stability of equilibrium (0,0) when , .

Let be the root of (17) and substitute it into (17), we have
Separating the real and imaginary parts, we have
Then we have

If (34) has no root or negative root, all the roots of (17) have negative real part. If it has positive roots, we know that the number of positive roots is finite. Denote them to be . From (33), we can get
Let , let , and let ; if , (16) has a pair of purely imaginary roots at , we can also prove that the purely imaginary roots are simple. Taking the derivative with respect to , we can get
Then, we have
Therefore, we have the following theorem.

Theorem 6. *Assume that the inequalities (H _{1})–(H_{3}) hold, then one has the following result.*(1)

*If , then all the roots of (17) have negative real parts. One can get that the equilibrium (0,0) of system (15) is locally asymptotically stable.*(2)

*If , and , then system (15) undergoes Hopf bifurcations at (0,0).*

#### 4. Direction and Stability of the Bifurcation

In this section, formula for determining the direction of Hopf bifurcation and the stability of bifurcation periodic solution of system (11) at , will be presented by employing the normal form method and center manifold theorem introduced by Hassard et al. in [12]. More precisely, we will compute the reduced system on the center manifold with the pair of conjugate complex, purely imaginary solution of the characteristic equation (17). By this reduction, we can determine the Hopf bifurcation direction, that is, to answer the question of whether the bifurcation branch of periodic solution exists locally for supercritical bifurcation or subcritical bifurcation.

Let , and let ; then is a critical value of Hopf bifurcation of system (11). With the translation , let , , we rewrite system (15) as follows: where . Furthermore, we can obtain the linear system of (38) as mentioned later while the nonlinear term is

Let , let , let , let , and let , where is a linear operator,

As is a one-parameter family of bounded linear operator in , by the Rise Representation Theorem, there exists a matrix whose components are bounded variation functions , such that , . Actually, we can take , where is Dirac delta function.

Next, we define where . We can rewrite (38) as where .

Denote , , , .

For , the adjoint operator of is defined as where is the transpose of .

We define the bilinear form where and .

We easily obtain that are eigenvalues of (17) by the translation . Then we have the following lemma.

Lemma 7. * is the eigenvector corresponding to and is the eigenvector of corresponding to . , , where
*

*Proof. *Letting be the eigenvector of corresponding to , we have
Calculating (48), we obtain
where is constant.

From (43), we can obtain
Therefore, one can show that

Letting be the eigenvector of corresponding to , based on (48), (50), we have
It is easy to obtain that and .

Now we compute as the following:
Then we have
On the other hand,
therefore, we have . This completes the proof.

In the remainder of this section, we use the same notation as in [12]. We first compute the center manifold at .

Let be the solution of (44) when , and define We have where On the center manifold, we have In fact, , are local coordinates for center manifold in the direction of , . Noting that is real if is real, we consider only real solutions in this paper. We rewrite (58) as follows: where From (56), we have Then we have It is easy to obtain We can also obtain Expanding (66) and comparing the coefficients, we obtain

We still need to compute and . Noticing that we rewrite as the following: where

As , we can get Comparing the coefficients with (69), we can get From (69), , where . Then we have From (43) and (71), we have where . Computing (73), we have where , .

, can be determined by setting in . As and are continuous on , then we have From (43), we obtain Substitute (74)–(77) to (78), we can obtain Therefore, we obtain Based on the previous analysis,we can see that each in (67) is determined by the parameters in system (11). Thus we can compute the following values: which determine the quantities of bifurcating periodic solutions in the center manifold at the critical value ; that is, we thus have the following theorem.

Theorem 8. *Using (81), one can compute the values of , , , . Therefore, we can answer the question of whether the bifurcation branch of periodic solution exists locally for supercritical bifurcation or subcritical bifurcation as the following.*(i)* determines the directions of the Hopf bifurcation: if (), then the Hopf bifurcation is supercritical (subcritical) and the bifurcating periodic solutions exists for .*(ii)* determines the stability of the bifurcating periodic solutions: the bifurcating periodic solutions are stable (unstable) if ().*(iii)* determines the period of the bifurcating periodic solutions: the period increases (decreases), if ().*

#### 5. Numerical Simulation Example

*Example 9. *Consider a single commodity market model as follows:
where denotes the price at time , and we choose , , , , , . Using Theorem 1, a quick computation revealed that .

If one choose , , from Lemma 3, we can calculate that , , and is the the critical value for Hopf bifurcation. It follows from Theorem 5 that if , is asymptotically stable. System (82) undergoes a Hopf bifurcation at . These conclusions are verified by the numerical simulation in Figures 1, 2, 3, and 4.

When , we can choose ; from (34) and (35), we can calculate that , ; from (37), we have . It follows from Theorem 6 that is asymptotically stable as far as and system (82) undergoes a Hopf bifurcation when . Furthermore, using Lemma 7, after simple computation, we have
Then we can obtain
It is easy to obtain
Therefore, the Hopf bifurcation of system (82) is supercritical, and the bifurcation periodic solutions are stable. This conclusions can be verified by the numerical simulation in Figures 5, 6, 7, and 8.

#### 6. Conclusion

The delay in production process has been considered in many economic models. However, few of them considered the delay of consumer consumption. In this paper, we first establish a class of economic models with two delays; the delay in production process and consumer consumption are both considered, and then the dynamics of this system have been investigated. Specifically, stability of equilibrium point and the existence of Hopf bifurcation, in great detail, are studied. We have derived some sufficient conditions to ensure that all the characteristic roots have negative real parts. We also show that a Hopf Bifurcation will occur once some parameters pass through critical values; that is, a family of periodic orbits bifurcates from the positive equilibrium point. At last, the direction of Hopf bifurcation and the stability of the bifurcating periodic orbits are discussed by applying the normal form theory and the center manifold theorem. Simulations show that the theoretically predicted values are in excellent agreement with the numerically observed behavior. Our results extend and complement some earlier publications.

#### Acknowledgments

The authors are extremely grateful to the editor and the anonymous reviewers for their constructive and valuable comments, which have contributed much to the improvement of this paper. This work was jointly supported by the National Natural Science Foundation of China under Grants nos. 11101053, 70921001, and 71171024, the Key Project of Chinese Ministry of Education under Grant no. 211118, the Excellent Youth Foundation of Educational Committee of Hunan Provincial no. 10B002, and the Scientific Research Funds of Hunan Provincial Science and Technology Department of China.

#### References

- M. L. Weidenbaum and S. C. Vogt, “Are economic forecasts any good?”
*Mathematical and Computer Modelling C*, vol. 11, pp. 1–5, 1988. View at Google Scholar · View at Scopus - A. Matsumoto and F. Szidarovszky, “Delay differential neoclassical growth model,”
*Journal of Economic Behavior and Organization*, vol. 78, no. 3, pp. 272–289, 2011. View at Publisher · View at Google Scholar · View at Scopus - M. Akio and S. Ferenc, “Nonlinear delay monopoly with bounded rationality,”
*Chaos, Solitons & Fractals*, vol. 45, no. 4, pp. 507–519, 2012. View at Publisher · View at Google Scholar - J. Bélair and M. C. Mackey, “Consumer memory and price fluctuations in commodity markets: an integrodifferential model,”
*Journal of Dynamics and Differential Equations*, vol. 1, no. 3, pp. 299–325, 1989. View at Publisher · View at Google Scholar · View at MathSciNet - G. Zhang, B. Chen, L. Zhu, and Y. Shen, “Hopf bifurcation for a differential-algebraic biological economic system with time delay,”
*Applied Mathematics and Computation*, vol. 218, no. 15, pp. 7717–7726, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - T. D. Howroyd and A. M. Russell, “Cournot oligopoly models with time delays,”
*Journal of Mathematical Economics*, vol. 13, no. 2, pp. 97–103, 1984. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - L. Li, “Several differential equation models in economic system,”
*Journal of the Graduate School of the Chinese Academy of Science*, vol. 20, no. 3, pp. 273–278, 2003 (Chinese). View at Google Scholar - J. Guckenheimer and P. Holmes,
*Non-linear Oscillations, Dynamical Systems and Bifurcation of Vector Fields*, Springer, New York, NY, USA, 1983. View at MathSciNet - P. Manfredi and L. Fanti, “Cycles in dynamic economic modelling,”
*Economic Modelling*, vol. 21, no. 3, pp. 573–594, 2004. View at Publisher · View at Google Scholar · View at Scopus - S. Ruan and J. Wei, “On the zeros of transcendental functions with applications to stability of delay differential equations with two delays,”
*Dynamics of Continuous, Discrete & Impulsive Systems A*, vol. 10, no. 6, pp. 863–874, 2003. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. Hale,
*Theory of Functional Differential Equations*, Springer, New York, NY, USA, 1977. View at MathSciNet - B. D. Hassard, N. D. Kazarinoff, and Y. H. Wan,
*Theory of Applications Hopf Bifurcation*, London Mathematical Society Lecture Note Series, Cambridge University Press, Cambridge, UK, 1981. View at MathSciNet