单峰映象(logistic map)是种二次的多项式映射(递推关系式),是一个由简单非线性方程产生混沌现象的经典范例。
绪论单峰映象(logistic map)是种二次的多项式映射(递推关系式),是一个由简单非线性方程产生混沌现象的经典范例。这种映射因生物学家Robert May在1976年发表的一篇论文而著名。单峰映象原本被Pierre François Verhulst用作一个人口学模型,后来应用在物种受到限制因素之下的数目。数学上可写成
其中
是介于0和1之间的数,表示在第n年的物种数目。
r是正整数,是根据繁殖和饿死率而得出的数。
单峰映象是根据以下两个现象:
当人口少时,繁殖增加的人口大致跟物种的原本总数目成正比;
环境资源有一“最大容量”,当人口接近最大容量时,人口可能会下降。
可是在特定初始条件及参数时,单峰映象的人口模型会出现负的人口数,较早期使用的Ricker模型也有混沌现象,但没有这种问题。
r=4的单峰映象是位元移动映射及参数μ=2的帐篷映射的非线性变换。
r的值对结果的影响变化参数r的值,其结果如下:
0和1之间:不论启始值数值为何,xn会越来越少,最后趋近于0。
1和2之间:不论启始值数值为何,xn会快速的趋近。
2和3之间:经过几次迭代,xn也会越来越接近,但一开始会在这个值左右振动,而收敛速率是线性的。
3:xn仍然会越来越接近,但收敛速率极为缓慢,不是线性的。
3和(约3.45)之间:针对几乎所有的初值,xn最后会在2个值之间持续的震荡,即 xn最后会是a,b,a,b...的变化,其数值和r有关。
3.45和大约3.54之间,针对几乎所有的初值,xn最后会在4个值之间持续的震荡。
约大于3.54: xn最后会在8个、16个、32个值……之间持续的震荡,至于r何时会令xn的值由n个到2n个,则和费根鲍姆常数有关。
约为3.5699:这样的振动消失,整个系统开始在混沌的状态之中。针对几乎所有的初值,都不会出现固定周期的震荡,初值再微小的变化,随着时间都会使结果产生明显的差异,这就是典型混沌的特性。
大于3.5699:整个系统在混沌的状态之中。不过,当中有些特定的r值还是使系统变成非混沌,有周期性的结果,这些区间称为“稳定岛”。例如当r约大于3.82,会出现3个值的周期,再大一点出现6个值及12个值的周期。
当r从大约3.5699到大约3.8284之间,系统混沌性质发展的现象有时会称为Pomeau–Manneville场景,其特征是周期性的震荡和非周期性的行为会穿插出现。此特征可以应用在半导体元件中。也有其他区域会使系统的周期为5个值,不管任意周期都存在某特定的r,使周期为指定值。1
大于4:针对几乎所有的初值,系统最后都会超过区间[0,1]并且发散。
对于任一个r值,最多只有一个稳定的极限环,若稳定极限环存在,几乎所有的点最后都会趋近极限。若某一个r值有一个稳定的极限环,可能也有无限个不同周期的不稳定极限。
这些情况可用分枝图表示,分枝图中的横轴是参数r的数值,纵轴中显示大部分初值下,稳态可能的x值,若最后数值会在2个值中震荡,分枝图上对应的数值就会有2点。若某特定r值下,已无法明确的看到有几个对应的点,当时系统可能已经在混沌状态下。
分枝图有自相似的特性。若将分枝图中 r=3.82的部分展开,只取三个分支中的一个。其图形会好像是原分枝图缩放及扭曲后的结果。针对所有非混沌的参数r都有此一特性。以此可以看出混沌和分形深入及明显的关系。
单峰映象和混沌和其他混沌系统比较,单峰映象较为简单,是一个说明混沌特性的很好的例子。简单来说,混沌就是对初始条件的高度灵敏度。r是在3.57及4之间的大部分数值都可以使单峰映象出现此一特性。对初始条件有高度灵敏度的常见原因是映射本身是对定义域的拉伸及折叠。单峰映象的二次差分方程可视为是对于区间(0,1)拉伸及折叠的过程。
利用二维及三维的相图可以看出一些单峰映象的特性。以r=4的单峰映象为例,二维相图为一抛物线,但是若用绘制三维相图,可看出进一步的结构,例如几个一开始很接近的点在迭代后开始发散.特别是位在斜率较大位置的点。
拉伸及折叠的结果使迭代的数列以指数形式发散(参照李亚普诺夫指数),可以用有混沌特性时,单峰映象的复杂及不可预测性说明这一点。事实上,数列的指数发散说明了混沌和不可预测性之间的关系:初值微小的误差在迭代过程中会以指数成长的方式增加。因此当对于初始状态的资讯中有微小的误差时.对未来状态的预测准确度也会随迭代次数增加而快速变差。
由于映象是限制在实数数线的一段范围内,其维度小于或等于1。依数值分析的结果,在r=3.5699456...时(刚开始混沌特性时),其关联维度为0.500 ± 0.005(Peter Grassberger,1983)、豪斯多夫维数大约是0.538,而分形维数为0.5170976...。2
有些混沌系统可对于其未来状态的可能性作准确的描述。若一个可能有混沌特性的动态系统存在吸引子,则存在一概率量测描述系统长期下,在吸引子各部分所花时间的比例。以r=4的单峰映象而言,启始状态在区间(0,1)中,而吸引子也在间(0,1)中,其概率量测对应参数a=0.5,b=0.5的Β分布,其不变测度为。不可预期性和随机不同,不过在一些情形下这二很类似,而且即使对单峰映象(或其他混沌系统)初值只有很少的资讯,仍然可以针对长期的分布作某种程度的预测。
相关条目混沌理论
刚性方程
李雅普诺夫稳定性
本词条内容贡献者为:
李嘉骞 - 博士 - 同济大学