#### Abstract

We develop a continuous time model with heterogeneous fundamentalists, imitators, and discrete time delays. We show that nonlinear dynamic phenomena, such as coexistence of attractors and local and global bifurcations, occur due to the existence of a time gap in the process of adjustment of market prices.

#### 1. Introduction

Starting from the seminal paper by Brock and Hommes [1], attention has been placed in building on models that try to mimic the behaviours of financial markets also in the absence of stochastic shocks. A mechanism that has proved to be fruitful in this direction has been the introduction of heterogeneous expectations of financial markets operators [2]. This literature has essentially been developed through the study of discrete time models [3]. For instance, Chiarella et al. [4] deepen the mathematical properties of a two-dimensional dynamic model with fundamentalists and chartists, while Tramontana et al. [5] develop a three-dimensional nonlinear dynamic model of internationally connected financial markets.

By pursuing the line of research of He et al. [6], He and Zheng [7], and He and Li [8], the aim of this paper is to inquire whether complex dynamic phenomena obtained in discrete time models also hold in continuous time models with* discrete* time delays. To this purpose, we use a model similar to that adopted by He and Westerhoff [9] and Naimzada and Ricchiuti [10]. We assume two experts that have different beliefs on the fundamental of an asset and choose their allocations by using the mean-variance criterion in every moment in time [2]. Financial operators are imitators and select allocations established by the two experts depending on performances obtained.

We show that this can generate nonlinear dynamics also in a continuous time framework with time delays. In particular, coexistence of attractors as well as local and global bifurcations can occur.

The rest of the paper is organised as follows. Section 2 sets up the model. Section 3 characterises local stability properties and local bifurcations of equilibria. Section 4 provides some numerical experiments to validate the theoretical results established in Section 3, while also showing the emergence of global bifurcations. Section 5 outlines the conclusions.

#### 2. The Model

We consider a continuous time version of the model developed by Naimzada and Ricchiuti [10] augmented with* discrete* time delays. Specifically, we set up a model with a risk-free asset, characterised by a perfectly elastic supply and an instantaneous rate of return , and a risky asset, with a price per share and a (stochastic) dividend process . There are two types of market operators (fundamentalists and imitators) and a market maker that behaves as a Walrasian auctioneer. We assume there exist two types (Type and Type ) of fundamentalists with different beliefs on the fundamental of the traded assets. They are myopic and behave on the basis of the mean-variance criterion at each moment in time. By considering (1) in Naimzada and Ricchiuti [10, page 173], we get the following continuous time version of wealth dynamics of the Type fundamentalist:
where , , and are the wealth, its time derivative, and the share of risky asset of fundamentalist at time , respectively.

The objective of Type fundamentalist is therefore the following: where is a parameter that measures the degree of risk aversion of both agents. By assuming that the variance is constant and is given by , the maximisation programme (2) gives the following market demands of the risky asset:

Similarly with Naimzada and Ricchiuti [10], we assume that both fundamentalists have the same correct expectations on dividend dynamics, that is, , but heterogeneous expectations on price per share dynamics, that is, . This means that every fundamentalist expects that the value of the risky asset tends to a level believed as being its fundamental. From (3), it follows that the share of risky asset of fundamentalist is given by where and . In addition, we set without loss of generality.

In this model, fundamentalists play the role of experts in the market and other agents follow fundamentalists’ choices depending on a mechanism that rewards the portfolio allocations based on the fundamental closer to the realised market price. Specifically, we assume the following adjustment rule: where is the share of agents that follow Type fundamentalist, is the time derivative of , and is the squared error related to expert . Then, this rule stimulates agents to adopt portfolio choices of the fundamentalist whose fundamental value deviates less from the price actually realised on the market.

In addition, we assume the existence of a Walrasian auctioneer that fixes the price of the risky asset according to the following mechanism: where is the speed of adjustment of prices over time. The dynamics of price defined in (6) determines the variation of the market price of the risky asset depending on the price and allocation choices of imitators (followers) made at . This is a quite realistic assumption with regard to the adjustment rule of prices, given that the time gap can be referred to small time intervals (e.g., minutes or hours) as those observed in actual financial markets.

In order to simplify notation, in what follows we omit the time dependence for variables and derivatives referred at time and use to indicate the state of a generic variable at time . Now, by using (4) and (5) to substitute into (6) to eliminate and , respectively, we get where .

#### 3. Existence of Equilibria and Local Bifurcations

In this section we perform the analysis of the delay differential equation defined in (7) (see, e.g., [11]).

Equilibria of (7) are obtained by setting and for all . By doing this, one finds that the nontrivial equilibria of (7) must satisfy the following equation: Consequently, (7) has three positive equilibria: , , and .

Let . Then (7) can be transformed into The characteristic equation of the linear part of (9) is given by where When there is no delay, that is, in (9), the characteristic equation becomes . Then, and are locally asymptotically stable (), while is unstable ().

Now, assume that in (9). We will investigate the location of the roots of the transcendental equation. First, it is immediate that (10) has no zero root. Next, we examine when this equation has pure imaginary roots.

Proposition 1. *Characteristic equation (10) has a pair of purely imaginary conjugate roots when , where
*

* Proof. *For , is a root of (10) if and only if
Separating the real and imaginary parts, we obtain
If or , then (14) yield and ; on the other hand, if , then and . The conclusion is immediate.

From (12), we note that there exists an inverse relationship between the value of (and thus of the speed of adjustment of prices over time ) and the bifurcation value of . This implies that the speed of adjustment of prices over time plays a destabilising role in the model (i.e., for a high value of the Hopf bifurcations occur for a low value of ), while the variance and the degree of risk aversion play a stabilising role (i.e., for high values of and the Hopf bifurcations occur for a high value of ).

Lemma 2. * is a simple purely imaginary root of (10), and all other roots satisfy for any integer .*

*Proof. *If is not simple, that is, , then one has , which is a contradiction. Let us assume that there exists a root such that and for some . From (14), we get . Therefore, the statement follows immediately.

Let be the roots of (10) close to that satisfy and , where and are defined by (12). The next result indicates that each crossing of the real part of characteristic roots at must be from left to right; that is, stability is lost at the smallest stability switch and it cannot be regarded later.

Lemma 3. *The following transversality condition holds.*

* Proof. *Differentiating both sides of (10) with respect to gives
Then, we have
which implies
concluding the proof.

Lemma 4. *Let or . If , then all roots of (10) have negative real parts; if , then all roots of (10) except have negative real parts; if for , then (10) has roots with positive real parts.*

* Proof. *The first part follows noting that the equilibrium is locally asymptotically stable when and so its stability can only be lost if eigenvalues cross the imaginary axis from left to right. Let be a root of (10) with . Then , , and thus, we derive , which is a contradiction the left-hand side term of this equation being a negative number. Since the rate of change of the real part of an eigenvalue with respect to when is positive, then the number of roots with positive real parts is increasing. Due to the fact above, the number of roots of the characteristic equation with positive real part will be constant for and equal to the number of eigenvalues with positive real parts when (i.e., zero being stable). For each subsequent interval , the number can be determined from the number in the previous interval and the number of roots with zero real part at .

Spectral properties in the previous lemma immediately lead to stability properties of the zero solution of (9) and, equivalently, of the positive equilibrium of (7).

Theorem 5. *Let be defined as in (12).*(1)*The positive equilibrium is unstable for all .*(2)*The positive equilibria and are locally asymptotically stable for and unstable for . Furthermore, (7) undergoes a Hopf bifurcation at when .*

*Remark 6. *Since when (7) collapses into a one-dimensional autonomous differential equation, then in the continuous time version of the model by Naimzada and Ricchiuti [10] no persistent oscillations can occur, and the dynamics are therefore monotonic and convergent towards one of the equilibria. This result stresses the importance of time delays in generating endogenous fluctuations in this kind of models.

#### 4. Numerical Simulations

In this section we provide some numerical simulations to validate the theoretical results on local bifurcations stated in Section 3, and we show the occurrence of global bifurcations as well. For this purpose, we fix the following parameter set, , , and , and let vary. For this parameter values the Hopf bifurcation occurs at . For (but sufficiently high), the dynamics of the model are oscillatory and convergent towards one of the equilibria depending on the initial conditions (as shown in Figure 1). Just after the Hopf bifurcation value there exist two attracting closed invariant curves (each of which with its basin of attraction), and then the dynamics of prices (and also the shares of imitators of Type and Type fundamentalists) show persistent oscillations (see Figures 2(a) and 2(b)) characterised by unique maximum and minimum points.

**(a)**

**(b)**

**(a)**

**(b)**

By increasing the value of , other phenomena are possible. In particular, we observe an increase in the number of local maxima and minima that resembles the sequence of period-doubling bifurcations shown by Naimzada and Ricchiuti [10] in a discrete time model (see Figures 3(a) and 3(b)). When increases further, a global bifurcation occurs and the two attractors merge each other. Then, a unique attractor is born (as shown in Figures 4(a) and 4(b)). At this stage, the dynamics of prices are characterised by large oscillations and possibly a high degree of unpredictability as suggested by the bifurcation diagram shown in Figure 5.

**(a)**

**(b)**

**(a)**

**(b)**

#### 5. Conclusion

The debate on whether it is better to adopt continuous time models or discrete time models to describe the behaviours of financial markets operators is still open, as pointed out by He and Li [8, page 974]: “The [discrete time] set up facilitates economic understanding of the role of heterogeneous expectations and mathematical analysis, it, however, faces a limitation when dealing with expectations formed from the lagged prices over different time horizons and a challenge to characterise the adaptive behaviour in a continuous-time.” Generally speaking, in the absence of stochastic shocks the continuous time framework tends to generate regular dynamics. Nevertheless, when one assumes that some economic processes (e.g., the adjustment of prices) react to changes occurred at a certain time gap, the dynamics generated by continuous time models tend to mimic those generated by discrete time models. This paper has confirmed the existence of nonlinear asset price dynamics in a continuous time version (with* discrete* time delays) of the model proposed by Naimzada and Ricchiuti [10], characterised by the existence of fundamentalists with heterogeneous expectations on the value of a risky asset.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publishing of this paper.