ACTIVIDAD EVALUATIVA EJE 3
MODELACION MATEMATICA
ESTUDIANTE:
ANGIE PAOLA GALLEGO MARIN
DOCENTE:
ANTONIO SALUSTIANO RODRIGUEZ PACHECO
FUNDACION UNIVERSITARIA DEL AREA
ANDINA
ECUACIONES DIFERENCIALES
INGENIERIA DE SISTEMAS
2023
Caso problema tarea eje 3 ED
Introducción
En muchas disciplinas del conocimiento, y por supuesto en Ingeniería, existen una gran
variedad de situaciones y problemas en los cuales se desea conocer un elemento variable
a partir de su coeficiente de variación, es decir, es necesario saber cómo varía dicho
elemento en función de una o varias variables. El propósito es determinar una función
desconocida mediante datos relacionados por una ecuación que contiene, por lo menos,
una de las derivadas de la función desconocida.
Los principios y leyes científicas que resultan de la experimentación se expresan en
ecuaciones matemáticas. Los Modelos Matemáticos son la idealización del problema físico
con el que nos encontramos y las ecuaciones diferenciales la forma de representarlos de
manera simplificada. Cada modelo matemático es una aproximación a la realidad del
problema físico, su aproximación y uso del modelo sólo depende de los criterios impuestos
a cada problema para su resolución. Si la intuición o la evidencia del experimento coinciden
con los resultados obtenidos por medio del modelo podremos determinar cuan útil es ese
modelo.
Las ecuaciones diferenciales ya formuladas deben ser resueltas, sujetas a condiciones
obtenidas del problema para determinar la incógnita o incógnitas involucradas, que en este
caso son funciones de variable real. Los procedimientos usados pueden producir una
solución exacta o soluciones aproximadas por medio del uso de cálculos numéricos con el
uso de las tecnologías de la comunicación y la información.
Situación problema resuelta. A continuación, se presenta una situación problema en la que
formulamos la ecuación diferencial, proponemos una ecuación diferencial y el método de
solución.
Una aplicación de las ecuaciones diferenciales: Desintegración de sustancias radiactivas.
Como sabemos, los núcleos de los átomos son composiciones combinadas de neutrones y
protones, pero por muchas razones, algunas de esas combinaciones son inestables hasta
llevar a la desintegración de los átomos. Esta condición se conoce como desintegración
radiactiva. La razón de cambio (o tasa de cambio) con que los núcleos de una sustancia se
desintegran es proporcional al número de núcleos de la sustancia restante en ciertotiempo
t. Si N es la cantidad de núcleos ya podemos sospechar (por la similitud de este modelo
matemático con las aplicaciones de los ejemplos anteriores), que la relación de
proporcionalidad nos queda (Figura 1):
𝑑𝑁
= 𝑘𝑁.
𝑑𝑡
Aunque el modelo matemático para la desintegración radiactiva que acabamos de deducir
es utilizado en el universo de las ciencias físicas, también se aplica en sistemas biológicos,
por ejemplo, para calcular cuánto puede ser la vida de un medicamento, es decir, una
medida de tendencia central que nos describa su vida media. En otras palabras, el tiempo
que tarda el organismo en eliminar la mitad de éste por los procesos biológicos de
metabolización o excreción.
Figura 2. Desintegración de sustancias radiactivas.
Fuente. http://biomodel.uah.es/tecnicas/isotopos/inicio.htm
Desintegración radiactiva1 del PB-209, isótopo radiactivo del plomo.
El PB-209, isótopo radiactivo del plomo, se desintegra con una rapidez proporcional a la
cantidad presente en algún instante y tiene una vida promedio de 3,3 horas2. Si inicialmente
se tiene un gramo de plomo, ¿cuánto tiempo debe transcurrir para que se desintegre un
80%?
Sea (𝑡) la cantidad de PB-209 (gramos) presente en el instante 𝑡 (horas).El
modelo matemático para la desintegración radiactiva es
𝑑𝑁
= −𝑘𝑁, 𝑘 > 0.
𝑑𝑡
1
Adaptación de un problema similar encontrado en http://ocw.uc3m.es/matematicas/ecuaciones-
diferenciales-ordinarias/pruebas-de-evaluacion-1/EDOexamsolene09.pdf
2
Esto significa que en cualquier momento hay el doble de isótopo que 3,3 horas después.
𝑑𝑁 = −𝑘𝑁 ⟹ 𝑑𝑁 = −𝑘𝑑𝑡 ⟹ 𝑙𝑛 𝑁 = 𝑘𝑡 + 𝑐 . Para
Resolviendo para 𝑁, tenemos 𝑑𝑡 𝑁
determinar la constante de integración 𝑐, es necesario establecer las condiciones iniciales,
es decir, para una población en determinado instante que tomamos como inicial, 𝑁 = 𝑁0
cuando 𝑡 = 0.
Volviendo al resultado preliminar 𝑙𝑛 𝑁 = 𝑘𝑡 + 𝑐 ⟹ 𝑙𝑛 𝑁0 = (0) + 𝑐 ⟹ 𝑐 = 𝑙𝑛 𝑁0 .
Entonces, 𝑙𝑛 𝑁 = 𝑘𝑡 + 𝑐 ⟹ 𝑙𝑛 𝑁 = 𝑘𝑡 + 𝑙𝑛 𝑁0 ⟹ 𝑒𝑙𝑛 𝑁 = (𝑘𝑡+𝑙𝑛 𝑁0 ) ⟹
𝑃 = 𝑒𝑘𝑡 ∙ 𝑒𝑙𝑛 𝑁0 .
En esta última expresión el término 𝑒𝑙𝑛 𝑁0 es una nueva constante 𝑃0, y llegamos a
𝑁 = 𝑁0𝑒𝑘𝑡.
Sabemos que en el instante inicial se tiene un gramo de plomo, es decir,
(0) = 1 ⟹ 1 = 𝑁0𝑒𝑘(0) ⟹ 1 = 𝑁0
Nos falta determinar 𝑘. Para hacerlo sabemos que la semivida del plomo es 3,3 horas, lo
que significa que después de 3,3 horas debe quedar la mitad del PB-209 del inicio, es decir,
(3,3) = 0,5. Entonces
1
1 (−3,3)𝑘 1 (−3,3)𝑘
1 = −3,3𝑘 ⟹ 𝑘 = 𝑙𝑛 2 𝑙𝑛 2
= ⟹ 𝑙𝑛 = 𝑙𝑛 ⟹ 𝑙𝑛 ⟹𝑘=
2 2 2 −3,3 3,3
𝑙𝑛 2
Al reemplazar los valores 1 = 𝑁0 y 𝑘= encontrados, la ecuación diferencial
3,3
particular nos queda
𝑙𝑛 2
− 𝑡
𝑁=𝑒 3,3
Con el modelo matemático encontrado, nos piden el tiempo que debe transcurrir para
que se desintegre un 80%, es decir, (𝑡) = 0,2 . Entonces
𝑙𝑛 2
− 𝑡 𝑙𝑛 2 𝑙𝑛 0,2 𝑙𝑛 20
0,2 = 𝑒 3,3 ⟹ 𝑙𝑛 0,2 = − 𝑡 ⟹ 𝑡 = 3,3 −𝑙𝑛 2 ⟹ 𝑡 = 3,3 ⟹ 𝑡 = 14,26.
3,3 𝑙𝑛 2
Es decir, deben pasar un poco más de 14 horas para que PB-209, isótopo radiactivo del
plomo se desintegre un 80%.
Problema propuesto situación 1
Una masa de m gramos cae verticalmente hacia abajo, bajo la influencia de la gravedad,
partiendo del reposo, siendo despreciable la resistencia del aire.
Usted debe realizar las siguientes etapas del problema:
1. Explicar las condiciones asociadas que describen el movimiento.
2. Formular la ecuación diferencial referente al problema.
3. Resolver la ecuación diferencial
Diagrama de fuerzas
Figura 2
Fuente. Propia
𝐴: 𝑝𝑜𝑠𝑖𝑐𝑖ó𝑛 𝑑𝑒 𝑙𝑎 𝑚𝑎𝑠𝑎 𝑚 𝑒𝑛 𝑒𝑙 𝑡𝑖𝑒𝑚𝑝𝑜 𝑡 = 0.
𝑃𝑖 ∶ 𝑙𝑎 𝑝𝑜𝑠𝑖𝑐𝑖ó𝑛 𝑑𝑒 𝑚 𝑒𝑛 𝑐𝑢𝑎𝑙𝑞𝑢𝑖𝑒𝑟 𝑡𝑖𝑒𝑚𝑝𝑜 𝑝𝑜𝑠𝑡𝑒𝑟𝑖𝑜𝑟 𝑡.
En cualquier problema de física que involucre cantidades vectoriales tales como fuerza,
desplazamiento, velocidad y aceleración, las cuales necesariamente requieren un
conocimiento de dirección, es conveniente establecer un sistema de coordenadas, junto
con la asignación de direcciones positivas y negativas. En este ejemplo observamos que la
variación se realiza respecto del eje x.
Teniendo en cuenta las condiciones iniciales, usted debe mostrar el procedimiento para
llegar a la ecuación diferencial de segundo orden
𝑑2𝑥
= 𝑔.
𝑑𝑡2
Sus soluciones son
𝑣 = 𝑔𝑡
1
𝑥= 𝑔𝑡2
2
Sea A en la figura la posición de la masa 𝑚 en el tiempo 𝑡 = 0, y sea 𝑃𝑖 la posición de 𝑚 en cualquier
tiempo posterior 𝑡. En cualquier problema de física que involucre cantidades vectoriales tales como
fuerza, desplazamiento, velocidad y aceleración, las cuales requieren necesariamente un
conocimiento de dirección, es se establece un sistema de coordenadas, junto con las direcciones
positivas y negativas en este ejemplo observamos que la variación se realiza respecto al eje x.
La velocidad instantánea es P es 𝑣 = 𝑑𝑥/𝑑𝑡. La aceleración instantánea en P es 𝑎 = 𝑑𝑣/𝑑𝑡 o 𝑎 =
𝑑2 𝑥/𝑑𝑡 2 . La fuerza que actúa es el peso, siendo su magnitud 𝑃 = 𝑚𝑔. Por la ley de Newton
tenemos:
𝑑𝑣 𝑑𝑣
𝑚 = 𝑚𝑔 𝑜 =𝑔
𝑑𝑡 𝑑𝑡
Puesto que la masa cae desde el reposo, vemos que 𝑣 = 0 cuando 𝑡 = 0, o en otras palabras 𝑣(0) =
0.
Nuestra formación matemática es el problema de valor inicial 𝑣(𝑡)
𝑑𝑣(𝑡)
= 𝑔 𝑣(0) = 0
𝑑𝑡
Aquí tenemos una ecuación de primer orden y su condición requerida. Otra manera de formular el
problema es escribir:
𝑑2 𝑥 𝑑2 𝑥
𝑚 2 = 𝑚𝑔 𝑜 =𝑔
𝑑𝑡 𝑑𝑡 2
En tal caso tenemos una ecuación de segundo orden en las variables 𝑥 𝑦 𝑡, y necesitamos dos
𝑑𝑥
condiciones para determinar 𝑥. Una de ellas es 𝑣 = 0 𝑜 𝑑𝑡 = 0 𝑒𝑛 𝑡 = 0. La segunda puede
obtenerse al notar que 𝑥 = 0 𝑒𝑛 𝑡 = 0 (puesto que conseguimos el origen de nuestro sistema de
coordenadas en A).
Su formulación matemática es:
𝑑2 𝑥 𝑑𝑥
= 𝑔 𝑥 = 0 𝑦 = 0 𝑒𝑛 𝑡 = 0
𝑑𝑡 2 𝑑𝑡
Cuando se establecen las ecuaciones diferenciales para describir algún fenómeno o ley, siempre las
acompañaremos de suficientes condiciones necesarias para la determinación de las arbitrarias en
solución general.
Solución.
𝑑𝑣
Empezamos con = 𝑔 (separación de las variables) obtenemos por integración 𝑣 = 𝑔𝑡 + 𝑐1 .
𝑑𝑡
𝑑𝑥
Puesto que 𝑣 = 0 cuando 𝑡 = 0, 𝑐 1 = 0, 𝑜 𝑣 = 𝑔𝑡, 𝑒𝑠𝑡𝑜 𝑒𝑠, 𝑑𝑡 = 𝑔𝑡.
1
Otra integración produce de la anterior ecuación 𝑥 = 𝑔𝑡 2 + 𝑐2 . Puesto que 𝑥 = 0 𝑒𝑛 𝑡 = 0.
2
1 𝑑2 𝑥
𝑐2 = 0 por tanto 𝑥 = 2 𝑔𝑡 2 . Podríamos haber llegado al mismo resultado al empezar con 𝑚 𝑑𝑡 2 =
𝑑2 𝑥
𝑚𝑔 o 𝑑𝑡 2
=𝑔
El signo nos indica que el objetivo se está movido en la en la dirección positiva, esto es, hacia abajo.
Se debería tener en cuenta que si hubiéramos tomado la dirección positiva hacia arriba la ecuación
𝑑𝑣 𝑑𝑣 𝑑2 𝑥
diferencial hubiera sido 𝑚 ( ) = −𝑚𝑔, 𝑒𝑠𝑡𝑜 𝑒𝑠, = −𝑔 𝑜 = −𝑔
𝑑𝑡 𝑑𝑡 𝑑𝑡 2
Esto conduciría a resultados equivalentes a los obtenidos.
Problema propuesto situación 2
Considere un circuito eléctrico consistente en una fuente de voltaje 𝐸 (batería o generador),
una resistencia 𝑅, y un inductor 𝐿, conectados en serie como se muestra en la figura
Fuente. Propia
Adoptamos una convención: la corriente fluye del lado positivo de la batería o generador
a través del circuito hacia el lado negativo.
Por la ley de Kirchhoff, la fuerza electromotriz, 𝐸, es igual a la caída de voltaje a través del
𝑑𝐼
inductor, 𝐿 𝑑𝑡 , más la caída de voltaje a través de la resistencia, 𝑅𝐼, tenemos como la
ecuación diferencial requerida para el circuito:
𝑑𝐼
𝐸=𝐿
𝑑𝑡
+ 𝑅𝐼
Usted debe realizar las siguientes etapas del problema:
1. Explicar las condiciones asociadas que describen el circuito.
2. Formular la ecuación diferencial referente al problema.
3. Resolver la ecuación diferencial.
Un generador con una 𝑓𝑒𝑚 se conecta en serie con una resistencia y un inductor. Si el interruptor K
se cierra en tiempo 𝑡 = 0, establezca una ecuación diferencial para la corriente y determine la
corriente en tiempo 𝑡.
Llamado a 𝐼 la corriente o intensidad de corriente que fluye según el primer circuito descrito,
dI
tenemos, por ley de Kirchhoff, 𝐸 = 𝐿𝑥 = dt + 𝑅𝐼
Puesto que el interruptor se cierra en 𝑡 = 0, debemos tener 𝐼 = 0 𝑒𝑛 𝑡 = 0.
Solución:
𝑑𝐼
La ecuación diferencial anterior 𝐸 = 𝑅𝐼 + 𝐿 𝑑𝑡 es una ecuación de primer orden lineal exacta,
𝑅
buscando un factor integrante obtenemos µ(𝑡) = 𝑒 2 𝑡 . multiplicado por ese factor la ecuación, da
𝑅 𝑅
𝑅 𝑅 𝑅 𝑅 𝑡 𝑅 𝑡
𝑡 𝑡 𝑡 𝑑𝐼 𝑡 𝑑(𝐼𝑒 2 ) 𝑡 𝐸𝑒 2
𝐸𝑒 2 = 𝑅𝐼𝑒 2 + 𝐿𝑒 2 , es decir 𝐸𝑒 2 = integrando 𝐼𝑒 2 = +c
𝑑𝑡 𝑑𝑡 10
Puesto que 𝐼 = 0 𝑒𝑛 𝑡 = 0, podemos con estas condiciones obtener la constante c.
