Feynman-Kac formulasi — Richard Feynman va Mark Kac nomi bilan ataladigan formula boʻlib, xususiy hosilali parabolik differensial tenglamalar hamda tasodifiy jarayonlarni oʻzaro bogʻlashga xizmat qiladi.

Xususan, ushbu formula xususiy hosilali differensial tenglamalarni tasodifiy hodisalar trayektoriyasi (Monte-Karlo metodi deb ham ataladi) yordamida yechish imkonini beradi. Shuningdek, tasodifiy jarayonning matematik kutilmasi, xususiy hosilali differensial tenglamaning yechimi sifatida hisoblanishi mumkin.

Bir oʻlchamli holat

Quyidagi differensial tenglamani koʻrib chiqamiz:

∂u∂t+μ(x,t)∂u∂x+12σ2(x,t)∂2u∂x2−V(x,t)u+f(x,t)=0(*){\displaystyle \frac{\partial u}{\partial t} + \mu(x,t) \frac{\partial u}{\partial x} + \tfrac{1}{2} \sigma^2(x,t) \frac{\partial^2 u}{\partial x^2} - V(x,t) u + f(x,t) = 0 \qquad (*)}

bu yerda u=u(x,t){\displaystyle u=u(x,t)} nomaʼlum funksiya, x∈R{\displaystyle x \in \mathbb{R}} va t∈[0,T]{\displaystyle t \in [0,T]} esa — mustaqil oʻzgaruvchilar, μ,σ,V,f{\displaystyle \mu, \sigma, V, f} — maʼlum funksiyalar. Feynman-Kac formulasiga binoan, (*) tenglamaning boshlangʻich shartlarga nisbatan yechimi

u(x,T)=ψ(x),{\displaystyle u(x,T)=\psi(x),}

shartli matematik kutilma koʻrinishida ifodalanishi mumkin:

u(x,t)=EQ[∫tTe−∫tsV(Xτ)dτf(Xs,s)ds+e−∫tTV(Xτ)dτψ(XT)|Xt=x],{\displaystyle u(x,t) = E^Q\left[ \int_t^T e^{- \int_t^s V(X_\tau)\, d\tau}f(X_s,s)ds + e^{-\int_t^T V(X_\tau)\, d\tau}\psi(X_T) \ \bigg| \ X_t=x \right],}

bu yerda Q{\displaystyle Q} — tasodifiy oʻlcham, Xt{\displaystyle X_t} esa Ito protsessi deb ataladi. Ushbu protsess quyidagi tasodifiy tenglama bilan ifodalanadi:

dXt=μ(Xt,t)dt+σ(Xt,t)dWtQ,(**){\displaystyle dX_t = \mu(X_t,t)\,dt + \sigma(X_t,t)\,dW_t^Q, \qquad (**)} bu yerda WtQ{\displaystyle W_t^Q} jarayon uchun boshlangʻich shart quyidagi koʻrinishda boʻladi:

X0=x{\displaystyle X_0=x}.

Koʻp oʻlchamli holat

Feynman-Kac formulasining koʻp oʻlchamli analogi ham mavjud boʻlib, oʻzgaruvchi quyidagi shartga boʻysunsagina uni qoʻllash mumkin: x=(x1,…,xn)∈Rn{\displaystyle x = (x_1, \ldots, x_n) \in \mathbb{R}^n}.

Bu holda (*) differensial tenglama quyidagi koʻrinishga keladi:

∂u∂t+∑i=1nμi(x,t)∂u∂xi+12∑i=1n∑j=1nγij(x,t)∂2u∂xi∂xj−V(x,t)u+f(x,t)=0{\displaystyle \frac{\partial u}{\partial t} + \sum_{i=1}^n \mu_i(x,t) \frac{\partial u}{\partial x_i} + \tfrac{1}{2} \sum_{i=1}^n\sum_{j=1}^n \gamma_{ij}(x,t) \frac{\partial^2 u}{\partial x_i\partial x_j} - V(x,t) u + f(x,t) = 0}

va n-oʻchamli Xt{\displaystyle X_t} tasodifiy jarayon quyidagi tenglama orqali ifodalanadi:

dXt=μ(Xt,t)dt+Σ(Xt,t)dWtQ,{\displaystyle dX_t = \mu(X_t,t)\,dt + \Sigma(X_t,t)\,dW_t^Q,}

bu yerda μ{\displaystyle \mu} — (μ1,…,μn){\displaystyle (\mu_1,\ldots,\mu_n)} elementlardan tuzilgan vektor-ustun, WtQ{\displaystyle W_t^Q} — n-oʻlchamli tasodifiy protsess, Σ=(σij){\displaystyle \Sigma = (\sigma_{ij})} — n-tartibli kvadat matritsa boʻlib, Γ=(γij){\displaystyle \Gamma = (\gamma_{ij})} matritsa bilan quyidagi formula orqali bogʻlangan:

Γ=Σ⋅Σ*,{\displaystyle \Gamma = \Sigma \cdot \Sigma^*,}

yulduzcha (*) belgisi transponirlashni bildiradi.