课堂用非线性动力学讲义第三部分

/dispbbs.asp?page=1&boardid=33&id=3262826
自然界中存在着大量复杂的几何形体（海岸线）等，它们具有以下特点： （1）形体是不规则的，其内外边界是不光滑的，但往往是“有规律的”粗糙。 （2）形体具有精细的结构，具有多重的甚至是无限重的尺度。形体没有一个统
依次放大其中的 x, y ∈ [0.55,0.7]× [0.15,0.21] ， [0.625,0.64]× [0.185,0.19] ，可以看出吸引 子的精细结构。
王跃方 《非线性动力学》大连理工大学 2010 年
24
0.21
0.190
0.20
0.189
0.19
0.188
0.18
0.187
0.17
0.16
0.186
0.15 0.55
0.60
0.65
0.70
0.185 0.625
0.630
0.635
0.640
其分形维数为 Dc = 1.2 左右。
王跃方 《非线性动力学》大连理工大学 2010 年
Newton 分形
Z n +1 = Z n − f ( Z n ) f ′( Z n ), n = 0,1,2,...
/fractal/Julia%20%E9%9B%86.htm
王跃方 《非线性动力学》大连理工大学 2010 年
11
6.4 规则分形的维数
前述的几种分形具有传统的几何学不能解释的独特特征，包括下面要讨论的维 数。在 Euclid 空间中，对维数的定义是这样的： 取一个长度为 l 的线段,把它放大两倍， 则放大后的长度为 2l ， 图形变化的倍数为 2，满足关系:
i =1 N
定义
Di = lim
I (r ) r →0 1 log r
在一般情况下， Di ≤ DC 这也是分形的特点之一。
王跃方 《非线性动力学》大连理工大学 2010 年
22
此外，还有关联维数，被定义为
d c = lim
a →0
ln ∑ Pi 2
i =1
N0
ln a
其中 a 为给定的距离， Pi 是相点落在第 i 个小盒子的概率。
1 1 例如，对于 Henon 映射的相平面，分别用 2 和 4 的 2 种方块去覆盖它，分别用
3.89，继续加密方块， 可以得到 Dc = 1.26 。 去了 95 块和 220 块，得到 Dc 分别为 6.57，
王跃方 《非线性动力学》大连理工大学 2010 年
20
再如中国的海岸线的分形维数。用不同的比例尺进行测量，得
正方体
王跃方 《非线性动力学》大连理工大学 2010 年
16
log 8 1 = D = 3 则 N = 8, β = 2 , 于是 s log 2
现在来看分形体。 （1）Cantor 集
1 log 2 = 0.6309 N = 2, β = , 于是∴ Ds = log 3 3
（2）Koch 雪花
Ultral-Fractal5.0 程序演示
王跃方 《非线性动力学》大连理工大学 2010 年
8
王跃方 《非线性动力学》大连理工大学 2010 年
9
b
王跃方 《非线性动力学》大连理工大学 2010 年
10
6.3 Julia 集
Z Z2 + c ，
Z ∈ C= , c constant ∈ C , p, q ∈ R ，Z 为坐标值。
王跃方 《非线性动力学》大连理工大学 2010 年
18
我们可以采用一个小方块（或圆片覆盖或填充）被测对象，通过统计覆盖所需 的小方块数来计算其维数，这种方法被称为盒子计数法（box 方法计算出的维数称为容量维数，用 Dc 表示。 设想用一根长度为 r 的尺子去测量线段，测到的数目与 r 成反比，即
有时候，为反映分形体在空间分布上的不均匀性，人们又提出了所谓的信息维 数的概念。其方法是先统计出分形体的某一部分落入第 i 个盒子的概率：
N N i (r ) Pi (r ) = Pi (r ) = 1 ∑ ， N (r ) i =1
来反映第 i 个盒子的填空程度。根据信息熵的定义
I (r ) = ∑ Pi log Pi (r )
2 Zn + c ( p, q )
q p c(p,q)=p+i*q
一般地说，这个格式给出以下的结果： （1） Z n
→ ∞ 当 n → ∞ （不稳定）
（2） Z n → 0 当 n → ∞ （渐近稳定） （3） Z n 有界，但 Z n
≠ 0 （稳定）
对屏幕上的所有像素点 ( p, q ) ，若它们的迭代结果不发散，则称这一点对应一个
王跃方 《非线性动力学》大连理工大学 2010 年
23
6.5 动力学中的分形现象
例1 Henon 映射中的吸引子
2 x n +1 = 1 + y n − 1.4 x n y n +1 = 0.3 y n
其相点在 R 中形成了一个吸引子，具有无标度性和自相似性，这样的吸引子具
2
有分形特征，被称为奇怪吸引子（strange attractor）
log 4 1 = 1.2618 N = 4, β = , ∴ Ds = log 3 3
王跃方 《非线性动力学》大连理工大学 2010 年
17
（3）Sierpinski 三角形
N = 3, β =
log 3 1 = 1.5849 , Ds = log 2 2
6.4 不规则分形的维数
对自然界中广泛存在的不规则的 （近似的或统计 意义上） 分形， 一般无法确知局部与整体的相似程度。
第六章
6.1 分形初步
例1 Weierstrass 曲线
∞
分形
大自然和科学现象的新几何描述
w( x) = ∑ a k cos(2π b k x)
k =0
处处连续而不可微 这对传统的基于光滑数学的空间描述方法提出了挑战。
王跃方 《非线性动力学》大连理工大学 2010 年
1
例2
Koch 曲线
n=0
n =1
（2）Koch 雪花
王跃方 《非线性动力学》大连理工大学 2010 年
14
放大倍数 L = 3, K =
4 log 4 × 3 = 4 ,∴ D f = = 1.2618 ， 1 < D f < 2 。 log 3 3
（3）Sierpinski 三角形 边长放大倍数 L = 2 ，面积放大倍数
log 3 3 = = 1.5849 ， 1 < D f < 2 。 D ) K = 3(4 × , f log 2 4
2
再取边长为 l 的立方体，把边长放大 2 倍，则放大后图形的体积为原来的 8 倍，
于是有: L3 = K ，其中当 L = 2 时， K = 8 。 由此看来，传统的 Euclid 空间的维数可以定义为：满足
王跃方 《非线性动力学》大连理工大学 2010 年
13
L
Df称为Hausdorff维数。
log N (r ) = 7.507 − 1.267 log r ，显然 Dc = 1.267 。据称，英国和挪威的海岸线的 Dc 分别为 1.3
和 1.52。 分维的大小反映了分形所占空间的程度，维数越大，空间被它占有的部分就越
王跃方 《非线性动力学》大连理工大学 2010 年
21
大，结构越致密。
3
3
3 →∞ 内边长 3 ⋅ 2
n
王跃方 《非线性动力学》大连理工大学 2010 年
3
例4
Cantor 集 ———— — — —
初始图 第一步 第二步 第三步
……
点的个数 → ∞ ，长度 → 0
例5
挪威的海岸线
王跃方 《非线性动力学》大连理工大学 2010 年
4
例 6 世界的真实面目：奇妙之旅--微观到宏观，从宏观到微观 /_ctrip/blog/item/1123193bcde742ea14cecb6f.html 视频： “十的力量”
< K , K 为某个阈值，则在该像
≥ K ，则在 ( p, q ) 描色 n 。
素点上着色，这给出了一个 M 集元素。否则若有 n = N 且 Z n 如此偏历整个荧屏，便画出了一幅 Mandelbrot 集。 Mandelbrot 的 Fortran 程序演示。
Z 0 = 0, (C ( p, q ) = ( p0 + p ⋅ ∆p ) + i (q0 − q ⋅ ∆q ))
下面介绍相似维。所谓的相似维的定义为
王跃方 《非线性动力学》大连理工大学 2010 年
15
Ds =
log N 1 log
β
其中 N 为组成某几何体的彼此相似的局部图形的数目， β 为各局部与整体的相似比。 例如：规则的正方形
log 4 1 = D 则 N = 4, β = 2 , 于是 s log 2 = 2
王跃方 《非线性动力学》大连理工大学 2010 年
6
6ห้องสมุดไป่ตู้2
Mandelbrot 集——二维世界中的分形几何体 1980 年 Mandelbrot 给出了 Mandelbrot 集。 其产生过程如下：
Z Z 2 + c ( p, q )
其中 Z , c ∈ C , p, q ∈ R 为屏幕的像素编号，令 Z 0 = 0 为给定的初值，由 此导出迭代格式 Z n= +1
L1 = K ,其中 L = 2 是线段长度的放大倍数， K = 2 是图形变化的倍数。
取一个长度（边长）为 l 的正方形，把边长放大 2 倍，则放大后正方形的面积为