Ir al contenido

Fórmula de Feynman-Kac

De Wikipedia, la enciclopedia libre

La fórmula de Feynman-Kac, llamada así por Richard Feynman y Mark Kac, establece un vínculo entre las ecuaciones parabólicas en derivadas parciales (PDE) y los procesos estocásticos. En 1947, cuando Kac y Feynman eran profesores de la Universidad de Cornell, Kac asistió a una conferencia de Feynman y comentó que los dos estaban trabajando en lo mismo desde diferentes direcciones.[1]​ Del intercambio surgió la fórmula de Feynman-Kac, que demostraba rigurosamente el caso real de la integral de caminos de Feynman. El caso complejo, que ocurre cuando se incluye el giro de una partícula, sigue siendo una problema abierto.[2]

La fórmula proporciona un método para resolver ciertas ecuaciones diferenciales parciales simulando los caminos aleatorios de un proceso estocástico. Recíprocamente, una clase importante de valores esperados de procesos aleatorios puede ser calculada por métodos deterministas.

Teorema

[editar]

Considérese la ecuación diferencial parcial

definida para todos y y sujeta a la condición terminal

donde son funciones conocidas, es un parámetro, y es la función incógnita. Entonces la fórmula de Feynman-Kac nos dice que la solución se puede escribir como una expectativa condicional

bajo la medida de probabilidad tal que es un proceso de Itô impulsado por la ecuación

con es un proceso de Wiener (también llamado movimiento browniano) bajo , y la condición inicial para es .

Demostración parcial

[editar]

Una demostración de que la fórmula anterior es una solución de la ecuación diferencial es larga, difícil y no se presenta aquí. Sin embargo, es razonablemente sencillo demostrar que, "si existe una solución", debe tener la forma anterior. Una demostración de este resultado menor sería la siguiente:

Sea la solución a la ecuación diferencial parcial anterior. Aplicación de la regla producto para procesos de Itô al proceso

se obtiene

A partir de:

el tercer término es y se puede eliminar. También tenemos que

Aplicando el lema de Itô a , se deduce que

El primer término contiene, entre paréntesis, la ecuación diferencial parcial anterior y, por lo tanto, es cero. Lo que queda es

Integrando esta ecuación de a , se concluye que

Al tomar expectativas, condicionado a <matemáticas>X_{t} = x</matemáticas>, y observar que el lado derecho es un Itô integral, que tiene expectativa cero,[3]​ de ello se deduce que

El resultado deseado se obtiene observando que

y finalmente

Comentarios

[editar]
  • La demostración anterior de que una solución debe tener la forma dada es esencialmente la de[4]​ con modificaciones para tener en cuenta .
  • La fórmula de expectativa anterior también es válida para las difusiones Itô N-dimensionales. La ecuación diferencial parcial correspondiente para se convierte en:[5] donde es decir, , donde denota la traspuesta de .
  • Esta expectativa puede ser aproximada usando el método de Monte Carlo o método cuasi-Monte Carlos.
  • Cuando fue publicado originalmente por Kac en 1949,[6]​ la fórmula de Feynman-Kac se presentó como una fórmula para determinar la distribución de ciertos funcionales de Wiener. Supongamos que deseamos encontrar el valor esperado de la función en el caso donde x(τ) es alguna realización de un proceso de difusión que comienza en x(0) = 0. La fórmula de Feynman-Kac dice que esta expectativa es equivalente a la integral de una solución a una ecuación de difusión. Específicamente, bajo las condiciones que , donde w(x, 0) = δ(x) y La fórmula de Feynman-Kac también puede interpretarse como un método para evaluar integrales funcionales de cierta forma. Si donde la integral se toma sobre todos paseo aleatorios, entonces donde w(x, t) es una solución a la Ecuación parabólica en derivadas parciales con condición inicial w(x, 0) = f(x).

Aplicaciones

[editar]

En finanzas cuantitativas, la fórmula de Feynman-Kac se utiliza para calcular eficientemente soluciones a la ecuación de Black-Scholes a opciones de precio en acciones[7]​ y bono cupón cero precios en modelos de estructura de términos afines.

En química cuántica, se usa para resolver la ecuación de Schrödinger con el método puro de difusión Monte Carlo.[8]

Véase también

[editar]

Referencias

[editar]
  1. Kac, Mark (1987). Enigmas of Chance: An Autobiography. University of California Press. pp. 115-16. ISBN 0-520-05986-7. 
  2. Glimm, James; Jaffe, Arthur (1987). Quantum Physics: A Functional Integral Point of View (2 edición). New York, NY: Springer. pp. 43-44. ISBN 978-0-387-96476-8. doi:10.1007/978-1-4612-4728-9. Consultado el 13 de abril de 2021. 
  3. Øksendal, Bernt (2003). «Theorem 3.2.1. iii)». En Springer-Verlag, ed. Ecuaciones diferenciales estocásticas. Una introducción con aplicaciones (en inglés) (6º edición). p. 30. ISBN 3540047581. 
  4. «PDE for Finance». 
  5. See Pham, Huyên (2009). Control estocástico de tiempo continuo y optimización con aplicaciones financieras. Springer-Verlag. ISBN 978-3-642-10044-4. 
  6. Kac, Mark (1949). «On Distributions of Certain Wiener Functionals». Transactions of the American Mathematical Society 65 (1): 1-13. JSTOR 1990512. doi:10.2307/1990512.  Este artículo está reimpreso en Baclawski, K.; Donsker, M. D., eds. (1979). Mark Kac: Probabilidad, teoría de números y física estadística, artículos seleccionados. Cambridge, Massachusetts: The MIT Press. pp. 268–280. ISBN 0-262-11067-9. 
  7. Paolo Brandimarte (6 de junio de 2013). «Chapter 1. Motivación». Numerical Methods in Finance and Economics: A MATLAB-Based Introduction. John Wiley & Sons. ISBN 978-1-118-62557-6. 
  8. Caffarel, Michel; Claverie, Pierre (15 de enero de 1988). «Development of a pure diffusion quantum Monte Carlo method using a full generalized Feynman–Kac formula. I. Formalismo». The Journal of Chemical Physics 88 (2): 1088-1099. Bibcode:88.1088C 1988JChPh.. 88.1088C. doi:10.1063/1.454227. 

Bibliografía

[editar]
  • Simon, Barry (1979). Academic Press, ed. Integración funcional y física cuántica. 
  • Hall, B. C. (2013). Springer, ed. Teoría cuántica para matemáticos.