1. Introduction
Diffusion processes are used as models in various applications, in particular in neuroscience to emulate the dynamics of the membrane potential of a neuron [
1]. Moreover, to take into account the spikes of the neuron, jump-diffusion processes have been proposed by Jahn et al. [
2] and Melanson and Longtin [
3], among others.
Now, diffusion and jump-diffusion processes both increase and decrease in any interval, however small. However, in some applications, it is not realistic to assume that the process can decrease or increase. For example, if represents the wear of a machine at time t, the stochastic process should increase with time.
One way to obtain a strictly increasing or decreasing process is to consider degenerate two-dimensional diffusion processes
defined by
where
is either a strictly positive or negative function and
is a standard Brownian motion. The functions
f and
v are such that
is a diffusion process. When
, the process
is called an integrated diffusion process. We can of course generalize the definition to the case when
is a jump-diffusion process.
The author has published a number of papers on integrated diffusion processes; see, for instance, Lefebvre [
4] for a recent one. Other papers on this topic include those by Lachal [
5], Makasu [
6], Metzler [
7], Caravelli et al. [
8] and Levy [
9].
Next, we define the first-passage time
where
C is a subset of
.
Let
be the joint probability density function (pdf) of the random vector
, with
. As is well known (see, for example, Cox and Miller [
10], p. 247), the function
satisfies the Kolmogorov backward equation
Moreover, the pdf
of the random variable
satisfies the same partial differential equation (PDE):
(subject to different boundary conditions). It follows that the moment-generating function of the random variable
, namely
where
, is a solution of the following PDE:
where
, etc. Furthermore, this equation is subject to the boundary condition
We now replace the diffusion process
by the jump-diffusion process defined by
where
is a Poisson process with rate
. The random variables
,
,
… are assumed to be independent and identically distributed (i.i.d.), and also independent of the Poisson process. We can state the following proposition.
Proposition 1. The function satisfies the integro-differential equation (dropping the dependence on α from the notation)for , where is the common density function of the s. As above, the equation is subject to the boundary condition (8). Proof. This result is obtained by generalizing the
infinitesimal generator of the jump-diffusion process
in Kou and Wang [
11] to the case when
and
are not necessarily constant functions. □
Remark 1. See also the remark after the proof of Proposition 2.
There are still few explicit solutions to first-passage problems for jump-diffusion processes. Kou and Wang [
11] obtained explicit formulae for the Laplace transform of the pdf of the first-passage time
to a constant boundary for a Wiener process with jumps having a double exponential distribution. This model was generalized or modified in various papers. In Chen et al. [
12],
was the first-exit time from a finite interval, while in Yin et al. [
13] the jumps were mixed-exponential random variables. Karnaukh [
14] considered the case when the parameters of the Wiener process depend on a finite Markov chain. In Lefebvre [
15], the jump sizes were assumed to be uniformly distributed, while in Abundo [
16] the jumps (positive and/or negative) were of a constant size and the boundaries were time-dependent. Because obtaining exact analytical solutions to such problems is difficult, some authors presented numerical techniques to obtain the quantities of interest; see, for example, Belkaid and Utzet [
17].
Di Crescenzo et al. [
18] computed bounds for first-crossing-time probabilities of a Wiener process with random jumps driven by a counting process. Fernández et al. [
19] proposed algorithms to compute double-barrier first-passage-time probabilities of a jump-diffusion process with an arbitrary jump size distribution.
D’Onofrio and Lanteri [
20] obtained numerical approximations for the density functions of first-passage times in the case of diffusion processes with state-dependent jumps. Finally, in Lefebvre [
21], the author was able to obtain exact solutions to optimal control problems for Wiener processes with exponential jumps.
In the current paper, explicit results will be presented for the first-passage time of a two-dimensional jump-diffusion process.
In the next section, the special case when the two-dimensional process
is
killed at the latest at the moment of the first jump, will be considered. We are also interested in the mean value of
, as well as in the probability that
will leave the continuation region through a given part of its boundary
. In
Section 3, the problem of maximizing or minimizing the time the controlled version of the process
spends in the continuation region
C will be treated. A particular problem will be solved explicitly in
Section 4. Finally, we will conclude with a few remarks in
Section 5.
2. Killed Processes
Assume that the random variables
are such that no overshoot is possible. That is, the degenerate two-dimensional jump-diffusion process
cannot jump over the boundary of the continuation region
C. Let
and
where
. We have the following corollaries.
Corollary 1. The function satisfies the integro-differential equation for , subject to the boundary condition Proof. It follows from the expansion of
into an infinite series (see Cox and Miller [
10]):
Notice that in our case,
will exist for any
because of Equation (18) below. □
Corollary 2. The probability is a solution of the integro-differential equation for . Moreover, the boundary condition is where .
In this paper, we consider the special case when the random variable
is such that the process
will leave the continuation region
C at the latest at the moment
of the first event of the Poisson process. Let
be the random variable that corresponds to
when
. We can write that
where
has an exponential distribution with parameter
. Furthermore, the sum in Equation (9) can be replaced by
, where
is the indicator variable of the event
, and the equation is valid for
. We can say that the process
is
killed at time
.
An application of the above problem is the following: as mentioned in
Section 1, a more realistic model for the wear of a machine is the degenerate two-dimensional process
defined in Equations (
1) and (2), when
is a strictly increasing function. Rishel [
22] proposed this model (in
n dimensions) in the context of an optimal control problem. If
denotes the remaining amount of deterioration that a device can undergo before it must be repaired or replaced, then
should be strictly negative instead. Moreover, the remaining lifetime of the device is the first-passage time to zero or to a level at which it is considered worn out.
Now, many electronic devices, especially mobile phones, are often replaced as soon as they break down, rather than being repaired. A mobile phone failure can be seen as a jump from the current value of to zero, so that the device is killed at the time of the jump. It is also possible that the device will be replaced before a failure occurs, due to normal wear and tear or because it has become obsolete. Thus, deterioration could also include the age of the device.
Similarly, in the case of humans, the downward jump to zero could occur during a massive heart attack or stroke.
Because we assume that
for any possible value
z of the random variable
Z, the integro-differential Equations (
10) and (
12) become, respectively, the partial differential equations
and
In the case of Equation (16), if
, then
whereas, we have
if
. If
belongs to
for some values of
z, and to
for other values of
z, then the integral in Equation (16) is replaced by
Solving integro-differential equations explicitly and exactly is not an easy task. In
Section 4, an example of a problem that we can indeed solve analytically will be presented. As above, the integro-differential equations will be reduced to PDE’s, and the method of similarity solutions will be used to transform these PDE’s into ordinary differential equations.
3. Optimal Control
In this section, we consider a controlled version of the two-dimensional process
:
where
is the control variable, which is assumed to be a continuous function, and
b is a non-zero constant. The aim is to find the value of the control that minimizes the expected value of the cost function
where
and
are constants. If the parameter
is positive (respectively negative), the optimizer must try to minimize (respectively maximize) the time spent by the controlled process in the continuation
C, taking the quadratic control costs into account. This type of problem is known as a
homing problem; see Whittle [
23] and/or [
24].
To solve the above problem, we can make use of dynamic programming. We define the
value function
That is,
is the expected cost (or reward, if it is negative) obtained by choosing the optimal value of the control variable in the interval
.
Proposition 2. The value function satisfies the second-order, non-linear partial integro-differential equationMoreover, we have the boundary condition Proof. First, thanks to Bellman’s principle of optimality, we can write that
Next, let
and
We have
Since
has a Poisson distribution with parameter
, we can write that
and
Hence,
Now, assuming that
is twice differentiable with respect to
x and to
y, making use of Taylor’s formula for functions of two variables, we obtain that
Furthermore, we have
and
, which implies that
Similarly, we find that
Indeed, by independence, we have
Let
. We compute
so that
Thus,
From Equation (31) and the above results, we deduce that
Dividing both sides of the above equation by
and letting
decrease to zero, we obtain the
dynamic programming equationFrom the preceding equation, we find that the optimal control
can be expressed as follows:
Substituting the optimal control into Equation (46), we obtain Equation (28).
Finally, the boundary condition (
29) follows at once from the definition of
, since
if
. □
Remark 2. Suppose that we set equal to zero in Equation (25) and that we replace the cost function defined in Equation (26) bywhere is the pdf of when the starting time is equal to . Then, since is actually a deterministic function, we may write thatProceeding as in the above proof, we find thatWe have for . Moreover, using the Leibniz integral rule,where we used the fact that because the two-dimensional process is time-invariant. Hence, setting equal to zero, we retrieve Equation (10). In the case of the killed processes considered in
Section 2, the integro-differential Equation (28) reduces to the non-linear PDE
The boundary condition remains the one in Equation (29).
In the next section, a particular problem will be treated. We will find the exact optimal control when the parameter is equal to zero, so that there are no jumps, and a numerical solution will be computed in the case when .
4. A Particular Problem
We consider the process
, defined by
That is,
is a standard Brownian motion with jumps. Moreover, we can write that
which implies that
Let
where
. Notice that
in the continuation region, so that
is strictly increasing with time.
Next, we define
Thus,
is a discrete random variable such that at time
the process will leave the continuation region, if it has not already done so. We can write that
where
is the Dirac delta function.
We deduce from Equation (19) that the moment-generating function of
satisfies the PDE
Based on this equation and the boundary conditions
if
or
, we look for a solution of the form
where
. This is an application of the method of similarity solutions, and
r is called the
similarity variable. For the method to work, both the equation and the boundary conditions must be expressed in terms of
r (after simplification). Here, we find that Equation (60) reduces to the
ordinary differential equation (ODE)
subject to the boundary conditions
, for
. With the help of the mathematical software program
Maple, we find that the general solution of the above equation can be written as
where
and
are Kummer functions. The constants
and
are uniquely determined from the boundary conditions
.
Since, as noted in
Section 2 (see Equation (18)),
, when
is large, the function
should be close to the moment-generating function of an exponential random variable with parameter
, namely
In
Figure 1, we present the functions
and
for
, when
,
,
and
. We observe that the two functions differ significantly. However, the two functions are very similar when
, as can be observed in
Figure 2. When
,
and
practically coincide for
.
Next, the function
satisfies the PDE (see Equation (20))
subject to
if
or
. Setting
, we obtain the ODE
with
, for
. We find that
The particular solution that satisfies the boundary conditions
is presented in
Figure 3.
Finally, let
This function is a solution of the PDE
Assuming that
, we obtain the ODE
whose general solution is
The solution that satisfies the boundary conditions
and
is shown in
Figure 4, when
and
.
To conclude this section, we will try to find the optimal control of the two-dimensional process
defined by
To do so, we must solve the non-linear second-order PDE
subject to
if
or
.
As above, we make use of the method of similarity solutions. We look for a solution of the form
. Equation (74) becomes
If
,
Maple is able to obtain the general solution of the preceding equation:
where
and
The constants
and
are determined by making use of the boundary conditions
.
When
,
Maple and
Mathematica are unable to provide an analytical expression for the solution of Equation (75). It is, however, not difficult to obtain a numerical solution for any choice of the parameters. For instance, if we choose
,
and
, we obtain the function
, as shown in
Figure 5, together with the function obtained when
. Finally, in
Figure 6, we present the corresponding optimal controls.
5. Conclusions
In this paper, we have considered degenerate two-dimensional jump-diffusion processes, defined in such a way that the first component of the vector is a strictly increasing or decreasing function with respect to time. This kind of process is more realistic than a one-dimensional diffusion or jump-diffusion process in many applications, especially when represents the age or wear of a certain device. We could generalize the model by incorporating more than one diffusion process . The diffusion processes could model the various variables that influence . For example, in the case of wear, important factors to consider are temperature, speed of use, etc.
In
Section 2, we obtained equations for functions defined in terms of a first-passage time for the processes
. Moreover, we treated an optimal control problem for these processes in
Section 3. Finally, a particular problem was solved explicitly in
Section 4.
As mentioned in
Section 1, there are few first-passage problems for one-dimensional jump-diffusion processes that have been solved exactly and explicitly so far. Here, we were able to find exact analytical expressions for quantities defined in terms of a first-passage time for a (degenerate) two-dimensional jump-diffusion process. Furthermore, in
Section 2, we saw that the processes considered in this paper could serve as models in real-life applications, such as the remaining amount of deterioration that a given device can undergo before it needs to be repaired or replaced.
In general, to solve this type of problem, it is necessary to find the solution of an integro-differential equation with partial derivatives. We considered the case when the process leaves the continuation region at the latest when the first event of the Poisson process occurs. In this case, the equation to be solved becomes a partial differential equation. Using the method of similarity solutions, it is sometimes possible to reduce this PDE to an ODE. It should also be possible to find a numerical solution to any particular problem.
As a follow-up to this work, we would like to find exact analytical solutions to problems where the equations to be solved are integro-differential equations; for example, by trying to transform the integro-differential equations into differential equations.