Abstract and Applied Analysis

Abstract and Applied Analysis / 2013 / Article
Special Issue

Dynamics of Delay Differential Equations with Its Applications

View this Special Issue

Research Article | Open Access

Volume 2013 |Article ID 962738 | https://doi.org/10.1155/2013/962738

Chuangxia Huang, Changlin Peng, Xiaohong Chen, Fenghua Wen, "Dynamics Analysis of a Class of Delayed Economic Model", Abstract and Applied Analysis, vol. 2013, Article ID 962738, 12 pages, 2013. https://doi.org/10.1155/2013/962738

Dynamics Analysis of a Class of Delayed Economic Model

Academic Editor: Zhichun Yang
Received14 Jan 2013
Accepted20 Feb 2013
Published04 Apr 2013


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 [27]. 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).(H1) The inequality hold: .(H2) The inequality hold: , where is a positive equilibrium point.(H3) 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 (H1) 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 (H2) hold: , then the equilibrium point (0,0) of system (15) is asymptotically stable.

Proof. As the inequality (H2) 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 (H3) 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 (H1)–(H3) 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 (H1)–(H3) 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 <