logo

イート・テイラー展開導出 📂確率微分方程式

イート・テイラー展開導出

定理1

dX(t)=f(Xt)dt+g(Xt)dWt,t[0,T] d X(t) = f \left( X_{t} \right) dt + g \left( X_{t} \right) d W_{t} \qquad , t \in [0, T] イートー過程が上のような自律確率微分方程式の解として与えられているとしよう。f,g:RRf,g : \mathbb{R} \to \mathbb{R}線形成長条件を満たすなら、つまりある定数KKに対して{f(Xt)K(1+Xt2)g(Xt)K(1+Xt2)\begin{cases} \left| f \left( X_{t} \right) \right| \le K \left( 1 + \left| X_{t} \right|^{2} \right) \\ \left| g \left( X_{t} \right) \right| \le K \left( 1 + \left| X_{t} \right|^{2} \right) \end{cases}が成り立ち、十分に多く微分可能なら以下が成り立つ。 Xt=X0+f(X0)0tds+g(X0)0tdWs+R X_{t} = X_{0} + f \left( X_{0} \right) \int_{0}^{t} ds + g \left( X_{0} \right) \int_{0}^{t} d W_{s} + R ここで、余剰remainderRRは以下の通りだ。LkL^{k}は導出過程で登場するオペレーターだ。 R=0tL0f(Xz)dzds+0tL1f(Xz)dWzds+0tL0g(Xz)dzdWs+0tL1g(Xz)dWzdWs R = \int_{0}^{t} L^{0} f \left( X_{z} \right) dz ds + \int_{0}^{t} L^{1} f \left( X_{z} \right) dW_{z} ds + \int_{0}^{t} L^{0} g \left( X_{z} \right) dz dW_{s} + \int_{0}^{t} L^{1} g \left( X_{z} \right) dW_{z} dW_{s}

説明

イトー・テイラー展開は、確率的テイラー公式stochastic Taylor formulaとも呼ばれる定理だ。式的には、積分の中にあったfXtf \circ X_{t}gXtg \circ X_{t}t=0t=0で評価されevaluated定数項として外に出てきて、それにより発生する誤差をRRでまとめたものと考えることができる。

導出

X(t)=X0+0tf(s)ds+0tg(s)dWs X (t) = X_{0} + \int_{0}^{t} f(s) ds + \int_{0}^{t} g(s) d W_{s} イートー過程の積分形をこのように考えてみよう。

イートーの公式: イートー過程{Xt}t0\left\{ X_{t} \right\}_{t \ge 0}が与えられているとしよう。 dXt=udt+vdWt d X_{t} = u dt + v d W_{t} 関数V(t,Xt)=VC2([0,)×R)V \left( t, X_{t} \right) = V \in C^{2} \left( [0,\infty) \times \mathbb{R} \right)に対してYt:=V(t,Xt)Y_{t} := V \left( t, X_{t} \right)とすると、{Yt}\left\{ Y_{t} \right\}もイートー過程で、次が成り立つ。 dYt=Vtdt+VxdXt+12Vxx(dXt)2=(Vt+Vxu+12Vxxv2)dt+VxvdWt \begin{align*} d Y_{t} =& V_{t} dt + V_{x} d X_{t} + {{ 1 } \over { 2 }} V_{xx} \left( d X_{t} \right)^{2} \\ =& \left( V_{t} + V_{x} u + {{ 1 } \over { 2 }} V_{xx} v^{2} \right) dt + V_{x} v d W_{t} \end{align*}

XtX_{t}に二次微分可能で連続な関数hC2(R)h \in C^{2} \left( \mathbb{R} \right)を取ると、イートーの公式により h(Xt)=h(X0)+0t[f(Xs)Xh(Xs)+12[g(Xs)]22X2h(Xs)]ds+0tg(Xs)Xh(Xs)dWs=h(X0)+0tL1h(Xs)ds+0tL1h(Xs)dWS \begin{align*} h \left( X_{t} \right) =& h \left( X_{0} \right) + \int_{0}^{t} \left[ f \left( X_{s} \right) {{ \partial } \over { \partial X }} h \left( X_{s} \right) + {{ 1 } \over { 2 }} \left[ g \left( X_{s} \right) \right]^{2} {{ \partial^{2} } \over { \partial X^{2} }} h \left( X_{s} \right) \right] ds \\ & + \int_{0}^{t} g \left( X_{s} \right) {{ \partial } \over { \partial X }} h \left( X_{s} \right) d W_{s} \\ =& h \left( X_{0} \right) + \int_{0}^{t} L^{1} h \left( X_{s} \right) ds + \int_{0}^{t} L^{1} h \left( X_{s} \right) d W_{S} \end{align*} ここで、L0L^{0}L1L^{1}は以下のように定義されるオペレーターだ。 L0:=fX+12g22X2L1:=gX \begin{align*} L^{0} &:= f {{ \partial } \over { \partial X }} + {{ 1 } \over { 2 }} g^{2} {{ \partial^{2} } \over { \partial X ^{2} }} \\ L^{1} &:= g {{ \partial } \over { \partial X }} \end{align*} これを形式的にh=fh = fh=gh = gに適用して、元の与えられていた積分形のイートー過程のf(Xt)f \left( X_{t} \right)g(Xt)g \left( X_{t} \right)に代入すると、次を得る。 X(t)=X0+0tf(s)ds+0tg(s)dWs=X0+0t(f(X0)+0sL1f(Xz)dz+0sL1f(Xz)dWz)ds+0t(g(X0)+0sL1g(Xz)ds+0sL1g(Xz)dWz)dWs=X0+0tf(X0)ds+0tg(X0)dWs+R=X0+f(X0)0tds+g(X0)0tdWs+R \begin{align*} X (t) =& X_{0} + \int_{0}^{t} f(s) ds + \int_{0}^{t} g(s) d W_{s} \\ =& X_{0} + \int_{0}^{t} \left( {\color{Red} f \left( X_{0} \right)} + \int_{0}^{s} L^{1} f \left( X_{z} \right) dz + \int_{0}^{s} L^{1} f \left( X_{z} \right) d W_{z} \right) ds \\ & + \int_{0}^{t} \left( {\color{Red} g \left( X_{0} \right)} + \int_{0}^{s} L^{1} g \left( X_{z} \right) ds + \int_{0}^{s} L^{1} g \left( X_{z} \right) d W_{z} \right) d W_{s} \\ =& X_{0} + \int_{0}^{t} f \left( X_{0} \right) ds + \int_{0}^{t} g \left( X_{0} \right) d W_{s} + R \\ =& X_{0} + f \left( X_{0} \right) \int_{0}^{t} ds + g \left( X_{0} \right) \int_{0}^{t} d W_{s} + R \end{align*}


  1. Panik. (2017). Stochastic Differential Equations: An Introduction with Applications in Population Dynamics Modeling: p215. ↩︎