Simulaciones de Ecuaciones Diferenciales de Primer Orden No Lineales
Simulaciones de Ecuaciones Diferenciales de Primer Orden No Lineales
Simulaciones de Ecuaciones Diferenciales de Primer Orden No Lineales
TEMA:
MODELADO CON ECUACIONES DIFERENCIALES
DE PRIMER ORDEN NO LINEALES MEDIANTE
MATHLAB - SIMULINK
INTEGRANTES:
BRYAN CHIMBORAZO
KEVIN GUADALUPE
MISHELLE VASCO
SEBASTIÁN VELÁSQUEZ
CURSO:
2do “A”
CARRERA:
“TELECOMUNICACIONES”
MATERIA
ECUACIONES DIFERENCIALES
PROFESOR:
Mgs. GERMAN MORENO
RIOBAMBA – ECUADOR
Contenido
Introducción……………………………………………………………………………...1
1
MODELACION MATEMATICA.
En general se entiende por “modelización matemática” el proceso por el cual se imita la realidad
en términos matemáticos. El objetivo evidentemente es bien explicar o comprender los
fenómenos naturales, bien encontrar respuestas a problemas técnicos o científicos. Toda
modelización lleva consigo un proceso de “idealización”. La realidad suele ser compleja y los
problemas reales habitualmente dependen de multitud de parámetros o variables, al mismo
tiempo que suelen estar inter-relacionados con otros procesos. El diseño de un modelo
matemático lleva aparejada la simplificación de muchos aspectos del problema real.
2
Pasos para modelar matemáticamente.
Identificación. Se trata de clarificar las preguntas que se intentan responder con el
modelo, formular el problema en palabras, documentar los datos relevantes e
identificar el mecanismo subyacente al problema real.
Suposiciones. El problema debe ser analizado para decidir los factores del mismo que
son importantes y aquellos que pueden ser ignorados. Con todo ello deben hacerse
suposiciones (o idealizaciones) lo más realistas posible.
Construcción. En este paso se “construye” el modelo, es decir, se traduce al lenguaje
matemático el problema (junto con las suposiciones anteriormente realizadas)
obteniéndose un conjunto de ecuaciones (o inecuaciones) después de haber
identificado las variables que deben intervenir en las mismas.
Análisis o Resolución. Se trata de la resolución del problema. Las soluciones consistirán
en general en funciones por medio de las cuales la o las variables dependientes se
expresarán en términos de la o las variables independientes. Por otro lado, se obtendra
información acerca de los parámetros que intervienen en el modelo.
Interpretación. En este paso, la solución matemática debe ser comparada con la
realidad para observar si se ajusta a lo conocido acerca del problema real. Se trata, en
definitiva, de interrumpir el proceso si se obtiene soluciones carentes de sentido real.
Validación. Una vez interpretada la solución, se comprueba numéricamente que
concuerda con los datos disponibles sobre el problema.
Implementación. Finalmente, se usa el modelo para describir el problema, se pueden
por tanto realizar predicciones sobre los valores de las variables. Es necesario prestar
atención al rango de validez del modelo.
3
Diferencias de entrada y salida.
Un gráfico muestra el conjunto de soluciones para una ecuación dada. En el caso de
ecuaciones lineales, el gráfico siempre será una línea. Por el contrario, una ecuación no
lineal puede parecerse a una parábola si es de grado 2, una forma de x curvada si es de
grado 3, o cualquier otro tipo de curva. Mientras que las ecuaciones lineales son siempre
rectas, las ecuaciones no lineales a menudo presentan curvas.
Excepciones.
Las relaciones lineales se pueden explicar mejor mediante ecuaciones lineales, donde el
aumento en una variable causa directamente el aumento o disminución de la otra. Por
ejemplo, el número de galletas que comes en un día podría tener un impacto directo en tu
peso como se ilustra mediante una ecuación lineal. Sin embargo, si se analiza la
distribución de células en mitosis, una ecuación exponencial no lineal encajaría mejor con
los datos.
EJEMPLOS Y SIMULACIONES.
1. Suponga que un estudiante es portador del virus de la gripe y regresa a un campus
aislado de 1 000 estudiantes. Si se supone que la razón con que se propaga el virus
no sólo a la cantidad x de estudiantes infectados sino también a la cantidad de
estudiantes no infectados, determine la cantidad de estudiantes infectados después
de 6 días si además se observa que después de cuatro días x(4) 50.
4
Simulación:
Grafica Solución:
5
2. La población P(t) de un suburbio de una gran ciudad en un instante cualquiera se rige
por:
En donde t se mide en meses. ¿Cuál es el valor límite de la población? ¿En qué momento será la
población igual a la mitad de su valor límite?
6
Simulación:
Grafica solución:
7
3.-Control de nivel del tanque mediante un ingreso y una salida de agua
Datos:
𝑏 = 0.09
𝑄1 = 0.05
Análisis
𝑑
𝑄1 (𝑡) − 𝑄2 (𝑡) = 𝐴 ℎ(𝑡)
𝑑𝑡
𝑑
𝑄1 (𝑡) − 𝐴2 𝑣2 = 𝐴 ℎ(𝑡)
𝑑𝑡
Ecuación de Bernoulli
𝑃1 𝑣1 2 𝑃2 𝑣2 2
+ +ℎ = + + ℎ2
𝛿𝑔 2𝑔 𝛿𝑔 2𝑔
𝑣2 2 = 2𝑔ℎ
𝒗𝟐 = √2𝑔𝒉
𝑑
𝑸𝟏 − 𝑸𝟐 = 𝐴 𝒉
𝑑𝑡
𝑑𝒉
𝑸𝟏 − 𝐴2 √2𝑔𝒉 = 𝐴
𝑑𝑡
𝑑𝒉
𝑸𝟏 − 𝑏√𝒉 = 𝐴
𝑑𝑡
𝑑𝑡 𝑑𝒉
=
𝐴 𝑸𝟏 − 𝑏√𝒉
𝑢 = 𝑸𝟏 − 𝑏√𝒉
√𝒉 = 𝑢 − 𝑸𝟏
8
𝑏𝑑ℎ
𝑑𝑢 = −
2√𝒉
2√𝒉
𝑑ℎ = − 𝑑𝑢
𝑏
2(𝑢 − 𝑸𝟏 )
𝑑ℎ = − 𝑑𝑢
𝑏
𝑑𝑡 𝑑𝒉
=
𝐴 𝑸𝟏 − 𝑏√𝒉
𝑑𝑡 2𝑸𝟏 − 2𝑢
= 𝑑𝑢
𝐴 𝑏𝑢
1 2𝑸𝟏 𝑑𝑢 2
∫ 𝑑𝑡 = ∫ − ∫ 𝑑𝑢
𝐴 𝑏 𝑢 𝑏
1 2𝑸𝟏 2
𝑡= ln|𝑢| − 𝑢 + 𝐶
𝐴 𝑏 𝑏
1 2𝑸𝟏 2
𝑡= ln|𝑸𝟏 − 𝑏√𝒉| − (𝑸𝟏 − 𝑏√𝒉) + 𝐶
𝐴 𝑏 𝑏
Para t=0 y h=0
2𝑸𝟏 2
0= ln|𝑸𝟏 | − 𝑸𝟏 + 𝐶
𝑏 𝑏
2
𝐶= 𝑸 (1 − ln|𝑸𝟏 |)
𝑏 𝟏
1 2𝑄1 2 2
𝒕= ln|𝑄1 − 𝑏√𝒉| − (𝑄1 − 𝑏√𝒉) + 𝑄1 (1 − ln|𝑄1 |)
𝐴 𝑏 𝑏 𝑏
9
Simulación:
Grafica solución:
10
Libro Dennis G. Zill (Problema 1)
4.- La cantidad N(t)N(t) de supermercados del país que están usando sistemas de revisión
computarizados se describe por el problema con valores iniciales
𝑑𝑁
= 𝑁(1 − 0.0005𝑁)
𝑑𝑡
𝑁(0) = 1
a) Use el concepto de sistema de fase para predecir cuantos supermercados se espera
que adopten el nuevo procedimiento en un periodo de tiempo largo. A mano, dibuje
una curva solución del problema con valores iniciales dados.
b) Resuelva el problema con valores iniciales y después utilice un programa de
graficación para comprobar y trazar la curva solución del inciso a) ¿Cuantas compañías
se espera adopten la nueva tecnología cuando t=10?
1 𝐴 𝐵
= +
𝑁(1 − 0.0005𝑁) 𝑁 1 − 0.0005𝑁
1 = 𝐴(1 − 0.0005𝑁) + 𝐵𝑁
1 = 𝐴 + (𝐵 − 0.0005𝐴)𝑁
𝐴=1
𝐵 = 0.0005
1 1 0.0005
= +
𝑁(1 + 0.0005𝑁) 𝑁 1 − 0.0005𝑁
1 0.0005
∫ 𝑑𝑁 + ∫ 𝑑𝑁 = ∫ 𝑑𝑡
𝑁 1 − 0.0005𝑁
ln|𝑁| − ln|1 − 0.0005𝑁| = 𝑡 + 𝐶
𝑁
ln | |=𝑡+𝐶
1 − 0.0005𝑁
𝑁 = 𝐶1 𝑒 𝑡 (1 − 0.0005𝑁)
𝐶1 𝑒 𝑡
𝑁=
1 + 0.0005𝐶1 𝑒𝑡
𝑁(0) = 1
𝐶1 𝑒 𝑡
𝑁=
1 + 0.0005𝐶1 𝑒𝑡
11
𝐶1 𝑒 0
1=
1 + 0.0005𝐶1 𝑒 0
𝐶1 − 0.0005𝐶1 = 1
𝐶1 = 1.00050025
1.00050025𝑒 𝑡
𝑁=
1 + 0.00050025𝑒 𝑡
𝟐𝟎𝟎𝟎
𝑵=
𝟏𝟗𝟗𝟗𝒆𝒕 + 𝟏
Simulación:
12
Grafica solución:
5.- Un reactor transforma plutonio 239 en uranio 238 que es relativamente estable
para uso industrial. Después de 15 años se determina que el 0.0043 por ciento de la
cantidad inicial A0 de plutonio se ha desintegrado. Determina la semivida1 de este
isótopo si la rapidez de desintegración es proporcional a la cantidad restante.
𝑑𝑥
= 𝜆𝑥
𝑑𝑡
𝑥(0) = 𝐴
𝑥(𝑡) = 𝐴𝑒 𝜆𝑡
0.0043
𝑥(15) = (1 − )𝐴
100
99.9957
𝐴𝑒 15𝜆 = 𝐴
100
13
1 99.9957
𝜆= ln ( )
15 100
Finalmente, la semivida de este isótopo es el valor t0 para el que se cumple la
condición x(t0) = A/2, por lo que
1
ln(
99.9957
)𝑡 𝐴
𝐴𝑒 15 100 =
2
𝑡 = 241790
Simulación:
Grafica solución:
14
6.- Cuando se combina os sustancias, A y B, se forman un compuesto C. La reacción entre
ambos es tal que, por cada gramo de A es usa 4 gramos de B. Se observa que a los 10 min se
han formado 30 gramos del producto C. Calcule la cantidad de C en función del tiempo si la
velocidad de la reacción es proporcional a las cantidades de Ay B que quedan y al principio hay
50 gramos de A y 32 gramos de B. ¿Qué cantidad de compuesto C hay a los 15 min? Interprete
la solución cuando t es infinito.
𝑥
𝐴
5
4𝑥
𝐵
5
𝑥 4
50 − 𝑦 32 − 𝑥
5 5
𝑑𝑥 𝑥 4𝑥
= 𝑘 (50 − ) (32 − )
𝑑𝑡 5 5
𝑑𝑥
= 𝑘(250 − 𝑥)(40 − 𝑥)
𝑑𝑡
1
𝑑𝑥 = 𝑘𝑑𝑡
(250 − 𝑥)(40 − 𝑥)
1 𝐴 𝐵
= +
(250 − 𝑥)(40 − 𝑥) 250 − 𝑥 40 − 𝑥
1 = 𝐴(40 − 𝑥) + 𝐵(250 − 𝑥)
𝐴 = −1/210
𝐵 = 1/210
1 1
− 210 𝑑𝑥 + 210 𝑑𝑥 = 𝑘𝑑𝑡
250 − 𝑥 40 − 𝑥
1 1 1 1
− ∫ 𝑑𝑁 + ∫ 𝑑𝑁 = 𝑘 ∫ 𝑑𝑡
210 250 − 𝑥 210 40 − 𝑥
250 − 𝑥
ln | | = 210𝑘𝑡 + 𝐶
40 − 𝑥
250 − 𝑥
= 𝐶𝑜 𝑒 210𝑘𝑡
40 − 𝑥
15
Cuando t=0 y x=0
24
𝐶𝑜 =
4
210𝑘 = 0.1258
250 − 𝑥
= 𝐶𝑜 𝑒 210𝑘𝑡 𝐷𝑒𝑠𝑝𝑒𝑗𝑜 𝑥 𝑦 𝑟𝑒𝑒𝑚𝑝𝑙𝑎𝑧𝑜 𝑣𝑎𝑙𝑜𝑟𝑒𝑠
40 − 𝑥
1 − 𝑒 −0.1258𝑡
𝑥 = 1000 ∗
25 − 4𝑒 −0.1258𝑡
𝑥 4
50 − 𝑦 32 − 𝑥
5 5
40 4
50 − = 42𝑔 𝑑𝑒 𝐴 𝑦 32 − ∗ 40 = 0𝑔 𝑑𝑒 𝐵
5 5
Simulación:
16
Grafica solución:
Recomendaciones:
Utilizar las versiones adecuadas de Matlab y Simulink para el óptimo
funcionamiento de acuerdo a las especificaciones de cada ordenador.
Ubicar de manera correcta cada una de los bloques de manera que describa los
procesos de modelaje realizado en cada ejercicio.
Una vez creada la simulación al momento de visualizar mediante Scope se deberá
ajustar la gráfica de manera correcta en función del tiempo para así poder
visualizar de una mejor manera la forma sentido y tendencia del resultado grafico
Conclusiones:
En esta práctica se utilizó simulink ya que es muy útil para realizar este tipo de
prácticas con la cual podemos resolver ecuaciones diferenciales, al realizar esta
práctica nos apoyamos en unos videos
Al realizar la práctica pudimos darnos cuenta de que teníamos un error en el
momento de poner los signos en la simulación debido a que habíamos colocado -
;+ cuando en realidad los signos debía de ser +,- ya que debía de ser igual a la
ecuación que estábamos realizando en la simulación entonces ya al cambiar los
signos como en la ecuación vimos que ya lo que nos mostraba estaba bien
17
Referencias
e.c.h. (s.f.). Ecuaciones de primer orden. Obtenido de
http://jacobi.fis.ucm.es/pparanda/EDNpdf/ppED1.pdf
18