科普
编辑
分享
反馈
  • 不变环面的破裂
  • 动手实践
  • 练习 1
  • 练习 2
  • 练习 3
  • 接下来学什么

不变环面的破裂

SciencePedia玻尔百科
定义

不变环面的破裂 是哈密顿动力系统中的一种基本现象,指系统在受到扰动时,共振机制导致其不变曲面失稳并最终瓦解。根据 KAM 定理,具有有理缠绕数的环面会在扰动下消失,而具有极不合理缠绕数(如黄金比例)的环面则能继续存在,形成了有序与混沌并存的复杂相空间。这一过程是解释行星轨道长期漂移(阿诺德扩散)、流体混沌混合以及原子级摩擦等多种物理现象的普适原理。

关键要点
  • KAM定理揭示,在微小扰动下,大多数具有“足够无理”卷绕数的不变环面会幸存,而有理环面则因共振而破裂。
  • 破裂的有理环面会转变为稳定的“岛屿”和“混沌之海”,随着扰动增强,混沌区域会通过奇里科夫重叠判据所描述的机制合并,导致大范围混沌。
  • 卷绕数为黄金分割比的“黄金环面”是最坚固的,它的破裂标志着系统从局部混沌向全局混沌转变的关键临界点。
  • 不变环面的破裂是一个普适现象,它统一了解释了从天体轨道、粒子束流稳定性到纳米尺度摩擦(超润滑)等多种跨学领域的物理过程。

引言

在物理学和数学的理想国度中,许多系统的运动可以用一种完美和谐的结构来描述:不变环面。这些高维的“甜甜圈”表面,代表着可积系统中绝对的秩序与可预测性,从行星轨道到微观粒子,它们的行为都如同精密的钟表一样周而复始。

然而,真实世界充满了不完美。一个微小的引力扰动,一个磁场中的瑕疵,都会对这个完美世界发起挑战。当现实的“微扰”之沙投入这架精密的宇宙钟表时,会发生什么?旧有的秩序是会轰然崩塌,还是会以一种更复杂的方式幸存下来?这正是动力系统理论中最深刻和迷人的问题之一。

本文将带领读者踏上一段探索之旅,揭示不变环面破裂的奥秘。在“核心概念”一章中,我们将学习KAM定理、共振和卷绕数等基本工具,理解秩序为何以及如何存续。接着,在“应用与跨学科连接”一章中,我们将看到这一理论如何在天体力学、粒子加速器物理、化学工程乃至纳米科学中展现其惊人的普适性。最后,通过一系列动手实践,你将有机会亲手模拟和验证这些迷人的动力学现象。

让我们首先从构建这个有序宇宙的基本砖石——不变环面及其核心属性开始。

核心概念

在我们深入探索之前,让我们先来玩一个思想游戏。想象一个完美无瑕的宇宙,一个由精确定律支配的钟表机构。在这个宇宙里,行星的轨道不是简单的椭圆,而是在更高维度的空间中,沿着一些光滑的、甜甜圈形状的表面运动。我们称这些表面为不变环面(invariant tori)。每个粒子,一旦出发,就永远被限制在它自己的“甜甜圈”赛道上,周而复始地进行着一种优美的、准周期的舞蹈。这个由无数嵌套的环面构成的和谐宇宙,就是物理学家所说的​可积系统(integrable system)。

节拍的重要性:卷绕数

每个环面上的运动都有其独特的“节拍”或“韵律”,我们用一个称为卷绕数(winding number),记作 ω\omegaω 的量来描述它。你可以把它想象成粒子在甜甜圈表面上,沿着短圈方向绕行的圈数与沿着长圈方向绕行圈数的平均比率。

这个数字至关重要,因为它决定了轨道的性质。如果 ω\omegaω 是一个有理数,比如 ω=3/5\omega = 3/5ω=3/5,这意味着粒子每绕长圈5次,就恰好绕短圈3次。它的轨迹最终会闭合,形成一条完美的周期性轨道。这就像一支节拍简单的华尔兹,舞步在固定的模式下不断重复。

然而,如果 ω\omegaω 是一个无理数,比如 ω=2−1\omega = \sqrt{2}-1ω=2​−1,那么粒子将永远不会回到它最初的相对位置。它的轨道会密密麻麻地缠绕在整个环面表面上,永不重复。这是一种更复杂的、非重复的韵律。

现实的闯入:微扰的力量

我们所处的真实宇宙,并非一个完美无瑕的钟表。小行星的轨道会受到木星引力的周期性“轻推”;用于约束等离子体的磁场线圈也难免存在微小的瑕疵。这些来自现实世界的不完美,在动力学中被称为微扰(perturbation)。一个微小但持续的扰动,就像在精密的钟表机构中投入一粒沙子。我们的问题是:这个由不变环面构成的和谐宇宙会发生什么?

直觉可能会告诉我们,微小的扰动只会导致微小的变化,也许环面的形状会略微扭曲,但整体结构依然存在。然而,大自然的选择远比这更加微妙和有趣。

共振的魔咒:脆弱的有理环面

想象一下推一个孩子荡秋千。如果你随意地推,时而顺势,时而逆势,秋千可能只是晃动几下。但如果你在秋千每次荡到最高点时都恰到好处地推一把,你的小力气就会累积起来,最终将秋千推得非常高。这就是共振(resonance)。

现在回到我们的环面。对于那些卷绕数为有理数(例如 ω=p/q\omega=p/qω=p/q)的环面,其上的运动是周期性的。这意味着,一个持续存在的微扰,总能在粒子运动到轨道上相同或相似的相位时,施加一个方向一致的“推力”。这个微小的推力会周期性地累积,就像推秋千一样,最终将粒子猛地推出原有的轨道,导致整个环面结构分崩离析。因此,具有简单节拍的环面,由于其周期性带来的共振效应,在微扰面前显得异常脆弱。

相反,那些拥有无理卷绕数的环面则幸运得多。由于其运动永不重复,微扰的“推力”作用在轨道上各种不同的相位上。有时是推动,有时是拉动,长远来看,这些作用力会相互抵消。就像在一个复杂的、没有固定节拍的舞蹈中,你无法通过有节奏的拍手来影响舞者的步伐。一个非常优雅的数学结果表明,对于特定的无理数运动,扰动力在长时间内的平均效应恰好为零!这正是无理环面能够抵御微扰的深层原因。

KAM定理:混沌边缘的秩序

20世纪中叶,三位杰出的数学家Andrey Kolmogorov、Vladimir Arnold和Jürgen Moser,发展出了一套惊人的理论,现在被称为​KAM定理。这个定理告诉我们,当一个可积系统受到微小扰动时,并非所有的环面都会被摧毁。

KAM定理的核心预测是:绝大多数的不变环面,特别是那些卷绕数“足够无理”的环面,将在微扰下幸存下来。它们只是被轻微地扭曲和变形,但依然将粒子的运动牢牢地束缚在自己的表面上。而那些卷绕数为有理数,或者“非常接近”有理数的环面,则会被摧毁。

这一定理揭示了一个深刻的图景:秩序与混沌可以并存。即使在引入了不完美之后,大量的规则运动结构依然存在。宇宙的钟表机构虽然出现了一些故障,但大部分零件仍在顽强地运转。

环面破碎之后:岛屿与混沌之海

那么,那些被摧毁的有理环面去了哪里?它们并没有凭空消失。取而代之的是一个更加复杂和美丽的结构:一串由更小的不变环面构成的稳定“​岛屿​”,被一片狭窄的“​混沌之海​”(或称随机层)所包围。

想象一下,原本一条平滑的环形赛道,在共振作用下破碎,变成了几个较小的环形赛道链(岛屿),而在这些小赛道之间,则是一片混乱的区域(混沌之海)。粒子如果落在岛屿上,它会继续进行规则的运动;但如果它不幸落入了混沌之海,它的轨迹将变得不可预测,仿佛失去了所有导航。

随着微扰强度的增加,这些因低阶有理数(如 1/2,2/3,3/21/2, 2/3, 3/21/2,2/3,3/2 等)破碎而产生的混沌之海会变得越来越宽。当来自不同共振区域的混沌之海扩张到相互接触时,它们就会合并成一片广阔的混沌海洋。这就是奇里科夫重叠判据(Chirikov overlap criterion)​所描述的情景。此时,粒子就可以在这片广阔的海洋中进行大范围的无序漫游,标志着大范围混沌(large-scale chaos)​的来临。

最后的秩序堡垒:黄金环面

在这个从有序到混沌的转变过程中,是否存在一个“最后的防线”?是否存在一个最坚固、最能抵抗微扰的环面?

答案是肯定的。这个最顽强的环面,其卷绕数必须是“最无理”的数——一个最难被有理数近似的数。借助连分数(continued fraction)​这一美妙的数学工具,我们可以精确地找到它。一个数越难被近似,其连分数展开中的系数就越小。那个所有系数都取最小正整数1的数,就是我们寻找的王者。这个数就是: ωgolden=11+11+11+…=5−12\omega_{golden} = \frac{1}{1+\frac{1}{1+\frac{1}{1+\dots}}} = \frac{\sqrt{5}-1}{2}ωgolden​=1+1+1+…1​1​1​=25​−1​ 是的,正是大名鼎鼎的​黄金分割比​!。

这个带有黄金分割比卷绕数的环面,被称为“最后的KAM环面”。它是在通往全局混沌的道路上,分隔开大片混沌区域的最后一道坚固屏障。当微扰强度大到足以摧毁这个“黄金环面”时,系统的大部分区域将陷入混沌的掌控之中。从行星的命运到核聚变实验的成败,这一从有序到混沌的转变,其背后隐藏的正是数字的算术性质与物理定律之间深刻而奇妙的联系。

应用与跨学科连接

在前面的章节中,我们已经深入探讨了可积系统中不变环面的精妙结构,以及它们在微小扰动下如何破裂,形成规则岛屿与混沌海洋并存的复杂画卷。这个过程,由深刻的KAM(Kolmogorov-Arnold-Moser)定理所描述,远非一个抽象的数学概念。它实际上是我们理解从行星轨道到原子摩擦等各种物理现象的一把钥匙。现在,让我们踏上一段旅程,去探索这一核心思想如何在众多科学和工程领域中开花结果,展现其令人惊叹的普适性和统一之美。

刚接触动力系统时,人们很容易陷入一种非黑即白的思维定式:一个系统要么是完全规则、可预测的(我们称之为“可积的”),就像一个理想的钟表;要么是完全混乱、不可预测的(我们称之为“遍历的”),像一锅沸水。人们可能会天真地假设,任何对理想系统的微小改动都会立即摧毁其所有秩序,使其堕入彻底的混沌。然而,大自然的真相要远比这更加丰富和迷人。KAM定理告诉我们,现实世界是一个“混合”的世界。在扰动之下,大部分规则运动的结构(不变环面)会顽强地存留下来,只是被略微扭曲;而在它们之间,狭窄的混沌区域会悄然浮现,如同大陆板块间的裂缝。正是这种有序与混沌的共存与交织,构成了我们将在下面看到的各种现象的基础。

从钟摆到银河:天体力学与粒子动力学中的回响

我们的旅程始于一个再熟悉不过的物体:单摆。一个简单的无扰动单摆,其运动是完美规则的。但在其运动的特定节点上施加一个微弱的周期性推力,情况就大为不同了。如果推力的频率与单摆自身的某个自然振动频率发生“共振”,这个推力的效应就会不断累积,最终足以打破原本稳定的轨道,在相空间中催生出一片小小的混沌区域。这个小小的桌面实验,揭示了一个宇宙级的原理。

现在,让我们将目光投向浩瀚的星空。一颗恒星在星系中的运动,远比简单的行星绕日模型复杂。它同时受到数千亿颗其他恒星引力的微弱扰动。在低能量状态下(可以类比于早期太阳系的宁静状态),许多恒星的运动是相当规则的,它们的轨道被限制在类似不变环面的结构上。然而,正如著名的Hénon-Heiles模型所揭示的,随着系统能量的增加,这些环面开始破裂,恒星的轨道变得混沌。这种混沌不是灾难,而是星系“演化”的一部分,它使得恒星能够在星系中混合,形成了我们今天观测到的复杂结构。从单摆的微小振动到星系的宏伟舞蹈,我们看到了同样的基本原理在起作用:共振导致环面破裂,从而在秩序中孕育出混沌。

如果我们将尺度再次缩小,深入到亚原子世界,同样的旋律仍在回响。在粒子加速器中,科学家们需要让质子或电子束在巨大的环形管道中稳定运行数十亿圈。每一次束流通过磁铁时,都会受到一次精确的“踢动”,这在数学上可以用一个被称为“标准映射”的模型来描述。这个“踢动”就是一种扰动。为了维持束流的稳定,我们必须确保粒子的轨道位于相空间中的“稳定岛”内——这些岛屿正是在扰动下幸存下来的KAM环面!对环面破裂机制的深刻理解,直接指导着工程师们如何设计磁场,以避免破坏性的共振,从而“驯服”这些以接近光速飞行的粒子。

工程世界的创造力:流体、化学与摩擦

不变环面的破裂不仅解释了自然现象,更启发了工程创新。想象一下,如何在没有湍流的微小通道(例如“芯片上的实验室”)中高效地混合两种液体?答案是:利用混沌!在平稳的流动中,流体粒子会沿着平滑的流线运动,这些流线就像是这个系统中的不变环面,彼此隔离,使得混合极其缓慢。然而,通过对流场施加一个微小的周期性扰动(比如一个行波),我们就可以在共振区域打破这些流线,创造出循环和拉伸的区域,即所谓的“混沌平流”。在这里,我们不是在避免混沌,而是在“按需设计”混沌,利用环面破裂的原理来达到一个非常实用的目的。

在化学工程领域,一个典型的连续搅拌釜反应器(CSTR)就如同一个不断搅拌的大锅。通过周期性地改变进入反应器的物料浓度和冷却剂温度,我们可以控制反应的进程。如果这两个控制信号的频率不成有理数比(即“非公度”的),系统就会进入一种准周期状态,其行为在相空间中表现为在一个二维环面上的运动。但是,如果我们将这些输入的扰动幅度(由参数 ε1\varepsilon_1ε1​ 和 ε2\varepsilon_2ε2​ 控制)调得太大,这个环面就会像一个被过度充气的轮胎一样“爆裂”,导致反应器的温度和产物浓度发生不可预测的混沌振荡。这种从准周期运动到混沌的转变,正是Ruelle-Takens-Newhouse理论所描述的经典路径之一。

也许最令人拍案叫绝的应用来自纳米科学。当一个原子链(例如石墨烯的一层)放置在另一个晶体表面上时,它能否自由滑动?日常经验告诉我们,物体间存在静摩擦力。但在原子尺度,情况变得微妙起来。Frenkel-Kontorova模型通过一个优雅的类比,将这个问题与标准映射直接联系起来。原子间的弹性力与衬底表面的周期性势场相互竞争。当原子链与衬底的晶格周期不成比例(非公度)且原子间相互作用足够强时,原子链的构型对应于一个光滑的KAM环面,整个链条可以无阻碍地滑动——这便是“超润滑”态,静摩擦力为零!然而,当衬底的“凹槽”势场相对于原子间的弹性变得更强时,这个环面就会破裂,发生所谓的“Aubry相变”。原子被钉扎在衬底的势阱中,滑动需要克服一个能量势垒,宏观上的静摩擦力也由此产生。一个抽象的环面破裂事件,竟是摩擦力起源的微观根源!

更深层次的连接:耗散、扩散与量子世界

通过上述例子,我们发现环面破裂似乎有两种不同的“后果”。在天体力学和粒子加速器这类能量守恒的“保守系统”中,环面破裂产生的是一个与幸存的规则岛屿共存的“混沌之海”。相空间的体积在演化中保持不变。然而,在微流体、化学反应器这类存在能量耗散的“耗散系统”中,环面破裂通常导向一个“奇异吸引子”——这是一个体积为零、具有分形维数的几何结构,附近的所有运动轨迹最终都会被“吸引”到它上面。尽管最终的混沌形态不同(“海洋”对“吸引子”),但其根本驱动力是相同的:共振环面的不稳定性。

故事还有一个更深邃的篇章。KAM定理的强大庇护主要适用于自由度较少的系统(如 N=2N=2N=2)。当系统的自由度增加到三个或更多时(例如,一个由三个耦合振子构成的系统,或者更真实的太阳系模型),情况变得更加复杂。幸存下来的KAM环面(现在是三维或更高维的“甜甜圈”)不再能像二维表面分割三维空间那样,完全隔离相空间的不同区域。它们之间存在着由破碎的共振环面构成的、错综复杂的“通道网络”,被称为“阿诺德网”。系统轨迹可以沿着这个网络极其缓慢地、近乎随机地漂移,从一个区域渗透到另一个区域。这种现象被称为“阿诺德扩散”。它是一种更隐蔽、更慢性的混沌形式,被认为是太阳系等复杂系统长期不稳定的根源之一。

最后,这场旅程将我们带到现代物理学的基石——量子力学。我们如何为一个经典系统找到其对应的量子能级?一个强大的半经典方法,即EBK量子化,其基础正是对经典运动所处的环面进行“量子化”。那么问题来了:如果一个经典系统是混沌的,它的环面已经破碎,我们该如何应用这个方法呢?答案是:我们无法应用。EBK量子化的整个框架都建立在存在良好定义的环面的前提之上。经典混沌的存在,从根本上改变了其对应量子系统的能谱结构和波函数特征。“量子混沌”这一整个领域,正是从这个深刻的问题中诞生的。

至此,我们的旅程暂告一段。从一个简单的单摆出发,我们穿越了星辰大海,深入了工程设备的核心,触摸了摩擦的本质,甚至窥见了量子世界的奥秘。所有这些看似毫不相关的现象,都被“不变环面的破裂”这一统一而优美的思想紧密地联系在一起。这正是物理学最激动人心之处:在纷繁复杂的表象之下,寻找那贯穿一切的、简洁而深刻的规律。

动手实践

练习 1

在复杂的动力系统中,我们常常需要借助数值工具来判断运动是规则的(在不变环面上)还是混沌的。最大李雅普诺夫指数 (λmax⁡\lambda_{\max}λmax​) 是诊断混沌的“黄金标准”,它量化了系统对初始条件的敏感性。这个练习将检验你解释计算结果以对轨道行为进行分类的能力,这是区分环面上规则运动和“混沌海”中无序运动的第一步。

问题​: 一位等离子体物理研究人员正在研究磁约束场中带电粒子运动的一个简化模型。在某个二维相空间投影中,粒子轨道的动力学由一组方程决定,这些方程依赖于单个可调控制参数 ϵ\epsilonϵ。对于较小的 ϵ\epsilonϵ 值,理论工作预测粒子轨道是准周期的,这意味着它们位于相空间中的一个不变环面上。当 ϵ\epsilonϵ 增加超过一个临界值时,这些环面预计会破裂,粒子的运动变得混沌,探索相空间中一个被称为“混沌海”的更大区域。

为了研究这一转变,该研究人员对几个不同的控制参数 ϵ\epsilonϵ 值,数值计算了最大李雅普诺夫指数 λmax\lambda_{max}λmax​。最大李雅普诺夫指数用于量化两条初始距离无穷小的轨道之间指数发散的平均速率。正的 λmax\lambda_{max}λmax​ 是混沌的一个标志,而对于环面上的规则准周期运动,λmax\lambda_{max}λmax​ 为零。由于数值模拟的有限时间和浮点运算,理论上为零的指数将被测量为一个在零附近波动的很小的值。计算出的指数单位是秒的倒数 (s−1s^{-1}s−1)。

下面显示了四个计算实验的结果:

实验ID控制参数 (ϵ\epsilonϵ)计算出的最大李雅普诺夫指数 (λmax\lambda_{max}λmax​),单位 s−1s^{-1}s−1
10.150.0030.0030.003
20.480.570.570.57
30.71−0.008-0.008−0.008
40.921.141.141.14

根据这些数据,以下哪个陈述正确地分类了实验中观察到的运动?

A. 实验1和3显示了不变环面上的运动;实验2和4显示了混沌运动。

B. 实验2和4显示了不变环面上的运动;实验1和3显示了混沌运动。

C. 实验1、2和3显示了不变环面上的运动;实验4显示了混沌运动。

D. 实验1显示了不变环面上的运动;实验2、3和4显示了混沌运动。

E. 所有实验都显示了混沌运动。

显示求解过程
练习 2

当一个可积系统受到微扰时,并非所有的环面都会以相同的方式被破坏。KAM 定理告诉我们,具有足够“无理”频率的环面会存活下来,而具有“有理”频率的环面则会破裂成一系列稳定的岛链。这项实践将深入探讨这一过程的核心,要求你分析计算出主共振岛的大小,这是微扰后出现的混合相空间中的一个关键特征。

问题​: 考虑一个单自由度哈密顿系统,其未扰动动力学由哈密顿量 H0(I)=12αI2H_0(I) = \frac{1}{2} \alpha I^2H0​(I)=21​αI2 描述,其中 III 和 θ\thetaθ 是正则作用量-角变量,α\alphaα 是一个正常数。该系统受到一个弱的时间周期性微扰,得到的总哈密顿量为: H(I,θ,t)=12αI2+ϵVcos⁡(kθ−mΩt)H(I, \theta, t) = \frac{1}{2} \alpha I^2 + \epsilon V \cos(k\theta - m\Omega t)H(I,θ,t)=21​αI2+ϵVcos(kθ−mΩt) 此处,ϵ\epsilonϵ 是一个微小的无量纲正常数,VVV 是一个具有能量单位的正常数,kkk 和 mmm 是正整数,Ω\OmegaΩ 是微扰的角频率。

这种微扰可以在相空间中导致共振岛的形成。当系统的未扰动频率 ω(I)=dH0dI\omega(I) = \frac{dH_0}{dI}ω(I)=dIdH0​​ 与微扰频率 Ω\OmegaΩ 可公度时,会在一个特定的作用量值 IrI_rIr​ 处发生主共振。对于由整数 kkk 和 mmm 定义的共振,在原始共振轨道附近会形成一系列被称为“岛”的稳定区域。这一系列岛的边界是一条分界线,它将岛内部的捕获运动与外部的非捕获运动分离开来。

确定该共振区域的全宽度 ΔI\Delta IΔI,其定义为作用量变量 III 沿分界线轨道的最大值与最小值之差。将您的答案表示为关于 α,ϵ,\alpha, \epsilon,α,ϵ, 和 VVV 的符号表达式。

显示求解过程
练习 3

现在,让我们将所有概念融会贯通。最后的这项实践是一个计算性的小项目,你将使用著名的模型系统——标准映射(Standard Map),来可视化环面破裂的整个过程。你将通过编程为不同的微扰强度 KKK 生成庞加莱截面,并使用你在第一个练习中学到的李雅普诺夫指数来量化系统从有序到混沌转变过程中不断增长的“混沌海”。

问题​: 你的任务是通过庞加莱截面为一保面积映射构建并分析相空间相图。该映射模拟了一个周期性受踢转子,并且是 Kolmogorov–Arnold–Moser (KAM) 定理的一个典范例子。考虑在二维环面上定义的标准映射,其通过映射 ΦK:(θn,pn)↦(θn+1,pn+1)\Phi_K : (\theta_n, p_n) \mapsto (\theta_{n+1}, p_{n+1})ΦK​:(θn​,pn​)↦(θn+1​,pn+1​) 的迭代给出:

pn+1=(pn+Ksin⁡θn) mod 2π,θn+1=(θn+pn+1) mod 2π,\begin{aligned} p_{n+1} &= \left(p_n + K \sin \theta_n\right) \bmod 2\pi,\\ \theta_{n+1} &= \left(\theta_n + p_{n+1}\right) \bmod 2\pi, \end{aligned}pn+1​θn+1​​=(pn​+Ksinθn​)mod2π,=(θn​+pn+1​)mod2π,​

其中 K≥0K \ge 0K≥0 是一个微扰参数,θ∈[0,2π)\theta \in [0,2\pi)θ∈[0,2π) 是一个角度,而 p∈[0,2π)p \in [0,2\pi)p∈[0,2π) 是一个类作用量变量。集合 {(θn,pn)}n≥0\{(\theta_n, p_n)\}_{n \ge 0}{(θn​,pn​)}n≥0​ 构成了庞加莱截面。所有角度必须以弧度为单位。

将映射 ΦK\Phi_KΦK​ 在相点 (θ0,p0)(\theta_0,p_0)(θ0​,p0​) 处的最大李雅普诺夫指数 λmax⁡\lambda_{\max}λmax​ 定义为以下极限:

λmax⁡(θ0,p0;K)=lim⁡N→∞1N ln⁡∥DΦKN(θ0,p0) v0∥,\lambda_{\max}(\theta_0,p_0;K) = \lim_{N \to \infty} \frac{1}{N}\,\ln \left\| D\Phi_K^N(\theta_0,p_0)\, v_0 \right\|,λmax​(θ0​,p0​;K)=N→∞lim​N1​ln​DΦKN​(θ0​,p0​)v0​​,

其中 DΦKND\Phi_K^NDΦKN​ 是该映射沿起始于 (θ0,p0)(\theta_0,p_0)(θ0​,p0​) 的轨道求值的 N 次复合的导数(雅可比矩阵),∥⋅∥\|\cdot\|∥⋅∥ 是欧几里得范数,v0≠0v_0 \ne 0v0​=0 是任意单位切向量。对于标准映射的单步雅可比矩阵 Jn=DΦK(θn,pn)J_n = D\Phi_K(\theta_n,p_n)Jn​=DΦK​(θn​,pn​),有:

Jn=(1+Kcos⁡θn1Kcos⁡θn1).J_n = \begin{pmatrix} 1 + K \cos \theta_n & 1\\ K \cos \theta_n & 1 \end{pmatrix}.Jn​=(1+Kcosθn​Kcosθn​​11​).

如果一条轨道的有限时间李雅普诺夫指数估计值 λ^max⁡\hat{\lambda}_{\max}λ^max​ (沿其轨道计算 N 次迭代得到) 满足 λ^max⁡>τ\hat{\lambda}_{\max} > \tauλ^max​>τ(其中 τ>0\tau > 0τ>0 是给定的阈值),则该轨道被认为是混沌的;否则,它被认为是规则的。你的任务是针对每个指定情况,使用上述定义数值近似 λ^max⁡\hat{\lambda}_{\max}λ^max​,并确定混沌初始条件的比例。

测试套件和规格:

  • 使用微扰参数集 KKK 等于四个值 K∈{ 0.0, 0.2, 0.98, 1.2 }K \in \{\,0.0,\,0.2,\,0.98,\,1.2\,\}K∈{0.0,0.2,0.98,1.2}。
  • 使用由 nθ=8n_\theta = 8nθ​=8 和 np=8n_p = 8np​=8 定义的均匀初始条件网格。具体而言,取
θ0(i)=2π inθ,i∈{0,1,…,7},p0(j)=2π jnp,j∈{0,1,…,7}.\theta_0(i) = \frac{2\pi\, i}{n_\theta},\quad i \in \{0,1,\ldots,7\}, \qquad p_0(j) = \frac{2\pi\, j}{n_p},\quad j \in \{0,1,\ldots,7\}.θ0​(i)=nθ​2πi​,i∈{0,1,…,7},p0​(j)=np​2πj​,j∈{0,1,…,7}.

这将在环面 [0,2π)×[0,2π)[0,2\pi)\times[0,2\pi)[0,2π)×[0,2π) 上产生 nθ×np=64n_\theta \times n_p = 64nθ​×np​=64 个初始条件 (θ0(i),p0(j))(\theta_0(i), p_0(j))(θ0​(i),p0​(j))。

  • 对于每个初始条件和每个 KKK 值,将映射迭代 N=3000N = 3000N=3000 步以近似 λ^max⁡\hat{\lambda}_{\max}λ^max​,并使用混沌分类阈值 τ=10−2\tau = 10^{-2}τ=10−2,即如果 λ^max⁡>10−2\hat{\lambda}_{\max} > 10^{-2}λ^max​>10−2,则将轨道分类为混沌,否则分类为规则。
  • 角度必须以弧度处理,所有计算都是无量纲的。

交付成果:

  • 对于测试套件中的每个 KKK 值,计算 64 个初始条件中混沌条件的比例。将这四个比例报告为四舍五入到三位小数的小数。

最终输出格式:

  • 你的程序应生成单行输出,其中包含结果,形式为由方括号括起来的、无空格的四个四舍五入分数的逗号分隔列表。例如,一个可接受的输出格式是
[0.000,0.031,0.734,0.953].[0.000,0.031,0.734,0.953].[0.000,0.031,0.734,0.953].
显示求解过程
接下来学什么
动力系统
尚未开始,立即阅读
KAM 定理
保守系统