En teoría de la probabilidad y estadística , los procesos de difusión son una clase de procesos de Markov de tiempo continuo con trayectorias de muestra casi seguramente continuas . Los procesos de difusión son estocásticos por naturaleza y, por lo tanto, se utilizan para modelar muchos sistemas estocásticos de la vida real. El movimiento browniano , el movimiento browniano reflejado y los procesos de Ornstein-Uhlenbeck son ejemplos de procesos de difusión. Se utilizan ampliamente en física estadística , análisis estadístico , teoría de la información , ciencia de datos , redes neuronales , finanzas y marketing .
Un modelo de trayectoria de difusión representa la trayectoria de una partícula incrustada en un fluido en movimiento y sometida a desplazamientos aleatorios debido a colisiones con otras partículas, lo que se conoce como movimiento browniano . La posición de la partícula es entonces aleatoria; su función de densidad de probabilidad en función del espacio y el tiempo se rige por una ecuación de convección-difusión .
Definición matemática Un proceso de difusión es un proceso de Markov con trayectorias de muestra continuas para el cual la ecuación directa de Kolmogorov es la ecuación de Fokker-Planck . [ 1 ]
Un proceso de difusión se define por las siguientes propiedades. Seaa i j ( incógnita , t ) {\displaystyle a^{ij}(x,t)} sean coeficientes uniformemente continuos yb i ( incógnita , t ) {\displaystyle b^{i}(x,t)} ser acotado, términos de deriva medibles de Borel. Hay una familia única de medidas de probabilidadPAG a ; b ξ , τ {\displaystyle \mathbb {P} _{a;b}^{\xi ,\tau }} (paraτ ≥ 0 {\displaystyle \tau \geq 0} ,ξ ∈ R d {\displaystyle \xi \in \mathbb {R} ^{d}} ) en el espacio canónicoΩ = do ( [ 0 , ∞ ) , R d ) {\displaystyle \Omega =C([0,\infty ),\mathbb {R} ^{d})} , con su Borelσ {\displaystyle \sigma } -álgebra, tal que:
1. (Condición inicial) El proceso comienza enξ {\displaystyle \xi } en ese momentoτ {\displaystyle \tau } :PAG a ; b ξ , τ [ ψ ∈ Ω : ψ ( t ) = ξ para 0 ≤ t ≤ τ ] = 1. {\displaystyle \mathbb {P} _{a;b}^{\xi ,\tau }[\psi \in \Omega :\psi (t)=\xi {\text{ para }}0\leq t\leq \tau ]=1.}
2. (Propiedad Martingala Local) Por cadaF ∈ do 2 , 1 ( R d × [ τ , ∞ ) ) {\displaystyle f\in C^{2,1}(\mathbb {R} ^{d}\times [\tau ,\infty ))} , el proceso
METRO t [ F ] = F ( ψ ( t ) , t ) − F ( ψ ( τ ) , τ ) − ∫ τ t ( L a ; b + ∂ ∂ s ) F ( ψ ( s ) , s ) d s {\displaystyle M_{t}^{[f]}=f(\psi (t),t)-f(\psi (\tau ),\tau )-\int _{\tau }^{t}{\bigl (}L_{a;b}+{\tfrac {\partial }{\partial s}}{\bigr )}f(\psi (s),s)\,ds} es una martingala local bajoPAG a ; b ξ , τ {\displaystyle \mathbb {P} _{a;b}^{\xi ,\tau }} parat ≥ τ {\displaystyle t\geq \tau } , conMETRO t [ F ] = 0 {\displaystyle M_{t}^{[f]}=0} parat ≤ τ {\displaystyle t\leq \tau } .
Esta familiaPAG a ; b ξ , τ {\displaystyle \mathbb {P} _{a;b}^{\xi ,\tau }} se llama elL a ; b {\displaystyle {\mathcal {L}}_{a;b}} -difusión.
Construcción SDE y generador infinitesimal Está claro que si tenemos unL a ; b {\displaystyle {\mathcal {L}}_{a;b}} -difusión, es decir( incógnita t ) t ≥ 0 {\displaystyle (X_{t})_{t\geq 0}} en( Ω , F , F t , PAG a ; b ξ , τ ) {\displaystyle (\Omega ,{\mathcal {F}},{\mathcal {F}}_{t},\mathbb {P} _{a;b}^{\xi ,\tau })} , entoncesincógnita t {\displaystyle X_{t}} satisface el SDEd incógnita t i = 1 2 ∑ k = 1 d σ k i ( incógnita t ) d B t k + b i ( incógnita t ) d t {\displaystyle dX_{t}^{i}={\frac {1}{2}}\,\sum _{k=1}^{d}\sigma _{k}^{i}(X_{t})\,dB_{t}^{k}+b^{i}(X_{t})\,dt} . Por el contrario, se puede construir esta difusión a partir de esa EDE sia i j ( incógnita , t ) = ∑ k σ i k ( incógnita , t ) σ j k ( incógnita , t ) {\displaystyle a^{ij}(x,t)=\sum _{k}\sigma _{i}^{k}(x,t)\,\sigma _{j}^{k}(x,t)} yσ i j ( incógnita , t ) {\displaystyle \sigma ^{ij}(x,t)} ,b i ( incógnita , t ) {\displaystyle b^{i}(x,t)} son Lipschitz continuas. Para ver esto, dejemosincógnita t {\displaystyle X_{t}} Resuelva la EDE comenzando enincógnita τ = ξ {\displaystyle X_{\tau }=\xi } . ParaF ∈ do 2 , 1 ( R d × [ τ , ∞ ) ) {\displaystyle f\in C^{2,1}(\mathbb {R} ^{d}\times [\tau ,\infty ))} , aplica la fórmula de Itô:d F ( incógnita t , t ) = ( ∂ F ∂ t + ∑ i = 1 d b i ∂ F ∂ incógnita i + v ∑ i , j = 1 d a i j ∂ 2 F ∂ incógnita i ∂ incógnita j ) d t + ∑ i , k = 1 d ∂ F ∂ incógnita i σ k i d B t k . {\displaystyle df(X_{t},t)={\bigl (}{\frac {\partial f}{\partial t}}+\sum _{i=1}^{d}b^{i}{\frac {\partial f}{\partial x_{i}}}+v\sum _{i,j=1}^{d}a^{ij}\,{\frac {\partial ^{2}f}{\partial x_{i}\partial x_{j}}}{\bigr )}\,dt+\sum _{i,k=1}^{d}{\frac {\partial f}{\partial x_{i}}}\,\sigma _{k}^{i}\,dB_{t}^{k}.} Reorganizar daF ( incógnita t , t ) − F ( incógnita τ , τ ) − ∫ τ t ( ∂ F ∂ s + L a ; b F ) d s = ∫ τ t ∑ i , k = 1 d ∂ F ∂ incógnita i σ k i d B s k , {\displaystyle f(X_{t},t)-f(X_{\tau },\tau )-\int _{\tau }^{t}{\bigl (}{\frac {\partial f}{\partial s}}+L_{a;b}f{\bigr )}\,ds=\int _{\tau }^{t}\sum _{i,k=1}^{d}{\frac {\partial f}{\partial x_{i}}}\,\sigma _{k}^{i}\,dB_{s}^{k},} cuyo lado derecho es una martingala local, que coincide con la propiedad de martingala local en la definición de difusión. La ley deincógnita t {\displaystyle X_{t}} definePAG a ; b ξ , τ {\displaystyle \mathbb {P} _{a;b}^{\xi ,\tau }} enΩ = do ( [ 0 , ∞ ) , R d ) {\displaystyle \Omega =C([0,\infty ),\mathbb {R} ^{d})} con la condición inicial correcta y la propiedad de martingala local. La unicidad se deriva de la continuidad de Lipschitz deσ , b {\displaystyle \sigma \!,\!b} . De hecho,L a ; b + ∂ ∂ s {\displaystyle L_{a;b}+{\tfrac {\partial }{\partial s}}} coincide con el generador infinitesimalA {\displaystyle {\mathcal {A}}} de este proceso. Siincógnita t {\displaystyle X_{t}} resuelve la EDE, luego paraF ( incógnita , t ) ∈ do 2 ( R d × R + ) {\displaystyle f(\mathbf {x} ,t)\in C^{2}(\mathbb {R} ^{d}\times \mathbb {R} ^{+})} , el generadorA {\displaystyle {\mathcal {A}}} esA F ( incógnita , t ) = ∑ i = 1 d b i ( incógnita , t ) ∂ F ∂ incógnita i + v ∑ i , j = 1 d a i j ( incógnita , t ) ∂ 2 F ∂ incógnita i ∂ incógnita j + ∂ F ∂ t . {\displaystyle {\mathcal {A}}f(\mathbf {x} ,t)=\sum _{i=1}^{d}b_{i}(\mathbf {x} ,t)\,{\frac {\partial f}{\partial x_{i}}}+v\sum _{i,j=1}^{d}a_{ij}(\mathbf {x} ,t)\,{\frac {\partial ^{2}f}{\partial x_{i}\partial x_{j}}}+{\frac {\partial f}{\partial t}}.}
Referencias ↑ "9. Procesos de difusión" (PDF) . Consultado el 10 de octubre de 2011 .