舌头排序
阿诺德舌头代表一系列锁相区间,在锁相区中运动是周期的。
每一个舌头都与一个有理数对应,若有理数用 p/q 表示,则对应的舌头记为 p/q 舌头。
这些舌头在Ω轴上的排列次序非常奇特。在一定分辨层次上,考虑任意
三个相邻的共振区,它们对应的有理数分别为 , p1
p1
, p2 和
p2
p3 ,则有如下关
p3
系:
也就是说,中间的有理数等于两边有理数的某种“加和”,这是一种特
1 1 1
殊的“加法”:分子相加,分母相加。看两组例子,对于 , , 和
4 3 2
2 1 2
, , 分别有“加法”关系
7 3 5
1 ⊕ 1 = 1 , 2 ⊕ 2 = 1 .
4 2 3 7 5 3
1 1
这两组数中都出现了 ,那么与 邻近的到底是什么数?说清这些将涉
3 3
及振动“等级”和纯数学中的一种奇怪“序列”。
舌头排序关系恰好对应于数论研究中的法里(John Farey)树。法里是英国的一位爱好广泛、知识渊博的“小人物”,因为他并不像其他人物那样身后留下大块的文章。据说他受好收藏奇石、矿物,对音乐、钱币、车轮、数学、天文都有涉猎,经常发表一些“豆腐干”文章。使其名垂青史的只是“法里序列”。我们先看“法里树”(Farey tree):
0 1
级: ,
1 1
第2 0 1 1
1
第3 0
1
第4 0
2 1
1 1 2 1
3 2 3 1
1 1 2 1
3 2 3 1
1
第5 0
4 3 5 2 5 3 4 1
1 1 2 1 3 2 3 1 4 3
1 5 4 7 3 8
5 2 5 3 4 1
5 7 2 7 5
, , , , ,
8 3 7 4 5 1
⋯⋯
0 1
除去两边的
和 不计,第 1 级有 0 个数,第 2 级有 1 个,第 3 级有 3
1 1
个,⋯,第 n 级有 2n-1-1 个。
设第 k 级总项数为 Sn,容易看出有递推关系 Sn+1=2Sn+1,其中 S1=1。所以 Sn=2n-1-1,其中 n 为正整数。法里树还可以用连分数(continued fraction)的形式表示:
第 1 级:[],[1]
第 2 级:[],[2],[1]
第 3 级:[],[3],[2],[1,2],[1]
第 4 级:[],[4],[3],[2,2],[2],[1,1,2],[1, 2],[1,3],[1]
⋯⋯
其中[]表示 0/1,「1]表示 1/1,[a,b,c,d,⋯]是连分数的简记形式,其含义为:
p =x + 1
a+ 1
b+ 1
1
c+ d+Λ
=[xo,a,b,c,d,⋯].
因为这时 p 不大于 q,所以 xo 为 0,在上面的法里树中,为记录方便略去了前面的 0。有理数具有有限的连分数展开式,无理数具有无限的连分数展开式。任何一个二次无理数都有一个从某一位后是循环的连分数展开式。这是 1770 年拉格朗日(J.L.Lagrange)证明的一个定理。所谓二次无理数是指形如
P ± D Q
的无理数,其中 P,Q,D 都是整数,D 是正的非完全平方数。比如 就
是二次无理数。设 x 为任意无理数,则 x 可以表示成 x=[xo,a,b,c,d,⋯的形式。1891 年胡尔维茨(A.Hurwitz)推进了狄利克雷(P.G.L.Dirichlet) 有理逼近定理,证明了一个更强的逼近定理:任何一个无理数 x 都有无穷多个满足不等式
x- < 1
q≥1
的有理逼近p / q,并且数 是最好可能的数,若用更大的数去代替
不等式都不再成立。历史上胡尔维茨证明此逼近定理没有用连分数,而是用法里序列的性质!
从逼近的观点看,连分数展开式“最简单的”数,却是最坏的数,因为它们极难逼近。设
ξ= 5-1 =[0,1,1,1,Λ ]=[0,1], 2
这里的ξ就是非常有名的黄金分割数,约等于 0.618。这个数是最难逼近的。在浑沌理论中,ξ所对应的无理 KAM 环面是“最坚韧的”、坚持最久的环面,当扰动增加时,这个数对应的黄金环面(golden torus)最难破坏,因而最后被破坏。一旦它也被破坏了,系统就进入全局完全浑沌。黄金分割数ξ可以用斐波那契(Fibonacci)序列近似:
ξ ≈ Fn ,
Fn+1
Fo =0,
F1 =1,
其中 Fn 是第 n 个斐波那契数,递推关系为Fn+1=Fn+Fn-1, n≥2.
这是一个常系数齐次线性差分方程。特征方程为r2─r─1=0,
求出 r,再根据初始条件,可以容易求得 Fn 的解析表达式,结果可能出乎意料,表达式中竟然含有无理数:
1 1+
5 n 1-
5 n
Fn=
2
-
.
2
实际上根据递推公式可以迅速求出 xn=Fn/Fn+1 的极限为黄金分割数0.618⋯。递推式两侧除以 Fn 有
1 =1+x x n
n-1 ,
当 n 趋于无穷大时,xn 与 xn-1 的极限相等,所以有1=x(1+x),
解此方程得
x= -1± ,
2
舍去负值有
x=ξ= 5-1 .
2
从逼近的角度看,存在无穷多个与ξ等价的无理数。一切与ξ等价的无理数连分数展开式在最后都有与ξ相同的循环节,即都有无数个 1 在循环。浑沌研究通过法里序列、连分数、无理数逼近理论等等与数学中的皇后
——数论发生了密切关联,特别是与丢番图逼近( Diophantine approximation)联系紧密。一位活跃于浑沌动力学领域的数学物理学家斯维丹诺维奇(P.Cvitanovic)说过:“我主要参考哈代(G.H.Hardy)和莱特(E.M. Wright)的著作。”这两位数论专家写的《数论导论》是数论方面的“圣经”,如今也成了浑沌学家的必读书。
“法里序列”不同于法里树。法里序列构成规则为:满足上面所述“加和”关系,并且第 n 行真分数由所有分母小于或等于 n 的真分数组成。可以从法里树中去掉分母大于 n 的分数直接得到法里序列。
F1: 0 ,1
1 1
F2 0 1 1
: , ,
1 2 1
F3 1 1 2 1
0
: , , , ,
1 3 2 3 1
F4 1 1 1 2 3 1
0
: , , , , , ,
1 4
F5 1
3 2 3
1 1 2
4 1
1 3 2
3 4 1
0
: , , , , , , , , , ,
1
F6 0
5 4 3 5
1 1 1 1
2 5 3
2 1 3
4 5 1
2 3 4 5
: , , , , , , , , , , , ,
1 6 5 4
1
1
⋯⋯
3 5 2
5 3 4 5 6
对于 F1,序列有 0 项;对于 F2,序列 1 项;对于 F2,序列有 3 项;对于 F4,序列有 5 项;对于 F5,序列有 9 项;对于 F6,序列有 11 项;对于 F7, 序列有 17 项。对于 F100,序列有 3043 项。任意给定一个整数 n,则 Fn 所具有的法里序列项数很难直接求得。不过,数论专家已猜测到项数趋于 3n2/π
2.对于 n=100,则 3n2/π2=3039.6355⋯⋯,很接近 3043。
以圆映射为例,横坐标记Ω,纵坐标记耦合强度 K(见图 5—6),这些舌头尖朝下。舌头排列顺序严格满足法里树关系。由于分母大的锁相(共振) 难以观察到,也可以用法里序列近似表示。
随着 K 的增大,舌头变宽。这就是说当耦合强度增加时,共 振区增加了, 当 K 增大到 1 时,共振区相重叠,系统出现浑沌运动。