Skip to content
BY 4.0 license Open Access Published by De Gruyter Open Access July 9, 2019

Razumikhin-type theorem on time-changed stochastic functional differential equations with Markovian switching

  • Xiaozhi Zhang EMAIL logo and Chenggui Yuan
From the journal Open Mathematics

Abstract

This work is mainly concerned with the exponential stability of time-changed stochastic functional differential equations with Markovian switching. By expanding the time-changed Itô formula and the Razumikhin theorem, we obtain the exponential stability results for the time-changed stochastic functional differential equations with Markovian switching. What’s more, we get many useful stability results by applying our new results to several important types of functional differential equations. Finally, an example is given to demonstrate the effectiveness of the main results.

MSC 2010: 34D20; 34K50

1 Introduction

The research for stochastic differential equations (SDEs) is a mature field, which plays an important role in modeling dynamic system considering uncertainty noise in many applied areas such as economics and finance, physics, engineering and so on. Many qualitative properties of the solution of stochastic functional differential equations (SFDEs) have been received much attention. In particular, the stability or asymptotic stability of SFDEs has been studied widely by more and more researchers ([1, 2, 3, 4, 5]).

Recently, Chlebak et al.[6] discussed sub-diffusion process and its associated fractional Fokker-Planck-Kolmogorov equations. The fractional partial differential equations are well known to be connected with limit process arising from continuous-time random walks. The limit process is time-changed Lev́y process, which is the first hitting time process of a stable subordinator (see [7, 8, 9] for details). The existence and stability of SDE with respect to time-changed Brownian motion recently have received much attention([10, 11]). Wu [12, 13] established the time-changed Itô formula of time-changed SDE, and then obtained the stability results. Subsequently, Nane and Ni [14] established the Itô formula for time-changed Lévy noise, then discussed the stability of the solution.

However, to the best of our knowledge, there are no results for the time-changed stochastic functional differential equations with Markovian switching published till now. Motivated strongly by the above, in this paper, we will study the stability of time-changed SFDEs with Markovian switching. By applying the time-changed Itô formula and Lyapunov function, we present the Razumikhin-type theorem([15, 16]) of the time-changed SFDEs with Markovian switching. More precisely, we consider the following SFDEs with Markovian switching driven by time-changed Brownian motions:

dx(t)=h(xt,t,Et,r(t))dt+f(xt,t,Et,r(t))dEt+g(xt,t,Et,r(t))dBEt (1.1)

on t ≥ 0 with {x(θ) : − τθ ≤ 0} = ξ CF0b ([−τ, 0];ℝn), where h, f, g are appropriately specified later.

In the remaining parts of this paper, further needed concepts and related background will be presented in Section 2. In Section 3, the exponential stability results of the time-changed SFDEs with Markovian switching will be given. Many useful types of results of stochastic delay differential equations and stochastic differential equations are presented in Section 4 and Section 5 respectively. Finally, an example is given to show the availability of the main results.

2 Preliminary

Throughout this paper, let (Ω, 𝔉, {𝔉}t≥0, P) be a complete probability space with the filtration {𝔉}t≥0 which satisfies the usual condition(i.e. {𝔉}t≥0 is right continuous and 𝔉 contains all the P-null sets in 𝔉). Let {U(t), t ≥ 0} be a right continuous with left limit (RCLL) increasing Lévy process that is called subordinator starting from 0. For a subordinator U(t), in particular, is a β -stable subordinator if it is a strictly increasing process denoted by Uβ(t) and characterized by Laplace transform

E[exp(sUβ(t))]=exp(tsβ),s>0,β(0,1).

For an adapted β -stable subordinator Uβ(t), define its generalized inverse as

Et:=Etβ=inf{s>0:Uβ(s)>t},

which means the first hitting time process. And Et is continuous since Uβ(t) is strictly increasing.

Let Bt be a standard Brownian motion independent on Et, define the following filtration as

Ft=s>tσ[Br:0rs]σ[Er:r0],

where σ1σ2 denotes the σ -algebra generated by the union of σ -algebras σ1 and σ2. It concludes that the time-changed Brownian motion BEt is a square integrable martingale with respect to the filtration {𝔉Et}t≥0. And its quadratic variation satisfies < BEt, BEt > = Et.([17])

Let r(t), t ≥ 0 be a right continuous Markov chain on the probability space taking values in a finite state space S = {1, 2, …, N} with generator Γ = (γij)N×N by

P{r(t+Δ)=j|r(t)=i}=rijΔ+o(Δ)ifij,1+rijΔ+o(Δ)ifi=j,

where Δ > 0, γij is the transition rate from i to j if ij and γii = ijγij . We assume that the Markov chain r(t) is independent on Brownian motion, it is well known that almost each sample path of r(t) is a right-continuous step function.

For the future use, we formulate the following generalized time-changed Itô formula.

Lemma 2.1

(The generalized time-changed Itô formula) Suppose Uβ(t) is a β -stable subordinator and Et is its associated inverse stable subordinator. Let x(t) be a 𝔉Et adapted process defined in (1.1). If V : ℝn × ℝ+ × ℝ+ × S → ℝ is a C2,1,1 (ℝn × ℝ+ × ℝ+ × S; ℝ) function, let

L1V(xt,t,Et,i)=Vt(x,t,Et,i)+Vx(x,t,Et,i)h(xt,t,Et,i)+j=1NγijV(x,t,Et,j)

and

L2V(xt,t,Et,i)=VEt(x,t,Et,i)+Vx(x,t,Et,i)f(xt,t,Et,i)+12trace[gTVxxg(xt,t,Et,i)],

then with probability one

V(x(t),t,Et,r(t))=V(x0,0,0,r(0))+0tL1V(xs,s,Es,r(s))ds+0tL2V(xs,s,Es,r(s))dEs+0tVx(x(s),s,Es,r(s))g(xs,s,Es,r(s))dBEs+0tR[V(x(s),s,Es,i0+h(r(s),l))V(x(s),s,Es,r(s))]μ(ds,dl),

where μ(ds, dl) = ν(ds, dl) − m(dl)ds is a martingale measure, ν(ds, dl) is a Poisson random measure with density dt × m(dl), in which m is the Lebesgue measure on ℝ.

Proof

Let y = [x, t1, t2]T = [x, t, Et]T, and G(y(t), r(t)) = V(xt, t, Et, r(t)). Based on the computation rules ([8]), we have

dtdt=dEtdEt=dtdEt=dtdBEt=dEtdBEt=0,dBEtdBEt=dEt.

Applying the multi-dimensional Itô formula([18]) to G(y(t), r(t)) yields that

G(y(t),r(t))=G(y(0),r(0))+0tGy(y(s),r(s))dy(s)+0t12dyTGyydy+0tj=1NγijG(y(s),j)ds+0tR[G(y(s),i0+h(r(s),l),x(s))G(y(s),r(s))]μ(ds,dl)=G(y(0),r(0))+0T[VxVt1Vt2]hdt+fdEt+gdBEtdt1dt2+0t12trace[gTVxxg]dEt+0tj=1NγijV(x(s),s,Es,j)ds+0tR[V(x(s),s,Es,i0+h(r(s),l))V(x(s),s,Es,r(s))]μ(ds,dl)=V(x0,0,0,r(0))+0tVx(x(s),s,Es,r(s))g(xs,s,Es,r(s))dBEs+0tVEs(x(s),s,Es,r(s))+Vxf(xs,s,Es,r(s))+12trace(gTVxxg)dEs+0tVt(x(s),s,Es,r(s))+Vxh(xs,s,Es,r(s))+j=1NγijV(x(s),s,Es,j)ds+0tR[V(x(s),s,Es,i0+h(r(s),l))V(x(s),s,Es,r(s))]μ(ds,dl).

This completes the proof. □

Corollary 2.1

Suppose Uβ(t) is a β -stable subordinator and Et is its associated inverse. Let x(t) be an 𝔉Et adapted process defined in (1.1). If V:ℝn × ℝ+ × ℝ+ × S → ℝ is a C2,1,1(ℝn × ℝ+ × ℝ+ × S; ℝ) function, then for any stopping time 0 ≤ t1t2 < ∞

EV(x(t2),t2,Et2,r(t2))=EV(x(t1),t1,Et1,r(t1))+Et1t2L1V(xs,s,Es,r(s))ds+Et1t2L2V(xs,s,Es,r(s))dEs

where L1 and L2 are defined in the lemma above.

In this paper, the following hypothesis is imposed on the coefficients h, f and g.

  1. Both h, f : ℝn × ℝ+ × ℝ+ × S → ℝn and g : ℝn × ℝ+ × ℝ+ × S → ℝn×m are Borel-measurable functions. They satisfy the Lipschitz condition. That is, there is L > 0 such that

    |h(ϕ1,t1,t2,i)h(ϕ2,t1,t2,i)||f(ϕ1,t1,t2,i)f(ϕ2,t1,t2,i)||g(ϕ1,t1,t2,i)g(ϕ2,t1,t2,i)|L||ϕ1ϕ2||

    for all t ≥ 0, iS and ϕ1, ϕ2C([ − τ, 0];ℝn).

  2. If x(t) is an RCLL and 𝔉Et-adapted process, then h(xt, t, Et, r(t)), f(xt, t, Et, r(t)), g(xt, t, Et, r(t)) ∈ 𝓛(𝔉Et), where 𝓛(𝔉Et) denotes the class of RCLL and 𝔉Et-adapted process.

3 Main results

In this section, we aim to establish the stability results of the system equation (1.1). Firstly, we have to guarantee the existence of the solution of the equation (1.1).

Lemma 3.1

Under the conditions of (H1) and (H2), for any initial data {x(θ) : − τθ ≤ 0} = ξ CF0b ([− τ, 0];ℝn), the equation (1.1) has a unique global solution.

Proof

Let T > 0 be arbitrary. It is known that ([18]) there is a sequence {τk}k≥0 of stopping times such that 0 < τ0 < τ1 < ⋯ < τk → ∞ and r(t) is constant on each interval [τk, τk+1), that is, for each k ≥ 0,

r(t)=r(τk),τkt<τk+1.

We first consider the equation on t ∈ [0, τ1T], it becomes

dx(t)=h(xt,t,Et,r(0))dt+f(xt,t,Et,r(0))dEt+g(xt,t,Et,r(0))dBEt

with initial data x0 = ξ CF0b ([−τ, 0]) has a unique solution on [−τ, τ1T]([4, 8]). Next, for t ∈ [τ1T, τ2T], the equation becomes

dx(t)=h(xt,t,Et,r(τ1T))dt+f(xt,t,Et,r(τ1T))dEt+g(xt,t,Et,r(τ1T))dBEt

with initial data xτ1T given above. Again we know the equation has a unique continuous solution on [τ1Tτ, τ2T]. Repeating the progress, we can see the equation has a unique solution x(t) on [−τ, T]. Since T is arbitrary, the existence and uniqueness have been proved. □

Now, let us consider the exponential stability of equation (1.1). We fix the Markov chain r(t) and let the initial data ξ vary in CF0b ([−τ, 0];ℝn). The solution of equation (1.1) is denoted as x(t; ξ) throughout this paper. Assume that h(0, t, Et, i) = 0, f(0, t, Et, i) = 0, g(0, t, Et, i) = 0, so the equation (1.1) have a trivial solution x(t; 0) = 0. Next, we establish a new Razumikhin theorem on p-th moment exponential stability for the time-changed SFDEs with Markovian switching.

Theorem 3.1

Let (H1) and (H2) hold. Let λ1, λ2, p, c1, c2, α be all positive numbers and q > 1. Assume that there exists a function V(x, t, Et, i) ∈ C2,1,1(Rn × [−τ, ∞) × [0, ∞) × S; R+) such that

c1|x|pV(x,t,Et,i)c2|x|p,(x,t,Et,i)Rn×[τ,)×[0,)×S (3.1)

and for all t > 0,

Emax1iNeαEtLjV(ϕ,t,Et,i)λjEmax1iNeαEtV(ϕ(0),t,Et,i)(j=1,2) (3.2)

provided ϕ = {ϕ(θ; − τθ ≤ 0)} satisfying

Emin1iNeαEt+θV(ϕ(θ),t+θ,Et+θ,i)qEmax1iNeαEtV(ϕ(0),t,Et,i) (3.3)

for allτθ ≤ 0. Then for all ξ CF0b ([−τ, 0], Rn)

E|x(t;ξ)|pc2c1E||ξ||peγt,t0, (3.4)

where γ = min{λ1, λ2, log(q)/τ}. In other words, the trivial solution of equation (1.1) is pth moment exponentially stable and the pth moment Lyapunov exponent is not greater thanγ.

Proof

For the initial data ξ CF0b ([−τ, 0], Rn) arbitrarily and we write x(t; ξ) = x(t) simply. Extend r(t) to [−τ, 0) by setting r(t) = r(0), and extend Et to [−τ, 0) by setting Et = E0. Let ε ∈ (0, γ) be arbitrary then set γ = γε . Define

U(t)=supτθ0Eeγ¯(t+θ+Et+θ)V(x(t+θ),t+θ,Et+θ,r(t+θ))fort0.

Since r(t) is right continuous, the fact that both Et and x(t) is continuous and E(supτst|x(s)|p) < ∞ for t ≥ 0, we can see 𝓔V(x(t), t, Et, r(t)) is right continuous on t ≥ − τ . Hence U(t) is well defined and right continuous. We claim that

D+U(t):=lim supl0+U(t+l)U(t)t0forallt0. (3.5)

To show this, we know that for each t ≥ 0, either U(t) > 𝓔[eγ(t+Et)V(x(t), t, Et, r(t))] or U(t) = 𝓔[eγ(t+Et)V(x(t), t, Et, r(t))].

  1. If U(t) > 𝓔[eγ(t+Et)V(x(t), t, Et, r(t))], it follows from the right continuity of 𝓔[eγ(t+Et)V(x(t), t, Et, r(t))] that for each l > 0 sufficiently small

    U(t)>E[eγ¯(t+l+Et+l)V(x(t+l),t+l,Et+l,r(t+l))]. (3.6)

    Noting that

    U(t+l)=supτθ0Eeγ¯(t+l+θ+Et+l+θ)V(x(t+l+θ),t+l+θ,Et+l+θ,r(t+l+θ))fort0,

    if l + θ > 0, by (3.6), we have

    Eeγ¯(t+l+θ+Et+l+θ)V(x(t+l+θ),t+l+θ,Et+l+θ,r(t+l+θ))U(t).

    Therefore, U(t + l) ≤ U(t). On the other hand, if l + θ ≤ 0, we set θ′ = l + θ, then

    U(t+l)=suplτθ0Eeγ¯(t+θ+Et+θ)V(x(t+θ),t+θ,Et+θ,r(t+θ))supτθ0Eeγ¯(t+θ+Et+θ)V(x(t+θ),t+θ,Et+θ,r(t+θ))=U(t).

    Therefore, for each t > 0, U(t + l) ≤ U(t) and D+ U(t) ≤ 0.

  2. If U(t) = 𝓔[eγ(t+Et)V(x(t), t, Et, r(t))], by the definition of U(t), one obtains that for −τθ ≤ 0,

    Eeγ¯(t+θ+Et+θ)V(x(t+θ),t+θ,Et+θ,r(t+θ))Eeγ¯(t+Et)V(x(t),t,Et,r(t)),

    it follows that

    Eeγ¯Et+θV(x(t+θ),t+θ,Et+θ,r(t+θ))eγ¯θEeγ¯EtV(x(t),t,Et,r(t))eγ¯τEeγ¯EtV(x(t),t,Et,r(t)).

    If 𝓔 [eγEtV(x(t), t, Et, r(t))] = 0, from (3.1) we can see that

    E[eγ¯Et+θc1|x(t+θ)|p]0,

    which yields that x(t + θ) = 0, −τθ ≤ 0. Since h(0, t, Et, i) = 0, f(0, t, Et, i) = 0 and g(0, t, Et, i) = 0 a.s. for all −τθ ≤ 0, one obtains that x(t + l) = 0 a.s. for all l > 0, hence U(t + l) = 0 and D+ U(t) = 0.

On the other hand, if 𝓔 [eγEtV(x(t), t, Et, r(t))] > 0, one can see that

Eeγ¯Et+θV(x(t+θ),t+θ,Et+θ,r(t+θ))<qEeγ¯EtV(x(t),t,Et,r(t))

for all −τθ ≤ 0 since eγτ < q. It follows from the condition (3.2) that

Emax1iNeγ¯EtLjV(ϕ,t,Et,i)<λjEmax1iNeγ¯EtV(ϕ(0),t,Et,i),j=1,2.

It means that

Eeγ¯EtLjV(xt,t,Et,r(t))<λjEeγ¯EtV(x(t),t,Et,r(t)),j=1,2,

then

Eeγ¯Et(γ¯V(x(t),t,Et,r(t))+LjV(xt,t,Et,r(t)))(λjγ¯)E[eγ¯EtV(x(t),t,Et,r(t))]<0.

By the right continuity of the process involved one can see that for all l > 0 sufficiently small,

Eeγ¯Es(γ¯V(x(s),s,Es,r(s))+LjV(xs,s,Es,r(s)))0,tst+l.

By the generalized time-changed Itô formula, we get that

Eeγ¯(t+l+Et+l)V(x(t+l),t+l,Et+l,r(t+l))Eeγ¯(t+Et)(V(x(t),t,Et,r(t)))=Ett+leγ¯(s+Es)[γ¯V(x(s),s,Es,r(s))+L1V(xs,s,Es,r(s))]ds+Ett+leγ¯(s+Es)[γ¯V(x(s),s,Es,r(s))+L2V(xs,s,Es,r(s))]dEs=tt+leγ¯sEeγ¯Es[γ¯V(x(s),s,Es,r(s))+L1V(xs,s,Es,r(s))]ds+tt+leγ¯sEeγ¯Es[γ¯V(x(s),s,Es,r(s))+L2V(xs,s,Es,r(s))]dEs0. (3.7)

Then U(t + l) ≤ U(t) for l > 0 sufficiently small.

Since

U(t+l)=supτθ0Eeγ¯(t+θ+l+Et+θ+l)V(x(t+l+θ),t+l+θ,Et+l+θ,r(t+l+θ)),

here we set θ′ = θ + l, if l + θ > 0, then 𝓔 [eγ(t+θ′+Et+θ)V(x(t+θ′), t + θ′, Et+θ, r(t + θ′))] ≤ U(t) from (3.7), otherwise, since U(t) = 𝓔[eγ(t+Et)V(x(t), t, Et, r(t))], then

Eeγ¯(t+θ+Et+θ)V(x(t+θ),t+θ,Et+θ,r(t+θ))U(t),

so, by the definition of supremum, U(t + l) = U(t) for l > 0 sufficiently small and D+ U(t) = 0. Therefore, the inequality (3.5) has been proved. It follows that

U(t)U(0),fort0.
Eeγ¯tc1|x|pEeγ¯(t+Et)V(x(t),t,Et,r(t))U(t)U(0)c2E||ξ||p

this means

E|x|pc2c1eγ¯tE||ξ||p=c2c1E||ξ||pe(γε)t.

Since ε is arbitrary, the required inequality (3.4) must hold. The proof is completed. □

4 Stochastic delay differential equations with Markovian switching

In this section, as a special case of equation (1.1), we consider the time-changed stochastic delay differential equation with Marking switching as follows,

dx(t)=H(x(t),x(tδ(t)),t,Et,r(t))dt+F(x(t),x(tδ(t)),t,Et,r(t))dEt+G(x(t),x(tδ(t)),t,Et,r(t))dBEt (4.1)

on t ≥ 0 with x0 = ξ CF0b ([−τ, 0];ℝn), where δ:ℝ+ → [0, τ] is Borel measure while

H,F:Rn×Rn×R+×R+×SRn

and

G:Rn×Rn×R+×R+×SRn×m.

We impose the following hypotheses:

  1. Both H, F : ℝn × ℝ+ × ℝ+ × S → ℝn and G : ℝn × ℝ+ × ℝ+ × S → ℝn×m are Borel-measurable functions. They satisfy the Lipschitz condition. That is, there is L > 0 such that

    |H(x,y,t1,t2,i)H(x¯,y¯,t1,t2,i)||F(x,y,t1,t2,i)F(x¯,y¯,t1,t2,i)||G(x,y,t1,t2,i)G(x¯,y¯,t1,t2,i)|L(|xx¯|+|yy¯|)

    for all t ≥ 0, iS and x, y, x, y ∈ ℝn.

  2. If x(t) is an RCLL and 𝓕Et-adapted process, then H(x(t), x(tδ(t)), t, Et, r(t)), F(x(t), x(tδ(t)), t, Et, r(t)), G(x(t), x(tδ(t)), t, Et, r(t)) ∈ 𝓛(𝔉Et), where 𝓛(𝔉Et) denotes the class of RCLL and 𝔉Et-adapted process.

If we define, for (ϕ, t, Ei, i) ∈ C([−τ, 0];ℝn) × ℝ+ × ℝ+ × S,

h(ϕ,t,Et,i)=H(ϕ(0),ϕ(δ(t)),t,Et,i),g(ϕ,t,Et,i)=G(ϕ(0),ϕ(δ(t)),t,Et,i),f(ϕ,t,Et,i)=F(ϕ(0),ϕ(δ(t)),t,Et,i),

then the equation (4.1) becomes the equation (1.1) and (H3) (H4) imply (H1) (H2). So, by Lemma 3.1, the equation (4.1) has a unique global solution which is again denoted by x(t; ξ). Furthermore, assume that H(0, 0, t, Et, i) = 0, F(0, 0, t, Et, i) = 0, G(0, 0, t, Et, i) = 0.

If VC2,1,1(ℝn × [−τ, ∞) × [0, ∞) × S;ℝ+), define L1V and L2V from ℝn × ℝn × ℝ+ × ℝ+ × S to ℝ respectively by

L1V(x,y,t,Et,i)=Vt(x,t,Et,i)+Vx(x,t,Et,i)H(x,y,t,Et,i)+j=1NγijV(x,t,Et,j),
L2V(x,y,t,Et,i)=VEt(x,t,Et,i)+Vx(x,t,Et,i)F(x,y,t,Et,i)+12trGTVxxG(x,y,t,Et,i).

Furthermore, we denote LFtp (Ω, ℝn) as the family of all 𝔉t-measurable ℝn-valued random variables X such that E|X|p < ∞. Meanwhile, we set

LjV(ϕ,t,Et,i)=LjV(ϕ(0),ϕ(δ(t)),t,Et,i),j=1,2

Theorem 4.1

Let (H3) and (H4) hold. Let λ1, λ2, p, c1, c2, α be all positive numbers and q > 1. Assume that there exists a function V(x, t, Et, i) ∈ C2,1,1(Rn × [−τ, ∞) × [0, ∞) × S; R+) such that

c1|x|pV(x,t,Et,i)c2|x|p,(x,t,Et,i)Rn×[τ,)×[0,)×S (4.2)

and for all t > 0,

Emax1iNeαEtLjV(X,Y,t,Et,i)λjEmax1iNeαEtV(X,t,Et,i)(j=1,2) (4.3)

provided X, Y LFtp (Ω, ℝn) satisfying

Emin1iNeαEt+θV(Y,tδ(t),Etδ(t),i)qEmax1iNeαEtV(X,t,Et,i) (4.4)

Then for all ξ CF0b ([−τ, 0], Rn)

E|x(t;ξ)|pc2c1E||ξ||peγt,t0, (4.5)

where γ = min{λ1, λ2, log(q)/τ}. In other words, the trivial solution of equation (4.1) is pth moment exponentially stable and the pth moment Lyapunov exponent is not greater thanγ.

Proof

Let ϕ = {ϕ(θ) : − τθ ≤ 0} ∈ LFtp ([−τ, 0], ℝn) satisfy (3.3). For X = ϕ(0), Y = ϕ(−δ(t)) ∈ LFtp (Ω, ℝn) satisfying

Emin1iNeαEt+θV(ϕ(δ(t)),tδ(t),Etδ(t),i)qEmax1iNeαEtV(ϕ(0),t,Et,i).

Then, from (4.3) we have

Emax1iNeαEtLjV(ϕ,t,Et,i)λjEmax1iNeαEtV(ϕ(0),t,Et,i)(j=1,2)

which is (3.2). Hence the conditions in Theorem 3.1 are satisfied and the conclusions follow. Applying the Theorem 3.1, the proof is completed. □

Theorem 4.2

Let (H3) and (H4) hold. Let p, c1, c2, α be all positive numbers and λ1j > λ2j ≥ 0, j = 1, 2. Assume that there exists a function V(x, t, Et, i) ∈ C2,1,1(Rn × [−τ, ∞) × [0, ∞) × S; R+) such that

c1|x|pV(x,t,Et,i)c2|x|p,(x,t,Et,i)Rn×[τ,)×[0,)×S (4.6)

and for all t > 0,

Emax1iNeαEtLjV(X,Y,t,Et,i)λ1jEmax1iNeαEtV(X,t,Et,i)+λ2jEmin1iNeαEt+θV(Y,tδ(t),Etδ(t),i)(j=1,2)

Then the trivial solution of equation (4.1) is pth moment exponentially stable and the pth moment Lyapunov exponent is not greater thanγ, where γ = min{λ1121, λ1222, log(q)/τ} with q > 1.

Proof

For t ≥ 0, q < λ1j/λ2j, j = 1, 2 and X, Y LFtp (Ω, ℝn) satisfying

Emin1iNeαEt+θV(Y,tδ(t),Etδ(t),i)qEmax1iNeαEtV(X,t,Et,i),

we can arrive that

Emax1iNeαEtLjV(X,Y,t,Et,i)λ1jEmax1iNeαEtV(X,t,Et,i)+λ2jEmin1iNeαEt+θV(Y,tδ(t),Etδ(t),i)(λ1jqλ2j)Emax1iNeαEtV(X,t,Et,i),

that is, (4.3) is satisfied with λj = λ1j2j, j = 1, 2. Then the conclusion follows form Theorem 4.1. □

5 Example

Let Et be generalized inverse of an β -stable subordinator Uβ(t). Let B(t) be a scalar Brownian motion and {r(t)} be a right-continuous Markov chain taking values in S = {1, 2} with generator Γ = {rij}2×2, here

γ11=γ12>0,γ21=γ22>0.

Assume that B(t) and r(t) are independent. Then let us consider the following one-dimensional linear stochastic differential equation with Markovian switching

dx(t)=ρ(r(t))x(t)dt+μ(r(t))x(tδ(t))dEt+σ(r(t))x(tδ(t))dBEt,t0 (5.1)

where

ρ(1)=1,ρ(2)=1;μ(1)=12,μ(2)=13;σ(1)=1,σ(2)=1.

The equation (5.1) can be regarded as the result of

dx(t)=x(t)dt12x(tδ(t))dEt+x(tδ(t))dBEt,t0 (5.2)

and

dx(t)=x(t)dt13x(tδ(t))dEt+x(tδ(t))dBEt,t0 (5.3)

switching to each other according to the movement of the Markovian chain r(t).

We define the function V : ℝ × ℝ+ × ℝ+ × S → ℝ+ by

V(x,t,Et,i)=ci|x|p

with ci = 1, c2 = c ∈ (0, 34 ). The operators have the following forms

L1V(x,t,Et,i)=(c1p)|x|p,i=1,(pc+44c)|x|p,i=2.
L2V(x,t,Et,i)=p(p1)2|x|p2|y|212p|x|p1|y|,i=1,cp(p1)2|x|p2|y|213cp|x|p1|y|,i=2.

Using the following inequality

aθb1θθa+(1θ)b,a,b>0,θ(0,1),

we can see that

L2V(x,t,Et,i)(p1)(p3)2|x|p+(p32)|y|p,i=1,c(p1)(p8)2|x|p+c(p4)|y|p,i=2.

Choose p = 2, 2 < c < 3, then

L1V(x,t,Et,i)=(c3)|x|pi=1,(42c)|x|p,i=2min{3c,2c4c}max{V(x,t,Et,1),V(x,t,Et,2)}.
L2V(x,t,Et,i)12|x|p+12|y|p,i=1,c3|x|p+2c3|y|p,i=212cmaxV(x,t,Et,1),V(x,t,Et,2)+23min{V(x,t,Et,1),V(x,t,Et,2)}.

By the Theorem 4.2 we conclude that the trivial solution of the equation (5.1) is pth moment exponentially stable.

6 Conclusions

The stochastic differential equations(SDEs) driven by time-changed Brownian motions is a new research area for recent years. In this paper, we have studied the exponential stability of the time-changed SDEs with Markovian switching, by expanding the time-changed Itô formula and the time-changed Razumikhin theorem. Our result generalizes that of SDEs in the literature. Due to the more construction of SDEs with time-change than the usual SDEs, our result is not a trivial generalization.

Acknowledgement

The Project is supported by the National Natural Science Foundation of China (11861040, 11771198, 11661053).

References

[1] Bakera C., Buckwarb E., Exponential stability in p-th mean of solutions, and of convergent Euler-type solutions, of stochastic delay differential equations. J. Comput. Appl. Math., 2005, 184, 404-42710.1016/j.cam.2005.01.018Search in Google Scholar

[2] Mao X., Stochastic functional differntial equations with Markovian switching, Funct. Differ. Equ., 1999, 6(3-4), 375-396Search in Google Scholar

[3] Mao X., Stochastic differential equations and applications, Horwood Publing, UK, 200710.1533/9780857099402Search in Google Scholar

[4] Protter P, Stochastic integration and differential equations, Springer-verlag, Berlin Heidelberg, 200510.1007/978-3-662-10061-5Search in Google Scholar

[5] Tan L., Jin W., Suo Y.Q., Stability in distribution of neutral stochastic functional differential equations, Statist. Probab. Lett., 2015, 107, 27-3610.1016/j.spl.2015.07.033Search in Google Scholar

[6] Chlebak L., Garmirian P., Wu Q., Sub-diffusion processes in Hilbert space and their associated stochastic differential equations and Fokker-Planck-Kolmogorov equations, arXiv:1610.00208v1Search in Google Scholar

[7] Bertoin J., Subordinators: examples and applications. In: Lectures on Probability Theory and Statistics, Berlin: Springer, 1999, 1-9110.1007/b72002Search in Google Scholar

[8] Kobayashi K., Stochastic calculus for a time-changed semimartingale and the associated stochastic differential equations. J. Theor. Probab., 2011, 24, 789-82010.1007/s10959-010-0320-9Search in Google Scholar

[9] Meerschaert M. M., Straka P., Inverse stable subordinators, Mathematical modelling of natural phenomena, 2013, 8(02), 1-1610.1051/mmnp/20138201Search in Google Scholar PubMed PubMed Central

[10] Hahn M., Kobayashi K., Umarov S., SDEs driven by a time-changed lévy process and their associated time-fractional order pseudo-differential equations, J. Theoret. Probab., 2012, 25(1), 262-27910.1007/s10959-010-0289-4Search in Google Scholar

[11] Jum E., Kobayashi K., A strong and weak approximation scheme for stochastic differential equations driven by a time-changed Brownian motion, Probab. Math. Statist., 2016, 36(2), 201-220Search in Google Scholar

[12] Wu Q., Stability analysis for a class of nonlinear time-changed systems, Cogent Math., 2016, 3, 1228273.10.1080/23311835.2016.1228273Search in Google Scholar

[13] Wu Q., Stability of stochastic differential equations with respect to time-changed Brownian motions, ArXiv: 1602.08160v1Search in Google Scholar

[14] Nane E., Ni Y., Stability of the solution of stochastic differential equation driven by time-changed lévy noise, P. Am. Math. Soc., 2017, 145, 3085-310410.1090/proc/13447Search in Google Scholar

[15] Mao X., Razumikhin-type theorems on exponential stability of stochastic functional differential equations, Stochastic Processes Appl., 1996, 65, 233-25010.1016/S0304-4149(96)00109-3Search in Google Scholar

[16] Zhu Q.X., Razumikhin-type theorem for stochastic functional differential equations with lévy noise and Markov switching, International J. Control, 2017, 90(8), 1703-171210.1080/00207179.2016.1219069Search in Google Scholar

[17] Magdziarz M., Path properties of subdiffusion - a martingale approach, Stoch. Models, 2010, 26(2), 256-27110.1080/15326341003756379Search in Google Scholar

[18] Skorohod A.V., Asymptotic methods in the theory of stochastic differential equations, Providence, RI: American Methmatical Society, 1989Search in Google Scholar

Received: 2018-07-21
Accepted: 2019-05-16
Published Online: 2019-07-09

© 2019 Zhang and Yuan, published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 Public License.

Downloaded on 27.4.2024 from https://www.degruyter.com/document/doi/10.1515/math-2019-0055/html
Scroll to top button