El objetivo de este trabajo es introducir el concepto de caos, mediante el análisis cualitativo de sistemas dinámicos discretos, cuyos modelos deterministas están dados por ecuaciones en diferencias no lineales. La variable independiente discreta es el tiempo y al cabo de cada intervalo de igual longitud la variable dependiente describira el estado del sistema. El análisis cualitativo se realiza geométricamente, mediante representaciones gráficas, sin aplicar métodos formales de resolución de la ecuación, usando algunas veces métodos iterativos

Bressan, Juan Carlos

Ferrazzi, Ana

Funes

Capítulo 2. Propuestas para la enseñanza de las matemáticas !Comité Latinoamericano de Matemática Educativa A. C. !647  Resumen. El objetivo de este trabajo es introducir el concepto de caos, mediante el análisis cualitativo de sistemas dinámicos discretos, cuyos modelos deterministas están dados por ecuaciones en diferencias no lineales. La variable independiente discreta es el tiempo y al cabo de cada intervalo de igual longitud la variable dependiente describirá el estado del sistema. El análisis cualitativo se realiza geométricamente, mediante representaciones gráficas, sin aplicar métodos formales de resolución de la ecuación, usando algunas veces métodos iterativos.  Palabras clave: caos, órbita, período, atractor, repulsor Abstract. The aim of this paper is to introduce the chaos concept. This is done through a qualitative analysis of discrete dynamical systems, whose deterministics models are given by non-linear differences equations. In each one, discrete time of same length is the independent variable and at the end of each interval the dependent variable will describe the state of the system. The qualitative analysis was conducted geometrically, using representations that let us to describe the change in state variables, whitout analitically resolution of equations but, using iterative methods if it is necessary. Key words: chaos, orbit, period, attracting, repelling !Introducción  Definiremos conceptos esenciales para el análisis cualitativo de una ecuación en diferencias ( )kk xfx =+1 , donde ( ) ( )fdomfimag ⊆ . Si ( )fdomx ∈0  la ecuación ( )kk xfx =+1  genera por iteración de f la sucesión ( ) ( );...;...;; 1010 kk xfxxfxx == + , que es la órbita de 0x , la cual se denotará { };...;...;;; 1100 +=Θ kkx xxxx . Diremos que x~  es punto fijo de ( )xf  cuando es abscisa de un punto de intersección de las gráficas de y=x con la de ( )xfy = , es decir, ( )xfx ~~ = . Un punto fijo x~  de ( )xf  es atractor, si existe algún entorno I de x~ , tal que para todo Ix ∈0 , su órbita converge a x~ , permaneciendo por tanto en dicho entorno. Un punto fijo x~  de ( )xf  es repulsor si existe algún entorno I de x~ , tal que para todo Ix ∈0  con xx ~0 ≠ , su órbita diverge de x~  saliendo por tanto de dicho entorno. El análisis cualitativo del sistema puede encararse: i) Graficando para algunos valores iniciales de 0x  los pares ( )kxk;  correspondientes a los primeros valores orbitales kx  de la órbita de 0x , con lo que tendremos una idea sobre el comportamiento de la sucesión correspondiente a ese 0x . ii) Mediante la representación para un determinado valor inicial 0x  de los pares ordenados ( )1; +kk xx  correspondientes a sucesivas iteraciones de f, obteniendo la línea de fases de 0x . La ANÁLISIS CUALITATIVO DE ECUACIONES EN DIFERENCIAS Y CAOS Juan Carlos Bressan y Ana Ferrazzi de Bressan Universidad de Buenos Aires  Universidad Argentina de la Empresa Argentina bressanjuancarlos22@gmail.com, anaferrazzibressan@gmail.com Acta Latinoamericana de Matemática Educativa 27 !Comité Latinoamericano de Matemática Educativa A. C. !648 reiteración del procedimiento para cada fDomx ∈0  genera la gráfica de f que conjuntamente con la bisectriz del 1ero y 3er cuadrantes permiten construir el diagrama de red. Para las funciones lineales ( ) baxxf +=  se sabe que un punto fijo x~  es atractor si 1<a , mientras que es repulsor si 1>a . Para cualquier función f continua y derivable en el intervalo de definición del sistema, el análisis del carácter de x~  puede encararse mediante la linealización de la función en x~ , válida en un intervalo infinitesimal alrededor de ese punto. Pero la pendiente de la función linealizada es ( )xf ~! ; luego aplicamos el siguiente criterio: Si ( ) 1~ <! xf , entonces x~  es atractor y si ( ) 1~ >! xf , entonces x~  es repulsor. Un modelo determinista caótico Consideremos la ecuación en diferencias kkk xxx 221 −=+ , asociada a la función ( ) ( )222 −=−= xxxxxf . Para su análisis cualitativo: (i) Obtenemos los puntos fijos, 0~1 =x  y 3~2 =x , resolviendo xxx 22 −= . (ii) Determinamos el comportamiento en cada x~ . Como ( ) 120 >=!f  y ( ) 143 >=!f , ambos puntos fijos son repulsores. (iii) Calculamos y representamos los primeros valores de la órbita para diversos 0x , por ejemplo si 1.00 =x  su órbita es { };...950559883.2;98760154.2;..;65906079.0;4161.0;19.0;1.01.00 −−=Θ =x . En la representación de los primeros valores de las órbitas de kkk xxx 221 −=+  se utilizan líneas de sostén para su mejor visualización. Esto puede verse en la Fig.1 donde representamos dichos valores orbitales para: 2,00 =x ; 3,00 =x ; 2,00 −=x ; 3,00 −=x ; 8,20 =x  y 7,20 =x  todos ellos pertenecientes al intervalo ( )3;1−=I . 1 2 3 4 5 6 7 8 9-1123kx( k)  Figura 1: Valores orbitales para distintas condiciones iniciales ( )3;10 −∈x   Capítulo 2. Propuestas para la enseñanza de las matemáticas !Comité Latinoamericano de Matemática Educativa A. C. !649 En la Fig. 1 podemos observar: 1º) Los puntos orbitales de cada órbita, están comprendidos en el intervalo ( )3;1−=I . 2º) Dividiendo, para cualquiera de las órbitas, el intervalo de valores orbitales en subintervalos de 5 décimas, aún con las pocas iteraciones efectuadas, son numerosos los subintervalos de esta órbita, “visitados” por otros valores orbitales de la misma órbita. Esta característica es la propiedad de mezcla. Para ilustrar otra característica de ciertas órbitas de este modelo, representamos los primeros puntos orbitales de 095.0−Θ  y de 0951.0−Θ  (ambos 0x  menores que 0~1 =x ). !1,5!1!0,500,511,522,533,51 4 7 10 13 16 19 22 25 28 31 34 37 40 43 46 49 52 55 58 61 64 67 70 73 76 79Series1Series2 Figura 2: Valores orbitales de órbitas de 095.0−Θ  y de 0951.0−Θ   La Fig. 2 muestra que la variación de un diezmilésimo en la condición inicial considerada, al avanzar el número de iteraciones, provoca variaciones inicialmente poco notorias en los valores orbitales, hasta producir valores que si bien fluctúan en el mismo intervalo, lo hacen con valores, que en ciertos intervalos, están totalmente diferenciados. Existe sensibilidad de las condiciones iniciales. Por ser 0~10 == xx  y 3~20 == xx  puntos fijos de la función, son de equilibrio del modelo, siendo sus órbitas periódicas { };...0..;;0;0;01 0~ =Θ =x ; { };...3..;;3;3;31 3~ =Θ =x . Además, 20 =x  (cero de f ); 10 =x  y 10 −=x  (abscisa y ordenada del vértice de f ), tienen por órbitas { }...;...;;;; 00022 =Θ ; { };...3;...;3;3;1;11 −=Θ  y { };...3;...;3;3;11 −=Θ− , que son eventualmente periódicas, pues al cabo de las primeras iteraciones se transforman en periódicas, en este caso de período 1. Nos resta ver el comportamiento de las órbitas cuando ( ) ( )+∞∪−∞−∈ ;31;0x . Puesto que ( ) ( )( ) ( )+∞=+∞∪−∞− ;3;31;f , para cualquier condición inicial 10 −<x  o bien 30 >x , será 3>kx  para 1≥k . Para el número creciente de iteraciones los valores orbitales también son crecientes, sin que existan intervalos de fluctuación para los puntos de sus trayectorias. Observamos en la Tabla 1 una divergencia constante de las órbitas a más infinito y que los valores Acta Latinoamericana de Matemática Educativa 27 !Comité Latinoamericano de Matemática Educativa A. C. !650 orbitales para 01,10 −=x  y 01,30 =x  son iguales; ello se debe a que ambos valores son simétricos con respecto a la abscisa del vértice de la parábola ( ) xxxf 22 −=  k 0 1 2 3 4 5 6 7 xk 3,01 3,0401 3,1620 3,6743 6,1518 25,5407 601,2457 360293.92 xk 3,009 3,0361 3,1456 3,6037 5,7793 21,8418 433,3796 186951,12 xk -1,01 3,0401 3,1620 3,6743 6,1518 25,5407 601,2457 360293,92 xk -1,02 3,0804 3,3281 4,4199 10,6956 93,0046 8463,85 71619827 Tabla 1: Valores orbitales de de kkk xxx 221 −=+  para 0093013 00 .;. == xx ; 021011 00 .;, −=−= xx  A continuación relacionaremos el comportamiento del modelo con la función compuesta ( ) ( ) ( )( )222))(( 22 −−−=== xxxxxffxfxg , cuyos puntos fijos son 0~1 =x , 3~2 =x , ( ) ( ) 52121~3 +=x  y ( ) ( ) 52121~4 −=x . Descartamos 1~x  y 2~x , por ser puntos fijos de f , quedando 3~x  y 4~x . Puesto que ( ) ( ) 102.4~~ 43 >≅"=" xgxg  resulta que 3~x  y 4~x  son puntos fijos repulsores de ( ) ( )xfxg 2= . Por otra parte, como ( ) 43 ~~ xxf =  y ( ) 34 ~~ xxf =  se obtienen los ciclos { }...;~;~...;;~;~ 434323~ xxxxx =Θ  y { }...;~;~...;;~;~ 343424~ xxxxx =Θ , que serán repulsores por ser 3~x  y 4~x  puntos fijos repulsores de ( ) ( )xfxg 2= . Calculamos, ( ) ( ) ( )[ ] ( ) ( )( )[ ]22222223 −−−−−−−= xxxxxxxxxf , cuyos puntos fijos son 87939.0;80194.0;3473.1;5321.2;2470.2;55496.0;3;0 −− , donde salvo 0 y 3 los restantes valores están calculados en forma aproximada. Esos valores dan lugar a { };...80194.0;55496.0;2470.2...;;80194.0;55496.0;2470.23 2470.2~ −−=Θ =x , { };...3473.1;5321.2;87939.0...;;3473.1;5321.2;87939.03 87939.0~ −−=Θ −=x , que son ciclos repulsores de orden 3. Al ver que existen puntos periódicos de período 3 en I, en virtud del teorema de Li y Yorke existen puntos periódicos de todos los períodos en I (ver pág 204 de Martín, Morán y Reyes). Esto asegura la densidad de los puntos periódicos en el intervalo I.   Capítulo 2. Propuestas para la enseñanza de las matemáticas !Comité Latinoamericano de Matemática Educativa A. C. !651 Diagramas de red Otra forma de reconocer el comportamiento de un modelo, es mediante la construcción del diagrama de red. Para ello, dibujamos la línea de fases ( ) xxxf 22 −=  de la ecuación recursiva kkk xxx 221 −=+  y la bisectriz del primero y tercer cuadrantes. Así, a partir de la condición inicial 0x  elegida, se dibuja una línea vertical hasta intersecar la línea de fases en el punto ( )( ) ( )1000 ;; xxxfx = . Desde esta intersección se traza una línea horizontal hasta intersecar la bisectriz en el punto ( ) ( )( ) ( )1100 xxxfxf ;; = , determinando así la nueva abscisa 1x  para calcular ( )1xf . A partir de este punto trazamos una vertical que intersecará la línea de fases en ( )( ) ( )2111 xxxfx ;; = . Este procedimiento puede reiterarse obteniendo el diagrama de red que permite conocer los valores aproximados de la órbita de 0x , leídos sobre cualquiera de los ejes, como muestran las Figs. 4 y 5  Figura 4: Diagrama de red para ( ) xxxf 22 −=  de la condición inicial 200 ,−=x .  Figura 5: Diagrama de red para ( ) xxxf 22 −=  de la condición inicial 200 ,=x .  Si bien las figuras anteriores no pueden asegurar que el modelo sea caótico en [ ]3;1− , nos permiten ver que ese carácter está directamente relacionado con la grafica de la función y sus puntos fijos. Por otra parte, sabemos que la función logística ( ) ( )xxxF −= 144  tiene un comportamiento caótico en [ ]1;0 , lo cual se prueba estableciendo una biyección con la gráfica de la función tienda de campaña dada por ( ) xxT 2=  si [ ]21;0∈x  y por ( ) ( )xxT −= 12  si ( ]1;21∈x , que tiene un comportamiento caótico en [ ]1;0 , como se demuestra en el libro de Martín, Morán y Reyes. De esta forma, si existe una transformación de semejanza de la gráfica de ( ) ( )2−= xxxf  en [ ]3;1−  a la de ( ) ( )xxxF −= 144  en [ ]1;0 , dejando invariante la bisectriz del Acta Latinoamericana de Matemática Educativa 27 !Comité Latinoamericano de Matemática Educativa A. C. !652 1ro. y 3er. cuadrantes, tendremos aseverada la caoticidad de ( )xf  en [ ]3;1− . Aseguramos tal semejanza puesto que por una simetría central en ( )0;0 , la grafica de ( )xf  en [ ]3;1− , se transforma en la de ( ) ( )xxxf +−= 21  en [ ]1;3− . Luego, mediante una traslación dada por ( )3;3V  obtenemos ( ) ( )xxxf −= 42  en [ ]4;0 . Finalmente por una homotecia de razón 41  resulta la gráfica de la función logística ( ) ( )xxxF −= 144  en [ ]1;0 . Los puntos fijos 0~1 =x  y 3~2 =x  de ( )xf  se transforman, respectivamente, en los puntos fijos 431 =x y 02 =x  de ( )xF4 . Además, cualquier [ ]3;10 −∈x  se transforma en ( )( ) [ ]1;0341 0 ∈− x . De esta forma ( ) ( )2−= xxxf  tendrá una órbita caótica para [ ]3;10 −∈x  si y sólo si ( ) ( )xxxF −= 144  tiene una órbita caótica para ( )( ) [ ]1;0341 0 ∈− x  Un modelo determinista no caótico Consideremos la ecuación en diferencias 221 ++−=+ kkk xxx , asociada a ( ) 22 ++−= xxxf . Sus puntos fijos 21 −=x~  y 22 =x~ , son repulsores pues ( ) 1122~1 >+=! xf  y ( ) 1122~2 >+−=" xf , generando los ciclos { }...;...;;;; 22221 2 −−−−=Θ−  y { }...;...;;;; 22221 2 =Θ  de período 1. Además, el punto notable 21+ , que es abscisa de una intersección de la gráfica de ( )xf  con la recta ( ) 2−=xr , genera el ciclo eventualmente periódico { }...;2...;;2;2;2121 −−−+=Θ + . -5 -4 -3 -2 -1 1 2 3 4-16-14-12-10-8-6-4-224x( k)x( k+1) Figura 6: 22 21 =−= xx~;~  puntos fijos y 21+  punto notable de ( ) 22 ++−= xxxf  -2.5 -2.0 -1.5 -1.0 -0.5 0.5 1.0 1.5 2.0 2.5 3.0 3.5-3-2-1123xx( k+1) Figura 7: 0~3 =x  y 2~4 =x  puntos fijos que consideraremos de ( ) ( )xfxg 2=  La función compuesta ( ) ( ) ( ) ( ) 222 2222 +++−+++−−== xxxxxfxg  tiene como puntos fijos 2~1 −=x , 22 =x~ , 0~3 =x  y 2~4 =x . Descartamos 1~x  y 2~x  por ser puntos fijos de ( )xf . Capítulo 2. Propuestas para la enseñanza de las matemáticas !Comité Latinoamericano de Matemática Educativa A. C. !653 Así, { };...2;0...;;2;0;2;020 =Θ  y { };...0;2...;;0;2;0;222 =Θ  son ciclos periódicos de orden 2 repulsores, pues ( ) ( ) 1320 >−="=" gg . Las figuras 6 y 7, exhiben los puntos fijos de ( )xf  y de ( )xf 2 , y además el punto notable de ( )xf . Determinamos ahora los puntos fijos de ( ) ( ) ( )( ) 237147144 23457823 +−−++−+−=== xxxxxxxxffxfxh : 802.0− ; 555.0 ; 247.2 ; 2− ; 2  donde los tres primeros, que forman un ciclo de período 3, están calculados con una aproximación de tres decimales y los dos restantes se descartan por ser puntos fijos de ( )xf . Como ( ) ;802.0247.2 −=f  ( ) 555.0802.0 =−f  y ( ) 247.2555.0 =f , es { };...555.0;802.0;247.2...;;555.0;802.0;247.23 247.2 −−=Θ . Observemos en la Fig. 8 que esos tres valores son raíces dobles de ( ) xxh = , de los que la bisectriz del 1ro. y 3er. cuadrantes es su recta tangente. Como la derivada ( )xh!  en cualquiera de esos puntos es 1, nos dificulta determinar el comportamiento del ciclo. Para ello, tomamos como condición inicial ( )21;20868.20 +−∈=x , valor aproximado de un punto fijo repulsor de ( )xf 4 ; de allí que la órbita comienza tomando valores aproximados del ciclo repulsor de orden 4 { };...90396.0;6602.1;26799.0;0868.24 0868.2~ −=Θ =x ; que, después de un período confuso de la iteración 16 a la 69, se estabiliza en el ciclo atractor de orden 3, como puede verse en la Fig. 9. El comportamiento atractor del ciclo de orden tres determina el carácter no caótico de la ecuación en diferencias 221 ++−=+ kkk xxx  en [ ]21;2 +−   !1!0,500,511,522,51 5 9 13 17 21 25 29 33 37 41 45 49 53 57 61 65 69 73 77 81 85 89 93 97 101105109x(k)k Figura 8: Gráfica de ( ) ( )xfxh 3= . Figura 9: Representación de la órbita de 0868.20 =x . Ahora, en la Tabla 2 veremos el comportamiento del modelo cuando [ ]21;20 +−∉x . k 0 1 2 3 4 5 xk 2,5 -1,75 -2,8125 -8,72265625 -82,8073883 -6937,87095 xk -1,7 -2,59 -7,2981 -58,5603636 -3485,87655 -12154819,2 !)!2! )!1! 1! 2! 3!)!2!2!x!(!!!k!)!x!(!!!k!+!1!)!Acta Latinoamericana de Matemática Educativa 27 !Comité Latinoamericano de Matemática Educativa A. C. !654 Tabla 2: Valores de la órbita de 221 ++−=+ kkk xxx  para 520 ,=x  y 710 ,−=x . Puesto que ( ) ( )( ) ( )−∞−=+∞+∪−∞− ;2;212;f , para cualquier condición inicial 20 −<x  o bien 210 +>x , será 2−<kx  para 1≥k . Para el número creciente de iteraciones los valores orbitales son monótonos decrecientes, resultando −∞=++∞→1lim kkx . Conclusiones En los dos ejemplos desarrollados los puntos fijos de la función ( )xf  son repulsores. En ( ) ( )2−= xxxf  existe un intervalo de caoticidad mientras que en ( ) 22 ++−= xxxf  no hay órbitas caóticas, pues existe un ciclo atractor de orden 3. La existencia de ciclos periódicos atractores implica que el modelo no presenta órbitas caóticas en el intervalo considerado. En algunos casos el carácter de los ciclos de orden 1 directamente permite afirmar la no caoticidad del modelo. Así, en el sistema modelado por ( ) kkk xxx 4331 +=+ , los puntos fijos: 211 −=x; 212 =x y 03 =x dan ciclos de orden 1, siendo repulsores los generados por 1x y 2x, mientras que es atractor el generado por 3x. Luego, todas las órbitas dadas por condiciones iniciales ( )21;210 −∈x  tienden a 03 =x; mientras que si ( )+∞∈ ;210x  su órbita tenderá a ∞+ , y si ( )21;0 −∞−∈x  tenderá a ∞− . Las Figs. 10 y 11 exhiben respectivamente los puntos fijos, así como las órbitas correspondientes a las diversas situaciones que muestran la no caoticidad del modelo. -1 1-11x( k)x( k+1) -4 -2 2 4 6 8 10 12 14 16 18 20-1.4-1.2-1.0-0.8-0.6-0.4-0.20.20.40.60.81.01.21.4kx( k) Figura 10. Puntos fijos de ( ) xxxf 43)(3 += : ( )211 −=x; ( )212 =x; 03 =x. Figura 11: Órbitas de la ecuación en diferencias ( ) kkk xxx 4331 +=+   Capítulo 2. Propuestas para la enseñanza de las matemáticas !Comité Latinoamericano de Matemática Educativa A. C. !655 Probar que un modelo es caótico en un intervalo I  requiere demostrar que: (i) Los puntos periódicos son densos en I ; (ii) Las órbitas son sensibles a las condiciones iniciales y (iii) Tiene la propiedad de mezcla, es decir, cualesquiera sea la división en subintervalos  nIII ,...,, 21  de I , para cualquier jI , existe algún jIx ∈0  tal que su órbita 0xΘ  tiene algún punto en cada uno de los subintervalos nIII ,...,, 21  de I . Esta propiedad por un resultado de M. Vellekoop y R. Berglund, cumple: Si I  es un intervalo finito o infinito de los números reales, la ecuación en diferencias ( )kk xfx =+1  con f  continua, es caótica en el intervalo I  si y sólo si tiene la propiedad de mezcla. Esta propiedad no es fácil de probar, de allí que buscamos el ejemplo del parágrafo 2, que permitió inferir las características de un modelo que tiene un intervalo de caoticidad, del cual por semejanza con la ecuación logística, queda indirectamente demostrado el comportamiento caótico del modelo en [ ]3;1−  Referencias bibliográficas Blanchard, P., Devaney, R. L. y Hall, G. R. (1999). Ecuaciones diferenciales. México: Ed. Thomson, Fernández Pérez, C., Vázquez Hernández, F. J. y Vegas Montaner, J. M. (2003). Ecuaciones diferenciales y en diferencias. Sistemas dinámicos. España: Ed. Thomson. Martín, M. A., Morán, M. y Reyes, M. (1998); Iniciación al caos. Sistemas dinámicos. España: Ed. Síntesis. 

Análisis cualitativo de ecuaciones en diferencias y caos

http://funes.uniandes.edu.co/12419/1/Bressan2014Analisis.pdf

Análisis cualitativo de ecuaciones en diferencias y caos

Abstract

Similar works

Full text

Available Versions

Funes