6.8. Fórmula de Itô

Uma outra propriedade fundamental da integral de Itô diz respeito a uma fórmula relacionada a mudanças de variáveis.

Considere um processo real {Xt}t0\{X_t\}_{t \geq 0} e uma função real contínua u:RR.u:\mathbb{R} \rightarrow \mathbb{R}. Podemos definir um outro processo {Yt}t0\{Y_t\}_{t \geq 0} através de

Yt=u(Xt). Y_t = u(X_t).

A questão é sobre a relação entre dYt\mathrm{d}Y_t e dXt.\mathrm{d}X_t. Quando uu e {Xt}t0\{X_t\}_{t \geq 0} são diferenciáveis, temos simplesmente

dYtdt=u(Xt)dXtdt. \frac{\mathrm{d}Y_t}{\mathrm{d}t} = u'(X_t)\frac{\mathrm{d}X_t}{\mathrm{d}t}.

Mesmo quando os caminhos amostrais de {Xt}t0\{X_t\}_{t \geq 0} são apenas de variação limitada, ainda temos {Yt}t0\{Y_t\}_{t \geq 0} de variação limitada e ainda temos que

dYt=u(Xt)dXt, \mathrm{d}Y_t = u'(X_t)\mathrm{d}X_t,

no sentido de integral de Riemann-Stieltjes, i.e.

abg(t) dYt=abu(Xt) dXt. \int_a^b g(t) \;\mathrm{d}Y_t = \int_a^b u'(X_t)\;\mathrm{d}X_t.

Mas e no caso em que XtX_t não é nem de variação limitada? Por exemplo, se quiseremos calcular Xt=sin(Wt),X_t = \sin(W_t), para o processo de Wiener, ou qualquer outra Xt=u(Wt)X_t = u(W_t)? Podemos começar por aí e calcular dXtdX_t para funções de processos de Wiener. Mas e se quiseremos, em cima disso, consider Yt=v(Xt)Y_t = v(X_t) e ainda Zt=w(Yt),Z_t = w(Y_t), etc.

Por conta disso, consideramos uma classe de processos que é, em um certo sentido, fechada para composições. E então calculamos a diferencial para processos nessa classe. Essa classe é a de processos de Itô, ou de difusão, ou seja, processos {Xt}t\{X_t\}_t satisfazendo uma equação da forma

dXt=At dt+Bt dWt, \mathrm{d}X_t = A_t \;\mathrm{d}t + B_t \;\mathrm{d}W_t,

para outros dados processos {At}t\{A_t\}_t e {Bt}t\{B_t\}_t adaptados a {Wt}t.\{W_t\}_t.

Para esses processos, teremos uma correção na fórmula de mudança de variáveis, nos levando à fórmula de Itô, que, nesse caso autônomo, tem a forma:

dYt=u(Xt) dXt+12Bt2u(Xt) dt. \mathrm{d}Y_t = u'(X_t) \;\mathrm{d}X_t + \frac{1}{2} B_t^2 u''(X_t)\;\mathrm{d}t.

Essa mesma análise pode ser feita no caso em que uu depende também da variável temporal, i.e. u:[0,)×RR.u:[0, \infty) \times \mathbb{R} \rightarrow \mathbb{R}. Nesse caso, temos

Yt=u(t,Xt). Y_t = u(t, X_t).

Denotamos o primeiro caso por autônomo e este por não autônomo. No caso não autônomo, a fórmula de Itô se torna

dYt=ut(t,Xt) dt+ux(t,Xt) dXt+12uxxBt2 dt. \mathrm{d}Y_t = u_t(t, X_t)\;\mathrm{d}t + u_x(t, X_t)\;\mathrm{d}X_t + \frac{1}{2} u_{xx} B_t^2\;\mathrm{d}t.

Observe que, em ambos os casos, se Bt=0,B_t = 0, para todo t,t, então {Xt}t0\{X_t\}_{t \geq 0} se reduz a um processo com caminhos diferenciáveis (ou de variação limitada) e recuperamos a fórmula clássica.

Ideia do termo de correção

O motivo da correção é simples. A ideia fundamental é que dWt\mathrm{d}W_t é de ordem dt1/2.\mathrm{d}t^{1/2}. Para obtermos uma expressão para dYt\mathrm{d}Y_t que seja até ordem dt,\mathrm{d}t, devemos manter o termo seguinte dWt2dt\mathrm{d}W_t^2 \sim \mathrm{d}t na expansão de Taylor, que no caso

Yt=u(Wt) Y_t = u(W_t)

significa

dYt=u(Wt) dWt+12u(Wt) dWt2u(Wt) dWt+12u(Wt) dt. \mathrm{d}Y_t = u'(W_t)\;\mathrm{d}W_t + \frac{1}{2}u''(W_t)\;\mathrm{d}W_t^2 \sim u'(W_t)\;\mathrm{d}W_t + \frac{1}{2}u''(W_t)\;\mathrm{d}t.

No caso Yt=u(Xt)Y_t = u(X_t) para um processo de difusão, temos, formalmente,

dYt=u(Xt) dXt+12u(Xt) dXt2=u(Xt)(At dt+Bt dWt)+12u(Xt)(At dt+Bt dWt)2=u(Xt)At dt+u(Xt)Bt dWt+12u(Xt)(At2 dt2+2AtBt dtdWt+Bt2 dWt2)=u(Xt)At dt+u(Xt)Bt dWt+12u(Xt)Bt2 dt. \begin{align*} \mathrm{d}Y_t & = u'(X_t)\;\mathrm{d}X_t + \frac{1}{2}u''(X_t)\;\mathrm{d}X_t^2 \\ & = u'(X_t)\left(A_t \;\mathrm{d}t + B_t\;\mathrm{d}W_t\right) + \frac{1}{2}u''(X_t)\left(A_t \;\mathrm{d}t + B_t\;\mathrm{d}W_t\right)^2 \\ & = u'(X_t)A_t \;\mathrm{d}t + u'(X_t) B_t\;\mathrm{d}W_t + \frac{1}{2}u''(X_t)\left(A_t^2 \;\mathrm{d}t^2 + 2A_tB_t\;\mathrm{d}t\mathrm{d}W_t + B_t^2\;\mathrm{d}W_t^2\right) \\ & = u'(X_t)A_t \;\mathrm{d}t + u'(X_t) B_t\;\mathrm{d}W_t + \frac{1}{2}u''(X_t) B_t^2\;\mathrm{d}t. \end{align*}

Acima, usamos que dt2\mathrm{d}t^2 é de ordem quadrática e que dtdWt\mathrm{d}t\mathrm{d}W_t é de ordem dt3/2,\mathrm{d}t^{3/2}, portanto são de ordems superiores e desprezados na aproximação de ordem linear, enquanto que dWt2dt.\mathrm{d}W_t^2 \sim \mathrm{d}t.

A mesma ideia funciona no caso não autônomo Yt=u(t,Xt).Y_t = u(t, X_t).

Mas essa ideia não constitui uma demonstração e ainda está muito vaga. Vamos tratar do detalhes, a seguir.

Fórmula de Itô no caso de um processo de Wiener

Para deixar as ideias claras, vamos considerar o caso mais simples em que Xt=WtX_t=W_t é o próprio processo de Wiener e u=u(x)u=u(x) é uma função duas vezes continuamente diferenciável definida para todo xR.x\in\mathbb{R}. Estamos interessados em obter uma fórmula para dYt,\mathrm{d}Y_t, onde {Yt}t0\{Y_t\}_{t\geq 0} é dado por

Yt=u(Wt). Y_t = u(W_t).

Considere uma malha temporal 0=t0<t1<<tn=T.0 = t_0 < t_1 < \cdots < t_n = T. Escrevemos

Yt=Y0+j=1n(YtjYtj1), Y_t = Y_0 + \sum_{j = 1}^n (Y_{t_j} - Y_{t_{j-1}}),

que é uma soma telescópica. Como uu é duas vezes continuamente diferenciável, temos a fórmula de expansão de Taylor com resto

u(b)=u(a)+abu(s) ds=u(a)+ab(u(a)+asu(ξ) dξ)ds=u(a)+u(a)(ba)+abξbu(ξ);ds dξ=u(a)+u(a)(ba)+abu(ξ)(bξ) dξ. \begin{align*} u(b) & = u(a) + \int_a^b u'(s) \;\mathrm{d}s \\ & = u(a) + \int_a^b \left(u'(a) + \int_a^s u''(\xi)\;\mathrm{d}\xi \right)\mathrm{d}s \\ & = u(a) + u'(a) (b - a) + \int_a^b \int_\xi^b u''(\xi);\mathrm{d}s \;\mathrm{d}\xi \\ & = u(a) + u'(a) (b - a) + \int_a^b u''(\xi)(b - \xi)\;\mathrm{d}\xi. \end{align*}

Assim, podemos escrever, fazendo a=Wtj1a = W_{t_{j-1}} e b=Wtjb=W_{t_j} em cada termo,

YtjYtj1=u(Wtj)u(Wtj1)=u(Wtj1)(WtjWtj1)+Wtj1Wtju(ξ)(Wtjξ) dξ=u(Wtj1)(WtjWtj1)+12u(Wtj1)(WtjWtj1)2+Wtj1Wtj(u(ξ)u(Wtj1))(Wtjξ) dξ. \begin{align*} Y_{t_j} - Y_{t_{j-1}} = u(W_{t_j}) - u(W_{t_{j-1}}) & = u'(W_{t_{j-1}})(W_{t_j} - W_{t_{j-1}}) + \int_{W_{t_{j-1}}}^{W_{t_j}} u''(\xi)(W_{t_j} - \xi) \;\mathrm{d}\xi \\ & = u'(W_{t_{j-1}})(W_{t_j} - W_{t_{j-1}}) \\ & \qquad + \frac{1}{2}u''(W_{t_{j-1}})(W_{t_j} - W_{t_{j-1}})^2 \\ & \qquad + \int_{W_{t_{j-1}}}^{W_{t_j}} (u''(\xi) - u''(W_{t_{j-1}}))(W_{t_j} - \xi) \;\mathrm{d}\xi. \end{align*}

Portanto,

YTY0=j=1nu(Wtj1)(WtjWtj1)+12j=1nu(Wtj1)(WtjWtj1)2+j=1nWtj1Wtj(u(ξ)u(Wtj1))(Wtjξ) dξ. \begin{align*} Y_T - Y_0 & = \sum_{j=1}^n u'(W_{t_{j-1}})(W_{t_j} - W_{t_{j-1}}) \\ & \qquad + \frac{1}{2} \sum_{j=1}^n u''(W_{t_{j-1}})(W_{t_j} - W_{t_{j-1}})^2 \\ & \qquad + \sum_{j=1}^n \int_{W_{t_{j-1}}}^{W_{t_j}} (u''(\xi) - u''(W_{t_{j-1}}))(W_{t_j} - \xi) \;\mathrm{d}\xi. \end{align*}

No limite de refinamento da malha, o primeiro termo converge para a integral de Itô de {u(Wt)}t0,\{u'(W_t)\}_{t \geq 0}, i.e.

j=1nu(Wtj1)(WtjWtj1)0Tu(Wt) dWt. \sum_{j=1}^n u'(W_{t_{j-1}})(W_{t_j} - W_{t_{j-1}}) \rightarrow \int_0^T u'(W_t)\;\mathrm{d}W_t.

No segundo termo, usamos que E[(WtjWtj1)2]=(tjtj1).\mathbb{E}[(W_{t_j} - W_{t_{j-1}})^2] = (t_j - t_{j-1}). É possível mostrar que

12j=1nu(Wtj1)(WtjWtj1)2120Tu(Wt) dt, \frac{1}{2}\sum_{j=1}^n u''(W_{t_{j-1}}) (W_{t_j} - W_{t_{j-1}})^2 \rightarrow \frac{1}{2}\int_0^T u''(W_t) \;\mathrm{d}t,

em média quadrática.

O último termo converge para zero em média quadrática, graças à continuidade da segunda derivada de u=u(x),u = u(x),

j=1nWtj1Wtj(u(ξ)u(Wtj1))(Wtjξ) dξ0. \sum_{j=1}^n \int_{W_{t_{j-1}}}^{W_{t_j}} (u''(\xi) - u''(W_{t_{j-1}}))(W_{t_j} - \xi) \;\mathrm{d}\xi \rightarrow 0.

Logo,

YTY0=0Tu(Wt) dWt+120Tu(Wt) dt. Y_T - Y_0 = \int_0^T u'(W_t)\;\mathrm{d}W_t + \frac{1}{2}\int_0^T u''(W_t) \;\mathrm{d}t.

Em forma diferencial, escrevemos que

dYt=u(Wt) dWt+12u(Wt) dt, \mathrm{d}Y_t = u'(W_t)\;\mathrm{d}W_t + \frac{1}{2}u''(W_t)\;\mathrm{d}t,

que é um caso particular da fórmula de Itô mencionada acima. Veremos agora a fórmula de Itô para processos de Itô.

Vamos ver as convergências acima em mais detalhes.

Convergência do segundo termo

Nesse caso, não há nada de especial em ser a segunda derivada. Basta mostrar que

j=1nHtj1(WtjWtj1)2120THt dt, \sum_{j=1}^n H_{t_{j-1}} (W_{t_j} - W_{t_{j-1}})^2 \rightarrow \frac{1}{2}\int_0^T H_t \;\mathrm{d}t,

para um processo qualquer {Ht}t\{H_t\}_t progessivamente mensurável e de quadrado integrável (como o é u(Wt)u''(W_t) quando uu é suave e {Wt}t\{W_t\}_t é o processo de Wiener).

O primeiro passo é aproximar, no sentido de média quadrática, o processo {Ht}t\{H_t\}_t por processos escada {Htm}t\{H_t^m\}_t adaptados a {Wt}t,\{W_t\}_t, conforme feito na parte de existência de integral de Itô. Cada HtmH_t^m é constante em uma partição 0=a0<a1<<am=T.0 = a_0 < a_1 < \ldots < a_m = T. Em cada subintervalo Ak=[ak1,ak),A_k = [a_{k-1}, a_k), temos índices

Jk={tjAk} J_k = \{t_j \in A_k\}

da partição {tj}j\{t_j\}_j usada no cálculo da integral. Assim, basta mostrar que

jJk(WtjWtj1)2akak1. \sum_{j\in J_k} (W_{t_j} - W_{t_{j-1}})^2 \rightarrow a_k - a_{k-1}.

Uma vez mostrado isso,

j=1nHtj1m(WtjWtj1)2=k=1nHak1jJk(WtjWtj1)2k=1nHak1m(akak1)=0THtm dt. \sum_{j=1}^n H_{t_{j-1}}^m (W_{t_j} - W_{t_{j-1}})^2 = \sum_{k=1}^n H_{a_{k-1}} \sum_{j\in J_k} (W_{t_j} - W_{t_{j-1}})^2 \rightarrow \sum_{k=1}^n H_{a_{k-1}}^m (a_k - a_{k-1}) = \int_0^T H_t^m \;\mathrm{d}t.

Agora, passando ao limite em m,m\rightarrow \infty, obtemos

j=1nHtj1(WtjWtj1)20THt dt. \sum_{j=1}^n H_{t_{j-1}} (W_{t_j} - W_{t_{j-1}})^2 \rightarrow \int_0^T H_t \;\mathrm{d}t.

Falta mostrar, então, que

jJk(WtjWtj1)2akak1. \sum_{j\in J_k} (W_{t_j} - W_{t_{j-1}})^2 \rightarrow a_k - a_{k-1}.

em média quadrática. Isso é feito diretamente. Defina τ=akak1\tau = a_k - a_{k-1} para facilitar a notação. Então

E[(jJk(WtjWtj1)2τ)2]=E[(jJk(WtjWtj1)2)2]+2τE[(jJk(WtjWtj1))]+τ2 \mathbb{E}\left[ \left( \sum_{j\in J_k} (W_{t_j} - W_{t_{j-1}})^2 - \tau \right)^2 \right] = \mathbb{E}\left[ \left( \sum_{j\in J_k} (W_{t_j} - W_{t_{j-1}})^2\right)^2 \right] + 2\tau \mathbb{E}\left[ \left( \sum_{j\in J_k} (W_{t_j} - W_{t_{j-1}})\right) \right] + \tau^2

Temos

E[(jJk(WtjWtj1)2)2]=i,jE[(WtiWti1)2(WtjWtj1)2]=ijE[(WtiWti1)2]E[(WtjWtj1)2]+jE[(WtjWtj1)4]=ijE[(WtiWti1)2]E[(WtjWtj1)2]+3jE[(WtjWtj1)2]2=ij(titi1)(tjtj1)+3j(tjtj1)2=i,j(titi1)(tjtj1)+2j(tjtj1)2=(i(titi1))(j(tjtj1))+2j(tjtj1)2=τ2+2j(tjtj1)2. \begin{align*} \mathbb{E}\left[ \left( \sum_{j\in J_k} (W_{t_j} - W_{t_{j-1}})^2\right)^2 \right] & = \sum_{i, j} \mathbb{E}\left[(W_{t_i} - W_{t_{i-1}})^2(W_{t_j} - W_{t_{j-1}})^2\right] \\ & = \sum_{i \neq j}\mathbb{E}\left[(W_{t_i} - W_{t_{i-1}})^2\right] \mathbb{E}\left[(W_{t_j} - W_{t_{j-1}})^2\right] + \sum_j \mathbb{E}\left[(W_{t_j} - W_{t_{j-1}})^4\right] \\ & = \sum_{i \neq j}\mathbb{E}\left[(W_{t_i} - W_{t_{i-1}})^2\right] \mathbb{E}\left[(W_{t_j} - W_{t_{j-1}})^2\right] + 3\sum_j \mathbb{E}\left[(W_{t_j} - W_{t_{j-1}})^2\right]^2 \\ & = \sum_{i \neq j} (t_i - t_{i-1})(t_j - t_{j-1}) + 3\sum_j (t_j - t_{j-1})^2 \\ & = \sum_{i, j} (t_i - t_{i-1})(t_j - t_{j-1}) + 2\sum_j (t_j - t_{j-1})^2 \\ & = \left(\sum_i (t_i - t_{i-1})\right)\left(\sum_j (t_j - t_{j-1})\right) + 2\sum_j (t_j - t_{j-1})^2 \\ & = \tau^2 + 2\sum_j (t_j - t_{j-1})^2. \end{align*}

Também temos

2τE[(jJk(WtjWtj1))]=2τjJkE[(WtjWtj1)2]=2τjJk(tjtj1)=2τ2. 2\tau \mathbb{E}\left[ \left( \sum_{j\in J_k} (W_{t_j} - W_{t_{j-1}})\right) \right] = 2\tau \sum_{j\in J_k}\mathbb{E}\left[(W_{t_j} - W_{t_{j-1}})^2 \right] = 2\tau \sum_{j\in J_k} (t_j - t_{j-1}) = 2\tau^2.

Assim, os termos com τ2\tau^2 se cancelam e sobra

E[(jJk(WtjWtj1)2τ)2]=2j(tjtj1)22maxjJk(tjtj1)j(tjtj1)=2τmaxjJk(tjtj1)0. \mathbb{E}\left[ \left( \sum_{j\in J_k} (W_{t_j} - W_{t_{j-1}})^2 - \tau \right)^2 \right] = 2\sum_j (t_j - t_{j-1})^2 \leq 2\max_{j\in J_k}(t_j - t_{j-1})\sum_j(t_j - t_{j-1}) = 2\tau \max_{j\in J_k}(t_j - t_{j-1}) \rightarrow 0.

Isso completa a demonstração da convergência em média quadrática do segundo termo.

Convergência do terceiro termo

Por densidade, basta mostrar no caso em que uu é três vezes continuamente diferenciável. Vamos assumir isso.

Quase todo caminho amostral de {Wt}t\{W_t\}_t é H\"older contínou com expoente θ<1/2\theta < 1/2 e é limitado em [0,T].[0, T]. Assim, fixado um caminho amostral, existe C>0C>0 tal que

Wt+τ(ω)Wt(ω)Cτθ |W_{t+\tau}(\omega) - W_t(\omega)| \leq C\tau ^\theta

Em particular, como W0=0,W_0 = 0,

Wt(ω)CTθ, |W_{t}(\omega)| \leq CT^\theta,

para todo t[0,T].t\in [0, T]. Além disso, assumindo uu'' continuamente diferenciável, dado ϵ>0,\epsilon > 0, existe K>0K > 0 tal que

w1,w2CTθu(w2)u(w1)Kw2w1. w_1, w_2 \leq CT^\theta \Rightarrow |u(w_2) - u(w_1)| \leq K|w_2 - w_1|.

Assim,

Wtj1Wtj(u(ξ)u(Wtj1))(Wtjξ) dξKWtjWtj13C2Ktjtj13θ. \left|\int_{W_{t_{j-1}}}^{W_{t_j}} (u''(\xi) - u''(W_{t_{j-1}}))(W_{t_j} - \xi) \;\mathrm{d}\xi\right| \leq K |W_{t_j} - W_{t_{j-1}}|^3 \leq C^2K |t_j - t_{j-1}|^{3\theta}.

Escolhendo 1/3<θ<1/2,1/3 < \theta < 1/2, temos 3θ>1,3\theta > 1, de modo que

j=1nWtj1Wtj(u(ξ)u(Wtj1))(Wtjξ) dξC2KTmaxjtjtj13θ10, \sum_{j=1}^n \int_{W_{t_{j-1}}}^{W_{t_j}} (u''(\xi) - u''(W_{t_{j-1}}))(W_{t_j} - \xi) \;\mathrm{d}\xi \leq C^2 K T \max_j |t_j - t_{j-1}|^{3\theta - 1} \rightarrow 0,

quando a malha é refinada. Ou seja, obtemos a convergência quase certamente, que implica na convergência em média quadrática.

Processo de Itô

Como dito acima, um processo de Itô é um processo {Xt}t0\{X_t\}_{t \geq 0} satisfazendo uma equação estocástica da forma

dXt=At dt+Bt dWt, \mathrm{d}X_t = A_t \;\mathrm{d}t + B_t \;\mathrm{d}W_t,

onde {Wt}t0\{W_t\}_{t \geq 0} é um processo de Wiener e {At}t0\{A_t\}_{t \geq 0} e {Bt}t0\{B_t\}_{t \geq 0} são não antecipativos em relação ao processo de Wiener.

Fórmula de Itô no caso autônomo

Considere, novamente, uma malha temporal 0=t0<t1<<tn=T0 = t_0 < t_1 < \cdots < t_n = T com maxj=1,,n{tjtj1}0.\max_{j=1, \ldots, n}\{t_j - t_{j-1}\} \rightarrow 0. Escrevemos

Yt=Y0+j=1n(YtjYtj1). Y_t = Y_0 + \sum_{j = 1}^n (Y_{t_j} - Y_{t_{j-1}}).

Assumindo, agora, que uu seja duas vezes continuamente diferenciável, podemos escrever

YtjYtj1=u(Xtj)u(Xtj1)=u(Xtj1)(XtjXtj1)+Xtj1Xtju(x)(Xtjx) dx=u(Xtj1)(XtjXtj1)+12u(Xtj1)(XtjXtj1)2+Xtj1Xtj(u(x)u(Xtj1))(Xtjx) dx. \begin{align*} Y_{t_j} - Y_{t_{j-1}} = u(X_{t_j}) - u(X_{t_{j-1}}) & = u'(X_{t_{j-1}})(X_{t_j} - X_{t_{j-1}}) + \int_{X_{t_{j-1}}}^{X_{t_j}} u''(x)(X_{t_j} - x) \;\mathrm{d}x \\ & = u'(X_{t_{j-1}})(X_{t_j} - X_{t_{j-1}}) \\ & \qquad + \frac{1}{2}u''(X_{t_{j-1}})(X_{t_j} - X_{t_{j-1}})^2 \\ & \qquad + \int_{X_{t_{j-1}}}^{X_{t_j}} (u''(x) - u''(X_{t_{j-1}}))(X_{t_j} - x) \;\mathrm{d}x. \end{align*}

Assim,

YTY0=j=1nu(Xtj1)(XtjXtj1)+12j=1nu(Xtj1)(XtjXtj1)2+j=1nXtj1Xtj(u(x)u(Xtj1))(Xtjx) dx. \begin{align*} Y_T - Y_0 & = \sum_{j=1}^n u'(X_{t_{j-1}})(X_{t_j} - X_{t_{j-1}}) \\ & \qquad + \frac{1}{2} \sum_{j=1}^n u''(X_{t_{j-1}})(X_{t_j} - X_{t_{j-1}})^2 \\ & \qquad + \sum_{j=1}^n \int_{X_{t_{j-1}}}^{X_{t_j}} (u''(x) - u''(X_{t_{j-1}}))(X_{t_j} - x) \;\mathrm{d}x. \end{align*}

No limite de refinamento da malha, o primeiro termo converge para a integral de Itô de {u(Xt)}t0\{u'(X_t)\}_{t \geq 0} com respeito {Xt}t0.\{X_t\}_{t \geq 0}. O último termo converge para zero, graças à continuidade da segunda derivada de u=u(x).u = u(x). Em relação ao segundo termo, temos

(XtjXtj1)2(Atj1(tjtj1)+Btj1(WtjWtj1))2Atj12(tjtj1)2+2Atj1Btj1(tjtj1)(WtjWtj1)+Btj12(WtjWtj1)2 \begin{align*} (X_{t_j} - X_{t_{j-1}})^2 & \approx \left(A_{t_{j-1}}(t_j - t_{j-1}) + B_{t_{j-1}}(W_{t_j} - W_{t_{j-1}}) \right)^2 \\ & \approx A_{t_{j-1}}^2(t_j - t_{j-1})^2 + 2A_{t_{j-1}}B_{t_{j-1}}(t_j - t_{j-1})(W_{t_j} - W_{t_{j-1}}) + B_{t_{j-1}}^2(W_{t_j} - W_{t_{j-1}})^2 \end{align*}

Ao somarmos e tomarmos o limite, os dois primeiros termos se anulam, pois contém um Δt\Delta t extra que converge para zero. Mais precisamente, no primeiro termo, vale, quase certamente,

12j=1nu(Xtj1)Atj12(tjtj1)212maxumaxjtjtj1j=1nAtj12(tjtj1)0×0TAt2 dt=0. \left|\frac{1}{2}\sum_{j=1}^n u''(X_{t_{j-1}}) A_{t_{j-1}}^2(t_j - t_{j-1})^2\right| \leq \frac{1}{2}\max |u''| \max_j|t_j - t_{j-1}| \sum_{j=1}^n A_{t_{j-1}}^2(t_j - t_{j-1}) \rightarrow 0 \times \int_0^T A_t^2 \;\mathrm{d}t = 0.

No segundo termo, em média quadrática, temos, como na isometria de Itô,

E[12j=1nu(Xtj1)2Atj1Btj1(tjtj1)(WtjWtj1)2]12maxumaxjtjtj1E[2j=1nAtj1Btj1(WtjWtj1)2]2maxu×0×0TE[At2Bt2] dt=0. \begin{align*} \mathbb{E}\left[\left|\frac{1}{2}\sum_{j=1}^n u''(X_{t_{j-1}}) 2A_{t_{j-1}}B_{t_{j-1}}(t_j - t_{j-1})(W_{t_j} - W_{t_{j-1}})\right|^2 \right] & \leq \frac{1}{2}\max |u''| \max_j|t_j - t_{j-1}| \mathbb{E}\left[\left| 2\sum_{j=1}^n A_{t_{j-1}}B_{t_{j-1}}(W_{t_j} - W_{t_{j-1}})\right|^2 \right] \\ & \rightarrow 2\max |u''| \times 0 \times \int_0^T \mathbb{E}[A_t^2 B_t^2] \;\mathrm{d}t = 0. \end{align*}

No último termo, usando que E[(WtjWtj1)2]=(tjtj1),\mathbb{E}[(W_{t_j} - W_{t_{j-1}})^2] = (t_j - t_{j-1}), obtemos

12j=1nu(Xtj1)Btj12(WtjWtj1)2120Tu(Xt)Bt2 dt. \frac{1}{2}\sum_{j=1}^n u''(X_{t_{j-1}}) B_{t_{j-1}}^2(W_{t_j} - W_{t_{j-1}})^2 \rightarrow \frac{1}{2}\int_0^T u''(X_t) B_t^2 \;\mathrm{d}t.

Informalmente, podemos escrever

(dX)2=At2 (dt)2+2AtBt dt dWt+Bt2(dWt)2=Bt2 dt, (\mathrm{d}X)^2 = A_t^2 \;(\mathrm{d}t)^2 + 2A_tB_t \;\mathrm{d}t\;\mathrm{d}W_t + B_t^2(\mathrm{d}W_t)^2 = B_t^2\;\mathrm{d}t,

com

(dt)2=0,dt dWt=0,(dWt)2=dt. (\mathrm{d}t)^2 = 0, \quad \mathrm{d}t\;\mathrm{d}W_t = 0, \quad (\mathrm{d}W_t)^2 = \mathrm{d}t.

De qualquer forma, obtemos, no limite,

YT=Y0+0Tu(Xt) dXt+120Tu(Xt)Bt dt. Y_T = Y_0 + \int_0^T u'(X_t)\;\mathrm{d}X_t + \frac{1}{2}\int_0^T u''(X_t) B_t \;\mathrm{d}t.

Ou seja, obtemos a fórmula de Itô no caso autônomo:

dYt=u(Xt) dXt+12Bt2u(Xt) dt. \mathrm{d}Y_t = u'(X_t) \;\mathrm{d}X_t + \frac{1}{2} B_t^2 u''(X_t)\;\mathrm{d}t.

Fórmula de Itô no caso não autônomo

Quando u=u(t,x)u = u(t, x) é duas vezes continuamente diferenciável e {Xt}t0\{X_t\}_{t\geq 0} é um processo de Itô satisfazendo

dXt=At dt+Bt dWt, \mathrm{d}X_t = A_t \;\mathrm{d}t + B_t \;\mathrm{d}W_t,

então o processo {Yt}t0\{Y_t\}_{t\geq 0} definido por Yt=u(t,Xt)Y_t = u(t, X_t) satisfaz a fórmula de Itô no caso não autônomo, a saber

dYt=ut(t,Xt) dt+ux(t,Xt) dXt+12uxxBt2 dt. \mathrm{d}Y_t = u_t(t, X_t)\;\mathrm{d}t + u_x(t, X_t)\;\mathrm{d}X_t + \frac{1}{2} u_{xx}B_t^2\;\mathrm{d}t.

Para lembrar dessa fórmula, pense primeiramente em uma expansão de segunda ordem

dYt=ut dt+ux dXt+12utt (dt)2+utx dt dXt+12uxx (dXt)2. \mathrm{d}Y_t = u_t\;\mathrm{d}t + u_x\;\mathrm{d}X_t + \frac{1}{2} u_{tt}\;(\mathrm{d}t)^2 + u_{tx} \;\mathrm{d}t \;\mathrm{d}X_t + \frac{1}{2} u_{xx}\;(\mathrm{d}X_t)^2.

Em seguida, eleva-se ao quadrado a equação para {Xt}t0\{X_t\}_{t \geq 0}:

(dXt)2=(At dt+Bt dWt)2=At2 (dt)2+2AtBt dt dWt+Bt2 (dWt)2. (\mathrm{d}X_t)^2 = \left(A_t \;\mathrm{d}t + B_t \;\mathrm{d}W_t\right)^2 = A_t^2\;(\mathrm{d}t)^2 + 2A_tB_t \;\mathrm{d}t\;\mathrm{d}W_t + B_t^2\;(\mathrm{d}W_t)^2.

Substituindo (dXt)2(\mathrm{d}X_t)^2 na expansão de segunda ordem de dYt,\mathrm{d}Y_t, obtemos

dYt=ut dt+ux dXt+12utt (dt)2+utx dt(At dt+Bt dWt)+12uxx(At2 (dt)2+2AtBt dt dWt+Bt2 (dWt)2). \begin{align*} \mathrm{d}Y_t & = u_t\;\mathrm{d}t + u_x\;\mathrm{d}X_t \\ & \qquad + \frac{1}{2} u_{tt}\;(\mathrm{d}t)^2 \\ & \qquad + u_{tx} \;\mathrm{d}t \left(A_t \;\mathrm{d}t + B_t \;\mathrm{d}W_t\right) \\ & \qquad + \frac{1}{2} u_{xx}\left(A_t^2\;(\mathrm{d}t)^2 + 2A_tB_t \;\mathrm{d}t\;\mathrm{d}W_t + B_t^2\;(\mathrm{d}W_t)^2\right). \end{align*}

Usando que, sob uma integração simples,

(dt)2=0,dt dWt=0,(dWt)2=dt, \begin{align*} (\mathrm{d}t)^2 & = 0, \\ \mathrm{d}t\;\mathrm{d}W_t & = 0, \\ (\mathrm{d}W_t)^2 & = \mathrm{d}t, \end{align*}

chega-se à fórmula de Itô

dYt=ut dt+ux dXt+12uxxBt2 dt. \mathrm{d}Y_t = u_t\;\mathrm{d}t + u_x\;\mathrm{d}X_t + \frac{1}{2} u_{xx} B_t^2 \;\mathrm{d}t.

Exercícios

  1. Mostre que

d(Wt2)=2Wt dWt+ dt. \mathrm{d}(W_t^2) = 2W_t\;\mathrm{d}W_t + \;\mathrm{d}t.
  1. Mostre que

d(tWt)=Wt dt+t dWt. \mathrm{d}(tW_t) = W_t\;\mathrm{d}t + t\;\mathrm{d}W_t.
  1. Mostre a "regra do produto" para o produto entre uma função determinística f=f(t)f=f(t) e um processo de Itô {Xt}t\{X_t\}_t:

d(f(t)Xt)=f(t)Xt dt+f(t) dXt. \mathrm{d}(f(t)X_t) = f'(t)X_t\;\mathrm{d}t + f(t)\;\mathrm{d}X_t.
  1. Considere dois processos de Wiener independentes {Wt1}t0\{W_t^1\}_{t\geq 0} e {Wt2}t0\{W_t^2\}_{t\geq 0} e considere uma função suave g=g(x1,x2),g=g(x_1, x_2), de R2\mathbb{R}^2 em R.\mathbb{R}. Obtenha a fórmula de Itô para o processo Xt=g(Wt1,Wt2),X_t = g(W_t^1, W_t^2), ou seja, mostre que

dXt=At dWt1+Bt dWt1+Ct dt \mathrm{d}X_t = A_t\;\mathrm{d}W_t^1 + B_t\;\mathrm{d}W_t^1 + C_t\;\mathrm{d}t

para processos apropriados {At}t0,\{A_t\}_{t\geq 0}, {Bt}t0\{B_t\}_{t\geq 0} e {Ct}t0.\{C_t\}_{t\geq 0}.

  1. (Regra de Itô para o produto) Considere dois processos de Itô {Xt}t0\{X_t\}_{t\geq 0} e {Yt}t0\{Y_t\}_{t\geq 0} satisfazendo

dXt=At dt+Bt dWtedYt=Ct dt+Dt dWt, \mathrm{d}X_t = A_t \;\mathrm{d}t + B_t\;\mathrm{d}W_t \quad \textrm{e} \quad \mathrm{d}Y_t = C_t \;\mathrm{d}t + D_t\;\mathrm{d}W_t,

respectivamente, em relação a um mesmo processo de Wiener {Wt}t0.\{W_t\}_{t\geq 0}. Argumente que

d(XtYt)=Yt dXt+Xt dYt+BtDt dt. \mathrm{d}(X_tY_t) = Y_t\;\mathrm{d}X_t + X_t\;\mathrm{d}Y_t + B_t D_t\;\mathrm{d}t.
Last modified: June 01, 2026. Built with Franklin.jl, using the Book Template.