Matemática Aplicada

Notas de Aula. Ecuaciones Diferenciales

Módulo III: EDOs

Dr. Leandro Giordano

ICO-UNGS Ecologı́a

Dr. Leandro Giordano (ICO-UNGS Ecologı́a) Matemática Aplicada 1 / 17

 Temario

1 Ámbito de Aplicación

2 Ejemplos de Aplicación

Dr. Leandro Giordano (ICO-UNGS Ecologı́a) Matemática Aplicada 2 / 17

 Temario

1 Ámbito de Aplicación

2 Ejemplos de Aplicación

Dr. Leandro Giordano (ICO-UNGS Ecologı́a) Matemática Aplicada 3 / 17

 Análisis de Sistemas Ambientales
Ecuaciones de conservación

Dr. Leandro Giordano (ICO-UNGS Ecologı́a) Matemática Aplicada 4 / 17

 Análisis de Sistemas Ambientales
Ecuaciones de conservación

Si se conocen I (t) y Q(t) ¿Cómo podrı́a hallarse x(t)?

Dr. Leandro Giordano (ICO-UNGS Ecologı́a) Matemática Aplicada 4 / 17

 Análisis de Sistemas Ambientales
Ecuaciones de conservación

Si se conocen I (t) y Q(t) ¿Cómo podrı́a hallarse x(t)?

Llave del problema: Ecuación de Conservación

dx
= I (t) − Q(t)
dt

Dr. Leandro Giordano (ICO-UNGS Ecologı́a) Matemática Aplicada 4 / 17

 Temario

1 Ámbito de Aplicación

2 Ejemplos de Aplicación

Dr. Leandro Giordano (ICO-UNGS Ecologı́a) Matemática Aplicada 5 / 17

 Concentración de sustancias contaminantes en un embalse

Dr. Leandro Giordano (ICO-UNGS Ecologı́a) Matemática Aplicada 6 / 17

 Concentración de sustancias contaminantes en un embalse

¿Cuáles son los flujos de entrada y salida?¿Cuál es la variable de estado?

Dr. Leandro Giordano (ICO-UNGS Ecologı́a) Matemática Aplicada 6 / 17

 Concentración de sustancias contaminantes en un embalse

Dr. Leandro Giordano (ICO-UNGS Ecologı́a) Matemática Aplicada 7 / 17

 Concentración de sustancias contaminantes en un embalse

Dr. Leandro Giordano (ICO-UNGS Ecologı́a) Matemática Aplicada 7 / 17

 Concentración de sustancias contaminantes en un embalse
Solución General

dQ
= rk + P − rc(t)
dt
dQ
= rk + P − rQ(t)/V
dt

Dr. Leandro Giordano (ICO-UNGS Ecologı́a) Matemática Aplicada 8 / 17

 Concentración de sustancias contaminantes en un embalse
Solución General

dQ
= rk + P − rc(t)
dt
dQ
= rk + P − rQ(t)/V
dt

dQ r
+ Q(t) = rk + P
dt V
Dr. Leandro Giordano (ICO-UNGS Ecologı́a) Matemática Aplicada 8 / 17
 Concentración de sustancias contaminantes en un embalse
Solución General

La ecuación

dQ r
+ Q(t) = rk + P
dt V

Es una EDO lineal de primer order, luego su solución general es (aplicando factor integrante):

R Z R
r /Vdt
Q(t)e = (rk + P)e r /Vdt dt
R Z R
Q(t)e r /V dt
= (rk + P) e r /V dt dt
Z
Q(t)e r /Vt = (rk + P) e r /Vt dt

V
Q(t)e r /Vt = (rk + P)e r /Vt + C
r
V (rk + P)
Q(t) = + Ce −r /Vt
r
Q(t) rk + P
= c(t) = + Ce −r /Vt
V r

Dr. Leandro Giordano (ICO-UNGS Ecologı́a) Matemática Aplicada 9 / 17

 Luego, si se sabe que c(t = 0) = c0 , para obtenerse la solución particular primero procedemos a
despejar el valor de la constante C:

(rk + P)
c0 = + Ce −r /V 0
r
(rk + P)
c0 = +C
r

Reordenando:

(rk + P)
C = c0 −
r

Dr. Leandro Giordano (ICO-UNGS Ecologı́a) Matemática Aplicada 10 / 17

 Luego, si se sabe que c(t = 0) = c0 , para obtenerse la solución particular primero procedemos a
despejar el valor de la constante C:

(rk + P)
c0 = + Ce −r /V 0
r
(rk + P)
c0 = +C
r

Reordenando:

(rk + P)
C = c0 −
r

Finalmente, reemplazando en la solución general, se obtiene:

(rk + P) (rk + P) −r /Vt

c(t) = + [c0 − ]e
r r

Dr. Leandro Giordano (ICO-UNGS Ecologı́a) Matemática Aplicada 10 / 17

 Luego, si se sabe que c(t = 0) = c0 , para obtenerse la solución particular primero procedemos a
despejar el valor de la constante C:

(rk + P)
c0 = + Ce −r /V 0
r
(rk + P)
c0 = +C
r

Reordenando:

(rk + P)
C = c0 −
r

Finalmente, reemplazando en la solución general, se obtiene:

(rk + P) (rk + P) −r /Vt

c(t) = + [c0 − ]e
r r

O lo que es lo mismo:

(rk + P) (rk + P) −r /Vt

c(t) = − e + c0 e −r /Vt
r r
(rk + P)
c(t) = (1 − e −r /Vt ) + c0 e −r /Vt
r

Dr. Leandro Giordano (ICO-UNGS Ecologı́a) Matemática Aplicada 10 / 17

 Proceso Demográfico

Supóngase que una población de una especie animal, en un sitio con

recursos ilimitados, presenta una tasa de reproducción de n individuos por
individuo y una tasa de mortalidad de m individuos por individuo, ambas
constantes en el tiempo. Sea P(t) la población en tiempo t ¿Cuál es la
ecuación diferencial que caracteriza el proceso?

Dr. Leandro Giordano (ICO-UNGS Ecologı́a) Matemática Aplicada 11 / 17

 Ası́, podemos plantear la siguiente ecuación de conservación

Dr. Leandro Giordano (ICO-UNGS Ecologı́a) Matemática Aplicada 12 / 17

 Ası́, podemos plantear la siguiente ecuación de conservación
dP
= nP(t) − mP(t)
dt
dP
= (n − m)P(t)
dt

Luego si k = n − m, se obtiene la ecuación del proceso Malthussiano, cuya resolución general la

vimos la clase pasada:
dP
= kP(t)
dt
dP
= kdt
P(t)
Z Z
dP
= kdt
P(t)
lnP(t) = kt + C
P(t) = e kt+C = e C e kt = Ce kt
Si para P(t = 0) tengo P0 ¿Cuál es la solución particular?
Dr. Leandro Giordano (ICO-UNGS Ecologı́a) Matemática Aplicada 12 / 17
 Para esto igualamos P(t) = P0 en el término correspondiente y evaluamos P(t) para t = 0 en el
restante:

P0 = Ce k0
C = P0

Y, de ahı́, la solución particular es

P(t) = P0 e kt

Dr. Leandro Giordano (ICO-UNGS Ecologı́a) Matemática Aplicada 13 / 17

 Modelos demográficos lineales

Finalmente, se nos puede postular que un sistema demográfico presenta entradas en función del
tiempo descritas por una función I (t) y salidas descritas por el producto de una función Q(t) y
P(t). Ası́ se obtiene:

Dr. Leandro Giordano (ICO-UNGS Ecologı́a) Matemática Aplicada 14 / 17

 Modelos demográficos lineales

Finalmente, se nos puede postular que un sistema demográfico presenta entradas en función del
tiempo descritas por una función I (t) y salidas descritas por el producto de una función Q(t) y
P(t). Ası́ se obtiene:

dP
= I (t) − Q(t)P(t)
dt

Dr. Leandro Giordano (ICO-UNGS Ecologı́a) Matemática Aplicada 14 / 17

 Modelos demográficos lineales

Finalmente, se nos puede postular que un sistema demográfico presenta entradas en función del
tiempo descritas por una función I (t) y salidas descritas por el producto de una función Q(t) y
P(t). Ası́ se obtiene:

dP
= I (t) − Q(t)P(t)
dt

Reordenado obtenemos:

dP
+ Q(t)P(t) = I (t)
dt

¿Qué tipo de EDO es y cómo se resuelve?

Dr. Leandro Giordano (ICO-UNGS Ecologı́a) Matemática Aplicada 14 / 17

 Por ejemplo si I (t) = I0 y Q(t) = k ¿cómo se procederı́a?

Dr. Leandro Giordano (ICO-UNGS Ecologı́a) Matemática Aplicada 15 / 17

 Por ejemplo si I (t) = I0 y Q(t) = k ¿cómo se procederı́a? Resuelva la EDO lineal y obtenga la
solución particular para P(t = 0) = P0

dP
= I0 − kP(t)
dt

dP
+ kP(t) = I0
dt Z
R R
kdt
P(t)e = I0 e kdt dt
Z
P(t)e kt = I0 e kt dt

e kt
P(t)e kt = I0 +C
k
I0
P(t) = + Ce −kt
k

Dr. Leandro Giordano (ICO-UNGS Ecologı́a) Matemática Aplicada 15 / 17

 Si se pide la solución particular para P(t = 0) = P0 :

I0
P0 = + Ce −k0
k
I0
P0 = + C
k
I0
C = P0 −
k
Finalmente:
I0 I0
P(t) = + (P0 − )e −kt
k k
I0
P(t) = ( − 1)e −kt + P0 e −kt
k

Dr. Leandro Giordano (ICO-UNGS Ecologı́a) Matemática Aplicada 16 / 17

 Supongamos que la cantidad de nuevos individuos en una población por
inmigración y natalidad queda dada mediante e −t/λ t y el decrecimiento
por emigración y mortalidad mediante kP(t)

¿Cuál es la ecuación diferencial (conservación) que describe el

proceso?

dP
= e −t/λ t − kP(t)
dt
dP
+ kP(t) = e −t/λ t
dt Z
R R
P(t)e kdt
= e −t/λ te kdt
dt
Z
P(t)e kt = e t(k−1/λ) tdt

Dr. Leandro Giordano (ICO-UNGS Ecologı́a) Matemática Aplicada 17 / 17