Suppose I have a model for the short rate r as (W(t) is standard Brownian motion)
r(t)=c+∫t0σ(s)2(t−s)ds+∫t0σ(s)dW(s)
I then want to find the dynamics of r, but how do I do that when the process itself contains integrals w.r.t Brownian motion? I get stuck when using Ito's formula and trying to calculate the integral ∂∂t∫t0σ(s)dW(s)
Answer
Let Xt=∫t0σ(s)dWs
denote a stochastic integral in the Itô sense. In that case one can write rt=f(t,Xt) where f:(t,x)→c+∫t0σ2(s)(t−s)ds+x
and use Itô's lemma to compute the differential drt=∂tf(t,Xt)dt+∂xf(t,Xt)dXt+∂xxf(t,Xt)d⟨X⟩t
where from (1) ∂tf(t,Xt)=∂t∫t0σ2(s)(t−s)ds=∫t0σ2(s)ds+1(σ2(t)(t−t))−0(σ2(0)(0−s))=∫t0σ2(s)ds
from Leibniz integral rule and ∂xf(t,Xt)=1,∂xxf(t,Xt)=0
along with, by definition of the Itô integral: dXt=σ(t)dWt,d⟨X⟩t=σ2(t)dt
such that finally: drt=(∫t0σ2(s)ds)dt+σ(t)dWt
And indeed by integrating this last equation from 0 to t one gets: rt−r0=∫t0(∫u0σ2(s)ds)du+∫t0σ(u)dWu
and noting that ∫t0∫u0σ2(s)dsdu=∫t0∫tsσ2(s)duds=∫t0σ2(s)(t−s)ds
by Fubini theorem, one gets rt=r0+∫t0σ2(s)(t−s)ds+∫t0σ(s)dWs
No comments:
Post a Comment