En matemáticas , algunos problemas de valores en la frontera se pueden resolver utilizando los métodos de análisis estocástico . Quizás el ejemplo más célebre sea la solución de Shizuo Kakutani en 1944 al problema de Dirichlet para el operador de Laplace utilizando el movimiento browniano . Sin embargo, resulta que para una gran clase de ecuaciones diferenciales parciales semielípticas de segundo orden, el problema de valor límite de Dirichlet asociado se puede resolver utilizando un proceso de Itō que resuelve una ecuación diferencial estocástica asociada .
Introducción: la solución de Kakutani al problema clásico de Dirichlet
Sea un dominio (un conjunto abierto y conexo ) en . Sea el operador de Laplace , sea una función acotada en la frontera y considere el problema:![{\displaystyle D}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\textstyle \mathbb {R} ^{n}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \Delta}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle\partial D}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\begin{casos}-\Delta u(x)=0,&x\in D\\\displaystyle {\lim _{y\to x}u(y)}=g(x),&x\ en \partial D\end{casos}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Se puede demostrar que si existe una solución, entonces el valor esperado de en el primer punto de salida (aleatorio) es para un movimiento browniano canónico que comienza en . Véase el teorema 3 en Kakutani 1944, p. 710.![{\displaystyle u}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle u(x)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle g(x)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle D}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle x}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
El problema de Dirichlet-Poisson
Sea un dominio in y sea un operador diferencial semielíptico on de la forma:![{\displaystyle D}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\textstyle \mathbb {R} ^{n}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle L}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\textstyle C^{2}(\mathbb {R} ^{n};\mathbb {R} )}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle L=\sum _{i=1}^{n}b_{i}(x){\frac {\partial }{\partial x_{i}}}+\sum _{i,j=1 }^{n}a_{ij}(x){\frac {\partial ^{2}}{\partial x_{i}\,\partial x_{j}}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
donde los coeficientes y son funciones continuas y todos los valores propios de la matriz no son negativos. Deja y . Considere el problema de Poisson :![{\ Displaystyle b_ {i}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \alpha (x)=a_{ij}(x)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\textstyle f\en C(D;\mathbb {R} )}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\textstyle g\in C(\partial D;\mathbb {R} )}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\begin{casos}-Lu(x)=f(x),&x\in D\\\displaystyle {\lim _{y\to x}u(y)}=g(x),&x \in \partial D\end{cases}}\quad {\mbox{(P1)}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
La idea del método estocástico para resolver este problema es la siguiente. Primero, se encuentra una difusión de Itō cuyo generador infinitesimal coincide con funciones soportadas de forma compacta . Por ejemplo, puede considerarse la solución de la ecuación diferencial estocástica:
![{\displaystyle A}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle C^{2}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle f:\mathbb {R} ^{n}\rightarrow \mathbb {R} }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle X}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \mathrm {d} X_{t}=b(X_{t})\,\mathrm {d} t+\sigma (X_{t})\,\mathrm {d} B_{t}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
donde es el movimiento browniano de n dimensiones, tiene componentes como los anteriores y el campo matricial se elige de modo que:![{\displaystyle B}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle b}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \sigma}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\frac {1}{2}}\sigma (x)\sigma (x)^{\top }=a(x),\quad \forall x\in \mathbb {R} ^{n} }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Para un punto , denotamos la ley del dato inicial dado , y denotamos la expectativa con respecto a . Denotemos la primera hora de salida de from .![{\displaystyle x\in \mathbb {R} ^{n}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \mathbb {P} ^{x}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle X}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle X_{0}=x}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \mathbb {E} ^{x}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \mathbb {P} ^{x}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \tau _ {D}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle X}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle D}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
En esta notación, la solución candidata para (P1) es:
![{\displaystyle u(x)=\mathbb {E} ^{x}\left[g{\big (}X_{\tau _ {D}}{\big )}\cdot \chi _{\{\tau _{D}<+\infty \}}\right]+\mathbb {E} ^{x}\left[\int _{0}^{\tau _{D}}f(X_{t})\ ,\mathrm {d} t\right]}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
siempre que sea una función acotada y que:![{\displaystyle g}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \mathbb {E} ^{x}\left[\int _{0}^{\tau _{D}}{\big |}f(X_{t}){\big |}\,\ mathrm {d} t\right]<+\infty }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Resulta que se requiere una condición más:
![{\displaystyle \mathbb {P} ^{x}{\big (}\tau _{D}<\infty {\big )}=1,\quad \forall x\in D}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Para todos , el proceso que comienza casi con seguridad sale en un tiempo finito. Bajo este supuesto, la solución candidata anterior se reduce a:![{\displaystyle x}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle X}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle D}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle u(x)=\mathbb {E} ^{x}\left[g{\big (}X_{\tau _ {D}}{\big )}\right]+\mathbb {E} ^ {x}\left[\int _{0}^{\tau _{D}}f(X_{t})\,\mathrm {d} t\right]}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
y resuelve (P1) en el sentido de que si denota el operador característico para (que concuerda con las funciones), entonces:![{\displaystyle {\mathcal {A}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle X}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle A}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle C^{2}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\begin{casos}-{\mathcal {A}}u(x)=f(x),&x\in D\\\displaystyle {\lim _{t\uparrow \tau _{D}} u(X_{t})}=g{\big (}X_{\tau _{D}}{\big )},&\mathbb {P} ^{x}{\mbox{-as,}}\ ;\forall x\in D\end{cases}}\quad {\mbox{(P2)}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Además, si satisface (P2) y existe una constante tal que, para todos :![{\textstyle v\in C^{2}(D;\mathbb {R} )}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle C}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle x\en D}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle |v(x)|\leq C\left(1+\mathbb {E} ^{x}\left[\int _{0}^{\tau _{D}}{\big |}g (X_{s}){\big |}\,\mathrm {d} s\right]\right)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
entonces .![{\displaystyle v=u}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Referencias
- Kakutani, Shizuo (1944). "Movimiento browniano bidimensional y funciones armónicas". Proc. Diablillo. Acad. Tokio . 20 (10): 706–714. doi : 10.3792/pia/1195572706 .
- Kakutani, Shizuo (1944). "Sobre los movimientos brownianos en el espacio n". Proc. Diablillo. Acad. Tokio . 20 (9): 648–652. doi : 10.3792/pia/1195572742 .
- Øksendal, Bernt K. (2003). Ecuaciones diferenciales estocásticas: una introducción con aplicaciones (Sexta ed.). Berlín: Springer. ISBN 3-540-04758-1.(Ver Sección 9)