Fórmula de Feynman–Kac

A fórmula de Feynman–Kac, que recebe este nome em homenagem ao físico norte-americano Richard Feynman e ao matemático polonês Mark Kac, estabelece uma ligação entre equações diferenciais parciais (EDPs) parabólicas e processos estocásticos. A fórmula oferece um método para resolver algumas EDPs pela simulação de caminhos aleatórios de um processo estocástico. Reciprocamente, uma importante classe de valores esperados de processos aleatórios pode ser computada por métodos determinísticos.

Fórmula

Considere a equação diferencial parcial

$\frac{\partial u}{\partial t}(x,t)+\mu (x,t)\frac{\partial u}{\partial x}(x,t)+\frac{1}{2}{\sigma }^2(x,t)\frac{{\partial }^{2}u}{\partial x^2}(x,t)-V(x,t)u(x,t)+f(x,t)=0,$

definida para todo $x$ em R e todo $t$ em $[0,T]$, sujeita à condição terminal

$u(x,T)=\psi (x),$

em que $\mu$, $\sigma$, $\psi$, $V$ e $f$ são funções conhecidas. $T$ é um parâmetro e $u:ℝ\times [0,T]\to ℝ$ é desconhecido. Então, a fórmula de Feynman–Kac nos diz que a solução pode ser escrita como um valor esperado condicional

$u(x,t)=E^Q[\underset{t}{\overset{T}{\int }}{e}^{-\underset{t}{\overset{r}{\int }}V(X_{\tau },\tau )d\tau }f(X_r,r)dr+e^{-\underset{t}{\overset{T}{\int }}V(X_{\tau },\tau )d\tau }\psi (X_T)|X_t=x]$

sob a medida de probabilidade $Q$, tal que $X$ é um processo de Itō dirigido pela equação

$dX=\mu (X,t)dt+\sigma (X,t)d{W}^Q,$

sendo $W^Q(t)$ um processo de Wiener (também chamado de movimento browniano) sob $Q$ e $x$ a condição inicial para $X(t)$.

Prova

Considere $u(x,t)$ a solução da EDP acima. Aplicando o lema de Itō ao processo

$Y(s)=e^{-\underset{t}{\overset{s}{\int }}V(X_{\tau },\tau )d\tau }u(X_s,s)+\underset{t}{\overset{s}{\int }}{e}^{-\underset{t}{\overset{r}{\int }}V(X_{\tau },\tau )d\tau }f(X_r,r)dr$

obtém-se

$\begin{array}{cc}dY=& d\left(e^{-\underset{t}{\overset{s}{\int }}V(X_{\tau },\tau )d\tau }\right)u(X_s,s)+e^{-\underset{t}{\overset{s}{\int }}V(X_{\tau },\tau )d\tau }du(X_s,s)\\ & +d\left(e^{-\underset{t}{\overset{s}{\int }}V(X_{\tau },\tau )d\tau }\right)du(X_s,s)+d(\underset{t}{\overset{s}{\int }}{e}^{-\underset{t}{\overset{r}{\int }}V(X_{\tau },\tau )d\tau }f(X_r,r)dr)\end{array}$

Já que

$d\left(e^{-\underset{t}{\overset{s}{\int }}V(X_{\tau },\tau )d\tau }\right)=-V(X_s,s)e^{-\underset{t}{\overset{s}{\int }}V(X_{\tau },\tau )d\tau }ds,$

o terceiro termo é $O(dtdu)$ e pode ser dispensado. Tem-se também que

$d(\underset{t}{\overset{s}{\int }}{e}^{-\underset{t}{\overset{r}{\int }}V(X_{\tau },\tau )d\tau }f(X_r,r)dr)=e^{-\underset{t}{\overset{s}{\int }}V(X_{\tau },\tau )d\tau }f(X_s,s)ds.$

Aplicado o lema de Itō novamente a $du(X_s,s)$, segue-se que

$\begin{array}{cc}dY=& e^{-\underset{t}{\overset{s}{\int }}V(X_{\tau },\tau )d\tau }(-V(X_s,s)u(X_s,s)+f(X_s,s)+\mu (X_s,s)\frac{\partial u}{\partial X}+\frac{\partial u}{\partial s}+\frac{1}{2}{\sigma }^2(X_s,s)\frac{{\partial }^{2}u}{\partial X^2})ds\\ & +e^{-\underset{t}{\overset{s}{\int }}V(X_{\tau },\tau )d\tau }\sigma (X,s)\frac{\partial u}{\partial X}dW.\end{array}$

O primeiro termo contém, entre parênteses, a EDP acima, sendo portanto zero. O que sobra é

$dY=e^{-\underset{t}{\overset{s}{\int }}V(X_{\tau },\tau )d\tau }\sigma (X,s)\frac{\partial u}{\partial X}dW.$

Integrando esta equação de $t$ a $T$, conclui-se que

$Y(T)-Y(t)=\underset{t}{\overset{T}{\int }}{e}^{-\underset{t}{\overset{s}{\int }}V(X_{\tau },\tau )d\tau }\sigma (X,s)\frac{\partial u}{\partial X}dW.$

Ao tomar valores esperados, condicionados em $X_t=x$, observando que o lado direito é uma integral de Itō, que tem valor esperado zero, segue-se que

$E[Y(T)\mid X_t=x]=E[Y(t)\mid X_t=x]=u(x,t).$

O resultado desejado é obtido observando que

$E[Y(T)\mid X_t=x]=E[e^{-\underset{t}{\overset{T}{\int }}V(X_{\tau },\tau )d\tau }u(X_T,T)+\underset{t}{\overset{T}{\int }}{e}^{-\underset{t}{\overset{r}{\int }}V(X_{\tau },\tau )d\tau }f(X_r,r)dr|X_t=x]$

e finalmente

$u(x,t)=E[e^{-\underset{t}{\overset{T}{\int }}V(X_{\tau },\tau )d\tau }\psi (X_T)+\underset{t}{\overset{T}{\int }}{e}^{-\underset{t}{\overset{s}{\int }}V(X_{\tau },\tau )d\tau }f(X_s,s)ds|X_t=x]$

Considerações

• A prova acima é essencialmente uma demonstração com modificações para computar $f(x,t)$.
• A fórmula de valor esperado acima também é válida para difusões de Itō de $n$ dimensões. A equação diferencial parcial correspondente para $u:{ℝ}^N\times [0,T]\to ℝ$ se torna:

$\frac{\partial u}{\partial t}+\sum _{i=1}^N{\mu }_i(x,t)\frac{\partial u}{\partial x_i}+\frac{1}{2}\sum _{i=1}^N\sum _{j=1}^N{\gamma }_{ij}(x,t)\frac{{\partial }^{2}u}{\partial x_i{x}_j}-r(x,t)u=f(x,t),$

em que,

${\gamma }_{ij}(x,t)=\sum _{k=1}^N{\sigma }_{ik}(x,t){\sigma }_{jk}(x,t),$

isto é, $\gamma =\sigma {\sigma }^{\prime }$, em que ${\sigma }^{\prime }$ denota a matriz transposta de $\sigma )$.

• O valor esperado pode ser então aproximado pelo método de Monte Carlo ou pelo método de quase-Monte Carlo.
• Quando foi publicado pela primeira vez por Kac em 1949, a fórmula de Feynman–Mac foi apresentada como uma fórmula para determinar a distribuição de certas funcionais de Wiener. Suponha que se deseja encontrar o valor esperado da função

$e^{-\underset{0}{\overset{t}{\int }}V(x(\tau ))d\tau }$

no caso em que $x(\tau )$ é alguma ocorrência de um processo de difusão iniciado em $x(0)=0$. A fórmula de Feynman-Kac diz que o valor esperado é equivalente à integral de uma solução para uma equação de difusão. Especificamente, sob as condições de que$uV(x)\ge 0$,

$E\left[e^{-u\underset{0}{\overset{t}{\int }}V(x(\tau ))d\tau }\right]=\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}w(x,t)dx$

em que $w(x,0)=\delta (x)$ e

$\frac{\partial w}{\partial t}=\frac{1}{2}\frac{{\partial }^{2}w}{\partial x^2}-uV(x)w.$

A fórmula de Feynman–Kac também pode ser interpretada como um método para avaliação de integrais funcionais de uma certa forma. Se

$I=\int f(x(0))e^{-u\underset{0}{\overset{t}{\int }}V(x(t))dt}g(x(t))Dx$

em que a integral é tomada sobre todos os passeios aleatórios, então,

$I=\int w(x,t)g(x)dx$

em que $w(x,t)$ é uma solução da equação diferencial parcial parabólica

$\frac{\partial w}{\partial t}=\frac{1}{2}\frac{{\partial }^{2}w}{\partial x^2}-uV(x)w$

com a condição inicial $w(x,0)=f(x)$.

Ver também

• Lema de Itō
• Teorema de Girsanov