Master Industriales. Tareas 3 y 4.
Ejercicio 1 (PVI para EDO. 3 puntos) En muchos casos el problema
de crecimiento de población, y(t) es la población en el instante t, se trata
con la ecuación (I) y ′ = k y, para una cierta constante k.
Pero en ciertos casos los recursos son limitados y hay que utilizar la ecuación
logı́stica, que se utiliza en diversos modelos de crecimiento de poblaciones,
propagación de enfermedades epidémicas, impacto de publicidad o difusión en
redes sociales entre otras. Constituye un refinamiento del modelo (I) cuando
al cabo de un tiempo, aparece la competencia entre algunos miembros, o falta
de recursos que impide que la población crezca por encima de una cierto
valor.
Básicamente podemos decir que el crecimiento es proporcional al tamaño de
la población en cada instante y a su diferencia con el tamaño máximo que
puede alcanzar. Si M es el tamaño máximo que puede alcanzar la población,
la ecuación quedarı́a
(II) y ′ = k y (M − y).
Con los datos M = 10000 y partiendo de una población de 100 habitantes.
Se pide
1. Utilizar el simbólico para resolver la ecuación (I) en un periodo entre
el instante inicial y 25 meses, sabiendo que al cabo de 12 meses la
población es de 1200 habitantes. Calcula la constante k y el instante
en el que se quedan sin recursos. [1 punto]
2. Resolver en simbólico la ecuación (II) en un periodo desde el inicio
hasta el mes 5, sabiendo que k = 10−3 . [0.5 punto]
3. Resolver el apartado anterior en numérico por uno de los métodos dado
en clase, tomando los pasos de discretización h = 0.01. Representar las
soluciones simbólica y numérica en una única gráfica. [0.5]
4. A partir de la solución analı́tica ver en qué momento la población al-
canza los 9000 habitantes. [0.5 puntos]
5. Utilizando la solución simbólica ver en qué momento la población al-
canza los 9000 habitantes.[0.5 puntos]
1
Ejercicio 2 (1 puntos) Sea dado el problema de transporte
′′
(
−u (x) + u′ (x) − u(x) = x22+1 ∀ x ∈ (0, 10)
u(0) = 0, u(10) = 1
que modela el flujo laminar de un fluido mediante los procesos de difusión,
convección y absorbción.
La solución u(x) del PVC representa la concentración de una especie quı́mica
transportada por el fluido a través de una corriente. Las condiciones de con-
torno son de tipo Dirichlet no homogéneas y definen el nivel de concentración
en los extremos del intervalo.
Calcular la solución numérica del problema en el intervalo espacial [0, 10]
utilizando diferencias finitas centradas para el término convectivo y un paso
de discretización h = 0.1. Representar gráficamente. Calcular los puntos
donde se encuentran el máximo y el mı́nimo de la concentración y los valores
de estos.