Diferencia entre revisiones de «Modelos epidemiológicos. Grupo C14»
(→. CASO DE POBLACIÓN DE RIESGO CONSTANTE) |
(→. CASO DE POBLACIÓN DE RIESGO CONSTANTE) |
||
| Línea 44: | Línea 44: | ||
Conocido el dato de población de riesgo y siendo éste constante,el problema queda simplificado a una única ecuación diferencial lineal y de coeficientes constantes.<br/> | Conocido el dato de población de riesgo y siendo éste constante,el problema queda simplificado a una única ecuación diferencial lineal y de coeficientes constantes.<br/> | ||
Para resolver dicha ecuación, se han utilizado los métodos de Euler y trapecio, dados los siguientes datos iniciales:<br/> | Para resolver dicha ecuación, se han utilizado los métodos de Euler y trapecio, dados los siguientes datos iniciales:<br/> | ||
| + | |||
<math> \left\{\begin{matrix} | <math> \left\{\begin{matrix} | ||
I_{0}=2000\\ | I_{0}=2000\\ | ||
| Línea 50: | Línea 51: | ||
c=0.01 | c=0.01 | ||
\end{matrix}\right.</math><br/> | \end{matrix}\right.</math><br/> | ||
| + | |||
Éste es el programa que resuelve la ecuación mediante el método de Euler:<br/> | Éste es el programa que resuelve la ecuación mediante el método de Euler:<br/> | ||
| + | |||
Y éste el que lo resuelve por el método del trapecio:<br/> | Y éste el que lo resuelve por el método del trapecio:<br/> | ||
===. <math> S=0 </math> === | ===. <math> S=0 </math> === | ||
Debido a que la ecuación lineal es de coeficientes constantes, es sencillo obtener la solución de manera analítica, resultando:<br/> | Debido a que la ecuación lineal es de coeficientes constantes, es sencillo obtener la solución de manera analítica, resultando:<br/> | ||
| + | |||
<math> I=2000e^{-0.31t}</math> | <math> I=2000e^{-0.31t}</math> | ||
===. <math> S=100 </math>=== | ===. <math> S=100 </math>=== | ||
Revisión del 17:11 4 mar 2015
| Trabajo realizado por estudiantes | |
|---|---|
| Título | Modelos epidemiológicos. Grupo C14 |
| Asignatura | Ecuaciones Diferenciales |
| Curso | Curso 2014-15 |
| Autores | Palacios Pintor, Pedro
Pontiveros Bermejo, Diego Reinoso Muñoz, Cristina Rojas Arranz, Almudena Torre Prado, Yago de la Vidal Sánchez, Nieves |
| Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura | |
En esta página discutiremos y analizaremos un modelo epidemiológico mediante métodos de cálculo numérico.
El modelo se basa en las siguientes hipótesis:
- El número de personas infectadas sólo se altera por el fallecimiento o cura de éstas, y se ve afectada por el número de contagios entre la población de riesgo.
- La tasa de individuos que pasan de la población de riesgo a estar infectados es proporcional a la interacción entre el número de individuos entre ambas clases.
Contenido
1 . INTERPRETACIÓN DEL MODELO
Matemáticamente el modelo viene definido por las siguientes ecuaciones:
[math] \left\{\begin{matrix} \frac{dS}{dt}=-aSI\\ \frac{dI}{dt}=aSI-bI-cI \end{matrix}\right. [/math]
Donde:
-La función [math] S [/math] indica la cantidad de personas en riesgo de ser infectadas e [math] I [/math] el número de individuos que ya han contraído la enfermedad.
- El valor de [math] a [/math] es un coeficiente relacionado con la probabilidad de contagio cuando se juntan individuos de la población de riesgo con la de infectados. El coeficiente [math] b [/math] es el porcentaje de infectados que sanan en una unidad de tiempo, en este caso días, y el [math] c [/math] el porcentaje de infectados que fallece cada día.
2 . CASO DE POBLACIÓN DE RIESGO CONSTANTE
Conocido el dato de población de riesgo y siendo éste constante,el problema queda simplificado a una única ecuación diferencial lineal y de coeficientes constantes.
Para resolver dicha ecuación, se han utilizado los métodos de Euler y trapecio, dados los siguientes datos iniciales:
[math] \left\{\begin{matrix}
I_{0}=2000\\
a=0.003\\
b=0.3\\
c=0.01
\end{matrix}\right.[/math]
Éste es el programa que resuelve la ecuación mediante el método de Euler:
Y éste el que lo resuelve por el método del trapecio:
2.1 . [math] S=0 [/math]
Debido a que la ecuación lineal es de coeficientes constantes, es sencillo obtener la solución de manera analítica, resultando:
[math] I=2000e^{-0.31t}[/math]