JOVANA
Library Glossary Getting Started Three Levels Fields How it works Mission
Join the mission
All guides

倍週期與通往混沌之路

一個調節旋鈕緩緩轉動,便能把系統從一個穩態帶到 2-週期、4-週期、8-週期——然後,在無窮多次倍週期被擠進一段有限範圍之後,邁入徹底的混沌。樸實的邏輯斯諦映射,在一條線上展示了這整段旅程。

從流到映射:每繞一圈看一眼

上一篇以一個霍普夫分岔作結,交給我們一個極限環——一個連續流會永遠繞行的持續振盪。為了理解*接下來*會發生什麼,當我們繼續轉動參數時,我們需要一副比完整軌跡更銳利的鏡片。訣竅是:別再連續地盯著迴圈看,而改為每繞一圈拍一張快照。挑一個軌道每一圈都會穿過的曲面,只按順序記下那些穿越點。這份頻閃式的紀錄就是[[poincare-map|龐加萊映射]]:它把 n 維中的連續流,化為 n-1 維中的離散規則,一個把每個穿越點送往下一個的函數。

這是一個貨真價實、而且威力強大的簡化。一個乾淨的極限環,變成映射的一個單一不動點——軌道每一圈都打在同一個穿越點上,於是快照從不移動。一個會晃動的振盪,每一圈回到稍微不同的位置、但兩圈後重複,便表現為一個週期-2 點:映射在兩個穿越點之間來回擺渡。流的長期生命,如今被編碼進一個函數的迭代之中,x -> F(x) -> F(F(x)) -> ...,而問系統會安頓到什麼,就變成問這些迭代落在哪裡。

邏輯斯諦映射:一條包辦一切的線

來認識這個著名的玩具。取一個介於 0 與 1 之間的數 x——把它想成今年的族群數量,佔棲地所能支撐之最大值的比例——並用[[logistic-map|邏輯斯諦映射]]算出明年的值,x_next = r x (1 - x)。那個單一旋鈕 r,落在 0 與 4 之間某處,是增長率。兩個因子在拉鋸:r x 在族群稀少時獎勵增長,而 (1 - x) 在擁擠浮現時加以節制。它看來不起眼——一條拋物線、一次乘法——然而迭代這一條線,就足以走完從平靜到混沌的整條路。

從溫和處起步。當 r 很小、低於 1 時,族群消亡:每個起始值都漂向不動點 x = 0。把 r 推過 1,則 x = 1 - 1/r 處一個新的不動點轉而變得穩定——族群年復一年安頓到單一的穩定水準,一個跨臨界分岔已把穩定性從一個不動點交給了另一個。至此一切尋常:一個參數、一個穩態,正是本級開篇那種分岔。有趣的部分,從這個穩態自身失去掌控時才開始。

一再地倍增

繼續把 r 往上轉。映射的一個不動點恰在該處斜率的大小小於 1 時才穩定——附近的點被拉進來。當 r 攀過 3,族群不動點處的那個斜率穿越 -1:該點變成一個排斥子,系統再也守不住單一的值。但它並不四散瓦解,而是開始*交替*,在兩個值之間高—低—高—低地跳,一個豐年接一個歉年,永遠如此。那個單一的穩態,誕下了一個穩定的週期-2 循環。這是一個倍週期分岔,也叫翻轉分岔,是上一篇霍普夫分岔的離散表親。

如今這道咒語自我重複。當 r 繼續上升,那兩個值各自都是*施行兩次*的映射 F(F(x)) 的不動點,而每一個又以同樣的方式失去穩定性,斜率穿越 -1。這個 2-循環翻轉成一個 4-循環:一種每隔四年才重複的四年模式。再推進,4-循環倍增為 8-循環,然後 16-循環、32-循環。這就是[[period-doubling-cascade|倍週期級聯]]——這條路上具決定性的那一段。系統的固有週期不斷倍增,1、2、4、8、16……每一次倍增,都由斜率在映射一個愈來愈多次迭代的版本中觸及 -1 所觸發。

  r in (1, 3)      ->  fixed point          period 1
  r ~ 3.000        ->  first doubling        period 2
  r ~ 3.449        ->  second doubling       period 4
  r ~ 3.544        ->  third doubling        period 8
  r ~ 3.564        ->  fourth doubling       period 16
     ...               windows shrink fast      ...
  r ~ 3.5699...    ->  accumulation point    chaos begins

  ratio of successive window lengths  ->  4.6692016...  (Feigenbaum delta)
邏輯斯諦映射的級聯:每一次倍增所需的 r 轉動愈來愈小,間距以一個近乎定值的比例縮小,收斂到費根鮑姆常數 4.669……

這裡是令人屏息之處。每一次接續的倍增,所需轉動旋鈕的幅度都比上一次*更小*——r 的窗口縮小,而且以一個近乎定值的比例縮小。一個窗口長度與下一個之比,收斂到一個普適的數,費根鮑姆常數 4.6692……,又因為窗口呈幾何級數縮小,無窮多次倍增被擠進 r 的一段有限範圍裡,全部在約 r = 3.5699 之前完成。過了這個累積點,週期實際上已變為無窮:軌道從不重複。我們經由一道無窮多級的、井然有序的階梯,抵達了混沌。

普適性:到處都是同一個數

一個族群玩具,為何該讓任何研究微分方程的人在意?因為費根鮑姆的發現是:這個數,4.669……,是普適的。它不是 r x (1 - x) 的怪癖。幾乎任何帶有單一光滑隆起的光滑映射——一條拋物線、一個正弦凸包、一百個彼此無關的公式——都會以*同一個*比例,倍週期地走向混沌。而同樣的級聯、同樣的常數,已在真實的連續系統中被測得:滴水的水龍頭、振盪的化學反應、受脅迫的電子電路、對流的流體。那些流的龐加萊映射,是喬裝過的單峰映射,於是它們原封不動地承繼了這條路徑。

對我們這門學科而言,這是一則深刻的訊息。我們在前面各級曾哀嘆,大多數常微分方程沒有封閉形式的解,連定性分析在非線性的世界裡也變得困難。普適性卻獻上一份意外的禮物:在混沌起始的附近,你那個特定方程的*定量*細節不再重要。級聯的形狀與常數 4.669……,僅由折疊與拉伸的幾何所決定,而非由你從哪一個物理系統出發所決定。一整類系統共享同一枚指紋,而你不必解出任何一個方程,便能辨認出通往混沌之路。

混沌是什麼——又不是什麼

過了累積點,把 x 的長期值對 r 畫出來,你便得到那張著名的[[bifurcation-diagram|分岔圖]]:一條線裂成兩條、再成四條、再成八條,然後模糊成密集陰影的混沌帶——卻又被一道道突然明亮、秩序井然的窗口貫穿,在那裡一個週期-3 的循環陡然重現,並再次從頭倍增。在混沌帶之內,軌道是[[deterministic-chaos|決定論混沌]],而「決定論」一詞要照字面理解:x_next = r x (1 - x) 裡沒有任何隨機性。只要給定確切的 x,整個未來就被固定了。沒有骰子、沒有雜訊、沒有任何隱藏之物。

然而它卻無法預測。調和這矛盾的關鍵是[[sensitive-dependence|對初始條件的敏感依賴]]:兩個相隔僅一髮之微的起始值,會隨著每次迭代被拉開,它們的間距平均每一步乘上一個固定的因子。這分離的速率,就是[[lyapunov-exponent|李雅普諾夫指數]];混沌恰恰是它為正的那種情形,於是任何誤差都呈指數增長。第十五位小數的捨入、那隻諺語裡的蝴蝶振翅——在不太多的步數之內,它便膨脹到主宰了答案。決定論在原則上固定了未來;敏感依賴卻使它在實務上無法得知,因為你永遠無法以無窮的精度指定當下。