APROX - NUMERlCA DE EDR

U NIVERSIDAD T ECNOLÓGICA DE P EREIRA
T ESIS DE M AESTRÍA
APROXIMACIÓN NUMÉRICA DE
ECUACIONES DIFERENCIALES CON
RETARDO POR EL MÉTODO DE
TRANSFORMACIÓN DIFERENCIAL
APLICADO A MODELOS BIOLÓGICOS
Autor: Director de tesis:

JHON FREDY TRUJILLO Dr. PEDRO PABLO
VALENCIA CÁRDENAS ALZATE
Trabajo presentado como requisito para optar

al título de MAGISTER EN ENSEÑANZA DE LA MATEMÁTICA
bajo la supervisión del
Grupo de Investigación en Ecuaciones Diferenciales No Lineales GEDNOL

Departamento de Matemáticas
Facultad de Ciencias Básicas
Universidad Tecnológica de Pereira
11 de agosto de 2016
I
Agradecimientos
Llegar a esta instancia en mis estudios de maestría ha sido posible obviamente a mi

esfuerzo y dedicación; sin embargo se requirió de alguna motivación externa y de algunos
acompañantes en el proceso que fueron tan importantes, y que posiblemente sin ellos hu-
biese sido imposible realizar este trabajo.
Iniciaré agradeciendo a mi familia, especialmente a mis padres por su incondicional

apoyo y constante ánimo en mi lucha por salir adelante. Ellos, junto con las ganas de triun-
far, fueron el motor que mantuvieron constante mi duro trasegar por este camino.
Agradezco a mi gran mentor, el Doctor Pedro Pablo Cárdenas Alzate por todo su gran
apoyo, por brindarme la posibilidad de aprender a su lado. Sus sugerencias, consejos y toda
la paciencia que tuvo conmigo en el desarrollo de este proyecto fue en definitiva la clave del
éxito, además de la gran motivación en el aprendizaje del tema estudiado.
Agradezco también a cada uno de los profesores que participaron en mi formación, en

especial al Doctor José Rodrigo González Granada , quien me motivó para entrar al fasci-
nante mundo de las ecuaciones diferenciales. Gracias por su paciencia y sus consejos en mi
proceso formativo.
Por último doy gracias también a cada uno de mis compañeros que cada fin de semana,
por teléfono o por redes sociales me ayudaron con mi formación y me motivaron a no decli-
nar. Ellos hicieron que este proceso tan complejo fuera además de interesante muy divertido.
Jhon Fredy Trujillo Valencia

II
Índice general
1. Introducción y generalidades 1
1.1. Resumen . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.2. Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.3. Objetivos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.3.1. General . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.3.2. Específicos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.4. Justificación . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.4.1. Planteamiento del problema . . . . . . . . . . . . . . . . . . . . . . . . 4
1.5. Metodología . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
2. Generalidades de las ecuaciones diferenciales con retardo 5

2.1. Ecuaciones Diferenciales con retardo (EDR) . . . . . . . . . . . . . . . . . . . . 5
2.2. Modelo matemático con retardo: Caso especial . . . . . . . . . . . . . . . . . . 6
2.3. Sobre la existencia y unicidad de soluciones . . . . . . . . . . . . . . . . . . . . 8
3. Método de Transformación Diferencial DTM 10

3.1. Antecedentes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
3.2. Generalidades del DTM . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
3.2.1. Transformación Diferencial . . . . . . . . . . . . . . . . . . . . . . . . . 11
4. DTM para ecuaciones diferenciales con retardo 16

4.1. Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
4.2. Extensión del DTM a ecuaciones con retardo . . . . . . . . . . . . . . . . . . . 16
4.3. Aplicaciones a ecuaciones con retardo . . . . . . . . . . . . . . . . . . . . . . . 22
5. Aplicación del DTM a algunos modelos Biológicos 30

5.1. Modelo para crecimiento de tumores en ratones . . . . . . . . . . . . . . . . . 30
5.1.1. Caso particular . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
5.2. Diseminación de infecciones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
5.2.1. Caso particular . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
5.3. Dinámica poblacional . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
5.3.1. Modelo de predicción . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
5.3.2. Caso particular 1. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
5.3.3. Caso particular 2. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
III
... a mis padres

1
Capítulo 1
Introducción y generalidades
1.1. Resumen
La interpretación aproximada de algunos fenómenos naturales ha llevado a introducir en
ciertos tipos de ecuaciones diferenciales cambios en la variable temporal llamados retardos,
lo cual hace que dichas ecuaciones y sus soluciones tengan un comportamiento más acorde
con la realidad. Estas ecuaciones, denominadas ecuaciones diferenciales con retardo requieren
de métodos complejos para su solución y en la mayoría de casos, tan sólo se logra realizar
una aproximación numérica.[11]
Una de las técnicas utilizadas en los últimos años para aproximar soluciones a modelos
de este tipo es el llamado método de transformación diferencial, más conocido por sus siglas en
inglés como DTM (Differential Transformation Method), método cuya estructura está deter-
minada por su definición y una serie de teoremas, los cuales son aplicados al problema en
cuestión. [9, 10]
En esta investigación mostraremos inicialmente un desarrollo teórico sobre el método de

transformación diferencial aplicado a ecuaciones diferenciales ordinarias con retardo en el
que se estudiarán las propiedades más importantes. Posteriormente se realizarán las demos-
traciones de los teoremas más relevantes en el DTM; se resolverán algunos modelos biológi-
cos que involucran ecuaciones diferenciales con retardo y ecuaciones integro-diferenciales.
Finalmente, se realizará la comparación numérica con otros métodos de aproximación y se
analizará la convergencia de dicho método en algunas soluciones.
1.2. Introducción
El Método de Transformación Diferencial tiene suma importancia en la resolución de
ecuaciones diferenciales no lineales (EDOs y EDPs) con valor inicial y en la frontera, con
la necesidad de conocer métodos que son generalmente de tipo semi-analítico numérico
utilizando soluciones en forma de serie de potencias. La transformación diferencial se apli-
ca a funciones de variable real, siendo así su argumento el orden k−ésimo de la derivada. [1]
De manera general, este método consiste en convertir una ecuación diferencial de domi-
nio real en una más simple (ecuación iterativa) de dominio natural. Dicha transformación se
realiza de igual manera que otra clase de transformadas, como es el caso de la transformada
de Laplace. Se emplea la definición y una serie de teoremas que se aplican dependiendo de
la característica de cada uno de los términos de la ecuación diferencial a transformar. En
Capítulo 1. Introducción y generalidades 2
esencia, este método tiene su razón de ser en las series de Taylor, ya que se busca llegar a
una serie infinita que represente la aproximación de la solución de la ecuación.
Muchos autores han investigado sobre el problema de encontrar soluciones a ecuacio-

nes diferenciales no lineales utilizando la transformación diferencial, realizando extensiones
y/o generalizaciones para su aplicación a la resolución de problemas físicos, biológicos y de
ciencias en general.[4, 6]
El método de transformación diferencial consiste básicamente en un proceso de tipo

semi-analítico numérico el cual se aplica a todos los miembros de la ecuación en estudio,
iterando con los valores de la variable k para así obtener los coeficientes Y (k) necesarios y
luego obtener y(x) por medio de su transformada inversa.
La teoría y aplicación de las ecuaciones diferenciales con retardo es un tema muy impor-
tante dentro de los sistemas dinámicos, la física y las matemáticas aplicadas, además un caso
particular de las ecuaciones diferenciales funcionales. Se ocupa de modelos en los cuales las
funciones en cada instante dependen no sólo de un tiempo t, sino también de los valores
para instantes anteriores, lo cual significa que la función incógnita y sus derivadas están
evaluadas en argumentos distintos. Algunas aplicaciones importantes se presentan en mo-
delos de dinámica poblacional, transferencia de calor, controles biológicos entre otros. Quizá
uno de los primeros en realizar modificaciones a una ecuación diferencial ordinaria y trans-
formarla en una ecuación con retardo fue G. E. Hutchinson. Como ejemplo interesante, en la
ecuación logística se asume que la tasa de natalidad de cierto organismo depende instantá-
neamente de los cambios de tamaño poblacional, sin embargo existen múltiples fenómenos
que afectan el proceso generando una demora.
Todos los acontecimientos de la naturaleza requieren la existencia de la variable tem-

poral para darse o completar un ciclo. Algunos de ellos toman tan poco tiempo que no
generan inconvenientes perceptibles, mientras que otros como los que comprometen pro-
cesos biológicos son de alta afectación en la toma de datos y en la medición científica. Esta
sutíl observación detectada tiempo atrás fue el motivo para que en 1948 Hutchinson, quien
trabajaba con dinámica poblacional modificara su ecuación introduciendo un retraso en la
tasa de crecimiento.
Las ecuaciones diferenciales con retardo, llamadas DDEs por sus siglas en inglés (De-
lay Differential Equations), aparecen en muchos campos como la tecnología de control, redes
de comunicación y gestión de la población biológica, por lo que han atraído considerable
atención en los últimos años, lo cual ha hecho que surjan cierto métodos para encontrar
soluciones aproximadas. El método de transformación diferencial es uno de ellos, fue pro-
puesto por primera vez por [20] en la solución de problemas lineales y no lineales de valores
iniciales en el análisis de circuitos eléctricos. Otros autores han enriquecido la literatura re-
lativa a este tema en los últimos años, entre ellos [7], quien resolvió el problema lineal y no
lineal con DTM de ecuaciones diferenciales con retardo y comparó dichos resultados con el
método de descomposición de Adomian. En [19], por ejemplo estudió algunas aplicaciones
del DTM bidimensional a ecuaciones diferenciales parciales.
En general, la principal ventaja del método es el hecho de que proporciona una aproxi-
mación analítica, en muchos casos una solución exacta en una serie rápidamente convergen-
te. Los resultados actuales tanto en la parte teórica como computacional muestran que las
ecuaciones diferenciales con retardo son un modelo revelador de las complejas dinámicas
presentes en la naturaleza, como es el caso de la ecuación depredador-presa de Lotka-Volterra.
A nivel de software se han realizado diversas simulaciones, la cuales permiten corroborar

la efectividad del método. Tal es el caso de Shampine - Thompson quienes desarrollaron el
dde23 basados en MATLAB para solucionar ecuaciones diferenciales con retardo. [12]
1.3. Objetivos
1.3.1. General
Aproximar la solución de ecuaciones diferenciales con retardo en modelos biológicos por
el método de transformación diferencial.
1.3.2. Específicos
Estudiar las propiedades de la transformación diferencial aplicados a ecuaciones dife-
renciales con retardo.
Desarrollar el método de transformación diferencial para ecuaciones diferenciales con

retardo.
Realizar simulaciones que evidencien la rapidez de convergencia del método.
Comparar las aproximaciones encontradas con DTM y los obtenidos por otros méto-
dos.
1.4. Justificación
El método de transformación diferencial ha demostrado una gran efectividad a la hora
de aproximar soluciones a cierto tipo de ecuaciones diferenciales. De igual manera ocurre
cuando este método se extiende a ecuaciones diferenciales con retardo. La gran dificultad
que se presenta al resolver este tipo de ecuaciones por métodos conocidos en la literatura
hacen del DTM una gran herramienta por su gran efectividad y rapidez de convergencia. [3]
Esta es una de las razones por las cuales es importante continuar con su estudio ya que es
una contribución a un conocimiento incipiente1 que promete un gran desarrollo. Este trabajo
pretende ser una herramienta de consulta, ya que describe el método y sus aplicaciones de
manera ordenada, explicado paso a paso, haciendo de este un gran aporte pedagógico a las
matemáticas, en especial a las ecuaciones diferenciales no lineales y sus aplicaciones.
1
Sobre todo en la comunidad académica Colombiana
1.4.1. Planteamiento del problema

La dificultad y en muchos casos la imposibilidad para encontrar soluciones explícitas a
ecuaciones diferenciales con retardo hace necesario el desarrollo de métodos que permitan
facilitar los procesos, entre ellos, los métodos semi-analítico numéricos.
Este trabajo trata, como aporte pedagógico, el desarrollar desde la parte teórica y práctica
el método de transformación diferencial para ecuaciones diferenciales con retardo aplicados
especialmente a modelos biológicos. La comparación con resultados obtenidos utilizando
otros métodos, permitirán buscar y establecer el grado de efectividad del DTM, haciendo
mejoras que conduzcan a otras posibles aplicaciones en diversos campos de las ciencias.
1.5. Metodología
Dado el objeto de estudio del DTM para ecuaciones con retardo y su aplicación a mode-
los biológicos , los procedimientos a seguir para la consecución de datos y elaboración del
documento final correspondieron principalmente a la ejecución de barridos bibliográficos.
El trabajo inició con la búsqueda de documentos electrónicos (vía Internet), en los cua-
les se indagaron diferentes aspectos como el contenido histórico, implementación numérica,
demostraciones diversas, generalizaciones y aplicaciones.
El objetivo fue tener una visión global del tema y así comprenderlo ampliamente. De ahí
se escribió el anteproyecto, el cual fue el punto de partida para la elaboración y puesta en
marcha de los diferentes objetivos planteados.
La segunda etapa se inició con el primer borrador del proyecto mientras se continuó con
la indagación de nuevos artículos que mostraron sus aplicaciones a ecuaciones con retardo.
Se avanzó con la elaboración de ejercicios evidenciando paso a paso el método a trabajar.
Luego se realizaron algunas demostraciones importantes de método.
La tercera etapa consistió en la consulta a investigadores, principalmente aquellos rela-

cionados con la biomatemática, buscando ajustar sus propias apreciaciones y conocimientos
del tema a los borradores construidos en las dos primeras etapas de la investigación. En
etapa del proyecto se hizo la revisión inicial del primer borrador avalada por el director de
tesis. Con esta aprobación y con las sugerencias dadas, se realizó el segundo borrador del
proyecto.
En la cuarta etapa del proyecto se realizaron los procesos de simulación con MatLab y
Mathematica, al igual que la elaboración del último documento. Se hizo otro tipo de revi-
siones por parte del director y se escribió el documento final con las últimas sugerencias
hechas.
5
Capítulo 2
Generalidades de las ecuaciones

diferenciales con retardo
2.1. Ecuaciones Diferenciales con retardo (EDR)

Cuando estudiamos modelos con sistemas dinámicos, representamos dicho modelo por
medio de una ecuación diferencial o un conjunto de ellas, específicamente si hablamos de un
sistema biológico el cual evoluciona con el tiempo t, denominamos así el estado a la variable
dependiente x (dependiente del tiempo t) caracterizando entonces el estado de un sistema
cualquiera en equilibrio.
Presentamos a continuación algunas definiciones esenciales en el estudio de las EDR.
Definición 2.1.1. (Espacio de estados). Denominamos el espacio de estados, representado por C al

conjunto C([−τ, 0], R) con τ ≥ 0.
Definición 2.1.2. (Variable de estado). Sea la función continua x : R → Rn , t real y C un espacio

de estados. Definimos la variable de estado (o estado del sistema) xt : [−τ, 0] → Rn del sistema en el
tiempo t como
xt (φ) = x(t + φ)
para todo φ ∈ [−τ, 0].
Definición 2.1.3. (EDR). Definimos ecuación diferencial con retardo a una ecuación de la forma
x0 (t) = f (t, x(t), x(t − r))
donde la función f es continua con parámetro r > 0 denominado el retardo de la ecuación.
Definición 2.1.4. Definimos problema de valor inicial de una EDR como

(
x0 (t) = f (t, x(t), x(t − r)), t ≥ γ
(2.1)
xγ = ϕ
donde γ ∈ R representa el tiempo inicial, ϕ ∈ C representa el estado del sistema respecto al tiempo γ
y la función f : R × C → Rn continua.
Definición 2.1.5. Definimos como una solución del problema de valor inicial (2.1) a la función
x : R → Rn la cual satisface precisamente (2.1).
Capítulo 2. Generalidades de las ecuaciones diferenciales con retardo 6
Es de anotar que en las ecuaciones diferenciales con retardo, el comportamiento de sus

soluciones (en un instante t) depende no sólo de la posición x(t) sino también del valor en
un instante anterior x(t − r).
Otro de los casos importantes de las EDR, específicamente el caso general, es el de las
ecuaciones funcionales, en las que el comportamiento de las soluciones en t depende de la
solución en el intervalo [t − r, t]. Este tipo de ecuaciones tienen la forma
x0 (t) = f (t, xt )
en la cual xt (s) = x(t + s) con xt : [−r, 0] → R.
Ahora bien, un problema de valor inicial para el caso general está dado por
(
x0 = f (t, xt )
(2.2)
xγ = ϕ
con xt (φ) := x(t + φ), φ ∈ [−τ, 0].
2.2. Modelo matemático con retardo: Caso especial

Supongamos la ecuación diferencial con retardo
x0 (t) = −αx(t − r) (2.3)
con α > 0. Podemos ver que esta ecuación es una EDR lineal cuyo retardo es r. Si tenemos
en cuenta la ecuación diferencial ordinaria (lineal)
dx
= −αx(t),
dt
podemos ver que el comportamiento de las soluciones varía respecto a los parámetros α y r,
ya que la ecuación característica asociada x + α = 0 es no lineal como sucede en el caso de
la EDO, sino que es una ecuación de tipo trascendente de la forma x + αe−rx = 0.
Por ejemplo, si tenemos la ecuación diferencial con retardo
x0 (t) = −0.2x(t − 1),
vemos entonces que si (α)r < e−1 , las soluciones convergerán a cero, tal y como pasa en el
caso de la EDO, concluyendo así que los retardos pequeños no influyen en la dinámica (ver
Figura 2.1).
Ahora bien, supongamos la ecuación con retardo
x0 (t) = −1.5x(t − 1),

1.8
1.6
1.4
1.2
0.8
0.6
0.4
0.2
0
−5 0 5 10 15 20
F IGURA 2.1: x0 (t) = −0.2x(t − 1).
por lo tanto, si los valores de (α)r varían entre e−1 y π2 , es decir, e−1 < (α)r < π2 , las soluciones
convergen de manera oscilatoria a cero, fenómeno característico de las EDR (ver figura 2.2).
1.5
0.5
−0.5
−1
−1.5
−2
−5 0 5 10 15 20 25 30 35 40
F IGURA 2.2: x0 (t) = −1.5x(t − 1).
Finalmente, supongamos la ecuación con retardo

0 T
x (t) = − x(t − 1)
2
Así, si (α)r ≥ π2 , tendremos que el comportamiento de las soluciones es periódico.

2.3. Sobre la existencia y unicidad de soluciones

Repasamos a continuación los elementos necesarios para garantizar la existencia y uni-
cidad de solución de las ecuaciones diferenciales con retardo.
Definición 2.3.1. Sean el conjunto acotado A = R × C, las constantes reales a, b y M, K > 0.
Decimos que la función f es Lipschitziana con constante K en A si
kf (t, φ) − f (t − ϕ)k ≤ Kkφ − ϕk, (2.4)
para todo t ∈ [a, b] y kφk, kϕk ≤ M .

Lema 2.3.1. Sea f : A → Rn continua que satisface (2.4). Entonces para todo t ∈ [a, b], existen
constantes M, L > 0 tal que
|f (t, ϕ)| ≤ L
con kϕk ≤ M .
Lema 2.3.2. (Gronwall). Sean las funciones ϕ(t), α(t) y β(t) continuas y no negativas en un in-
tervalo I. Sean C ≥ 0 y t0 , t ∈ I. Si
Z t

ϕ(t) ≤ C + (α(s)ϕ(s) + β(s))ds
(2.5)
t0
entonces Z t R
t α(s)ds
ϕ(t) ≤ C + β(s)ds e
t0
(2.6)
t0
Teorema 2.3.1. (Existencia y unicidad de solución). Sea f una función continua la cual satisface
(2.4), entonces existe una solución única del problema (2.1) en el intervalo [γ − τ, γ + B], con B =
mı́n{τ, ML
}, γ ∈ R y ϕ ∈ C tales que kϕk ≤ M . Además, si K es la constante Lipschitziana, entonces

máx x(β, ϕ) − x(β, ψ) ≤ kϕ − ψkeKB , (2.7)
γ−τ ≤β≤γ+B
donde kϕk, kψk ≤ M .

Nota 2.3.1. Entre las ecuaciones diferenciales ordinarias y las ecuaciones diferenciales con retardo
podemos encontrar una diferencia importante respecto a la naturaleza del valor inicial, es decir, para
las EDO, la condición inicial no es más que un escalar (o vector), y en el caso de las EDR la condición
inicial es un conjunto de valores de x sobre un intervalo dado, es decir, de naturaleza funcional de la
forma x0 : [−τ, 0] → R.
Nota 2.3.2. En el siguiente ejemplo podemos ver que si consideramos otra condición inicial, entonces
las soluciones del nuevo PVI podrían intersectarse con las soluciones retardadas, es decir, del problema
anterior.
Ejemplo 1.
Sean los problemas de valor inicial
(
x0 (t) = − π2 x(t − 1)
(2.8)
x0 (t) = 0, ∀ t ∈ [−1, 0]
y (
y 0 (t) = − π2 y(t − 1)
(2.9)
y0 (t) = sin π2 t , ∀ t ∈ [−1, 0]

Podemos ver así que en el problema (2.8) la solución es x ≡ 0 (función nula); mientras que en el
problema (2.9) tenemos que para todo t ∈ [0, 1]
y 0 (t) = − π2 sin π2 (t − 1) ,

donde integrando obtenemos

Z t
− π2 π

y(t) = sin 2
(u − 1) du + y0 (0)
0
Z t
π π

= 2
cos 2
u du
0
= sin π2 t ,

∀ t ∈ [0, 1]
Este resultado puede generalizarse diciendo que para todo m ≥ 1,

Z t
π π
π

y(t) = y(m − 1) + cos 2
u du = sin 2
t
2 m−1
donde vemos que ambas soluciones se intersectan para todo t = kπ con k ∈ N y t ∈ [m − 1, m] (ver
figura 2.3).
1.5
0.5
−0.5
−1
−1.5
−2 0 2 4 6 8 10
F IGURA 2.3: Intersección con las soluciones retardadas

10
Capítulo 3
Método de Transformación Diferencial

DTM
3.1. Antecedentes
El método de transformación diferencial es prácticamente desconocido en los países
occidentales y se ha utilizado en décadas recientes en países como Irán, Iraq, Turquía, Ara-
bia Saudita y demás repúblicas islámicas para la resolución de ecuaciones diferenciales no
lineales (ordinarias o parciales), como método numérico y cuasi analítico en problemas de
ciencias e ingeniería. Ha aportado en los últimos años velocidad en la convergencia, aho-
rrando tiempo de cómputo en la consecución de la solución en series de potencias para
problemas que por los métodos estrictamente analíticos son en gran medida difíciles o im-
posibles de resolver.
Fue utilizado inicialmente por Zhou en la solución de problemas en circuitos eléctricos

y, dada su rápida convergencia, se hizo extensivo a toda una gama de problemas de física e
ingeniería, los cuales requerían de un método que permitiera minimizar tiempos de cálculo
y en algunos casos, encontrar la solución exacta al problema mismo a partir de las series de
potencias que se arrojaban como soluciones.
El método de transformación diferencial nace como una nueva alternativa para la re-
solución de las ecuaciones diferenciales no lineales, las cuales son la representación de la
mayoría de los fenómenos físicos que requieren de interpretación y solución matemática.
La solución en series infinitas y a través de los distintos métodos numéricos para ciertas
ecuaciones diferenciales brinda una opción importante en eventos en los cuales hallar su
solución por métodos analíticos típicos se hace difícil y en algunos casos, imposible.
La transformación diferencial, al igual que todo un grueso de transformadas (Fourier, La-

place, Mellin, etc.), es un recurso importante para la solución de problemas con valor inicial
y con valor en la frontera para ecuaciones diferenciales no lineales ordinarias y parciales,
aportando soluciones en series infinitas, muchas de las cuales pueden ser representación
de funciones explicitas conocidas [13]. Su aplicabilidad se evidencia en el uso del método
de transformación diferencial para encontrar soluciones, aplicando la transformación y su
transformación inversa, a través de iteraciones consecutivas que conducen a visualizar el
comportamiento de dichas soluciones.
Capítulo 3. Método de Transformación Diferencial DTM 11
3.2. Generalidades del DTM

3.2.1. Transformación Diferencial
Sea y(x) una función continuamente diferenciable en un intervalo de la forma (x0 −r, x0 +
r) con r > 0. Definimos la transformación diferencial de y(x) escrita Y (k) como
1 dk y(x)

Y (k) = (3.1)
k! dxk x=x0
El valor de la función y(x) es obtenido a partir de la transformación diferencial inversa de

Y (k)
X∞
y(x) = (x − x0 )k Y (k) (3.2)
k=0
Podemos observar que al sustituir Y (k) en la transformada diferencial inversa obtenemos

∞
dk y(x)

X 1 k
y(x) = (x − x0 ) , (3.3)
k=0
k! dxk x=x0
es decir, la solución corresponde a la serie de Taylor en x0 de la función y(x). Ahora bien, si

asumimos que x0 = 0, tenemos que
1 dk y(x)

Y (k) = (3.4)
k! dxk x=0
y
∞
dk y(x)

X 1 k
y(x) = (x) (3.5)
k=0
k! dxk x=0
La transformación diferencial cumple con propiedades de linealidad al igual que otras

transformaciones conocidas en la literatura. En los siguientes teoremas se enuncian y mues-
tran algunas de estas propiedades, al igual que otras correspondientes a operaciones de
producto, derivación y transformaciones propias de algunas funciones especiales.
Sean las funciones f (x), g(x) y h(x) (continuamente diferenciables) con transformadas
diferenciales F (x), G(x) y H(x) respectivamente. Enunciamos los siguientes teoremas y pre-
sentamos sus respectivas pruebas. [2, 5, 14]
Teorema 3.2.1. Si f (x) = g(x) ± h(x), entonces F (k) = G(k) ± H(k).

Demostración.
1 dk

F (k) = (g(x) ± h(x))
k! dxk x=x0
k k

1 d d
= g(x) ± k h(x)
k! dxk dx x=x
k k0
1 d g(x) 1 d
= ± h(x)
k! dxk x=x0 k! dxk x=x0
= G(k) ± H(k)
Teorema 3.2.2. Si f (x) = αg(x) con α ∈ R, entonces F (k) = αG(k).

Demostración.
1 dk

F (k) = (αg(x))
k! dxk x=x0
k
1 d
=α g(x)
k! dxk x=x0
= αG(k)
dy
Teorema 3.2.3. Si f (x) = , entonces F (k) = (k + 1)Y (k + 1).
dx
Demostración.
1 dk dy

F (k) =
k! dxk dx x=x0
1 dk+1 y

=
k! dxk+1 x=x0
k+1
1 d y
= (k + 1)
(k + 1)! dxk+1 x=x0
= (k + 1)Y (k + 1)
dn y
Teorema 3.2.4. Si f (x) = n , entonces F (k) = (k + 1)(k + 2) · · · (k + n)Y (k + n)
dx
Demostración.
1 dk dn y

F (k) =
k! dxk dxn x=x0
1 dk+n y

=
k! dxk+n x=x0
k+n
1 d y
= (k + 1)(k + 2) · · · (k + n)
(k + n)! dxk+n x=x0
= (k + 1)(k + 2) · · · (k + n)Y (k + n)
Pk
Teorema 3.2.5. Si f (x) = g(x) h(x), entonces F (k) = r=0 G(k)H(k − r).
Demostración.
1 dk

F (k) = g(x) h(x)
k! dxk x=x0
" k #
1 X k dk−L g(x) dL h(x)
=
k! L=0 L dxk−L dxL
x=x0
k
dk−L g(x) dL h(x)
X 1
k!
=
L=0
k! (k − L)!L! dxk−L dxL x=x0
k
" #
1 dL h(x) dk−L g(x)

X 1
=
L=0
L! dxL x=x0 (k − L)! dxk−L x=x0
k
X
= H(L)G(k − L)
L=0
Teorema 3.2.6. Si f (x) = xn , entonces F (k) = δ(k − n), donde δ corresponde a la función delta de
Dirac.
Demostración.
1 dk (xn )

F (k) =
k! dxk x=x0
1
n(n − 1) · · · (n − (k − 1))xn−k x=x0

=
k!
1
Ahora, si por ejemplo x0 = 1, tenemos que F (k) = [n(n − 1) · · · (n − (k − 1))]. Si k = n,
k!
6 n, entonces F (k) = 0. Por consiguiente, F (k) =
entonces F (k) = 1. Finalmente, si k =
δ(k − n).
αk
Teorema 3.2.7. Si f (x) = eαx , entonces F (k) = .
k!
Demostración.
1 dk eαx

F (k) =
k! dxk x=x0 =0
1
αk eαx x=x0 =0

=
k!
αk
=
k!
con α ∈ R.
(−1)k+1
Teorema 3.2.8. Si f (x) = ln(x), entonces F (k) = .
k
Demostración.
1 dk ln(x)

F (k) =
k! dxk x=1
1 (−1)k+1 (k − 1)!

=
k! xk x=1
1 (−1)k+1
= (k − 1)!
k(k − 1)! (xk ) x=1
(−1)k+1 (−1)k+1
= =
k(x)k k
Teorema 3.2.9. Si f (x) = 1, entonces F (k) = δ(k).

Demostración.
1 dk (1)

F (k) =
k! dxk x=x0
(
0, k > 0
=
1, k = 0
por lo tanto, F (k) = δ(k).
Teorema 3.2.10. (Cárdenas). Si f (x) = xm y(x), entonces

(
0, k<m
F (k) =
Y (k − m), k ≥ m
para todo m ∈ N.
Demostración. La prueba la realizaremos utilizando inducción. Asumamos inicialmente que
es verdad para m, entonces probaremos que es verdad para m + 1. Usando DTM al producto
tenemos:
f (x) = xm+1 y(x) = x (xm y(x))
y así obtenemos
k
X
F (k) = Y1 (k1 )Y2 (k − k1 )
k1 =0
donde y1 (x) = x y y2 (x) = xm y(x). Por lo tanto

k
X
F (k) = δ(k1 − 1)Y2 (k − k1 )
k1 =0
Cuando k = 0, la expresión se reduce a
δ(−1)Y2 (0) = 0
Ahora, si k > 0, entonces todos los términos de la suma son iguales a cero excepto k1 = 1,
obteniendo así
(
0, k−1<m
F [xm+1 y(x)](x) = F [xm y(x)](k − 1) =
Y (k − 1 − m), k − 1 ≥ m
16
Capítulo 4
DTM para ecuaciones diferenciales con

retardo
4.1. Introducción
Consideremos inicialmente la ecuación diferencial con retardo de orden n
y (n) (x) = f (x, y(x), y(α1 (x)), . . . , y(αm (x))) (4.1)
con x ∈ [0, b] sujeta a las condiciones
y(x), y 0 (x), . . . , y (n−1) (x) = (h1 (x), h2 (x), . . . , hn (x))

donde x ∈ [a, 0]. Acá, a es el mínimo de los valores αi (x) ≤ x (funciones retardo) para todo
x ∈ [0, b] y todo i ∈ [1, m]. Asumamos que las funciones retardo αi (x) y gi (x) son suficiente-
mente suaves.
En la siguiente sección hacemos una extensión de lo tratado en el capítulo anterior, es

decir, la extensión del método de transformación diferencial aplicado a ecuaciones diferen-
ciales con retardo.
4.2. Extensión del DTM a ecuaciones con retardo

Usando el concepto central de la transformación diferencial de una función f (x), for-
mulamos los siguientes teoremas que son centrales para el desarrollo de esta investigación.
[14]
Teorema 4.2.1. Asumamos que F (k) y Y (k) son las transformadas diferenciales de las funciones
f (x) y y(x) respectivamente y α, β ∈ (0, 1). Entonces si f (x) = y(αx) tenemos que F (k) = αk Y (k).
Demostración. Por definición de DTM, tenemos
dn f (x) dn y(αx)
=
dxn dxn
Capítulo 4. DTM para ecuaciones diferenciales con retardo 17
Por tanto, aplicando los teoremas del DTM llegamos a
dk
k!F (k) = αk y(x̃), con αx = x̃
dx̃k
= αk k!Y (k)
= αk Y (k), k = 0, 1, 2, . . .
x 1
Corolario 4.2.1. Si f (x) = y con α ≥ 1, entonces F (k) = Y (k).
α αk
Demostración.
dk dk x
f (x) = y
dxk dxk α
1 dk
k!F (k) = k k y(x̃), con αx = x̃
α dx̃
1
= k k!Y (k)
α
1
F (k) = k Y (k), k = 0, 1, 2, . . .
α
Teorema 4.2.2. Si f (x) = y1 (αx) y2 (βx), con α, β ≥ 1, entonces

k
X
F (k) = αl β k−l Y1 (l)Y2 (k − l)
l=0
Demostración. De antemano, recordemos la fórmula de Leibnitz para la derivada n-ésima de

un producto, es decir,
n
(n)
X n (n−k) (k)
(uv) = u v (4.2)
k=0
k
Ahora bien, sean Y1 (k) y Y2 (k) las transformadas diferenciales de las funciones y1 y y2 res-
pectivamente, entonces
dk dk
f (x) = [y1 (αx) y2 (βx)]
dxk dxk
k l
X k d dk−l
k!F (k) = l
(y 1 (αx)) k−l
(y2 (βx))
l=0
l dx dx
k
X k l dl dk−l
= α l y1 (x̃)β k−l k−l y2 (x̂), con x̃ = αx, x̂ = βx
l=0
l dx̃ dx̂
k
X k!
αl l!Y1 (l) β k−l (k − l)!Y2 (k − l)

=
l=0
(k − l)!l!
Xk
F (k) = αl β k−l Y1 (l)Y2 (k − l)
l=0
x
x
Corolario 4.2.2. Si f (x) = y1 y2 con α, β ≥ 1, entonces
α β
k
X 1 1
F (k) = Y1 (l)Y2 (k − l)
l=0
αl β k−l
Demostración.
dk dk

x x
f (x) = y 1 y 2
dxk dxk α β
k
X k 1 dl 1 dk−l

= l dx̃l
y 1 (x̃) y (x̂), con x̃ = αx , x̂ =
k−l dx̂k−l 2
x
β
l=0
l α β
k
X k! 1 1
= l
l!Y1 (l) k−l (k − l)!Y2 (k − l)
l=0
(k − l)!l! α β
k
X 1 1
k!F (k) = k! l
Y1 (l) k−l Y2 (k − l)
l=0
α β
k
X 1 1
F (k) = Y1 (l)Y2 (k − l)
l=0
αl β k−l
dn
Teorema 4.2.3. Si f (x) = y(αx), entonces
dxn
(k + n)! k+n
F (k) = α Y (k + n)
k!
Demostración.
dk dk
y (n) (αx)

k
f (x) = k
dx dx
dk+n
= αk+n k+n y(x̃), con αx = x̃
dx̃
k!F (k) = αk+n (k + n)!Y (k + n)
(k + n)!
F (k) = αk+n Y (k + n)
k!
(k + n)! k+n
= α Y (k + n)
k!
dn
Corolario 4.2.3. Si f (x) = y(x − α) con α > 0, entonces
dxn
M
(k + n)! X i−k−n i
F (k) = (−α) Y (i)
k! i=k+n
k+n
para M → ∞.
La demostración de este corolario utiliza un argumento similar al del teorema inmedia-

tamente anterior.
Z mx
1
Teorema 4.2.4. Si f (x) = y(αs) ds, entonces F (k) = mk αk−1 Y (k − 1).
0 k
Demostración.
Z mx
dk dk

f (x) = k y(αs) ds
dxk dx 0
Z mx
dk−1 d

= k−1 y(αs) ds
dx dx 0
dk−1
= k−1 [my(αmx)]
dx
dk−1
k!F (k) = m(αm)k−1 k−1 y(x̃), con x̃ = αmx
dx̃
k−1 k−1
= m(α) m (k − 1)!Y (k − 1)
mk (α)k−1 (k − 1)!Y (k − 1)
=
k(k − 1)!
1
= mk αk−1 Y (k − 1)
k
Z mx
Corolario 4.2.4. Si f (x) = y1 (αs)y2 (βs) ds, entonces
0
k−1
1 X k l k−l−1
F (k) = m αβ Y1 (l)Y2 (k − l − 1)
k l=0
Demostración.
Z mx
dk dk
f (x) = k y1 (αs)y2 (βs) ds
dxk dx 0
Z mx
dk−1 d

= k−1 y1 (αs)y2 (βs) ds
dx dx 0
dk−1
= k−1 [my1 (αmx)y2 (βmx)]
dx
dk−1
= m k−1 [y1 (αmx)y2 (βmx)]
dx
k−1
dl dk−l−1

X k−1
=m (αm)l l y1 (x̃)(βm)k−l−1 k−l−1 y2 (x̂)
l=0
l dx̃ dx̂
k−1
X k−1
=m (αm)l l!Y1 (l)(βm)k−l−1 (k − l − 1)!Y2 (k − l − 1)
l=0
l
k−1
X (k − 1)!
k!F (k) = m (αm)l l!Y1 (l)(βm)k−l−1 (k − l − 1)!Y2 (k − l − 1)
l=0
l!(k − 1 − l)!
k−1
1 X
= mk αl β k−l−1 Y1 (l)Y2 (k − l − 1)
k l=0
Teorema 4.2.5. Si f (x) = g(x − α) con α > 0, entonces

N
i−k i
X
F (k) = (−α) G(i),
i=k
k
con N → ∞.
Demostración. Utilizando los teoremas para la transformada diferencial de un producto de
funciones tenemos
∞
X
f (x) = G(k)(x − x0 − α)k
k=0
∞ k
X X
k−i k
= G(k) (−1) (x − x0 )i αk−i
k=0 i=0
i
∞ Xk
X
k−i k
= (−1) (x − x0 )i αk−i G(k)
k=0 i=0
i
∞ X
∞
X
k−i k k−i
i
= (−1) (x − x0 ) α G(k)
i=0 k=i
i
∞ ∞
k−i k
X X
i
= (x − x0 ) (−1) αk−i G(k)
i=0 k=i
i
∞ ∞
i−k i
X X
k
= (x − x0 ) (−α) G(i)
k=0 i=k
k
Teorema 4.2.6. (Cárdenas-Trujillo). Si f (x) = g 2 (x − α) con α > 0, entonces

k
" N N #
X X i X h1
F (k) = (−α)i−l G(i) · (−α)h1 −k+l G(h1 )
l=0 i=l
l h =k−l
k − l
1
Para la prueba de este teorema debemos tener presente que podemos escribir la función
f (x) como f (x) = g(x − α) · g(x − α); posteriormente utilizar la transformada diferencial de
un producto, es decir,
X k
F (k) = G(l)H(k − l)
l=0
donde F es la transformada de la función f (x) y G, H son las transformadas de la función

g(x − α).
Demostración. Empezamos escribiendo como se mencionó antes que
f (x) = g 2 (x − α) = g(x − α) · g(x − α)
Por lo tanto, aplicando la fórmula de Leibnitz (4.2) tenemos

k l
dk dk X k d dk−l
f (x) = [g(x − α) · g(x − α)] = g(x − α) g(x − α)
dxk dxk l=0
l dx l dx k−l
Ahora bien, tomando x − α = x∗ obtenemos

k l k
X k d dk−l X k dl ∗ d
k−l
l
g(x − α) k−l
g(x − α) = ∗ l
g(x ) ∗ k−l
g(x∗ )
l=0
l dx dx l=0
l dx dx
donde utilizando la definición de transformación diferencial llegamos a

k k
X k dl ∗ d
k−l
∗
X k
∗l
g(x ) ∗ k−l
g(x ) = l!G(l)(k − l)!G(k − l)
l=0
l dx dx l=0
l
k
X
= k!G(l)G(k − l) (4.3)
l=0
Aplicando el teorema 4.2.5 , las expresiones G(l) y G(k − l) en (4.3) se transforman en

N
i−l i
X
G(l) = (−α) G(i) (4.4)
i=l
l
y
N
X
h1 −k+l h1
G(k − l) = (−α) G(h1 ) (4.5)
h1 =k−l
k−l
Finalmente, reemplazando (4.4) y (4.5) en (4.3) se tiene el resultado final
dk dk
f (x) = [g(x − α) · g(x − α)]
dxk dxk
Xk X k
= k!G(l)G(k − l) = k! G(l)G(k − l)
l=0 l=0
k
" N N #
X X i X h1
F (k) = (−α)i−l G(i) · (−α)h1 −k+l G(h1 )
l=0 i=l
l h =k−l
k − l
1
Corolario 4.2.5. Si f (x) = y1 (αx)y2 (x − β), entonces

k N
X X
l h1 −k+l h1
F (k) = α (−β) Y1 (l)Y2 (h1 )
l=0 h1 =k−l
k−l
Corolario 4.2.6. Si f (x) = y1 (x − α)y2 (x − β), entonces

k
" N N #
X X i X h1
F (k) = (−α)i−l G(i) · (−β)h1 −k+l G(h1 )
l=0 i=l
l h =k−l
k−l
1
Para la prueba de estos dos últimos corolarios, utilizamos el teorema 4.2.6 donde el es-
quema de la demostración es similar a dicho teorema.
4.3. Aplicaciones a ecuaciones con retardo

Ejemplo 1.
Supongamos queremos resolver el siguiente problema de valor inicial
(
y 00 (x) − 43 y(x) − y( x2 ) + x2 − 2 = 0
y(0) = 0, y 0 (0) = 0
para todo x ∈ [0, 1]. Es fácil ver que la función y(x) = x2 es solución del problema en cues-
tión. Ahora, aplicando los teoremas del DTM para ecuaciones con retardo tenemos: [15]
(k + 2)! 1
Y (k + 2) = 34 Y (k) + k Y (k) − δ(k − 2) + 2δ(k − 0)
k! 2

k! 3 1
Y (k + 2) = + Y (k) − δ(k − 2) + 2δ(k) (4.6)
(k + 2)! 4 2k
De las condiciones iniciales tenemos que para y(0) = 0, Y (0) = 0 y para y 0 (0) = 0, Y (1) = 0.
Utilizando la ecuación de recurrencia (4.6) tenemos que:
Si k = 0, entonces

0! 3 1
Y (2) = + Y (0) − δ(−2) + 2δ(0)
2! 4 1
Y(2)=1
Si k = 1, entonces

1! 3 1
Y (3) = + Y (1) − δ(−1) + 2δ(1)
3! 4 2
Y(3)=0
Si k = 2, entonces

2! 3 1 1
Y (4) = + Y (2) − δ(0) + 2δ(2) = (1(1) − 1)
4! 4 4 12
Y(4)=0
Si k = 3, entonces

3! 3 1
Y (5) = + Y (3) − δ(1) + 2δ(3)
5! 4 8
Y(5)=0
Si k = 4, entonces

4! 3 1
Y (6) = + Y (4) − δ(2) + 2δ(4)
6! 4 16
Y(6)=0
y así sucesivamente. Por lo tanto, reemplazando en la ecuación (3.2) tenemos:
y(x) = Y (0)x0 + Y (1)x1 + Y (2)x2 + Y (3)x3 + . . .

= 0x0 + 0x1 + 1x2 + 0x3 + . . .
= x2
la cual corresponde efectivamente a la solución exacta.

Ejemplo 2.
Supongamos queremos resolver el siguiente problema de valor inicial [16]
(
y 000 (x) − 2y 2 x2 + 1 = 0

y(0) = 0, y 0 (0) = 1, y 00 (0) = 0
para todo x ∈ [0, 1]. Podemos verificar sin problema que la función y(x) = sin x es solución
exacta de la ecuación diferencial. Además, la ecuación puede ser escrita como
y 000 (x) − 2y x2 y x2 + 1 = 0,

esto, con el fin de aplicar los teoremas del DTM. En efecto, tenemos que
k
(k + 3)! X1 1
Y (k + 3) = −δ(k − 0) + 2 Y (l)Y (k − l)
k! l=0
2l 2k−l
k
X 1 1
= −δ(k) + 2 Y (l)Y (k − l)
l=0
2l 2k−l
k
X 1
= −δ(k) + 2 k
Y (l)Y (k − l)
l=0
2
Así llegamos a la ecuación de recurrencia para la ecuación con retardo

" k
#
k! X 1
Y (k + 3) = −δ(k) + 2 Y (l)Y (k − l) (4.7)
(k + 3)! l=0
2k
De las condiciones iniciales tenemos que para y(0) = 0, Y (0) = 0; para y 0 (0) = 1, Y (1) = 1 y
para y 00 (0) = 0, Y (2) = 0. Utilizando la ecuación (4.7) tenemos que:
Si k = 0, entonces

0! 1
Y (3) = −δ(0) + 2 0 Y (0)Y (0)
3! 2
1
Y (3) = −
3!
Si k = 1, entonces

1! 1 1
Y (4) = −δ(1 − 0) + 2 1 Y (0)Y (1) + 1 Y (1)Y (0)
4! 2 2
Y (4) = 0
Si k = 2, entonces

2! 1 1 1 2! 1
Y (5) = −δ(2) + 2 2 Y (0)Y (2) + 2 Y (1)Y (1) + 2 Y (2)Y (0) = 2 2
5! 2 2 2 5! 2
1
Y (5) =
5!
Si k = 3, entonces

3! 1 1 1 1
Y (6) = −δ(3) + 2 3 Y (0)Y (3) + 3 Y (1)Y (2) + 3 Y (2)Y (1) + 3 Y (3)Y (0)
6! 2 2 2 2
3!
= (0 + 2(0)) = Y(6)=0
6!
Si k = 4, entonces

4! 1 1 1 1
Y (7) = −δ(4) + 2 4 Y (0)Y (4) + 4 Y (1)Y (3) + 4 Y (2)Y (2) + 4 Y (3)Y (1)
7! 2 2 2 2

1 4! 1 1 −4!4
+ 4 Y (4)Y (0) = 2 − 4 − 4 =
2 7! 2 3! 2 3! 7!24 3!
1
Y (7) = −
7!
y así sucesivamente. Por lo tanto, utilizando la ecuación (3.2) llegamos a
y(x) = Y (0)x0 + Y (1)x1 + Y (2)x2 + Y (3)x3 + . . .

1 1
= 0x0 + 1x1 + 0x2 − x3 + 0x4 + x5 + . . .
3! 5!
1 3 1 5 1 7
= x − x + x − x + ...
3! 5! 7!
la cual corresponde con la solución exacta mencionada inicialmente, es decir, y(x) = sin x.
Este resultado podemos escribirlo en forma compacta como

0, k = 2m
Y (k) = m
(−1)
 , k = 2m + 1
(2m + 1)!
para todo m ∈ N.
Ahora bien, es importante resaltar que el problema inmediatamente anterior puede escri-
birse como un sistema de ecuaciones diferenciales realizando sustituciones adecuadas. Por
tanto, mostramos a continuación la aplicación del método de transformación diferencial a
sistemas de ecuaciones diferenciales con retardo. [17]
Ejemplo 3.
Sea el problema de valor inicial (ver ejemplo 2)
(
y 000 (x) − 2y 2 x2 + 1 = 0

y(0) = 0, y 0 (0) = 1, y 00 (0) = 0
para todo x ∈ [0, 1]. Por medio de las sustituciones y1 = y, y2 = y 0 y y3 = y 00 podemos escribir
la ecuación diferencial con retardo como el sistema

0
y1 = y2

y20 = y3

 0
y3 = 2y12 ( x2 ) − 1
donde y10 = y 0 = y2 , y20 = y 00 = y3 y y30 = y 000 = −1 + 2y12 ( x2 ). Acá, las condiciones iniciales son
y1 (0) = y3 (0) = 0 y y2 (0) = 1. Es fácil ver que como y(x) = sin x es la solución exacta al PVI,
entonces y1 (x) = sin(x), y2 (x) = cos(x) y y3 (x) = − sin(x).
Ahora podemos aplicar el DTM al sistema, obteniendo respectivamente
(k + 1)!
Y1 (k + 1) = Y2 (k) (4.8)
k!
(k + 1)!
Y2 (k + 1) = Y3 (k) (4.9)
k!
k
(k + 1)! X 1 1
Y3 (k + 1) = −δ(k − 0) + 2 Y1 (l)Y1 (k − l) (4.10)
k! l=0
2l 2k−l
De nuevo, las condiciones iniciales nos arrojan que Y1 (0) = Y3 (0) = 0 y Y2 (0) = 1. Entonces,
iterando las ecuaciones (4.8), (4.9) y (4.10) tenemos
Si k = 0,
0!
Y1 (1) = Y2 (0) = 1
1!
0!
Y2 (1) = Y3 (0) = 0
1!
0! 1
Y3 (1) = −1 + 2 0 Y1 (0)Y1 (0) = -1
1! 2
Si k = 1,
1!
Y1 (2) = Y2 (1) = 0
2!
1! 1
Y2 (2) = Y3 (1) = −
2! 2!

1! 1 1
Y3 (2) = 0 + 2 1 Y1 (0)Y1 (1) + 1 Y1 (1)Y1 (0) = 0
2! 2 2
Si k = 2,
2! −1
Y1 (3) = Y2 (2) =
3! 3!
2!
Y2 (3) = Y3 (2) = 0
3!
2! 1 1 1
Y3 (3) = 0 + 2 2 Y1 (0)Y1 (2) + 2 Y1 (1)Y1 (1) + 2 Y1 (2)Y1 (0)
3! 2 2 2

2! 1 2! 1 1
= 0+2 2 = =
3! 2 3! 2 3!
y1 (x) = Y1 (0)x0 + Y1 (1)x1 + Y1 (2)x2 + Y1 (3)x3 + . . .

1 1
= 0x0 + 1x1 + 0x2 − x3 + 0x4 + x5 + . . .
3! 5!
1 3 1 5 1 7
= x − x + x − x + ...
3! 5! 7!
≈ sin(x)
y2 (x) = Y2 (0)x0 + Y2 (1)x1 + Y2 (2)x2 + Y2 (3)x3 + . . .

1
= 1x0 + 0x1 − x2 + 0x3 + . . .
2!
1 2 1 4
= 1 − x + x − ...
2! 4!
≈ cos(x)
y3 (x) = Y3 (0)x0 + Y3 (1)x1 + Y3 (2)x2 + Y3 (3)x3 + . . .

1
= 0x0 + (−1)x1 + 0x2 + x3 + 0x4 + . . .
3!
1 3 1 5 1 7
= −x + x − x + x + . . .
3! 5! 7!
≈ − sin(x)
que son en efecto las soluciones mencionadas inicialmente.
Ejemplo 4.
Consideremos a continuación el siguiente problema integro-diferencial con retardo pro-
porcional [18]  Z x

 0 x
1 3
y(t)y 2t dt = 0

y (x) − y 2 − 2
0
(4.11)
y(0) = y 0 (0) = 1

Utilizando el corolario 4.2.4 tenemos

" k−1 #
k k−l−1
(k + 1)! 1 1 1X 1 1
Y (k + 1) − k Y (k) − (1)l Y (l)Y (k − l − 1) = 0
k! 2 2 k l=0 3 2
o lo que es lo mismo
" k−1 k k−l−1
#
1 1 1 X 1 1
Y (k + 1) = k
Y (k) + Y (l)Y (k − l − 1) (4.12)
k+1 2 2(k) l=0 3 2
De las condiciones iniciales y(0) = y 0 (0) = 1 tenemos Y(0)=1 y Y(1)=1 . Ahora bien,
iterando la fórmula de recurrencia (4.12) tenemos:
Si k = 1, entonces
" 0 ! #
1 1 1 1 1
Y (2) = Y (1) + Y (0)Y (0)
2 2 2 3 2

1 1 1 1
= Y (1) + Y (0)Y (0) =
4 4 3 3
Si k = 2, entonces
" 2 1 2 0 !#
1 1 1 1 1 1 1
Y (3) = Y (2) + Y (0)Y (1) + Y (1)Y (0)
3 4 4 3 2 3 2

1 1 1 1 1 1
= + + =
3 12 4 18 9 24
Si k = 3, entonces
" 3 2 3 1 3 0 !#
1 1 1 1 1 1 1 1 1
Y (4) = Y (3) + Y (0)Y (2) + Y (1)Y (1) + Y (2)Y (0)
4 8 6 3 2 3 2 3 2

1 1 1 1 1 1 169
= + + + =
4 192 6 324 54 81 62208
Si k = 4, entonces
" 4 3 4 2 4 1
1 1 1 1 1 1 1 1 1
Y (5) = 4
Y (4) + Y (0)Y (3) + Y (1)Y (2) + Y (2)Y (1)
5 2 8 3 2 3 2 3 2
4 0 !#
1 1
+ Y (3)Y (0)
3 2

1 1 169 1 1 1 1 1 125
= + + + + =
5 16 62208 8 15552 972 486 1944 995328
y(x) = Y (0)x0 + Y (1)x1 + Y (2)x2 + Y (3)x3 + . . .

1 1
= 1x0 + 1x1 + 0x2 − x3 + 0x4 + x5 + . . .
3! 5!
1 2 1 3 169 4 125 5
=1+x+ x + x + x + x + ...
3 24 62208 995328
30
Capítulo 5
Aplicación del DTM a algunos modelos

Biológicos
A continuación estudiaremos la aplicación del método de transformación diferencial a

modelos de tipo biomatemáticos en los cuales se presentan retardos de diferente tipo.
5.1. Modelo para crecimiento de tumores en ratones

El crecimiento de tumores de tipo ascities Ehrlich2 en ratones se modela por medio del
siguiente problema de valor inicial con retardo

y 0 (t) = ry(t − α) 1 − y(t − α)

C (5.1)
y(t) = y0 (t) ≥ 0

para todo t ∈ [−α, 0]. Acá, y(t) se relaciona con la cantidad de células (concentración) en el
ratón; r es la tasa de proporcionalidad (neta) de reproducción de células del tumor; C es la
capacidad de almacenamiento y α es el retardo que muestra la duración de un ciclo de la
multiplicidad de células.
5.1.1. Caso particular

Consideremos a continuación el siguiente caso particular del modelo (5.1)

y 0 (t) = 0.1y t t

1−y
2 2 (5.2)
y(0) = 0.1

Podemos ver que la ecuación diferencial asociada a este problema puede escribirse de
manera tal que los teoremas del DTM para ecuaciones con retardo sean de fácil aplicación,
es decir,
t t t
y(t) = 0.1y − 0.1y y
2 2 2
Ahora, aplicando DTM a esta última expresión llegamos a la ecuación de recurrencia
k
(k + 1)! 1 X1 1
Y (k + 1) = 0.1 k Y (k) − 0.1 Y (l)Y (k − l)
k! 2 l=0
2l 2k−l
2
Paul Ehrlich (1854-1915). Médico y Bacteriólogo alemán, ganador del premio Nobel de Medicina en 1908.
Capítulo 5. Aplicación del DTM a algunos modelos Biológicos 31
" k
#
k! 0.1 X 1
Y (k + 1) = Y (k) − 0.1 Y (l)Y (k − l) (5.3)
(k + 1)! 2k l=0
2k
De la condición inicial y(0) = 0.1 tenemos que Y (0) = 0.1. Ahora, iterando la ecuación (5.3)
obtenemos:
Si k = 0, entonces

0! 0.1 1
Y (1) = Y (0) − 0.1 0 Y (0)Y (0)
1! 20 2
= 1 [0.1Y (0) − 0.1Y (0)Y (0)] = (0.1)(0.1) − (0.1)3 = 0.009
Si k = 1, entonces

1! 0.1 1 1
Y (2) = Y (1) − 0.1 1 Y (0)Y (1) + 1 Y (1)Y (0)
2! 21 2 2

1 1 1
= 0.05Y (1) − 0.1 Y (0)Y (1) + Y (1)Y (0)
2 2 2
1
= [0.05(0.009) − 0.1 (0.5(0.1)(0.009) + 0.5(0.009)(0.1))] = 0.00018
2
Si k = 2, entonces

2! 0.1 0.1
Y (3) = Y (2) − 2 (Y (0)Y (2) + Y (1)Y (1) + Y (2)Y (0))
3! 22 2

1 0.1 0.1
= Y (2) − (Y (0)Y (2) + Y (1)Y (1) + Y (2)Y (0))
3 4 4

1 0.1 2

= 0.025(0.00018) − 0.1(0.00018) + (0.009) + (0.00018)(0.1)
3 4
= −0.00333416
Si k = 3, entonces

3! 0.1 0.1
Y (4) = Y (3) − 3 (Y (0)Y (3) + Y (1)Y (2) + Y (2)Y (1) + Y (3)Y (0))
4! 23 2

1 0.1
= 0.0125(−0.00333416) − (0.1(−0.00333416) + (0.009)(0.00018) + (0.00018)(0.009)
4 8

+(−0.00333416)(0.1)) = −0.00000928326
y(x) = Y (0)x0 + Y (1)x1 + Y (2)x2 + Y (3)x3 + Y (4)x4 + . . .

= 0.1 + 0.009x + 0.00018x2 − 0.00333416x3 − 0.00000928326x4 + . . .
Figure 3. Example 4.2 of Oberle and Pesch.

1.2
0.8
0.6
0.4
0.2
−0.2
−0.4
λ = Solución exacta
−0.6
λ =DTM
−0.8
0 0.5 1 1.5 2 2.5
F IGURA 5.1: Modelo para crecimiento de tumores.
5.2. Diseminación de infecciones

Según el modelo de Kermack-McKendrick, el número de nuevas infecciones por unidad
de tiempo es proporcional al producto y1 (t)y2 (t), donde y1 (t) es la fracción susceptible de la
población en estudio; y2 (t) es la fracción infectada y y3 (t) representa la fracción inmunizada.
Ahora, asumiendo que el nuevo número de personas inmunizadas es proporcional al de

las infectadas, obtenemos el modelo

0
y1 (t) = −y1 (t)y2 (t)

y20 (t) = y1 (t)y2 (t) (5.4)

 0
y3 (t) = y2 (t)
Acá, por comodidad suponemos que todas las constantes de proporcionalidad son iguales
a 1. De igual forma podemos asumir que la población inmunizada vuelve a ser susceptible
después de un periodo fijo de tiempo que llamaremos α1 . Así pues, de igual forma podemos
introducir un periodo de incubación α2 obteniendo el modelo

0
y1 (t) = −y1 (t)y2 (t − α1 ) + y2 (t − α2 )

y20 (t) = y1 (t)y2 (t − α1 ) − y2 (t) (5.5)

 0
y3 (t) = y2 (t) − y2 (t − α2 )
5.2.1. Caso particular

Consideremos a continuación un caso especial del modelo (5.5)

0
y1 (t) = −y1 (t)y2 (t − 1) + y2 (t − 10)

y20 (t) = y1 (t)y2 (t − 1) − y2 (t) (5.6)

 0
y3 (t) = y2 (t) − y2 (t − 10)
sujeto a las condiciones iniciales y1 (0) = 0, y2 (0) = 1, y3 (0) = 5. Aplicando los teoremas
del DTM para ecuaciones con retardo vistos en el capítulo anterior tenemos:
k N N
(k + 1)! X X
l h1 −k+l h1 X
l−k l
Y1 (k + 1) = − (1) (−1) Y1 (l)Y2 (h1 ) + (−10) Y2 (l)
k! l=0 h1 =k−l
k−l l=k
k
k N
(k + 1)! X X
l h1 −k+l h1
Y2 (k + 1) = (1) (−1) Y1 (l)Y2 (h1 ) − Y2 (k)
k! l=0 h1 =k−l
k − l
N
(k + 1)! X
l−k l
Y3 (k + 1) = Y2 (k) − (−10) Y2 (l)
k! l=k
k
Después de simplificar obtenemos el sistema de ecuaciones de recurrencia:
" k N N #
1 X X h1
X l
Y1 (k + 1) = − (−1)h1 −k+l Y1 (l)Y2 (h1 ) + (−10)l−k Y2 (l)
k+1 l=0 h1 =k−l
k − l l=k
k
" k N #
1 X X
h1 −k+l h1
Y2 (k + 1) = (−1) Y1 (l)Y2 (h1 ) − Y2 (k)
k + 1 l=0 h =k−l k−l
1
" N #
1 X l
Y3 (k + 1) = Y2 (k) − (−10)l−k Y2 (l) (5.7)
k+1 l=k
k
En la figura 5.2 podemos ver la eficacia del método de transformación diferencial aplica-
do al sistema de ecuaciones diferenciales con retardo (5.6).
5.3. Dinámica poblacional

A continuación presentamos el modelo en el que la población de cierta especie dada por
la función y(x) donde su tasa de crecimiento (creciente o decreciente) y está denotada por α.
Así, si la población se incrementa tendremos entonces una escasez de alimento al igual que
espacio en el hábitat.
En general, el comportamiento de este tipo de fenómenos es modelado por el siguiente

problema de valor inicial
(
y 0 (x) = α(β − y(x))y(x), x ≥ 0
(5.8)
y(x0 ) = y0
Sol.Exacta
5
Sol.Exacta
Sol.Exacta
DTM
4 DTM
0
0 1 2 3 4 5 6 7 8 9 10
F IGURA 5.2: Diseminación de infecciones.
5.3.1. Modelo de predicción

Asumiendo que la tasa de crecimiento es dependiente de la población en las generaciones
precedentes, obtenemos el modelo con retardo
y 0 (x) = α(β − y(x − τ ))y(x) (5.9)
Ahora, haciendo el cambio de variable z(x) = ατ y(τ x) tenemos que
z 0 (x) = y 0 (τ x)(ατ 2 )
De acuerdo a y(τ x) vemos que
y 0 (τ x) = α(β − y(τ (x − 1)))y(τ x)
z 0 (x) = (ατ β − ατ y(τ (x − 1)))ατ y(τ x)

= (ατ β − z(x − 1))z(x)
Finalmente, retornando el cambio de variable de z por y y ατ β por β llegamos a
y 0 (x) = (β − y(x − 1))y(x) (5.10)
5.3.2. Caso particular 1.

Como caso particular de (5.10) tenemos el siguiente modelo

y 0 (t) = 3.5y(t) 1 − y(t − 0.74)

19 (5.11)
y(0) = 19.001

Entonces, para aplicar el DTM reescribimos la ecuación diferencial con retardo en la forma
y 0 (t) = 3.5y(t) − 0.184211y(t)y(t − 0.74)
Aplicando a esta ecuación los teoremas del DTM obtenemos:

k N !
(k + 1)! X X h1
Y (k + 1) = 3.5Y (k) − 0.184211 (1)l (−0.74)h1 −k+l Y (l)Y (h1 )
k! l=0 h =k−l
k−l
1
donde después de algunas simplificaciones obtenemos la ecuación de recurrencia

" k N #
1 X X h1
Y (k + 1) = 3.5Y (k) − 0.184211 (−0.74)h1 −k+l Y (l)Y (h1 ) (5.12)
k+1 l=0 h =k−l
k−l
1
Acá, de la condición inicial y(0) = 19.001 tenemos el dato Y (0) = 19.001.
20
DTM vs. Sol. Exacta
19.5
19
y(t)
18.5
18
17.5
0 5 10 15
tiempo t
F IGURA 5.3: Dinámica poblacional.
5.3.3. Caso particular 2.

Consideremos ahora otro modelo particular de la forma (5.10), es decir,
(
y 0 (t) = −3y(t − 1) (1 + y(t))
(5.13)
y(0) = 0.1
Como en el caso anterior, podemos escribir esta ecuación de forma tal que los teoremas
del DTM sean aplicables sin problema alguno, es decir, el problema (5.13) adopta la forma
y 0 (t) = −3y(t − 1) − 3y(t)y(t − 1)
Ahora, aplicando DTM a esta última ecuación obtenemos

N k N
(k + 1)! X
l−k l
X X
l h1 −k+l h1
Y (k + 1) = −3 (−1) Y (l) − 3 (1) (−1) Y (l)Y (h1 )
k! l=k
k l=0 h =k−l
k−l
1
donde después de algunas simplificaciones llegamos a la ecuación de recurrencia
" N k N #
−3 X l−k l
X X
h1 −k+l h1
Y (k + 1) = (−1) Y (l) + (−1) Y (l)Y (h1 ) (5.14)
k + 1 l=k k l=0 h =k−l
k−l
1
Acá, la condición inicial y(0) = 0.1 implica que Y (0) = 0.1. Después de iterar la ecuación de
recurrencia (5.14), obtenemos una excelente aproximación con respecto a la original evalua-
da con técnicas computacionales. (Ver figura 5.4)
1.2
0.8
0.6
0.4
0.2
−0.2
−0.4
−0.6 λ = 3 (Solución exacta)

λ =DTM
−0.8
0 0.5 1 1.5 2 2.5
F IGURA 5.4: Dinámica poblacional.

37
Conclusiones
En este trabajo mostramos las propiedades mas relevantes del método de transformación
diferencial utilizando su definición. Presentamos las pruebas de cada uno de los teoremas
del DTM para ecuaciones diferenciales ordinarias sin retardo, al igual que las propiedades
a partir de la transformada inversa con la expansión en series infinitas de Taylor, su linea-
lidad, la relación existente de la transformada diferencial de una función con las derivadas
de orden n, al igual que las funciones trascendentes de mayor aplicación.
En la aplicación del DTM a ecuaciones diferenciales con retardo utilizamos inicialmen-

te ejemplos elementales que permitieran ilustrar y comparar su solución con la solución
exacta hallada por métodos analíticos, visualizando así el comportamiento en cuanto a la
convergencia de las soluciones. Aplicamos posteriormente el DTM para problemas con más
alto grado de complejidad, como los problemas de valor inicias de tipo integro-diferencial
con diferentes tipos de retardo. Generalizamos el método para poder aplicarlo a sistemas
de ecuaciones diferenciales con retardo, mostrando nuevamente la rápida convergencia en
comparación con su solución exacta.
Durante esta investigación, propusimos el teorema 4.2.6 y los corolarios 4.2.5 y 4.2.6,
fundamentales a la hora de aplicar el DTM a ecuaciones con retardo en el cual la presencia
de términos cuadráticos del tipo g 2 (x − α) hacen que los teoremas conocidos no sean de fácil
aplicación. Presentamos igualmente la prueba de dicho teorema y con base a ello se puede
fácilmente demostrar los corolarios mencionados.
Utilizando Matlab y Mathematica, pudimos mostrar que la aplicación del DTM a pro-
blemas biológicos en los cuales había presencia de diferentes tipo de retardo es muy eficien-
te; esto se pudo comprobar con los resultados obtenidos por diferentes investigadores en
el campo de la Biomatemática. Las ecuaciones de recurrencia mostraron que los teoremas
expuestos en este trabajo hacen del DTM una herramienta semi-analítica de muy rápida
convergencia en comparación con las aproximaciones hechas con otro tipo de técnicas. Esto
puede verse en el último capítulo donde mostramos la aplicación a diferentes tipo de mode-
los biológicos.
El método de transformación diferencial induce a la investigación futura a problemas que

presenten ecuaciones con retardos mas generales (sistemas de ecuaciones diferenciales con
retardo) y sobre todo a temas como la estabilidad de algunas de las soluciones obtenidas, al
igual que el diseño de programas de computación propios para el DTM. Estos se presentan
en el panorama del mismo como problemas abiertos, así como su relación con otros métodos
tales como el método de Adomian, la perturbación homotópica, etc.
38
Bibliografía
[1] Arboleda, A. (2012). Resolución de Ecuaciones Diferenciales No Lineales por el Método de

Transformación Diferencial, Tesis de Maestría en Enseñanza de la Matemática, Universi-
dad Tecnológica de Pereira, Colombia.
[2] Cárdenas A, P.P. (2013).The Zhou’s method for solving the white-dwarfs equation. Applied
Mathematics., 10C(4):28-32.
[3] Cárdenas A, P.P.(2014). An Iterative Method for Solving Two Special Cases of Lane-Emden
Type Equation. American Journal of Computational Mathematics., 5(4): 243-253.
[4] Bani, M. and Amundsen, D. (2014). Oscillatory traveling waves for a population diffusion
model with two age and nonlocality induced by maturation delay. Applied Mathematics., 34:
309-324.
[5] Abazari, R. and Kilicman, A. (2014). Application of Differential Transform Method on Non-
linear Integro-Differential Equations with Proportional Delay. Neural Computing and Ap-
plications., 24: 391-397.
[6] Biazar, J., Islam, M. and Islam, M. (2012). Differential Transform Method for Special Systems
of Integral Equations. Journal of King Saud University (Science)., 24: 211-214.
[7] Liu, B., Zhou, X. and Du, Q. (2015). Differential Transform Method for Some Delay Differen-
tial Equations. Applied Mathematics., 6:585-593.
[8] Raslan, K., and Zain, F. (2013). Comparison study between diferential transform method and
adomian decomposition method for some delay differential equations. Int. J. Phys. Sci., 17:744-
2013.
[9] Fadhel, F.S., and Mohammed, J.(2011). Extend diferential transform methods for solving
diferential equations with multiple delay., Applied Mathematics., 24(03):10.
[10] Mirzaee, F., and Latifi, L.(2011). Numerical Solution of Delay Differential Equations by Diffe-
rential Transform Method., Mathematics Department, Malayer University, Malayer, Iran.
Vol 20, No. 78/2.
[11] Liz, E.(2006). Sobre ecuaciones diferenciales con retraso, dinámica de poblaciones y números
primos., Materials Matematics, page 24.
[12] Shampine, L., and Thompson, S.(2000). Solving Delay Differential Equations with dde23.
Mathematics Department, Southern Methodist University. Dallas, TX 75275.
[13] Kumar, A., and Arora, R. (2013). Solutions of the coupled system of Burgers’ equations and
coupled Klein-Gordon equation by RDT Method., International Journal of Advances in Ap-
plied Mathematics and Mechanics, 1(2): 103-115.
BIBLIOGRAFÍA 39
[14] Ayaz, F. (2004). Solutions of the system of differential equations by differential transform met-
hod., Applied Mathematics and Computation. 147:547-567.
[15] Blanco, L., and Estrella E. (2013). Solving delay differential systems with history functions
by the Adomian decomposition method., Appl. Math. Comput. 218:5994-6011.
[16] Smarda, Z., and Diblik, J.(2013). Extension of the differential transformation method to non-
linear differential equations and integro-differential equations with proportional delays., Ad-
vances in Difference Equations. 69:1-12.
[17] Valli, D., and Kurths, J. (2014). Synchronization in coupled Ikeda delay systems., The Euro-
pean Physical Journal. 223:1465-1479.
[18] Karako, F., and Bereketoglu, H. (2009). Solutions of Delay Differential Equations by Using
Differential Transform Method., International Journal of Computer Mathematics. 86:914-
923.
[19] Candan, T., and Dahiya, R.(2004). Oscillation behavior of n-th order neutral differential equa-
tions with continuous delay., Journal of Mathematical Analysis and Applications. 290:
105-112.
[20] Zhou, J.(1986). Differential transformation and its application for electrical circuits., Huar-
jung University Press, Wuuhahn, China.

APROX - NUMERlCA DE EDR

Cargado por

Copyright:

Formatos disponibles

APROX - NUMERlCA DE EDR

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

APROX - NUMERlCA DE EDR

Cargado por

Copyright:

Formatos disponibles

U NIVERSIDAD T ECNOLÓGICA DE P EREIRA

Autor: Director de tesis:

Trabajo presentado como requisito para optar

Grupo de Investigación en Ecuaciones Diferenciales No Lineales GEDNOL

Llegar a esta instancia en mis estudios de maestría ha sido posible obviamente a mi

Iniciaré agradeciendo a mi familia, especialmente a mis padres por su incondicional

Agradezco también a cada uno de los profesores que participaron en mi formación, en

Jhon Fredy Trujillo Valencia

2. Generalidades de las ecuaciones diferenciales con retardo 5

3. Método de Transformación Diferencial DTM 10

4. DTM para ecuaciones diferenciales con retardo 16

5. Aplicación del DTM a algunos modelos Biológicos 30

... a mis padres

En esta investigación mostraremos inicialmente un desarrollo teórico sobre el método de

Muchos autores han investigado sobre el problema de encontrar soluciones a ecuacio-

El método de transformación diferencial consiste básicamente en un proceso de tipo

Todos los acontecimientos de la naturaleza requieren la existencia de la variable tem-

A nivel de software se han realizado diversas simulaciones, la cuales permiten corroborar

Desarrollar el método de transformación diferencial para ecuaciones diferenciales con

Realizar simulaciones que evidencien la rapidez de convergencia del método.

1.4.1. Planteamiento del problema

La tercera etapa consistió en la consulta a investigadores, principalmente aquellos rela-

Generalidades de las ecuaciones

2.1. Ecuaciones Diferenciales con retardo (EDR)

Presentamos a continuación algunas definiciones esenciales en el estudio de las EDR.

Definición 2.1.1. (Espacio de estados). Denominamos el espacio de estados, representado por C al

Definición 2.1.2. (Variable de estado). Sea la función continua x : R → Rn , t real y C un espacio

x0 (t) = f (t, x(t), x(t − r))

donde la función f es continua con parámetro r > 0 denominado el retardo de la ecuación.

Definición 2.1.4. Definimos problema de valor inicial de una EDR como

Es de anotar que en las ecuaciones diferenciales con retardo, el comportamiento de sus

en la cual xt (s) = x(t + s) con xt : [−r, 0] → R.

con xt (φ) := x(t + φ), φ ∈ [−τ, 0].

2.2. Modelo matemático con retardo: Caso especial

x0 (t) = −αx(t − r) (2.3)

Por ejemplo, si tenemos la ecuación diferencial con retardo

x0 (t) = −0.2x(t − 1),

Ahora bien, supongamos la ecuación con retardo

x0 (t) = −1.5x(t − 1),

F IGURA 2.1: x0 (t) = −0.2x(t − 1).

F IGURA 2.2: x0 (t) = −1.5x(t − 1).

Finalmente, supongamos la ecuación con retardo

Así, si (α)r ≥ π2 , tendremos que el comportamiento de las soluciones es periódico.

2.3. Sobre la existencia y unicidad de soluciones

kf (t, φ) − f (t − ϕ)k ≤ Kkφ − ϕk, (2.4)

para todo t ∈ [a, b] y kφk, kϕk ≤ M .

donde kϕk, kψk ≤ M .

donde integrando obtenemos

Este resultado puede generalizarse diciendo que para todo m ≥ 1,

F IGURA 2.3: Intersección con las soluciones retardadas

Método de Transformación Diferencial

Fue utilizado inicialmente por Zhou en la solución de problemas en circuitos eléctricos

La transformación diferencial, al igual que todo un grueso de transformadas (Fourier, La-

3.2. Generalidades del DTM

El valor de la función y(x) es obtenido a partir de la transformación diferencial inversa de

Podemos observar que al sustituir Y (k) en la transformada diferencial inversa obtenemos

es decir, la solución corresponde a la serie de Taylor en x0 de la función y(x). Ahora bien, si

La transformación diferencial cumple con propiedades de linealidad al igual que otras