We investigate a class of stochastic partial differential equations with Markovian switching. By using the Euler-Maruyama scheme both in time and in space of mild solutions, we derive sufficient conditions for the existence and uniqueness of the stationary distributions of numerical solutions. Finally, one example is given to illustrate the theory.
1. Introduction
The theory of numerical solutions of stochastic partial differential equations (SPDEs) has been well developed by many authors [1–5]. In [2], Debussche considered the error of the Euler scheme for the nonlinear stochastic partial differential equations by using Malliavin calculus. Gyöngy and Millet [3] discussed the convergence rate of space time approximations for stochastic evolution equations. Shardlow [5] investigated the numerical methods of the mild solutions for stochastic parabolic PDEs derived by space-time white noise by applying finite difference approach.
On the other hand, the parameters of SPDEs may experience abrupt changes caused by phenomena such as component failures or repairs, changing subsystem interconnections, and abrupt environmental disturbances [6–9], and the continuous-time Markov chains have been used to model these parameter jumps. An important equation is a class of SPDEs with Markovian switching
(1)dX(t)=[AX(t)+f(X(t),r(t))]dt+g(X(t),r(t))dW(t),t≥0.
Here the state vector has two components X(t) and r(t), the first one is normally referred to as the state while the second one is regarded as the mode. In its operation, the system will switch from one mode to another one in a random way, and the switching among the modes is governed by the Markov chain r(t).
Since only a few SPDEs with Markovian switching have explicit formulae, numerical (approximate) schemes of SPDEs with Markovian switching are becoming more and more popular. In this paper, we will study the stationary distribution of numerical solutions of SPDEs with Markovian switching. Bao et al. [10] investigated the stability in distribution of mild solutions to SPDEs. Bao and Yuan [11] discussed the numerical approximation of stationary distribution for SPDEs. For the stationary distribution of numerical solutions of stochastic differential equations in finite-dimensional space, Mao et al. [12] utilized the Euler-Maruyama scheme with variable step size to obtain the stationary distribution and they also proved that the probability measures induced by the numerical solutions converge weakly to the stationary distribution of the true solution. But since the mild solutions of SPDEs with Markovian switching do not have stochastic differential, a significant consequence of this fact is that the Itô formula cannot be used for mild solutions of SPDEs with Markovian switching directly. Consequently, we generalize the stationary distribution of numerical solutions of the finite dimensional stochastic differential equations with Markovian switching to that of infinite dimensional cases.
Motived by [11–13], we will show in this paper that the mild solutions of SPDE with Markovian switching (1) have a unique stationary distribution for sufficiently small step size. So this paper is organised as follows: in Section 2, we give necessary notations and define Euler-Maruyama scheme of mild solutions. In Section 3, we give some lemmas and the main result in this paper. Finally, we will give an example to illustrate the theory in Section 4.
2. Statements of Problem
Throughout this paper, unless otherwise specified, we let (Ω,ℱ,{ℱt}t≥0,ℙ) be complete probability space with a filtration {ℱt}t≥0 satisfying the usual conditions (i.e., it is increasing and right continuous while ℱ0 contains all ℙ-null sets). Let (H,〈·,·〉H,∥·∥H) be a real separable Hilbert space and W(t) an H-valued cylindrical Brownian motion (Wiener process) defined on the probability space. Let IG be the indicator function of a set G. Denote by (ℒ(H),∥·∥) and (ℒHS(H),∥·∥HS) the family of bounded linear operators and Hilbert-Schmidt operator from H into H, respectively. Let r(t), t≥0, be a right-continuous Markov chain on the probability space taking values in a finite state space 𝕊={1,2,…,N} with the generator Γ=(γij)N×N given by
(2)ℙ{r(t+δ)=j∣r(t)=i}={γijδ+o(δ)ifi≠j,1+γijδ+o(δ)ifi=j,
where δ>0. Here γij>0 is the transition rate from i to j if i≠j while
(3)γii=-∑j≠iγij.
We assume that the Markov chain r(·) is independent of the Brownian motion W(·). It is well known that almost every sample path of r(·) is a right-continuous step function with finite number of simple jumps in any finite subinterval of ℝ+:=[0,+∞).
Consider SPDEs with Markovian switching on H(4)dX(t)=[AX(t)+f(X(t),r(t))]dt+g(X(t),r(t))dW(t),t≥0,
with initial value X(0)=x∈H and r(0)=i∈𝕊. Here f:H×𝕊→H, g:H×𝕊→ℒHS(H). Throughout the paper, we impose the following assumptions.
(A,𝒟(A)) is a self-adjoint operator on H generating a C0-semigroup {eAt}t≥0, such that ∥eAt∥≤e-αt for some α>0. In this case, -A has discrete spectrum 0<ρ1≤ρ2≤⋯≤limi→∞ρi=∞ with corresponding eigenbasis {ei}i≥1 of H.
Both f and g are globally Lipschitz continuous. That is, there exists a constant L>0 such that
(5)∥f(x,j)-f(y,j)∥H2∨∥g(x,j)-g(y,j)∥HS2≤L∥x-y∥H2,∀x,y∈H,j∈𝕊;
There exist μ>0 and λj>0, (j=1,2,…,N) such that
(6)2λj〈x-y,f(x,j)-f(y,j)〉H+λj∥g(x,j)-g(y,j)∥HS2+∑l=1Nγjlλl∥x-y∥H2≤-μ∥x-y∥H2,∀x,y∈H,j∈𝕊.
It is well known (see [1, 8]) that under (A1)–(A3), (4) has a unique mild solution X(t) on t≥0. That is, for any X(0)=x∈H and r(0)=i∈𝕊, there exists a unique H-valued adapted process X(t) such that(7)X(t)=etAx+∫0te(t-s)Af(X(s),r(s))ds+∫0te(t-s)Ag(X(s),r(s))dW(s).
Moreover, the pair Z(t)=(X(t),r(t)) is a time-homogeneous Markov process.
Remark 1.
We observe that (A2) implies the following linear growth conditions:
(8)|f(x,j)|H2∨∥g(x,j)∥HS2≤L-(1+∥x∥H2),∀x∈H,j∈𝕊,
where L-=2maxj∈𝕊(L∨|f(0,j)|H2∨∥g(0,j)∥HS2).
Remark 2.
We also establish another property from (A3): (9)2λj〈x,f(x,j)〉H+λj∥g(x,j)∥HS2+∑l=1Nγjlλl∥x∥H2≤2λj〈x,f(x,j)-f(0,j)〉H+λj∥g(x,j)-g(0,j)∥HS2+∑l=1Nγjlλl∥x∥H2+2λj〈x,f(0,j)〉H+2λj〈g(x,j)-g(0,j),g(0,j)〉HS+λj∥g(0,j)∥HS2≤-μ∥x∥H2+μ4∥x∥H2+4λj2∥f(0,j)∥Hμ+μ4L∥g(x,j)-g(0,j)∥HS2+4Lλj2μ∥g(0,j)∥HS2+λj∥g(0,j)∥HS2≤-μ∥x∥H2+μ4∥x∥H2+μ4∥x∥H2+4λj2∥f(0,j)∥Hμ+4Lλj2μ∥g(0,j)∥HS2+λj∥g(0,j)∥HS2≤-μ2∥x∥H2+α1,∀x∈H,j∈𝕊,
where α1:=maxj∈𝕊[(4λj2∥f(0,j)∥H2/μ)+(4Lλj2/μ)∥g(0,j)∥HS2+λj∥g(0,j)∥HS2] and 〈T,S〉HS:=∑i=1∞〈Tei,Sej〉H for S,T∈ℒHS(H).
Denote by Zx,i(t)=(Xx,i(t),ri(t)) the mild solution of (4) starting from (x,i)∈H×𝕊. For any subset A∈𝔅(H),B⊂𝕊, let ℙt((x,i),A×B) be the probability measure induced by Zx,i(t),t≥0. Namely,
(10)ℙt((x,i),A×B)=ℙ(Zx,i∈A×B),
where 𝔅(H) is the family of the Borel subset of H.
Denote by 𝒫(H×𝕊) the family by all probability measures on H×𝕊. For P1,P2∈𝒫(H×𝕊), define the metric d𝕃 as follows:
(11)d𝕃(P1,P2)=supφ∈𝕃|∑j=1N∫Hφ(u,j)P1(du,j)-∑j=1N∫φ(u,j)P2(du,j)|,
where 𝕃={φ:H×𝕊→ℝ:|φ(u,j)-φ(v,l)|≤∥u-v∥H+|j-l|, and |φ(u,j)|≤1, for u,v∈K,j,l∈𝕊}.
Remark 3.
It is known that the weak convergence of probability measures is a metric concept with respect to classes of test function. In other words, a sequence of probability measures {Pk}k≥1 of 𝒫(H×𝕊) converges weakly to a probability measure P0∈𝒫(H×𝕊) if and only if limk→∞d𝕃(Pk,P0)=0.
Definition 4.
The mild solution Z(t)=(X(t),r(t)) of (4) is said to have a stationary distribution π(·×·)∈𝒫(H×𝕊) if the probability measure ℙt((x,i),(·×·)) converges weakly to π(·×·) as t→∞ for every i∈𝕊, and every x∈U, a bounded subset of H, that is,
(12)limt→∞d𝕃(ℙt(x,i),π(·×·))=limt→∞(supφ∈𝕃|∑j=1N𝔼φ(Zx,i(t))=limt→∞==limt→∞=-∑j=1N∫Hφ(u,j)π(du,j)|)=0.
By Theorem 3.1 in [10] and Theorem 3.1 in [14], we have the following.
Theorem 5.
Under (A1)–(A3), the Markov process Z(t) has a unique stationary distribution π(·×·)∈𝒫(H×𝕊).
For any n≥1, let πn:H→Hn:=Span{e1,e2,…,en} be the orthogonal projection. Consider SPDEs with Markovian switching on Hn,
(13)dXn(t)=[AnXn(t)+fn(Xn(t),r(t))]dt+gn(Xn(t),r(t))dW(t),
with initial data Xn(0)=πnx=∑i=1n〈x,ei〉Hei, x∈H. Here An=πnA, fn=πnf, gn=πng.
Therefore, we can observe that
(14)Anx=Ax,etAnx=etAx,〈x,fn〉H=〈x,f〉H,〈x,gn〉H=〈x,g〉H,∀x∈Hn,
By the property of the projection operator and (A2), we have
(15)∥An(x-y)∥H2∨∥fn(x,j)-fn(y,j)∥H2∨∥gn(x,j)-gn(y,j)∥HS2≤λn2∥(x-y)∥H2∨∥f(x,j)-f(y,j)∥H2∨∥g(x,j)-g(y,j)∥HS2≤(λn2∨L)∥(x-y)∥H2,∀x,y∈Hn,j∈𝕊.
Hence, (13) admits a unique strong solution {Xn(t)}t≥0 on Hn (see [8]).
We now introduce an Euler-Maruyama based computational method. The method makes use of the following lemma (see [15]).
Lemma 6.
Given Δ>0, then {r(kΔ),k=0,1,2,…} is a discrete Markov chain with the one-step transition probability matrix
(16)P(Δ)=(Pi,j(Δ))N×N=eΔΓ.
Given a fixed step size Δ>0 and the one-step transition probability matrix P(Δ) in (16), the discrete Markov chain {r(kΔ),k=0,1,2,…} can be simulated as follows: let r(0)=i0, and compute a pseudorandom number ξ1 from the uniform (0,1) distribution.
Define
(17)r(Δ)={i,i∈𝕊-{N}suchthat∑j=1i-1Pr(0),j(Δ)≤ξ1<∑j=1iPr(0),j(Δ),N,∑j=1N-1Pr(0),j(Δ)≤ξ1,
where we set ∑j=10Pr(0),j(Δ)=0 as usual. Having computed r(0),r(Δ),…,r(kΔ), we can compute r((k+1)Δ) by drawing a uniform (0,1) pseudorandom number ξk+1 and setting
(18)r((k+1)Δ)={i,i∈𝕊-{N}suchthat∑j=1i-1Pr(kΔ),j(Δ)≤ξk+1<∑j=1iPr(kΔ),j(Δ),N,∑j=1N-1Pr(kΔ),j(Δ)≤ξk+1.
The procedure can be carried out repeatedly to obtain more trajectories.
We now define the Euler-Maruyama approximation for (13). For a stepsize Δ∈(0,1), the discrete approximation Y-n(kΔ)≈Xn(kΔ), is formed by simulating from Y-n(0)=πnx,r(0)=r0, and
(19)Y-n((k+1)Δ)=eΔAn{Y-n(kΔ)+fn(Y-n(kΔ),r(kΔ))Δ+gn(Y-n(kΔ),r(kΔ))ΔWk},
where ΔWk=W((k+1)Δ)-W(kΔ)).
To carry out our analysis conveniently, we give the continuous Euler-Maruyama approximation solution which is defined by
(20)Yn(t)=etAnπnx+∫0te(t-⌊s⌋)Anfn(Yn(⌊s⌋),r(⌊s⌋))ds+∫0te(t-⌊s⌋)Angn(Yn(⌊s⌋),r(⌊s⌋))dW(s)=etAnπnx+∫0te(t-⌊s⌋)Afn(Yn(⌊s⌋),r(⌊s⌋))ds+∫0te(t-⌊s⌋)Agn(Yn(⌊s⌋),r(⌊s⌋))dW(s),
where ⌊t⌋=[t/Δ]Δ and [t/Δ] denotes the integer part of t/Δ and Yn(0)=Y-n(0)=πnx, and Yn(kΔ)=Y-n(kΔ).
It is obvious that Yn(t) coincides with the discrete approximation solution at the gridpoints. For any Borel set A∈𝔅(Hn),x∈Hn,i,j∈𝕊, let Z-n(kΔ)=(Y-n(kΔ),r(kΔ)),
(21)ℙn,Δ((x,i),A×{j}):=ℙ(Z-n(Δ)∈A×{j}∣Z-n(0)=(x,i)),ℙkn,Δ((x,i),A×{j}):=ℙ(Z-n(kΔ)∈A×{j}∣Z-n(0)=(x,i)).
Following the argument of Theorem 5 in [13], we have the following.
Lemma 7.
{Z-n(kΔ)}k≥0 is a homogeneous Markov process with the transition probability kernel ℙn,Δ((x,i),A×{j}).
To highlight the initial value, we will use notation {Z-n,(x,i)(kΔ)}.
Definition 8.
For a given stepsize Δ>0, {Z-n,(x,i)(kΔ)}k≥0 is said to have a stationary distribution {πn,Δ(·×·)}∈𝒫(Hn×𝕊) if the k-step transition probability kernel ℙkn,Δ((x,i),·×·) converges weakly to πn,Δ(·×·) as k→∞, for every (x,i)∈Hn×𝕊, that is,(22)limk→∞d𝕃(Pkn,Δ((x,i),·×·),πn,Δ(·×·))=0.
We will establish our result of this paper in Section 3.
Theorem 9.
Under (A1)–(A3), for a given stepsize Δ>0, and arbitrary x∈Hn, i∈𝕊, {Z-n,(x,i)(kΔ)}k≥0 has a unique stationary distribution πn,Δ(·×·)∈𝒫(Hn×𝕊).
3. Stationary in Distribution of Numerical Solutions
In this section, we shall present some useful lemmas and prove Theorem 9. In what follows, C>0 is a generic constant whose values may change from line to line.
For any initial value (x,i), let Yn,x,i(t) be the continuous Euler-Maruyama solution of (20) and starting from (x,i)∈H×𝕊. Let Xx,i(t) be the mild solution of (4) and starting from (x,i)∈H×𝕊.
Lemma 10.
Under (A1)–(A3), then
(23)𝔼∥Yn,x,i(t)-Yn,x,i(⌊t⌋)∥H2≤3(ρn2+2L-)Δ(1+𝔼∥Yn(⌊t⌋)∥H2).
Proof.
Write Yn,x,i(t)=Yn(t), Yn,x,i(⌊t⌋)=Yn(⌊t⌋). From (20), we have
(24)Yn(⌊t⌋)=e⌊t⌋Aπnx+∫0⌊t⌋e(⌊t⌋-⌊s⌋)Afn(Yn(⌊s⌋),r(⌊s⌋))ds+∫0⌊t⌋e(⌊t⌋-⌊s⌋)Agn(Yn(⌊s⌋),r(⌊s⌋))dW(s).
Thus,
(25)Yn(t)-Yn(⌊t⌋)=e(t-⌊t⌋)A×(e⌊t⌋Aπnx×+×××+∫0⌊t⌋e(⌊t⌋-⌊s⌋)A××××++×××fn(Yn(⌊s⌋),r(⌊s⌋))ds+∫0⌊t⌋e(⌊t⌋-⌊s⌋)A×gn(Yn(⌊s⌋),r(⌊s⌋))dW(s)e⌊t⌋A)-Yn(⌊t⌋)+∫⌊t⌋te(t-⌊s⌋)Afn(Yn(⌊s⌋),r(⌊s⌋))ds+∫⌊t⌋te(t-⌊s⌋)Agn(Yn(⌊s⌋),r(⌊s⌋))dW(s)=(e(t-⌊t⌋)A-1)Yn(⌊t⌋)+∫⌊t⌋te(t-⌊s⌋)Afn(Yn(⌊s⌋),r(⌊s⌋))ds+∫⌊t⌋te(t-⌊s⌋)Agn(Yn(⌊s⌋),r(⌊s⌋))dW(s).
Then, by the Hölder inequality and the Itô isometry, we obtain
(26)𝔼∥Yn(t)-Yn(⌊t⌋)∥H2≤3{∫⌊t⌋t𝔼∥(e(t-⌊t⌋)A-1)Yn(⌊t⌋)∥H2+𝔼∫⌊t⌋t∥fn(Yn(⌊s⌋),r(⌊s⌋))∥H2ds+𝔼∫⌊t⌋t∥gn(Yn(⌊s⌋),r(⌊s⌋))∥HS2ds}.
From (A1), we have
(27)𝔼∥(e(t-⌊t⌋)A-1)Yn(⌊t⌋)∥H2=𝔼∥∑i=1n(e-ρi(t-⌊t⌋)-1)〈Yn(⌊t⌋),ei〉Hei∥H2≤(1-e-ρn(t-⌊t⌋))2𝔼∥Yn(⌊t⌋)∥H2≤ρn2Δ2𝔼∥Yn(⌊t⌋)∥H2,
here we use the fundamental inequality 1-e-a≤a, a>0. And, by (8), it follows that
(28)𝔼∫⌊t⌋t∥fn(Yn(⌊s⌋),r(⌊s⌋))∥H2ds+𝔼∫⌊t⌋t∥gn(Yn(⌊s⌋),r(⌊s⌋))∥HS2ds≤2L-Δ(1+𝔼∥Yn(⌊t⌋)∥H2).
Substituting (27) and (28) into (26), the desired assertion (23) follows.
Lemma 11.
Under (A1)–(A3), if Δ<min{1,1/3(ρn2+2L-),((4αp+μ)/(8q+4qρn2L-+4qL-+24qr^+6qL(ρn2+2L-)))2}, then there is a constant C>0 that depends on the initial value x but is independent of Δ, such that the continuous Euler-Maruyama solution of (20) has
(29)supt≥0𝔼∥Yn,x,i(t)∥H≤C,
where q=max1≤i≤Nλi, p=min1≤i≤Nλi.
Proof.
Write Yn,x,i(t)=Yn(t), ri(kΔ)=r(kΔ). From (20), we have the following differential form:
(30)dYn(t)={AYn(t)+e(t-⌊t⌋)Afn(Yn(⌊t⌋),r(⌊t⌋))}dt+e(t-⌊t⌋)Agn(Yn(⌊t⌋),r(⌊t⌋))dW(t),
with Yn(0)=πnx.
Let V(x,i)=λi∥x∥H2. By the generalised Itô formula, for any θ>0, we derive from (30) that
(31)eθt𝔼(λr(t)∥Yn(t)∥H2)≤λi∥x∥H2+𝔼∫0teθsλr(s)×{∥gn(Yn(⌊s⌋),r(⌊s⌋))∥HS2θ∥Yn(s)∥H2+2〈Yn(s),AYn(s)〉H+2〈Yn(s),e(s-⌊s⌋)Afn(Yn(⌊s⌋),r(⌊s⌋))〉H+∥gn(Yn(⌊s⌋),r(⌊s⌋))∥HS2}ds+𝔼∫0teθs∑l=1Nγr(s)lλl∥Yn(s)∥H2ds≤q∥x∥H2+θq𝔼∫0teθs∥Yn(s)∥H2ds-2αp𝔼∫0teθs∥Yn(s)∥H2ds+𝔼∫0teθsλr(s)×{2〈Yn(s),e(s-⌊s⌋)Af(Yn(⌊s⌋),r(⌊s⌋))〉H+∥g(Yn(⌊s⌋),r(⌊s⌋))∥HS2〈Yn(s),e(s-⌊s⌋)Af(Yn(⌊s⌋),r(⌊s⌋))〉H}ds+𝔼∫0teθs∑l=1Nγr(s)lλl∥Yn(s)∥H2ds.
By the fundamental transformation, we obtain that
(32)〈Yn(t),e(t-⌊t⌋)Af(Yn(⌊t⌋),r(⌊t⌋))〉H=〈Yn(t),f(Yn(t),r(t))〉H+〈Yn(t),(e(t-⌊t⌋)A-1)f(Yn(t),r(t))〉H+〈Yn(t),e(t-⌊t⌋)A(f(Yn(⌊t⌋),r(⌊t⌋))Yn(t)(t)(t),e(t-⌊t⌋)Ae(t-⌊t⌋)A-f(Yn(t),r(t)))e(t-⌊t⌋)A〉H.
By Höld inequality, we have
(33)∥g(Yn(⌊t⌋),r(⌊t⌋))∥HS2=∥(g(Yn(t),r(t))-g(Yn(⌊t⌋),r(⌊t⌋)))g(Yn(t),r(t))-(g(Yn(t),r(t))-g(Yn(⌊t⌋),r(⌊t⌋)))∥HS2≤(1+Δ1/2)∥g(Yn(t),r(t))∥HS2+(1+Δ-1/2)×∥(g(Yn(t),r(t))-g(Yn(⌊t⌋),r(⌊t⌋)))∥HS2.
Then, from (31), we have
(34)eθt𝔼(λr(t)∥Yn(t)∥H2)≤q∥x∥H2+θq𝔼∫0teθs∥Yn(s)∥H2ds-2αp𝔼∫0teθs∥Yn(s)∥H2ds+𝔼∫0teθsλr(s){∥g(Yn(s),r(s))∥HS22〈Yn(s),f(Yn(s),r(s))〉H++𝔼∫0teθsλr(s)+∥g(Yn(s),r(s))∥HS2}ds+𝔼∫0teθs∑l=1Nγr(s)lλl∥Yn(s)∥H2ds+𝔼∫0teθsλr(s)×{∥g(Yn(s),r(s))-g(Yn(⌊s⌋),r(⌊s⌋)))∥HS22〈Yn(s),(e(s-⌊s⌋)A-1)f(Yn(s),r(s))〉H+2〈Yn(s),e(s-⌊s⌋)A×(f(Yn(⌊s⌋),r(⌊s⌋))-f(Yn(s),r(s)))e(s-⌊s⌋)A〉H+Δ1/2∥g(Yn(s),r(s))∥HS2+(1+Δ-1/2)×∥(g(Yn(s),r(s))-g(Yn(⌊s⌋),r(⌊s⌋)))∥HS2}ds≤q∥x∥H2+θq𝔼∫0teθs∥Yn(s)∥H2ds-2αp𝔼∫0teθs∥Yn(s)∥H2ds-μ2𝔼∫0teθs∥Y(s)∥H2ds+α1∫0teθsds+𝔼∫0teθsλr(s)×{∥g(Yn(s),r(s))-g(Yn(⌊s⌋),r(⌊s⌋)))∥HS22〈Yn(s),(e(s-⌊s⌋)A-1)f(Yn(s),r(s))〉H×+2〈Yn(s),e(s-⌊s⌋)A(f(Yn(⌊s⌋),r(⌊s⌋))Yn(s)Yn(s),e(s-⌊s⌋)AYn(s),e(s-⌊s⌋)A-f(Yn(s),r(s)))e(s-⌊s⌋)A〉H+Δ1/2∥g(Yn(s),r(s))∥HS2+(1+Δ-1/2)×∥(g(Yn(s),r(s))-g(Yn(⌊s⌋),r(⌊s⌋)))∥HS2}ds:=J1(t)+J2(t)+J3(t)+J4(t).
By the elemental inequality: 2ab≤(a2/κ)+κb2, a,b∈ℝ, κ>0, and (8), (27), we obtain that, for Δ<1,
(35)J2(t)≤𝔼∫0teθsλr(s){f(Yn(s),r(s))∥H2Δ1/2∥Yn(s)∥H2≤𝔼+∫0teθsλr(s)+Δ-1/2∥(e(s-⌊s⌋)A-1)≤𝔼+∫0teθsλr(s)×f(Yn(s),r(s))∥H2}ds≤𝔼∫0teθsλr(s)Δ1/2∥Yn(s)∥H2ds+𝔼∫0teθsλr(s)Δ-1/2ρn2Δ2L-(1+∥Yn(s)∥H2)ds≤𝔼∫0teθsλr(s){(Δ1/2+Δ1/2ρn2L-)≤𝔼+∫0teθsλr(s)×∥Yn(s)∥H2+Δ1/2ρn2L-}ds≤q𝔼∫0teθs{Δ1/2(1+ρn2L-)∥Yn(s)∥H2+Δ1/2ρn2L-}ds.
By (A2) and (8), we have
(36)∥(f(Yn(⌊t⌋),r(⌊t⌋))-f(Yn(t),r(t)))∥H2≤2∥(f(Yn(⌊t⌋),r(⌊t⌋))-f(Yn(⌊t⌋),r(t)))∥H2+2∥(f(Yn(⌊t⌋),r(t))-f(Yn(t),r(t)))∥H2≤8L-(1+∥Yn(⌊t⌋)∥H2)I{r(t)≠r(⌊t⌋)}+2L∥Yn(t)-Yn(⌊t⌋)∥H2.
Similarly, we have(37)∥g(Yn(t),r(t)))-g(Yn(⌊t⌋),r(⌊t⌋))∥HS2≤8L-(1+∥Yn(⌊t⌋)∥H2)I{r(t)≠r(⌊t⌋)}+2L∥Yn(t)-Yn(⌊t⌋)∥H2.
Thus, we obtain from (36) that
(38)J3(t)≤2𝔼∫0teθsλr(s)〈Yn(s),e(s-⌊s⌋)A+≤2𝔼∫0teθsλr(s)×(f(Yn(⌊s⌋),r(⌊s⌋))≤+2𝔼2𝔼∫0teθsλr(s)-f(Yn(s),r(s)))e(s-⌊s⌋)A〉Hds≤𝔼∫0teθsλr(s){Δ1/2∥Yn(s)∥H2ds+Δ-1/2∥e(s-⌊s⌋)A∥2≤2𝔼∫0teθsλr(s)×∥f(Yn(⌊s⌋),r(⌊s⌋))≤2𝔼2𝔼∫0teθsλr(s)-f(Yn(s),r(s))∥H2∥e(s-⌊s⌋)A∥2}ds≤q𝔼∫0teθs{Δ1/2∥Yn(s)∥H2+Δ-1/28L-+≤q𝔼∫0teθs×(1+∥Yn(⌊s⌋)∥H2)I{r(s)≠r(⌊s⌋)}+≤q𝔼∫0teθs+2Δ-1/2L∥Yn(s)-Yn(⌊s⌋)∥H2}ds.
By Markov property, we compute
(39)𝔼[(1+∥Y(⌊t⌋)∥H2)I{r(t)≠r(⌊t⌋)}]=𝔼(𝔼[(1+∥Y(⌊t⌋)∥H2)I{r(t)≠r(⌊t⌋)}∣r(⌊t⌋)])=𝔼(𝔼[(1+∥Y(⌊t⌋)∥H2)∣r(⌊t⌋)])×𝔼[I{r(t)≠r(⌊t⌋)}∣r(⌊t⌋)]=𝔼(1+∥Y(⌊t⌋)∥H2)∑i∈𝕊I{r(⌊t⌋)=i}ℙ(r(t)≠i∣r(⌊t⌋)=i)=𝔼(1+∥Y(⌊t⌋)∥H2)∑i∈𝕊I{r(⌊t⌋)=i}×∑j≠i(γij(t-⌊t⌋)+o(t-⌊t⌋))=𝔼(1+∥Y(⌊t⌋)∥H2)(maxi∈𝕊(-γii)Δ+o(Δ))∑i∈𝕊I{r(⌊t⌋)=i}≤γ^Δ𝔼(1+∥Y(⌊t⌋)∥H2),
where γ^=N[1+max1≤i≤N(-γii)]. Substituting (39) into(38) gives
(40)J3(t)≤q𝔼∫0teθs{Δ1/2∥Yn(s)∥H2≤q𝔼+∫0teθs+2Δ-1/2L∥Yn(s)-Yn(⌊s⌋)∥H2}ds+q∫0t8eθsΔ1/2γ^L-𝔼(1+∥Y(⌊s⌋)∥H2)ds.
Furthermore, due to (37) and (39), we have
(41)J4(t)=𝔼∫0teθsλr(s)×{Δ1/2∥g(Yn(s),r(s))∥HS2+(1+Δ-1/2)∥g(Yn(s),r(s)))-g(Yn(⌊s⌋),r(⌊s⌋))g(Yn(s),r(s)))∥HS2}ds≤q𝔼∫0teθsΔ1/2L-(1+∥Yn(s)∥H2)ds+q(1+Δ-1/2)×∫0teθs8L-(1+∥Yn(⌊s⌋)∥H2)I{r(s)≠r(⌊s⌋)}ds+2Lq(1+Δ-1/2)∫0teθs∥Yn(s)-Yn(⌊s⌋)∥H2ds≤qΔ1/2L-𝔼∫0teθs(1+∥Yn(s)∥H2)ds+16qγ^Δ1/2L-∫0teθs(1+∥Yn(s)∥H2)ds+2Lq(1+Δ-1/2)∫0teθs∥Yn(s)-Yn(⌊s⌋)∥H2ds.
On the other hand, by Lemma 10, when 3(ρn2+2L-)Δ≤1, we have
(42)𝔼∥Yn(t)∥H2≤2𝔼∥Yn(t)-Yn(⌊t⌋)∥H2+2𝔼∥Yn(⌊t⌋)∥H2≤6(ρn2+2L-)Δ(1+𝔼∥Yn(⌊t⌋)∥H2)+2𝔼∥Yn(⌊t⌋)∥H2≤4𝔼∥Yn(⌊t⌋)∥H2+2.
Putting (35), (40), and (41) into (34), we have
(43)eθt𝔼(λr(t)∥Yn(t)∥H2)≤q∥x∥H2+∫0teθs[α1+qρn2Δ1/2L-≤q∥x∥H2+∫0teθs+24qΔ1/2γ^L-+qΔ1/2L-]ds+𝔼∫0teθs[qθ-2αp-μ2+qΔ1/2(2+ρn2L-)+qΔ1/2L-]×∥Yn(s)∥H2ds+24qΔ1/2γ^L-𝔼×∫0teθs∥Yn(⌊s⌋)∥H2ds+(4qΔ-1/2L+2qL)𝔼∫0teθs∥Yn(s)-Yn(⌊s⌋)∥H2ds.
By Lemma 10 and the inequality (42), we obtain that
(44)eθt𝔼(λr(t)∥Yn(t)∥H2)≤q∥x∥H2+∫0teθs[Δ1/2α1+2qθ-4αp-μ≤q∥x∥H2++∫0teθs+3qρn2Δ1/2L-+24qΔ1/2γ^L-≤q∥x∥H2++∫0teθs+3qΔ1/2L-+4qΔ1/2]ds+∫0teθs[4qθ-8αp-2μ+4qΔ1/2(2+ρn2L-)++∫0teθs+4qΔ1/2L-+24qΔ1/2γ^L-]∥Yn(⌊s⌋)∥H2ds+6qLΔ1/2(ρn2+2L-)𝔼∫0teθs(1+∥Yn(⌊s⌋)∥H2)ds≤q∥x∥H2+∫0teθs[α1+2qθ-4αp-μ+3qρn2Δ1/2L-≤q∥x∥H2++∫0teθs+24qΔ1/2γ^L-+3qΔ1/2L-+4qΔ1/2≤q∥x∥H2++∫0teθs+6qLΔ1/2(ρn2+2L-)]ds+∫0teθs[4qθ-8αp-2μ+4qΔ1/2(2+ρn2L-)+∫0t+eθs+4qΔ1/2L-+24qΔ1/2γ^L-+∫0t+eθs+6qLΔ1/2(ρn2+2L-)]∥Yn(⌊s⌋)∥H2ds.
Let θ=(4αp+μ)/4q, for Δ<((4αp+μ)/(8q+4qρn2L-+4qL-+24qr^+6qL(ρn2+2L-)))2, then
(45)peθt𝔼(∥Y(t)∥H2)≤q∥x∥H2+∫0teθs[α1+4αp+μ2]ds.
That is,
(46)supt≥0𝔼(∥Y(t)∥H2)≤C.
Lemma 12.
Let (A1)–(A3) hold. If Δ<min{1,1/18(ρn2+2L),((2αp+μ)/(4q+2qL+2qρn2L+12qLr^))2}, then
(47)limt→∞𝔼∥Yn,x,i(t)-Yn,y,i(t)∥H2=0uniformlyforx,y∈U,
where U is a bounded subset of Hn.
Proof.
Write Yn,x,i(t)=Yx(t), Yn,y,i(t)=Yy(t), ri(kΔ)=r(kΔ). From (20), it is easy to show that
(48)(Yx(t)-Yy(t))-(Yx(⌊t⌋)-Yy(⌊t⌋))=(Yx(t)-Yx(⌊t⌋))-(Yy(t)-Yy(⌊t⌋))=(e(t-⌊t⌋)A-1)(Yx(⌊t⌋)-Yy(⌊t⌋))+∫⌊t⌋te(t-⌊s⌋)A(fn(Yx(⌊s⌋),r(⌊s⌋))+∫⌊t⌋te(t-⌊s⌋)A+-fn(Yy(⌊s⌋),r(⌊s⌋)))ds+∫⌊t⌋te(t-⌊s⌋)A((Yy(⌊s⌋),r(⌊s⌋))gn(Yx(⌊s⌋),r(⌊s⌋))+∫⌊t⌋te(t-⌊s⌋)A+-gn(Yy(⌊s⌋),r(⌊s⌋)))dW(s).
By using the argument of Lemma 10, we derive that, if Δ<1,
(49)𝔼∥(Yx(t)-Yy(t))-(Yx(⌊t⌋)-Yy(⌊t⌋))∥H2≤3(ρn2+2L)Δ𝔼∥Yx(⌊t⌋)-Yy(⌊t⌋)∥H2,(50)𝔼∥(Yx(t)-Yy(t))∥H2=𝔼∥(Yx(t)-Yy(t))-(Yx(⌊t⌋)-Yy(⌊t⌋))+(Yx(⌊t⌋)-Yy(⌊t⌋))∥H2≤(1+2)𝔼∥(Yx(t)-Yy(t))+≤(1+2)𝔼-(Yx(⌊t⌋)-Yy(⌊t⌋))∥H2+(1+12)𝔼∥Yx(⌊t⌋)-Yy(⌊t⌋)∥H2≤9(ρn2+2L)Δ𝔼∥(Yx(⌊t⌋)-Yy(⌊t⌋))∥H2+1.5𝔼∥(Yx(⌊t⌋)-Yy(⌊t⌋))∥H2.
If Δ<1/18(ρn2+2L), then(51)𝔼∥(Yx(t)-Yy(t))∥H2≤2𝔼∥(Yx(⌊t⌋)-Yy(⌊t⌋))∥H2.
Using (30) and the generalised Itô formula, for any θ>0, we have
(52)eθt𝔼(λr(t)∥Yx(t)-Yy(t)∥H2)≤λi∥x-y∥H2+𝔼∫0teθsλr(s){∥gn(Yx(⌊s⌋),r(⌊s⌋))-gn(Yy(⌊s⌋),r(⌊s⌋))∥HS2θ∥Yx(s)-Yy(s)∥H2𝔼+𝔼∫0teθsλr(s)+2〈Yx(s)-Yy(s),eθs𝔼+𝔼∫0teθsλr(s)AYx(s)-Yy(s)〉H𝔼+𝔼∫0teθsλr(s)+2〈(fn(Yx(⌊s⌋),r(⌊s⌋))-fn(Yy(⌊s⌋),r(⌊s⌋)))Yx(s)-Yy(s),e(s-⌊s⌋)A++𝔼∫0teθsλr(s)+2×((Yy(⌊s⌋),r(⌊s⌋))fn(Yx(⌊s⌋),r(⌊s⌋))++𝔼+𝔼∫0teθsλr(s)+2-fn(Yy(⌊s⌋),r(⌊s⌋)))〉H𝔼+𝔼∫0teθsλr(s)+∥gn(Yx(⌊s⌋),r(⌊s⌋))+𝔼+𝔼∫0teθsλ+r(s)-gn(Yy(⌊s⌋),r(⌊s⌋))∥HS2}ds+𝔼∫0teθs∑l=1Nγr(s)lλl∥Yx(s)-Yy(s)∥H2ds≤q∥x-y∥H2+qθ𝔼∫0teθs∥Yx(s)-Yy(s)∥H2ds-2αp𝔼∫0teθs∥Yx(s)-Yy(s)∥H2ds+𝔼∫0teθsλr(s){∥g(Yx(⌊s⌋),r(⌊s⌋))-g(Yy(⌊s⌋),r(⌊s⌋))∥HS22〈Yx(s)-Yy(s),e(s-⌊s⌋)A++𝔼∫0teθsλr(s)×(f(Yx(⌊s⌋),r(⌊s⌋))++𝔼∫0teθsλr(s)-f(Yy(⌊s⌋),r(⌊s⌋)))e(s-⌊s⌋)A〉H++𝔼∫0teθsλr(s)+∥g(Yx(⌊s⌋),r(⌊s⌋))+++𝔼∫0teθsλr(s)+-g(Yy(⌊s⌋),r(⌊s⌋))∥HS2}ds+𝔼∫0teθs∑l=1Nγr(s)lλl∥Yx(s)-Yy(s)∥H2ds.
By the fundamental transformation, we obtain that
(53)〈Yx(t)-Yy(t),e(t-⌊t⌋)A×(f(Yx(⌊t⌋),r(⌊t⌋))-f(Yy(⌊t⌋),r(⌊t⌋)))e(t-⌊t⌋)A〉H=〈Yx(t)-Yy(t),f(Yx(t),r(t))-f(Yy(t),r(t))〉H+〈Yx(t)-Yy(t),(e(t-⌊t⌋)A-1)×(f(Yx(t),r(t))-f(Yy(t),r(t)))e(t-⌊t⌋)A〉H+〈Yx(t)-Yy(t),e(t-⌊t⌋)A×(f(Yx(⌊t⌋),r(⌊t⌋))-f(Yy(⌊t⌋),r(⌊t⌋)))r++(⌊t⌋)-(f(Yx(t),r(t))-f(Yy(t),r(t)))e(t-⌊t⌋)A〉H.
By the Höld inequality, we have
(54)∥g(Yx(⌊t⌋),r(⌊t⌋))-g(Yy(⌊t⌋),r(⌊t⌋))∥HS2≤(1+Δ1/2)∥g(Yx(t),r(t))-g(Yy(t),r(t))∥HS2+(1+Δ-1/2)∥(g(Yx(t),r(t))-g(Yy(t),r(t)))-(g(Yx(⌊t⌋),r(⌊t⌋))-g(Yy(⌊t⌋),r(⌊t⌋)))∥HS2.
Then, from (52) and (A3), we have
(55)eθt𝔼(λr(t)∥Yx(t)-Yy(t)∥H2)≤q∥x-y∥H2+(qθ-2αp-μ)𝔼×∫0teθs∥Yx(s)-Yy(s)∥H2ds+2𝔼∫0teθsλr(s)×〈Yx(s)-Yy(s),(e(s-⌊s⌋)A-1)×(f(Yx(s),r(s))-f(Yy(s),r(s)))e(s-⌊s⌋)A〉Hds+2𝔼∫0teθsλr(s)×〈(f(Yx(s),r(s))-f(Yy(s),r(s)))Yx(s)-Yy(s),e(s-⌊s⌋)A×(f(Yx(⌊s⌋),r(⌊s⌋))-f(Yy(⌊s⌋),r(⌊s⌋)))-(f(Yx(s),r(s))-f(Yy(s),r(s)))〉Hds+𝔼∫0teθsλr(s){,r(⌊s⌋))Yy)∥HS2Δ1/2∥g(Yx(s),r(s))+𝔼++𝔼𝔼∫0teθsλr(s)-g(Yy(s),r(s))∥HS2++𝔼∫0teθsλr(s)+(1+Δ-1/2)++𝔼∫0teθsλr(s)×∥(g(Yx(s),r(s))+++𝔼𝔼𝔼𝔼∫0teθsλr(s)-g(Yy(s),r(s)))++𝔼𝔼∫0teθsλr(s)-(g(Yx(⌊s⌋),r(⌊s⌋))++𝔼𝔼+𝔼∫0teθsλr(s)-g(Yy(⌊s⌋),𝔼𝔼++++𝔼∫0teθsλr(s)λr(s)r(⌊s⌋)Yy)Yy)∥HS2}ds:=G1(t)+G2(t)+G3(t)+G4(t).
By (A2) and (27), we have, for Δ<1,
(56)G2(t)≤𝔼∫0teθsλr(s){(f(Yx(s),r(s))-Yy(s),r(s))Δ1/2∥Yx(s)-Yy(s)∥H2≤𝔼∫0teθsλr(s)+Δ-1/2∥(e(s-⌊s⌋)A-1)≤𝔼∫0teθsλr(s)+Δ-1/2×(f(Yx(s),r(s))≤𝔼∫0teθsλr(s)+Δ-1/2+-Yy(s),r(s))∥H2}ds≤𝔼∫0teθsλr(s)(Δ1/2+Δ3/2ρn2L)∥Yx(s)-Yy(s)∥H2ds≤q𝔼∫0teθsΔ1/2(1+ρn2L)∥Yx(s)-Yy(s)∥H2ds.
It is easy to show that
(57)G3(t)≤q𝔼∫0teθs{-(f(Yx(s),r(s))-Yy(s),r(s))∥H2Δ1/2∥Yx(s)-Yy(s)∥H2≤q𝔼∫0teθs+Δ-1/2∥(e(s-⌊s⌋)A)∥2≤q𝔼∫0teθs×∥f(Yx(⌊s⌋),r(⌊s⌋))≤q𝔼∫0teθs+-f(Yy(⌊s⌋),r(⌊s⌋))≤q𝔼∫0teθs+-(f(Yx(s),r(s))≤q𝔼∫0teθs++-Yy(s),r(s))∥H2}ds≤q𝔼∫0teθsΔ1/2∥Yx(s)-Yy(s)∥H2ds+G-3(t).
By (39), we have
(58)G-3(t)≤2qΔ-1/2𝔼∫0teθs[∥(f(Yx(s),r(s))-Yy(s),r(s))∥H2∥f(Yx(⌊s⌋),r(⌊s⌋))≤2qΔ-1/2𝔼∫0teθs-f(Yy(⌊s⌋),r(⌊s⌋))∥H2≤2qΔ-1/2𝔼∫0teθs+∥(f(Yx(s),r(s))≤2qΔ-1/2𝔼∫0teθs+-Yy(s),r(s))∥H2]×I{r(s)≠r(⌊s⌋)}ds≤4qΔ-1/2L𝔼∫0teθs∥Yx(⌊s⌋)-Yy(⌊s⌋)∥H2×I{r(s)≠r(⌊s⌋)}ds≤4qLγ^Δ1/2𝔼∫0teθs∥Yx(⌊s⌋)-Yy(⌊s⌋)∥H2ds.
Therefore, we obtain that
(59)G3(t)≤qΔ1/2𝔼∫0teθs∥Yx(s)-Yy(s)∥H2ds+4qLγ^Δ1/2𝔼∫0teθs∥Yx(⌊s⌋)-Yy(⌊s⌋)∥H2ds.
On the other hand, using the similar argument of (58), we have
(60)G4(t)≤q𝔼∫0teθs{Δ1/2L∥Yx(s)-Yy(s)∥H2G4(t)≤q𝔼∫0teθs+(1+Δ-1/2)4Lγ^ΔG4(t)≤q𝔼∫0teθs×∥Yx(⌊s⌋)-Yy(⌊s⌋)∥H2}ds.
Hence, we have
(61)peθt𝔼(∥Yx(t)-Yy(t)∥H2)≤q∥x-y∥H2+𝔼∫0teθs[qθ-2αp-μ+q(1+ρn2L)Δ1/2+qΔ1/2+qLΔ1/2]∥Yx(s)-Yy(s)∥H2ds+𝔼∫0teθs[4qLγ^Δ1/2+(Δ1/2+Δ)4qLγ^]×∥Yx(⌊s⌋)-Yy(⌊s⌋)∥H2ds.
By (50), we obtain that
(62)peθt𝔼(∥Yx(t)-Yy(t)∥H2)≤q∥x-y∥H2+𝔼∫0teθs[Δ1/22qθ-4αp-2μ+≤q∥x-y∥H2+𝔼∫0teθs+4qΔ1/2+2qLΔ1/2+≤q∥x-y∥H2+𝔼∫0teθs+2qρn2LΔ1/2+12qLγ^Δ1/2]×∥Yx(⌊s⌋)-Yy(⌊s⌋)∥H2ds.
Let θ=(2αp+μ)/2q, for Δ<((2αp+μ)/(4q+2qL+2qρn2L+12qLr^))2, then the desired assertion (47) follows.
We can now easily prove our main result.
Proof of Theorem 9.
Since Hn is finite-dimensional, by Lemma 3.1 in [12], we have
(63)limk→∞d𝕃(ℙkn,Δ((x,i),·×·),ℙkn,Δ((y,i),·×·))=0,
uniformly in x,y∈Hn,i,j∈𝕊.
By Lemma 7, there exists πn,Δ(·×·)∈𝒫(Hn×𝕊), such that
(64)limk→∞d𝕃(ℙkn,Δ((0,1),·×·),πn,Δ(·×·))=0.
By the triangle inequality (63) and (64), we have
(65)limk→∞d𝕃(ℙkn,Δ((x,i),·×·),πn,Δ(·×·))≤limk→∞d𝕃(ℙkn,Δ((x,i),·×·),ℙkn,Δ((0,1),·×·))+limk→∞d𝕃(ℙkn,Δ((0,1),·×·),πn,Δ(·×·))=0.
4. Corollary and Example
In this section, we give a criterion based M-matrices which can be verified easily in applications.
For each j∈𝕊, there exists a pair of constants βj and δj such that, for x,y∈H,
(66)〈x-y,f(x,j)-f(y,j)〉H≤βj∥x-y∥H2,∥g(x,j)-g(y,j)∥HS2≤δj∥x-y∥H2.
Moreover, 𝒜:=-diag(2β1+δ1,…,2βN+δN)-Γ is a nonsingular M-matrix [8].
Corollary 13.
Under (A1), (A2), and (A4), for a given stepsize Δ>0, and arbitrary x∈Hn,i∈𝕊, {Z-n,(x,i)(kΔ)}k≥0 has a unique stationary distribution πn,Δ(·×·)∈𝒫(Hn×𝕊).
Proof.
In fact, we only need to prove that (A3) holds. By (A4), there exists (λ1,λ2,…,λN)T>0, such that (q1,q2,…,qN)T=𝒜(λ1,λ2,…,λN)T>0.
Set μ=min1≤j≤Nqj, by (66), we have
(67)2λj〈x-y,f(x,j)-f(y,j)〉H+λj∥g(x,j)-g(y,j)∥H2+∑l=1Nγjlλl∥x-y∥H2≤2λjβj∥x-y∥H2+δjλj∥x-y∥H2+∑l=1Nγjlλl∥x-y∥H2=(2λjβj+δjλj+∑l=1Nγjlλl)∥x-y∥H2=-qj∥x-y∥H2≤-μ∥x-y∥H2
In the following, we give an example to illustrate the Corollary 13.
We take H=L2(0,π) and A=∂2/∂ξ2 with domain 𝒟(A)=H2(0,π)∩H01(0,π), then A is a self-adjoint negative operator. For the eigenbasis ek(ξ)=(2/π)1/2sin(kξ), ξ∈[0,π], Aek=-k2ek, k∈ℕ. It is easy to show that
(69)∥etAx∥H2=∑i=1∞e-2k2t〈x,ei〉H2≤e-2t∑i=1∞〈x,ei〉H2.
This further gives that
(70)∥etA∥≤e-t,
where α=1, thus (A1) holds.
Let W(t) be a scalar Brownian motion, let r(t) be a continuous-time Markov chain values in 𝕊=1,2, with the generator
(71)Γ=(-221-1),B(1)=B1=(-0.3-0.1-0.2-0.2),B(2)=B2=(-0.4-0.2-0.3-0.2).
Then λmax(B1TB1)=0.1706, λmax(B2TB2)=0.3286.
Moreover, g satisfies
(72)∥g(x,j)-g(y,j)∥HS2≤δj∥x-y∥H2,
where δ1=0.1, δ2=0.06.
Defining f(x,j)=B(j)x, then
(73)∥f(x,j)-f(y,j)∥HS2∨∥g(x,j)-g(y,j)∥HS2≤(λmax(BjTBj)∨δj)∥x-y∥H2<0.33∥x-y∥H2,〈x-y,f(x,j)-f(y,j)〉H≤12〈x-y,(BjT+Bj)(x-y)〉H≤12λmax(BjT+Bj)∥x-y∥H2.
It is easy to compute
(74)β1=12λmax(B1T+B1)=-0.0919,β2=12λmax(B2T+B2)=-0.03075.
So the matrix 𝒜 becomes
(75)𝒜=diag(0.0838,0.0015)-Γ=(2.0838-2-11.0015).
It is easy to see that 𝒜 is a nonsingular M-matrix. Thus, (A4) holds. By Corollary 13, we can conclude that (68) has a unique stationary distribution πn,Δ(·×·).
Acknowledgment
This work is partially supported by the National Natural Science Foundation of China under Grants 61134012 and 11271146.
Da PratoG.ZabczykJ.1992Cambridge, UKCambridge University Press10.1017/CBO9780511666223MR1207136DebusscheA.Weak approximation of stochastic partial differential equations: the nonlinear case2011802738911710.1090/S0025-5718-2010-02395-6MR2728973ZBL1217.65011GyöngyI.MilletA.Rate of convergence of space time approximations for stochastic evolution equations2009301296410.1007/s11118-008-9105-5MR2465711ZBL1168.60025JentzenA.KloedenP. E.2011CBMS-NSF Regional Conference Series in Applied Mathematics10.1137/1.9781611972016MR2856611ShardlowT.Numerical methods for stochastic parabolic PDEs1999201-212114510.1080/01630569908816884MR1683281ZBL0919.65100AthansM.Command and control (C2) theory: a challenge to control science19873242862932-s2.0-0023324578HighamD. J.MaoX.YuanC.Preserving exponential mean-square stability in the simulation of hybrid stochastic differential equations2007108229532510.1007/s00211-007-0113-yMR2358006ZBL1137.65007MaoX.YuanC.2006London, UKImperial CollegeMR2256095MaritonM.1990Marcel DekkerBaoJ.HouZ.YuanC.Stability in distribution of mild solutions to stochastic partial differential equations201013862169218010.1090/S0002-9939-10-10230-5MR2596056ZBL1195.60087BaoJ.YuanC.Numerical approximation of stationary distribution for SPDEshttp://arxiv.org/abs/1303.1600MaoX.YuanC.YinG.Numerical method for stationary distribution of stochastic differential equations with Markovian switching2005174112710.1016/j.cam.2004.03.016MR2102646ZBL1066.65016YuanC.MaoX.Stationary distributions of Euler-Maruyama-type stochastic difference equations with Markovian switching and their convergence2005111294810.1080/10236190412331314150MR2112804ZBL1076.65012YuanC.MaoX.Asymptotic stability in distribution of stochastic differential equations with Markovian switching2003103227729110.1016/S0304-4149(02)00230-2MR1950767ZBL1075.60541AndersonW. J.1991New York, NY, USASpringer10.1007/978-1-4612-3038-0MR1118840ZBL0926.91029