domingo, 26 de abril de 2015

ecuaciones diferenciales



Din�mica de poblaciones



Aqu� est�n algunas preguntas naturales relacionadas con los problemas de la poblaci�n:
  • �Cu�l ser� la poblaci�n de cierto pa�s dentro de diez a�os?
  • �C�mo estamos protegiendo los recursos contra la extinci�n?
Puede ser dicho m�s sobre el problema pero, en esta peque�a revisi�n que no los discutiremos detalladamente. Para ilustrar el uso de ecuaciones diferenciales con respecto a este problema consideramos el modelo matem�tico m�s f�cil ofrecido para gobernar la din�mica de la poblaci�n de cierta especie. Com�nmente se llama el modelo exponencial , es decir, el �ndice del cambio de la poblaci�n es proporcional a la poblaci�n existente. Es decir si (t) mide a poblaci�n, tenemos
displaymath58,
donde  k es el �ndice de crecimiento. Es bastante f�cil ver que si > 0, tenemos crecimiento, y si < 0, disminuir� la poblaci�n. �sta es una ecuaci�n lineal con soluciona
displaymath66,
donde  tex2html_wrap_inline68  est� la poblaci�n inicial, es decir.  tex2html_wrap_inline70  . Por lo tanto, concluimos lo siguiente:
  • si > 0, la poblaci�n crece y contin�a ampli�ndose al infinito, es decir,displaymath74
  • si < 0,  la poblaci�n disminuye y tiende a 0, o sea, se extingue.
Claramente, el primer caso, > 0, no es real, con lo que el modelo no es bueno. Deben considerarse ciertas limitaciones ambientales. La complicaci�n es que el crecimiento de la poblaci�n es limitado eventualmente por un cierto factor, generalmente por los recursos esenciales. Cuando una poblaci�n est� lejos de sus l�mites del crecimiento puede crecer exponencialmente. Sin embargo, al acercarnos al l�mite, el modelo no representa la realidad.
Otro modelo fue propuesto para remediar este defecto del modelo exponencial. Se llama el modelo log�stico (tambi�n llamado modelo de Verhulst-Verhulst-Pearl ). La ecuaci�n diferencial para este modelo es
displaymath80,
donde tenemos un tama�o M limitador para la poblaci�n (tambi�n llamada la capacidad de carga). Claramente, cuando es peque�o comparado a M , la ecuaci�n reduce la exponencial. Para solucionar esta ecuaci�n reconocemos una ecuaci�n no lineal que sea separable. Las soluciones constantes son = 0 y . Las soluciones no-constantes pueden ser obtenidas separando las variables
displaymath92,
e integrando
displaymath94
As�
displaymath96,
lo cu�l da
displaymath98
Las manipulaciones algebraicas f�ciles dan
displaymath100
donde nos queda una constante C. Solucionando para , conseguimos
displaymath106
Si consideramos la condici�n inicial  tex2html_wrap_inline70  (se asume que  tex2html_wrap_inline68  no es igual a 0 o M), conseguimos
displaymath116,
lo cu�l, substituido una vez en la expresi�n para (t) y simplificado, nos da
displaymath120
Es f�cil ver que
displaymath122






Din�mica de poblaciones: Ejemplo 2


Ejemplo: Sea (t)  la poblaci�n de cierta especie animal. Asumamos que (t) satisface la ecuaci�n log�stica del crecimiento
displaymath23

1.
�Es la ecuaci�n diferencial separable?
2.
�Es la ecuaci�n diferencial aut�noma?
3.
�Es la ecuaci�n diferencial lineal?
4.
Sin solucionar la ecuaci�n diferencial, d� un bosquejo del gr�fico de P (t).
5.
�Cu�l es el comportamiento a largo plazo de la poblaci�n P(t)?
6.
Demuestre que la soluci�n es de la formadisplaymath27
Encuentre .
Pista: Utilice la condici�n inicial y el resultado de 5.
7.
�D�nde tiene la soluci�n un punto de la inflexi�n?
Pista: Esto se puede hacer sin usar la respuesta de 6.
8.
�Se puede destacar algo sobre el �ndice de crecimiento de la poblaci�n (t) en el punto de la inflexi�n (encontrado en 7)?

Respuestas:

1.
Esta ecuaci�n diferencial es aut�noma, es decir la variable t falta. Por lo tanto, esta ecuaci�n es de hecho separable.
2.
�La respuesta es s�! (v�ase 1.) Cada ecuaci�n diferencial aut�noma es separable.
3.
La ecuaci�n no es lineal debido a la presencia de  tex2html_wrap_inline92  .
4.
El gr�fico de la soluci�n (t) es como sigue:

 
5.
Claramente, debido a la condici�n inicial (v�ase el gr�fico de la soluci�n abajo), tenemosdisplaymath96,
200 que son la capacidad de carga.
6.
Solucionemos esta ecuaci�n (utilice la t�cnica para solucionar ecuaciones separables). Primero, buscamos las soluciones constantes (los puntos del equilibrio o puntos cr�ticos). Tenemosdisplaymath98
Entonces, las soluciones no-constantes pueden ser generadas separando las variables
displaymath100,
e integrando
displaymath102.
conseguimos
displaymath104,
o cu�l da
displaymath106
Por lo tanto, tenemos
displaymath108.
Las manipulaciones algebraicas f�ciles dan
displaymath110,
donde  tex2html_wrap_inline112  .
Por lo tanto, todas las soluciones son
displaymath114
donde C es un par�metro.

Observaci�n: Podemos reescribir las soluciones no-constantes como
displaymath118,
donde est�n dos par�metros a y B. Si utilizamos las condiciones
displaymath124
podremos conseguir la soluci�n deseada. De hecho, tenemos
displaymath126
As�,  tex2html_wrap_inline128  y por lo tanto  tex2html_wrap_inline130  .
7.
Del gr�fico (v�ase 4.), el gr�fico de P (t) tiene un punto de la inflexi�n entre 0 y 200. Encontr�moslo. Puesto que debemos tener ' ' (t) = 0, conseguimosdisplaymath136,
donde hemos utilizado que   tex2html_wrap_inline138  . Puesto que nuestra soluci�n no es una de las dos soluciones constantes nos dejan solamente con la ecuaci�n
displaymath140
Esto se simplifica
displaymath142,
lo cu�l da = 100 (a medio camino entre 0 y 200).

8.
El �ndice de crecimiento en = 100 es m�ximo.

 

No hay comentarios:

Publicar un comentario