Dinamik sistemler · x → r x (1 − x)

Lojistik dönüşüm

Lojistik dönüşüm tek satırlık x → r x (1 − x) kuralıdır: r arttıkça uzun vadeli davranış tek değerden ikiye, dörde ve kaosa geçer.

Bu görselleştirmede lojistik dönüşüm

Lojistik dönüşüm, 0≤r≤40 \le r \le 4 parametresiyle [0,1][0, 1] aralığını kendisine götüren f(x)=rx(1−x)f(x) = r x (1 - x) dönüşümüdür. x0x_0 başlangıcından xn+1=f(xn)x_{n+1} = f(x_n) işlemini tekrar tekrar uygulamak bir dizi verir; bu sayfa dizinin uzun vadeli davranışının rr ile nasıl değiştiğini inceler.

xn+1=r xn (1−xn),xn∈[0,1],0≤r≤4x_{n+1} = r\,x_n\,(1 - x_n),\qquad x_n \in [0, 1],\qquad 0 \le r \le 4

Görseller ne gösteriyor?

Çatallanma diyagramı, büyük görsel: yatayda rr, düşeyde xx. Her sütun merkezindeki rr değerini alır, x0=1/πx_0 = 1/\pi ile başlar, ilk TT adımı atlar ve sonraki KK adımın nereye düştüğünü kaydeder; daha parlak, daha sık demektir. TT ve KK sağ üst köşede yazılır ve yakınlaştıkça büyür. Dolayısıyla sütun yörüngenin tüm geçmişi değil, tek başlangıcın uzun vadeli örneklemidir. “Tipik”, (uzunluk anlamında) hemen hemen her başlangıç demektir; 1/π1/\pi onları temsil eder, ancak her rr için tipik olduğunun kanıtı yoktur.

Diyagramın altında aynı rr eksenindeki λ şeridi, λ(r)=lim⁡1n∑ln⁡∣f′(xi)∣\lambda(r) = \lim \frac1n \sum \ln|f'(x_i)| Lyapunov üssünün sayısal tahminidir (aynı tipik başlangıç, 1000 adım atlanır, 4000 adımın ortalaması alınır). İki başlangıç arasındaki sonsuz küçük farkın ortalama üstel büyüme hızıdır: 0 altında farklar ortalamada küçülür, 0 üstünde büyür. Süper kararlı parametrede (12\tfrac12 çekici döngü üzerindedir) λ=−∞\lambda = -\infty olur; sayısal bulunan parametrede periyot etiketli aşağı okla çizilir.

Panelin başındaki örümcek ağı grafiği: y=f(x)y = f(x) eğrisi (veya f2f^2, f3f^3, f4f^4) ve y=xy = x köşegeni. x0x_0 noktasından eğriye yukarı, köşegene yatay: bir adım. Dolu noktalar çekici, boş noktalar itici sabit noktalardır.

Tek tuşla açılan zaman serisi: nn'ye göre xnx_n; sağında, bu yörüngenin 10510^5 adımının 40 kutunun her birindeki payı. r=4r = 4 için koyu çizgi kuramsal olasılıktır.

Neden böyle görünüyor?

x∗=1−1/rx^* = 1 - 1/r sabit noktasında eğim f′(x∗)=2−rf'(x^*) = 2 - r olur. 1<r<31 < r < 3 için ∣f′(x∗)∣<1|f'(x^*)| < 1 sağlanır ve komşularını çeker. rr, 3'ü geçince f′(x∗)f'(x^*), −1-1 değerini geçer: çekicilik kaybolur; yanında 4+2r−r24 + 2r - r^2 çarpanlı r+1±(r−3)(r+1)2r\frac{r + 1 \pm \sqrt{(r - 3)(r + 1)}}{2r} 2 periyotlu döngü belirir. Çift, f2f^2 eğrisinin köşegenle yeni kesiştiği yerdir; görmek için örümcek ağını f2f^2 yapın. r=1+6r = 1 + \sqrt6 değerinde 2 periyotlu döngünün çarpanı −1-1 olur ve o da bölünür.

Katlanmalar giderek hızlanır ve r∞≈3.5699r_\infty \approx 3.5699 noktasında birikir. Ardışık aralıkların oranı Feigenbaum sabiti δ≈4.6692\delta \approx 4.6692 değerine yaklaşır; yalnızca bu dönüşümde değil, ikinci dereceden maksimuma sahip bütün bir tek tepeli dönüşüm sınıfında aynıdır (Feigenbaum 1978'de buldu; sonra Lanford bilgisayar desteğiyle kanıtladı).

Kaos içindeki pencere: rr, 1+8≈3.82841 + \sqrt8 \approx 3.8284 değerini geçince f3f^3 köşegene değer; biri çekici, diğeri itici iki 3 periyotlu döngü doğar (eyer-düğüm çatallanması, çarpan +1+1). Pencerenin içinde kendi periyot katlanması dizisi bulunur.

r=4r = 4 için f(sin⁡2πu)=sin⁡2(2πu)f(\sin^2 \pi u) = \sin^2(2\pi u): x=sin⁡2(πu)x = \sin^2(\pi u) koymak, u↦2u mod 1u \mapsto 2u \bmod 1 ikiye katlama dönüşümünü ff üzerine taşır (yarı eşleniklik: uu ve 1−u1 - u aynı xx değerini verir). Buradan tipik başlangıcın uzun vadeli yoğunluğu 1πx(1−x)\frac{1}{\pi\sqrt{x(1 - x)}} ve λ=ln⁡2\lambda = \ln 2 tam olarak elde edilir.

Kaotik görünen bölgenin her yanında periyodik pencereler bulunur. Çekici döngü varsa tipik başlangıçlar ona yaklaşır; yoksa çizilen sütun tüm başlangıçların sınıflandırması değil, belirtilen başlangıç değerinin uzun vadeli örneklemidir.

Yanlış yorumlamaya açık noktalar

Tarihçe

“Lojistik” adı Verhulst'un 1838 nüfus artışı diferansiyel denkleminden gelir. May 1976'da Nature'da bu ayrık dönüşümü kullanarak basit modellerin çok karmaşık davranabileceğini gösterdi [1]. Ulam ve von Neumann r=4r = 4 için değişmez yoğunluğu 1947'de zaten belirtmişti [2]. Feigenbaum 1978'de periyot katlanması dizisinin evrensel sabitlerini buldu [3]; Singer en fazla bir çekici döngü için ölçüt verdi [4].

İlgili kavramlar

Kaynaklar

  1. R. M. May, “Simple mathematical models with very complicated dynamics”, Nature 261 (1976) 459–467. doi:10.1038/261459a0
  2. S. M. Ulam, J. von Neumann, “On combination of stochastic and deterministic processes” (abstract), Bulletin of the American Mathematical Society 53 (1947) 1120.
  3. M. J. Feigenbaum, “Quantitative universality for a class of nonlinear transformations”, Journal of Statistical Physics 19 (1978) 25–52. doi:10.1007/BF01020332
  4. D. Singer, “Stable orbits and bifurcation of maps of the interval”, SIAM Journal on Applied Mathematics 35 (1978) 260–267.
  5. O. E. Lanford III, “A computer-assisted proof of the Feigenbaum conjectures”, Bulletin of the American Mathematical Society (N.S.) 6 (1982) 427–434.
  6. M. V. Jakobson, “Absolutely continuous invariant measures for one-parameter families of one-dimensional maps”, Communications in Mathematical Physics 81 (1981) 39–88.
  7. J. Graczyk, G. Świątek, “Generic hyperbolicity in the logistic family”, Annals of Mathematics 146 (1997) 1–52.
  8. M. Lyubich, “Dynamics of quadratic polynomials, I–II”, Acta Mathematica 178 (1997) 185–297.

Ek okumalar: Vikipedi: Logistic map (İngilizce); MacTutor Matematik Tarihi: Robert May (İngilizce); MacTutor Matematik Tarihi: Mitchell Feigenbaum (İngilizce).