Diffusion processes

[22/100%]
Show
LaTeX

Feller's diffusion approximation to the branching process. Let D={D(t):t≥0}D = \left\{ D(t): t \geq 0\right\} be a diffusion process with instantaneous mean a(t,x)=αxa(t, x) = \alpha x and instantaneous variance b(t,x)=βxb(t, x) = \beta x where α\alpha and β\beta are positive constants. Let D(0)=dD(0) = d. Show that the moment generating function of D(t)D(t) is

M(t,θ)=exp⁡{2αdθeαtβθ(1−eαt)+2α} M(t, \theta ) = \exp \left\{ \frac{2 \alpha d \theta e^{\alpha t}}{\beta \theta \left(1-e^{\alpha t}\right)+2 \alpha }\right\}

Find the mean and variance of D(t)D(t), and show that P(D(t)=0)→e−2dα/β\mathbb {P}\left(D(t) = 0\right) \rightarrow e^{-2 d \alpha / \beta } as t→∞t \rightarrow \infty.

?

Let DD be an Ornstein-Uhlenbeck process with D(0)=0D(0) = 0, and place reflecting barriers at −c-c and dd where c,d>0c, d > 0. Find the limiting distribution of DD as t→∞t \rightarrow \infty.

?

Let DD be a Wiener process with drift mm, and suppose that D(0)=0D(0) = 0. Place absorbing barriers at the points x=−ax = -a and x=bx = b where aa and bb are positive real numbers. Show that the probability pap_{a} that the process is absorbed at −a-a is given by

pa=e2mb−1e2m(a+b)−1 p_{a} = \frac{e^{2 m b}-1}{e^{2 m(a+b)}-1}
?

Let X={X(t):t≥0}X = \left\{ X(t): t \geq 0\right\} be a simple birth-death process with parameters λn=nλ\lambda_{n} = n \lambda and μn=nμ\mu_{n} = n \mu. Suggest a diffusion approximation to XX.

?

Let DD be a diffusion with instantaneous mean and variance a(t,x)a(t, x) and b(t,x)b(t, x), and let M(t,θ)=E[eθD(t)]M(t, \theta ) = \mathbb {E}\left[e^{\theta D(t)}\right], the moment generating function of D(t)D(t). Use the forward diffusion equation to derive Bartlett's equation:

∂M∂t=θa(t,∂∂θ)M+12θ2b(t,∂∂θ)M \frac{\partial M}{\partial t} = \theta a\left(t, \frac{\partial }{\partial \theta }\right) M+\frac{1}{2} \theta ^{2} b\left(t, \frac{\partial }{\partial \theta }\right) M

where we interpret

g(t,∂∂θ)M=∑nγn(t)∂nM∂θn g\left(t, \frac{\partial }{\partial \theta }\right) M = \sum _{n} \gamma _{n}(t) \frac{\partial ^{n} M}{\partial \theta ^{n}}

if g(t,x)=∑n=0∞γn(t)xng(t, x) = \sum_{n = 0}^{\infty } \gamma_{n}(t) x^{n}.

?

Write down Bartlett's equation in the case of the Wiener process DD having drift mm and instantaneous variance 11, and solve it subject to the boundary condition D(0)=0D(0) = 0.

?

Write down Bartlett's equation in the case of an Ornstein-Uhlenbeck process DD having instantaneous mean a(t,x)=−xa(t, x) = -x and variance b(t,x)=1b(t, x) = 1, and solve it subject to the boundary condition D(0)=0D(0) = 0.

?

Let W1(t),W2(t),W3(t)W_{1}(t), W_{2}(t), W_{3}(t) be independent Wiener processes. The positive R=R(t)R = R(t) such that R2=W12+W22+W32R^{2} = W_{1}^{2}+W_{2}^{2}+W_{3}^{2} is the three-dimensional Bessel process. Show that RR is a Markov process. Is this result true in a general number nn of dimensions?

?

If WW is a Wiener process and the function g:R→Rg: \mathbb {R} \rightarrow \mathbb {R} is continuous and strictly monotone, show that g(W)g(W) is a continuous Markov process.

?

Exponential martingale, geometric Brownian motion. Let WW be a standard Wiener process and define S(t)=eat+bW(t)S(t) = e^{a t+b W(t)}. Show that:

?
(a)

SS is a Markov process,

(b)

SS is a martingale (with respect to the filtration generated by WW) if and only if a+12b2=0a+\frac{1}{2} b^{2} = 0, and in this case⁡E[S(t)]=1\operatorname {case} \mathbb {E}\left[S(t)\right] = 1.

Find the transition density for the Markov process of Exercise (13.3.9a).

?

Feller's diffusion approximation to the branching process, Example (13.3.12) revisited.

Let X={X(t):t≥0}X = \left\{ X(t): t \geq 0\right\} be the diffusion model of Example (13.3.12) with X(0)=x>0X(0) = x > 0, and with instantaneous mean and variance axa x and bxb x where a<0a < 0 and b>0b > 0. Let Y=∫0TX(u)duY = \int_{0}^{T} X(u) d u be the integrated size of the accumulated population up to the extinction time T=inf⁡{t≥0:X(t)=0}T = \inf \left\{ t \geq 0: X(t) = 0\right\}. Show that the moment generating function m(x,ψ)=Ex[e−ψY]m(x, \psi ) = \mathbb {E}_{x}\left[e^{-\psi Y}\right] satisfies

12bd2mdx2+admdx−ψm=0,ψ≥0 \frac{1}{2} b \frac{d^{2} m}{d x^{2}}+a \frac{d m}{d x}-\psi m = 0, \quad \psi \geq 0

Deduce that YY has a 'first-passage' distribution of the type of Exercise (13.4.2), and specify which.

?

Let DD be a standard Wiener process with drift mm starting from D(0)=d>0D(0) = d > 0, and suppose that there is a reflecting barrier at the origin. Show that the density function fr(t,y)f^{\mathrm{r}}(t, y) of D(t)D(t) satisfies fr(t,y)→0f^{\mathrm{r}}(t, y) \rightarrow 0 as t→∞t \rightarrow \infty if m≥0m \geq 0, whereas fr(t,y)→2∣m∣e−2∣m∣yf^{\mathrm{r}}(t, y) \rightarrow 2\left|m\right| e^{-2\left|m\right| y} for y>0y > 0, as t→∞t \rightarrow \infty if m<0m < 0.

?

Let WW be a standard Wiener process started at W(0)=x∈(0,a)W(0) = x \in (0, a), with absorbing barriers at 0 and aa. Let CC be the process WW conditioned on the event that WW is absorbed at aa. Show that the instantaneous mean and variance of CC, when C(t)=c∈(0,a)C(t) = c \in (0, a), are 1/c1 / c and 11, respectively.

?

Let TT be the first passage time by CC to the point aa. Show that Ex[T]=13(a2−x2)\mathbb {E}_{x}\left[T\right] = \frac{1}{3}\left(a^{2}-x^{2}\right) for x∈(0,a)x \in (0, a).

?

Show that Lemma 8.10 remains true for any dd - dimensional Feller process XtX_{t} with continuous paths and generator L=∑j,k=1dajk(x)∂j∂k+∑j=1dbj(x)∂jL = \sum_{j, k = 1}^{d} a_{j k}(x) \partial_{j} \partial_{k}+\sum_{j = 1}^{d} b_{j}(x) \partial_{j} such that a11(x)⩾a0>0a_{11}(x) \geqslant a_{0} > 0 and ∣b1(x)∣⩽b0<∞\left|b_{1}(x)\right| \leqslant b_{0} < \infty.

?

Let (A,D(A))(A, \mathfrak {D}(A)) be the generator of a diffusion process in the sense of Definition 23.1 and denote by a,ba, b the diffusion and drift coefficients. Show that a∈C(Rd,Rd×d)a \in \mathcal{C}\left(\mathbb {R}^{d}, \mathbb {R}^{d \times d}\right) and b∈C(Rd,Rd)b \in \mathcal{C}\left(\mathbb {R}^{d}, \mathbb {R}^{d}\right).

?

Show that under the assumptions of Proposition 23.5 we can interchange integration and differentiation: ∂2∂xj∂xk∫p(t,x,y)u(y)dy=∫∂2∂xj∂xkp(t,x,y)u(y)dy\frac{\partial^{2}}{\partial x_{j} \partial x_{k}} \int p(t, x, y) u(y) d y = \int \frac{\partial^{2}}{\partial x_{j} \partial x_{k}} p(t, x, y) u(y) d y and that the resulting function is in C∞(Rd)\mathcal{C}_{\infty }\left(\mathbb {R}^{d}\right).

?

Let L=L(x,Dx)=∑ijaij(x)∂i∂j+∑jbj(x)∂j+c(x)L = L\left(x, D_{x}\right) = \sum_{i j} a_{i j}(x) \partial_{i} \partial_{j}+\sum_{j} b_{j}(x) \partial_{j}+c(x) be a second order differential operator. Determine its formal adjoint L∗(y,Dy)L^{*}\left(y, D_{y}\right), i.e. the linear operator satisfying ⟨Lu,ϕ⟩L2=⟨u,L∗ϕ⟩L2\left\langle L u, \phi \right\rangle_{L^{2}} = \left\langle u, L^{*} \phi \right\rangle_{L^{2}} for all u,ϕ∈Cc∞(Rd)u, \phi \in \mathcal{C}_{c}^{\infty }\left(\mathbb {R}^{d}\right). What is the formal adjoint of the operator L(x,D)+∂∂tL(x, D)+\frac{\partial }{\partial t} ?

?

Complete the proof of Proposition 23.6 (Kolmogorov's forward equation).

?

Let (Bt)t⩾0\left(B_{t}\right)_{t \geqslant 0} be a BM1\mathrm{BM}^{1}. Show that Xt:=(Bt,∫0tBsds)X_{t}: = \left(B_{t}, \int_{0}^{t} B_{s} d s\right) is a diffusion process and determine its generator.

?

Let (Xt)t⩾0\left(X_{t}\right)_{t \geqslant 0} be a diffusion process with the infinitesimal generator L=L(x,D)=A∣Cc∞L = L(x, D) = \left.A\right|_{\mathcal{C}_{c}^{\infty }} as in (23.1). Write Mtu=u(Xt)−u(X0)−∫0tLu(Xr)drM_{t}^{u} = u\left(X_{t}\right) - u\left(X_{0}\right) - \int_{0}^{t} L u\left(X_{r}\right) d r.

?
(a)

Show that Mu=(Mtu)t⩾0M^{u} = \left(M_{t}^{u}\right)_{t \geqslant 0} is an L2L^{2} martingale for all u∈Cc∞(Rd)u \in \mathcal{C}_{c}^{\infty }\left(\mathbb {R}^{d}\right).

(b)

Show that Mf=(Mtf)t⩾0M^{f} = \left(M_{t}^{f}\right)_{t \geqslant 0} is a local martingale for all f∈C2(Rd)f \in \mathcal{C}^{2}\left(\mathbb {R}^{d}\right).

(c)

Denote by ⟨Mu,Mϕ⟩t,u,ϕ∈Cc∞(Rd)\left\langle M^{u}, M^{\phi }\right\rangle_{t}, u, \phi \in \mathcal{C}_{c}^{\infty }\left(\mathbb {R}^{d}\right), the quadratic covariation of the L2L^{2} martingales MuM^{u} and MϕM^{\phi }, cf. (15.8) and Theorem 15.15.b). Show that

⟨Mu,Mϕ⟩t=∫0t(L(uϕ)−uLϕ−ϕLu)(Xs)ds=∫0t∇u(Xs)⋅a(Xs)∇ϕ(Xs)ds. \begin{aligned} \left\langle M^{u}, M^{\phi }\right\rangle _{t} & = \int _{0}^{t}(L(u \phi ) - u L \phi - \phi L u)\left(X_{s}\right) d s \\ & = \int _{0}^{t} \nabla u\left(X_{s}\right) \cdot a\left(X_{s}\right) \nabla \phi \left(X_{s}\right) d s . \end{aligned}

(x⋅y\left(x \cdot y\right. denotes the Euclidean scalar product and ∇=(∂1,…,∂d)⊤.)\left.\nabla = \left(\partial_{1}, \ldots , \partial_{d}\right)^{\top }.\right)

Remark. The operator Γ(u,ϕ):=L(uϕ)−uLϕ−ϕLu\Gamma (u, \phi ): = L(u \phi ) - u L \phi - \phi L u is called carré - du - champ or mean - field operator.

References

[GrSt20] Grimmett G.R., Stirzaker D.R. Probability and Random Processes. 4th ed., Oxford University Press, 2020.
[Sch21] Schilling R.L. Brownian Motion: A Guide to Random Processes and Stochastic Calculus. 3rd ed., De Gruyter, 2021.