Articulo de referencia

Ecuación diferencial con retardo

En matemáticas , las ecuaciones diferenciales con retardo ( EDR ) son un tipo de ecuación diferencial en la que la derivada de la función desconocida en un momento dado se expre...

En matemáticas , las ecuaciones diferenciales con retardo ( EDR ) son un tipo de ecuación diferencial en la que la derivada de la función desconocida en un momento dado se expresa en términos de los valores de la función en momentos anteriores. Las EDR también se denominan sistemas con retardo temporal , sistemas con efecto residual o tiempo muerto, sistemas hereditarios, ecuaciones con argumento desviado o ecuaciones diferenciales en diferencias. Pertenecen a la clase de sistemas con un estado funcional , es decir, ecuaciones diferenciales parciales (EDP) que son de dimensión infinita, a diferencia de las ecuaciones diferenciales ordinarias (EDO) que tienen un vector de estado de dimensión finita. Cuatro puntos pueden dar una posible explicación de la popularidad de las EDR: [ 1 ]

  1. El efecto posterior es un problema aplicado: es bien sabido que, junto con las crecientes expectativas de rendimiento dinámico, los ingenieros necesitan que sus modelos se comporten de manera más similar al proceso real. Muchos procesos incluyen fenómenos de efecto posterior en su dinámica interna. Además, los actuadores , sensores y redes de comunicación que ahora participan en bucles de control de retroalimentación introducen dichos retrasos. Finalmente, además de los retrasos reales, los retardos temporales se utilizan con frecuencia para simplificar modelos de orden muy alto. Por lo tanto, el interés por las ecuaciones diferenciales con retardo (EDR) sigue creciendo en todas las áreas científicas y, especialmente, en la ingeniería de control.
  2. Los sistemas con retardo siguen siendo resistentes a muchos controladores clásicos : cabría pensar que la solución más sencilla consistiría en sustituirlos por aproximaciones de dimensión finita. Desafortunadamente, ignorar efectos que se representan adecuadamente mediante ecuaciones diferenciales con retardo no es una alternativa general: en el mejor de los casos (retardos constantes y conocidos), conlleva el mismo grado de complejidad en el diseño del control. En los peores casos (retardos variables en el tiempo, por ejemplo), puede resultar desastroso en términos de estabilidad y oscilaciones.
  3. La introducción voluntaria de retardos puede beneficiar al sistema de control . [ 2 ]
  4. A pesar de su complejidad, las EDD a menudo aparecen como modelos simples de dimensión infinita en el área muy compleja de las ecuaciones diferenciales parciales (EDP).

Una forma general de la ecuación diferencial con retardo de tiempo paraincógnita(t)Rnorte{\displaystyle x(t)\in \mathbb {R} ^{n}}es ddtincógnita(t)=F(t,incógnita(t),incógnitat),{\displaystyle {\frac {d}{dt}}x(t)=f(t,x(t),x_{t}),} dóndeincógnitat={incógnita(τ):τt}{\displaystyle x_{t}=\{x(\tau ):\tau \leq t\}}representa la trayectoria de la solución en el pasado. En esta ecuación,F{\displaystyle f}es un operador funcional deR×Rnorte×do1(R,Rnorte){\displaystyle \mathbb {R} \times \mathbb {R} ^{n}\times C^{1}(\mathbb {R} ,\mathbb {R} ^{n})}aRnorte.{\displaystyle \mathbb {R} ^{n}.}

Ejemplos

  • Retraso continuoddtincógnita(t)=F(t,incógnita(t),0incógnita(t+τ)dμ(τ)){\displaystyle {\frac {d}{dt}}x(t)=f{\left(t,x(t),\int _{-\infty }^{0}x(t+\tau )\,d\mu (\tau )\right)}}
  • Retardo discretoddtincógnita(t)=F(t,incógnita(t),incógnita(tτ1),,incógnita(tτmetro)){\displaystyle {\frac {d}{dt}}x(t)=f(t,x(t),x(t-\tau _{1}),\dots ,x(t-\tau _{m}))}paraτ1>>τmetro0.{\displaystyle \tau _{1}>\dots >\tau _{m}\geq 0.}
  • Lineal con retardos discretosddtincógnita(t)=A0incógnita(t)+A1incógnita(tτ1)++Ametroincógnita(tτmetro){\displaystyle {\frac {d}{dt}}x(t)=A_{0}x(t)+A_{1}x(t-\tau _{1})+\dots +A_{m}x(t-\tau _{m})}dóndeA0,,AmetroRnorte×norte{\displaystyle A_{0},\dotsc ,A_{m}\in \mathbb {R} ^{n\times n}}.
  • Ecuación del pantógrafoddtincógnita(t)=aincógnita(t)+bincógnita(λt),{\displaystyle {\frac {d}{dt}}x(t)=ax(t)+bx(\lambda t),}donde a , b y λ son constantes y 0 < λ < 1. Esta ecuación y algunas formas más generales reciben su nombre de los pantógrafos de los trenes. [ 3 ] [ 4 ]

Resolución de ecuaciones diferenciales con retardo

Las ecuaciones diferenciales con retardo (EDR) se resuelven principalmente de forma gradual mediante un principio denominado método de pasos. Por ejemplo, consideremos la EDR con un único retardo. ddtincógnita(t)=F(incógnita(t),incógnita(tτ)){\displaystyle {\frac {d}{dt}}x(t)=f(x(t),x(t-\tau ))}

con la condición inicial dadaϕ:[τ,0]Rnorte{\displaystyle \phi \colon [-\tau ,0]\to \mathbb {R} ^{n}}. Entonces la solución en el intervalo[0,τ]{\displaystyle [0,\tau ]}es dado porψ(t){\displaystyle \psi (t)}que es la solución al problema de valor inicial no homogéneoddtψ(t)=F(ψ(t),ϕ(tτ)),{\displaystyle {\frac {d}{dt}}\psi (t)=f(\psi (t),\phi (t-\tau )),} conψ(0)=ϕ(0){\displaystyle \psi (0)=\phi (0)}Esto puede continuarse para los intervalos sucesivos utilizando la solución del intervalo anterior como término no homogéneo. En la práctica, el problema de valor inicial se suele resolver numéricamente.

Ejemplo

SuponerF(incógnita(t),incógnita(tτ))=aincógnita(tτ){\displaystyle f(x(t),x(t-\tau ))=ax(t-\tau )}yϕ(t)=1{\displaystyle \phi (t)=1}. Entonces el problema de valor inicial se puede resolver con integración,

incógnita(t)=incógnita(0)+s=0tddtincógnita(s)ds=1+as=0tϕ(sτ)ds,{\displaystyle x(t)=x(0)+\int _{s=0}^{t}{\frac {d}{dt}}x(s)\,ds=1+a\int _{s=0}^{t}\phi (s-\tau )\,ds,}

es decir,incógnita(t)=at+1{\displaystyle x(t)=at+1}donde la condición inicial viene dada porincógnita(0)=ϕ(0)=1{\displaystyle x(0)=\phi (0)=1}. De manera similar, para el intervalo t[τ,2τ]{\displaystyle t\in [\tau ,2\tau ]}nos integramos y ajustamos a la condición inicial,

incógnita(t)=incógnita(τ)+s=τtddtincógnita(s)ds=(aτ+1)+as=τt(a(sτ)+1)ds=(aτ+1)+as=0tτ(as+1)ds,{\displaystyle {\begin{aligned}x(t)=x(\tau )+\int _{s=\tau }^{t}{\frac {d}{dt}}x(s)\,ds&=(a\tau +1)+a\int _{s=\tau }^{t}\left(a(s-\tau )+1\right)ds\\&=(a\tau +1)+a\int _{s=0}^{t-\tau }\left(as+1\right)ds,\end{aligned}}}

es decir,incógnita(t)=(aτ+1)+a(tτ)(12a(tτ)+1).{\textstyle x(t)=(a\tau +1)+a(t-\tau )\left({\frac {1}{2}}{a(t-\tau )}+1\right).}

Reducción a EDO

En algunos casos, las ecuaciones diferenciales pueden representarse en un formato que se parece a las ecuaciones diferenciales con retardo .

  • Ejemplo 1 Consideremos una ecuaciónddtincógnita(t)=F(t,incógnita(t),0incógnita(t+τ)miλτdτ).{\displaystyle {\frac {d}{dt}}x(t)=f\left(t,x(t),\int _{-\infty }^{0}x(t+\tau )e^{\lambda \tau }\,d\tau \right).}Introduciry(t)=0incógnita(t+τ)miλτdτ{\displaystyle y(t)=\int _{-\infty }^{0}x(t+\tau )e^{\lambda \tau }\,d\tau }para obtener un sistema de EDOddtincógnita(t)=F(t,incógnita,y),ddty(t)=incógnitaλy.{\displaystyle {\frac {d}{dt}}x(t)=f(t,x,y),\quad {\frac {d}{dt}}y(t)=x-\lambda y.}
  • Ejemplo 2 Una ecuaciónddtincógnita(t)=F(t,incógnita(t),0incógnita(t+τ)porque(ατ+β)dτ){\displaystyle {\frac {d}{dt}}x(t)=f\left(t,x(t),\int _{-\infty }^{0}x(t+\tau )\cos(\alpha \tau +\beta )\,d\tau \right)}es equivalente addtincógnita(t)=F(t,incógnita,y),ddty(t)=porque(β)incógnita+αz,ddtz(t)=pecado(β)incógnitaαy,{\displaystyle {\frac {d}{dt}}x(t)=f(t,x,y),\quad {\frac {d}{dt}}y(t)=\cos(\beta )x+\alpha z,\quad {\frac {d}{dt}}z(t)=\sin(\beta )x-\alpha y,}dóndey=0incógnita(t+τ)porque(ατ+β)dτ,z=0incógnita(t+τ)pecado(ατ+β)dτ.{\displaystyle y=\int _{-\infty }^{0}x(t+\tau )\cos(\alpha \tau +\beta )\,d\tau ,\quad z=\int _{-\infty }^{0}x(t+\tau )\sin(\alpha \tau +\beta )\,d\tau .}

La ecuación característica

De forma similar a las EDO , muchas propiedades de las EDD lineales pueden caracterizarse y analizarse utilizando la ecuación característica . [ 5 ] La ecuación característica asociada a la EDD lineal con retardos discretos ddtincógnita(t)=A0incógnita(t)+A1incógnita(tτ1)++Ametroincógnita(tτmetro){\displaystyle {\frac {d}{dt}}x(t)=A_{0}x(t)+A_{1}x(t-\tau _{1})+\dots +A_{m}x(t-\tau _{m})}es el polinomio exponencial dado pordet(λI+A0+A1miτ1λ++Ametromiτmetroλ)=0.{\displaystyle \det(-\lambda I+A_{0}+A_{1}e^{-\tau _{1}\lambda }+\dotsb +A_{m}e^{-\tau _{m}\lambda })=0.}

Las raíces λ de la ecuación característica se denominan raíces características o autovalores, y el conjunto de soluciones se conoce a menudo como espectro . Debido a la exponencial en la ecuación característica, la EDO tiene, a diferencia del caso de la EDO, un número infinito de autovalores, lo que hace que el análisis espectral sea más complejo. Sin embargo, el espectro posee algunas propiedades que pueden aprovecharse en el análisis. Por ejemplo, aunque existe un número infinito de autovalores, solo hay un número finito de autovalores en cualquier franja vertical del plano complejo. [ 6 ]

Esta ecuación característica es un problema de valores propios no lineal y existen muchos métodos para calcular el espectro numéricamente. [ 7 ] [ 8 ] En algunas situaciones especiales es posible resolver la ecuación característica explícitamente. Consideremos, por ejemplo, la siguiente EDD: ddtincógnita(t)=incógnita(t1).{\displaystyle {\frac {d}{dt}}x(t)=-x(t-1).} La ecuación característica es λmiλ=0.{\displaystyle -\lambda -e^{-\lambda }=0.} Esta ecuación tiene un número infinito de soluciones para λ complejo . Están dadas por λ=Wk(1),{\displaystyle \lambda =W_{k}(-1),} donde W k es la k -ésima rama de la función W de Lambert , por lo tanto: incógnita(t)=incógnita(0)miWk(1)t.{\displaystyle x(t)=x(0)\,e^{W_{k}(-1)\cdot t}.}

Otro ejemplo

La siguiente DDE: [ 9 ]ddt(t)=2(2t+1)2(2t1).{\displaystyle {\frac {d}{dt}}u(t)=2u(2t+1)-2u(2t-1).}

Tener como solución enR{\displaystyle \mathbb {R} }la función: [ 10 ](t)={F(t+1),|t|<10,|t|1{\displaystyle u(t)={\begin{cases}F(t+1),&|t|<1\\0,&|t|\geq 1\end{cases}}}conF(t){\displaystyle F(t)}la función de Fabius , conocida como función de Rvachëv up . [ 11 ]

Aplicaciones

Véase también

Referencias

  1. Richard, Jean-Pierre (2003). "Sistemas de retardo de tiempo: una visión general de algunos avances recientes y problemas abiertos". Automatica . 39 (10): 1667– 1694. doi : 10.1016/S0005-1098(03)00167-5 .
  2. Lavaei, Javad; Sojoudi, Somayeh ; Murray, Richard M. (2010). "Implementación simple basada en retardo de controladores de tiempo continuo" . Actas de la Conferencia Americana de Control de 2010. págs. 5781–5788 . doi : 10.1109/ACC.2010.5530439 . ISBN  978-1-4244-7427-1. S2CID 1200900 . 
  3. Griebel, Thomas (2017-01-01). "La ecuación del pantógrafo en el cálculo cuántico" . Tesis de maestría .
  4. Ockendon, John Richard; Tayler, AB; Temple, George Frederick James (1971-05-04). "La dinámica de un sistema de recolección de corriente para una locomotora eléctrica" . Actas de la Royal Society de Londres. Serie A, Ciencias Matemáticas y Físicas . 322 (1551): 447– 468. Bibcode : 1971RSPSA.322..447O . doi : 10.1098/rspa.1971.0078 . S2CID 110981464 . 
  5. Michiels, Wim; Niculescu, Silviu-Iulian (2007). Estabilidad y estabilización de sistemas con retardo de tiempo . Avances en diseño y control. Sociedad de Matemáticas Industriales y Aplicadas. pp. 3–32 . doi : 10.1137/1.9780898718645 . ISBN  978-0-89871-632-0.
  6. Michiels, Wim; Niculescu, Silviu-Iulian (2007). Estabilidad y estabilización de sistemas con retardo de tiempo . Avances en diseño y control. Sociedad de Matemáticas Industriales y Aplicadas. pág. 9. doi : 10.1137/1.9780898718645 . ISBN  978-0-89871-632-0.
  7. Michiels, Wim; Niculescu, Silviu-Iulian (2007). Estabilidad y estabilización de sistemas con retardo de tiempo . Avances en diseño y control. Sociedad de Matemáticas Industriales y Aplicadas. pp. 33–56 . doi : 10.1137/1.9780898718645 . ISBN  978-0-89871-632-0.
  8. Appeltans, Pieter; Michiels, Wim (2023-04-29). "Análisis y diseño de controladores de sistemas con retardo de tiempo utilizando TDS-CONTROL. Un tutorial y manual". arXiv : 2305.00341 [ math.OC ].
  9. Juan Arias de Reyna (2017). "Aritmética de la función de Fabius". arXiv : 1702.06487 [ math.NT ].
  10. "A288163 - Oeis" .
  11. "A288163 - Oeis" .
  12. Makroglou, Athena; Li, Jiaxu; Kuang, Yang (1 de marzo de 2006). "Modelos matemáticos y herramientas de software para el sistema regulador glucosa-insulina y la diabetes: una visión general" . Matemáticas Numéricas Aplicadas . Artículos seleccionados, Tercera Conferencia Internacional sobre Soluciones Numéricas de Ecuaciones de Volterra y Retardo. 56 (3): 559– 573. doi : 10.1016/j.apnum.2005.04.023 . ISSN 0168-9274 . 
  13. Salpeter, Edwin E.; Salpeter, Shelley R. (1998-02-15). "Modelo matemático para la epidemiología de la tuberculosis, con estimaciones del número reproductivo y la función de retardo de la infección" . American Journal of Epidemiology . 147 (4): 398– 406. doi : 10.1093/oxfordjournals.aje.a009463 . ISSN 0002-9262 . PMID 9508108 .  
  14. Kajiwara, Tsuyoshi; Sasaki, Toru; Takeuchi, Yasuhiro (2012-08-01). "Construcción de funcionales de Lyapunov para ecuaciones diferenciales con retardo en virología y epidemiología" . Nonlinear Analysis: Real World Applications . 13 (4): 1802– 1826. doi : 10.1016/j.nonrwa.2011.12.011 . ISSN 1468-1218 . 
  15. Gopalsamy, K. (1992). Estabilidad y oscilaciones en ecuaciones diferenciales con retardo de la dinámica de poblaciones . Matemáticas y sus aplicaciones. Dordrecht, Países Bajos: Kluwer Academic Publishers. doi : 10.1007/978-94-015-7920-9 . ISBN 978-0792315940.
  16. Kuang, Y. (1993). Ecuaciones diferenciales con retardo y aplicaciones en dinámica de poblaciones . Matemáticas en ciencia e ingeniería. San Diego, CA: Academic Press. ISBN 978-0080960029.
  17. López, Álvaro G. (2020-09-01). "Sobre un origen electrodinámico de las fluctuaciones cuánticas" . Dinámica no lineal . 102 (1): 621– 634. arXiv : 2001.07392 . Bibcode : 2020NonDy.102..621L . doi : 10.1007/s11071-020-05928-5 . ISSN 1573-269X . S2CID 210838940 .  

Lecturas adicionales

  • Bellen, Alfredo; Zennaro, Marino (2003). Métodos numéricos para ecuaciones diferenciales con retardo . Matemáticas numéricas y computación científica. Oxford, Reino Unido: Oxford University Press. ISBN 978-0198506546.
  • Bellman, Richard; Cooke, Kenneth L. (1963). Ecuaciones diferenciales en diferencias (PDF) . Matemáticas en ciencia e ingeniería. Nueva York, NY: Academic Press. ISBN 978-0120848508.{{cite book}}: Incompatibilidad de ISBN/Fecha ( ayuda )
  • Briat, Corentin (2015). Sistemas lineales con parámetros variables y retardo de tiempo: análisis, observación, filtrado y control . Avances en retardos y dinámica. Heidelberg, Alemania: Springer-Verlag. ISBN 978-3662440490.
  • Driver, Rodney D. (1977). Ecuaciones diferenciales ordinarias y con retardo . Ciencias matemáticas aplicadas. Vol.  20. Nueva York, NY: Springer-Verlag. doi : 10.1007/978-1-4684-9467-9 . ISBN 978-0387902319.
  • Erneux, Thomas (2009). Ecuaciones diferenciales con retardo aplicadas . Encuestas y tutoriales en ciencias matemáticas aplicadas. Vol.  3. Nueva York, NY: Springer Science+Business Media. doi : 10.1007/978-0-387-74372-1 . ISBN 978-0387743714.
  • Skip Thompson (ed.). "Ecuaciones diferenciales con retardo" . Scholarpedia .