Semana 02 Teoría

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 27

Ecuaciones

diferenciales
Ecuaciones diferenciales de
primer orden: Parte I
Semana 02: Teoría
Profesores del curso:
Hermes Pantoja Carhuavilca
Sergio Quispe Rodríguez
Patricia Reynoso Quispe
Cristina Navarro Flores
Orlando Galarza Gerónimo
César Barraza Bernaola
César Vergaray Albujar

Profesores: Utec-Ciencias
Índice
1 Ecuaciones diferenciales
exactas
2 Modelos de ecuaciones
diferenciales de 1er orden
Objetivos

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales August 22, 2023 2 / 26


ECUACIONES
DIFERENCIALES
EXACTAS 1
Logros
?
?
Identifica si una EDO es exacta.
Resuelve una EDO exacta.
Identifica modelos de problemas reales expresados como EDOs.
Resuelve modelos de problemas reales expresados como EDOs usando el
método de variables separables.

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales August 22, 2023 4 / 26


Conceptos preliminares
Si z = f (x, y ) es una función de dos variables, entonces su diferencial está dado
por la expresión

∂f ∂f
dz = d(f (x, y )) = dx + dy
∂x ∂y
Si en particular f (x, y ) = C, donde C es una constante, entonces

∂f ∂f
dx + dy = 0 (1)
∂x ∂y
Dada una ecuación diferencial de la forma

M(x, y )dx + N(x, y )dy = 0 (2)

Podemos hacer una comparación entre las ecuaciones (1) y (2).

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales August 22, 2023 5 / 26


Ecuaciones exactas
Dada la ecuación diferencial
M(x, y )dx + N(x, y )dy = 0 (3)
Se puede argumentar que, de forma general, la solución está definida implícitamente
por f (x, y ) = C, siendo C una constante. Calculando el diferencial de esta última
expresión, obtenemos:
∂f ∂f
dx + dy = 0. (4)
∂x ∂y
Por lo tanto, si se cumple que

∂f (x, y )
M(x, y ) = (5)
∂x
∂f (x, y )
N(x, y ) = (6)
∂y
entonces la expresión f (x, y ) = C define la solución de (3).
Universidad de Ingeniería y Tecnología Ecuaciones diferenciales August 22, 2023 6 / 26
Derivando (5) y (6) respecto de y y x, respectivamente

∂M ∂2f
=
∂y ∂y ∂x
∂N ∂2f
=
∂x ∂x∂y
∂M ∂N
⇒ = ¡Condición importante!
∂y ∂x

Algunos comentarios
La condición importante es una consecuencia de suponer que la solución de
la ecuación (3) es de la forma f (x, y ) = C.
¿Será posible aplicar este método a cualquier ecuación de la forma dada por
(3)?
Dada una ecuación como (3), ¿Cómo hallamos la función f (x, y )?

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales August 22, 2023 7 / 26


Definición
La ecuación diferencial
M(x, y )dx + N(x, y )dy = 0 (7)
es exacta si se satisface que

∂M ∂N
= ← Condición importante
∂y ∂x

Solución:
Si encontramos una función f (x, y ), que cumple con

∂f (x, y ) ∂f (x, y )
M(x, y ) = , N(x, y ) = . (8)
∂x ∂y

Entonces, la solución de la ecuación (7) está dada por f (x, y ) = C.

Siempre que la ecuación sea exacta, es posible hallar f (x, y ) a partir de las ecuaciones (8).

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales August 22, 2023 8 / 26


Ejemplo
Resolver la siguiente ED: yx + 6x dx + (ln x − 2) dy = 0


Solución:
1
Paso 1: Derivadas parciales: ∂M ∂N
∂y = ∂x = x . Luego, la ED es Exacta.
Paso 2: Escogemos una de las funciones, M o N. En este caso escogeremos M y la
integramos en función a x:
Z 
y 
f (x, y ) = + 6x dx = y ln x + 3x 2 + h(y )
x
Paso 3: Recordamos que h(y ) es constante para x. Luego, derivamos esta función
f (x, y ) en función a y .
∂f (x, y )
= ln x − 2 = N
∂y
Paso 4: Recordamos que esta derivada parcial en función a y es igual a N.
Z
ln x + h′ (y ) = ln x − 2 ⇒ h(y ) = −2dy = −2y

Luego : f (x, y ) = y ln x + 3x 2 − 2y = C

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales August 22, 2023 9 / 26


Ejercicio
Resolver la siguiente ED:
dy xy 2 − cos x sin x
=
dx y (1 − x 2 )
Solución:
Paso 1: Levamos a su forma equivalente e identificamos M(x, y ) y N(x, y )

(cos(x) sin(x) − xy 2 )dx + y (1 − x 2 )dy = 0


⇒ M(x, y ) = cos(x) sin(x) − xy 2
⇒ N(x, y ) = y (1 − x 2 )
Paso 2: Verificamos si es exacta
∂M ∂
= (cos(x) sin(x) − xy 2 ) = −2xy
∂y ∂y
∂N ∂ 
y (1 − x 2 ) = −2xy

=
∂x ∂x
La ecuación es exacta.
Universidad de Ingeniería y Tecnología Ecuaciones diferenciales August 22, 2023 10 / 26
? Paso 3: Debido a que la ecuación es exacta, existe una función f (x, y ) tal que
∂f (x, y )
= M(x, y ) = cos(x) sin(x) − xy 2 (9)
∂x
∂f (x, y )
= N(x, y ) = y (1 − x 2 ) (10)
∂y
Paso 4: Integramos la ecuación (9) respecto de x:
Z
1 1
f (x, y ) = (cos(x) sin(x) − xy 2 )dx + g(y ) = sin2 (x) − x 2 y 2 + g(y ) (11)
2 2

Paso 5: Derivamos la ecuación (11) respecto de y y usamos la ecuación (10):


 
∂f ∂ 1 2 1 2 2 dg(y ) dg(y )
= sin (x) − x y + g(y ) = −x 2 y + = y −x 2 y ⇒ =y
∂y ∂y 2 2 dy dy

Paso 6: Hallamos g(y ) integrando respecto de y :


y2
g(y ) = (12)
2
Universidad de Ingeniería y Tecnología Ecuaciones diferenciales August 22, 2023 11 / 26
? Paso 7: Reemplazando (12) en (11) obtenemos f (x, y ):

1 1 y2
f (x, y ) = sin2 (x) − x 2 y 2 +
2 2 2
Finalmente, la solución del problema es

1 1 y2
sin2 (x) − x 2 y 2 + =C
2 2 2
Paso 8: Usamos la condición inicial
1 1 22
sin2 (0) − (0)2 (2)2 + =C ⇒ C=2
2 2 2
Resultando la solución
1 1 y2
sin2 (x) − x 2 y 2 + =2
2 2 2

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales August 22, 2023 12 / 26


Dinámica Poblacional y Decaimiento Radiactivo
El modelo de dinámica poblacional propuesto por el inglés Tomás Malthus indica
que la tasa de crecimiento de una población (que no se ve afectada por algún
medio externo) es proporcional a la población existente en un determinado
momento (bacterias que crecen en una placa de Petri, por ejemplo).
El modelo de decaimiento radiactivo supone que la velocidad dA/dt a la que
se desintegran los núcleos de una sustancia es proporcional a la cantidad
(más exactamente, al número de núcleos) A(t) de la sustancia que queda en el
tiempo t:

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales August 22, 2023 13 / 26


2
MODELOS DE
ECUACIONES
DIFERENCIALES
DE 1er ORDEN
Dinámica Poblacional y Decaimiento Radiactivo
En las gráfica de la izquierda se ve el crecimiento de una población determinada
sin ninguna restricción ni agente externo que afecte su crecimiento.
En la grafica de la derecha se ve la cantidad de material radiactivo en el
transcurso del tiempo.

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales August 22, 2023 15 / 26


Ejercicios
1 La población de un pueblo crece con una razón proporcional a dicha población en
el instante t. La población inicial de 500 aumenta 15% en 10 años. ¿Cuál será la
población pasados 30 años? ¿Qué tan rápido está creciendo la población en t = 30?.
Solución:
Paso 1: Construir la ecuación diferencial
dP
= kP
dt
Paso 2: Indentificar la condición inicial y el dato adicional
P(0) = 500 P(10) = 500 + 15%(500) = 575
Paso 3: Separar las variables
dP
= kdt
P
Paso 4: Aplicar integrales
Z Z
1
dP = kdt ⇒ P = Aekt Continuar ...
P
Universidad de Ingeniería y Tecnología Ecuaciones diferenciales August 22, 2023 16 / 26
2 El isótopo radioactivo del plomo Pb-209 decae con una razón proporcional a la cantidad
presente en el tiempo “t” y tiene una vida media de 3.3 horas. Si al principio había 1
gramo de plomo. ¿Cuánto tiempo debe transcurrir para que decaiga 90%?
Solución:
Paso 1: Construir la ecuación diferencial: dX dt = kX
Paso 2: Indentificar la condición inicial y el dato adicional: X (0) = 1, X (3.3) = 0.5
Paso 3: Separar las variables y aplicar integrales:
Z Z
1
dX = kdt ⇒ ln |X | = kt + C1
X
⇒ eln|X | = ekt+C1 ⇒ |X | = ekt+C1
⇒ X = Aekt

Paso 5: Usar las condiciones del paso 2:

X (0) = Ae0 = 1 ⇒ A=1


1
X (3.3) = e3.3k = 0.5 ⇒ k= ln(0.5) = −0.21
3.3

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales August 22, 2023 17 / 26


⇒ X (t) = e−0.21t

? Paso 6: Se busca el tiempo “t” tal que la masa X (t) de plomo en ese intante se
el 10% de la cantidad inicial (decae el 90%).

0.1 = e−0.21t ⇒ ln(0.1) = ln e−0.21t




ln(0.1) = −0.21t ⇒ t = 10.96

Finalmente, debe transcurrir aproximadamente 11 horas para que decaiga el 90%


del plomo inicial.

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales August 22, 2023 18 / 26


3 Datado mediante C-14. Este procedimiento se basa en el hecho de que el isótopo
C-14 se va desintegrando de un cuerpo desde el día de su fallecimiento. El modelo
de decaimiento radioactivo que se basa en el modelo de Malthus nos indica que la
tasa de desitegración del C-14 es proporcional a la cantidad de C-14 presente en ese
instante y se sabe que la vida media del C-14 es igual a 5730 años.
Si se encuentra un hueso fosilizado que contiene 0.1% de su cantidad original de C-14,
determine la edad del fósil.
Solución:
Paso 1: Identificar el modelo. De acuerdo con el enunciado y el resultado del
ejercico anterior, la masa del C-14 está dado por

A(t) = A0 ekt , donde k es la constante de decaimiento

Paso 2: Determinar el valor de k . A partir del dato de vida media del C-14,

1
A0 = A0 e5730k
2  
1 1 1
⇒ = e5730k ⇒ k= ln
2 5730 2

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales August 22, 2023 19 / 26


? Paso 3: Se busca el tiempo “t” tal que la masa A(t) del C-14 en ese instante sea
el 0.1% de la cantidad inicial A0 :
0.1 ln(0.5)
A0 = A0 e 5730 t
100  
5730 0.1
⇒t = ln = 57104
ln(0.5) 100

Finalmente, se obtiene que la edad del fósil es de 57104 años, aproximadamente.

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales August 22, 2023 20 / 26


Ley de enfriamiento/calentamiento de Newton
?
La velocidad a la que cambia la temperatura de un cuerpo es proporcional a la
diferencia entre la temperatura del cuerpo y la temperatura del medio circundante,
la llamada temperatura ambiente.
A
Si T (t) representa la temperatura de un cuerpo en el instante t, Tm la temperatura
del ambiente y dT /dt la razón con la cual la temperatura del cuerpo cambia,
entonces

dT
= k (T − Tm ) (13)
dt
Para analizar:
En cualquier caso, enfriamiento o calentamiento, si Tm es una constante se tiene
que k < 0.

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales August 22, 2023 21 / 26


Ejemplos
Un cuerpo que tiene una temperatura de 70◦ F es depositado (en el instante t = 0) en un
lugar donde la temperatura se mantiene a 40◦ F . Después de 3 minutos, la temperatura del
cuerpo ha disminuido a 60◦ F .
¿Cuál es la temperatura del cuerpo después de 5 minutos?
¿Cuánto tiempo pasará para que el cuerpo tenga 50◦ F de temperatura?
Solución:
Paso 1: Se plantea la ecuación diferencial que rigue este fenómeno.
dT
= k (T − Tm )
dt
Paso 2: Indentificar la condición inicial y los datos adicionales.

T (0) = 70◦ F , Tm = 40◦ F , T (3) = 60◦ F

Paso 3: Resolver la ecuación del paso 1. Sugerencia: Use separación de variables.

T (t) = Cekt + Tm

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales August 22, 2023 22 / 26


? Paso 4: Utilizar los datos del paso 2.

T (t) = Cekt + 40
⇒ T (0) = C + 40 = 70 ⇒ C = 30

Además
T (3) = 30e3k + 40 = 60
 
3k 1 20
⇒ 30e = 20 ⇒ k = ln = −0.1352
3 30
Por lo tanto, la temperatura del cuerpo está dado por:

T (t) = 30e−0.1252t + 40

Paso 5: Responda las preguntas.

T (5) = 55.26◦ F

T (t) = 50 ⇒ t = 8.13 min

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales August 22, 2023 23 / 26


Solución

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales August 22, 2023 24 / 26


Conclusiones
?
?
1 Es importante determinar si una EDO es lineal o no y si lo es, usando el
factor integrante, se convierte en una EDO de variables separables y es
mucho más sencilla su solución.
2 Existen varios modelos de problemas reales que pueden ser expresados
como ecuaciones diferenciales
3 Los modelos vistos hasta ahora, como el crecimiento poblacional o el
enfriamiento de Newton, pueden ser resueltos usando separación de
variables.

Universidad de Ingeniería y Tecnología Ecuaciones diferenciales August 22, 2023 25 / 26


Gracias

También podría gustarte