Featured image of post 碎形維度與豪斯多夫測度:超越整數壁壘的小數維度與自我相似性的科學

碎形維度與豪斯多夫測度:超越整數壁壘的小數維度與自我相似性的科學

曼德布洛特集合、海岸線悖論,以及介於一維與二維之間的分數維度。幾何學與測度論所揭示的自然秩序。

前言:重新思考維度的概念

我們日常經驗中的空間,通常被認知為三維的歐幾里得空間。紙上的線是一維的,平面是二維的,而立體是三維的。這是自古希臘歐幾里得以來,數千年間作為人類空間認知基礎的堅定直覺。然而,當我們觀察自然界中複雜的形狀時,這種「整數維度」的典範便面臨了決定性的極限。雲朵不是球體,山脈不是圓錐,海岸線不是圓弧。樹木的分支、血管的網絡、閃電的軌跡等,自然界中可見的許多形態,都具有本質上不同於平滑歐幾里得幾何學對象的「粗糙度(roughness)」。

為了解決以數學來描述這種自然界複雜性的問題,「碎形幾何學(Fractal Geometry)」作為一種新語言誕生了。支撐其理論根基的,是源自實分析與測度論深淵的概念:「豪斯多夫測度(Hausdorff measure)」以及隨之而來的「豪斯多夫維度(Hausdorff dimension)」。本文將從直覺的幾何學到嚴謹的測度論,徹底解說碎形維度是如何被定義、計算,並應用於理解自然現象。


第一章:歐幾里得幾何學的極限與「自然的粗糙度」

本華·曼德布洛特的提問『英國的海岸線有多長?』

有一個象徵著碎形幾何學開端的著名問題。「英國的海岸線有多長?(How Long Is the Coast of Britain?)」這是本華·曼德布洛特(Benoit Mandelbrot)於1967年發表在科學期刊《Science》上的論文標題。

乍看之下,這個問題似乎只是個測量問題。但其中卻潛藏著深刻的悖論。假設為了測量海岸線的長度,我們使用一把固定長度的尺(例如長度 $\eta = 100 \text{ km}$)來近似海岸線。如果將尺的長度變小($\eta = 10 \text{ km}, 1 \text{ km}, 1 \text{ m}, \dots$),測量出的總長度會如何變化呢?如果是平滑的曲線(如圓或拋物線等),隨著尺變小,長度會收斂到一個固定的有限值。這是古典曲線長度(弧長)的定義。

然而,在真實的海岸線上卻並非如此。尺越小,那些原本隱藏在尺之間被忽略的小海灣和岩石的凹凸就會被測量進去,使得海岸線的長度無限增加。也就是說,在測量尺度 $\eta$ 趨近於 $0$ 的極限下,海岸線長度 $L(\eta)$ 會發散至無限大。

理查森效應

這個現象由氣象學家路易斯·弗萊·理查森(Lewis Fry Richardson)憑經驗發現。理查森測量了各國國界線和海岸線在不同尺度下的長度,發現測量尺度 $\eta$ 與測得長度 $L(\eta)$ 之間存在著以下的冪定律:

$$ L(\eta) \propto \eta^{1-D} $$

在這裡,$D$ 是一個常數,海岸線越複雜,$D$ 的值就越大。理查森本人將這個 $D$ 視為經驗常數,但曼德布洛特給予了它深遠的數學解釋。也就是說,他認為這個 $D$ 正好代表了目標物的「維度」。

在平滑一維曲線的情況下,$D=1$,因此 $L(\eta) \propto \eta^0 = 1$,長度會收斂至固定值。但是,像英國海岸線這種極度複雜的邊界,$D \approx 1.25$,因此 $1 - D = -0.25 < 0$,隨著 $\eta \to 0$,$L(\eta) \to \infty$。這個大於 $1$ 且小於 $2$ 的實數維度,正是「碎形維度」最初的萌芽。


第二章:自我相似性與相似維度

碎形(fractal)一詞源自拉丁語「fractus(破碎的、斷裂的)」,由曼德布洛特創造。碎形最基本的特徵之一就是「自我相似性(self-similarity)」。這指的是將整體放大後,其中包含了與整體相同結構的性質。

利用這種自我相似性,可以導出「相似維度(Similarity Dimension)」這種直覺的維度定義。

相似維度 $D$ 的直覺推導

讓我們思考一下平滑歐幾里得圖形的性質。

  • 將一維的線段縮小為 $1/r$ 的尺度時,要構成原本的線段,需要 $r^1$ 個縮小的線段。
  • 將二維的正方形各邊縮小為 $1/r$ 時,要構成原本的正方形,需要 $r^2$ 個小正方形。
  • 將三維的立方體各邊縮小為 $1/r$ 時,要構成原本的立方體,需要 $r^3$ 個小立方體。
$$ N = r^d $$$$ \log N = d \log r $$

解出維度 $d$,可得以下定義:

$$ d = \frac{\log N}{\log r} $$

將此定義擴展至不具備整數維度的自我相似圖形,就是「相似維度」。

$$ D_s = \frac{\log N}{\log(1/r)} $$

這裡的 $r$ 是縮小率($0 < r < 1$),$N$ 是使用該縮小圖形來完全覆蓋原圖形所需的數量。(當 $r$ 作為縮小率時,分母會是 $\log(1/r)$。在前面的例子中 $r$ 是倍率,因此請注意符號的定義)。

康托爾集合(Cantor Set)

由格奧爾格·康托爾(Georg Cantor)於1883年引入的這集合,是測度論中最重要的反例之一。 其建構方法如下:

  1. 從區間 $[0, 1]$ 開始(步驟0)。
  2. 移除中間的 1/3,即 $(1/3, 2/3)$(步驟1:區間變為 $[0, 1/3] \cup [2/3, 1]$)。
  3. 移除剩下每個區間中間的 1/3。
  4. 無限重複此過程。
$$ D = \frac{\log 2}{\log 3} \approx 0.6309 $$

這是一個大於0維(點)且小於1維(線)的集合。令人驚訝的是,這個集合的勒貝格測度(長度)為 $0$,卻包含了不可數無限個點。

科赫曲線(Koch Curve)

這是由海里格·馮·科赫(Helge von Koch)於1904年提出,連續但在任何地方皆不可微的曲線。

  1. 將線段三等分。
  2. 將中間的區間替換為以該區間為底邊的正三角形的兩邊。
  3. 對所有的線段重複此過程。
$$ D = \frac{\log 4}{\log 3} \approx 1.2618 $$

謝爾賓斯基墊片(Sierpinski Gasket)

$$ D = \frac{\log 3}{\log 2} \approx 1.5849 $$

面積(二維勒貝格測度)為0,但一維長度卻是無限大。


第三章:豪斯多夫外測度與豪斯多夫維度的嚴謹定義

相似維度很直覺且易於計算,但只能適用於具備「嚴格自我相似性」的圖形。為了決定自然界的碎形或數學上複雜集合(失去自我相似性的集合)的維度,需要一個基於實分析與測度論、嚴謹且普遍的維度定義。這就是「豪斯多夫維度(Hausdorff Dimension)」。

1918年,費利克斯·豪斯多夫(Felix Hausdorff)擴展了卡拉西奧多里(Carathéodory)的測度論方法,為任意非負實數 $d$ 定義了 $d$-維外測度。

$\delta$-覆蓋($\delta$-cover)

$$ E \subset \bigcup_{i=1}^\infty U_i \quad \text{且} \quad \operatorname{diam}(U_i) \leq \delta $$

則稱之為 $E$ 的 $\delta$-覆蓋。在此,$\operatorname{diam}(U_i)$ 是 $U_i$ 的直徑(上確界距離)$\sup_{x,y \in U_i} \|x - y\|$。

豪斯多夫外測度 $\mathcal{H}^d(E)$

固定一個非負實數 $d \geq 0$。對於 $E$ 的任意 $\delta$-覆蓋 $\{U_i\}$,考慮其各自直徑的 $d$ 次方之和,並取其下確界。

$$ \mathcal{H}_\delta^d(E) = \inf \left\{ \sum_{i=1}^\infty (\operatorname{diam} U_i)^d \mathrel{\Big|} \{U_i\} \text{ 是 } E \text{ 的 } \delta\text{-覆蓋} \right\} $$

當 $\delta$ 逐漸變小,覆蓋的條件變嚴格,因此下確界的集合變窄,$\mathcal{H}_\delta^d(E)$ 呈現單調遞增(非遞減)。因此,$\delta \to 0$ 的極限存在(包含 $\infty$)。

$$ \mathcal{H}^d(E) = \lim_{\delta \to 0} \mathcal{H}_\delta^d(E) = \sup_{\delta > 0} \mathcal{H}_\delta^d(E) $$

這個 $\mathcal{H}^d(E)$ 被稱為 $d$ 維的豪斯多夫測度。從測度論的角度來說,這是一個具備波萊爾正則性(Borel regularity)的外測度(滿足卡拉西奧多里條件),並成為波萊爾集族上滿足可數可加性的真實測度。

在整數維度 $d = n$ 的情況下,$\mathcal{H}^n(E)$ 與一般的 $n$ 維勒貝格測度只有常數倍的差異(若調整正規化常數則會完全一致)。

作為跳躍臨界值的豪斯多夫維度 $\dim_H(E)$

$$ \sum (\operatorname{diam} U_i)^s = \sum (\operatorname{diam} U_i)^{s-d} (\operatorname{diam} U_i)^d \leq \delta^{s-d} \sum (\operatorname{diam} U_i)^d $$

當 $\delta \to 0$ 時,$\delta^{s-d} \to 0$,因此 $\mathcal{H}^s(E) = 0$。 反過來說,若 $\mathcal{H}^s(E) > 0$,則對於任意 $d < s$,$\mathcal{H}^d(E) = \infty$。

也就是說,當 $d$ 從 $0$ 開始增加時,$\mathcal{H}^d(E)$ 在某個臨界點之前會始終為 $\infty$,一旦超過那個臨界點,就會始終為 $0$,這展現出一種極端的「跳躍」。這個臨界的 $d$ 值就被定義為 豪斯多夫維度(Hausdorff Dimension)。

$$ \dim_H(E) = \inf \{ d \geq 0 \mid \mathcal{H}^d(E) = 0 \} = \sup \{ d \geq 0 \mid \mathcal{H}^d(E) = \infty \} $$

這個定義有著令人驚嘆的優雅:即使目標集合 $E$ 不具備自我相似性,甚至是任何病態的集合,都能為所有度量空間的子集嚴格且唯一地確定其維度。康托爾集合和科赫曲線的豪斯多夫維度,與前述的相似維度完全一致。


第四章:計盒維度(容量維度)、資訊維度與堆積維度

豪斯多夫維度在數學上是最洗鍊的概念,但不適合數值計算或實驗數據分析(因為必須從無限的覆蓋模式中求出下確界並取極限)。因此,在應用數學與物理學中,使用了更容易計算的碎形維度定義。

計盒維度(容量維度, Box-counting Dimension)

用邊長為 $\varepsilon$ 的網格(格子)來分割空間,將與目標集合 $E$ 相交的盒子(網格點)數量記為 $N(\varepsilon)$。此時,計盒維度 $\dim_B(E)$ 的定義如下:

$$ \dim_B(E) = \lim_{\varepsilon \to 0} \frac{\log N(\varepsilon)}{-\log \varepsilon} $$

這個定義極具實用性,也是在影像分析等領域估算碎形維度演算法(覆蓋法)的基礎。不過,它也有數學上的缺點。例如,有理數集合 $\mathbb{Q} \cap [0,1]$ 的計盒維度是 $1$,但因為是可數集合,其豪斯多夫維度為 $0$。一般而言,$\dim_H(E) \leq \dim_B(E)$ 必定成立。

資訊維度(Information Dimension)與廣義維度

當碎形集合具有不均勻分布時,單純計算盒子數量是不夠的。假設每個盒子 $i$ 所包含的測度(機率)為 $P_i$,利用夏農熵(Shannon entropy) $I(\varepsilon) = - \sum P_i \log P_i$ ,可以定義資訊維度 $D_1$。

$$ D_1 = \lim_{\varepsilon \to 0} \frac{\sum P_i \log P_i}{\log \varepsilon} $$

此外,在基於阿爾弗雷德·雷尼(Alfréd Rényi)廣義熵的「多重碎形(Multifractal)」理論中,進一步發展出了廣義維度(Rényi維度) $D_q$ 的概念。

堆積維度(Packing Dimension)

由 Tricot 在1980年代引入的堆積維度 $\dim_P(E)$,是與豪斯多夫維度對偶的概念。豪斯多夫維度採取「覆蓋集合」的途徑,而堆積維度則採取「在集合內填入(堆積)球體」的途徑。 嚴格來說,$\dim_H(E) \leq \dim_P(E) \leq \dim_{\overline{B}}(E)$(上計盒維度)的關係成立,在機率論的集合分析中是非常強大的工具。


第五章:曼德布洛特集合與朱利亞集合的複動力系統

談論碎形幾何學時,絕對無法避開複動力系統(Complex Dynamics)的世界。特別是由複平面上極為簡單的二次映射所產生的「曼德布洛特集合(Mandelbrot set)」,被認為是數學史上最複雜美麗的圖形之一。

複二次映射 $z_{n+1} = z_n^2 + c$

考慮以複數 $c \in \mathbb{C}$ 為參數的動力系統。從初始值 $z_0 = 0$ 出發,用以下的遞迴方程式產生數列 $\{z_n\}$。

$$ z_{n+1} = z_n^2 + c $$

當 $n \to \infty$ 時,不會使數列發散而保持有界的參數 $c$ 之集合,我們稱之為 曼德布洛特集合 $\mathcal{M}$。

$$ \mathcal{M} = \left\{ c \in \mathbb{C} \mathrel{\Big|} \sup_{n} |z_n| < \infty, \text{ 其中 } z_0 = 0 \right\} $$

另一方面,固定 $c$ 而改變初始值 $z_0$ 時,能讓數列保持有界的初始值之集合(的邊界),稱為 朱利亞集合(Julia set)。曼德布洛特集合可以說發揮著無限多個朱利亞集合型錄(連通性參數空間)的功能。

關於邊界豪斯多夫維度的宍倉定理

曼德布洛特集合的邊界 $\partial \mathcal{M}$,具有超乎想像複雜的碎形結構。無論怎麼放大,都會不斷出現被無數細絲相連、無限縮小的微小曼德布洛特集合複製品(迷你曼德布洛特)。

這條邊界線的「複雜度」在數學上究竟有多高呢?1998年,日本數學家宍倉光廣(Mitsuhiro Shishikura)證明了複動力系統中具有紀念碑意義的定理。

$$ \dim_H(\partial \mathcal{M}) = 2 $$

儘管只是平面(二維)上單純的邊界「線(一維性的事物)」,其豪斯多夫維度卻達到了空間維度本身的 $2$,這個事實意味著 $\partial \mathcal{M}$ 在複平面上擁有為了極限填滿空間而產生的扭曲、彎折與無數微細結構。然而,其二維勒貝格測度(面積)是否為正,仍是現代數學中未解的超級難題。


第六章:物理學與自然界中的碎形

碎形幾何學與豪斯多夫維度超越了純數學的框架,對物理學、生物學、宇宙論等所有自然科學領域產生了破壞性的深遠影響。自然界似乎比起歐幾里得幾何學,更偏好選擇碎形幾何學。

亂流(Turbulence)與流體力學

流體力學中最難解的現象「亂流」,具有碎形結構。根據能量級聯理論(由理查森和柯爾莫哥洛夫等人提出),亂流中的大漩渦會崩解為較小的漩渦,然後小漩渦再繼續崩解,以自我相似的方式重複這個過程。計算發生能量耗散區域(耗散結構)的碎形維度,已成為數學上解析納維-斯托克斯方程式的途徑之一。

布朗運動軌跡 $D=2$

$$ \dim_H(\text{Brownian path}) = 2 \quad \text{almost surely} $$

這表示,雖然是由一維時間參數所產生的曲線,但它密集地探索空間的程度,甚至具有了二維空間面積般的延展。

星系的大尺度結構(宇宙論)

抬頭仰望夜空,星辰似乎是隨機散佈的;但如果以三維方式映射宇宙大尺度下星系的分布(如史隆數位巡天等),就會浮現出由絲狀的超星系團與巨大的空洞(Void)所構成的「宇宙大尺度結構」。分析此物質分布的相關函數,顯示在某個尺度下具有碎形維度 $D \approx 1.2$ 到 $2.0$ 的自我相似性。重力造成的物質自我組織化,產生了碎形。

肺泡與血管網的最佳傳輸網絡

在生物學領域,碎形同樣普遍。人類的肺(支氣管的分支結構)、心血管系統、大腦的神經網路等,都具有碎形結構。 為什麼自然淘汰選擇了碎形呢?因為這是為了「在有限體積(空間)內,塞入無限表面積」的最佳解。支氣管透過碎形分支,在保持肺部體積固定的同時,使進行氧氣交換的表面積最大化,更將血液輸送到身體每個角落細胞的能量損失降到最低。生命的優化機制,完全符合了數學上碎形維度的法則。

結論:維度的連續性與新自然觀

歐幾里得幾何學賦予我們的「整數維度」,是人類大腦為了簡化並理解世界而極為有用的近似模型。但是,隨著測度論的發展與曼德布洛特的直覺所誕生的「碎形」與「豪斯多夫維度」,證明了維度不必然是 $0, 1, 2, 3$ 這種離散的值,也可以作為實數的連續體而存在。

豪斯多夫維度是將潛藏在自然界中的「粗糙度」、「無限細節」、「混沌中的秩序」予以量化的終極尺度。從海岸線、樹木、閃電、宇宙結構,一直到我們自身的身體結構,碎形可以說是宇宙普遍的設計語言。

作為數學抽象極致的測度論(豪斯多夫外測度),竟能如此精確地描述物理世界的實態,這個事實讓我們對存在於數學與自然科學之間神祕的對應關係留下了強烈的印象。碎形幾何學,從根本上改變了我們看待世界的方式。

comments powered by Disqus