Aproximación de Poisson Compuesta para Ecuaciones Estocásticas de Volterra con Núcleos Singulares

Información Básica

• ID del Artículo: 2510.27323
• Título: Compound Poisson Approximation for Stochastic Volterra Equations with Singular Kernels
• Autores: Xicheng Zhang (Universidad Tecnológica de Beijing), Yuanlong Zhao (Universidad de Wuhan)
• Clasificación: math.PR (Teoría de Probabilidades)
• Fecha de Presentación: 31 de octubre de 2025
• Enlace del Artículo: https://arxiv.org/abs/2510.27323

Resumen

Este artículo establece la convergencia fuerte de ecuaciones diferenciales estocásticas (EDEs) y ecuaciones estocásticas de tipo Volterra bajo aproximación de procesos de Poisson compuestos, derivando tasas de convergencia explícitas. En comparación con el método clásico de Euler-Maruyama, la ventaja clave de la aproximación de Poisson compuesta es que no requiere continuidad del coeficiente de deriva en la variable temporal, e incluso puede manejar problemas de singularidad. Experimentos numéricos verifican la estabilidad del método.

Antecedentes de Investigación y Motivación

1. Problema Central

Este artículo estudia el problema de aproximación numérica de ecuaciones diferenciales estocásticas y ecuaciones estocásticas de Volterra, con énfasis especial en casos donde los coeficientes presentan singularidades temporales.

2. Importancia del Problema

• Significado Teórico: Las ecuaciones estocásticas de Volterra ocupan un lugar importante en el análisis estocástico, particularmente en la descripción de sistemas estocásticos con efectos de memoria
• Valor Aplicado: Ecuaciones estocásticas impulsadas por movimiento Browniano fraccional, ecuaciones estocásticas de Navier-Stokes y otros modelos importantes pueden reducirse a este tipo de problemas
• Desafío Numérico: Cuando las funciones núcleo presentan singularidades, los métodos numéricos tradicionales enfrentan problemas de convergencia y estabilidad

3. Limitaciones de Métodos Existentes

• Método Euler-Maruyama: Requiere continuidad Hölder de los coeficientes en la variable temporal, incapaz de manejar casos discontinuos o singulares
• Métodos EM Mejorados: Aunque relajan ciertas condiciones (como la condición one-sided Lipschitz), aún requieren cierta regularidad temporal
• Métodos para Núcleos Singulares: La literatura existente se enfoca principalmente en núcleos no singulares; la teoría para manejar núcleos singulares no es suficientemente completa

4. Motivación de la Investigación

• Desarrollar métodos numéricos capaces de manejar singularidades temporales
• Establecer estimaciones de tasas de convergencia explícitas
• Proporcionar herramientas numéricas confiables para procesos estocásticos importantes como el movimiento Browniano fraccional

Contribuciones Principales

1. Establecimiento de Teoría de Convergencia Fuerte para Aproximación de Poisson Compuesta: Primera demostración sistemática de convergencia fuerte de EDEs y ecuaciones estocásticas de Volterra bajo aproximación de Poisson compuesta
2. Derivación de Tasas de Convergencia Explícitas:
   • Para EDEs estándar: tasa de convergencia $\varepsilon^{\gamma \wedge \frac{\beta}{2}}$
   • Para ecuaciones estocásticas de Volterra: tasa de convergencia $\varepsilon^{\gamma/(2(2+\gamma))}$
3. Manejo Innovador de Coeficientes Singulares: El método no requiere continuidad temporal del coeficiente de deriva y puede manejar singularidades temporales
4. Provisión de Marco Teórico Completo: Incluye estimaciones de momentos, aplicación de desigualdades de Gronwall, lemas técnicos y sistema completo de pruebas
5. Presentación de Ejemplos de Aplicación Específicos: Verificación detallada de condiciones de hipótesis para movimiento Browniano fraccional (parámetro de Hurst $H \in (0,1/2) \cup (1/2,1)$)
6. Verificación Numérica: Dos experimentos numéricos verifican la efectividad y estabilidad del método

Explicación Detallada del Método

Definición de Tareas

Problema 1: EDE Estándar Considérese la ecuación diferencial estocástica en $\mathbb{R}^d$: $X_t = X_0 + \int_0^t \sigma(s, X_s) dW_s + \int_0^t b(s, X_s) ds$

donde $\sigma: \mathbb{R}_+ \times \mathbb{R}^d \to \mathbb{R}^d \otimes \mathbb{R}^m$, $b: \mathbb{R}_+ \times \mathbb{R}^d \to \mathbb{R}^d$, y $W_t$ es un movimiento Browniano estándar $m$-dimensional.

Problema 2: Ecuación Estocástica de Volterra$Y_t = Y_0 + \int_0^t \sigma(t, s, Y_s) dW_s + \int_0^t b(t, s, Y_s) ds$

La característica clave es que los coeficientes dependen del tiempo actual $t$ y de la variable de integración $s$, reflejando efectos de memoria.

Arquitectura del Modelo

1. Construcción de Proceso de Poisson Compuesto

Proceso de Poisson Base: Sea $(T_k)_{k \in \mathbb{N}}$ una secuencia de variables aleatorias exponenciales independientes e idénticamente distribuidas con parámetro 1: $P(T_k \geq t) = e^{-t}, \quad t \geq 0$

Defínanse los tiempos de salto: $S_0 = 0$, $S_n = S_{n-1} + T_n$

Proceso de Poisson estándar: $N_t = \max\{n: S_n \leq t\}$

Proceso de Poisson Escalado: Para $\varepsilon > 0$ dado, defínase: $N_t^\varepsilon := \varepsilon N_{t/\varepsilon}, \quad \tilde{N}_t^\varepsilon := N_t^\varepsilon - t$

donde $N_t^\varepsilon$ es un proceso de Poisson con tamaños de salto $\varepsilon$ e intensidad $1/\varepsilon$.

Proceso de Poisson Compuesto: $W_{N_t^\varepsilon} = \sum_{k=1}^{N_{t/\varepsilon}} (W_{k\varepsilon} - W_{(k-1)\varepsilon})$

Este es un proceso de Poisson compuesto con medida de intensidad $\nu_\varepsilon \sim N(0, \varepsilon I_d)$.

2. Esquema Numérico

Aproximación de Poisson Compuesta para EDE: $X_t^\varepsilon = X_0 + \int_0^t \sigma(s, X_{s-}^\varepsilon) dW_{N_s^\varepsilon} + \int_0^t b(s, X_{s-}^\varepsilon) dN_s^\varepsilon$

Forma completamente discreta: $X_t^\varepsilon = X_0 + \sum_{k=1}^{N_{t/\varepsilon}} [\sigma(S_k^\varepsilon, X_{S_{k-1}^\varepsilon}^\varepsilon)(W_{k\varepsilon} - W_{(k-1)\varepsilon}) + \varepsilon b(S_k^\varepsilon, X_{S_{k-1}^\varepsilon}^\varepsilon)]$

Aproximación de Poisson Compuesta para Ecuación Estocástica de Volterra: $Y_t^\varepsilon = Y_0 + \int_0^t \sigma(t, s, Y_{s-}^\varepsilon) dW_{N_s^\varepsilon} + \int_0^t b(t, s, Y_{s-}^\varepsilon) dN_s^\varepsilon$

Forma completamente discreta: $Y_t^\varepsilon = Y_0 + \sum_{k=1}^{N_{t/\varepsilon}} [\sigma(t, S_k^\varepsilon, Y_{S_{k-1}^\varepsilon}^\varepsilon)(W_{k\varepsilon} - W_{(k-1)\varepsilon}) + \varepsilon b(t, S_k^\varepsilon, Y_{S_{k-1}^\varepsilon}^\varepsilon)]$

Puntos de Innovación Técnica

1. Discretización Temporal Aleatoria

A diferencia de las mallas de tiempo determinísticas, se utilizan los tiempos de salto del proceso de Poisson como puntos de discretización, constituyendo una estrategia de discretización aleatoria con ventajas clave:

• Introduce aleatoriedad de forma natural, independiente de los incrementos del movimiento Browniano
• No requiere continuidad temporal de los coeficientes
• Puede manejar singularidades temporales

2. Estrategia de Descomposición de Error

Para la EDE, el error se descompone en cuatro términos: $X_t - X_t^\varepsilon = I_1(t) + I_2(t) + I_3(t) + I_4(t)$

donde:

• $I_1(t)$: Diferencia entre la trayectoria verdadera y su valor en puntos de tiempo aleatorio
• $I_2(t)$: Diferencia entre movimiento Browniano y proceso de Poisson compuesto (término más crítico)
• $I_3(t)$: Diferencia de coeficientes en diferentes trayectorias
• $I_4(t)$: Error de discretización del término de deriva

3. Lemas Técnicos Clave

Lema 2.1 (Estimaciones de Momentos para Tiempos de Salto de Poisson): Para $\alpha \in (0,1]$, $\beta \geq 0$: $\mathbb{E}|r^\alpha - (S_k^\varepsilon)^\alpha|^\beta \lesssim (k\varepsilon)^{\alpha\beta} k^{-\beta/2}, \quad r \in [(k-1)\varepsilon, k\varepsilon]$

Esta estimación es más refinada que usar directamente $|a^\alpha - b^\alpha| \leq |a-b|^\alpha$.

Lema 2.3 (Estimación de Desviación del Proceso de Poisson): Para cualquier $p > 0$: $\mathbb{E}|N_t^\varepsilon - t|^p \leq C\varepsilon^{p/2}(t^{p/2} \vee t)$

4. Tratamiento Refinado del Término $I_2(t)$

Este es el término más difícil, requiriendo:

• Representación mediante medida aleatoria de Poisson compensada
• Desigualdad máxima de Doob
• Hipótesis de continuidad Hölder en variable temporal $(H_\sigma^t)$
• Estimaciones de suma refinadas

Para ecuaciones de Volterra, el tratamiento de $I_2$ es aún más complejo, requiriendo introducir un parámetro de discretización intermedio $\delta$ y descomponerlo en cuatro subtérminos.

Configuración Experimental

Conjunto de Datos (Escenarios de Experimentos Numéricos)

Experimento 1: EDE Lineal con Deriva de Dependencia Temporal Singular

$dX_t = \mu(t)X_t dt + \sigma(t)X_t dW_t$

donde: