Featured image of post 分形维数与豪斯多夫测度:超越整数壁垒的小数维度与自相似性科学

分形维数与豪斯多夫测度:超越整数壁垒的小数维度与自相似性科学

曼德博集合、海岸线悖论,以及介于一维与二维之间的分数维度。几何学与测度论揭示的自然秩序。

引言:重新思考维度的概念

我们日常体验的空间被认知为三维的欧几里得空间。纸上的线是一维的,平面是二维的,而立体是三维的。自古希腊欧几里得以来,这数千年来一直是我们人类空间认知基础的坚定直觉。然而,当我们观察自然界中复杂的形状时,这种“整数维度”的范式面临着决定性的局限。云不是球体,山不是圆锥,海岸线也不是圆弧。树木的分支、血管的网络、闪电的轨迹等,自然界中可见的许多形态,都具有与平滑的欧几里得几何学对象本质上不同的“粗糙度(roughness)”。

为了用数学来描述自然界的这种复杂性,“分形几何学”作为一种新的语言应运而生。支撑其理论基础的,是诞生于实分析与测度论深渊的概念——“豪斯多夫测度”及其附带的“豪斯多夫维数”。本文将从直观的几何学到严格的测度论,全面深入地讲解分形维数是如何定义、计算,并应用于理解自然现象的。


第一章:欧几里得几何学的局限与“自然的粗糙度”

本华·曼德博之问《英国的海岸线有多长?》

有一个著名的疑问标志着分形几何学的开端。那就是本华·曼德博(Benoit Mandelbrot)于1967年在科学杂志《Science》上发表的论文标题:“英国的海岸线有多长?(How Long Is the Coast of Britain?)”。

乍看之下,这个问题似乎只是个单纯的测量问题。然而,其中潜藏着深刻的悖论。为了测量海岸线的长度,假设我们使用一定长度的尺子(例如长度 $\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)

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)

这是海里格·冯·科赫在1904年构想出的一种连续但处处不可微的曲线。

  1. 将线段分为3等分。
  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年,费利克斯·豪斯多夫扩展了卡拉西奥多里的测度论方法,定义了针对任意非负实数 $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$ 维的豪斯多夫测度。从测度论的角度来说,这是一个具有博雷尔完全性的外测度(满足卡拉西奥多里条件),并在博雷尔集族上满足可数可加性,成为一个真正的测度。

在整数维度 $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$,则利用香农熵 $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)扩展熵的“多重分形”理论中,这发展成了广义维数(雷尼维数)$D_q$ 的概念。

计箱维数(Packing Dimension)

在1980年代由 Tricot 引入的计箱维数 $\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