MODELIZACION DE DATOS EXPERIMENTALES EJERCICIOS -ECUACION LOGISTICA / SISTEMA PRESA - DEPREDADOR Alvaro Guerrero Bonilla Ejercicio 1 :

Para t(0) hay P toneladas de peces ; y su evolucion esta dada por: dP/dt= rP(1-P/k) ; a) Si en un instante t(0) P=10, y se comienza a pescar a razon de 5 Ton/mes, en los 3 primeros meses y a razon de 1 ton/mes el resto de meses de la temporada: Cual es la poblecion al cabo de 15 meses? Cual es el valor minimo de la poblacion? A partir del mes 3 , tiende la poblacion a estacionarse? Hallar dicho valor de equilibrio Solucion: creamos la funcion captura para estarla a la ecuacion original: Utilizando matlab: >> edit captura function y=captura(t) y=5*(t<=3)+1(t>3); Utilizando dfield7:

Damos condiciones inicales : Options->Keyboard Input ->

Para t=0 ; P(0)=10

1 DE15

MODELIZACION DE DATOS EXPERIMENTALES EJERCICIOS -ECUACION LOGISTICA / SISTEMA PRESA - DEPREDADOR Alvaro Guerrero Bonilla
Obtenemos entonces:

Averiguamos la poblacion en el mes 15 , al calcular la curva en el intervalo 0t15

Obtenemos la matriz de resultados: en matlab: >>M=[dfdata1.t;dfdata1.P]; >>M A partir de t cercano a 4 -> P=0 Por lo tanto la poblacion tiende a la extincion. calculo del valor minimo de la poblacion : [a,b]=min(dfdata2.P) Ans => a=-19,5428 ; b=169 => para hallar el mes : >>dfdata.t(b) = > ans = 6,2965 , maximo valor de t para el que habia entregado valores de P

2 DE15

MODELIZACION DE DATOS EXPERIMENTALES EJERCICIOS -ECUACION LOGISTICA / SISTEMA PRESA - DEPREDADOR Alvaro Guerrero Bonilla Ejercicio 2 :
Supongamos una condicion de captura constante : Capt2(t)=k : Hallar equlibrios Al variar k en el intervalo(o,5) , cambian los equilibrios? DE ser asi , para que valor de k se produce cambio en el numero de equilibrios? Que tipo de bifurcacion de equilibrio tiene lugar? Realizar retrato de fase de cada uno de los casos. Debemos realizar el retrato de fase para cada valor de H(t)=k a)H(t)=0

Para H(t)=0 Equilibrios P(1-P/10)=0 PE1=0 PE2=10

Nota : Las curvas se dibujan dando click en cada area del diagrama Para visualizar : Options->lineas de phase ; solution dirtection (both) Repetimos el mismo procedimiento midificando H(t) : H(t)=1 Para H(t)=1 Equilibrios P(1-P/10)-1=0 P-P^2/10-1=0 PE1= 8,872983346 PE2 1,127016654

3 DE15

MODELIZACION DE DATOS EXPERIMENTALES EJERCICIOS -ECUACION LOGISTICA / SISTEMA PRESA - DEPREDADOR Alvaro Guerrero Bonilla
H(t)=2 Para H(t)=2 Equilibrios P(1-P/10)-2=0 P-P^2/10-2=0 PE1= 7,236067977 PE2 2,763932023

H(t)=3 Para H(t)=3 Equilibrios P(1-P/10)-3=0 P-P^2/10-3=0 PE1= #NUM! PE2 #NUM! No tiene puntos de equilibrio El maximo valor de H(t) Para el que se presenta un PE corresponde a : 1-4/10*H(t)=0 => H(t)=10/4=2,5 H(t)=2,5 Para H(t)=3 Equilibrios P(1-P/10)-3=0 P-P^2/10-3=0 PE1= 5 PE2 5

Ejercicio 3 :
Sea la ecuacion: y''-*y'+y+y^3=0 a) Variando en el intervalo(-1,1) ; encontrar y clasificar los distintos tipos de equilibrio: x'=y ;=>x'=0=y y'=*y-x-y3 ;=> Para y=0 => x=0 4 DE15

MODELIZACION DE DATOS EXPERIMENTALES EJERCICIOS -ECUACION LOGISTICA / SISTEMA PRESA - DEPREDADOR Alvaro Guerrero Bonilla
Debemos verificar las estabilidades para los valores de entre (-1,1): Utilizando pplane7

Solutions-> keyboard input-> damos condiciones iniciales x(0)=0, Y(0)=0

Para obtener total claridad sobre las caracteristicas del equilibrio Solutions->Find an equilibrium point->

Confirmamos que (0,0) es una espiral atractiva=> estable

5 DE15

MODELIZACION DE DATOS EXPERIMENTALES EJERCICIOS -ECUACION LOGISTICA / SISTEMA PRESA - DEPREDADOR Alvaro Guerrero Bonilla

Dando click sobre la grafica logramos la represen tacion de las curvas Se observa que no hay orbitas periodicas

Repetimos los mismos pasos para =a= 0 Y con Keyboard input ingresamos condiciones iniciales (0,0) Ahora, Solutions->Find equilibrium point , encontramos PE y dibujamos la linealizacion:

; => (0,0)=Centro en equilibrio

Ahora para =a= 0 : Evaluando para a=1

6 DE15

MODELIZACION DE DATOS EXPERIMENTALES EJERCICIOS -ECUACION LOGISTICA / SISTEMA PRESA - DEPREDADOR Alvaro Guerrero Bonilla
hallamos el equilibrio

Confirmamos que (0,0) es una espiral repulsiva = inestable

Ejercicio 4 :
Sea x'=x(2-x)-axy ; a>0,5 y'=y(3-y)-2xy Resolver las ecuaciones para diferentes valores de a>0,5 Solucion: 1 - Buscamos los puntos de equilibrio , que corresponden a las condiciones x'=0 y Y'=0; en condicion independiente de a : a) Trabajamos con la ecuacion x'=x(2-x)-axy=0 => PE1 : x'=0 para (x=0,y=0) PE1=(0,0) PE2 : si y =0 =>x(2-x)=0 => x=2 => PE2=(2,0) Tenemos entonces los puntos de equilibrio(0,0) y (2,0) PE3 : Si x=y =>2y-y^2-ay^2=0; =>y^2(1+a)=2Y; =>y=2/(1+a) =>Linea de ceros =>H318 (x=2/(1+a), y=2/(1+a)) PE3=(2/(1+a), 2/(1+a) b)Trabajando con la ecuacion : y'=y(3-y)-2xy=0 PE4: Cuando x=0 => y(3-y)=0 : y=3 => PE4=(0, 3) PE5: Cuando y=0 => X=0 = PE1 PE6: Si y=x=>3Y-Y^2-2Y^2=0; =>3Y^2=3Y =>Y=1 ; (X=1, Y=1) PE6=(1,1) Comprobamos los puntos de equilibrio utilizando pplane7: Variando a entre (0,5 , 1) :

Iniciamos para a=0,5

7 DE15

MODELIZACION DE DATOS EXPERIMENTALES EJERCICIOS -ECUACION LOGISTICA / SISTEMA PRESA - DEPREDADOR Alvaro Guerrero Bonilla
(0,3) Punto silla

(2,0)=Punto atractivo = estable

(0,0)= Punto repulsivo = Inestable

Ahora para a=1 ; (0,3)=Punto atractivo = estable

(1,1)=Punto silla

(2,0)=Punto atractivo = estable

(0,0)= Punto repulsivo = Inestable

Comprobamos los 4 puntos encontrados numerica . mente para valores superiores a a=0,65 Aparece un foco atractivo adicional generandose una bifurcacion en el modelo

8 DE15

MODELIZACION DE DATOS EXPERIMENTALES EJERCICIOS -ECUACION LOGISTICA / SISTEMA PRESA - DEPREDADOR Alvaro Guerrero Bonilla
EJERCICIO 5: Demostracion cualitativa de los puntos de equilibrio en el sistema presa - depredador , y su apliacion para determinacion de las funciones sobre una base de mediciones . SISTEMAS PRESA DEPREDADOR Objetivos : Mediante la utilizacin de mtodos numricos y cualitativos, encontrar los puntos de equilibrio, anlisis de estabilidad alrededor de dichos puntos, y demostracin que el sistema tiene soluciones periodicas. EJ.: Sean X(t) y Y(t) , las poblaciones de presas y depredadores respectivamente. Entonces: Para las presas: 1. La velocidad con que varia la poblacin de presas x es proporcional a su poblacin existente en un momento t, =>dx/dt=a1x(t) 2. La velocidad con que varia la poblacion de presas x es proporcional al numero de encuentros con los depredadores y ,causando su disminucin =>dx/dt= - a2x(t)y(t) Entonces: dX(t)/dt =a1X(t)-a2X(t)Y(t) De forma similar , para los depredadores su poblacin variara con base en: 3. El numero de depredadores existentes en un momento t , disminuidos por su tasa de mortalidad =>dy/dt= -b1y(t) 4. El numero de encuentros depredadores presa , ser causa de incremento en la poblacin :=> dy/dt=b2X(t)Y(t) , Combiando ambos efectos : dy/dt= -b1y(t)+b2x(t)y(t) ANALISIS CUALITATIVO DEL MODELO: Encontrar puntos de equilibrio: Los puntos de equilibrio se dan en la condicin : X(t)=0 , Y(t)=0. => X(t)= dX(t)/dt =a1X(t)-a2X(t)Y(t) =0 ; x(a1-a2Y)=0 => Y(t)= dY(t)/dt= -b1Y(t)+ b2X(t)Y(t)=0 ; Y(-b1+b2X)=0 ; Soluciones : (x(t),Y(t)) : PE1(0,0) ; PE2(b1/b2,a1/a2) Para estas soluciones , dados unos valores iniciales para presas y depredadores, sus nmeros permanecern constantes, por lo tanto estas soluciones son llamadas de punto de equilibrio o estacionarias.

En el plano x,y :

Punto de estabilidad a1/a2 b1/b2 9 DE15

MODELIZACION DE DATOS EXPERIMENTALES EJERCICIOS -ECUACION LOGISTICA / SISTEMA PRESA - DEPREDADOR Alvaro Guerrero Bonilla
Para la obtencin de la solucin general , dividimos la ecuacin 1 entre la 2 , obteniendo: dX(t)/dt =a1X(t)-a2X(t)Y(t) (1) dy/dt= -b1y(t)+b2x(t)y(t) (2) => dx/dy=( a1x-a2xy)/(-b1y+b2xy) => dx/dy= x(a1-a2y)/y(-b1+b2x) => (-b1+b2x)dx/x=(a1-a2y)dy/y => integrando : -b1logx+b2x=a1logy-a2y +K ; (K = constante arbitraria) b2x+a2y-(b1logx+a1logy)=K , => b2x+a2y-log(x^(b1)*y^(a1))=K K depende de las condiciones iniciales , => K=b2*x(0)+ a2*y(0) log(x(0)^b1*y(0)^a1) LINEALIZACION DEL SISTEMA : Definimos las funciones: F(x,y) = a1x-a2xy ; G(x,y)=-b1y+b2xy Calculamos las derivadas parciales : dF/dX= a1-a2y ; dF/dY= -a2x ; dG/dX=b2y ; dG/dY=-b1 + b2x , => Representados en la matriz jacobiana : dF/dX dF/dY dG/dX dG/dY Ahora , sustituimos los valores (x,y) PE1 y PE2 encontrados , en la matriz Jacobiana para hallar las soluciones : Para PE1 (x=0, y=0) => J(0,0)= a1 0 0 -b1 La solucin es de la forma: {a1,-b1} =>{(1,0),(0,1)} Para las funciones (t) tenemos : X(t)= c1*exp(a1t) ; Y(t)= c2*exp(-b1t). Para el punto PE2 (b1/b2, a1/a2) => J(0,0) = 0 - a2b1/b2 b2a1/b2 0

TABLA Indice captura conejos y linces compaa Hudson Bay 1900-1920 Ao Presa Depredador 1 30 4 2 47,2 6,1 3 70,2 9,8 4 77,4 35,2 5 36,3 59,4 6 20,6 41,7 7 18,1 19 8 21,4 13

10 DE15

MODELIZACION DE DATOS EXPERIMENTALES EJERCICIOS -ECUACION LOGISTICA / SISTEMA PRESA - DEPREDADOR Alvaro Guerrero Bonilla
9 10 11 12 13 14 15 16 17 18 19 20 21 22 25,4 27,1 40,3 57 76,6 52,3 19,5 11,2 7,6 14,6 16,2 24,7 8,3 9,1 7,4 8 12,3 19,5 45,7 51,1 29,7 15,8 9,7 10,1 8,6

1. Valores iniciales: X(0)=30 ; Y(0)= 4 , (datos de la tabla) 2. Valores que satisfagan las ecuaciones : X(t)=b1/b2 , Y(t)= a1/a2 Tomamos los valores de las medias( tabla) => b1/b2=34,08 ; a1/a2=20,17 Vamos a obtener los valores bases a1 y b1: a1 : Buscamos 2 aos consecutivos donde tengamos un alto incremento en el numero de conejos , con una baja poblacin de linces : Ao 11 =27,1 conejos Ao 12 = 40,3 conejos De acuerdo con la ecuacin X(t)=X(t-1)*exp(a1) => a1 = Ln(X(t)/X(t-1)) => 76,6= 57 * expa1 => a1 = Ln(76,6/57) =0,3968 B1. Buscamos un punto con poblacin bajo de conejos , que represente un rpido descenso den los linces : Ao 7: 19 linces Ao 6: 41,7 linces Y(t)=Y(t-1)*exp(-b1) => b1= -Ln(Y(t)/Y(t-1)) ; b1=-Ln(19/41,7) =0,7860 ; =>a1/a2 =20,17=> a2= a1/20,17 =0,3968/20,17 =0,0197= a2 =>b1/b2=34,08 =>b2=b1/34,08=0,786/34,08=0,023=b2 Tendremos entonces la siguiente representacin de las ecuaciones : X(t)= a1x(t)-a2x(t)y(t) = 0,4x-0,02xy ; x(0)=30 Y(t)= -b1y(t)+b2x(t)y(t)= -0,8y+0,023xy ; y(0)=4

EJERCICIOS CAPITULO 4 - MODELOS DEPREDADOR- PRESA MODELOS Y SISTEMAS DINAMICOS Ejercicio 1 :

11 DE15

MODELIZACION DE DATOS EXPERIMENTALES EJERCICIOS -ECUACION LOGISTICA / SISTEMA PRESA - DEPREDADOR Alvaro Guerrero Bonilla
Demostrar que enel punto de equilibrio no trivial , el centro existente para =0, se transforma en foco inestable si 0<<1. Dibujar diagrama de fases a) Identificamos los puntos de equilibrio , que corresponderan a los valores x,y para los cuales dx/dt =0 y dy/dt=0: a.1)En la ecuacion dx/dt=x(1-y-x) => tenemos las siguientes condiciones : x=0 -> y=0 : PE1(0,0); Para y=x => x(1-x)=0 ; x=1/ => La recta x=1/(1+) es una linea de ceros = PE2 Si y=0; => x(1-X)=0 =>x=1/ => PE3(1/,0) a.2)En la ecuacion dy/dt=-y(1-x-y), obtenemos: Para y=0 , entonces x=0 ; PE1(0,0); Para y=x=> -y(1-y-y)=0; =>y=y^2(1+)=> y=1/(1+); es una linea de ceros =PE2 Para x=0 => -y(1-y)=0 =>-y+y^2=0 ; y=1/ =>P4=(0,1/) 1/ 1/(1+)=(xe2 =(xe2,ye4) , ye2) Evolucion del nodo 0 #DIV/0! 1 Nodo silla 0 2 0,66666667 Nodo silla 0 1 0,5 Punto de equilibrio

0 0,5 1

PE1

12 DE15

MODELIZACION DE DATOS EXPERIMENTALES EJERCICIOS -ECUACION LOGISTICA / SISTEMA PRESA - DEPREDADOR Alvaro Guerrero Bonilla Ejercicio 2 :
Un modelo depredador - presa para describir la iinteraccion de herbivorso H y plancton P es: dP/dt=r*P(K-P-(BH/(C+P))) ; dH/dt= DH(P/(C+P)-AH) donde : r,K,A,B,C,H son constantes positivas 2.1) describa brevemente las supocisones ecologicas del modelo La ecuacion dp/dt=rP(K-P-(BH/C+P)) , se ajusta al modelo dx/xt = x(r(1-x/K)-ky/(x+D)) , lo que garantiza que la ecuacion de ajuste al modelo Lotka - Volterra : a-Es dependiente del numero de depredadores ; bCumple con la condicion de ser =0 , para la condicion x=0 quie hace la funcion dp/dt=0. c- Aparenta incluir algun grado de saturacion para el cntrol establecido por los depredadores. En la ecuacion dH/dt , se evidencia que el crecimiento de la poblacion de depredadores se condiciona por el numero de las presas. Las condiciones indicadas no se alteran al ser r,K,A,B,C,H , constantes positivas. 2.2)Demostrar que existe un estado de equilibrio positivo (po,ho), para cada a>0 , k>0 Solucion: Primero buscamos las condiciones para las cuales se cumple: dp/dr=0 = p(k-p-h/(1+p)) =0 ,y ,dh/dr=dh(p/(1+p)-ah) Trabajando con la primera ecuacion: p(k-p-h/(1+p)) =0 Solucion 1: p=0 ; de esta solucion se deriva Un punto Q1=(0,0) Solucion 2 :Para h=p => pk-p^2-p^2/(1+p)=0; => (pk-p^2)(1+p)=p^2 => (k-p)(1+p)=p; => La curva (p,(k-p)(1+p)) representa una curva de ceros Sobre esta curva se encontraran los puntos(P1 y P2) correspondientes a las soluciones para: p=[(k-2)+/RAIZ((2-k)^2+4k)]/2 De las cuales P1 debe corresponder al cuadrante 1 y P2 al cuadrante 3 Trabajando con la segunda ecuacion: dh/dr = dh(p/(1+p)-ah) Solucion 4 : h=0 => Corresponbde a la solucion 1 Solucion 5: p=h=> dh(h/(1+h) - ah)=0; => dh^2/(1+h)-adh^2=0; => dh^2/((1+h)=adh^2 1/(1+h)=a; => 1=a+ah=> ah=-a => h=-1 => La recta (-1,h) representa otro curva de ceros Por lo tanto Se puede identificar qu el punto Q1(-1,0) es un punto de equilibro Comprobacion cuantitativa: Sean k=1, a=1, d=1 x y x1 -2,5 -5,25 -1 -2,2 -3,84 -1 -2 -3 -1 -1,5 -1,25 -1 -1 0 -1 -0,5 0,75 -1 0 1 -1 0,5 0,75 -1 0,618033989 0,61803399 -1 dp/dr=p(k-p-h/(1+p))=0 0 [Link] 13 DE15

MODELIZACION DE DATOS EXPERIMENTALES EJERCICIOS -ECUACION LOGISTICA / SISTEMA PRESA - DEPREDADOR Alvaro Guerrero Bonilla
1 0 -1 1,5 -1,25 -1 2 -3 -1 2,2 -3,84 -1 2,5 -5,25 -1 Soluciones a la ecuacion p=[(k-2)+/-RAIZ((2-k)^2+4k)]/2: k P1 p2 1 0,61803399 -1,61803399 => P1=(0.61803399, 0.61803399) es una solucion en el cuadrante 1 hacemos variando k a 0,5 x y x1 -2,5 -4,5 -1 -2,2 -3,24 -1 -2 -2,5 -1 -1,5 -1 -1 -1 0 -1 -0,5 0,5 -1 0 0,5 -1 0,280776406 0,28077641 -1 dp/dr=p(k-p-h/(1+p))=0 0 [Link] 0,5 0 -1 1 -1 -1 1,5 -2,5 -1 2 -4,5 -1 2,2 -5,44 -1 2,5 -7 -1 Soluciones a la ecuacion p=[(k-2)+/-RAIZ((2-k)^2+4k)]/2: k P1 p2 0,5 0,28077641 -1,78077641 => P1=(0.28077641,0.28077641) es una solucion en el cuadrante 1 2.3 -Dibujar curvas de pendientes verticales y horizontales. Observar los cambios cualitativos que se producen con la variacion de k Condiciones: k=1, a=1,d=1

Al variar el valor k desde 0 hasta 1 , observamos como pasamos de una condicion sin soluciones en el cuadrante 1(cuadrante de maximo interes), hasta la mayor area bajo la curva para el rango en estudio.

14 DE15

MODELIZACION DE DATOS EXPERIMENTALES EJERCICIOS -ECUACION LOGISTICA / SISTEMA PRESA - DEPREDADOR Alvaro Guerrero Bonilla
2.4 - Estrategia para determinar la estabilidad del estado de equilibrio : Se deben considerar las trayectorias alrededor de P1. Evaluando la matriza jacobiana de la funcion en le punto P1 => Si los autovalores de la matriz tienen parte no nula y ademas tiene valor negativo , se puede asegurar la estabilidad del punto Si las orbitas espirales se acercan al p1 = Estable , si se alejan = inestable.

15 DE15