Problem 18.13 Medium Free 1 variation

Mean of the Vasicek Rate

Let WtW_t be a standard Brownian motion and let RtR_t satisfy the Vasicek SDE

dRt=(a−bRt) dt+σ dWt,R0=r,dR_t = (a - bR_t)\,dt + \sigma\,dW_t, \qquad R_0 = r,

for constants a,b,σ>0a,b,\sigma > 0 and r>0r>0. Write E[RT]\mathbb{E}[R_T] as an explicit function of a,b,r,Ta,b,r,T (note that it does not depend on σ\sigma). Evaluate this function for T=10T=10, a=0.2a=0.2, b=0.1b=0.1, r=1r=1: the result has the form A+BeCA + Be^{C} for integers A,B,CA,B,C. Find A+B+CA+B+C.

Work it out, then compare with the answer.

SolutionShow the worked solution

Key idea. The drift is linear in RtR_t, so the integrating factor ebte^{bt} removes RtR_t from the drift. The noise term then has deterministic integrand and mean zero.

By Itô’s product rule (the factor ebte^{bt} has finite variation, so there is no cross-variation term),

d(ebtRt)=bebtRt dt+ebt dRt=aebt dt+σebt dWt.d\big(e^{bt}R_t\big) = be^{bt}R_t\,dt + e^{bt}\,dR_t = ae^{bt}\,dt + \sigma e^{bt}\,dW_t .

Integrating from 00 to TT gives

RT=re−bT+ab(1−e−bT)+σ∫0Te−b(T−s) dWs.R_T = re^{-bT} + \frac{a}{b}\big(1-e^{-bT}\big) + \sigma\int_0^T e^{-b(T-s)}\,dW_s .

The last integral has a deterministic, square-integrable integrand. It is therefore a centred Gaussian random variable, with mean 00. Hence

E[RT]=re−bT+ab(1−e−bT),\mathbb{E}[R_T] = re^{-bT} + \frac{a}{b}\big(1-e^{-bT}\big),

which does not involve σ\sigma. Equivalently, m(t)=E[Rt]m(t)=\mathbb{E}[R_t] solves m′=a−bmm'=a-bm with m(0)=rm(0)=r.

With T=10T=10, a=0.2a=0.2, b=0.1b=0.1 and r=1r=1, we have bT=1bT=1 and a/b=2a/b=2. So

E[R10]=e−1+2(1−e−1)=2−e−1.\mathbb{E}[R_{10}] = e^{-1} + 2\big(1-e^{-1}\big) = 2 - e^{-1}.

This gives A=2A=2, B=−1B=-1, C=−1C=-1. The representation is unique because ee is transcendental. Therefore

E[RT]=re−bT+ab(1−e−bT),A+B+C=0\boxed{\mathbb{E}[R_T] = re^{-bT} + \tfrac{a}{b}\big(1-e^{-bT}\big),\qquad A+B+C = 0}

Variations

The same idea with a twist. Each solution gives only the changed step.

Variation 18.13a

Let WtW_t be a standard Brownian motion and let XtX_t satisfy the CIR-type SDE

dXt=κ(θ−Xt) dt−σXt dWt,X0=x>0,dX_t = \kappa(\theta - X_t)\,dt - \sigma\sqrt{X_t}\,dW_t, \qquad X_0 = x > 0,

for constants κ,θ,σ>0\kappa,\theta,\sigma>0. Assume 2κθ>σ22\kappa\theta > \sigma^2, which guarantees Xt>0X_t>0 almost surely, so Xt\sqrt{X_t} is well defined. Write E[XT]\mathbb{E}[X_T] as an explicit function of x,κ,θ,Tx,\kappa,\theta,T (it does not depend on σ\sigma). Evaluate this function for T=10T=10, κ=0.2\kappa=0.2, θ=2\theta=2, x=5x=5.

Work it out, then compare with the answer.

Solution

The same integrating factor gives

d(eκtXt)=κθeκt dt−σeκtXt dWt.d\big(e^{\kappa t}X_t\big) = \kappa\theta e^{\kappa t}\,dt - \sigma e^{\kappa t}\sqrt{X_t}\,dW_t .

Now the stochastic integral has a random integrand, so we must check that it has mean zero.

Step 1: E[Xt]\mathbb{E}[X_t] is finite. Localise with stopping times τn↑∞\tau_n\uparrow\infty at which the stochastic integral is stopped. This gives

E[eκ(t∧τn)Xt∧τn]≤x+θ(eκt−1).\mathbb{E}\big[e^{\kappa(t\wedge\tau_n)}X_{t\wedge\tau_n}\big] \le x+\theta(e^{\kappa t}-1).

Since X≥0X\ge 0, Fatou’s lemma yields sup⁡t≤TE[Xt]<∞\sup_{t\le T}\mathbb{E}[X_t]<\infty.

Step 2: the integral is a true martingale. By Step 1, E∫0Te2κsXs ds<∞\mathbb{E}\int_0^T e^{2\kappa s}X_s\,ds<\infty, so the integral is a true martingale with mean zero.

Taking expectations, m(t)=E[Xt]m(t)=\mathbb{E}[X_t] solves m′=κ(θ−m)m'=\kappa(\theta-m) with m(0)=xm(0)=x. Hence

E[XT]=θ+(x−θ)e−κT.\mathbb{E}[X_T]=\theta+(x-\theta)e^{-\kappa T}.

With κT=2\kappa T=2, θ=2\theta=2 and x=5x=5, this is 2+3e−2≈2.4062+3e^{-2}\approx 2.406.

E[XT]=θ+(x−θ)e−κT,E[X10]=2+3e−2\boxed{\mathbb{E}[X_T]=\theta+(x-\theta)e^{-\kappa T},\qquad \mathbb{E}[X_{10}]=2+3e^{-2}}