Positivity of solution of nonhomogeneous stochastic differential equation with non-Lipschitz diffusion
We give a sufficient condition on coefficients of a nonhomogeneous stochastic differential equation with non-Lipschitz diffusion for a solution starting from arbitrary nonrandom positive point to stay positive. Some examples of application of the condition mentioned above are considered.
Збережено в:
Дата: | 2008 |
---|---|
Автори: | , |
Формат: | Стаття |
Мова: | English |
Опубліковано: |
Інститут математики НАН України
2008
|
Онлайн доступ: | http://dspace.nbuv.gov.ua/handle/123456789/4570 |
Теги: |
Додати тег
Немає тегів, Будьте першим, хто поставить тег для цього запису!
|
Назва журналу: | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
Цитувати: | Positivity of solution of nonhomogeneous stochastic differential equation with non-Lipschitz diffusion / Y. Mishura, S. Posashkova // Theory of Stochastic Processes. — 2008. — Т. 14 (30), № 3-4. — С. 77-88. — Бібліогр.: 6 назв.— англ. |
Репозитарії
Digital Library of Periodicals of National Academy of Sciences of Ukraineid |
irk-123456789-4570 |
---|---|
record_format |
dspace |
spelling |
irk-123456789-45702009-12-08T12:00:31Z Positivity of solution of nonhomogeneous stochastic differential equation with non-Lipschitz diffusion Mishura, Y. Posashkova, S. We give a sufficient condition on coefficients of a nonhomogeneous stochastic differential equation with non-Lipschitz diffusion for a solution starting from arbitrary nonrandom positive point to stay positive. Some examples of application of the condition mentioned above are considered. 2008 Article Positivity of solution of nonhomogeneous stochastic differential equation with non-Lipschitz diffusion / Y. Mishura, S. Posashkova // Theory of Stochastic Processes. — 2008. — Т. 14 (30), № 3-4. — С. 77-88. — Бібліогр.: 6 назв.— англ. 0321-3900 http://dspace.nbuv.gov.ua/handle/123456789/4570 en Інститут математики НАН України |
institution |
Digital Library of Periodicals of National Academy of Sciences of Ukraine |
collection |
DSpace DC |
language |
English |
description |
We give a sufficient condition on coefficients of a nonhomogeneous stochastic differential equation with non-Lipschitz diffusion for a solution starting from arbitrary nonrandom positive point to stay positive. Some examples of application of the condition mentioned above are considered. |
format |
Article |
author |
Mishura, Y. Posashkova, S. |
spellingShingle |
Mishura, Y. Posashkova, S. Positivity of solution of nonhomogeneous stochastic differential equation with non-Lipschitz diffusion |
author_facet |
Mishura, Y. Posashkova, S. |
author_sort |
Mishura, Y. |
title |
Positivity of solution of nonhomogeneous stochastic differential equation with non-Lipschitz diffusion |
title_short |
Positivity of solution of nonhomogeneous stochastic differential equation with non-Lipschitz diffusion |
title_full |
Positivity of solution of nonhomogeneous stochastic differential equation with non-Lipschitz diffusion |
title_fullStr |
Positivity of solution of nonhomogeneous stochastic differential equation with non-Lipschitz diffusion |
title_full_unstemmed |
Positivity of solution of nonhomogeneous stochastic differential equation with non-Lipschitz diffusion |
title_sort |
positivity of solution of nonhomogeneous stochastic differential equation with non-lipschitz diffusion |
publisher |
Інститут математики НАН України |
publishDate |
2008 |
url |
http://dspace.nbuv.gov.ua/handle/123456789/4570 |
citation_txt |
Positivity of solution of nonhomogeneous stochastic differential equation with non-Lipschitz diffusion / Y. Mishura, S. Posashkova // Theory of Stochastic Processes. — 2008. — Т. 14 (30), № 3-4. — С. 77-88. — Бібліогр.: 6 назв.— англ. |
work_keys_str_mv |
AT mishuray positivityofsolutionofnonhomogeneousstochasticdifferentialequationwithnonlipschitzdiffusion AT posashkovas positivityofsolutionofnonhomogeneousstochasticdifferentialequationwithnonlipschitzdiffusion |
first_indexed |
2025-07-02T07:46:53Z |
last_indexed |
2025-07-02T07:46:53Z |
_version_ |
1836520480399425536 |
fulltext |
Theory of Stochastic Processes
Vol.14 (30), no.3-4, 2008, pp.77-88
YULIYA MISHURA AND SVITLANA POSASHKOVA
POSITIVITY OF SOLUTION OF
NONHOMOGENEOUS STOCHASTIC
DIFFERENTIAL EQUATION WITH
NON-LIPSCHITZ DIFFUSION
We give a sufficient condition on coefficients of a nonhomogeneous
stochastic differential equation with non-Lipschitz diffusion for a so-
lution starting from arbitrary nonrandom positive point to stay posi-
tive. Some examples of application of the condition mentioned above
are considered.
1. Introduction
The main goal of this paper is to investigate the question about positivity
of solution of nonhomogeneous stochastic differential equation with non-
Lipschitz diffusion. Such stochastic differential equations arise in modelling
asset prices and interest rates on financial markets.
For example, Cox-Ingersoll-Ross interest rate model has the form:
rt = r0 +
∫ t
0
a(b− rs)ds+
∫ t
0
σ
√
rsdWs, t ≥ 0,
where r0, a, σ are real positive constants. It is easy to see that for rt > b the
drift is negative and for rt < b it is positive, so the solution of this equation
is mean-reverting.
Positivity of solutions of the stochastic differential equation with homo-
geneous coefficients of the form
X(t) = X0 +
∫ t
0
b(X(s))ds+
∫ t
0
σ(X(s))dW (s), t ≥ 0 (1)
was studied in [5], Chapter VI. This property is quite important for pro-
cesses, which model interest rate dynamics on financial market, because the
interest rate must be positive.
2000 Mathematics Subject Classifications. 60H10.
Key words and phrases. Cox-Ingersoll-Ross model, stochastic differential equation,
positivity of solutions.
77
78 YULIA MISHURA AND SVITLANA POSASHKOVA
The paper is organized as follows. Section 2 is devoted to nonhomoge-
neous stochastic differential equations with non-Lipschitz diffusion and con-
tains the sufficient condition of positivity of a solution starting from a non-
random positive value. In the paper [2] such sufficient condition was proved
in a particular case, where the diffusion was of the form σ(t, x) = σ(t)
√
x.
Section 3 contains some examples of application of the sufficient condi-
tion mentioned above, in particular we consider a nonhomogeneous version
of Cox-Ingersoll-Ross model.
2. Positivity of solution of nonhomogeneous stochastic
differential equation with non-Lipschitz diffusion
Consider a stochastic differential equation
X(t) = X0 +
∫ t
0
b(s,X(s))ds+
∫ t
0
σ(s,X(s))dW (s), t ≥ 0, (2)
where the initial value X0 > 0 is nonrandom, the coefficients b, σ : R+ ×
R → R are measurable, {W (t), t ≥ 0} is a Wiener process with respect to
filtration {Ft, t ≥ 0} on a probability space (Ω,F , P ).
Assume that the coefficients of this equation satisfy the following Ya-
mada conditions (Y1)–(Y4) (see e.g. [5,6]):
(Y1) the functions b and σ are jointly continuous.
(Y2) the coefficients grow at most linearly in x
|b(t, x)| + |σ(t, x)| ≤ C(1 + |x|), t ≥ 0, x ∈ R.
(Y3) the drift is Lipschitz continuous
|b(t, x) − b(t, y)| ≤ C |x− y| , t ≥ 0, x, y ∈ R.
(Y4) there exists such an increasing function ρ : R+ → R+, that∫
0+
ρ−2(u)du = +∞, and
|σ(t, x) − σ(t, y)| ≤ Cρ(|x− y|), t ≥ 0, x, y ∈ R.
Remark 2.1. An example of such function ρ is ρ(x) = xα, 1
2
≤ α ≤ 1, that
is why this condition is a kind of Hölder condition for σ. Other examples
are: ρ(x) = x1/2(ln 1/x)1/2, ρ(x) = x1/2(ln 1/x)1/2(ln ln 1/x)1/2 etc.
Definition 2.1. We say that the pair (X,W ) is a strong solution of the
equation (2), if X is the process adapted to the filtration
{F̄W
t
}
generated
by process W .
POSITIVITY OF SOLUTION 79
Definition 2.2. The equation (2) has the property of uniqueness of tra-
jectories if any two solutions (X,W ) and (X̃,W ), adapted to the same
filtration, satisfy the equality
P (∀t > 0 : X(t) = X̃(t)) = 1.
Theorem 2.1.([6]) Under conditions (Y1)–(Y4) the equation (2) has the
unique (relatively to its trajectories) solution from the class L2(Ω×[0, T ],F⊗
B([0, T ]), P × λ) for all T > 0, there λ is Lebesgue measure on R+.
Assume that for all x > 0 and t > 0 σ2(t, x) > 0, and that for all z > 0
the following inequality holds:
inf
x∈[z,+∞)
inf
t>0
σ2(t, x) = d = d(z) > 0. (3)
Let also for some ε > 0 there exist a positive continuous function A(t), t > 0,
such that for all t > 0 and x ∈ (0, ε)
2b(t, x)
σ2(t, x)
≥ A(t)
x
. (4)
The main result of the paper is the following one.
Theorem 2.2. Under the assumptions (Y1)–(Y4), (3)-(4) and the condi-
tion
∀t > 0 A(t) > 1 (equivalently ∀T > 0 inf
(0,T ]×(0,ε)
2b(t, x)x
σ2(t, x)
> 1)
the trajectories of the process {X(t), t ≥ 0} will be positive with probability 1.
Remark 2.2. For homogeneous Cox-Ingersoll-Ross model described by the
stochastic differential equation
drt = a(b− rt)dt+ σ
√
rtdWt
with a initial condition r0 > 0, the necessary and sufficient condition pro-
viding positivity of the solution is ab ≥ σ2
2
(see [5]), and in a nonhomo-
geneous case, if b(t, x) = a(t)(b(t) − x), σ(t, x) = σ(t)
√
x, then the coef-
ficient A(t) is A(t) = 2a(t)b(t)
σ2(t)
. So we obtain the condition 2a(t)b(t)
σ2(t)
> 1 or
a(t)b(t) > σ2(t)
2
, that is a generalization of the sufficient condition for the ho-
mogeneous case (a corresponding sufficient condition for nonhomogeneous
Cox-Ingersoll-Ross model also is proved in [2])
80 YULIA MISHURA AND SVITLANA POSASHKOVA
Proof. At first we prove that the solution of the equation (2) is a semimartin-
gale. Indeed, under condition of linear growth the process {Y (t), t ≥ 0} of
the form
Y (t) =
∫ t
0
σ(s,X(s))dW (s)
is a local martingale, because∫ t
0
σ2(s,Xs)ds ≤
∫ t
0
C2(1+Xs)
2ds ≤ 2C2t+2C2
∫ t
0
X2
sds < +∞ P −a.s.,
where the last inequality is a consequence of Theorem 2.1. So,
P{
∫ t
0
σ2(s, Ys)ds < +∞} = 1
means that {Y (t), t ≥ 0} is a local martingale. Now let us show that
Z(t) =
∫ t
0
b(s,X(s))ds, t ≥ 0 is the process of bounded variation. Indeed,
from the condition of linear growth (Y2) we can see that∫ t
0
|b(s,Xs)|ds ≤
∫ t
0
C(1 + |Xs|)ds = Ct+ C
∫ t
0
|Xs|ds < +∞,
where the last inequality follows by Theorem 2.1.
So, the process {X(t), t ≥ 0} is a sum of a local martingale and a process
of bounded variation, that is why it is a semimartingale,
It is sufficient to prove that the trajectories of the process X are positive
on any [0, T ]. So, let T > 0 be fixed.
Let ε > 0 be the constant from the condition (4). For x ∈ (ε,+∞) by
the condition (Y 2) there exists a constant K > 0 such that
|b(t, x)| ≤ K(1 + x).
We denote
p := K
(
1 +
1
ε
)
,
and then obtain the following inequality
∀x ∈ (ε,+∞) |b(t, x)| ≤ px.
We define
V (x) =
{
lnx, if 0 < x < ε,
1
ε
exp{−pε2
2d
} ∫ x
ε
exp{pu2
2d
}du+ ln ε, if ε ≤ x <∞.
It is easy to see that the function V (x) is continuously differentiable on
(0,+∞) and V ′(x) has the form
V ′(x) =
{
1
x
, if 0 < x < ε,
1
ε
exp{−pε2
2d
} exp{px2
2d
}, if ε ≤ x <∞.
POSITIVITY OF SOLUTION 81
We consider the differential operator
L =
1
2
σ2(t, x)
∂2
∂x2
+ b(t, x)
∂
∂x
.
Under condition (4) we have that for all x ∈ (0, ε) the following holds
LV =
1
2
σ2(t, x)(− 1
x2
) + b(t, x)
1
x
=
1
2x
σ2(t, x)(
2b(t, x)
σ2(t, x)
− 1
x
) ≥
≥ 1
2x
σ2(t, x)(
A(t)
x
− 1
x
) =
1
2x
σ2(t, x)
A(t) − 1
x
> 0.
For x ∈ (ε,+∞) the following inequality takes place
LV =
1
2
σ2(t, x) exp{−pε
2
2d
}px
dε
exp{px
2
2d
} + b(t, x)
1
ε
exp{−pε
2
2d
} exp{px
2
2d
} ≥
≥ d
px
dε
exp{−pε
2
2d
} exp{px
2
2d
} − px
1
ε
exp{−pε
2
2d
} exp{px
2
2d
} = 0.
Thus LV ≥ 0 for all x > 0, x
= ε.
Let 0 < m < X0 < M be fixed constants. We consider the random
variable
τm,M = inf{t : X(t) = m X(t) = M}
We use the following generalization of Ito’s formula from the article [3]. Let
X = X(t), t ≥ 0 be a continuous semimartingale and let f : R+ → R+ be a
continuous function of bounded variation.
We define lfs (X) as the local time of the process X on the curve f of the
form:
lfs (X) = P − lim
δ↓0
1
2δ
∫ s
0
I(f(r) − δ < Xr < f(r) + δ)d < X,X >r .
We make a remark that in our case the curve is given by the equation
f(s) = ε, because the second derivative of the function V has a break
only at this point. The process {X(t), t ≥ 0} is a semimartingale from
the proof above. That is why for any semimartingale, also for the process
{X(t), t ≥ 0}, the local time exists at any point, it follows from Theorem
5.52 on the page 186 in [1]. So, lfs (X) = lεs(X) exists.
We have as a consequence of Theorem 2.1 from the paper [3] that if we
set C = [0, ε] and D = [ε,+∞) (and taking into account that the function
V does not depend on t and V ′(x) is continuous on interval [0,+∞) in our
case), we obtain
V (X(τm,M ∧ T )) = V (X0) +
∫ τm,M∧T
0
V ′(X(s))dX(s)+
82 YULIA MISHURA AND SVITLANA POSASHKOVA
+
1
2
∫ τm,M∧T
0
V ′′(X(s))I(X(s)
= ε)d < X,X >s= V (X0)+
+
∫ τm,M∧T
0
V ′(X(s))b(s,X(s))ds+
∫ τm,M∧T
0
V ′(X(s))σ(s,X(s))dW (s)+
+
1
2
∫ τm,M∧T
0
V ′′(X(s))I(X(s)
= ε)σ2(s,X(s))ds,
where we have used the equality
d < X,X >s= σ2(s,X(s))ds.
Functions b(t, x) and σ(t, x) satisfy Yamada conditions. The function
b(t, x) is bounded for x ∈ [m,M ] and separated from 0 and from ∞; σ(t, ·)
is uniformly at Hölder continuous. Then it follows from Theorem 3.2.1 and
Corollary 3.2.2 in [4] that the transition probability function P (s, x; t,Γ) for
the process {X(t), t ≥ 0} exists and has the density
P (s, x; t,Γ) =
∫
Γ
p(s, x; t, y)dy,
where p(s, x; t, y), 0 ≤ s < t, x, y ∈ R is a positive function, which is
continuous in all variables. We set s := 0, x := X0 and obtain that the
distribution function of the process {X(t), t ≥ 0} has the form
P1(t,Γ) =
∫
Γ
p1(t, y)dy.
Thus, the process {X(t), t ≥ 0} has a density, so
P (X(t) = ε) = 0, ∀t > 0.
Further,
0 ≤ |E
∫ τm,M∧T
0
(V ′(X(s))b(s,X(s))I(X(s) = ε)ds| ≤
≤ |V ′(ε)| max
s∈[0,T ]
|b(s, ε)|
∫ T
0
P (X(s) = ε)ds = 0,
so,
E
∫ τm,M∧T
0
(V ′(X(s))b(s,X(s))I(X(s) = ε)ds = 0.
Take an expectation of both sides
V (X(τm,M ∧ T )) = V (X0) +
∫ τm,M∧T
0
V ′(X(s))b(s,X(s))ds+
POSITIVITY OF SOLUTION 83
+
∫ τm,M∧T
0
V ′(X(s))σ(s,X(s))dW (s)+
+
1
2
∫ τm,M∧T
0
V ′′(X(s))I(X(s)
= ε)σ2(s,X(s))ds,
and receive:
E(V (X(τm,M ∧ T ))) = V (X0) + E
∫ τm,M∧T
0
V ′(X(s))b(s,X(s))ds+
+E
∫ τm,M∧T
0
1
2
V ′′(X(s))I(X(s)
= ε)σ2(s,X(s))ds =
= V (X0) + E
∫ τm,M∧T
0
V ′(X(s))b(s,X(s))I(X(s) = ε)ds+
+E
∫ τm,M∧T
0
(V ′(X(s))b(s,X(s)) +
1
2
V ′′(X(s))σ2(s,X(s))I(X(s)
= ε)ds =
= V (X0) + E
∫ τm,M∧T
0
(V ′(X(s))b(s,X(s))+
+
1
2
V ′′(X(s))σ2(s,X(s))I(X(s)
= ε)ds.
The function V (x) is bounded from above on [0,M ], so we have
E(V (X(τm,M ∧ T ))) ≤ max
x∈[0,M ]
V (x)P (τ = τM ∧ T ) + V (m)P (τ = τm) ≤
≤ max
x∈[0,M ]
V (x) + V (m)P (τ = τm),
where τM (τm) is the first moment that the boundary m (boundary M) is
reached and τ = τm ∧ τM ∧ T .
Thus,
max
x∈[0,M ]
V (x) + V (m)P (τ = τm) ≥ V (X0)+
+E(
∫ τm,M∧T
0
(V ′(X(s))b(s,X(s)) +
1
2
V ′′(X(s))σ2(s,X(s)))I(X(s)
= ε)ds)
= V (X0) + E(
∫ τm,M∧T
0
LV (X(s))I(X(s)
= ε)ds) ≥ 0.
Let m tends to 0 and we obtain from above that
V (m) → −∞, m→ 0,
so the left-hand side of the inequality tends to −∞ and at the same time
the right-hand side of the inequality is nonnegative. Thus,
P (τ = τ0) = 0. (5)
84 YULIA MISHURA AND SVITLANA POSASHKOVA
So, for any fixed M and T the equality(5) holds. The proof follows when
M , T tend to +∞. �
4. Examples
Consider examples of application of Theorem 2.2 to some stochastic dif-
ferential equations.
Example 3.1.
Consider a stochastic differential equation
X(t) = X0 +
∫ t
0
(p(s) − q(s)X(s))ds+
∫ t
0
c(t)(X(s))αdW (s), t ≥ 0, (6)
where the initial value X0 > 0 is nonrandom, the functions p(t), q(t) and
c(t) are positive, continuous and inft>0 c(t) = d > 0, a constant 1
2
≤ α < 1.
Then in our notations
b(t, x) = p(t) − q(t)x,
σ2(t, x) = c2(t)x2α.
Thus, the left-hand part of inequality (4) can be rewritten in the form
2b(t, x)
σ2(t, x)
=
2(p(t) − q(t)x)
c2(t)x2α
=
2p(t)
c2(t)x2α
− 2q(t)
c2(t)
x1−2α ≥ 2p(t)
c2(t)x2α
Then, we put in (4) ε = 1, A(t) = 2p(t)
c2(t)
. For 1
2
≤ α < 1 we have that
1
x2α−1
A(t)
x
≥ A(t)
x
,
because 2α− 1 > 0 and x ∈ (0, ε) = (0, 1). So, by Theorem 2.2 we receive,
under the condition A(t) > 1 for all t > 0, which in our case has the form
∀t > 0 : p(t) >
c2(t)
2
,
that the trajectories of the process {X(t), t ≥ 0} are positive with proba-
bility 1.
Example 3.2.
Consider a stochastic differential equation
X(t) = X0 +
∫ t
0
(p(s) − q(s)X(s))ds+
POSITIVITY OF SOLUTION 85
+
∫ t
0
√
(X(s) ∨ 0)((c(t) −X(s)) ∨ 0)dW (s), t ≥ 0, (7)
where the initial value 0 < X0 < c(0) is nonrandom, the functions p(t), q(t)
and c(t) are positive, continuous and c(t) is also nondecreasing. Assume
that
p(t) < q(t)c(t) ∀t > 0.
Then we have in our notations
b(t, x) = p(t) − q(t)x,
σ(t, x) =
√
(x ∨ 0)((c(t) − x) ∨ 0).
Let us prove that such stochastic differential equation has a solution. It is
sufficient to verify the fulfillment of Yamada conditions on the coefficients.
We can consider the interval t ∈ [0, T ] and then tend T to +∞.
(Y1) From the initial condition on functions p(t), q(t) and c(t) we can
see that b(t, x) and σ(t, x) are jointly continuous.
(Y2)
|b(t, x)| + |σ(t, x)| ≤ p(t) + q(t)|x| + |x| + c(t) + |x| = p(t) + c(t)+
+|x|(1 + q(t)) ≤ max
t∈[0,T ]
(p(t) + c(t), 1 + q(t))(1 + |x|)
(Y3)
|b(t, x) − b(t, y)| = |q(t)(x− y)| ≤ max
t∈[0,T ]
q(t)|x− y|
(Y4)
|σ(t, x) − σ(t, y)| = |
√
x(c(t) − x) ∨ 0 −
√
y(c(t) − y) ∨ 0|
Let the function g : R+ × R × R → R be defined as
g(t, x, y) = |
√
x(c(t) − x) ∨ 0 −
√
y(c(t) − y) ∨ 0|.
We assume, that t ∈ [0, T ], and consider the next cases.
1) If (x < 0 or x > c(t)) and (y < 0 or y > c(t)) then
g(t, x, y) = 0 ≤
√
c(t)
√
|x− y|.
2) If 0 ≤ x ≤ c(t) and y < 0 then
g(t, x, y) = |
√
x(c(t) − x)| ≤ |
√
(x− y)(c(t) − x)| ≤ |
√
(x− y)c(t)| =
=
√
c(t)
√
|x− y|.
86 YULIA MISHURA AND SVITLANA POSASHKOVA
Note that under the symmetry the same estimation for g(t, x, y) holds when
x < 0 and 0 ≤ y ≤ c(t).
3) If 0 ≤ x ≤ c(t) and y > c(t) then
g(t, x, y) = |
√
x(c(t) − x)| ≤ |
√
x(y − x)| ≤ |
√
c(t)(y − x)| =
=
√
c(t)
√
|x− y|.
Note that under the symmetry the same estimation for g(t, x, y) holds when
x > c(t) and 0 ≤ y ≤ c(t).
4) If 0 ≤ x ≤ c(t) and 0 ≤ y ≤ c(t) then
g(t, x, y) = |
√
x(c(t) − x)−
√
y(c(t) − y)| ≤
√
|x(c(t) − x) − y(c(t) − y)| =√
|x(c(t) − x) − y(c(t) − x) + y(c(t) − x) − y(c(t) − y)| =
=
√
|(x− y)(c(t) − x) + y(c(t) − x− c(t) + y)| =
=
√
|(x− y)(c(t) − x− y)| ≤
√
c(t)|x− y| =
√
c(t)
√
|x− y|.
Thus, in all cases we have that
g(t, x, y) ≤
√
c(t)
√
|x− y|.
That is why
max
t∈[0,T ]
g(t, x, y) ≤ ( max
t∈[0,T ]
√
c(t))
√
|x− y|.
We obtain that
|σ(t, x) − σ(t, y)| = |
√
x(c(t) − x) ∨ 0 −
√
y(c(t) − y) ∨ 0| ≤
≤ ( max
t∈[0,T ]
√
c(t))
√
|x− y|.
Then ρ(|x− y|) =
√|x− y|, ρ(x) =
√
x increases and∫
0+
ρ−2(u)du =
∫
0+
1
u
du = +∞.
It means that Yamada’s conditions are fulfilled and the equation has the
unique solution.
Let us prove that X(t) ≥ 0 a.s. for all t ≥ 0. Define τ1 = inf{t : X(t) =
−δ}, where δ > 0 is some fixed constant. Assume that P{τ1 < ∞}. So, a
constant r < τ1 exists such that
X(t) < 0 t ∈ (r, τ1) a.s.
However, in this case
dX(t) = (p(t) − q(t)X(t))dt > 0.
POSITIVITY OF SOLUTION 87
So, the function t→ X(t) increases on (r, τ1) which is impossible.
Now we prove that X(t) ≤ c(t) a.s. for all t ≥ 0. Define τ2 = inf{t :
c(t) − X(t) = −δ}, where δ > 0 is some fixed constant. Assume that
P{τ2 <∞}. So, a constant R < τ2 exists such that
−δ < c(t) −X(t) < −δ + θ t ∈ (R, τ2) a.s.,
where θ is such constant that p(t)−q(t)c(t)
q(t)
+ θ < 0, t ∈ (R, τ2). Then
c(t) + δ − θ < X(t) < c(t) + δ t ∈ (R, τ2).
But
dX(t) = (p(t) − q(t)X(t))dt < (p(t) − q(t)(c(t) + δ − θ))dt =
= ((p(t) − q(t)c(t) + q(t)θ) − q(t)δ)dt < 0.
We obtain that dX(t) < 0, t ∈ (R, τ2), thus, the function t → X(t)
decreases on (r, τ1), but this is impossible (because the function c(t) − x(t)
must decrease and the function c(t) is nondecreasing). From this we obtain
that the initial stochastic differential equation can be rewritten in the form
X(t) = X0+
∫ t
0
(p(s)−q(s)X(s))ds+
∫ t
0
√
X(s)(c(t) −X(s))dW (s), t ≥ 0.
We know from above that b(t, x) = p(t) − q(t)x, σ2(t, x) = x(c(t) − x)
Thus, the inequality (4) can be rewritten in the form
2b(t, x)
σ2(t, x)
=
2(p(t) − q(t)x)
x(c(t) − x)
=
2p(t)
x(c(t) − x)
− 2q(t)
c(t) − x
.
As x ∈ (0, ε), where ε < c(0), then we can estimate the above expression as
2p(t)
x(c(t) − x)
− 2q(t)
c(t) − x
≥ 2p(t)
xc(t)
− 2q(t)
c(t) − ε
≥ 2p(t)
xc(t)
Denote A(t) = 2p(t)
c(t)
. Thus, by Theorem 2.2 the condition A(t) > 1 must
be satisfied for all t > 0, and in our case it has the form
∀t > 0 : p(t) >
c(t)
2
.
It is sufficient for trajectories of the process {X(t), t ≥ 0} be positive a.s.
So, if
c(t)
2
< p(t) < q(t)c(t) ∀t > 0,
then trajectories of the process {X(t), t ≥ 0} given by the stochastic differ-
ential equation (7) are positive with probability 1.
88 YULIA MISHURA AND SVITLANA POSASHKOVA
4. Conclusion
We declare the sufficient condition on coefficients which provides a.s.
positivity of the trajectories of the solution of the stochastic differential
equation with nonhomogeneous coefficients and non-Lipschitz diffusion. The
result of this paper is applied to some stochastic differential equations, in
particular for nonhomogeneous Cox-Ingersoll-Ross model.
References
1. Jacod J., Calcul Stochastique Et Problemes De Martingales, Berlin: Sprin-
ger, (1979), 539 p. (French)
2. Mishura Y. S., Posashkova S. V., Shevchenko G. M., Properties of solu-
tions of stochastic differential equation with nonhomogeneous coefficients
and non-Lipschitz diffusion, Theory of Probability and Math.Statistics,
79, (2008), submitted. (Ukrainian)
3. Peskir G., A Change-of-Variable Formula with Local Time on Curves, J.
Theoretical Probability, 3, (2005), 499–535.
4. Stroock D. W., Varadhan S. R.§., Multidimensional Diffusion Processes,
Berlin Heidelberg New York Springer-Verlag, (1979), 340 p.
5. Watanabe S., Ikeda N., Stochastic differential equations and diffusion pro-
cesses, M.: Science, (1986), 448 p.
6. Yamada T., Sur l’approximation des solutions d’équations différentielles
stochastiques, Z. Wahr. Verw. Geb., 36, (1976), 153–164.
Department of Probability Theory and Mathematical Statistics,
Kyiv National Taras Shevchenko University, Kyiv, Ukraine
E-mail address: myus@univ.kiev.ua, revan1988@gmail.com
|