SciELO - Scientific Electronic Library Online

 
vol.41 issue2A first look at knot theory and Khovanov homology author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

  • On index processCited by Google
  • Have no similar articlesSimilars in SciELO
  • On index processSimilars in Google

Share


Revista Integración

Print version ISSN 0120-419XOn-line version ISSN 2145-8472

Integración - UIS vol.41 no.2 Bucaramanga July/Dec. 2023  Epub Aug 09, 2023

https://doi.org/10.18273/revint.v41n2-2023004 

Artículos originales

Aspectos dinámicos del acoplamiento skew de la familia logística

Dynamical aspects of skew coupling of the logistic family

Osvaldo Osuna a  

Cristian Jesús Rojas-Milla b  

a Universidad Michoacana, Instituto de Física y Matemáticas, Morelia, México. osvaldo.osuna@umich.mx

b Universidad del atlántico, Departamento de matemática, Barranquilla, Colombia. E-mail: cristianrojas@mail.uniatlantico.edu.co.


Resumen.

En este trabajo el objetivo fundamental es estudiar algunos aspectos de la evolución asintótica de las órbitas obtenidas por iterar el endomorfismo F µ,c (x, y)= (f µ (x), f µ (y) + c(xy)), donde μ > 1 es el parámetro de la familia logística: f µ (x) = µx(1 x) y (0, 1), y ε ϵ (0;1), es el parámetro de acoplamiento. Este mapa biparamétrico es un híbrido entre dos clásicos en la teoría de los sistemas dinámicos, por un lado el paradigmático mapa cuadrático y por el otro el acoplamiento skew. El resultado principal será construir de manera detallada para determinados parámetros (μ, ϵ ) un conjunto compacto invariante. Así como estudiar el comportamiento asintótico dentro de dicho compacto y su complemento, para ello de obtiene una descripción del comportamiento de las preimágenes de zonas en ℝ2 que juegan un rol importante en la comprensión de la dinámica del acoplamiento.

MSC2010:

37C05, 37N35, 37C70.

Palabras clave: Mapa logístico; acoplamiento skew; conjunto compacto invariante

Abstract.

Our main aim is to study some aspects of the asymptotic evolution of the orbits obtained by iterating the endomorphism F µ,c (x, y) = (f µ (x), f µ (y) + c(xy)), where > 1 is the parameter of the logistic family: f (x) = x(1 x) and ε ϵ (0;1); is the coupling parameter. This biparametric map is hybrid between two classics in the theory of dynamical systems, the paradigmatic quadratic map and the skew coupling. The main result will be to construct in detail for certain parameters (μ, ϵ ) an invariant compact set. As well as a study of the asymptotic behavior within this compact and its complement, for this we obtain a description of the behavior of the preimages of zones in ℝ2 that play an important role in understanding the dynamics of the coupling.

Keywords: Logistic map; skew coupling; invariant compact set

1. Introducción y antecedentes

Un número considerable de modelos matemáticos en dinámica discreta, incluso provenientes de diferentes procesos de discretización de ecuaciones diferenciales, conducen a una clase especial de sistemas dinámicos, los llamados skew products (por su terminología en inglés). Estos fueron definidos por Anzai [4] en 1951 para producir cierta clase de contraejemplos en teoría ergódica; tales sistemas aparecen de manera natural en el estudio de automorfismos de nilvariedades, ver [10]. Esta clase de sistemas dinámicos en forma general se definen de la siguiente manera: dados dos espacios topológicos X y Y , un endomorfismo F: X × YX × Y es llamado skew product si existen aplicaciones f: XX y g: X × YY tales que

El estudio de las transformaciones skew product tiene una amplia literatura matemática, por ejemplo: En [1], [2] se estudian propiedades ergódicas los skew productos, en particular, se establecen la fórmula de Abramov- Roklhin y así como cálculos de la entropía para skew productos. En [3] se analizan atractores extraños de tales mapas. Puede encontrarse diversas propiedades de bifurcación de mapas skew en [5]. En [6] se presentan conexiones de productos skew con la dinámica de mapas en el intervalo. Se presentan propiedades dinámicas a través de derivada Schwarziana en [7]. En [11] y [14], se analizan conjuntos invariantes relevantes de dichos mapas.

Dada una función f: ℝ → ℝ, digamos continua, denominamos acoplamiento skew de f al mapa F: ℝ2 → ℝ2 definido, para cada x; y 2 ℝ, por

Destacamos que esta clase particular de skew product ha sido considerado como modelo relevante para la transmisión de información a través de señales caóticas, la cual va acompañada frecuentemente de la noción de sincronización de sistemas; ver por ejemplo [8] y [9]. En este trabajo estudiaremos aspectos dinámicos del acoplamiento skew para la familia logística F μ, ϵ: ℝ2 → ℝ2 dado por

Donde μ > 1 y ε ϵ (0;1):

2. Puntos críticos e imagen crítica

Como el mapa anterior F μ, ϵ es cuadrático, entonces es no invertible y en consecuencia el conjunto J 0 de sus puntos críticos, el conjunto J 1 = F μ, ϵ (J 0) de sus valores críticos (imagen crítica), así como su preimagen F μ, ϵ1(J 1); tienen un papel importante en la descripción de la dinámica de F μ, ϵ. De hecho, J 1 divide el plano en regiones abiertas cada una con un número constante de preimágenes, mientras que F μ, ϵ1(J 1) divide ℝ2 en regiones abiertas que se mapean difeomórficamente sobre las regiones determinadas por J 1; más adelante precisaremos estos divisiones, antes recordamos que J 0 = f(x; y) ϵ2: det D (x;y) F μ, ϵ = 0g, donde D (x;y) F μ, ϵ denota la diferencial (o matriz jacobiana) de F F μ, ϵ en el punto (x; y). Dado que

J 0 es la unión de las rectas , ambas se

intersectan en el punto crítico siendo que . Sigue por tanto que ; de hecho

donde , para cualquier t ϵ ℝ. En consecuencia,

donde esto es, es un rayo vertical con punto inicial en es una parábola horizontal hacia la izquierda con vértice en . Las curvas que definen a J 1 dividen el plano en tres regiones abiertas: Z 0 ;Z 2 y Z 4. Estas regiones abiertas son definidas por:

cada una de estas están caracterizada por el constante número de preimágenes que tienen sus puntos; por ejemplo: en Z 0 no hay preimagen, en efecto, por definición (x 0 ; y 0) Є Z 0 siempre que x 0 >μ4. Así, si la preimagen de Z 0 fuera no vacía existiría (x; y) Є R 2 tal que f μ (x; y) = (x 0 ; y 0). Luego la ecuación cuadrática μx(1 x) = x 0 tendría solución real, por lo que su discriminante debe cumplir μ 2 - 4μx0 ≥ 0; así x 0 tendría que ser menor que μ4 lo cual contradice que esta en Z 0. De manera análoga se tiene que los elementos en Z 2 tienen dos preimagenes y los de Z 4 tienen cuatro, ver Figura 1.

Figura 1. Imagen crítica

Observe que si y sólo si z = c μ ; además, cada punto en F μ,Є (ι1) tiene exactamente dos preimágenes, estas están dispuestas en 1 simétricamente respecto de c μ ; así que la restricción de F μ,Є a ι1 es 2 a 1. Con respecto a ι2 la situación es diferente; en primer lugar, es simple ver que es inyectiva, sin embargo, las preimágenes de puntos en F μ,Є (ι2 ) requieren algo de atención. La rama superior de la parábol con tiene como preimagen 1 a 1 al rayo no obstante, cada punto en la rama inferior tiene tres preimágenes, con una de las cuales se recorre el rayo , las otras dos describen un par de curvas (del tipo ) simétricas respecto al rayo El hecho de tener tres zonas con número diferentes de preimágenes implica, entre otras cosas, que las dinámicas del acoplamiento skew de F μ,Є y el producto (Є = 0) no son topológicamente conjugadas. Por otra parte, las curvas que definen dividen el plano en seis zonas: Z i 2 (i = 1; 2) y Z j 4 (j = 1;2;3; 4); cada zona Z i 2 (resp. Z j 4 ) es mapeada de forma difeomórfica sobre Z 2 (resp. Z 4); ver Figura 2.

Figura 2. Ilustración de F μ,Є -1(J 1)

3. Puntos fijos y diagrama de bifurcación

Otros elementos importantes de considerar en la descripción de la dinámica de F ; son sus puntos fijos y la variación de su naturaleza hiperbólica. Es fácil ver que los puntos fijos de F μ,Є son: O denotará al origen del plano

Las respectivas matrices jacobianas en esos puntos son:

El siguiente cuadro refleja los autovalores y autoespacios relativos a cada punto fijo.

Punto fijo Autovalores Autoespacios asociados
O λ 1 = µ λ 2 = µc E 1 = {(x, y) : x = y} E 2 = {(x, y) : x = 0}
P µ,c λ 1 = µ λ 2 = −µ + c + 2
2−−2
Q µ,c λ 1 = 2 − µ λ 2 = 2 − µc E 1 = {(x, y) : x = y} E 2 = {(x, y) : x = 0}
R µ,c λ 1 = 2 − µ λ 2 = µ + c

En base a esta información se construye, en términos de los parámetros μ y Є, un diagrama de bifurcación para cada punto fijo. Por ejemplo, para el origen y P μ,Є tenemos en primer lugar, que estos puntos coinciden cuando μ = 1 + Є; de hecho, en la región 1 < μ < 1 + Є el origen O es silla hiperbólica y P μ,Є repulsor; cuando 1 + Є < μ < 3 + Є, O es repulsor y P μ,Є silla hiperbólica; mientras que si μ > 3 + Є ambos puntos son repulsores.

Dado que la fibra 0 x ℝ, es decir el eje y es invariante, la dinámica en ella es dada por la función h μ,Є (t) = f μ(t) Єt, la cual mediante el cambio de variables x = μ μ-ϵty de parámetro ρ= μ - Є es conjugada a . Similares análisis y comentarios se hacen para los restantes puntos fijos.

  • Los puntos fijos O y P μ,Є se asocian para producir una bifurcación transcrítica cuando μ= 1 + Є; como se establece en [13] [Capítulo 3], ver Figura 3 donde aparecen las líneas de bifurcación.

Figura 3 Diagrama de bifurcación para O y P μ,Є  

Para O siempre vale λ 1 > λ 2. Este punto fijo es silla hiperbólica cuando µ < 1 + Є y repulsor cuando µ > 1 + Є

● Cuando µ < 1 + Є el punto P µ,c es repulsor, en el intervalo 1 + Є < µ < 3 + Є es silla hiperbólica y vuelve a ser repulsor cuando µ > 1 + Є. Cuando µ = 3 + Є, λ 2 = -1 y se inicia una cascada de duplicación de período (ver [13] [Capítulo 3 ]), esta bifurcación se despliega sobre el conjunto {(x, y) : x = 0}, donde la dinámica es gobernada por h 0(t) = (µЄ)tµt 2. Sobre el autoespacio E 2 opera la menor expansión excepto en el caso µ < 1 +2.

●El punto fijo Q µ,Є es atractor si µ < 3 − Є, es silla hiperbólica cuando 3 − Є < µ < 3 y es repulsor si µ > 3. Cuando µ = 3− Є se tiene λ 2 = −1 y comienza un cascada de duplicación de período la cual se manifiesta en la fibra invariante y es gobernada por la dinámica de . Mientras que si µ = 3, entonces λ 1 = -1 y se da inicio a otra cascada de duplicación de período pero esta vez sobre la diagonal, cuya dinámica es regida por f µ .

● Para R µ, Є siempre se tiene λ 2 > |λ 1| y λ2 > 1. Si µ < 3, este punto fijo es una silla hiperbólica; y es repulsor cuando µ > 3. El autovalor λ1 alcanza el valor de -1 cuando µ = 3, ello marca el inicio de una cascada de duplicación de período esta vez sobre una curva invariante cuya dinámica es regida por f µ ; esto es mostrado en [12]. Ver Figura 4.

Figura 4 Diagrama de duplicación de periodo sobre una curva 

4. Infinito como atractor

El cuarto ingrediente en nuestro trabajo es la propiedad de atracción en el infinito. En general, un endomorfismo F: ℝ n → ℝ n admite a 1 como atractor si existe una vecindad w de 1 (complemento de un compacto) tal que F(W ) W y para cada z Є W su órbita positiva es no acotada: ||F n (z)k|| →+∞ cuando n →+∞. La cuenca de atracción de 1 es el conjunto B de todos los puntos cuyas órbita es no acotada; es muy simple mostrar que B es un conjunto abierto, de hecho . Se llama cuenca inmediata de 1 a la mayor componente conexa en B , se denota por B 0 .

Note que si W es como arriba, entonces un punto z tiene órbita acotada si z Є ℝ n / W . Destacamos que el compacto ℝ n / W no tiene porque ser F-invariante, de hecho puede contener componentes de la cuenca de atracción de 1. Tal es el caso de la familia logística f μ; es bien conocido que W = (-∞,0) (1,+ ∞) es f μ -invariante y todo punto en W tiene órbita no acotada para cualquier > 1. Sin embargo, para todo μ > 4, el compacto [0;1] tiene infinitas componentes de la cuenca de atracción de + ∞; es más, los puntos en [0;1] con órbita acotada definen un conjunto de Cantor con medida de Lebesgue nula. Note que (-∞,0) [ (1,+ ∞) es la cuenca inmediata de 1 para f μ.

Suponga que F: ℝ n → ℝ n tiene a ∞ como atractor y que W es como arriba, entonces Λ= ℝ n B es el conjunto de puntos cuya órbita (positiva) es acotada; se adopta la nomenclatura de la dinámica compleja y llamamos a conjunto de Julia de F. Si para cada subconjunto A de ℝ n hacemos A c = ℝ n A (complemento de A), entonces

No es difícil mostrar que Λ es completamente invariante; es decir, F 1(Λ) = Λ. Por tanto en Λ hay dinámica, posiblemente no trivial. Así que en general se desea conocer tan-to la dinámica de F restricta a Λ como la geometría de este conjunto. Otro asunto de gran interés es conocer la frontera de Λ; esta separaría (en cierto sentido) la cuenca de atracción en 1 con las cuencas de atracción de otros atractores, que obviamente están en Λ; todo esto es de manera general, y muchas veces (dependiendo de cómo es F) sus respuestas o bien son triviales, o en ocasiones difíciles de formular.

En nuestro análisis de la dinámica de F μ,Є como en (1.3), la meta es explicar cómo encontrar un conjunto compacto K2, tal que W = ℝ2 K se comporte como lo discutido en el párrafo anterior. Una aproximación para la búsqueda de este compacto es la siguiente. Sabemos que F μ,Є mapea la fibra {x}X ℝ en la fibra {f μ (x)} X ℝ, la fibra 0 X ℝ es fija, la fibra {1}X ℝ es mapeada en 0 X ℝ y toda fibra {f μ (x)} X ℝ está a la izquierda de 0 X ℝ cualquiera sea x [0;1]. Por otra parte, también es cierto que la imagen por F μ,Є de la horizontal {y = β} es la parábola

Con base a esto afirmamos que existen β1 ≤ 0 y β2 ≥ 1 tales que K = [0;1] [β1 ; β2] es un compacto invariante con órbitas acotadas. Más aún, los valores de β1 y β2 pueden ser determinados de manera que: h μ,Є 1) = h μ,Є 2) y (0, β1) Є P (μ, Є, β1); observe que (0;1) Kc luego P (μ, Є, β1)(0, β1) K c . Una ilustración gráfica de lo afirmado es dado en la Figura 5.

Figura 5 Ilustración de las preimágenes de las horizontales por los valores β1 y β2 deseados. 

5. Sobre las componentes del mapa Skew

En esta sección analizamos la dinámica de las componentes del mapa skew, explorando la existencia de un compacto invariante. Para simplificar las cuentas para cada x Є ℝ hacemos

obviamente F μ,Є (x,y) = (f μ (x); h x (y)), para todo x; y Є ℝ. Resaltamos que para todo n ≥ 1 y x; y 2 ℝ se tiene

donde .

Ahora profundizaremos en el entendimiento de h x , primero notemos que h xx(y) = μ-Є-2μy, por lo que el único punto crítico de y el valor máximo de h x es También es claro que un punto es fijo para h x si y sólo si es raíz de (μ - Є -1)y - μy 2 + Єx = 0. Por tanto h x tiene puntos fijos si y sólo si en tal caso, estos puntos fijos son dados por

obviamente α(x) es la menor raíz de (μ - Є - 1) γ - μγ2 + Єx = 0. También es claro que si x es un punto fijo de f μ, entonces (x;α(x)) es el punto fijo de F μ,Є en la fibra {x} x ℝ con menor ordenada; más precisamente

Denotemos por b(x) la otra preimagen de α(x) según h x ; es decir

en otras palabras

claramente α(x) + b(x) = 2d, α (x) ≤ 0 y b(x) > 0. Por otra parte, es simple verificar que

Así, el gráfico de h x tiene sus puntos fijos como lo indica la Figura 6, es decir: (α(x); a(x)) en el t/ercer cuadrante, (a^(x); a^(x)) en el primer cuadrante y (α(x); b(x) en el cuarto cuadrante.

Figura 6 Gráfica de h x , recordemos que a^ (x) representa un punto fijo de h x  

Note que α(x) siempre es repulsor y a^(x) puede transitar por una bifurcación “flip” en efecto tal bifurcación se caracteriza por la perdida de hiperbolicidad de el punto fijo como consecuencia también se pierde estabilidad y periodicidad del sistema, anticipando la llegada de una una cascada de duplicación de período, ver [13] [Capítulo 3] para más detalles.

Para cada x Є [0;1], el mapa h x se escribe como f ρ (x) , en forma más precisa se tiene

Lema 5.1. Dado x Є [0;1] fijo, el mapa h x es topológicamente conjugado a un elemento f ρ (x) de la familia logística.

Demostración. Para x Є [0; 1] fijo, escribimos , consideremos la transformación afín, con α y β por determinar, su inversa es , hagamos , entonces se cumplen las igualdades

Tomando entonces g se escribe como g(y) = y 2+k, donde , el término sobrante.

Ahora tomamos una segunda transformación afín , tenemos que , obtenemos las igualdades

Haciendo δ2 δ + k = 0, en particular tomando la raíz positiva y γ : = , entonces resulta Eligiendo ρ (x) = , entonces por transitividad de la conjugación, también h x es conjugada a (x) , lo cual concluye la prueba.

Esquemáticamente ilustramos el lema anterior donde la parábola de h x se ve como lo ilustra la Figura 7.

Figura 7 Gráfica de hx, en los nuevos x’ - y’ . 

Así, la conjugación lleva la parábola h x en la logística cuya gráfica es dada en la Figura 8.

Figura 8 Gráfica de f ρ. Recuerde que ρ-1ρ representa un punto fijo de f ρ

Note que los intervalos I x = [a(x); b(x)] son estrictamente crecientes: siempre que x 2 > x 1; esto es, a(x 2) < a(x 1) y b(x 2) > b(x 1) para todo x 2 > x 1; también es claro que

Adicionalmente, como {x} x ℝ es mapeado por , para todo t [0;1], las órbitas acotadas de F μ,Є viven en la banda [0;1] x ℝ; es esta la razón por la que se considera solamente x Є [0;1] al analizar los puntos con órbitas acotadas. Note adicionalmente que si x es cualquier punto en [0;1], entonces en el intervalo I x hay puntos cuya h x -órbita (positiva) se sale de I x si y sólo si μ > 4α (x); esto sigue de la conocida propiedad que en [0;1] hay puntos cuya f ρ -órbita abandona este intervalo si y sólo si ρ > 4. Dicho de otra manera: el valor máximo M x de la función h x excede a b(x) si y sólo si (b(x) a(x)) μ> 4. Como

se tiene que

en particular, M x > b(x) para cada x Є [0;1] si y sólo si μ > 4 + Є.

6. Detección de un buen compacto invariante

Este apartado de nuestro estudio tiene como meta detectar un compacto invariante Q tal que la dinámica del acoplamiento skew F μ,Є fuera de Q sea trivial, es decir todo punto fuera de Q tenga ω-límite vacío este es el resultado principal de nuestra investigación, el cual enunciamos a continuación

Teorema 6.1. El sistema (1.3) admite un conjunto compacto invariante Q tal que su complemento es positivamente invariante y todo punto en Q c tiene conjunto ω -límite vacío.

Demostración. Recordemos que dado un número real β, la imagen de la recta horizontal y = β es la parábola P (μ,Є,β) de ecuación:

Al hacer uso de la función h x (h x (y) = (y) -Єy + Єx), la expresión de arriba es:

Note que cuando esta parábola es la de valores críticos

cuyas ramas inferior y superior son dadas, respectivamente, por

En general, las ramas superior e inferior de la parábola P (μ,Є,β) son dadas por

en particular los valores de y en la expresión anterior cuando x = 0 son

Deseamos un valor de β de forma que h 1(β) = β, así β es a(1) ó a^ (1); tomamos el menor de ellos, α(1), y definimos el compacto Q dado por el rectángulo [0; 1]x[α(1);b(1)]. De lo discutido se desprende lo siguiente:

1.QP(c, µ, α(1)) = QP(c, µ, b(1)) = {(0, α(1))}.

2. Dado β [α (1), b(1)], si (0, y β ) es el punto de intersección entre la rama superior de P(Є, µ, β) y Q. Si β < α(1), entonces y fi < fi; por lo que la ordenada de E µ,c (x, fi), para cualquier x, es menor que β.

Sean de lo anterior se deduce que de hecho, también permite concluir que los puntos en Q - tienen órbita no acotada. En realidad esto es una situación particular de la propiedad general que se enuncia en el siguiente resultado que puede encontrarse en [12], Lema 2.1.

Lema 6.2. Sean E: ℝ n n continua, W n tal que E(W ) W. Si existe una función continua L: W A tal que la diferencia órbitalL = L °F - L sea definida (positiva o negativa), entonces los puntos en W tienen órbita no acotada.

Note que Λ es el conjunto de puntos cuyas órbitas son acotadas; observe también que Λ es completamente invariante; esto es, F μ,Є −1 (Λ) = Λ.

La vecindad Q arriba considerada no es la única con la propiedad descrita. Otra forma de obtener una tal vecindad es la siguiente. Primero considere cualquier es simple verificar que la preimagen, según F μ,Є, de la recta definida por y = x + t es la circunferencia C t (ver Figura 9) cuya ecuación es

Figura 9. Circunferencia C t

Claramente todas estas circunferencias tienen el mismo centro en el punto c = (c 1 ; c 2), donde .

Definamos para cada x; y Є ℝ; note que la curva de nivel t; i.e. g -1(t), es justamente ; además, el único punto crítico de g es el centro de las circunferencias C t , adicionalmente:

así que es una partición de observe que mientras mayor sea el valor del nivel, menor es el radio de la respectiva circunferencia. Por otra parte, el lugar geométrico donde ocurren las tangencias de las circunferencias C t con rectas y = -x + δ es dado por la recta de ecuación y = x - ϵμ; tales puntos le llamaremos puntos de tangencia antidiagonal, los cuales se ilustran en rojo en la Figura 9. En consecuencia los puntos de tangencias antidiagonal de C t son aquellos que resuelven

podemos ver que estos puntos son: ; donde

Con estas observaciones note que la recta antidiagonal tangente a C t en p es la dada por y = x + δ(t) donde

Observe que los puntos fijos de t → δ(t) son aquellos valores t tales que la recta que genera a C t (y = x + t) coincide con la recta antidiagonal tangente (y = x + δ (t)) a C t en el punto p. Es simple mostrar que t → δ (t) tiene a lo más dos puntos fijos y siempre uno de ellos es negativo. Si β es el menor de tales puntos fijos, entonces

Con esta selección del punto fijo β siguen inmediatamente:

  • La recta y = x + β no solo genera a C β como preimagen sino que también es tangente antidiagonal por el respectivo punto p.

  • Para cada t < β, la recta y = x + t (que genera a C t como preimagen) está por debajo de y = x + δ(t) (tangente antidiagonal a C t ) y C β int C t .

  • Mientras que para cada t > β, la recta y = x + t está por encima de y = x + β(t) y en este caso C t int C β.

Para el β seleccionado denotamos por W el semiespacio definido por todos los (x; y) tales que y < x + β. Es simple mostrar que F μ,Є (W ) W y al considerar la función L dada por L(x; y) = x + y se tiene ΔL < 0 en W . De esta forma, al usar el Lema 6.2 se deduce que ext C β  B , y por tanto la validez de (6.6) con Γ = C β. La Figura 10 ilustra lo anterior.

Esto nos da otra construcción de un tal Q y concluye la prueba.

Figura 10 Significado geométrico de δ(t). 

En lo que sigue analizamos algunas propiedades de los conjuntos invariantes construidos.

7. Posiciones relativas

Fijemos Q = [0;1] [a(1); b(1)] como la primera aproximación compacta para describir el conjunto de puntos con órbitas acotadas (complemento de la cuenca B del atractor en el infinito). En adelante este conjunto es denotado por Λμ,Є pues depende de los parámetros Є y μ ya sabemos . Dado que queremos analizar preimágenes de Q, es importante conocer las posiciones relativas de Q y el conjunto de valores críticos J 1: unión del rayo y la parábola crítica . Para esto debe tenerse en cuenta:

1. La ordenada b(1) de los vértices superiores de Q.

2. Las ordenadas de los puntos de intersección de con el eje x = 0, estos son: los mismos corresponden con los valores máximos de las funciones h 0 y h 1 respectivamente.

3. El vértice de la parábola .

Para facilitar, posiblemente, las notaciones y cuentas hagamos

Una cuenta directa muestra, con que:

Una primera observación obvia: cada uno de los segundos factores en las expresiones anteriores son positivas para todo valor μ > 1 y ε Є (0 1); basta observar que δ> 0 y que ; de donde

Es fácil ver que (7.2) se cumple para todo µ ≥ 4 + c y c ∈ (0, 1). Tampoco es difícil verificar que D < 0 para cada c ∈ (0, 1) y todo 1 < µ ≤ 3; así mismo, D < 0 si µ = 4 − c, cualquiera sea 0 < c < 1. Sin embargo, para cada µ ≤ 4 cercano a 4, siempre hay al gún c ∈ (0, 1) de forma que D > 0. Para mostrar esto fijemos cualquier µ > 1, dado que y sigue inmediatamente que si sólo si . Observe que esta desigualdad es válida para µ = 4, luego por continuidad existe δ > 0 tal que para todo 0 < P < δ la desigualdad continua siendo válida para µ = 4 − P. De esta forma, para todo P ∈ (0, δ) y cada ε ∈ (0, 1) se cumple de acá que para cada ρ Є (0,δ) podemos elegir ε ∈ (0; 1) de manera que (7.2) valga para este ∈ y μ = 4-ρ

En cuanto a D 1, no es difícil ver que D 1 < 0 para todo μ≤ 4 + ε y cada ε ∈ (0,1). De hecho, dado que se tiene cualesquiera sean μ > 1 y ε ∈ (0,1); y como la afirmación sigue. También se muestra que D 1 > 0 para todo c (0, 1) y µ 5+c. Además, cuando µ = 5 se tiene que D 1 varía de signo con c; de hecho D 1 > 0 si c es próximo de 0, y es negativo si c es cercano a 1.

Con relación a D 2 dado que

se puede demostrar sin dificultad:

a) Para todo µ ≥ 4 y c ∈ (0, 1), D 2 > 0.

b) Para cada se tiene D 2 < 0.

c) Cuando µ varía en el intervalo cambia de signo.

Las siguientes gráficas en Figura 11 evidencian numéricamente lo anterior.

Figura 11 Posiciones relativas 

8. Comentarios sobre preimágenes

Consideremos cualquier punto (u, v) Єℝ2 que admita preimágen; es decir, que exista al menos un punto (x, y) tal que E µ,c (x, y) = (u, v); esto equivale a f µ (x) = u y k x (y) = v; recuerde que h x : ℝ → ℝ es definida, para cada x ∈ ℝ, fijo por h x (t) = (µЄ)tµt 2 + Єx, cualquiera sea t ∈ ℝ; recuerde además que

Es claro que f μ(x) = u admite solución en x si y sólo si u ≤μ4; en tal caso

definen las soluciones de f μ(x) = u. Por otra parte, para cada uno de los valores u-dependientes de arriba, x_(u) y x +(u), la ecuación h x (y) = v admite soluciones si y sólo si v ≤ max h x ; es decir

en cualquiera de estos casos, las soluciones en y de h x (y) = v son descritas por

Las desigualdades en (8.2) determinan lugares geométricos distinguidos e importantes para el análisis y descripción de la dinámica de F μ,Є. La primera de estas desigualdades corresponde a la clausura de la región abierta Z 4 cuya frontera es la rama inferior de la parábola y el rayo ; la segunda desigualdad define justamente el epigrafo de denotamos por Z 2 el interior de este epigrafo, ver Figura 12.

Analicemos el comportamiento de los puntos en cada una de las zonas arriba descritas.

● Caso A. . Obviamente e inmediato verificar que u . Note que en estas condiciones y claramente , de manera que (8.3) tiene sentido tanto para x = x _(u) como para x = x+(u); por lo cual F -1 μ,Є (u,v) contiene a los más cuatro puntos:

Figura 12 Ilustración del conjunto de valores críticos J 1 = P (μ,Є, 12 (1-ϵμ)) {(μ4,t) :t  μ2+ϵ24μ} junto a las clausuras de las zonas Z 2 y Z 4

Siendo que

Aprovechamos las expresiones en (8.1), (8.4) y (8.5) para hacer algunas precisiones sobre el caso que estamos tratando: A1.A1.

1. Si , entonces los pares ordenados en (8.4) son coincidentes con el punto crítico ; es decir

Por otra parte, si , entonces además de , se tiene

en consecuencia todo punto en el rayo tiene exactamente dos preimágenes, estas se ubican simétricamente respecto del punto crítico . Observe que en cuanto v recorre , las preimágenes de (u;v) recorren toda la línea crítica

2. Si (u;v) está en la rama inferior de , excepto el valor crítico entonces se muestra:

Por tanto F μ,Є -1 (u; v) contiene exactamente tres puntos:

Haciendo uso de las fórmulas en (8.5) y el hecho que (u; v) está en la rama inferior de se deduce por un lado que

y además, al recorrer (u;v) la rama inferior de y el par de puntos recorren, respectiva y unívocamente, el rayo:

y la parábola de ecuación

Ver Figura 13 para una ilustración geométrica.

Figura 13 Ilustración de la ubicación las preimagen de un punto (u,v) en la rama inferior de la parábola Pμ,ϵ, 12 1-ϵμ   

3. Supongamos ahora que (u,v) está en la región abierta Z 4 cuya frontera es la rama inferior de y el rayo . Esto significa que

Atendiendo a estas condiciones y a las fórmulas en (8.1) y (8.5) se tiene

de donde se desprenden las siguientes conclusiones:

A3.1 Los puntos se disponen simétricamente con respecto al rayo crítico .

A3.2 Los puntos se ubican al lado derecho de ι1 y son simétricos con respecto de ι2 pero fuera del epigrafo de la parábola ι2, (8.7).

La ubicación de estos cuatro puntos se distribuye conforme indica la Figura 14.

● Caso B. Esto indica que (u;v) está en el epigrafo de . Dado que hemos analizado el caso en que (u,v) esté en la rama inferior de , supondremos que (u,v) está por encima de esa rama, es decir: ; luego no es difícil ver que:

por lo que de (8.3) sigue que a lo más hay dos puntos (x,y) en F μ,Є -1 (u,v):

B1.B1.

1. Supongamos que el punto (u,v) está sobre la rama superior de , esto es: . En este caso , por tanto F μ,Є -1 (u,v) es un único punto: ; además, al recorrer (u,v) esa ramade recorre unívocamente el rayo .

2. Si (u,v) está en el interior del epigrafo de , entonces

De acá entonces que F μ,Є -1 (u,v) contenga exactamente dos puntos:

también se deduce que estos puntos se encuentran en el interior del epigrafo de la parábolae y se disponen simétricamente con respecto al rayo ; ver Figura 14.

Figura 14 Ubicación de los cuatro puntos en F μ,Є -1 (u; v) cuando (u,v) está en la región abierta Z 4 (zona en verde), y ubicación de las dos preimágenes cuando (u,v) está en Z 2 (zona en gris). 

En este trabajo hemos usado varias técnicas que para determinados parámetros ( ; ) nos permiten detectar un conjunto compacto invariante para (1.3), con lo cuál tenemos herramientas para estudiar zonas ricas en dinámica. En particular, se tienen zonas donde las órbitas son acotadas y otras donde no son acotadas.

Agradecimientos:

Los autores expresan su agradecimiento al prof. Neptalí Romero de la UCLA en Venezuela por su estimulo e interesantes discusiones. También agradecemos a los árbitros anónimos por sus observaciones que ayudaron a mejorar el presente trabajo.

Referencias

[1] L. M. Abramov and V. A. Rohlin, Vestnik Leningrad. Univ. 17(7) (1962), 5-13; English transl., Amer. Math. Soc. Transl. 48(2) (1965), 255-265. [ Links ]

[2] R. L. Adler, A note on the entropy of skew product transformations, Proc. Amer. Math. Soc. 14 (1963), 665-669. [ Links ]

[3] L. Alsedà and M. Misiurewicz, Skew Product Attractors and Concavity. Proc. Amer. Math. Soc. 143 (2015), 703-716. [ Links ]

[4] Hirotada Anzai. Ergodic Skew Product Transformations on the Torus. Osaka Math. J. 3(1), (1951) 83-99. [ Links ]

[5] L. J. Díaz, S. Esteves and J. Rocha. Skew product cycles with rich dynamics: From totally non-hyperbolic dynamics to fully prevalent hyperbolicity. Dyn. Syst. 31(1). (2016) 1-40. [ Links ]

[6] L. S. Efremova. Small C1-smooth perturbations of skew products and the partial integrability property. Appl. Math. Nonlinear Sci. 5(2) (2020) 317-328. [ Links ]

[7] S. Fadaei, G. Keller and F. H. Ghane. Invariant graphs for chaotically driven maps. Nonlinearity. 31. (2018) 5329-5349. [ Links ]

[8] M. Hasler. Synchronization of chaotic systems and transmission of information. Internat. J. Bifur. Chaos Appl. Sci. Engrg. 8(4) (1998) 647-659. [ Links ]

[9] M. Hasler and Y. L. Maistrenko. An introduction to synchronization of chaotic systems: Coupled skew tent maps. IEEE Trans. Circuits Syst. 144 (1997) 856 - 866. [ Links ]

[10] W. Parry. Ergodic properties of affine transformations and flows on nilmanifolds. Amer. J. Math. 91, (1969) 757-771. [ Links ]

[11] Jana Rodríguez Hertz and Raúl Ures. On the Three-Legged Accesibility Property. New Trends in One-Dimensional Dynamics. María José Pacífico and Pablo Guarino (eds.) Springer Pro-ceedings in Mathematics & Statistics. 285 (2019) 239-248. [ Links ]

[12] N. Romero, A. Rovella and R. Vivas. Invariant Manifolds and Synchronization for Coupled Logistic Mappings. Int. J. Pure Appl. Math. Sci. 4(1) (2007) 39-57. [ Links ]

[13] Stephen Wiggins. Introduction to Applied Nonlinear Dynamical Systems and Chaos. Springer, (1990). [ Links ]

[14] D. Lind and B. Marcus (1995) An Introduction to Symbolic Dynamics and Coding. Cambridge University Press, Cambridge. [ Links ]

Para citar este artículo: O. Osuna y C.J. Rojas-Milla, Aspectos dinámicos del acoplamiento skew de la familia logística, Rev. Integr. Temas Mat., 42 (2023), No. 2, 125-145. doi: 10.18273/revint.v41n2-2023004

Recibido: 04 de Julio de 2023; Aprobado: 13 de Diciembre de 2023

Creative Commons License Este es un artículo publicado en acceso abierto bajo una licencia Creative Commons