走向混沌的道路

走向混沌的道路第一页，共五十页，编辑于2023年，星期一第五节

保守系统中的不规则运动1．可积与不可积系统2．扰动与KAM定律3.有理环面破裂与同（异）宿结构4.阿诺德扩散5.标准映射第二页，共五十页，编辑于2023年，星期一保守系统1．可积与不可积系统

分析力学里人们常用广义坐标q和广义动量p来表示系统的变量，它们是系统哈密顿量正则共轭变量。第一项为动能，第二项为势能。如果系统仅受势能力作用，则系统机械能守恒，称为保守系统。实际系统还受耗散力作用，系统能量就不守恒，因此称为耗散系统（如：阻尼单摆、范德玻耳方程、洛伦兹方程、平方映射等。)

取相空间某区域的全部状态为初始状态，区域形状将会因各代表点的运动速度可能不同而发生变化。耗散系统因能量散失而有吸引子存在，相空间内所有轨线都要收缩到吸引子上，初始一定大小的相空间在运动中逐渐减小，在t→∞时趋向于零。保守系统的里不存在吸引域，也就不会有吸引子，相空间是守恒的。第三页，共五十页，编辑于2023年，星期一保守系统耗散系统的相空间在运动中逐渐减小，在t→∞时趋向于零保守系统的相空间是守恒的1．可积与不可积系统第四页，共五十页，编辑于2023年，星期一可积系统定义一类力学系统哈密顿系统的正则方程

例：无阻尼的自由单摆它的广义坐标为摆角，q=θ，广义动量为单摆的其哈密顿量为：1．可积与不可积系统第五页，共五十页，编辑于2023年，星期一可积系统由正则方程得无阻尼单摆方程：两边同乘得全微分方程

积分

积分常数C称为运动积分

说明自由单摆是一个完全可积系统，运动状态由能量决定1．可积与不可积系统第六页，共五十页，编辑于2023年，星期一可积系统

不是所有的哈密顿系统都可以化成全微分形式，因此它们是不可积的（或者说没有运动积分）。一般将哈密顿系统分两类：完全可积的与不可积的。实际上完全可积系统是极少数，绝大多数系统是不可积。对一个N自由度的保守系统，其哈密顿为N对广义动量

p1,p2,…,pN与广义坐标q1,q2,…,qN

的函数，运动方程为:如果完全可积，要求有N个独立全微分方程，即要有N个独立的运动常数

C1．可积与不可积系统第七页，共五十页，编辑于2023年，星期一

由正则共轭变量(p,q)可变换出另一对正则共轭变量，称作用-角度变量。对单自由度系统

母函数

用了变量以后，运动方程为积分得

上式说明I与系统能量E一样，也是运动积分，I=I(H)可以写为I(E)。第二式定义了系统运动频率：，它不是常数，与I有关，是非线性的。的一个重要特性是周期性，在一个周期内S的变化为：

作用-角度变量1．可积与不可积系统第八页，共五十页，编辑于2023年，星期一简谐振子以单位质量的简谐振子为例，它的哈密顿函数为：由得由哈密顿函数可得采用作用变量与角变量定义，得1．可积与不可积系统第九页，共五十页，编辑于2023年，星期一环面运动

采用作用变量与角变量之后，一个保守系统可以用相空间内的环面上的运动来表示。对N自由度系统，一切周期的或准周期的有界运动，是在2N维环面上的具有N个频率的运动。对2自由度系统，有两组作用-角度变量：与，相空间是四维的。给出一组的同心圆，是环绕这组圆的转角，给出另一组同心圆。如能量E恒定，相空间中的运动由四维降低到三维，运动限制在三维能面上。在为常数的二维环面上，有：1．可积与不可积系统第十页，共五十页，编辑于2023年，星期一环面运动

当系统两个特征频率比

w=ω1/ω2为有理数，运动是周期的，轨线经有限次的绕环后闭合，环面由无数条周期轨道组成，称为有理环面。如w为无理数，称为无理环面，环面由无数条准周期轨线组成，每一条轨线都随时间一圈一圈地覆盖整个环面。如在=常数处设置一平面(庞加莱截面)，则截面上那些环面曲线与截面交点构成了一个圆。二维环面运动1．可积与不可积系统第十一页，共五十页，编辑于2023年，星期一扰动2．扰动与KAM定律

设系统未受扰动时运动是可积的，哈密量为，受扰动时的哈密量为：ε为无量纲参数，它决定扰动的强度。如周期为T的扰动，则扰动可为将展开成级数式中n为扰动频率，n,m

某整数。

代入运动方程积分得第十二页，共五十页，编辑于2023年，星期一

将零级近似代入运动方程解得一级近似解系统的运动频率与I有关，当在某个值Ir上出现扰动频率v与系统ω(Ir)频率间的共振时，则：或因为ω(Ir)与Ir有关，因此此时的共振称为非线性共振。当发生非线性共振时，一级近似解中分式的分母等于零，得到发散得结果，这就是著名的小分母发散问题。扰动2．扰动与KAM定律第十三页，共五十页，编辑于2023年，星期一扰动将对系统产生两种不同的影响：（1）非线性共振：一个很小的扰动可将导致有理环面发生重大改变。这相应于相空间的共振有理环面。（2）非共振情况：相应于无理环面。这是动力学的一个久未解决的问题。1954年，前苏联数学家哥尔摩格洛夫(Kolmogorov)提出了一个环面不变定理，后来为阿诺德(Arnold)所证明，美国数学家莫瑟(Moser)在某些条件下也证明了该定理。现在常称为KAM定理。该定理考虑一个近可积系统，即对一完全可积系统施加了一个很小的完全不可积扰动。KAM定理：如果扰动很小，大多数非共振的不变环面并不消失，只是发生一些微小的变形。满足KAM定理的绝大多数轨道，其运动仍然限制在N维环面上，环面上的运动仍然是准周期的。这些未被破坏的环称为KAM环。

KAM定理2．扰动与KAM定律第十四页，共五十页，编辑于2023年，星期一3.有理环面破裂与同(异)宿结构

扰动的一级近似解没有回答非线性共振对有理环面的影响。数学上可用庞加莱截面解决。

未受扰动时，在给定能面中取q2=常数的庞加莱截面，轨线与截面的交点在I1=常数的圆上。轨线相继两次穿越截面的时间间隔：

q1

每次的改变量为庞加莱截面上点的运动为一二维映射，称为扭转映射：

受扰动时（ε≠0），扭转映射T0变成Tε，略去下标后有式中f与g由扰动项确定。扭转映射第十五页，共五十页，编辑于2023年，星期一3.有理环面破裂与同(异)宿结构

研究绕卷数的有理环面在扰动下的变化。记有理环面为Go

，它为一个由映射Tom的不动点组成的圆。两条不变曲线G+

与G-

，它们分别位于Go

的两边。在Tem的作用下，在每个q=常数的圆半径上，存在着转动角度刚好2p的点，它们只有径向移动而没有转动，这些点连结起来构成了在Tem作用下的曲线Ge

，Tom的映像闭合曲线Tom(Ge)。曲线Ge与Tom(Ge)两者保围的面积相等且相交，共有2m个交点。扭转映射对有理环面作用第十六页，共五十页，编辑于2023年，星期一3.有理环面破裂与同(异)宿结构

Tom的作用：Go圆上的点刚好转动2p,曲线G+上的点转动小于2p,曲线G-的点就转动大于2p。因此Go上的点不动，G+的点顺时转动，G-的点逆时转动。

Tem的作用：如扰动eV很小，Tem的作用不会改变与圆上点转动情况，顺时针或反时针转动均不会变化，只是在每个q=常数的圆半径上，存在着转动角度刚好2p的点，它们只有径向移动而没有转动，将这些点连结起来，就构成了在Tem作用下的曲线Ge

。

此外,还有Tom的映像闭合曲线Tom(Ge)。由于保守系统，曲线Ge与Tom(Ge)两者不仅保围的面积相等，且相交，共有2m个交点。根据相交点附近点移动的走向，可判定其中一半是椭圆点，另一半是双曲点，它们相间地分布。

椭圆不动点附近有较小有理环面，是一些区域较小的规则运动。扰动也将使它们受到破坏，情况与上面相类似。扰动产生更高一级的椭圆不动点及围绕它们更较小一级的规则运动区。如此的破坏过程还会继续发展下去，以至产生规则与不规则运动交织在一起的无穷堪套的自相似结构。扭转映射对有理环面作用第十七页，共五十页，编辑于2023年，星期一3.有理环面破裂与同(异)宿结构

无阻尼单摆，或负线性恢复力杜芳方程相图上，均有四条流线通过界轨线上的双曲不动点，其中两条流向双曲点，另两条背离双曲点。

数学上称这些流线为不变曲线或流形(manifold)。与相图上的真正相轨线不同，现在这些双曲点出现在环面截面上，它们由截面上的点构成。

两条背离双曲点的流线称为稳定流形ws，两条流向双曲点的流线称为不稳定流形wu。若沿着一条稳定流形ws从双曲不动点O出发，将会连接到的一条不稳定流形wu进入双曲点O。若前后两个是不同的点,称为异宿点，当前后两个是同点,称为同宿点。同(异)宿点第十八页，共五十页，编辑于2023年，星期一3.有理环面破裂与同(异)宿结构(1)

系统受到小幅周期性扰动，对同宿点，代表点沿着稳定流形离开双曲点后，不会连接到不稳定流形而进入双曲点，一般是与不稳定流形相交(异宿点情况相似)；

(2)

一旦发生相交，交点上的流形与同宿点或异宿点特点相同：即两条流入流线与两条流出流线，但它们不是不动点。同(异)宿结构过双曲点是不变曲线，线上的点经过多次映射不会跑出该线。扰动的影响是：同宿结构第十九页，共五十页，编辑于2023年，星期一3.有理环面破裂与同(异)宿结构同(异)宿结构异宿结构(3)由于映射是连续的，在双曲点O外会产生一系列新同宿点，且要反复作用无限次数才能沿着接近到双曲不动点O。

(4)在到达双曲不动点O前，流线与交叉产生的同(异)宿点越来越密，总共有无限多个点。由于是保守系统，相继两次交叉所包围的面积是定值，于是越趋近双曲不动点O，在新宿点会越来越密同时，振荡幅度越来越大。扰动的影响是：第二十页，共五十页，编辑于2023年，星期一4.阿诺德扩散

若不可积扰动足够小，KAM无理环面仍然可以保持，有理环面发生破裂，产生出一系列新椭圆不动点与双曲不动点。新生椭圆不动点在扰动连续作用下，又产生出更小一级的椭圆不动点与双曲不动点；在双曲不动点附近，通过不动点的稳定与不稳定流形构成复杂的异宿结构。它们是二维环面上的非规则运动区。KAM环面将环面上规则与非规则的运动区域分隔开来。在整个环面上共存了规则的与非规则的运动。扰动对二维环面运动影响

第二十一页，共五十页，编辑于2023年，星期一4.阿诺德扩散

高维KAM环面能否成为等能面的边界，对不满足KAM定理的导致不规则运动的少数轨道(即不稳定轨道)起限制作用。N个自由度系统具有2N维的相空间和2N－1维等能面。

N维环面成为等能面边界的条件：N≥2N－2。即只有N≤2系统的环面才有可能把等能面包围起来或分割成几个部分。N≥3系统不会满足这样条件。高维相空间里KAM环面不会被等能面隔离，不稳定轨道会在各个KAM环面之间来回穿行，并逐步扩散到整个等能面上去，这种现象被称为阿诺德扩散。

高维空间的阿诺德扩散第二十二页，共五十页，编辑于2023年，星期一5.

标准映射

与耗散系统中两系统间的同步与锁模相对应，单自由度保守系统受周期性外力作用时产生不规则运动。系统哈密顿为：运动方程为：选取时间系列：t0,t1,t2,…，运动方程退化为离散映射：该映射给出I

和q

在前后两个时间点上关系，又可写成：离散映射T为扰动的特征时间第二十三页，共五十页，编辑于2023年，星期一5.

标准映射稳定性条件设扰动势只与广义坐标q

有关：

二维映射

二维映射稳定性条件为：本征值方程本征值方程的解第二十四页，共五十页，编辑于2023年，星期一5.

标准映射采用无量纲作用量略去常数相位因子由不动点的条件即：标准映射可得由此可得两个奇点

与耗散系统圆映射情况有点相似，在保守系统中，标准映射给出的复杂运动行为与参数K取值有关。对不同参数K值进行数值计算，并画出I-q平面上的相。

标准映射第二十五页，共五十页，编辑于2023年，星期一第二十六页，共五十页，编辑于2023年，星期一5.

标准映射随机海K=0.5时，在周期1与周期2区之间，除了一些横贯左右的点线外，还存在一些光滑曲线。点线是破裂了的有理线，光滑曲线是没有破裂的KAM不变线。说明随机性受KAM不变线约束而存在于局部区域内；

K=1.0时，周期1与周期2区之间的光滑曲线消失，成了弥散点，KAM线破裂了。随KAM不变线破裂，局部随机区逐步向全局扩散，汇成广泛弥散开的大海—随机海(Stochasticsea)。这是一种全局性的混沌形态。第二十七页，共五十页，编辑于2023年，星期一5.

标准映射标准映射相图第二十八页，共五十页，编辑于2023年，星期一第六节电子混沌电路1.外激励的非线性LC谐振电路2.微分方程的模拟电子电路3.实际动力体系的电子模拟电路

第二十九页，共五十页，编辑于2023年，星期一

如果在通常的LC谐振电路中，使用非线性电阻、电容或电感等一些非线性电子元件，如果这些元件的数值（电阻值、电容量或电感量）是外加电压或频率的函数，就构成一种非线性LC

谐振电路，在合适的驱动电压作用下，将会呈现非线性动力学系统的混沌特性。这里以单结晶体管作为非线性电阻，变容二极管作为非线性电容为例，介绍一下外激励非线性LC

电路进入混沌的基本情况。

1.外激励非线性LC谐振电路非线性LC谐振电路

第三十页，共五十页，编辑于2023年，星期一

这里将一只单结晶体管联结为一个非线性电阻。单结晶体管也称双基极二极管，它有两个基极b1

和b2和一个发射极e

。e对基极b1和b2就是一个PN

结，具有单向导电性。两个基极之间的电阻是半导体的内电阻。当b2-

b1

间加上电压ubb，b2-

b1

间的发射极电位决定于内电阻上的分压比h：

当发射极外加电压高于e点电位的二极管压降ud(0.6V)时，ue>hubb+ud,e-

b1导通。在混沌电路中，通常将基极b2和发射极联在一起，成了二端元件。当ue<hubb+ud时,单结晶体管的伏安特性呈正电阻，电路中电流为b2-

b1

电流，当ue>hubb+ud时，e-

b1导通，随电流增大，两端电压下降，成负电阻。1.外激励非线性LC谐振电路单结晶体管伏安特性第三十一页，共五十页，编辑于2023年，星期一单结晶体管混沌电路

1.外激励非线性LC谐振电路

单结晶体管可选BT33D，电路如下。这是串联的强迫振荡电路，非线性电阻(R2)即单结晶体管，电感L和电容C，电阻R1、R3都是线性元件。电阻R3（也是电流取样电阻）和电源构成了单结晶体管的偏置电路。L＝30mH，C＝0.03μF，R1＝125Ω，R3＝33.8kΩ，Ec＝28V

第三十二页，共五十页，编辑于2023年，星期一电路中的混沌状态

设驱动为：以f与u0为该电路的两个可变参数，研究取不同的数值时系统的运动状态：

研究发现存在以下一些工作状态：

①驱动电压的幅度u0＝0时，电路处于单稳状态；

②固定驱动电压幅度，改变驱动信号的频率f，当逐步接近电路的固有频率时，电路出现突变的锁频状态；

③固定驱动电压的频率f，将驱动电压的幅度u0由小到大增加，在u0的一定变化范围内出现倍周期分岔与混沌运动。但在不同频率f下，出现倍周期分岔与混沌的u0值范围不同，倍周期分岔与混沌的过程也不一样；

④当f远离f0时，随输入信号幅度值逐步升高，电路出现倍周期分岔－混沌－反倍周期的过程；

⑤当u0值进一步升高时，倍周期分岔与混沌现象消失，电路表现出一般非线性电路所共有的畸变波形的特征。1.外激励非线性LC谐振电路第三十三页，共五十页，编辑于2023年，星期一电路中的混沌状态

设驱电压动为：

以f与u0为该电路的两个可变参数，研究取不同的数值时系统的运动状态：

固定驱动电压的频率f，将驱动电压的幅度u0由小到大增加，在u0的一定变化范围内出现倍周期分岔与混沌运动。

1.外激励非线性LC谐振电路第三十四页，共五十页，编辑于2023年，星期一

这是一个广泛研究过的电路。通常它由电阻R，电感L与一只硅二极管D阻成。硅二极管在反向偏置时其极间电容即为回路电容，且其容量随电压变化而变化，但与所加电压是非线性关系。正向偏置时为一直流电压源。图中右侧的场效应管电路是提高输入阻抗以便用示波进行测量。电路中，R＝100W，L＝0.25H。1.外激励非线性LC

谐振电路二极管－电感混沌电路

第三十五页，共五十页，编辑于2023年，星期一二极管－电感混沌电路

实验时，电路上加上正弦电压，并使外加电压频率与由电感和二极管反向电容组成的回路的固有频率相接近，即接近共振状态。逐步改变输入电压的大小，同时用示波器检测二极管上的电压波形。当输入电压很小时，输出电压为高度相等的整流波形。随着外加电压的增加，二极管上出现一串高度不等的整流波形电压。这时二极管电压经历了倍周期分岔，并最终进入混沌。根据实验测波形可以计算费根鲍姆常数等各项常数，结果如表。

11―15db13.2db功率谱中平均峰高比6.3±0.36.619噪声指数2.4±0.12.50…费根鲍姆第二常数4.26±0.14.66920…费根鲍姆第一常数实验值理论值各项数值1.外激励非线性LC

谐振电路第三十六页，共五十页，编辑于2023年，星期一2.微分方程模拟电子电路

能产生混沌行为的典型非线性常微分方程是由三个独立变量的一阶微分方程，如洛伦兹方程或罗斯勒方程。在研究动力系统的相图时，常将一个二阶微分方程化为两个一阶方程。同样地，可以通过数学上变换，可将三个一阶微分方程化为一个三阶微分方程来。因此可将产生混沌行为非线性微分方程写成一个三阶方程：式中为各项系数，为一非线性函数，系数为包括0在内的各种正负常数，因此从一般方程可变换出许多形式的非线性方程。但不是所有的方程都具有混沌行为。三阶方程的系统有三个李雅普诺夫指数，具有混沌行为要求至少有一个是正值。一个系统的李雅普诺夫指数之和代表了相体积沿轨道的平均变化速率，如果三个指数之和小于零，，系统相体积在运动中会逐渐减小，是耗散系统；如等于零，，相体积不变，此为保守系统。三阶非线性微分方程第三十七页，共五十页，编辑于2023年，星期一非线性函数电路2.微分方程模拟电子电路非线性函数f(x)

在电路上可用二极管及线性放大器的适当的联接来实现。第三十八页，共五十页，编辑于2023年，星期一

f(x)=|x|

的电路

这种情况的微分方程可以如下构成：式中只有一个控制参数A。由图可见，该电路以为电压驱动，用了三个串接的反向积分器以产生，和x等三个量，并一定比例将三个信号及一个由电池产生的直流电压相加起来，再反馈到第一个积分器的输入端。

2.微分方程模拟电子电路第三十九页，共五十页，编辑于2023年，星期一

f(x)=|x|的电路2.微分方程模拟电子电路

为了得到合适的参数A取值范围，可计算方程解x(t)

对A的分岔图。取方程右边的正负号为负号，初始条件为，得的分岔图。

可见当参数A由0.8向0.64085逼近时，系统以倍周期分岔进入混沌。A值小于0.64085后，出现类似于平方映射那样的带有大大小小窗口的典型的混沌带，而在附近出现周期3窗口，是不稳定的。方程在这里的长时间演化是对无穷大发散的。第四十页，共五十页，编辑于2023年，星期一2.微分方程模拟电子电路

f(x)=|x|的电路

f(x)=|x|混沌电路的三个李雅普诺夫指数为:可见三个指数之和小于零，这是耗散系统。所以该电路存在奇怪吸引子。第四十一页，共五十页，编辑于2023年，星期一单二极管非线函数电路

这是用单二极管完成的非线性函数为R(x)微分方程：该方程的三个李雅普诺夫指数分别为(0.042,0,0.342)，也是耗散系统。电路上与上面f(x)=|x|电路基本相同，主要差别在以阻容积分电路代替了第二个有源积分器。电路中各个电阻值均为1kΩ。类似地可以用非线性函数D(x)组成的微分方程，与R(x)电路的差别在只在二极管正反向不同。

2.微分方程模拟电子电路第四十二页，共五十页，编辑于2023年，星期一跃变非线性函数电路

跃变非线性函数利用了运算放大器的内在非线性特性，即理想放大器开环特性。当输入电压过零时，理想放大器输出从负饱和值跃变到正饱和值。一个简单的以跃变非线性函数的微分方程为：式中不同正负号有不同的李雅普诺夫指数，取正号时为，[0.152,0,0.652]，负号时为[0.601,0,1.101]。2.微分方程模拟电子电路取方程中正号时电路。图中的电阻R的选取是使运算放大器正向饱和电流为1mA，其它电阻均为1kΩ，电容的单位为微法。第四十三页，共五十页，编辑于2023年，星期一

对方程取负号，并设x项的系数是可调参数B，则方程变为：此图是对不同参数B时，在平面内的吸引子形式。2.微分方程模拟电子电路跃变非线性函数电路

对方程的不同正负号，奇怪吸引子有不相的形式。当取正号时，奇怪吸引子是一种单折带形式，有点象罗斯勒吸引子；当取负号时，奇怪吸引子是一种双旋结构，类似于洛仑兹吸引子。第四十四页，共五十页，编辑于2023年，星期一

这是一个具有非线性电阻的混沌电路电路，是由美籍华人蔡少棠首先发起研究的。图中电阻R是非线