# Numerical Study on Stochastic Diabetes Mellitus Model with Additive Noise

Authors: Zhifang Zhang; Qingyi Zhan; Xiangdong Xie
Journal: Computational and Mathematical Methods in Medicine
 (2019)

---

## Abstract

This article focuses on the numerical analysis and simulation of the stochastic diabetes mellitus model with additive noise. The existence and uniqueness theorem of the solution under some appropriate assumptions is established. And, the mean square stability and convergence of numerical solutions are proposed, too. The practical use of these theorems is demonstrated in the numerical computations of the stochastic diabetes mellitus model and the value for the forecast of the tendency of diabetes mellitus in a given time.


---

## Body

## 1. Introduction

At present, with the development of the society and the increment of the economy, diabetes mellitus is becoming more and more popular in the world. In fact, diabetes mellitus is the name given to a group of different conditions in which there is not the right amount of insulin to stabilize the amount of sugar in the body. As we know, there exist two forms of diabetes mellitus. Type I diabetes mellitus, which depends on insulin, most often occurs in young people, while type II diabetes mellitus, which does not depend on insulin, usually develops in the aged.It is well known that there are many finished works which utilize two methods to investigate the diabetes mellitus in the view of mathematical model. One pays attention to microscopic action of the nosogenesis of diabetes mellitus and forms many mathematical models such as ordinary differential equations and partial differential equation (refer [1] and the references therein). Another takes account of the macroscopic case of the size of population of diabetes mellitus in a given time. This leads to insight into the control of diabetes mellitus, and with the current alarming increase in the incidence of the disease of a given considered region, this area has gained increased interest and importance.However, due to many uncertainties or random influences, where these uncertainties come from the influence of diet, physical activity level, and the age dynamical distribution of population, noises should be taken into account. The existing deterministic mathematical models of diabetes mellitus [2] need to be revised so that it can simulate the fact more really. Therefore, we expand it to the case of stochastic differential equations (SDEs), whose applications describe many natural phenomena in meteorology, biology, and so on [3, 4]. As far as we know, till now, there has been little investigation of the diabetes mellitus mathematical model in the view of SDEs in the literature. Stochastic numerical analysis is still an interesting method of studying epidemic disease tendency of diabetes mellitus.The main motivations of this work are twofold. On one side of the coin, the classical results about the deterministic mathematical model of diabetes mellitus are the base of this research. A variety of mathematical models have been used for different aspects of diabetes mellitus, and many important results, which can reveal the facts of diabetes mellitus, are obtained (refer [2, 5–9] and the references therein). On the other hand, it has been attracted by some random phenomena which often appear in the population dynamics of diabetes mellitus. These need numerical simulations to direct the control policy of diabetes mellitus. Furthermore, it is the fact that there exists our earlier work [10, 11] on stability analysis and numerical simulation of SDEs. For example, there are some work on numerical analysis of SDEs [11–15] and numerical simulation of SDEs [10]. These carry out the foundation of numerical analysis [16].In this work, we first prove the existence and uniqueness theorem of the solution of the stochastic diabetes mellitus model under some assumptions. Then, the mean square stability and convergence are proposed. And, numerical examples are shown to illustrate the possibility of the stochastic mathematical model of diabetes mellitus and the value in the forecast of the tendency of diabetes mellitus in a given time. These results show that, under some appropriate conditions, SDEs can simulate the epidemic disease tendency of diabetes mellitus more accurately, whose value can be well estimated by the numerical approximative solution.A more detailed outline of this paper is as follows. Section2 shows some relevant concepts and norms which will be utilized later. Section 3 is devoted to the theoretical analysis of the stochastic diabetes mellitus model, that is, the existence of the solution, mean square stability, and convergence. Section 4 presents numerical experiments of the stochastic diabetes mellitus model in some given areas. Illustrative numerical results for the main theorem are included. Section 5 is addressed to the conclusions of this article.


## 2. Preliminaries

### 2.1. Generated Differential Equation Model

Here, we consider a stochastic model which can describe the dynamic behaviour of diabetes mellitus.LetΩ,ℱ,ℙ be a canonical Wiener space, ℱtt∈ℝ+ be its natural normal filtration, and Wtt∈ℝ+ be a standard one-dimensional Brownian motion defined on the space Ω,ℱ,ℙ. We assume that Ω≔ω∈Cℝ+,ℝ:ω0=0 endowed with the compact-open topology. In the realization, Wtω≔ωt, where ω⋅∈Ω, which means that the elements of Ω can be identified with the paths of the Wiener process. Based on the conclusions of deterministic ordinary differential equations about diabetes mellitus in [1, 2], we consider a class of Itô SDEs in the form of(1)dCt=−λt+θtCt+λtNtdt+g1tdWt,dNt=It−νt+δtCt−μtNtdt+g2tdWt,where t∈ℝ+; the quantity of diabetes mellitus which has complications in a special research region at time t is written as CtCt∈ℝ+, and C0=C0; Nt is on behalf of the scale of the population which has diabetes mellitus in a special research region at time t, namely, Nt=Ct+Dt. Here, DtDt∈ℝ+ presents the quantity of diabetes mellitus which has no complications in a special research region at time t; the morbidity of diabetes mellitus in a special research region at time t is represented as It; μt stands for the mortality rate, the chance of a diabetes mellitus person who is developing a complication is written as λt, the proportion whose complications are mended is shown as γt, the parameter νt presents the rate at which diabetic patients with complication become severely disable, the parameter δt shows the mortality rate due to complications, θt≔γt+μt+νt+δt denotes the sum of the above parameters, and g1tandg2t are functions with respect to t, which denote the uncertain influences [2].


### 2.2. Basic Notations and Assumptions

We make use of the following notations which is similar to [11]:(i)
LetL2Ω,ℙ be the space of all square-integrable random variables x:Ω⟶ℝd.(ii)
The norm of a random variablex=x1,x2,…,xd∈L2Ω,ℙ is defined as(2)x2=∫Ωx1ω2+x2ω2+,…,+xdω2dℙω1/2<∞.(iii)
The norm of a stochastic processxt,ω is defined as xt,ω2=supt∈ℝ+xtω2<∞, where xtω∈L2Ω,ℙ and t∈ℝ+.(iv)
We define the norm of random matrixAL2Ω,ℙ=EA21/2, where A is a random matrix and ⋅ is the operator norm.(v)
Unless otherwise stated, the norms⋅2 and ⋅L2Ω,ℙ are usually denoted as ⋅ in sequels.In this paper, we also make the following assumptions which are used for the theoretical analysis [11].Hypothesis 2.1.
(i)
The initial valuesC0 and D0 are bounded; that is,(3)maxC0,D0≤K1,
forK1>0
(ii)
Assume that the functionλ:ℝ⟶ℝ is continuous, measurable function and the function θ:ℝ⟶ℝ is continuous, too
(iii)
The functionsλ and θ are globally bounded with respect to t. That is, there exists a positive constant J such that
(4)
max
t
θ
t
,
λ
t
≤
J
,
holds fort∈ℝ+
(iv)
The functionsg1t and g2t are globally bounded. That is, there exists a constant K2>0 such that
(5)
max
t
g
1
t
,
g
2
t
≤
K
2
,
holds fort∈ℝ+


### 2.3. Equivalent Form

SDE (1) can be rewritten in the matrix-vector form as follows:(6)dX=bX,tdt+BX,tdWt,where(7)Xt=CtNt,X0=C0N0,BX,t=g1tg2t,(8)bX,t=−λt+θtCt+λtNtIt−νt+δtCt−μtNt.We define(9)θ:−∞,+∞×Ω⟶Ω,θtωs=ωt+s−ωt.and Δ≔s,t∈ℝ2,s≤t. By the conclusions in [3], SDE (1) generates a stochastic flow φ:Δ×ℝ2×Ω⟶ℝ2 when the solution of SDE (1) exists uniquely, which is usually written as φs,t,x0,ω≔φs,t,ωx0 on the metric dynamical systems Ω,ℱ,ℙ,θt. The stochastic flow φ is given by(10)φs,t,ωX0=X0+∫stbφs,r,ωX0,rdr+∫stBφs,r,ωX0,rdWr,t≥s.


## 2.1. Generated Differential Equation Model

Here, we consider a stochastic model which can describe the dynamic behaviour of diabetes mellitus.LetΩ,ℱ,ℙ be a canonical Wiener space, ℱtt∈ℝ+ be its natural normal filtration, and Wtt∈ℝ+ be a standard one-dimensional Brownian motion defined on the space Ω,ℱ,ℙ. We assume that Ω≔ω∈Cℝ+,ℝ:ω0=0 endowed with the compact-open topology. In the realization, Wtω≔ωt, where ω⋅∈Ω, which means that the elements of Ω can be identified with the paths of the Wiener process. Based on the conclusions of deterministic ordinary differential equations about diabetes mellitus in [1, 2], we consider a class of Itô SDEs in the form of(1)dCt=−λt+θtCt+λtNtdt+g1tdWt,dNt=It−νt+δtCt−μtNtdt+g2tdWt,where t∈ℝ+; the quantity of diabetes mellitus which has complications in a special research region at time t is written as CtCt∈ℝ+, and C0=C0; Nt is on behalf of the scale of the population which has diabetes mellitus in a special research region at time t, namely, Nt=Ct+Dt. Here, DtDt∈ℝ+ presents the quantity of diabetes mellitus which has no complications in a special research region at time t; the morbidity of diabetes mellitus in a special research region at time t is represented as It; μt stands for the mortality rate, the chance of a diabetes mellitus person who is developing a complication is written as λt, the proportion whose complications are mended is shown as γt, the parameter νt presents the rate at which diabetic patients with complication become severely disable, the parameter δt shows the mortality rate due to complications, θt≔γt+μt+νt+δt denotes the sum of the above parameters, and g1tandg2t are functions with respect to t, which denote the uncertain influences [2].


## 2.2. Basic Notations and Assumptions

We make use of the following notations which is similar to [11]:(i)
LetL2Ω,ℙ be the space of all square-integrable random variables x:Ω⟶ℝd.(ii)
The norm of a random variablex=x1,x2,…,xd∈L2Ω,ℙ is defined as(2)x2=∫Ωx1ω2+x2ω2+,…,+xdω2dℙω1/2<∞.(iii)
The norm of a stochastic processxt,ω is defined as xt,ω2=supt∈ℝ+xtω2<∞, where xtω∈L2Ω,ℙ and t∈ℝ+.(iv)
We define the norm of random matrixAL2Ω,ℙ=EA21/2, where A is a random matrix and ⋅ is the operator norm.(v)
Unless otherwise stated, the norms⋅2 and ⋅L2Ω,ℙ are usually denoted as ⋅ in sequels.In this paper, we also make the following assumptions which are used for the theoretical analysis [11].Hypothesis 2.1.
(i)
The initial valuesC0 and D0 are bounded; that is,(3)maxC0,D0≤K1,
forK1>0
(ii)
Assume that the functionλ:ℝ⟶ℝ is continuous, measurable function and the function θ:ℝ⟶ℝ is continuous, too
(iii)
The functionsλ and θ are globally bounded with respect to t. That is, there exists a positive constant J such that
(4)
max
t
θ
t
,
λ
t
≤
J
,
holds fort∈ℝ+
(iv)
The functionsg1t and g2t are globally bounded. That is, there exists a constant K2>0 such that
(5)
max
t
g
1
t
,
g
2
t
≤
K
2
,
holds fort∈ℝ+


## 2.3. Equivalent Form

SDE (1) can be rewritten in the matrix-vector form as follows:(6)dX=bX,tdt+BX,tdWt,where(7)Xt=CtNt,X0=C0N0,BX,t=g1tg2t,(8)bX,t=−λt+θtCt+λtNtIt−νt+δtCt−μtNt.We define(9)θ:−∞,+∞×Ω⟶Ω,θtωs=ωt+s−ωt.and Δ≔s,t∈ℝ2,s≤t. By the conclusions in [3], SDE (1) generates a stochastic flow φ:Δ×ℝ2×Ω⟶ℝ2 when the solution of SDE (1) exists uniquely, which is usually written as φs,t,x0,ω≔φs,t,ωx0 on the metric dynamical systems Ω,ℱ,ℙ,θt. The stochastic flow φ is given by(10)φs,t,ωX0=X0+∫stbφs,r,ωX0,rdr+∫stBφs,r,ωX0,rdWr,t≥s.


## 3. Theoretical Results

### 3.1. Existence of Equation (1)’s Solutions

The following result guarantees the existence of solutions for SDEs and is a direct consequence of Theorem 3.2.4 in [17].Lemma 3.1.
Suppose thatb:ℝ2×0,T⟶ℝ2 and B:ℝ2×0,T⟶M2 are continuous and satisfy the following conditions for some constant L:(i)
bX,t−bX^,t≤LX−X^,BX,t−BX^,t≤LX−X^ for all 0≤t≤T and X,X^∈ℝ2(ii)
bX,t≤L1+X,BX,t≤L1+X for all 0≤t≤T and X∈ℝ2(iii)
LetX0 be any ℝ2-valued random variable such that EX02<∞ and X0 is independent of ℱ0
Then, there exists a unique solutionX∈L2Ω,ℙ of the stochastic differential equation:(11)dX=bX,tdt+BX,tdWt,0≤t≤T,X0=X0.
By the conclusions of Lemmas 3.1, we obtain the following theorem:Theorem 3.2.
Suppose that SDE (1) satisfies Hypothesis 2.1 and the initial conditions are given in Section 2.1.
Then, SDE (1) has a uniqueness solution Xt=Ct,Nt for all 0≤t≤T.Proof.
In order to utilize Lemma 3.1 to this problem, we only need to check that the conditions of this theorem satisfy its three hypotheses.
First and foremost, Hypothesis (iii) obviously holds.
Secondly, by the assumptions of SDE (1), we obtain that(12)bX,t−bX^,t=−λt+θtCt−C^t+λtNt−N^t−νt+δtCt−C^t−μtNt−N^t=−λt+θt,λt−νt+δt,−μtCt−C^tNt−N^t≤−λt+θt,λt−νt+δt,−μt⋅X−X^≤At⋅X−X^,where(13)At=−λt+θt,λt−νt+δt,−μt.
It follows from the definition of random matrix and Hypothesis 2.1(ii) that we obtain(14)bX,t−bX^,t≤L⋅X−X^,where(15)A≤maxt∈0,TEλt+θt2+λ2t+νt+δt2+μ2t1/2≤7J≔L.
By the similar way, we can prove that(16)BX,t−BX^,t=g1tCt−C^tg2tCt−C^t=g1tg2t⋅X−X^.
It follows from the definition of random matrix and Hypothesis 2.1 (iii) that we obtain(17)BX,t−BX^,t≤K2X−X^.
This completes the check of the first hypothesis. Last but not least, it follows from Hypothesis 2.1 (iv) that(18)bX,t=−λt+θtCt+λtNt−νt+δtCt−μtNt=−λt+θt,λt−νt+δt,−μtCtNt≤A⋅X≤LX≤L1+X.
By the same way, we can obtain that(19)BX,t=g1tg2t≤K2≤L≤L1+X.
This completes the check of the second hypothesis.
Therefore, the conclusion of Theorem 3.2 follows from Lemma 3.1. The proof is finished.


### 3.2. Mean-Square Asymptotical Stability

In this section, we investigate the mean-square uniformly asymptotic stability of the solutionφs,t,ωX0 of SDE (1). The pullback method is a powerful tool to the proof of uniformly asymptotic stability. To be precise, let us introduce some related definition [18].Definition 3.1.
The solutionφs,t,ωX0 of SDE (1) is said to be mean-square asymptotically stable if, for any given ϵ>0, every other solution φs,t,ωX^0 of SDE (1) satisfies(20)limt⟶+∞φs,t,ωX0−φs,t,ωX^0=0,for any bounded ℱs-measurable bounded initial values X0 and X^0, respectively, where X0−X^0<ε.Theorem 3.3.
Assume that for any initial valuesX0 and X^0∈L2Ω,ℙ, the coefficients of SDE (1) satisfy Theorem 3.2, then the solution φt−τ,t,θ−τωX0 of SDE (1) is mean-square asymptotically stable.Proof.
First and foremost, letφt−τ,t,θ−τωX^0 be another solution of SDE (1) and ϵ>0 be an arbitrary constant. If X0−X^0≤ϵ, it follows from (9) and the method which is used to estimate [10] that(21)Eφt−τ,t,θ−τωX0−φt−τ,t,θ−τωX^02≤I1+I2+I3,(22)I1=3EX0−X^02,I2=3E∫t−τtbφt−τ,t,θ−τωX0,r−bφt−τ,t,θ−τωX^0,rdr2,(23)I3=3E∫t−τtBφt−τ,t,θ−τωX0,r−Bφt−τ,t,θ−τωX^0,rdWr2.
It follows from the Cauchy–Schwarz inequality and the Lipschitz condition of the functionb that we have(24)I2≤3τ∫t−τtEbφt−τ,t,θ−τωX0,r−bφt−τ,t,θ−τωX^0,r2dr≤3τ⋅L2∫t−τtEφt−τ,t,θ−τωX0−φt−τ,t,θ−τωX^02dr.
The Itô isometry and the global Lipschitz condition of the functionB imply that(25)I3≤3∫t−τtEbφt−τ,t,θ−τωX0,r−Bφt−τ,t,θ−τωX^0,r2dr≤3L2∫t−τtEφt−τ,t,θ−τωX0−φt−τ,t,θ−τωX^02dr.
By the Gronwall inequality, there exists a numberM1 such that(26)φt−τ,t,θ−τωX0−φt−τ,t,θ−τωX^0≤M1,where(27)M1=3X0−X^03τ+1L2expτ.
Therefore, by the fact thatM1⟶0 as τ⟶−∞, we obtain that(28)limt⟶+∞φt−τ,t,θ−τωX0−φt−τ,t,θ−τωX^0=0.
Then, by Definition 3.1, it is mean-square asymptotically stable.
This completes the proof.


### 3.3. Mean-Square Convergence

The finite time interval0,t is divided into N subintervals with the length Δt. The exact solution of SDE (1) in 0,t has the form(29)Xt,ω=X0+∫0NΔtbX0,rdr+∫0NΔtBX0,rdWr.The Euler–Maruyama scheme is applied to SDE (6), and we have the following form(30)Xk+1=Xk+bXk,tkΔtk+BXk,tkΔWk,where tk=kΔt,k=0,1,2,…,N and ΔWk=Wk+1−Wk.The Milstein scheme is applied to SDE (6), and we have the following form(31)Xk+1=Xk+bXk,tkΔtk+BXk,tkΔWk+12BXk,tk∂BXk,tk∂XΔWk2−Δtk.The following result shows that the numerical approximationXk to the solution of SDE (1) is mean-square convergent to the exact solution of SDE (1) under some conditions.Theorem 3.4.
Assume that, for any initial valueX0∈L2Ω,ℙ, the coefficients of SDE (1) satisfy Theorem 3.3; then, the numerical approximation Xk to the solution of SDE (1) by Euler–Maruyama scheme and Milstein scheme is mean-square convergent, and the convergence order is 0.5.Proof.
We are interested in the mean square convergence to zero of the error(32)ek=Xk−Xtk,ω,where Xtk,ω denotes the theoretical solution of SDE (1) at the time tk. From the expression of Xt,ω, we obtain(33)Xtk+1,ω=Xtk,ω+∫kΔtk+1ΔtbXtk,ω,rdr+∫kΔtk+1ΔtBXtk,ω,tkdWr.
Then, it implies that(34)EXk+1−Xtk+1,ω2≤I4+I5+I6,
where(35)I4=3EXk−Xtk,ω2,I5=3EbXk,tkΔtk−∫kΔtk+1ΔtbXtk,ω,rdr2,(36)I6=3EBXk,tkΔwk−∫kΔtk+1ΔtBXtk,ω,rdWr2.
We notice from the Cauchy–Schwarz inequality and the global Lipschitz condition of functionb that we can obtain(37)I5≤3E∫kΔtk+1ΔtbXk,rdr−∫kΔtk+1ΔtbXtk,ω,rdr2≤3ΔtkL2∫kΔtk+1ΔtEXk−Xtk,ω2dr.
Then, the Itô isometry and the global Lipschitz condition of the functionB imply that(38)I6≤3E∫kΔtk+1ΔtBXk,rdWr−∫kΔtk+1ΔtBXtk,ω,rdWr2≤3L2∫kΔtk+1ΔtEXk−Xtk,ω2dr.
It follows from the Gronwall inequality that there exists a numberM2 such that(39)EXk+1−Xtk+1,ω2≤M2,where(40)M2=3EXk−Xtk,ω2⋅3Δtk+1L2expΔtk.
By the fact thatM2 tends to zero as Δtk⟶0, that is,(41)limΔtk⟶0M2=0,we obtain(42)limΔtk⟶0ek+1=limN⟶+∞Xk+1−Xtk+1,ω=0.
Therefore, it is mean-square convergent, and the convergence order is 0.5. We have established the theorem.


## 3.1. Existence of Equation (1)’s Solutions

The following result guarantees the existence of solutions for SDEs and is a direct consequence of Theorem 3.2.4 in [17].Lemma 3.1.
Suppose thatb:ℝ2×0,T⟶ℝ2 and B:ℝ2×0,T⟶M2 are continuous and satisfy the following conditions for some constant L:(i)
bX,t−bX^,t≤LX−X^,BX,t−BX^,t≤LX−X^ for all 0≤t≤T and X,X^∈ℝ2(ii)
bX,t≤L1+X,BX,t≤L1+X for all 0≤t≤T and X∈ℝ2(iii)
LetX0 be any ℝ2-valued random variable such that EX02<∞ and X0 is independent of ℱ0
Then, there exists a unique solutionX∈L2Ω,ℙ of the stochastic differential equation:(11)dX=bX,tdt+BX,tdWt,0≤t≤T,X0=X0.
By the conclusions of Lemmas 3.1, we obtain the following theorem:Theorem 3.2.
Suppose that SDE (1) satisfies Hypothesis 2.1 and the initial conditions are given in Section 2.1.
Then, SDE (1) has a uniqueness solution Xt=Ct,Nt for all 0≤t≤T.Proof.
In order to utilize Lemma 3.1 to this problem, we only need to check that the conditions of this theorem satisfy its three hypotheses.
First and foremost, Hypothesis (iii) obviously holds.
Secondly, by the assumptions of SDE (1), we obtain that(12)bX,t−bX^,t=−λt+θtCt−C^t+λtNt−N^t−νt+δtCt−C^t−μtNt−N^t=−λt+θt,λt−νt+δt,−μtCt−C^tNt−N^t≤−λt+θt,λt−νt+δt,−μt⋅X−X^≤At⋅X−X^,where(13)At=−λt+θt,λt−νt+δt,−μt.
It follows from the definition of random matrix and Hypothesis 2.1(ii) that we obtain(14)bX,t−bX^,t≤L⋅X−X^,where(15)A≤maxt∈0,TEλt+θt2+λ2t+νt+δt2+μ2t1/2≤7J≔L.
By the similar way, we can prove that(16)BX,t−BX^,t=g1tCt−C^tg2tCt−C^t=g1tg2t⋅X−X^.
It follows from the definition of random matrix and Hypothesis 2.1 (iii) that we obtain(17)BX,t−BX^,t≤K2X−X^.
This completes the check of the first hypothesis. Last but not least, it follows from Hypothesis 2.1 (iv) that(18)bX,t=−λt+θtCt+λtNt−νt+δtCt−μtNt=−λt+θt,λt−νt+δt,−μtCtNt≤A⋅X≤LX≤L1+X.
By the same way, we can obtain that(19)BX,t=g1tg2t≤K2≤L≤L1+X.
This completes the check of the second hypothesis.
Therefore, the conclusion of Theorem 3.2 follows from Lemma 3.1. The proof is finished.


## 3.2. Mean-Square Asymptotical Stability

In this section, we investigate the mean-square uniformly asymptotic stability of the solutionφs,t,ωX0 of SDE (1). The pullback method is a powerful tool to the proof of uniformly asymptotic stability. To be precise, let us introduce some related definition [18].Definition 3.1.
The solutionφs,t,ωX0 of SDE (1) is said to be mean-square asymptotically stable if, for any given ϵ>0, every other solution φs,t,ωX^0 of SDE (1) satisfies(20)limt⟶+∞φs,t,ωX0−φs,t,ωX^0=0,for any bounded ℱs-measurable bounded initial values X0 and X^0, respectively, where X0−X^0<ε.Theorem 3.3.
Assume that for any initial valuesX0 and X^0∈L2Ω,ℙ, the coefficients of SDE (1) satisfy Theorem 3.2, then the solution φt−τ,t,θ−τωX0 of SDE (1) is mean-square asymptotically stable.Proof.
First and foremost, letφt−τ,t,θ−τωX^0 be another solution of SDE (1) and ϵ>0 be an arbitrary constant. If X0−X^0≤ϵ, it follows from (9) and the method which is used to estimate [10] that(21)Eφt−τ,t,θ−τωX0−φt−τ,t,θ−τωX^02≤I1+I2+I3,(22)I1=3EX0−X^02,I2=3E∫t−τtbφt−τ,t,θ−τωX0,r−bφt−τ,t,θ−τωX^0,rdr2,(23)I3=3E∫t−τtBφt−τ,t,θ−τωX0,r−Bφt−τ,t,θ−τωX^0,rdWr2.
It follows from the Cauchy–Schwarz inequality and the Lipschitz condition of the functionb that we have(24)I2≤3τ∫t−τtEbφt−τ,t,θ−τωX0,r−bφt−τ,t,θ−τωX^0,r2dr≤3τ⋅L2∫t−τtEφt−τ,t,θ−τωX0−φt−τ,t,θ−τωX^02dr.
The Itô isometry and the global Lipschitz condition of the functionB imply that(25)I3≤3∫t−τtEbφt−τ,t,θ−τωX0,r−Bφt−τ,t,θ−τωX^0,r2dr≤3L2∫t−τtEφt−τ,t,θ−τωX0−φt−τ,t,θ−τωX^02dr.
By the Gronwall inequality, there exists a numberM1 such that(26)φt−τ,t,θ−τωX0−φt−τ,t,θ−τωX^0≤M1,where(27)M1=3X0−X^03τ+1L2expτ.
Therefore, by the fact thatM1⟶0 as τ⟶−∞, we obtain that(28)limt⟶+∞φt−τ,t,θ−τωX0−φt−τ,t,θ−τωX^0=0.
Then, by Definition 3.1, it is mean-square asymptotically stable.
This completes the proof.


## 3.3. Mean-Square Convergence

The finite time interval0,t is divided into N subintervals with the length Δt. The exact solution of SDE (1) in 0,t has the form(29)Xt,ω=X0+∫0NΔtbX0,rdr+∫0NΔtBX0,rdWr.The Euler–Maruyama scheme is applied to SDE (6), and we have the following form(30)Xk+1=Xk+bXk,tkΔtk+BXk,tkΔWk,where tk=kΔt,k=0,1,2,…,N and ΔWk=Wk+1−Wk.The Milstein scheme is applied to SDE (6), and we have the following form(31)Xk+1=Xk+bXk,tkΔtk+BXk,tkΔWk+12BXk,tk∂BXk,tk∂XΔWk2−Δtk.The following result shows that the numerical approximationXk to the solution of SDE (1) is mean-square convergent to the exact solution of SDE (1) under some conditions.Theorem 3.4.
Assume that, for any initial valueX0∈L2Ω,ℙ, the coefficients of SDE (1) satisfy Theorem 3.3; then, the numerical approximation Xk to the solution of SDE (1) by Euler–Maruyama scheme and Milstein scheme is mean-square convergent, and the convergence order is 0.5.Proof.
We are interested in the mean square convergence to zero of the error(32)ek=Xk−Xtk,ω,where Xtk,ω denotes the theoretical solution of SDE (1) at the time tk. From the expression of Xt,ω, we obtain(33)Xtk+1,ω=Xtk,ω+∫kΔtk+1ΔtbXtk,ω,rdr+∫kΔtk+1ΔtBXtk,ω,tkdWr.
Then, it implies that(34)EXk+1−Xtk+1,ω2≤I4+I5+I6,
where(35)I4=3EXk−Xtk,ω2,I5=3EbXk,tkΔtk−∫kΔtk+1ΔtbXtk,ω,rdr2,(36)I6=3EBXk,tkΔwk−∫kΔtk+1ΔtBXtk,ω,rdWr2.
We notice from the Cauchy–Schwarz inequality and the global Lipschitz condition of functionb that we can obtain(37)I5≤3E∫kΔtk+1ΔtbXk,rdr−∫kΔtk+1ΔtbXtk,ω,rdr2≤3ΔtkL2∫kΔtk+1ΔtEXk−Xtk,ω2dr.
Then, the Itô isometry and the global Lipschitz condition of the functionB imply that(38)I6≤3E∫kΔtk+1ΔtBXk,rdWr−∫kΔtk+1ΔtBXtk,ω,rdWr2≤3L2∫kΔtk+1ΔtEXk−Xtk,ω2dr.
It follows from the Gronwall inequality that there exists a numberM2 such that(39)EXk+1−Xtk+1,ω2≤M2,where(40)M2=3EXk−Xtk,ω2⋅3Δtk+1L2expΔtk.
By the fact thatM2 tends to zero as Δtk⟶0, that is,(41)limΔtk⟶0M2=0,we obtain(42)limΔtk⟶0ek+1=limN⟶+∞Xk+1−Xtk+1,ω=0.
Therefore, it is mean-square convergent, and the convergence order is 0.5. We have established the theorem.


## 4. Numerical Experiments

### 4.1. Experimental Preparation

Based on the sampling statistical data from Fujian Province, PR China, shown in Table1, we consider the following stochastic differential equations of diabetes mellitus:(43)dCt=−0.03tCt+0.02tNtdt+K3sintdWt,dNt=0.05t−0.007tCt−0.002tNtdt+K3costdWt.That is,(44)λt=0.02t,θt=0.01t,It=0.05t,νt+δt=0.007t,μt=0.002t,(45)g1t=K3sint,g2t=K3cost,which are obtained by the fitting method in a one-dimensional space of real numbers. And, due to the randomness of sampling data and the periodic property of the considered noise, the sine and cosine functions are used to control the intensity of the added noise.Table 1
Sampling statistical data from Fujian Province, PR China, in five years (2012–2016).
Years
λ
t
θ
t
I
t
ν
t
+
δ
t
μ
t
2012
0.0211
0.0089
0.0573
0.00526
0.00191
2013
0.0401
0.0101
0.1010
0.00611
0.00081
2014
0.0612
0.0293
0.1536
0.00650
0.00118
2015
0.0812
0.0387
0.2006
0.00706
0.00157
2016
0.0991
0.0591
0.2436
0.00664
0.00198It follows from Theorem 3.2 that there exist solutions of SDE (43). As shown in [10], in order to obtain the Brownian trajectory, we can construct the positive time path and reflect it against point zero. We select the meshes as follows [14, 19]:(46)t=500,Δt=0.01,N+1=501.Brownian trajectories are generated by the following method:(47)W0=0,Wi+1Δt=WiΔt+ψi+1,where(48)ψi=N0,Δt,i=1,2,…,N+1.Utilizing Theorem 3.2 and the one-step numerical scheme (EM scheme [19]) to solve SDE (43) with the initial value X0, we obtain(49)Ck+1t=Ckt+−0.03tkCkt+0.02tkNktΔtk+K3sintkΔWk,Nk+1t=Nkt+0.05tk−0.007tkCkt−0.002tkNktΔtk+K3costkΔWk.Then, we obtain a numerically computed solution of SDE (43).And, we obtain the graphs for numerical approximations to the solutions in the time interval0,500 as shown in Figure 1.Figure 1
Numerical solutionsCt and Nt with corresponding different starting points X0=0.65 and X0=0.1 and K3=5.0.
(a)
(b)As we see, there exist random phenomena with different starting pointsX0=0.65 and X0=0.1 at time t=0. And, it shows the fact that the numerical results can match the reality very well, where the reality means the real statistic data of diabetes mellitus in Fujian Province, PR China.


### 4.2. Numerical Results

To begin with, in order to check the convergence of numerical approximations, we plot the curves from different starting points at the timet=0 in the same graph. In Figures 2 and 3, their starting points X0,Y0 are 0.70,0.10 and 0.65,0.11, respectively. As time goes on, the trajectories tend to be more and more close and the difference of two numerical solutions becomes zero, too. This presents the fact that whatever starting points we select, two numerical solutions arrive at the same trajectory when time goes forward. That is to say, the solution of SDE (49) is a stochastic process which is different for every ω∈Ω. These confirm the fact that the numerical results are close to the epidemic disease tendency of diabetes mellitus in Fujian Province, PR China.Figure 2
Convergence of the numerical solutionCt with different starting points.
(a)
(b)Figure 3
Convergence of the numerical solutionNt with different starting points.
(a)
(b)Secondly, to check the stability of the true solution, we plot the curves from different starting points at timet=0 in the same graph in a long time, such as 40 years. As we see from Figures 4 and 5, the starting points X0,Y0 are 0.70,0.10, 0.65,0.11, and 0.79,0.09, respectively. As time progresses, the trajectories become close in a given small neighbourhood of one orbit. This also reflects the fact that if the starting points we choose are in a given region, as we move forward in time, numerical solutions arrive at the tolerant neighbourhood of one orbit which depend on different ω∈Ω. In other words, the solution of SDE (49) is not sensitive to the change of the initial value.Figure 4
Stability of the numerical solutionCt in 40 years.Figure 5
Stability of the numerical solutionNt in 40 years.


### 4.3. Control Policy and Its Practical Use

As we can see, the graph ofNt shows that the total number of diabetes mellitus oscillates at the beginning and increases later. However, the graph of Ct shows that the number of diabetes mellitus with complication decreases at first and then increases, too. Moreover, Figures 4 and 5 also present the fact that Nt and Ct both can exist random periodic phenomenon with the period of approximately 10 years. Therefore, efforts should be taken to move the situation out of the current endemic case. That is, we need to employ proper measure, such as nutrition supplement, gene therapy, physical activity and health education, and the number of incidence of diabetes mellitus will not exacerbate at least.


## 4.1. Experimental Preparation

Based on the sampling statistical data from Fujian Province, PR China, shown in Table1, we consider the following stochastic differential equations of diabetes mellitus:(43)dCt=−0.03tCt+0.02tNtdt+K3sintdWt,dNt=0.05t−0.007tCt−0.002tNtdt+K3costdWt.That is,(44)λt=0.02t,θt=0.01t,It=0.05t,νt+δt=0.007t,μt=0.002t,(45)g1t=K3sint,g2t=K3cost,which are obtained by the fitting method in a one-dimensional space of real numbers. And, due to the randomness of sampling data and the periodic property of the considered noise, the sine and cosine functions are used to control the intensity of the added noise.Table 1
Sampling statistical data from Fujian Province, PR China, in five years (2012–2016).
Years
λ
t
θ
t
I
t
ν
t
+
δ
t
μ
t
2012
0.0211
0.0089
0.0573
0.00526
0.00191
2013
0.0401
0.0101
0.1010
0.00611
0.00081
2014
0.0612
0.0293
0.1536
0.00650
0.00118
2015
0.0812
0.0387
0.2006
0.00706
0.00157
2016
0.0991
0.0591
0.2436
0.00664
0.00198It follows from Theorem 3.2 that there exist solutions of SDE (43). As shown in [10], in order to obtain the Brownian trajectory, we can construct the positive time path and reflect it against point zero. We select the meshes as follows [14, 19]:(46)t=500,Δt=0.01,N+1=501.Brownian trajectories are generated by the following method:(47)W0=0,Wi+1Δt=WiΔt+ψi+1,where(48)ψi=N0,Δt,i=1,2,…,N+1.Utilizing Theorem 3.2 and the one-step numerical scheme (EM scheme [19]) to solve SDE (43) with the initial value X0, we obtain(49)Ck+1t=Ckt+−0.03tkCkt+0.02tkNktΔtk+K3sintkΔWk,Nk+1t=Nkt+0.05tk−0.007tkCkt−0.002tkNktΔtk+K3costkΔWk.Then, we obtain a numerically computed solution of SDE (43).And, we obtain the graphs for numerical approximations to the solutions in the time interval0,500 as shown in Figure 1.Figure 1
Numerical solutionsCt and Nt with corresponding different starting points X0=0.65 and X0=0.1 and K3=5.0.
(a)
(b)As we see, there exist random phenomena with different starting pointsX0=0.65 and X0=0.1 at time t=0. And, it shows the fact that the numerical results can match the reality very well, where the reality means the real statistic data of diabetes mellitus in Fujian Province, PR China.


## 4.2. Numerical Results

To begin with, in order to check the convergence of numerical approximations, we plot the curves from different starting points at the timet=0 in the same graph. In Figures 2 and 3, their starting points X0,Y0 are 0.70,0.10 and 0.65,0.11, respectively. As time goes on, the trajectories tend to be more and more close and the difference of two numerical solutions becomes zero, too. This presents the fact that whatever starting points we select, two numerical solutions arrive at the same trajectory when time goes forward. That is to say, the solution of SDE (49) is a stochastic process which is different for every ω∈Ω. These confirm the fact that the numerical results are close to the epidemic disease tendency of diabetes mellitus in Fujian Province, PR China.Figure 2
Convergence of the numerical solutionCt with different starting points.
(a)
(b)Figure 3
Convergence of the numerical solutionNt with different starting points.
(a)
(b)Secondly, to check the stability of the true solution, we plot the curves from different starting points at timet=0 in the same graph in a long time, such as 40 years. As we see from Figures 4 and 5, the starting points X0,Y0 are 0.70,0.10, 0.65,0.11, and 0.79,0.09, respectively. As time progresses, the trajectories become close in a given small neighbourhood of one orbit. This also reflects the fact that if the starting points we choose are in a given region, as we move forward in time, numerical solutions arrive at the tolerant neighbourhood of one orbit which depend on different ω∈Ω. In other words, the solution of SDE (49) is not sensitive to the change of the initial value.Figure 4
Stability of the numerical solutionCt in 40 years.Figure 5
Stability of the numerical solutionNt in 40 years.


## 4.3. Control Policy and Its Practical Use

As we can see, the graph ofNt shows that the total number of diabetes mellitus oscillates at the beginning and increases later. However, the graph of Ct shows that the number of diabetes mellitus with complication decreases at first and then increases, too. Moreover, Figures 4 and 5 also present the fact that Nt and Ct both can exist random periodic phenomenon with the period of approximately 10 years. Therefore, efforts should be taken to move the situation out of the current endemic case. That is, we need to employ proper measure, such as nutrition supplement, gene therapy, physical activity and health education, and the number of incidence of diabetes mellitus will not exacerbate at least.


## 5. Conclusion

The main result of this article is the numerical simulation of stochastic diabetes mellitus model. The results show that the methods are effective and the numerical results can match the results of theoretical analysis and reality. Although some progresses are made, more simple and practical models and methods will be shown in our future work.



---

*Source: 5409180-2019-06-10.xml*