Frontiers journals are at the top of citation and impact metrics

Original Research ARTICLE

Front. Appl. Math. Stat., 29 September 2016 |

A Heterogeneous Agent Model of Asset Price with Three Time Delays

  • 1Department of Economics, International Center for further Development of Dynamic Economic Research, Chuo University, Hachioji, Japan
  • 2Department of Applied Mathematics, University of Pécs, Pécs, Hungary

This paper considers a continuous-time heterogeneous agent model of a financial market with one risky asset, two types of agents (i.e., the fundamentalists and the chartists), and three time delays. The chartist's demand is determined through a nonlinear function of the difference between the current price and a weighted moving average of the delayed prices whereas the fundamentalist's demand is governed by the difference between the current price and the fundamental value. The asset price dynamics is described by a nonlinear delay differential equation. Two main results are analytically and numerically shown:

(i) a single delay destabilizes the market price and generates cyclic oscillations around the equilibrium;

(i) under multiple delays, stability loss and gain repeatedly occur as the length of the delay increases.

1. Introduction

Persistent volatility is a prominent characteristic feature of financial markets. It is, however, well-known that the efficient market hypothesis cannot account for the discrepancy between the observed market price and the fundamental value of the asset. Recently, heterogeneous agent models (HAMs) have been developed to explain a wide variety of financial market behavior such as temporary bubbles, sudden market crashes and price resistance in discrete-time as well as in continuous-time framework. In his survey of recent developments of HAMs, Hommes [1] discusses that nonlinear discrete-time HAMs can generate various dynamics ranging from cyclic fluctuations to chaotic behavior. Among others there is Chiarella et al. [2] who propose a discrete-time HAM with a moving average (MA) rule having various memory lengths. Their main finding is that the length of the MA rule can be a source of complicated dynamics in destabilized financial markets in deterministic and stochastic processes of difference equations. On the other hand, continuous-time HAMs are also examined in various ways in which dynamics are described by ordinary or delay differential equations. They have a long history since Zeeman [3] and Beja and Goldman [4]. Moreover, Chiarella [5], which is a development of Beja and Goldman [4], shows that the market price tends to a stable limit cycle under a nonlinear demand function of the risky asset when the equilibrium is unstable. More recently, He and Zheng [6] reconstruct the discrete-time model of Chiarella et al. [2] in a continuous-time framework in which the expected price is formed with a moving average of the past (delay) prices. Their main result concerns a double edge effect on the stability caused by the length of the memory or delay: an increase in delay can destabilize the market price and also stabilize it. In their model as well as in subsequent studies on continuous-time HAMs such as He and Li [7] and Xu et al. [8], this interesting phenomenon is obtained under the assumption that infinitely many past price data are available at no charge of cost. Needless to remember the well-known line in economics, “there is no such thing as a free lunch,” it is to go too far to get necessary information for nothing. At least two ways are possible to render this extreme assumption to more realistic one, the first is to introduce information cost associated with collecting the price data and the second is to limit availability of the past prices. In the present paper, we take the second way and then reconsider the price dynamics when the chartists forecast an expected price, using only a limited information on past prices.

Following the framework of Dibeh [9], this paper constructs a dynamic HAM of a speculative asset with three delays (i.e., three past prices) and investigates the effect of time delays on the asset price dynamics. Dibeh [9] conducts mainly with numerical simulations. Analytical developments of a one-delay version and a two-delay version of Dibeh's model are already presented by Qu and Wei [10] and Matsumoto and Szidarovszky [11], respectively. The model we analyze is a continuation of these preceding studies. One of our main findings is that stability switches from stability to instability as well as from instability to stability can occur in three delay framework. That is, the double edge effect can also be observed in continuous-time HAMs with a finite number of delay prices.

The structure of the paper is as follows. Section Model constructs a HAM model with fundamentalists and chartists. The stability switching curve is analytically derived, which divides the delay space into stability and instability regions. Section Dynamics has three subsections, in each of which we conduct a stability and bifurcation analysis under different circumstances. Section Concluding Remarks concludes the paper with further research directions.

2. Model

We consider an asset pricing model with one risky asset and two traders, fundamentalists and chartists. Let p(t) be the price of the risky asset at time t. Since the fundamentalists believe that p(t) eventually converges to f, the fundamental (i.e., equilibrium) price of the asset, they will sell the asset if p(t) > f and buy it if p(t) < f. The simplest form of the demand function of the fundamentalist is


where m ∈ (0, 1) is the fraction of the fundamentalists in the market. On the other hand, the chartists base their decisions of market participation on the price trend of the asset. Their demand function is


where s(t) is a weighted average of the past price trends formulated at time t,

s(t)=i = 1nαi[p(tτi  1)p(tτi)]

where τi ≥ 0 is a time delay with τ0 = 0 and αi ≥ 0 denotes a weight of a price change satisfying i=1nαi=1. The demand function g(s) is assumed to be hyperbolic tangent as in Chiarella [5] and Dibeh [9],


The market demand is the sum of demands of the fundamentalists and chartists,


The price dynamics follows Dibeh's formulation in which the growth rate of price is determined by the market demand,


that is written as a nonlinear delay differential equation,

p˙(t)=(1m)p(t)tanh[s(t)]+mp(t)[fp(t)].    (1)

Notice that s(t) can be rewritten as



p¯n(t)=j = 1nηjp(tτj)


ηj=αjαj + 1α1 for 1jn1 and ηn=αnα1(=1j = 1n  1ηj).

Since p-n(t) is considered to be the weighted (moving) average of the past n prices, the chartists believe that the price will rise (fall) when the current price is above (below) the average of the n delay prices. If infinitely many price data is available with finite or infinite memory length, the moving average can be presented by continuously distributed time delay such as

pξ(t)=t  ξtη(s)p(s)ds

where 0 < ξ ≤ ∞ is the memory length and η(s) is a weighting function. This form is used in He and Zheng [6].

Concerning the specification of s(t), Qu and Wei [10] examine the case of n = 1 in which the trend includes one delay price with α1 = 1


It is shown first that there exists a threshold value of the delay at which stability is lost and second that the stationary point can be stable for smaller values and bifurcates to a limit cycle for larger values. Matsumoto and Szidarovszky [11] consider the case of n = 2 in which two delay prices are used to formulate the trend,






With two delay prices, the coefficients are positive and add up to unity in interpolation (i.e., α1 > α2) whereas one coefficient is negative, the other greater than unity and the sum is also unity in extrapolation (i.e., α1 < α2). It is demonstrated that the two delay model has a double edge effect with which stability losses and gains can repeatedly occur.

In this study, we move one more step forward and draw attention to the case of n = 3 for which the average trend is





η1=α1α2α1,η2=α2α3α1 and η3=α3α1.

For the sake of convenience, we adopt the following form,


The unique stationary point of dynamic Equation (1) is identical with the equilibrium price,

pe=f=p(t)=p(tτj) for all t0 and j=1,2,3.

It is clear that at the stationary point there is no price trend, se = 0. Linearizing the nonlinear dynamic equation in a neighborhood of the stationary point yields

p˙δ(t)=αpδ(t)+β1pδ(tτ1)+β2pδ(tτ2)β3pδ(tτ3)    (2)

where pδ(t)=p(t)-pe and the coefficients are defined as


With condition i = 13αi=1, it can be verified that


The locus of α3 = 0 divides the first quadrant of the (α1, α2) plane into two parts, upper and lower right triangles, and condition α3 ≥ 0 eliminates the upper one. The lower triangle is further divided into six subparts by the three loci of α2 = α1, α3 = α2, and α3 = α1. Those divisions are depicted in Figure 1 in which the relations among the magnitudes of α1, α2, and α3 are determined in the following way,

α1>α2>α3 in region I,α3>α2>α1 in region IV,α1>α3>α2 in region II,α2>α3>α1 in region V,α3>α1>α2 in region III,α2>α1>α3 in region VI.

Figure 1. Divisions of the (α1, α2) plane with α3 ≥ 0.

Region I is colored in gray and we will limit our analysis to it below. The line of α2 = (2 − 3m)/2(1 − m) − α1 is shown to be dotted and downward sloping1. It is also checked that β1 = 0, β2 = 0, and β3 = 0 hold, respectively, on the loci of α2 = α1, α3 = α2, and α3 = 0. On these boundaries the three-delay Equation (2) is reduced to one of three two-delay equations, depending on which βi value becomes zero. Matsumoto and Szidarovszky [11] consider a two delay equation that corresponds to the one with β3 = 0. Their analytical method can be applied to the other two equations as well. In this study, we will confine our attention to the regions in which βi > 0. Although, the parametric region I is a smaller part of the whole feasible region, the three delay model restricted to region I has more freedom concening the parametric choice than the two delay model that is defined only on the line α2 = 1 − α1.

2.1. Construction of Stability Switching Curves

Substituting an exponential solution pδ(t)=eλtu into the linearized Equation (2) yields the corresponding characteristic equation

λαβ1eλτ1β2eλτ2+β3eλτ3=0    (3)

Dividing the left hand side of Equation (3) by λ − α and denote the result by a(λ),

a(λ)=1+a1(λ)eλτ1+a2(λ)eλτ2+a3(λ)eλτ3=0    (4)



Since λ = 0 is not a solution of Equation (3), a pair of pure imaginary roots must exist if stability switch occurs. We thus assume λ = , ω > 0 is a solution of Equation (3). For λ = iω,

a1(iω)=αβ1+iβ1ωα2+ω2 and |a1(iω)|=|β1|α2+ω2,
a2(iω)=αβ2+iβ2ωα2+ω2 and |a2(iω)|=|β2|α2+ω2


a3(iω)=αβ3+iβ3ωα2+ω2 and |a1(iω)|=β3α2+ω2.

Each term of a() may be viewed as a vector in the complex plane. Hence solving a() = 0 analytically is equivalent to constructing a quadrangle geometrically as shown in Figure 2 in which these four vectors are arranged from head to tail.


Figure 2. |1| and |aj(iω)| for j = 1, 2, 3 form a quadrangle.

The determination of the stability switching surface in the (τ1, τ2, τ3) plane is challenging2. In order to make the problem manageable, we reduce the three-delay equation to a two-delay equation by fixing the value of τ3 at a certain positive level and then construct the stability switching curve in the (τ1, τ2) plane in which we can apply the method used earlier in Matsumoto and Szidarovszky [11], and Gu et al. [14]. Notice that vector 1 is connected to vector a3(iω)e-iωτ3 in Figure 2. If the sum of these vectors is denoted by ad(, τ3),


then a(λ) with λ = can be rewritten as


As already shown in Figure 2, a() = 0 means that these three vectors in a() must form a triangle having the dotted base with interior angles, θ1 and θ2. Sufficient and necessary conditions for forming a triangle are given by

(i)fd(ω)=|a1(iω)|+|a2(iω)||ad(iω,τ3)|0,(ii)f1(ω)=|a2(iω)|+|ad(iω,τ3)||a1(iω)|0,(iii)f2(ω)=|ad(iω,τ3)|+|a1(iω)||a2(iω)|0.    (5)

Each inequality condition implies that the length of any segment of the triangle is not greater than the sum of the lengths of the remaining two segments. Define the crossing frequency set Ω of all ω > 0 such that a() = 0 holds for at least one delay combination (τ1, τ2, τ3) ≥ 0. For ω = 0 in a() leads to

a(iω)|ω = 0=1α(α+β1+β2β3)<0.

This inequality implies 0 ∉ Ω.

By the law of cosine, we can determine the values of θ1 and θ2 of the triangle in Figure 2,




Vertices A and B maybe located above or below the horizontal axis and the slope of segment AB can be positive or negative. So we have eight possibilities to construct a quadrangle. However, a simple geometric consideration shows that there are only two different possibilities:




Solving these equations for τ1 and τ2 yields the threshold values of the delays,

τ1±(ω,k)=1ω[arg(a1(iω))arg(ad(iω,τ3))+(2k1)π     ±θ1(ω)]    (6)


τ2(ω,n)=1ω[arg(a2(iω))arg(ad(iω,τ3))+(2n1)π    θ2(ω)].    (7)

The stability switching curves consist of two sets of parametric segments,

L1(k,n)={τ1+(ω,k),τ2(ω,n)} for k,n=0,1,2,...    (8)


L2(k,n)={τ1(ω,k),τ2+(ω,n)} for k,n=0,1,2,...    (9)

where the parameter ω runs through the crossing frequency set Ω.

Substituting τ1 = τ2 = τ3 = 0 into Equation (3) presents


This inequality implies that the steady state with no delays is stable. Next, assuming that τ1 = τ2 = 0, we examine the existence of a threshold value of τ3 at which the system loses stability. Without loss of generality, assuming λ = , ω > 0 and then substituting it into Equation (3) with τ1 = τ2 = 0 reduce the characteristic equation to


We denote α + β1 + β2 by Δ where


and break down the characteristic equation into the real and imaginary parts,

Δ+β3cosωτ3=0ωβ3sinωτ3=0.    (10)

Moving Δ and ω to the right hand side and adding the squared equations give








Notice that the dotted curve in Figure 1 is described by α2 = φ(α1). It is clear first that

β3+Δ0 according to φ(α1)α2

and second that

β3+Δ<0 if m2/3.

This leads to the following result.

Proposition 1. Given τ1 = τ2 = 0, the solution of delay system (1) is stable for any τ3 ≥ 0 if the fundamentalists dominate over the chartists in the sense that m ≥ 2/3.

If β3 + Δ > 0, then there is a ω such as ω-=β32-Δ2>0. In this case, loss of stability for λ=iω- can be shown in the following way. The characteristic equation with τ1 = τ2 = 0 can be written as

β3eλτ3=β1+β2+αλ.    (11)

Taking λ as a function of τ3 and differentiating it with respect to τ3 yield


where the right hand side of Equation (11) are used. Substituting λ = and taking the real part of the resulting expression give

d[Reλ]dτ3|λ = iω=Re[λ(β1+β2+αλ)1τ3(β1+β2+αλ)]=Re[ω2+iω(β1+β2)1τ3(β1+β2+α)iω1+τ3(β1+β2+α)iω1+τ3(β1+β2+α)iω]=ω2[1τ3(β1+β2+α)]2+ω2>0.

This inequality implies that all roots cross the imaginary axis at from left to right as τ3 increases so stability is lost. Correspondingly, solving the first equation of Equation (10) for τ3 presents the threshold value

τ¯3(α1,α2,m)=1ω¯cos1(Δβ3)>0    (12)

and the system is stable for τ3<τ-3 and unstable otherwise3. We summarize the results obtained so far.

Proposition 2. Given τ1 = τ2 = 0, the solution of delay system (1) is stable for τ3<τ-3(α1,α2,m), loses stability for τ3=τ-3(α1,α2,m) and becomes unstable for all τ3>τ-3(α1,α2,m) if the following condition holds,


whereas it is always stable for any τ3 > 0 otherwise.

It is to be noticed that with τ1 = τ2 = 0, the three delay model is reduced to a one delay model and thus this stability result is essentially the same as the one of Qu and Wei [10]. In the next section we will proceed to the case of τi > 0 for i = 1, 2, 3.

3. Dynamics

Central to this study is the problem of dynamics associated with three fixed delays. This problem is complex and it seems unlikely that an analytical approach to nonlinear delay differential Equation (1) would generate fruitful results on dynamics as it might be difficult to solve the equation. In order to understand global dynamics with delays better, it is useful to take a numerical approach. To this end, we make several simplifying assumptions and then numerically examine dynamics of the asset price.

Assumption 1. f = 5, m = 0.4 and τ3 = 5.

Assumption 2. αi for i = 1, 2, 3 satisfy α1 > α2 > α3 > 0.

Assumption 3. τi for i = 1, 2, 3 satisfy τ1 < τ2 < τ3.

Since the parameter values in Assumption 1 are determined only for computational convenience, the main results to be obtained below hold for any other values, provided that m < 2/3. Assumption 2 implies that the more recent trend has the more weight. This is not necessary but simplifies the analysis. Under m = 0.4 of Assumption 1, region I is located above the line α2 = φ(α1) which is the dotted downward-sloping line passing through point (1/3, 1/3) in Figure 1. Due to Proposition 2, the stationary point is always stable in region I where α2 > φ(α1) if τ1 = τ2 = 0. Thus it is easy to see how increasing values of these delays affect stability of the steady state. Assumption 3 imposes a natural ordering among the three delays, that is, an older price has a longer delay.

The rest of this section is divided into three subsections. We will look at the τ1-effect on dynamics caused by a change in τ1 in Section Delay Effect I: τ1-effect, explore the τ2-effect in Section Delay Effect II: τ2-effect, and focus on considerations of the τ3-effect in Section Delay Effect III: τ3-effect.

3.1. Delay Effect I: τ1-effect

We start with the τ1-effect and present some numerical examples to see how different values of τ1 affect dynamics of Equation (1). We first check the triangle conditions. Subtracting the third equation from the second equation Equation in (5) yields

f1(ω)f2(ω)=6f(1m)α2+ω2(α213)    (13)

implying that the relative location of the curves of f1(ω) and f2(ω) depends on whether the value of α2 is larger or less than 1/3. To simplify the analysis, we assume that the value of α2 is kept constant at 1/3;

Assumption 4. α2 = 1/3

Assumption 2 implies ηi > 0 for all i, that is, the weights of the delay prices to calculate the moving average are positive. Their relative magnitudes depend on the selected values of αi and thus there are many possibilities. Among them, we will limit analysis to only two cases, η1 = η2 > η3 if α1 < 1/2 and η3 > η1 = η2 if α1 < 1/2 since Assumption 4 makes η1 = η2 always4. With Assumption 4, f1(ω) = f2(ω) holds regardless of the value of α1. Further, fi(ω) > 0, i = 1, 2, for ω > 0 which will be numerically confirmed below. Consequently, it is enough to verify only the location and the shapes of the fd(ω) curve. We consider dynamics over an interval [α10,2/3] of α1. The left hand side extreme value, α100.356 is numerically determined so as to make the maximum of fd(ω) equal to zero. Thus for α1<α10, fd(ω) < 0 for all ω > 0, implying that one of the triangle conditions is always violated. The other extreme value, α1 = 2/3, together with α2 = 1/3 makes α3 = 0 via α2 = 1 − α1. For α1 > 2/3, the nonnegative condition α3 ≥ 0 is violated. In the following numerical analysis, α1 is increased from 6/15 to 9/15 with an increment of 1/15. The resultant four fd(ω) curves are colored in black in Figure 3 where the red curve with α1 = 1/2 and the blue curves with those extreme values are also depicted. The corresponding f1(ω) = f2(ω) curves are illustrated in the dotted curves with the same color and are seen to be positive for all ω ≥ 0. In the gray region below the lower blue curve, no stability switch occurs and thus the steady state is stable. We do not consider the τ1-effect there. On the other hand, the gray region above the upper blue curve is not feasible and thus eliminated from consideration. Figure 3 also illustrates various shapes of the fd(ω) curves and shows that increasing α1 shifts the curve upward and removes its unevenness. Indeed, the lower blue curve with α1=α10 located at the bottom of the yellow region is high-wave shaped while the higher blue curve with α1 = 2/3 at the top is monotonically downward-sloping. In addition, as the red curve passes through the origin, we then see that the curves below the red curve have two intersections with the horizontal axis and the curves above have only one intersection.


Figure 3. The fd(ω) curves with various values of α1.

To be more specific, we pick up α1 = 6/15 and construct a stability switching curve. The corresponding fd(ω) curve is the black one just above the lower blue curve in Figure 3 and intersects the horizontal axis twice at

ω¯10.438 and ω¯20.582.

The triangle conditions in Equation (5) hold for ωΩ=[ω-1,ω-2]. We now ready to derive the stability switching curves, L1(k, n) and L2(k, n) defined in (8) and (9). Taking k = n = 1, we obtain the points of τ1±(ω,1) and τ2(ω,1) by changing values of ω with increment 0.0002 inside the interval [ω-1,ω-2] and then plot these points in the (τ1, τ2) plane. In Figure 4A, segment L2(1, 1) is illustrated as a real curve in the region above the diagonal and segment L1(1, 1) is a dotted curve in the region below5. Two segments form an elliptical-shaped closed curve having the starting point S=(τ1+(ω-1,1),τ2-(ω-1,1)) and the ending point E=(τ1-(ω-2,1),τ2+(ω-2,1)) on the diagonal where, for i = 1, 2,

τi±(ω¯1,1)=τs8.478 and τi±(ω¯2,1)=τe2.789.

Figure 4. Dynamics for α1 = 6/15 and α2 = 1/3: (A) Stability region, (B) Bifurcation diagram.

Since the delays have two constraints, τi ≤ 5 for i = 1, 2 and τ1 < τ2 due to Assumptions 1 and 3, the feasible region should be above the diagonal line and subject to τ2 ≤ 5. It is colored in yellow and further sub-divided into two sub-regions by the L2(1, 1) segment that intersects two lines, one is the diagonal at (τe, τe) and the other is the horizontal line at τ2 = 56.

It is clear from the yellow region that stability is preserved for τ1 and τ2 such as 0<τ1<τe and τ1<τ2<τe. When the varying pair (τ1, τ2) with τ2>τe crosses the L2(1, 1) segment, we have the stability switching as the real part of at least one eigenvalue turns to be positive and then the stationary point loses stability. To see this, we perform a simulation by increasing the value of τ1 along the dotted horizontal line at τ2 = 4 that crosses the L2(1, 1) curve at point (τ1c,4) with τ1c1.6407. The numerical results are summarized in Figure 4B in which a bifurcation diagram with respect to τ1 is depicted. For τ1τ1c, a pair of (τ1, 4) stays within the stability region. So the stationary point p* = f is stable and the corresponding part of the bifurcation diagram is a horizontal line at p(t) = f. Once τ1>τ1c, the pair is in the instability region above the L2(1, 1) segment. The stationary point, therefore, loses stability and bifurcates to a limit cycle having two extreme values (i.e., maximum and minimum). The existence of the limit cycle is confirmed by the Hopf bifurcation theorem. Figure 4B further implies that the limit cycle becomes larger as τ1 gets larger and the stability is never regained for τ1 ≤ 48. The τ1-effect is summarized as follows:

Proposition 3. Stability of dynamic Equation (1) with respect to τ1 depends on the selected value τ-2 of τ2;

If τ¯2τe,then the steady state is stable for τ1<τ¯2;Ifτ¯2>τe,then there is the threshold value τ1c such that thesteady state is stable for τ1<τ1c and loses stability forτ1τ1c.

where τ1c is the τ1-point of the intersection of the switching curve and the horizontal line at τ-2.

Taking various values of α1 from the interval [α10,2/3] and constructing the stability switching curves, we detect the τ1-effect more. In Figure 5A, the switching curves with various values of α1 are illustrated and their end-points on the diagonal are denoted by the green dots. The upper most blue curve has α1=α10 and the lower most blue curve has α1 = 2/3 while the curve shifts downward as α1 increases from α10 to 2/3. As far as the curve is single-valued in τ2, we have essentially the same result on the τ1-effect as in the case of τ1 = 6/15: fixing the value of τ2 at some level and increasing the value of τ1, the stability is preserved for τ1 below the threshold value τ1c and it is lost for τ1>τ1c. As seen in Figure 5A, this threshold value becomes smaller as the value of α1 becomes larger9, a larger value of α1 reinforces the destabilizing τ1-effect.


Figure 5. Delay destabilizing and stabilizing effects: (A) Stability switching curves; (B) Bifurcation diagram.

3.2. Delay Effect II: τ2-effect

Taking τ1 fixed, we are next concerned with the τ2-effect, how changes in the value of τ2 affect dynamics. We first return to Figure 4A in which both values of α1 and α2 are fixed at 6/15 and 1/3, respectively. The downward-sloping shape of the switching curve in the yellow region indicates that changing the value of τ2 with constant τ1 may generate similar effect to the τ1-effect. Let τ1m denote the τ1-value of the intersection of the stability switching curve and the horizontal line at τ2 = 510. If τ1 is selected to be less than τ1m, then the steady state is stable for any τ2 that should be above the diagonal and ≤5. If τ1m<τ1<τe, then the steady state is stable for τ2<τ2c and bifurcates to a limit cycle for τ2>τ2c where τ2c is defined in the same way as τ1c and depends on the selected value of τ1. One more case exists and is specific to the τ2-effect, namely, if τ1>τe, any feasible pair of τ1 and τ2 is in the unstable region so the steady state is unstable.

We now turn attention again to Figure 5A and consider the τ2-effect when α1 is increased but α2 is still fixed at 1/3. It is already seen that the stationary point is stable in the region left to the stability switching curve and increasing the value of α1 shifts the switching curve leftward. The stability switching curve with α = α0 is not depicted as it is located outside Figure 5A. The shift of the curve implies that a larger value of α1 reinforces the destabilizing τ2-effect in the sense that it makes the stability region smaller. In addition, we have interesting dynamic phenomenon when the switching curve with α1 being close to 2/3 has forward- and backward-bending segments. This shape suggests the stability loss and gain with respect to τ2. If the initial point of (τ1, τ2) is selected on the diagonal such that the stationary point is stable, then the vertical line standing at this value of τ1 denoted as τ10=0.85 crosses the stability switching curve three times at

τ2a1.014,τ2b1.687 and τ2c3.359

as illustrated in Figure 5A. That is, the stability is lost at the first intersection and a limit cycle emerges. The cycle is getting larger and then smaller to merge the steady state at the value of the second intersection at which the stability is regained. Further increasing τ2 crosses the stability switching curve again at which the stability is lost again. The corresponding bifurcation diagram is given in Figure 5B in which the stability losses and gains are shown. This τ2-effect is summarized as follows:

Proposition 4. With a value of τ1 being close to 2/3, the increasing delay τ2 has a double edge effect in a sense that it can destabilize the market price to generate a limit cycle, then stabilize it and later destabilize it again.

3.3. Delay Effect III: τ3-effect

We now look into the τ3-effect, the effect caused by changing the value of τ3. To this end, α1 = (6 − 0.1)/12 is taken and dynamics with three different values of τ3, τ3 = 5, τ3 = 6, and τ3 = 7, are considered. The corresponding switching curves are illustrated in different colors in Figure 6A (that is, a blue curve with τ3 = 5, red curves with τ3 = 6, and black curves with τ3 = 7). As τ3 increases, it is observed first that the feasible regions of τ1 and τ2 increase since the upper bound of τ2 becomes larger and second that the downward-sloping switching curve shifts upward, which implies an enlargement of the stability region. It is further observed that island-shaped switching curves emerge for τ3 = 6 and τ3 = 7 and the steady state is unstable inside them. To detect the reason why the island-shaped region is born under larger values of τ3, we draw the green f1(ω) = f2(ω) curve and the red fd(ω) curve for τ3 = 7 in Figure 6B. It is seen that f1(ω) = f2(ω) > 0 for ω > 0 and the larger τ3 value makes the second hump of the fd(ω) curve high enough to cross over the horizontal line. As a result, the red curve intersects the horizontal line four times at


Figure 6. τ3-effects on switching curves: (A) Stability switching curves; (B) Triangle conditions.

In consequence, the triangle conditions (5) hold in the two intervals, [ω1, ω2] and [ω3, ω4]. It is verified that the switching curve defined on [ω1, ω2] is downward-sloping and the one on [ω3, ω4] forms an island shape. Both curves are colored in black. The existence of the smaller red islands indicates that the fd(ω) curve also intersects the horizontal axis four times even when τ3 = 6. A critical value of τ3 for the birth of the island is somewhere between τ3 = 5 and τ3 = 6, for which the maximum of the second hump of fd(ω) becomes zero.

Proposition 5. As far as τ3 is relatively small, increasing the value of τ3 stabilizes the steady state by shifting the downward-sloping switching curve upward whereas for a relatively larger value, increasing value of τ3 additionally generates an island-shaped switching curve within which the steady state is unstable.

To see the τ3-effect on dynamics, we perform two simulations. In the first simulation, we increase the value of τ1 along the horizontal line at τ2=τ20(4.356). This line successively crosses the stability switching curves, the blue curve at τ1a(0.972), the red curve at τ1b(1.372) and the black curve at τ1c(=2). Although, these bifurcating values are denoted as black dots, they are not labeled in Figure 6A to avoid confusion. The corresponding bifurcation diagrams are illustrated in Figure 7A in which the black curve is depicted first, the red curve is then put on it and finally the blue curve is further placed upon. Each diagram has a qualitatively similar shape as the one shown in Figure 4B, that is, the steady state loses stability at the bifurcation value and a limit cycle emerges for a larger value of τ1. Thus larger τ3 value does not alter qualitative aspects of the τ1-effect, however, does affect its quantitative aspects because the bifurcation value of τ1 becomes larger as τ3 increases. Increasing value of τ3 enhances the stabilizing τ1-effect by delaying the loss of stability.


Figure 7. Delay effects on dynamics: (A) τ1-effect; (B) τ2-effect.

In the second simulation, we shift the emphasis away from the τ1-effect to the τ2-effect. Fixing τ3 = 7 and starting at point (2, 2) on the diagonal, we increase the value of τ2 from 2 to 7 along the vertical real line at τ1=τ10(=2) in Figure 6A. The line crosses the lower black island-shaped curve at τ2a(2.563), the black downward sloping curve at τ2b(4.356) and then upper black circle-shaped curve at τ2c(6.283). These threshold values are not labelled in Figure 6A by the same reason as before. The simulation results are given in Figure 7B in which stability gain and loss occur in the following way:

(i) we initially have a unstable steady state (in other words, a stable limit cycle) when the initial point is selected on the diagonal in Figure 6A11, that is inside the half black circle. The real part of one eigenvalue is positive.

(ii) the limit cycle becomes smaller as τ2 increases from 2 to τ2a and merges with the steady state for τ2=τ2a at which instability is switched to stability. The positive real part turns to be negative when the delay τ2 crosses the imaginary axis from left to right.

(iii) for τ2a<τ2<τ2b, the model has a stable steady state and stability loss occurs for τ2=τ2b. The real parts of at least one eigenvalue turn to be positive again when the delay τ2 crosses the imaginary axis.

(iv) for τ2>τ2b, Figure 7B suggests the existence of stable limit cycles and the instability of the steady state. Although further increasing τ2 intersects the island switching curve at point (2,τ2c) at which the real part of another eigenvalue changes sign. However no stability switch occurs, since there is already at least one other eigenvalue with positive real part.

Proposition 6. Increasing the value of τ3 does not affect essentially dynamics with respect to τ1 and can affect the dynamics with respect to τ2 as it can generate multiple stability switching.

4. Concluding Remarks

This paper constructs a heterogeneous agent model with three fixed delays and considers its dynamic behavior both analytically and numerically. The dependence of the delay effects on the changes in the lengths of the delays is also studied. The common result is that a longer delay can destabilize the market and can give rise to cyclic oscillations around the equilibrium (i.e., fundamental) price. This result clarifies the instability condition and thus complements the numerical study of Dibeh [9]. In addition, it is found that under multiple delays, stability loss and gain repeatedly occur as the length of a delay increases. He and Zheng [6] observed this phenomenon in a financial market model with continuously distributed time delay with finite memory length that involved infinitely many past price observations and called it the double edge effect. Our model is based on one, two or three past observations. It is thus shown in this paper that the same effect can occur even under finite delay information. In our future research the model of this paper will be extended to involve two risky assets and the stability and instability of the equilibrium will be examined in the extended models.

Author Contributions

AM: Performs simulations; FS: Constructs the model.

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 authors thank two referees for very helpful comments and highly appreciate the finanical supports from the MEXT-Supported Program for the Strategic Research Foundation at Private Universities 2013-2017, the Japan Society for the Promotion of Science (Grant-in-Aid for Scientific Research (C), 25380238, 26380316 and 16k03556) and Chuo University (Joint Research Grant). The usual disclaimers apply.


1. ^m = 0.4 is assumed and we will refer to this line later.

2. ^It is possible to have a switching surface in the 3D space. See Almodaresi and Bozorg [12] and Gu and Naghnaeian [13].

3. ^Solving the second equation of Equation (10) for τ3 gives the same result in a different form.

4. ^The remaining cases will be examined in a subsequent paper.

5. ^The segment Li(k, n) shifts upward if n increases and rightward if k increases. For (k, n)/2 and (k, n) = 0, the switching segments are defined but located ouside of the region with τi ≤ 5 for i = 1, 2, so they are not depicted in Figure 4A.

6. ^The τ1-value of the intersection is obtained in the following way: first solving τ2+(ω,1)=5 for ω yields a solution ωm ≃ 0.560 and then substituting it into τ1-(ω,1) gives τ1m0.887. So the black segment depicted in the yellow region is defined for ω ∈ [ωm, ω2].

7. ^It is possible to obtain this value by following the same procedure with which we obtained the value of τ1m just above.

8. ^Mathematically, it might be possible to regain stability for τ1 > 4. However, economically it is not the case as τ1 ≤ τ2 ≤ 5 is imposed by Assumptions 1 and 3.

9. ^This is not the case of the non-monotonics curve with α1 = 9/15.

10. ^τ1m” is not labeled in Figure 4A to avoid figure congestion.

11. ^More precisely, we select a constant initial function, g(t) = 2 for t ≤ 0.


1. Hommes C. Heterogeneous agent models in economics and finance. In: Tesfatsion L, Judd KL, editors. Handbook of Computational Economics, Agent-based Computational Economics, II. North-Holland: Amsterdam (2006). p. 1109–86.

Google Scholar

2. Chiarella C, He XZ, Hommes C. A dynamic analysis of moving average rules. J Econ Dyn Control (2006) 30:1729–53. doi: 10.1016/j.jedc.2005.08.014

CrossRef Full Text | Google Scholar

3. Zeeman E. On the unstable behavior of stock exchanges. J Math Econ. (1974) 1:39–49.

Google Scholar

4. Beja A, Goldman M. On the dynamic behavior of prices in disequilibrium. J Finance (1980) 35:235–48. doi: 10.2307/2327380

CrossRef Full Text | Google Scholar

5. Chiarella C. The dynamics of speculative behavior. Ann Oper Res. (1992) 37:101–23.

Google Scholar

6. He XZ, Zheng M. Dynamics of moving average rules in a continuous-time financial market model. J Econ Behav Organ. (2010) 76:615–34. doi: 10.1016/j.jebo.2010.08.005

CrossRef Full Text | Google Scholar

7. He XZ, Li K. Heterogeneous beliefs and adaptive behavior in a continuous-time asset price model. J Econ Dyn Control (2012) 36:973–87. doi: 10.1016/j.jedc.2012.02.002

CrossRef Full Text | Google Scholar

8. Xu X, Liu J, Guo L, Xu Z. Oscillatory dynamics in a continuous-time delay asset price model with dynamical fundamental price. Comput Econ. (2015) 45:517–29. doi: 10.1007/s10614-014-9428-9

CrossRef Full Text | Google Scholar

9. Dibeh G. Speculative dynamics in a time-delay model of asset prices. Physica A (2005) 355:199–208. doi: 10.1016/j.physa.2005.02.084

CrossRef Full Text | Google Scholar

10. Qu Y, Wei J. Global hopf bifurcation analysis for a time-delayed model of asset prices. Discrete Dyn Nat Soc. (2010) 2010:432821. doi: 10.1155/2010/432821

CrossRef Full Text | Google Scholar

11. Matsumoto A, Szidarovszky F. Heterogeneous Agent Model of Asset Price with Time Delays. IERCU Discussion Paper #247, Institute of Economic Research, Chuo University (2015).

12. Almodaresi A, Bozorg M. Stability crossing surfaces for linear time-delay systems with three delays. Int J Control (2009). 82:2304–10. doi: 10.1080/00207170903037739

CrossRef Full Text | Google Scholar

13. Gu K. Naghnaeian M. On stability crossing set for systems with three delays. IEEE Trans Automat Contr. (2011) 56:11–26. doi: 10.1109/TAC.2010.2050162

CrossRef Full Text | Google Scholar

14. Gu K, Niculescu S, Chen J. On stability crossing curves for general systems with two delays. J Math Anal Appl. (2005) 311:231–53. doi: 10.1016/j.jmaa.2005.02.034

CrossRef Full Text | Google Scholar

Keywords: heterogeneous agents model, three time delays, stability switching curves, delay effect, bifurcation

Citation: Matsumoto A and Szidarovszky F (2016) A Heterogeneous Agent Model of Asset Price with Three Time Delays. Front. Appl. Math. Stat. 2:15. doi: 10.3389/fams.2016.00015

Received: 24 June 2016; Accepted: 14 September 2016;
Published: 29 September 2016.

Edited by:

Iryna Sushko, Institute of Mathematics (NAS Ukraine), Ukraine

Reviewed by:

Davide Radi, University Carlo Cattaneo, Italy
Luca Guerrini, Marche Polytechnic University, Italy

Copyright © 2016 Matsumoto and Szidarovszky. 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) or licensor 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: Akio Matsumoto,