2025-11-22T07:28:16.184474

Existence of Invariant Measures for Delay Equations with Stochastic Negative Feedback

Bosch, van Gaans, Lunel
We provide sufficient conditions for the existence of invariant probability measures for generic stochastic differential equations with finite time delay. Applications include the Mackey--Glass equations and Nicholson's blowflies equation, each perturbed by a (small) multiplicative noise term. Solutions to these stochastic negative feedback systems persist globally and all solutions are bounded above in probability. It turns out that the occurrence of finite time blowups and boundedness in probability of solutions and solution segments are closely related. A non-trivial invariant measure is shown to exist if and only if there is at least one initial condition for which the solution remains bounded away from zero in probability. The noise driving the dynamical system is allowed to be an integrable Lévy process.
academic

Existencia de Medidas Invariantes para Ecuaciones con Retardo y Retroalimentación Negativa Estocástica

Información Básica

  • ID del Artículo: 2501.00141
  • Título: Existence of Invariant Measures for Delay Equations with Stochastic Negative Feedback
  • Autores: M. van den Bosch (Universidad de Leiden), O.W. van Gaans (Universidad de Leiden), S.M. Verduyn Lunel (Universidad de Utrecht)
  • Clasificación: math.DS (Sistemas Dinámicos), math.PR (Teoría de Probabilidad)
  • Fecha de Publicación: 3 de enero de 2025
  • Enlace del Artículo: https://arxiv.org/abs/2501.00141

Resumen

Este artículo proporciona condiciones suficientes para la existencia de medidas de probabilidad invariantes para ecuaciones diferenciales estocásticas generales con retardo temporal finito. Las aplicaciones incluyen la ecuación de Mackey-Glass y la ecuación de Nicholson's blowflies perturbadas por términos de ruido multiplicativo (pequeño). Las soluciones de estos sistemas estocásticos de retroalimentación negativa existen globalmente y están acotadas en sentido probabilístico. Se encuentra que la ocurrencia de explosión en tiempo finito está estrechamente relacionada con la acotación probabilística de las soluciones y sus segmentos de trayectoria. Existen medidas invariantes no triviales si y solo si existe al menos una condición inicial tal que la solución se mantiene alejada de cero en sentido probabilístico. El ruido que impulsa el sistema dinámico puede ser un proceso de Lévy integrable.

Antecedentes de Investigación y Motivación

Descripción del Problema

Este artículo estudia la existencia de medidas invariantes para ecuaciones diferenciales con retardo que poseen retroalimentación negativa estocástica. Los objetos de estudio central son ecuaciones diferenciales con retardo deterministas de la forma:

x'(t) = -γ(t)x(t) + r(t)f(x(t-τ))

y sus versiones con perturbaciones estocásticas.

Importancia de la Investigación

  1. Significado Biológico: La ecuación de Mackey-Glass se utiliza para modelar la concentración de glóbulos blancos, la ecuación de Nicholson's blowflies describe la dinámica de poblaciones de insectos, estos modelos tienen aplicaciones importantes en fisiología y ecología
  2. Desafíos Matemáticos: El comportamiento caótico de ecuaciones diferenciales con retardo aún no ha sido demostrado rigurosamente, la búsqueda de medidas invariantes proporciona nuevas vías para comprender el comportamiento a largo plazo de estos sistemas
  3. Perturbaciones Estocásticas: Los sistemas reales inevitablemente están sujetos a ruido, el estudio de versiones estocásticas es más cercano a las aplicaciones prácticas

Limitaciones de Métodos Existentes

  • Para la ecuación de Mackey-Glass y la ecuación de Nicholson's blowflies, los fenómenos caóticos aún no han sido demostrados matemáticamente
  • Los resultados existentes sobre atractores globales solo se aplican a regiones de parámetros que permiten soluciones estacionarias y periódicas
  • Falta un enfoque sistemático para tratar medidas invariantes de ecuaciones diferenciales estocásticas con retardo general

Contribuciones Principales

  1. Establecimiento de un Marco Teórico General: Se proporcionan condiciones suficientes para la existencia de medidas invariantes para ecuaciones diferenciales estocásticas con retardo temporal finito
  2. Demostración de la Existencia Global de Soluciones: Todas las soluciones existen globalmente y están acotadas en sentido probabilístico
  3. Caracterización de Condiciones Necesarias y Suficientes para la Existencia de Medidas Invariantes: Existen medidas invariantes no triviales si y solo si existe al menos una condición inicial tal que la solución se mantiene alejada de cero en sentido probabilístico
  4. Extensión de Tipos de Ruido: Se permite que el ruido sea un proceso de Lévy integrable, no limitado únicamente al movimiento browniano
  5. Aplicación a Modelos Clásicos: Se demuestra la existencia de medidas invariantes para versiones estocásticas de la ecuación de Mackey-Glass y la ecuación de Nicholson's blowflies

Explicación Detallada de Métodos

Definición de la Tarea

Se estudian ecuaciones diferenciales estocásticas con retardo autónomas:

dX(t) = a(X_t) dt + b(X_{t-}) dM(t), t ≥ 0
X(u) = Φ(u), u ∈ [-τ, 0]

donde M(t) es una semimartingala, a, b son funcionales localmente Lipschitz, y X_t denota el proceso de segmento de la solución en el tiempo t.

Arquitectura del Método Principal

1. Técnica de Cambio de Variables

Mediante transformación logarítmica Y(t) = log X(t), la ecuación original se transforma en:

dY(t) = [-γ(t) + r(t)e^{-Y(t)}f(e^{Y(t-τ)})] dt + a(Y_t,t) dt + b(Y_t,t) dW(t)

Esta transformación preserva la positividad de las soluciones, facilitando el análisis.

2. Método de Estimación de Trayectorias

  • Estimación de Cota Superior: Utilizando la propiedad de retroalimentación negativa, cuando la solución es suficientemente grande, el término de deriva se vuelve negativo, previniendo que la solución tienda a infinito
  • Estimación de Cota Inferior: Cuando f(0) > 0, se utilizan técnicas de desigualdades integrales para obtener cotas inferiores de la solución

3. Análisis de Acotación Probabilística

Se establecen conexiones entre tres conceptos clave:

  • Existencia global
  • Acotación probabilística de soluciones
  • Acotación probabilística de segmentos de trayectoria

4. Condiciones de Compacidad

Mediante el teorema de Arzelà-Ascoli y la topología de Skorokhod, se establecen condiciones suficientes para la compacidad de segmentos de solución:

  • Condición (i): Coeficientes de ruido acotados
  • Condición (ii): El supremo del término de deriva está acotado en sentido probabilístico
  • Condición (iii): La norma del supremo de segmentos de solución está acotada en sentido probabilístico

Puntos de Innovación Técnica

1. Técnica de Supremo en Tiempo Inverso

Para procesos con deriva negativa impulsados por Lévy:

Y(t) = -∫₀ᵗ a(s) ds + ∫₀ᵗ b(s) dL(s)

Se demuestra que el proceso de supremo en tiempo inverso (sup_{0≤θ≤t}(Y(t)-Y(θ)))_{t≥0} está acotado en sentido probabilístico.

2. Extensión del Método de Krylov-Bogoliubov

Se extiende el teorema de existencia clásico de Krylov-Bogoliubov al espacio de funciones continuas por la derecha con límites por la izquierda D-τ,0, manejando el caso de ruido con saltos.

3. Aplicación de Características de Semimartingalas

Se utilizan las características locales de semimartingalas (B^X, C^X, ν^X) para demostrar la compacidad del proceso de segmento de solución, aplicable a tipos de ruido más amplios.

Configuración Experimental

Verificación Teórica

Este trabajo es principalmente teórico, verificando la corrección de los resultados mediante demostración matemática, en lugar de experimentos numéricos.

Modelos de Aplicación

  1. Ecuación de Mackey-Glass:
    x'(t) = -γx(t) + r·x(t-1)/(1+x(t-1)^p), p > 1
    
  2. Ecuación de Nicholson's blowflies:
    x'(t) = -γx(t) + r·x(t-1)e^{-px(t-1)}, p > 0
    

Tipos de Ruido

  • Movimiento browniano
  • Procesos de Lévy con saltos acotados
  • Procesos de Lévy integrables

Resultados Principales

Teorema 1.1 (Resultado Principal)

Supóngase que f: ℝ → ℝ es localmente Lipschitz continua, no negativa en (0,∞) y acotada superiormente, inf_{t≥0} γ(t) > 0, sup_{t≥0} r(t) < ∞. Entonces la solución de la ecuación diferencial estocástica con retardo:

dX(t) = [-γ(t)X(t) + r(t)f(X(t-τ))] dt + σX(t)c(X_t) dW(t)

es única, existe globalmente, y para casi todas las condiciones iniciales F₀-medibles no negativas está acotada en sentido probabilístico.

Además, cuando γ y r son constantes:

  • Si f(0) > 0, entonces casi todas las soluciones se mantienen alejadas de cero en sentido probabilístico, existiendo una distribución estacionaria μ con soporte en (0,∞)
  • Si f(0) = 0, entonces existen dos distribuciones estacionarias distintas: la medida de Dirac δ₀ y una distribución estacionaria μ con soporte en (0,∞) (siempre que exista al menos una condición inicial tal que la solución se mantenga alejada de cero en sentido probabilístico)

Corolario 1.2 (Ecuación de Mackey-Glass)

Para la ecuación de Mackey-Glass determinista, cuando p > 1, inf_{t∈-τ,0} φ(t) > 0, lim inf_{t→∞} r(t)/γ(t) > 1, la solución se mantiene alejada de cero y acotada. Cuando γ y r son constantes y r > γ, la ecuación posee una medida invariante no trivial en C-1,0.

Trabajo Relacionado

Teoría de Ecuaciones Diferenciales con Retardo

  • La teoría clásica de Hale & Verduyn Lunel proporciona el marco fundamental para este trabajo
  • El trabajo relacionado con la ecuación de Wright proporciona precedentes para perturbaciones estocásticas
  • La investigación sobre existencia de atractores globales proporciona referencias para el análisis de comportamiento asintótico

Ecuaciones Diferenciales Estocásticas

  • Teoría de ecuaciones diferenciales estocásticas de Mao
  • Teoría de integración de semimartingalas de Protter
  • Trabajo de Da Prato & Zabczyk sobre ergodicidad de sistemas de dimensión infinita

Teoría de Medidas Invariantes

  • Teorema de existencia clásico de Krylov-Bogoliubov
  • Investigación de propiedades de Feller para sistemas con retardo finito
  • Teoría de compacidad y ergodicidad

Conclusiones y Discusión

Conclusiones Principales

  1. Se establece una teoría general para la existencia de medidas invariantes en ecuaciones de retardo con retroalimentación negativa estocástica
  2. Se demuestra la existencia global de soluciones y su acotación probabilística
  3. Se caracterizan las condiciones necesarias y suficientes para la existencia de medidas invariantes no triviales
  4. Se proporciona una base matemática rigurosa para modelos clásicos de biomatemática

Limitaciones

  1. Caso f(0) = 0: Cuando f(0) = 0, la demostración de que las soluciones se mantienen alejadas de cero en sentido probabilístico sigue siendo un problema abierto
  2. Unicidad: La unicidad de la medida invariante no ha sido demostrada
  3. Aspecto Computacional: Falta un método efectivo para calcular numéricamente medidas invariantes
  4. Generalización de Dimensión Superior: La extensión del método a sistemas requiere trabajo adicional

Direcciones Futuras

  1. Investigar estimaciones de cota inferior de soluciones cuando f(0) = 0
  2. Explorar unicidad, regularidad y continuidad absoluta de medidas invariantes
  3. Desarrollar métodos numéricos para calcular medidas invariantes
  4. Extender a sistemas de retardo multidimensionales
  5. Investigar estabilidad asintótica y velocidad de convergencia

Evaluación Profunda

Fortalezas

  1. Innovación Teórica: La extensión exitosa del método de Krylov-Bogoliubov a ecuaciones diferenciales estocásticas con retardo demuestra una dificultad técnica muy alta
  2. Sistematicidad del Método: Se establece una cadena teórica completa desde la existencia local hasta la existencia global, desde la acotación probabilística hasta la compacidad, y luego a la existencia de medidas invariantes
  3. Valor de Aplicación: Proporciona una base teórica rigurosa para modelos importantes de biomatemática
  4. Profundidad Técnica: La aplicación de técnicas de supremo en tiempo inverso y características de semimartingalas refleja un nivel técnico superior

Deficiencias

  1. Completitud: El tratamiento del caso f(0) = 0 no es suficientemente completo, representando una deficiencia teórica importante
  2. Practicidad: Falta verificación numérica y métodos computacionales, dificultando la verificación de resultados teóricos en la práctica
  3. Generalidad: El método se aplica principalmente a ecuaciones escalares, la extensión a sistemas requiere trabajo adicional

Impacto

  1. Contribución Académica: Realiza una contribución importante a la teoría de ecuaciones diferenciales estocásticas con retardo, llenando un vacío teórico en el campo
  2. Perspectivas de Aplicación: Proporciona nuevas herramientas para el análisis de modelos estocásticos en biología, ecología y otros campos
  3. Valor de Métodos: Las técnicas desarrolladas tienen valor de referencia importante para investigaciones relacionadas

Escenarios Aplicables

Los métodos de este artículo son aplicables a:

  • Modelado de sistemas biológicos con mecanismos de retroalimentación negativa
  • Análisis de sistemas dinámicos de retardo perturbados por ruido
  • Sistemas estocásticos donde se requiere investigar el comportamiento a largo plazo
  • Ecuaciones diferenciales con retardo que poseen ruido multiplicativo

Referencias

El artículo cita 88 referencias importantes que abarcan múltiples campos incluyendo ecuaciones diferenciales con retardo, análisis estocástico, teoría ergódica y otros, reflejando una base teórica sólida y una perspectiva académica amplia.