Logistic映射

**單峰映射**（英語：Logistic map）是種二次多項式的映射（遞迴關係式），是一個由簡單非線性方程式式產生混沌現象的經典範例。這種映射因生物學家 Robert May 在 1976 年發表的一篇論文而著名，一定程度上是離散時間的族群 / 人口模型，類似於皮埃爾・弗朗索瓦・韋呂勒的邏輯斯諦函數。。單峰映射實質上是邏輯斯諦函數的差分方程式，其數學表達為：

　
: $x_{n+1}=rx_{n}(1-x_{n})$

其中

 * $x_{n}$ 是介於 0 和 1 之間的數，表示當前族群數量 / 人口數量與環境承載力的比值。
* $r$ 是正整數，是根據繁殖和死亡率而得出的數。

單峰映射的方程式旨在描述以下兩個現象：

1. 當族群數量 / 人口少時，繁殖增加的個體數大致跟族群原本的總數目成正比；
2. 高密度導致的死亡，環境資源有其承載力（最大容量），當種接近最大容量時，增長率下降的速度與環境承載力減去當前族群數量的差成正比。

參數 r 通常取[0, 4]區間內的值，因此 xn 在[0, 1]上保持有界。_r_ = 4 的情景是 dyadic 映射及參數 μ=2 的帳篷映射的非線性變換。當 _r_ > 4 時，族群數量會出現負值（該問題在更早前的同樣表現出混沌動態的 Ricker 模型中不會出現）。也可在[−2, 0]的區間內取 r 值，該情形下 xn 有界，處於[−0.5, 1.5]之間。

## r 的值對結果的影響

變化參數 $r$ 的值，其結果如下：

* 0 和 1 之間：不論起始值數值為何， $x_{n}$ 會越來越少，最後趨近於 0。
* 1 和 2 之間：不論起始值數值為何， $x_{n}$ 會快速的趨近 ${\frac {r-1}{r}}$ 。
* 2 和 3 之間：經過幾次迭代， $x_{n}$ 也會越來越接近 ${\frac {r-1}{r}}$ ，但一開始會在這個值左右振動，而收斂速率是線性的。
* 3： $x_{n}$ 仍然會越來越接近 ${\frac {r-1}{r}}$ ，但收斂速率極為緩慢，不是線性的。
* 3 和 $1+{\sqrt {6}}$ （約 3.45）之間：針對幾乎所有的初值， $x_{n}$ 最後會在 2 個值之間持續的震盪，即 $x_{n}$ 最後會是 a,b,a,b... 的變化，其數值和 $r$ 有關。
* 3.45 和大約 3.54 之間，針對幾乎所有的初值， $x_{n}$ 最後會在 4 個值之間持續的震盪。
* 約大於 3.54： $x_{n}$ 最後會在 8 個、16 個、32 個值…… 之間持續的震盪，至於 $r$ 何時會令 $x_{n}$ 的值由 n 個到 2n 個，則和費根鮑姆常數 $\delta =4.669...$ 有關。
* 約為 3.5699：這樣的振動消失，整個系統開始在混沌的狀態之中。針對幾乎所有的初值，都不會出現固定週期的震盪，初值再微小的變化，隨著時間都會使結果產生明顯的差異，這就是典型混沌的特性。
* 大於 3.5699：整個系統在混沌的狀態之中。不過，當中有些特定的 $r$ 值還是使系統變成非混沌，有週期性的結果，這些區間稱為「穩定島」。例如當 $r$ 約大於 3.82，會出現 3 個值的週期，再大一點出現 6 個值及 12 個值的週期。
* 當 $r$ 從大約 3.5699 到大約 3.8284 之間，系統混沌性質發展的現象有時會稱為 Pomeau–Manneville 場景，其特徵是週期性的震盪和非週期性的行為會穿插出現。此特徵可以應用在半導體元件中。也有其他區域會使系統的週期為 5 個值，不管任意週期都存在某特定的 $r$ ，使週期為指定值。
* 大於 4：針對幾乎所有的初值，系統最後都會超過區間 [0,1] 並且發散。

對於任一個 $r$ 值，最多只有一個穩定的極限環，若穩定極限環存在，幾乎所有的點最後都會趨近極限。若某一個 $r$ 值有一個穩定的極限環，可能也有無限個不同週期的不穩定極限。

這些情況可用分枝圖表示，分枝圖中的橫軸是參數 $r$ 的數值，縱軸中顯示大部份初值下，穩態可能的 $x$ 值，若最後數值會在 2 個值中震盪，分枝圖上對應的數值就會有 2 點。若某特定 $r$ 值下，已無法明確的看到有幾個對應的點，當時系統可能已經在混沌狀態下。

分枝圖有自相似的特性。若將分枝圖中 $r=3.82$ 的部份展開，只取三個分支中的一個。其圖形會好像是原分枝圖縮放及扭曲後的結果。針對所有非混沌的參數 $r$ 都有此一特性。以此可以看出混沌和碎形深入及明顯的關係。

## 單峰映射和混沌

和其他混沌系統比較，單峰映射較為簡單，是一個說明混沌特性的很好的例子。簡單來說，混沌就是對初始條件的高度靈敏度。 $r$ 是在 3.57 及 4 之間的大部份數值都可以使單峰映射出現此一特性。對初始條件有高度靈敏度的常見原因是映射本身是對定義域的拉伸及摺疊。單峰映射的二次差分方程式可視為是對於區間 (0,1) 拉伸及摺疊的過程。

利用二維及三維的相圖可以看出一些單峰映射的特性。以 $r=4$ 的單峰映射為例，二維相圖為一拋物線，但是若用 $(x_{n},x_{n+1},x_{n+2})$ 繪製三維相圖，可看出進一步的結構，例如幾個一開始很接近的點在迭代後開始發散．特別是位在斜率較大位置的點。

拉伸及摺疊的結果使迭代的數列以指數形式發散（參照李亞普諾夫指數），可以用有混沌特性時，單峰映射的複雜及不可預測性說明這一點。事實上，數列的指數發散說明了混沌和不可預測性之間的關係：初值微小的誤差在迭代過程中會以指數成長的方式增加。因此當對於初始狀態的資訊中有微小的誤差時．對未來狀態的預測準確度也會隨迭代次數增加而快速變差。

由於映象是限制在實數數線的一段範圍內，其維度小於或等於 1。依數值分析的結果，在 r=3.5699456... 時（剛開始混沌特性時），其關聯維度為 0.500 ± 0.005（Peter Grassberger,1983）、郝斯多夫維數大約是 0.538，而碎形維數為 0.5170976...。

有些混沌系統可對於其未來狀態的可能性作準確的描述。若一個可能有混沌特性的動態系統存在吸引子，則存在一機率量測描述系統長期下，在吸引子各部份所花時間的比例。以 $r=4$ 的單峰映射而言，起始狀態在區間 (0,1) 中，而吸引子也在間 (0,1) 中，其機率量測對應參數 $a=0.5,b=0.5$ 的 Β 分布，其不變測度為 $\pi ^{-1}x^{-1/2}(1-x)^{-1/2}$ 。不可預期性和隨機不同，不過在一些情形下這二很類似，而且即使對單峰映射（或其他混沌系統）初值只有很少的資訊，仍然可以針對長期的分布作某種程度的預測。

## 部份情形下的解

在 $r=4$ 及 $r=2$ 的特殊情形下，單峰映射有解析解。不過大部份情形下的通解只能以統計的方式預測。 $r=4$ 的解為

　
: $x_{n+1}=\sin ^{2}(2^{n}\theta \pi )$

其中初始條件參數 $\theta $ 是由 $\theta ={\tfrac {1}{\pi }}\sin ^{-1}(x_{0}^{1/2})$ 求得。針對有理數的 $\theta $ ，有限次數的迭代後 $x_{n}$ 就會變成一個週期性的數列。不過幾乎所有的 $\theta $ 都是無理數，此時 $x_{n}$ 不會重複，因此沒有週期解。此解可以清楚的看出混沌的二個重要特徵：拉伸及摺疊。係數 2n 表示拉伸的指數成長，因此造成蝴蝶效應，也就是對初始值的高度靈敏性，而解中包括正弦函數的平方，使解摺疊在 [0, 1] 的範圍內。

 $r=2$ 的解為：

$$
x_{n}={\frac {1}{2}}-{\frac {1}{2}}(1-2x_{0})^{2^{n}}
$$

對於 $x_{0}\in [0,1)$ 。此解沒有混沌的特性。由於對於不包括不穩定固定點 0 在內的 $x_{0}$ ，當 _n_ 趨近無限大時 $(1-2x_{0})^{2^{n}}$ 會趨近於零，因此 $x_{n}$ 會趨近穩定的固定點 ${\tfrac {1}{2}}.$ 。

## _r_ = 4 時找任意週期的循環

_r_ = 4 時，幾乎所有的初值都會使單峰映射出現混沌特性，不過也存在無限個初值會使單峰映射最後呈週期性變化，而且對於所有正整數，都存在一初值使單峰映射的週期為該正整數。可以利用單峰映射和位元位移映射之間的關係來找出任何週期的循環。若 _x_ 依照單峰映射 $x_{n+1}=4x_{n}(1-x_{n})\,$ 而 _y_ 依照位元位移映射

　
: $y_{n+1}={\begin{cases}2y_{n}&0\leq y_{n}<0.5\\2y_{n}-1&0.5\leq y_{n}<1,\end{cases}}$

則二個變數的關係如下：

　
: $x_{n}=\sin ^{2}(2\pi y_{n})$ .

位元位移映射其名稱是因為當 _y_ 以二進制表示時，映射會將二進制的數字左移一位。例如若數字是二進制的循環小數，循環節為 001，則位元位移映射的序列為 001001001... →010010010... →100100100... →001001001...，為週期為 3 的循環，循環節為 010, 011, 100, 101, 110 時也會有類似情形，這些循環小數都可以轉換為對應的分數，上例若以分數表示為：1/7 → 2/7 → 4/7 → 1/7。轉換到 r=4 的單峰映射後，為 611260467... → .950484434... → .188255099... → .611260467...。其他週期為 3 的循環也可以轉換為單峰映射。依相同方式也可以找出在
位元位移映射下，任意週期的循環，再轉換為單峰映射。

不過幾乎所有在區間 [0, 1) 的數字都是無理數，而初始值為無理數的位元位移映射沒有循環的特性，因此對應的單峰映射也沒有循環的特性。

## 相關條目

* 混沌理論
* 曼德博集合，邏輯斯蒂映射可以看作是其圖像的一部分。
* 剛性方程式
* 李雅普諾夫穩定性

## 參考資料

### 教科書

* Sprott, Julien Clinton. Chaos and Time-Series Analysis. Oxford University Press. 2003. ISBN 0-19-850840-9.
* Strogatz, Steven. Nonlinear Dynamics and Chaos. Perseus Publishing. 2000. ISBN 0-7382-0453-6.
* Tufillaro, Nicholas; Tyler Abbott, Jeremiah Reilly. An experimental approach to nonlinear dynamics and chaos. Addison-Wesley New York. 1992. ISBN 0-201-55441-0.

## 外部連結

* 曼德博集合和單峰映射的聯繫
* The Chaos Hypertextbook（頁面存檔備份，存於網際網路檔案館）
* Java Applet
* Logistic Map Simulation（頁面存檔備份，存於網際網路檔案館）. A Java applet simulating the Logistic Map by Yuval Baror.
* Logistic Map. Contains an interactive computer simulation of the logistic map.
* The Chaos Hypertextbook（頁面存檔備份，存於網際網路檔案館）. An introductory primer on chaos and fractals.
* Interactive Logistic map with iteration and bifurcation diagrams in Java.
* Interactive Logistic map showing fixed points.
* Macintosh Quadratic Map Program （頁面存檔備份，存於網際網路檔案館）
* The transition to Chaos and the Feigenbaum constant（頁面存檔備份，存於網際網路檔案館）- JAVA applet
* The Logistic Map and Chaos（頁面存檔備份，存於網際網路檔案館） by Elmer G. Wiens
* Complexity & Chaos (audiobook) by Roger White. Chapter 5 covers the Logistic Equation.
* "History of iterated maps（頁面存檔備份，存於網際網路檔案館）," in _一種新科學_ by 史蒂芬・沃爾夫勒姆. Champaign, IL: Wolfram Media, p. 918, 2002.
* Discrete Logistic Equation（頁面存檔備份，存於網際網路檔案館） by Marek Bodnar after work by Phil Ramsden, Wolfram 演示項目.
* Multiplicative coupling of 2 logistic maps（頁面存檔備份，存於網際網路檔案館） by C. Pellicer-Lostao and R. Lopez-Ruiz after work by Ed Pegg Jr, Wolfram 演示項目.
* Using SAGE to investigate the discrete logistic equation（頁面存檔備份，存於網際網路檔案館）