Proceso de Wiener - Wiener process

Una realización única de un proceso Wiener unidimensional
Una sola realización de un proceso de Wiener tridimensional

En matemáticas , el proceso de Wiener es un proceso estocástico de tiempo continuo de valor real llamado en honor al matemático estadounidense Norbert Wiener por sus investigaciones sobre las propiedades matemáticas del movimiento browniano unidimensional. A menudo también se le llama movimiento browniano debido a su conexión histórica con el proceso físico del mismo nombre observado originalmente por el botánico escocés Robert Brown . Es uno de los procesos de Lévy más conocidos ( procesos estocásticos de càdlàg con incrementos independientes estacionarios ) y ocurre con frecuencia en matemáticas puras y aplicadas , economía , finanzas cuantitativas , biología evolutiva y física .

El proceso de Wiener juega un papel importante tanto en matemáticas puras como aplicadas. En matemáticas puras, el proceso de Wiener dio lugar al estudio de las martingalas de tiempo continuo . Es un proceso clave en términos del cual se pueden describir procesos estocásticos más complicados. Como tal, juega un papel vital en el cálculo estocástico , los procesos de difusión e incluso la teoría del potencial . Es el proceso impulsor de la evolución de Schramm-Loewner . En matemáticas aplicadas , el proceso de Wiener se utiliza para representar la integral de un proceso gaussiano de ruido blanco y, por lo tanto, es útil como modelo de ruido en ingeniería electrónica (ver ruido browniano ), errores de instrumentos en la teoría de filtrado y perturbaciones en la teoría de control .

El proceso de Wiener tiene aplicaciones en todas las ciencias matemáticas. En física se utiliza para estudiar el movimiento browniano , la difusión de partículas diminutas suspendidas en un fluido y otros tipos de difusión a través de las ecuaciones de Fokker-Planck y Langevin . También forma la base para la rigurosa formulación integral de trayectoria de la mecánica cuántica (mediante la fórmula de Feynman-Kac , una solución a la ecuación de Schrödinger se puede representar en términos del proceso de Wiener) y el estudio de la inflación eterna en cosmología física . También es prominente en la teoría matemática de las finanzas , en particular el modelo de fijación de precios de opciones de Black-Scholes .

Caracterizaciones del proceso Wiener

El proceso Wiener se caracteriza por las siguientes propiedades:

  1. tiene incrementos independientes : para todos los incrementos futuros son independientes de los valores pasados ,
  2. tiene incrementos gaussianos: normalmente se distribuye con media y varianza ,
  3. tiene caminos continuos: es continuo en .

Que el proceso tenga incrementos independientes significa que si 0 ≤ s 1 < t 1s 2 < t 2 entonces W t 1  -  W s 1 y W t 2  -  W s 2 son variables aleatorias independientes, y la condición similar se cumple para n incrementos.

Una caracterización alternativa del proceso de Wiener es la denominada caracterización de Lévy que dice que el proceso de Wiener es casi con seguridad una martingala continua con W 0 = 0 y variación cuadrática [ W t , W t ] = t (lo que significa que W t 2  -  t es también una martingala).

Una tercera caracterización es que el proceso de Wiener tiene una representación espectral como una serie sinusoidal cuyos coeficientes son N (0, 1) variables aleatorias independientes . Esta representación se puede obtener mediante el teorema de Karhunen-Loève .

Otra caracterización de un proceso de Wiener es la integral definida (desde el tiempo cero al tiempo t ) de un proceso gaussiano correlacionado delta ("blanco") de media cero, varianza unitaria .

El proceso de Wiener se puede construir como el límite de escala de una caminata aleatoria u otros procesos estocásticos de tiempo discreto con incrementos independientes estacionarios. Esto se conoce como teorema de Donsker . Al igual que la caminata aleatoria, el proceso de Wiener es recurrente en una o dos dimensiones (lo que significa que casi con seguridad regresa a cualquier vecindario fijo del origen infinitamente a menudo) mientras que no es recurrente en las dimensiones tres y superiores. A diferencia del paseo aleatorio, es invariante en escala , lo que significa que

es un proceso de Wiener para cualquier constante α distinta de cero. La medida de Wiener es la ley de probabilidad en el espacio de funciones continuas g , con g (0) = 0, inducida por el proceso de Wiener. Una integral basada en la medida de Wiener se puede llamar integral de Wiener .

Proceso de salchicha como límite de caminata aleatoria

Vamos a ser iid variables aleatorias con media 0 y varianza 1. Para cada n , definir un proceso estocástico de tiempo continuo

Esta es una función de paso aleatorio. Los incrementos de son independientes porque son independientes. Para n grande , está cerca del teorema del límite central. El teorema de Donsker afirma que a medida , se acerca a un proceso de Wiener, lo que explica la ubicuidad del movimiento browniano.

Propiedades de un proceso de Wiener unidimensional

Propiedades básicas

La función de densidad de probabilidad incondicional , que sigue una distribución normal con media = 0 y varianza = t , en un tiempo fijo t :

La expectativa es cero:

La varianza , usando la fórmula computacional, es t :

Estos resultados se derivan inmediatamente de la definición de que los incrementos tienen una distribución normal , centrada en cero. Por lo tanto

Covarianza y correlación

La covarianza y correlación (donde ):

Estos resultados se derivan de la definición de que los incrementos que no se superponen son independientes, de los cuales solo se utiliza la propiedad de que no están correlacionados. Supongamos eso .

Sustituyendo

llegamos a:

Dado que y son independientes,

Por lo tanto

Un corolario útil para la simulación es que podemos escribir, para t 1 < t 2 :

donde Z es una variable normal estándar independiente.

Representación de Wiener

Wiener (1923) también dio una representación de un camino browniano en términos de una serie aleatoria de Fourier . Si son variables gaussianas independientes con media cero y varianza uno, entonces

y

representar un movimiento browniano en . El proceso escalado

es un movimiento browniano encendido (cf. teorema de Karhunen-Loève ).

Corriendo máximo

La distribución conjunta del máximo corriente

y W t es

Para obtener la distribución incondicional de , integre sobre −∞ < wm  :

la función de densidad de probabilidad de una distribución semi-normal . La expectativa es

Si en un momento el proceso de Wiener tiene un valor conocido , es posible calcular la distribución de probabilidad condicional del máximo en el intervalo (cf. Distribución de probabilidad de los puntos extremos de un proceso estocástico de Wiener ). La función de distribución de probabilidad acumulada del valor máximo, condicionada por el valor conocido , es:

Auto-semejanza

Una demostración de la escala browniana, que muestra la disminución de c . Tenga en cuenta que las características promedio de la función no cambian al hacer zoom, y tenga en cuenta que aumenta cuadráticamente más rápido horizontalmente que verticalmente.

Escala browniana

Para cada c > 0, el proceso es otro proceso de Wiener.

Inversión del tiempo

El proceso para 0 ≤ t ≤ 1 se distribuye como W t para 0 ≤ t ≤ 1.

Inversión de tiempo

El proceso es otro proceso de Wiener.

Una clase de martingalas brownianas

Si un polinomio p ( x , t ) satisface la PDE

luego el proceso estocástico

es una martingala .

Ejemplo: es una martingala, que muestra que la variación cuadrática de W en [0, t ] es igual a t . De ello se deduce que el tiempo esperado de la primera salida de W de (- c , c ) es igual ac 2 .

De manera más general, para cada polinomio p ( x , t ), el siguiente proceso estocástico es una martingala:

donde a es el polinomio

Ejemplo: el proceso

es una martingala, lo que muestra que la variación cuadrática de la martingala en [0, t ] es igual a

Acerca de las funciones p ( xa , t ) más generales que los polinomios, ver martingalas locales .

Algunas propiedades de las rutas de muestra

El conjunto de todas las funciones w con estas propiedades es de la medida Wiener completa. Es decir, una ruta (función de muestra) del proceso de Wiener tiene todas estas propiedades casi con seguridad.

Propiedades cualitativas

  • Para cada ε> 0, la función w toma valores (estrictamente) positivos y (estrictamente) negativos en (0, ε).
  • La función w es continua en todas partes pero no diferenciable en ninguna parte (como la función de Weierstrass ).
  • Los puntos del máximo local de la función w son un conjunto contable denso; los valores máximos son diferentes por pares; cada máximo local es agudo en el siguiente sentido: si w tiene un máximo local en t entonces
Lo mismo ocurre con los mínimos locales.
  • La función w no tiene puntos de incremento local, es decir, ningún t > 0 satisface lo siguiente para algunos ε en (0, t ): primero, w ( s ) ≤ w ( t ) para todos los s en ( t - ε, t ), y en segundo lugar, w ( s ) ≥ w ( t ) para todos los s en ( t , t + ε). (El aumento local es una condición más débil que w aumenta en ( t - ε, t + ε).) Lo mismo es válido para la disminución local.
  • La función w es de variación ilimitada en cada intervalo.
  • La variación cuadrática de w sobre [0, t] es t.
  • Los ceros de la función w son un conjunto perfecto denso en ninguna parte de medida de Lebesgue 0 y dimensión de Hausdorff 1/2 (por lo tanto, incontable).

Propiedades cuantitativas

Ley del logaritmo iterado
Módulo de continuidad

Módulo de continuidad local:

Módulo de continuidad global (Lévy):

Hora local

La imagen de la medida de Lebesgue en [0, t ] debajo del mapa w (la medida de avance ) tiene una densidad L t (·). Por lo tanto,

para una amplia clase de funciones f (a saber: todas las funciones continuas; todas las funciones integrables localmente; todas las funciones mensurables no negativas). La densidad L t es (más exactamente, puede y se elegirá que sea) continua. El número L t ( x ) se denomina hora local en x de w en [0, t ]. Es estrictamente positivo para todo x del intervalo ( a , b ) donde a y b son el menor y el mayor valor de w en [0, t ], respectivamente. (Para x fuera de este intervalo la hora local evidentemente se desvanece.) Trata como una función de dos variables x y t , el tiempo local es todavía continua. Tratada como una función de t (mientras que x es fija), la hora local es una función singular que corresponde a una medida no atómica en el conjunto de ceros de w .

Estas propiedades de continuidad son bastante no triviales. Tenga en cuenta que la hora local también se puede definir (como la densidad de la medida de avance) para una función suave. Entonces, sin embargo, la densidad es discontinua, a menos que la función dada sea monótona. En otras palabras, existe un conflicto entre el buen comportamiento de una función y el buen comportamiento de su hora local. En este sentido, la continuidad del tiempo local del proceso de Wiener es otra manifestación de la no suavidad de la trayectoria.

Tasa de información

La tasa de información del proceso de Wiener con respecto a la distancia de error al cuadrado, es decir, su función cuadrática de tasa-distorsión , viene dada por

Por lo tanto, es imposible codificar usando un código binario de menos de bits y recuperarlo con un error cuadrático medio esperado menor que . Por otro lado, para cualquiera , existe un código binario suficientemente grande de no más de elementos distintos, de modo que el error cuadrático medio esperado en la recuperación de este código es como máximo .

En muchos casos, es imposible codificar el proceso de Wiener sin muestrearlo primero. Cuando el proceso de Wiener se muestrea a intervalos antes de aplicar un código binario para representar estas muestras, la compensación óptima entre la tasa de código y el error cuadrático medio esperado (al estimar el proceso de Wiener de tiempo continuo) sigue la representación paramétrica

donde y . En particular, es el error cuadrático medio asociado solo con la operación de muestreo (sin codificación).

Procesos relacionados

Procesos Wiener con deriva ( azul ) y sin deriva ( rojo ).
Procesos 2D Wiener con deriva ( azul ) y sin deriva ( rojo ).
El generador de un movimiento browniano es ½ veces el operador de Laplace-Beltrami . La imagen de arriba es del movimiento browniano en una variedad especial: la superficie de una esfera.

El proceso estocástico definido por

se denomina proceso de Wiener con deriva μ y varianza infinitesimal σ 2 . Estos procesos agotan los procesos continuos de Lévy .

A grandes rasgos, aparecen dos procesos aleatorios en el intervalo de tiempo [0, 1] al condicionar el proceso de Wiener a desaparecer en ambos extremos de [0,1]. Sin más condicionamiento, el proceso toma valores tanto positivos como negativos en [0, 1] y se denomina puente browniano . Condicionado también para permanecer positivo en (0, 1), el proceso se llama excursión browniana . En ambos casos un tratamiento riguroso implica un procedimiento limitante, ya que la fórmula P ( A | B ) = P ( AB ) / P ( B ) no se aplica cuando P ( B ) = 0.

Se puede escribir un movimiento browniano geométrico

Es un proceso estocástico que se utiliza para modelar procesos que nunca pueden tomar valores negativos, como el valor de las acciones.

El proceso estocástico

se distribuye como el proceso de Ornstein-Uhlenbeck con parámetros , , y .

El tiempo de alcanzar un solo punto x > 0 por el proceso de Wiener es una variable aleatoria con la distribución de Lévy . La familia de estas variables aleatorias (indexadas por todos los números positivos x ) es una modificación continua a la izquierda de un proceso de Lévy . La modificación continua a la derecha de este proceso viene dada por los tiempos de la primera salida de los intervalos cerrados [0, x ].

El tiempo local L = ( L x t ) xR , t ≥ 0 de un movimiento browniano describe el tiempo que el proceso pasa en el punto x . Formalmente

donde δ es la función delta de Dirac . El comportamiento de la hora local se caracteriza por los teoremas de Ray-Knight .

Martingalas brownianas

Sea A un evento relacionado con el proceso de Wiener (más formalmente: un conjunto, medible con respecto a la medida de Wiener, en el espacio de funciones), y X t la probabilidad condicional de A dado el proceso de Wiener en el intervalo de tiempo [0 , t ] (más formalmente: la medida de Wiener del conjunto de trayectorias cuya concatenación con la trayectoria parcial dada en [0, t ] pertenece a A ). Entonces el proceso X t es una martingala continua. Su propiedad de martingala se deriva inmediatamente de las definiciones, pero su continuidad es un hecho muy especial: un caso especial de un teorema general que establece que todas las martingalas brownianas son continuas. Una martingala browniana es, por definición, una martingala adaptada a la filtración browniana; y la filtración browniana es, por definición, la filtración generada por el proceso Wiener.

Movimiento browniano integrado

La integral de tiempo del proceso de Wiener

se llama movimiento browniano integrado o proceso de Wiener integrado . Surge en muchas aplicaciones y puede ser demostrado tener la distribución N (0, t 3 /3), calculado utilizando el hecho de que la covarianza del proceso de Wiener es .

Para el caso general del proceso definido por

Entonces, para ,

De hecho, siempre es una variable aleatoria normal de media cero. Esto permite la simulación de lo dado tomando

donde Z es una variable normal estándar y

El caso de corresponde a . Todos estos resultados pueden verse como consecuencias directas de la isometría de Itô . El proceso de Wiener integrado n veces es una variable normal de media cero con varianza . Esto viene dado por la fórmula de Cauchy para la integración repetida .

Cambio de tiempo

Cada martingala continua (comenzando en el origen) es un proceso de Wiener cambiado en el tiempo.

Ejemplo: 2 W t = V (4 t ) donde V es otro proceso de Wiener (diferente de W pero distribuido como W ).

Ejemplo. donde y V es otro proceso de Wiener.

En general, si M es una martingala continua, entonces donde A ( t ) es la variación cuadrática de M en [0, t ], y V es un proceso de Wiener.

Corolario. (Ver también los teoremas de convergencia de la martingala de Doob ) Sea M t una martingala continua, y

Entonces solo son posibles los siguientes dos casos:

otros casos (como, por ejemplo, etc.) son de probabilidad 0.  

Especialmente, una martingala continua no negativa tiene un límite finito (como t → ∞) casi con seguridad.

Todo lo que se indica (en esta subsección) para martingalas vale también para martingalas locales .

Cambio de medida

Una amplia clase de semimartingales continuos (especialmente, de procesos de difusión ) está relacionada con el proceso de Wiener a través de una combinación de cambio de tiempo y cambio de medida .

Usando este hecho, las propiedades cualitativas indicadas anteriormente para el proceso de Wiener se pueden generalizar a una amplia clase de semimartingalas continuas.

Proceso de Wiener de valor complejo

El proceso de Wiener de valores complejos puede definirse como un proceso aleatorio de valores complejos de la forma donde y son procesos de Wiener independientes (de valor real).

Auto-semejanza

Escala browniana, inversión de tiempo, inversión de tiempo: lo mismo que en el caso de valores reales.

Invarianza de rotación: para cada número complejo de manera que el proceso sea ​​otro proceso de Wiener de valor complejo.

Cambio de tiempo

Si es una función completa, entonces el proceso es un proceso de Wiener de valor complejo modificado en el tiempo.

Ejemplo: donde

y es otro proceso de Wiener de valor complejo.

En contraste con el caso de valor real, una martingala de valor complejo generalmente no es un proceso de Wiener de valor complejo cambiado en el tiempo. Por ejemplo, la martingala no lo es (aquí ya son procesos de Wiener independientes, como antes).

Ver también

Notas

Referencias

enlaces externos