逻辑斯谛映射

单峰映射(英语:Logistic map)是种二次多项式映射递回关系式),是一个由简单非线性方程式产生混沌现象的经典范例。这种映射因生物学家Robert May在1976年发表的一篇论文而著名,[1]一定程度上是离散时间的种群/人口模型,类似于皮埃尔·弗朗索瓦·韦吕勒英语Pierre François Verhulst逻辑斯谛函数[2]。单峰映射实质上是逻辑斯谛函数的差分方程,其数学表达为:

其中

  • 是介于0和1之间的数,表示当前种群数量/人口数量与环境承载力的比值。
  • 是正整数,是根据繁殖和死亡率而得出的数。

单峰映射的方程旨在描述以下两个现象:

  1. 当种群数量/人口少时,繁殖增加的个体数大致跟种群原本的总数目成正比;
  2. 高密度导致的死亡,环境资源有其承载力(最大容量),当种接近最大容量时,增长率下降的速度与环境承载力减去当前种群数量的差成正比。

参数r通常取[0, 4]区间内的值,因此xn[0, 1]上保持有界。r = 4的情景是dyadic映射英语bit-shift map及参数μ=2的帐篷映射的非线性变换。当r > 4时,种群数量会出现负值(该问题在更早前的同样表现出混沌动态的Ricker模型英语Ricker model中不会出现)。也可在[−2, 0]的区间内取r值,该情形下xn有界,处于[−0.5, 1.5]之间。[3]

r的值对结果的影响

 

变化参数 的值,其结果如下:

  • 0和1之间:不论起始值数值为何, 会越来越少,最后趋近于0。
  • 1和2之间:不论起始值数值为何, 会快速的趋近 
  • 2和3之间:经过几次迭代, 也会越来越接近 ,但一开始会在这个值左右振动,而收敛速率是线性的。
  • 3: 仍然会越来越接近 ,但收敛速率极为缓慢,不是线性的。
  • 3和 (约3.45)之间:针对几乎所有的初值, 最后会在2个值之间持续的震荡,即 最后会是a,b,a,b...的变化,其数值和 有关。
  • 3.45和大约3.54之间,针对几乎所有的初值, 最后会在4个值之间持续的震荡。
  • 约大于3.54: 最后会在8个、16个、32个值……之间持续的震荡,至于 何时会令 的值由n个到2n个,则和费根鲍姆常数 有关。
  • 约为3.5699:这样的振动消失,整个系统开始在混沌的状态之中。针对几乎所有的初值,都不会出现固定周期的震荡,初值再微小的变化,随著时间都会使结果产生明显的差异,这就是典型混沌的特性。
  • 大于3.5699:整个系统在混沌的状态之中。不过,当中有些特定的 值还是使系统变成非混沌,有周期性的结果,这些区间称为“稳定岛”。例如当 约大于3.82,会出现3个值的周期,再大一点出现6个值及12个值的周期。
  •  从大约3.5699到大约3.8284之间,系统混沌性质发展的现象有时会称为Pomeau–Manneville场景英语Pomeau–Manneville scenario,其特征是周期性的震荡和非周期性的行为会穿插出现。此特征可以应用在半导体元件中[4]。也有其他区域会使系统的周期为5个值,不管任意周期都存在某特定的 ,使周期为指定值。
  • 大于4:针对几乎所有的初值,系统最后都会超过区间[0,1]并且发散。

对于任一个 值,最多只有一个稳定的极限环,若稳定极限环存在,几乎所有的点最后都会趋近极限[5]:13。若某一个 值有一个稳定的极限环,可能也有无限个不同周期的不稳定极限。

这些情况可用分枝图英语Bifurcation diagram表示,分枝图中的横轴是参数 的数值,纵轴中显示大部份初值下,稳态可能的 值,若最后数值会在2个值中震荡,分枝图上对应的数值就会有2点。若某特定 值下,已无法明确的看到有几个对应的点,当时系统可能已经在混沌状态下。

分枝图有自相似的特性。若将分枝图中 的部份展开,只取三个分支中的一个。其图形会好像是原分枝图缩放及扭曲后的结果。针对所有非混沌的参数 都有此一特性。以此可以看出混沌和分形深入及明显的关系。

 

单峰映射和混沌

 
r=4的单峰映射,二维及三维的相图,其中可看出其拉伸及折叠的特性
 
x[n+1]-x[n] vs x[n] 动画
 
单峰映射的蛛网图英语cobweb diagram,在超过3.57的大部份r值可看出其混沌的特性

和其他混沌系统比较,单峰映射较为简单,是一个说明混沌特性的很好的例子。简单来说,混沌就是对初始条件的高度灵敏度。 是在3.57及4之间的大部份数值都可以使单峰映射出现此一特性。对初始条件有高度灵敏度的常见原因是映射本身是对定义域的拉伸及折叠。单峰映射的二次差分方程可视为是对于区间(0,1)拉伸及折叠的过程。

利用二维及三维的相图可以看出一些单峰映射的特性。以 的单峰映射为例,二维相图为一抛物线,但是若用 绘制三维相图,可看出进一步的结构,例如几个一开始很接近的点在迭代后开始发散.特别是位在斜率较大位置的点。

拉伸及折叠的结果使迭代的数列以指数形式发散(参照李亚普诺夫指数),可以用有混沌特性时,单峰映射的复杂及不可预测性说明这一点。事实上,数列的指数发散说明了混沌和不可预测性之间的关系:初值微小的误差在迭代过程中会以指数成长的方式增加。因此当对于初始状态的资讯中有微小的误差时.对未来状态的预测准确度也会随迭代次数增加而快速变差。

由于映象是限制在实数数线的一段范围内,其维度小于或等于1。依数值分析的结果,在r=3.5699456...时(刚开始混沌特性时),其关联维度英语correlation dimension为0.500 ± 0.005[6](Peter Grassberger,1983)、豪斯多夫维数大约是0.538[7],而分形维数为0.5170976...[6]

 
r=3.5的单峰映射,在头三次迭代后的函数

有些混沌系统可对于其未来状态的可能性作准确的描述。若一个可能有混沌特性的动态系统存在吸引子,则存在一机率量测英语Probability measure描述系统长期下,在吸引子各部份所花时间的比例。以 的单峰映射而言,起始状态在区间(0,1)中,而吸引子也在间(0,1)中,其机率量测对应参数 Β分布[8],其不变测度英语invariant measure 。不可预期性和随机不同,不过在一些情形下这二很类似,而且即使对单峰映射(或其他混沌系统)初值只有很少的资讯,仍然可以针对长期的分布作某种程度的预测。

部份情形下的解

  的特殊情形下,单峰映射有解析解[9]。不过大部份情形下的通解只能以统计的方式预测[10] 的解为[9][11]

 

其中初始条件参数 是由 求得。针对有理数的 ,有限次数的迭代后 就会变成一个周期性的数列。不过几乎所有的 都是无理数,此时 不会重复,因此没有周期解。此解可以清楚的看出混沌的二个重要特征:拉伸及折叠。系数2n表示拉伸的指数成长,因此造成蝴蝶效应,也就是对初始值的高度灵敏性,而解中包括正弦函数的平方,使解折叠在[0, 1]的范围内。

 的解为:

 

对于 。此解没有混沌的特性。由于对于不包括不稳定固定点0在内的 ,当n趋近无限大时 会趋近于零,因此 会趋近稳定的固定点 

r = 4时找任意周期的循环

r = 4时,几乎所有的初值都会使单峰映射出现混沌特性,不过也存在无限个初值会使单峰映射最后呈周期性变化,而且对于所有正整数,都存在一初值使单峰映射的周期为该正整数。可以利用单峰映射和位元位移映射英语bit-shift map之间的关系来找出任何周期的循环。若x依照单峰映射 y依照位元位移映射

 

则二个变数的关系如下:

 .

位元位移映射其名称是因为当y以二进制表示时,映射会将二进制的数字左移一位。例如若数字是二进制的循环小数,循环节为001,则位元位移映射的序列为001001001... →010010010... →100100100... →001001001...,为周期为3的循环,循环节为010, 011, 100, 101, 110 时也会有类似情形,这些循环小数都可以转换为对应的分数,上例若以分数表示为:1/7 → 2/7 → 4/7 → 1/7。转换到r=4的单峰映射后,为611260467... → .950484434... → .188255099... → .611260467...。其他周期为3的循环也可以转换为单峰映射。依相同方式也可以找出在 位元位移映射下,任意周期的循环,再转换为单峰映射。

不过几乎所有在区间[0, 1)的数字都是无理数,而初始值为无理数的位元位移映射没有循环的特性,因此对应的单峰映射也没有循环的特性。

相关条目

参考资料

  1. ^ May, Robert M. Simple mathematical models with very complicated dynamics. Nature. 1976, 261 (5560): 459–467. Bibcode:1976Natur.261..459M. PMID 934280. S2CID 2243371. doi:10.1038/261459a0. hdl:10338.dmlcz/104555 . 
  2. ^ "埃里克·韦斯坦因. Logistic Equation. MathWorld. 
  3. ^ Tsuchiya, Takashi; Yamagishi, Daisuke. The Complete Bifurcation Diagram for the Logistic Map. Z. Naturforsch. February 11, 1997, 52a: 513–516. 
  4. ^ Carson Jeffries; Jose Perez. Observation of a Pomeau–Manneville intermittent route to chaos in a nonlinear oscillator. Physical Review A. 1982, 26 (4): 2117–2122. Bibcode:1982PhRvA..26.2117J. doi:10.1103/PhysRevA.26.2117. 
  5. ^ Collet, Pierre, and Jean-Pierre Eckmann英语Jean-Pierre Eckmann, Iterated Maps on the Interval as Dynamical Systems, Birkhauser, 1980.
  6. ^ 6.0 6.1 Peter Grassberger and I. Procaccia. Measuring the strangeness of strange attractors. Physica D. 1983, 9 (1–2): 189–208. Bibcode:1983PhyD....9..189G. doi:10.1016/0167-2789(83)90298-1. 
  7. ^ Peter Grassberger. On the Hausdorff dimension of fractal attractors. Journal of Statistical Physics. 1981, 26 (1): 173–179. Bibcode:1981JSP....26..173G. doi:10.1007/BF01106792. 
  8. ^ Jakobson, M.,"Absolutely continuous invariant measures for one-parameter families of one-dimensional maps," Communications in Mathematical Physics 81, 1981, 39-88.
  9. ^ 9.0 9.1 Schröder, Ernst. Über iterierte Funktionen. Math. Ann. 1870, 3 (2): 296–322. doi:10.1007/BF01443992. 
  10. ^ Little, M.; Heesch, D. Chaotic root-finding for a small class of polynomials (PDF). Journal of Difference Equations and Applications. 2004, 10 (11): 949–953. doi:10.1080/10236190412331285351. [永久失效链接]
  11. ^ Lorenz, Edward (1964), "The problem of deducing the climate from the governing equations," Tellus 16 (February): 1-11.
  12. ^ Veritasium. 这个方程会改变你的世界观. bilibili. 事件发生在 8:37. 2022-02-02 [2022-10-25]. (原始内容存档于2022-10-26) (英语). 

教科书

外部链接