Articulo de referencia

Función de Feigenbaum

En el estudio de sistemas dinámicos, el término función de Feigenbaum se ha utilizado para describir dos funciones diferentes introducidas por el físico Mitchell Feigenbaum : [ ...

En el estudio de sistemas dinámicos, el término función de Feigenbaum se ha utilizado para describir dos funciones diferentes introducidas por el físico Mitchell Feigenbaum : [ 1 ]

Idea

Ruta de duplicación de período hacia el caos

En el mapa logístico,

tenemos una funciónFr(incógnita)=rincógnita(1incógnita){\displaystyle f_{r}(x)=rx(1-x)}y queremos estudiar qué sucede cuando iteramos el mapa muchas veces. El mapa podría caer en un punto fijo , un ciclo fijo o el caos. Cuando el mapa cae en un ciclo fijo estable de longitudnorte{\displaystyle n}, encontraríamos que la gráfica deFrnorte{\displaystyle f_{r}^{n}}y el gráfico deincógnitaincógnita{\displaystyle x\mapsto x}se cruza ennorte{\displaystyle n}puntos y la pendiente de la gráfica deFrnorte{\displaystyle f_{r}^{n}}está delimitado en(1,+1){\displaystyle (-1,+1)}en esas intersecciones.

Por ejemplo, cuandor=3.0{\displaystyle r=3.0}, tenemos una única intersección, con pendiente limitada en(1,+1){\displaystyle (-1,+1)}, lo que indica que se trata de un único punto fijo estable.

Comor{\displaystyle r}aumenta a más allár=3.0{\displaystyle r=3.0}, el punto de intersección se divide en dos, lo que es una duplicación de período. Por ejemplo, cuandor=3.4{\displaystyle r=3.4}Hay tres puntos de intersección, siendo el del medio inestable y los otros dos estables.

Comor{\displaystyle r}aprochesr=3.45{\displaystyle r=3.45}, se produce otra duplicación de período de la misma manera. Las duplicaciones de período ocurren cada vez con mayor frecuencia, hasta que en cierto puntor3.56994567{\displaystyle r\approx 3.56994567}Las duplicaciones de período se vuelven infinitas y el mapa se vuelve caótico. Esta es la ruta de duplicación de período hacia el caos .

Relación entreincógnitanorte+2{\displaystyle x_{n+2}}yincógnitanorte{\displaystyle x_{n}}cuandoa=2.7{\displaystyle a=2.7}Antes de que ocurra la bifurcación de duplicación de período, la órbita converge al punto fijo.incógnitaF2{\displaystyle x_{f2}}.
Relación entreincógnitanorte+2{\displaystyle x_{n+2}}yincógnitanorte{\displaystyle x_{n}}cuandoa=3{\displaystyle a=3}La pendiente tangente en el punto fijoincógnitaF2{\displaystyle x_{f2}}. es exactamente 1, y se produce una bifurcación de duplicación de período.
Relación entreincógnitanorte+2{\displaystyle x_{n+2}}yincógnitanorte{\displaystyle x_{n}}cuandoa=3.3{\displaystyle a=3.3}El punto fijoincógnitaF2{\displaystyle x_{f2}}se vuelve inestable, dividiéndose en un ciclo estable periódico-2.
Cuandor=3.0{\displaystyle r=3.0}, tenemos una única intersección, con pendiente exactamente+1{\displaystyle +1}, lo que indica que está a punto de experimentar una duplicación de período.
Cuandor=3.4{\displaystyle r=3.4}Hay tres puntos de intersección, siendo el del medio inestable y los otros dos estables.
Cuandor=3.45{\displaystyle r=3.45}, hay tres puntos de intersección, siendo el del medio inestable y los otros dos con pendiente exactamente+1{\displaystyle +1}, lo que indica que está a punto de experimentar otra duplicación de período.
Cuandor3.56994567{\displaystyle r\approx 3.56994567}, hay infinitas intersecciones, y hemos llegado al caos a través de la ruta de duplicación de período .

Límite de escala

Aproximación al límite de escala comor{\displaystyle r}aprochesr=3.5699{\displaystyle r^{*}=3.5699\cdots }desde abajo.
En el punto del caosr=3.5699{\displaystyle r^{*}=3.5699\cdots }, mientras repetimos las duplicaciones de períodoFr1,Fr2,Fr4,Fr8,Fr16,{\displaystyle f_{r^{*}}^{1},f_{r^{*}}^{2},f_{r^{*}}^{4},f_{r^{*}}^{8},f_{r^{*}}^{16},\dots }Los gráficos parecen asemejarse entre sí, excepto que están encogidos hacia el centro y rotados 180 grados, convergiendo en un fractal.

Al observar las imágenes, se puede notar que en el punto del caosr=3.5699{\displaystyle r^{*}=3.5699\cdots }, la curva deFr{\displaystyle f_{r^{*}}^{\infty }}parece un fractal. Además, a medida que repetimos las duplicaciones de períodoFr1,Fr2,Fr4,Fr8,Fr16,{\displaystyle f_{r^{*}}^{1},f_{r^{*}}^{2},f_{r^{*}}^{4},f_{r^{*}}^{8},f_{r^{*}}^{16},\dots }Los gráficos parecen similares entre sí, excepto que están reducidos hacia el centro y rotados 180 grados.

Esto nos sugiere un límite de escala: si duplicamos repetidamente la función, entonces debemos escalarla porα{\displaystyle \alpha }para una cierta constanteα{\displaystyle \alpha }:F(incógnita)αF(F(incógnita/α)){\displaystyle f(x)\mapsto -\alpha f(f(-x/\alpha ))}entonces, en el límite, terminaríamos con una funcióngramo{\displaystyle g}que satisfacegramo(incógnita)=αgramo(gramo(incógnita/α)){\displaystyle g(x)=-\alpha g(g(-x/\alpha ))}Además, a medida que los intervalos de duplicación de período se vuelven cada vez más cortos, la relación entre dos intervalos de duplicación de período converge a un límite, la primera constante de Feigenbaum.δ=4.6692016{\displaystyle \delta =4.6692016\cdots }.

Para valores incorrectos del factor de escalaα{\displaystyle \alpha }, el mapa no converge a un límite, pero cuandoα=2.5029{\displaystyle \alpha =2.5029\dots }, converge.
En el punto del caosr=3.5699{\displaystyle r^{*}=3.5699\cdots }, mientras repetimos la iteración de la ecuación funcionalF(incógnita)αF(F(incógnita/α)){\displaystyle f(x)\mapsto -\alpha f(f(-x/\alpha ))}conα=2.5029{\displaystyle \alpha =2.5029\dots }, encontramos que el mapa sí converge a un límite.

La constanteα{\displaystyle \alpha }se puede encontrar numéricamente probando muchos valores posibles. Para los valores incorrectos, el mapa no converge a un límite, pero cuando esα=2.5029{\displaystyle \alpha =2.5029\dots }, converge. Esta es la segunda constante de Feigenbaum.

Régimen caótico

En el régimen caótico,Fr{\displaystyle f_{r}^{\infty }}, el límite de las iteraciones del mapa, se convierte en bandas oscuras caóticas intercaladas con bandas brillantes no caóticas.

En el régimen caótico,Fr{\displaystyle f_{r}^{\infty }}, el límite de las iteraciones del mapa, se convierte en bandas oscuras caóticas intercaladas con bandas brillantes no caóticas.

Otros límites de escala

Cuandor{\displaystyle r}aprochesr3.8494344{\displaystyle r\approx 3.8494344}, tenemos otro enfoque de duplicación de período para el caos, pero esta vez con períodos 3, 6, 12, ... Esto nuevamente tiene las mismas constantes de Feigenbaumδ,α{\displaystyle \delta ,\alpha }. El límite deF(incógnita)αF(F(incógnita/α)){\textstyle f(x)\mapsto -\alpha f(f(-x/\alpha ))}También es la misma función. Este es un ejemplo de universalidad .

Mapa logístico que se aproxima al límite de escalamiento del caos de duplicación de período.r=3.84943{\displaystyle r^{*}=3.84943\dots }desde abajo. En el límite, esto tiene la misma forma que la der=3.5699{\displaystyle r^{*}=3.5699\cdots }, puesto que todas las rutas de duplicación de período hacia el caos son las mismas (universalidad).

También podemos considerar la ruta de triplicación de períodos hacia el caos eligiendo una secuencia der1,r2,{\displaystyle r_{1},r_{2},\dots }de tal manera quernorte{\displaystyle r_{n}}es el valor más bajo del período-3norte{\displaystyle 3^{n}}ventana del diagrama de bifurcación . Por ejemplo, tenemosr1=3.8284,r2=3.85361,{\displaystyle r_{1}=3.8284,r_{2}=3.85361,\dots }, con el límiter=3.854077963{\displaystyle r_{\infty }=3.854077963\dots }. Este tiene un par diferente de constantes de Feigenbaum.δ=55.26,α=9.277{\displaystyle \delta =55.26\dots ,\alpha =9.277\dots }. [ 2 ] YFr{\displaystyle f_{r}^{\infty }}converge al punto fijo aF(incógnita)αF(F(F(incógnita/α))){\displaystyle f(x)\mapsto -\alpha f(f(f(-x/\alpha )))}Como otro ejemplo, period-4-pling tiene un par de constantes de Feigenbaum distintas de las de period-doubling, aunque period-4-pling se alcanza mediante dos duplicaciones de periodo. En detalle, definar1,r2,{\displaystyle r_{1},r_{2},\dots }de tal manera quernorte{\displaystyle r_{n}}es el valor más bajo del período-4norte{\displaystyle 4^{n}}ventana del diagrama de bifurcación. Entonces tenemosr1=3.960102,r2=3.9615554,{\displaystyle r_{1}=3.960102,r_{2}=3.9615554,\dots }, con el límiter=3.96155658717{\displaystyle r_{\infty }=3.96155658717\dots }. Este tiene un par diferente de constantes de Feigenbaum.δ=981.6,α=38,82{\displaystyle \delta =981.6\dots ,\alpha =38.82\dots }.

En general, cada ruta de multiplicación de periodo hacia el caos tiene su propio par de constantes de Feigenbaum. De hecho, normalmente hay más de una. Por ejemplo, para el proceso de multiplicación de periodo 7, existen al menos 9 pares diferentes de constantes de Feigenbaum. [ 2 ]

Generalmente,3δ2α2{\textstyle 3\delta \approx 2\alpha ^{2}}y la relación se vuelve exacta a medida que ambos números aumentan hasta el infinito:límiteδ/α2=2/3{\displaystyle \lim \delta /\alpha ^{2}=2/3}.

ecuación funcional de Feigenbaum-Cvitanović

Esta ecuación funcional surge en el estudio de mapas unidimensionales que, en función de un parámetro, experimentan una cascada de duplicación de período. Descubierta por Mitchell Feigenbaum y Predrag Cvitanović , [ 3 ] la ecuación es la expresión matemática de la universalidad de la duplicación de período. Especifica una función g y un parámetro α mediante la relación

gramo(incógnita)=αgramo(gramo(incógnita/α)){\displaystyle g(x)=-\alpha g(g(-x/\alpha ))}

con las condiciones iniciales{gramo(0)=1,gramo(0)=0,gramo(0)<0.{\displaystyle {\begin{cases}g(0)=1,\\g'(0)=0,\\g''(0)<0.\end{cases}}}Para una forma particular de solución con una dependencia cuadrática de la solución cerca de x = 0, α = 2,5029... es una de las constantes de Feigenbaum .

La serie de poder degramo{\displaystyle g}es aproximadamente [ 4 ]gramo(incógnita)=11.52763incógnita2+0,104815incógnita4+0,026705incógnita6+O(incógnita8){\displaystyle g(x)=1-1.52763x^{2}+0.104815x^{4}+0.026705x^{6}+O(x^{8})}

Renormalización

La función de Feigenbaum se puede derivar mediante un argumento de renormalización . [ 5 ]

La función de Feigenbaum satisface [ 6 ]gramo(incógnita)=límitenorte1F(2norte)(0)F(2norte)(incógnitaF(2norte)(0)){\displaystyle g(x)=\lim _{n\to \infty }{\frac {1}{F^{\left(2^{n}\right)}(0)}}F^{\left(2^{n}\right)}\left(xF^{\left(2^{n}\right)}(0)\right)}para cualquier mapa en la línea realF{\displaystyle F}al comienzo del caos.

Función de escala

La función de escala de Feigenbaum proporciona una descripción completa del atractor del mapa logístico al final de la cascada de duplicación de periodo. El atractor es un conjunto de Cantor y, al igual que el conjunto de Cantor del tercio medio, puede cubrirse con un conjunto finito de segmentos, todos mayores que un tamaño mínimo d n . Para un d n fijo, el conjunto de segmentos forma una cobertura Δ n del atractor. La razón de los segmentos de dos coberturas consecutivas, Δ n y Δ n +1, puede aproximarse a una función σ , la función de escala de Feigenbaum.

Véase también

Notas

  1. Feigenbaum, MJ (1976) "Universalidad en la dinámica discreta compleja", Informe anual de la División Teórica de Los Alamos 1975-1976
  2. 1 2 Delbourgo, R.; Hart, W.; Kenny, BG (1985-01-01). "Dependencia de las constantes universales del período de multiplicación en mapas no lineales" . Physical Review A. 31 ( 1): 514– 516. Bibcode : 1985PhRvA..31..514D . doi : 10.1103/PhysRevA.31.514 . ISSN 0556-2791 . PMID 9895509 .  
  3. La nota a pie de página en la página 46 de Feigenbaum (1978) afirma: "Esta ecuación exacta fue descubierta por P. Cvitanović durante una discusión y en colaboración con el autor".
  4. Iii, Oscar E. Lanford (mayo de 1982). "Una demostración asistida por computadora de las conjeturas de Feigenbaum" . Boletín (Nueva Serie) de la Sociedad Matemática Americana . 6 (3): 427– 434. doi : 10.1090/S0273-0979-1982-15008-X . ISSN 0273-0979 . 
  5. Feldman, David P. (2019). Caos y sistemas dinámicos . Princeton. ISBN 978-0-691-18939-0OCLC 1103440222 .​ {{cite book}}: CS1 mantenimiento: falta el editor de ubicación ( enlace )
  6. Weisstein, Eric W. "Función de Feigenbaum" . mathworld.wolfram.com . Consultado el 7 de mayo de 2023 .

Bibliografía

  • Feigenbaum, M. (1978). "Universalidad cuantitativa para una clase de transformaciones no lineales". Journal of Statistical Physics . 19 (1): 25– 52. Bibcode : 1978JSP....19...25F . CiteSeerX 10.1.1.418.9339 . doi : 10.1007 /BF01020332 . MR 0501179. S2CID 124498882 .   
  • Feigenbaum, M. (1979). "Las propiedades métricas universales de las transformaciones no lineales". Journal of Statistical Physics . 21 ( 6): 669– 706. Bibcode : 1979JSP....21..669F . CiteSeerX 10.1.1.418.7733 . doi : 10.1007/BF01107909 . MR 0555919. S2CID 17956295 .   
  • Feigenbaum, Mitchell J. (1980). "La transición al comportamiento aperiódico en sistemas turbulentos" . Communications in Mathematical Physics . 77 (1): 65– 86. Bibcode : 1980CMaPh..77...65F . doi : 10.1007/BF01205039 . S2CID 18314876 . 
  • Epstein, H.; Lascoux, J. (1981). "Propiedades de analiticidad de la función de Feigenbaum" . Commun. Math. Phys . 81 (3): 437– 453. Bibcode : 1981CMaPh..81..437E . doi : 10.1007/BF01209078 . S2CID 119924349 . 
  • Feigenbaum, Mitchell J. (1983). "Comportamiento universal en sistemas no lineales". Physica . 7D ( 1– 3): 16– 39. Bibcode : 1983PhyD....7...16F . doi : 10.1016/0167-2789(83)90112-4 .Bound as Order in Chaos, Actas de la Conferencia Internacional sobre Orden y Caos celebrada en el Centro de Estudios No Lineales, Los Alamos, Nuevo México 87545, EE. UU., del 24 al 28 de mayo de 1982 , Eds. David Campbell, Harvey Rose; North-Holland Amsterdam ISBN 0-444-86727-9.
  • Lanford III, Oscar E. (1982). "Una demostración asistida por computadora de las conjeturas de Feigenbaum" . Bull. Am. Math. Soc . 6 (3): 427– 434. doi : 10.1090/S0273-0979-1982-15008-X . MR 0648529 . 
  • Campanino, M.; Epstein, H.; Ruelle, D. (1982). "Sobre la ecuación funcional de Feigenbaum"gramogramo(λincógnita)+λgramo(incógnita)=0{\displaystyle g\circ g(\lambda x)+\lambda g(x)=0}". Topología . 21 (2): 125– 129. doi : 10.1016/0040-9383(82)90001-5 . MR 0641996 . 
  • Lanford III, Oscar E. (1984). "Una prueba más breve de la existencia del punto fijo de Feigenbaum". Commun. Math. Phys . 96 (4): 521– 538. Bibcode : 1984CMaPh..96..521L . CiteSeerX 10.1.1.434.1465 . doi : 10.1007/BF01212533 . S2CID 121613330 .  
  • Epstein, H. (1986). "Nuevas pruebas de la existencia de las funciones de Feigenbaum" (PDF) . Commun. Math. Phys . 106 (3): 395– 426. Bibcode : 1986CMaPh.106..395E . doi : 10.1007/BF01207254 . S2CID 119901937 . 
  • Eckmann, Jean-Pierre ; Wittwer, Peter (1987). "Una demostración completa de las conjeturas de Feigenbaum". J. Stat. Phys . 46 (3/4): 455. Bibcode : 1987JSP ....46..455E . doi : 10.1007/BF01013368 . MR 0883539. S2CID 121353606 .  
  • Stephenson, John; Wang, Yong (1991). "Relaciones entre las soluciones de la ecuación de Feigenbaum" . Appl. Math. Lett . 4 (3): 37– 39. doi : 10.1016/0893-9659(91)90031-P . MR 1101871 . 
  • Stephenson, John; Wang, Yong (1991). "Relaciones entre autofunciones asociadas con soluciones de la ecuación de Feigenbaum" . Appl. Math. Lett . 4 (3): 53– 56. doi : 10.1016/0893-9659(91)90035-T . MR 1101875 . 
  • Briggs, Keith (1991). "Un cálculo preciso de las constantes de Feigenbaum" . Math. Comp . 57 (195): 435– 439. Bibcode : 1991MaCom..57..435B . doi : 10.1090/S0025-5718-1991-1079009-6 . MR 1079009 . 
  • Tsygvintsev, Alexei V.; Mestel, Ben D.; Obaldestin, Andrew H. (2002). «Fracciones continuas y soluciones de la ecuación de Feigenbaum-Cvitanović» . Cuentas Rendus de la Academia de Ciencias, Serie I. 334 (8): 683– 688. doi : 10.1016/S1631-073X(02)02330-0 .
  • Mathar, Richard J. (2010). "Representación en serie de Chebyshev de la función de duplicación de período de Feigenbaum". arXiv : 1008.4608 [ math.DS ].
  • Varin, VP (2011). "Propiedades espectrales del operador de duplicación de período" . KIAM Preprint . 9 . arXiv : 1202.4672 .
  • Weisstein, Eric W. "Función Feigenbaum" . MundoMatemático .