9.2. Equação de Fokker-Planck para equações diferenciais estocásticas

Considere, agora, uma equação diferencial estocástica

dXt=f(t,Xt) dt+g(t,Xt) dWt, \mathrm{d}X_t = f(t, X_t)\;\mathrm{d}t + g(t, X_t)\;\mathrm{d}W_t,

com uma condição initial

Xtt=0=X0. \left. X_t \right|_{t = 0} = X_0.

Suponhamos, para todos os efeitos, que exista uma solução definida pelo menos em um intervalo de tempo [0,T].[0, T]. Suponhamos, ainda, que cada XtX_t tenha uma função densidade de probabilidade p(t,x)p(t, x) bem definida. A questão que queremos investigar é sobre a evolução dessa função.

Vamos considerar funções testes Φ(t,x)\Phi(t, x) que se anulem nos extremos do intervalo temporal, i.e. Φ(T,x)=Φ(0,x)=0,\Phi(T, x) = \Phi(0, x) = 0, e que decaiam suficientemente rápido quando x.|x|\rightarrow \infty.

Temos, pela fórmula de Itô, que

Φ(t,XT)=Φ(0,X0)+0T(Φt(t,Xt)+Φx(t,Xt)f(t,Xt)+12Φxx(t,Xt)g(t,Xt)2) dt+0TΦx(t,Xt)g(t,Xt) dWt. \Phi(t, X_T) = \Phi(0, X_0) + \int_0^T \left(\Phi_t(t, X_t) + \Phi_x(t, X_t)f(t, X_t) + \frac{1}{2}\Phi_{xx}(t, X_t)g(t, X_t)^2\right)\;\mathrm{d}t + \int_0^T \Phi_x(t, X_t)g(t, X_t)\;\mathrm{d}W_t.

Como Φ\Phi se anula em t=0t=0 e t=T,t=T, obtemos

0T(Φt(t,Xt)+Φx(t,Xt)f(t,Xt)+12Φxx(t,Xt)g(t,Xt)2) dt+0TΦx(t,Xt)g(t,Xt) dWt=0. \int_0^T \left(\Phi_t(t, X_t) + \Phi_x(t, X_t)f(t, X_t) + \frac{1}{2}\Phi_{xx}(t, X_t)g(t, X_t)^2\right)\;\mathrm{d}t + \int_0^T \Phi_x(t, X_t)g(t, X_t)\;\mathrm{d}W_t = 0.

Tomando o valor esperado, chegamos a

0TE[Φt(t,Xt)+Φx(t,Xt)f(t,Xt)+12Φxx(t,Xt)g(t,Xt)2] dt=0. \int_0^T \mathbb{E}\left[\Phi_t(t, X_t) + \Phi_x(t, X_t)f(t, X_t) + \frac{1}{2}\Phi_{xx}(t, X_t)g(t, X_t)^2\right]\;\mathrm{d}t = 0.

Em termos da função densidade de probabilidade,

0T(Φt(t,x)+Φx(t,x)f(t,x)+12Φxx(t,x)g(t,x)2)p(t,x) dx dt=0. \int_0^T \int_{-\infty}^{\infty} \left(\Phi_t(t, x) + \Phi_x(t, x)f(t, x) + \frac{1}{2}\Phi_{xx}(t, x)g(t, x)^2\right)p(t, x)\;\mathrm{d}x\;\mathrm{d}t = 0.

Integrando por partes, em t,t, a primeira integral, temos

0TΦt(t,x)p(t,x) dx dt=0TΦt(t,x)p(t,x) dt dx=(Φ(T,x)Φ(T,0)0TΦ(t,x)pt(t,x) dt) dx. \begin{align*} \int_0^T \int_{-\infty}^{\infty} \Phi_t(t, x) p(t, x)\;\mathrm{d}x\;\mathrm{d}t & = \int_{-\infty}^{\infty} \int_0^T \Phi_t(t, x) p(t, x)\;\mathrm{d}t \;\mathrm{d}x \\ & = \int_{-\infty}^{\infty} \left( \Phi(T, x) - \Phi(T, 0) - \int_0^T \Phi(t, x) p_t(t, x)\;\mathrm{d}t \right) \;\mathrm{d}x. \end{align*}

Novamente, como Φ(t,x)\Phi(t, x) se anula em t=0t=0 e t=T,t=T, isso nos dá

0TΦt(t,x)p(t,x) dx dt=0TΦ(t,x)pt(t,x) dt dx. \int_0^T \int_{-\infty}^{\infty} \Phi_t(t, x) p(t, x)\;\mathrm{d}x\;\mathrm{d}t = -\int_{-\infty}^{\infty} \int_0^T \Phi(t, x) \frac{\partial p}{\partial t}(t, x)\;\mathrm{d}t \;\mathrm{d}x.

Integrando as outras duas integrais por partes em x,x, temos

0T(Φx(t,x)f(t,x)+12Φxx(t,x)g(t,x)2)p(t,x) dx dt=0T(Φ(t,x)x(f(t,x)p(t,x))+12Φ(t,x)2x2(g(t,x)2p(t,x))) dx dt. \begin{align*} \int_0^T \int_{-\infty}^{\infty} & \left(\Phi_x(t, x)f(t, x) + \frac{1}{2}\Phi_{xx}(t, x)g(t, x)^2\right)p(t, x)\;\mathrm{d}x\;\mathrm{d}t \\ & = \int_0^T \int_{-\infty}^{\infty} \left(-\Phi(t, x)\frac{\partial}{\partial x}(f(t, x)p(t, x)) + \frac{1}{2}\Phi(t, x)\frac{\partial^2}{\partial x^2}(g(t, x)^2 p(t, x))\right)\;\mathrm{d}x\;\mathrm{d}t. \end{align*}

Juntando os termos, chegamos a

0T(pt(t,x)x(f(t,x)p(t,x))+122x2(g(t,x)2p(t,x)))Φ(t,x) dx dt. \int_0^T \int_{-\infty}^{\infty} \left(-\frac{\partial p}{\partial t}(t, x) -\frac{\partial}{\partial x}(f(t, x)p(t, x)) + \frac{1}{2}\frac{\partial^2}{\partial x^2}(g(t, x)^2 p(t, x))\right)\Phi(t, x)\;\mathrm{d}x\;\mathrm{d}t.

Como a função teste é arbitrária, chegamos na equação de Fokker-Planck para a evolução da distribuição de probabilidades

pt(t,x)+x(f(t,x)p(t,x))122x2(g(t,x)2p(t,x))=0. \frac{\partial p}{\partial t}(t, x) + \frac{\partial}{\partial x}(f(t, x)p(t, x)) - \frac{1}{2}\frac{\partial^2}{\partial x^2}(g(t, x)^2 p(t, x)) = 0.

Essa também é conhecida como equações progressivas de Kolmogorov para equações diferenciais estocásticas, ou forward Kolmogorov diffusion equation. As equações de Kolmogorov (progressivas ou regressivas - forward or backward) são definidas em contextos mais gerais de processos de Markov com tempo contínuo e espaço discreto ou contínuo.

Exercícios

  1. Deduza a equação de Fokker-Planck no caso multi-dimensional

dXt=f(t,Xt) dt+G(t,Xt) dWt, \mathrm{d}\mathbf{X}_t = f(t, \mathbf{X}_t)\;\mathrm{d}t + G(t, \mathbf{X}_t)\;\mathrm{d}\mathbf{W}_t,

em que Xt\mathbf{X}_t é um processo com valores em Rd,\mathbb{R}^d, f:R×RdRd,f:\mathbb{R} \times \mathbb{R}^d \rightarrow \mathbb{R}^d, G:R×RdRd×k,G:\mathbb{R} \times \mathbb{R}^d \rightarrow \mathbb{R}^{d\times k}, e {Wt}t0\{\mathbf{W}_t\}_{t\geq 0} é um processo com valores em Rk\mathbb{R}^k composto de kk processos de Wiener independentes, com d,kNd, k \in\mathbb{N} arbitrários.

Last modified: June 01, 2026. Built with Franklin.jl, using the Book Template.