Damped jump-telegraph processes
We study a one-dimensional Markov modulated random walk with jumps. It is assumed that amplitudes of jumps as well as a chosen velocity regime are random and depend on a time spent by the process at a previous state of the underlying Markov process.
Equations for the distribution and equations for its moments are derived. We characterise the martingale distributions in terms of observable proportions between jump and velocity regimes.
keywords:inhomogeneous jump-telegraph process, Volterra equation, martingale measure
Msc:primary 60J27; secondary 60J75, 60K99
Telegraph processes with different switchings and velocity regimes are studied recently in connection with possibility of different applications such as, for instance, queuing theory (see Zacks (2004), Stadje and Zacks (2004)) and mathematical biology (see Hadeler (1999)). Special attention is devoted to financial applications (see Ratanov (2007), López and Ratanov (2012)). In the latter case, an arbitrage reasoning demands the presence of jumps.
The motions with deterministic jumps are studied in detail, see the formal expressions of the transition densities in Ratanov (2007), Di Crescenzo and Martinucci (2013). Such a model is developed for the option pricing problem, which is based on the risk-neutral approach, see Ratanov (2007). If the jump amplitudes are random, the case is less known. The telegraph processes of this type are studied earlier only under the assumption of mutual independence of jump values and jump amplitudes, see Stadje and Zacks (2004) and Di Crescenzo and Martinucci (2013). Similar setting were used for the purposes of financial applications, López and Ratanov (2012).
We present here a jump-telegraph process when an amplitude of the next jump depends on the (random) time spent by the process at the previous state. This approach is of special interest for the economical and the financial applications, everywhere when the comportment of process relates with friction and memory.
Assume that the particle moves with random (and variable) velocities performing jumps of random amplitude whenever the velocity is changed. More precisely, the actual velocity regime and the amplitude of the next jump are defined as (alternated) functions of the time spent by the particle at the previous state. We assume also that the time intervals between the subsequent state changes have sufficiently arbitrary alternated distributions. It creates an effect of damping process where a friction is generated by means of memory.
2 Generalised jump-telegraph processes: distribution
Let be a probability space. Consider two continuous-time Markov processes . The subscript indicates the initial state, (with probability 1). Assume that are left-continuos a. s.
Let be a Markov flow of switching times. The increments are independent and possess alternated distributions (with the distribution functions , , the survival functions and the densities ). We assume that , i. e. the state process is started at the switching instant. The distributions of and depend on the initial state . For brevity, we will not always indicate this dependence.
Consider a particle moving on with two alternated velocity regimes and . These velocities are described by two continuous functions . At each instant the particle takes the velocity regime , where is the (random) time spent by the particle at the previous state. We define a pair of the (generalised) telegraph processes driven by variable velocities as follows,
The integral is named the integrated telegraph process.
Let be a counting process. Notice that, and and .
The integrated telegraph process can be interpreted as the sum of random number of random variables. If , then
if , then the integrated telegraph process is expressed as
Notice that Simplifying we write instead of .
In the same manner we define the jump component. Let and be a pair of deterministic continuous (or, at least, boundary measurable) functions. Consider telegraph processes (2.1) based on instead of ,
An integrated jump process is defined in the form of compound Poisson process by the integral
The amplitude of a jump depends on the time spent by the particle at the current state.
It describes the particle which moves, alternating the velocity regimes at random times , starting from the origin at the velocity regime . Each velocity reversal is accompanied by jumps of random amplitudes, is the current particle’s position.
Conditioning on the first velocity reversal, notice that
Here denotes the equality in distribution. At each of two equalities the first term represents the movement without velocity reversal; the second one is the sum of three terms: the path till the first reversal, the jump value and the movement which is initiated after the first reversal.
The distribution of is separated into the singular and absolutely continuous parts.
The singular part of the distribution corresponds to the movement without any velocity reversals; let be the respective conditional distribution, if the initial state is fixed: for any Borel set we set
We denote the corresponding expectation by . On the space of (continuous) test-functions consider the linear functional (generalised function), It is easy to see that
The generalised function
can be viewed as the distribution “density”. Here is the Dirac measure (of unit mass) at point .
The absolutely continuous part of the distribution of is characterised by the densities
corresponds to the absolutely continuous part of distribution of .
Conditioning on the first velocity reversal, similarly to (2.6) we obtain the following equations,
(if the inner integrals are understood in the sense of the theory of generalised functions). Summing up in (2.8) we get the system of integral equations for (complete) distribution densities,
Here and are defined by (2.7).
Notice that , if (say, ). Define the functions , : for ,
Equations (2.8) have the following solution:
where the displacements are defined as the sum of alternating jumps, where if is odd, and if k is even.
Summing up we obtain the solution of (2.9):
where and are the modified Bessel functions.
3 Generalised jump-telegraph processes: moments
Therefore the expectations follow the equations of Volterra type:
Integrating by parts at the latter integral we have
which gives the following simplification for functions :
Here we denote .
Equations for variances can be derived similarly:
Let be a locally bounded measurable function. Assume that
Then the expectations
exist, and they satisfy the system
Functions satisfy the equations
where are generated by the preceding moments,
System (3.8) possesses a unique solution, see e.g. Linz (1985). Under appropriate assumptions the solution can be found explicitly. Consider the following example. Let the distributions of interarrival times are exponential:
In this particular case system (3.8) is solved by
where . Here we use the matrix notations ,
To check it, notice that system (3.8) is equivalent to ODE with zero initial condition:
where We get this equation by differentiating in (3.8) with subsequent integration by parts. Clearly, the unique solution is
Integrating by parts in (3.11) we obtain
Now, the desired representation (3.10) follows from
Let be the filtration, generated by .
The proof can be done by computing the conditional expectation for . Indeed,
According to the Markov property applied to the processes and we have
where and are copies of and respectively, independent of . Therefore,
Here denotes the integrated jump-telegraph process, which is initiated from the state , and is based on and . The latter expectation is equal to zero, , if and only if (4.1) holds.∎
Notice that if (4.1) holds, then the direction of jump should be opposite to the (mean) velocity value.
If the jump-telegraph processes and defined by (2.5) are martingales, then
Moreover, and are martingales, if and only if the distribution densities of interarrival times satisfy the following integral relations:
The latter equality is equivalent to (4.4).
Notice that by definition Hence, condition (4.3) is fulfilled. ∎
In this framework various particular cases of the martingale distributions and the corresponding distributions of interarrival times can be presented by applying Corollary 3. Consider the following examples.
Exponential distribution. Assume that functions and are proportional:
Relations (4.4) mean that the integrated jump-telegraph process is the martingale if the distributions of interarrival times are exponential:
Identities (4.6) can be written in detail as follows. The (observable) parameters of the model, i. e. the regimes of velocities and the regimes of jumps , satisfy the equations
with some positive constants and . These equations help to compute the switching intensities and by using the (observable) proportion between velocity and jump values. On the other hand, if mean velocity regimes are given, and , from these equations we can conclude that small jumps occur with high frequency, and big jumps are rare. The direction of jump should be opposite to the velocity sign, see also Remark 2.
In the framework of (2.5) we assume that the Markov flow of switching times has interarrival intervals which are exponentially distributed with alternated constant intensities . Let the velocity regimes and jump amplitudes are given, and they are proportional as in (4.6),
The martingale measure for exists and it is unique.
According to the Girsanov Theorem, see Ratanov (2007), we apply Radon-Nikodym derivative of the form