最新二章分岔与奇怪吸引子
非线性电路与混沌讲解
混沌运动的主要特征
初值敏感性
长期行为的不可预见性
（例：洛仑兹的天气模型）
下一页
看两条轨道是如何分道扬镳的 （初值分别为0.506，0.506127）
倍周期分岔
解释倍周期分岔现象，我们从混沌描述中最重要的的一维 非线性迭代方程式入手。这类方程中最有典型意义的是虫口方程。
xn1 xn (1 xn )
生命的整体性，包括： （1）时间上整体，即生命活动的高度有序性；
（2）空间上整体，即生命结构和物质相互作用、相互影响形成网 络状整体。如形态上，人体结构和功能混沌调节机制（如免疫网 络调控、心脏、肺、肠的分形结构学原理，心电的混沌产生与心 脏普肯野氏纤维分形分布联系等等。在物质信号相互作用的动力 学研究上，其相互作用不仅仅是激活、失活或抑制、促进这一简 单的关系，而应包括复杂的数学过程，这种复杂的数学过程应是 非线性的，很可能符合混沌原理。这一现象至今尚未有人意识到 这一点，但将来肯定会的！当然，想彻底了解生命意义，这条路 还很长很长。Biblioteka 实验原理GL
C2
C1
R
R 有源非线性负阻元件, G 电导，C1和C2 电容。
各区域的作用：黄色区： 产生振荡，蓝色区： 移相， 粉色 区：有源非线性元件。 下一页
～
有源非线性元件的电压、电流特性：
I
0
V
上图是一个非线性负阻的电压电流特性曲线，它的 实现方法有许多种，本实验使用的是Kennedy于1993年提 出的方法。它采用了两个运算放大器和六个配置电阻来 实现的。由于我们主要研究的是元件的外部效应，即电 路两端的电压和流过其电流的关系，因此我们可以把上 述元件看成是一个黑匣子即—有源非线性负阻。 下一页
在数学上把天气(气候)预报问题提成初值问题，即用 动力学的方法进行预报，从认识论上讲就是把大气看成是确 定论的系统，这在较短的时间尺度内是行得通的，而在时间 较长的时候却是有问题的，主要是大气运动是非线性、强迫 和耗散的。由这三大特点，可以得到一幅这样的图像：误差 是随着时间呈指数增加的，初始场的作用随着时间是衰减的， 必须考虑能量的补充和耗散。Lorenz发现了“蝴蝶效应”， 指的就是初始场微小的不确定性的指数放大。这就提出了确 定论预报的可预报性问题，中期数值天气预报逐日预报的可 预报时限大约是两三周左右的时间。 也就是说进行长期预报是不可能的。
非线性物理学概论
非线性科学概要为《非线性物理概论》一书写的序言汪秉宏上一世纪初量子力学和相对论的发现，因为提出了突破人们传统思维的新概念，将人类的世界观推进到超越经典的领域，而被公认为是物理学或更确切地说是科学的两次革命。
牛顿创立的经典力学被发现并不始终是正确的。
当深入到微观尺度（<10-8cm），应该取代为量子力学，当物体的速度接近于光速（～10 10cm/s），则相对论是正确的。
非线性科学作为科学的一个新分支，如同量子力学和相对论一样，也将我们引向全新的思想，给予我们惊人的结果。
非线性科学的诞生，进一步宣布了牛顿的经典决定论的局限性。
它指出，即使是通常的宏观尺度和一般物体的运动速度，经典决定论也不适用于非线性系统的混沌轨道的行为分析。
非线性科学涵盖各种各样尺度的系统，涉及以任意速率运动的对象，这一事实丝毫不降低这一新学科的创新性，恰恰相反，刚好说明它具有广泛的应用性。
从这一点来看，其实非线性科学的诞生和发展更有资格被称为科学的一场革命。
非线性科学，目前有六个主要研究领域，即：混沌、分形、模式形成、孤立子、元胞自动机，和复杂系统。
而构筑多种多样学科的共同主题乃是所研究系统的非线性。
一个系统，如果其输出不与其输入成正比，则它是非线性的。
例如一个介电晶体，当其输出光强不再与输入光强成正比，就成为非线性介电晶体。
例如弹簧，当其位移变得很大时，胡克定律就失效，弹簧变为非线性振子。
又例如单摆，仅当其角位移很小时，行为才是线性的。
实际上，自然科学或社会科学中的几乎所有已知系统，当输入足够大时，都是非线性的。
因此，非线性系统远比线性系统多得多，客观世界本来就是非线性的，线性只是一种近似。
任何系统在线性区和非线性区的行为之间存在显着的定性上的差别。
例如单摆的振荡周期在线性区不依赖于振幅，但在非线性区，单摆的振荡周期是随振幅而变的。
从数学上看，非线性系统的特征是迭加原理不再成立。
迭加原理是指描述系统的方程的两个解之和仍为其解。
分岔图做法[1]
![分岔图做法[1]](https://img.taocdn.com/s3/m/f8c583b4172ded630b1cb664.png)
>>混沌研究总结篇------一、分岔图（系统）先打个提纲，这几天把自己混沌相关知识研究学习内容总结一下。
首先简绍几个基本概念：一、自治系统一个n阶自治的连续动态系统可以表示为可以理解为对于自治的连续系统，上相量场f是不依赖于时间t的。
二、非自治系统一个n阶非自治的连续动态系统可以表示为可以理解为对于非自治的连续系统，向量场f不仅依赖于状态变量x，而且依赖于时间t，如Duffing振子。
三、庞加莱映射庞加莱映射是一个传统的用来离散化连续系统的方法。
庞加莱映射可以用(n-1)阶的离散映射来取代n阶的连续系统。
庞加莱映射的用处正在于减小系统的阶数，并且在连续系统和离散系统之间建立了一座桥梁。
对于n阶自治系统，其对应的解对就着轨迹。
当选择作为一个(n-1)维的超平面，这样轨迹将穿越超平面。
难点主要是超平面的选取，使其对应的解穿越超平面，就可以得到一个领域内的庞加莱映射。
对于n阶非自治系统，若其外加强迫力的最小周期是T,j最终的庞加莱映射可以定义为相应的轨道P(xk)是对某个轨迹每隔T时刻采样一次获得，这种操作和每隔T时刻的频闪观测仪的行为很相似。
所以要想得到一个系统的庞加莱映射，这段话一定要好好理解，当真真知道这中间说的含义，庞加莱映射这么画其实也已经知道国。
四、分岔图分岔图的横坐标是一个变化的参数，纵坐标是你要求的某一个量的随着各参数的变化情况，而poincare则是我们选取横坐标上的某参数的某一个具体值时截面图，只不过poincare截面的选取其实可以是任意的。
下面主要研究的混沌系统有：Logistic、Henon、Lorenz、Duffing、Rossler、Chen、混沌电机模型等系统系统先说Chen系统，因为和课题有一定的关系，而且自己以后起家也得从Chen系统入手。
系统方程如下：dx/dt=a*(y-x)dy/dt=(c-a)*x+c*y-x*zdz/dt=x*y-b*z就是对此方程中不同参数a、b、c下对系统画分岔图，研究混沌系统(1)给定a、c，画b关于系统的分岔图结果如下图所示CODE：function fenchatuchenclc;clearXA=35;XC=28;Z=[];for XB=linspace(2,,100);options = odeset('RelTol',1e-6,'AbsTol',[1e-4 1e-4 1e-5]);[T,X]=ode45('chen',[0,50],[-5 0 5],options,XA,XB,XC);n=length(X);for k=round(n/2):nif abs(X(k,1))<1Z=[Z,XB+abs(X(k,2))*i];endendendfigureplot(Z,'.','markersize',1)title('chen映射分岔图')xlabel('b'),ylabel('|x| where x=0')这组代码不完全是自己的,现在见解其中一些方法在进行自己系统的绘制，这个程序的具体原理我会在后面给出来的。
分岔与奇怪吸引子
d2t
dt
对于平衡点 I1 邻域有:
I(t)I0exp (t)
引进参数作用 量I 与角度量q
x 2I cosq
d d It 2 Ico 2qs1 2Isi2qn
I0 是初始对 I1 的偏离小量。作用量I
随时间指数增长, I1是不稳定解, 为不 稳定焦点。
相位求平均
/ C2
dI dt
I2 C
1. 切分岔
数学模型
利用方程： 由
dx x2
dt
得d平x/d衡t 点0
x0
(a)当μ＜0时,解 x0 为虚数，因此不存在奇点， (b)当μ＞0时出现两个奇点, x0 ，
说明上述方程的解在 x0=0 处发生了分裂。 μ＞0 两个奇点的稳定性
在解 x0 附近取一点，计算它与平衡点距离随
时间变化。设距离：x x0
向不动点B
.
2 .平方映射的不动点
不动点的稳定性
非线性动力学核心问题之一就是研究系统的稳定性问题。
上述计算可见，当μ<3时迭代走向不动点，当μ>3迭代值出现持续振荡，
说明迭代在μ= 3附近发生了变化，稳定不动点变得不稳定了。
如一维映射 xn+1f具(有,x不n)动点，即有解
x f(,x)
设 n 为对不动点的偏离量，需继续迭代，有： xn + 1f(,xn)
.
2 .平方映射的不动点
平方映射的不动点
通过作图或数值计算表明，计算可以得到一个不变的终值，它被称为映 射的不动点。一个映射的不动点就是xi与xi+1相同时的数值，它不再因继续 迭代而发生变化。对平方映射，不动点为：
解此方程得：
xi xi(1xi)
近代物理实验思考题答案
近代物理实验思考题答案一、夫兰克—赫兹实验1解释曲线I p －V G2形成的原因答；充汞的夫兰克-赫兹管，其阴极K 被灯丝H 加热，发射电子。
电子在K 和栅极G 之间被加速电压KG U 加速而获得能量，并与汞原子碰撞，栅极与板极A 之间加反向拒斥电压GA U ，只有穿过栅极后仍有较大动能的电子，才能克服拒斥电场作用，到达板极形成板流A I 。
2实验中，取不同的减速电压V p 时，曲线I p －V G2应有何变化？为什么？答；减速电压增大时，在相同的条件下到达极板的电子所需的动能就越大，一些在较小的拒斥电压下能到达极板的电子在拒斥电压升高后就不能到达极板了。
总的来说到达极板的电子数减小，因此极板电流减小。
3实验中，取不同的灯丝电压V f 时，曲线I p －V G2应有何变化？为什么？答；灯丝电压变大导致灯丝实际功率变大，灯丝的温度升高，从而在其他参数不变得情况下，单位时间到达极板的电子数增加，从而极板电流增大。
灯丝电压不能过高或过低。
因为灯丝电压的高低，确定了阴极的工作温度，按照热电子发射的规律，影响阴极热电子的发射能力。
灯丝电位低，阴极的发射电子的能力减小，使得在碰撞区与汞原子相碰撞的电子减少，从而使板极A 所检测到的电流减小，给检测带来困难，从而致使A GK I U 曲线的分辨率下降；灯丝电压高，按照上面的分析，灯丝电压的提高能提高电流的分辨率。
但灯丝电压高, 致使阴极的热电子发射能力增加，同时电子的初速增大，引起逃逸电子增多，相邻峰、谷值的差值却减小了。
二、塞曼效应1、什么叫塞曼效应，磁场为何可使谱线分裂？答；若光源放在足够强的磁场中时，原来的一条光谱线分裂成几条光谱线，分裂的谱线成分是偏振的，分裂的条数随能级的类别而不同。
后人称此现象为塞曼效应。
原子中电子的轨道磁矩和自旋磁矩合成为原子的总磁矩。
总磁矩在磁场中受到力矩的作用而绕磁场方向旋进从而可以使谱线分离2、叙述各光学器件在实验中各起什么作用?答；略3、如何判断F-P 标准具已调好？答；实验时当眼睛上下左右移动时候，圆环无吞吐现象时说明F-P 标准具的两反射面平行了。
分岔ppt课件
2 .平方映射的不动点
<1时走向不动点 A
当参数<1时，抛物线高度较低，与迭代线只有一个交点A。这时不管
初值如何，迭代最终趋于原点，原点是唯一的不动点。图b是随迭代次 数 n 的变化曲线，这是最终衰变到零的指数衰变曲线。在生态上，虽 然初始有一定的种群数量，但受到环境的制约最终走向了灭绝。
2 .平方映射的不动点
准备:
1. 建立坐标系 x n1 ~ x n
2. 作条抛物线：
x n1 x n (1 x n )
3. 作对角线，称恒等线
x n1 x n
通过它做投影。
1.平方映射
作图计算
平方映射 x n1 x n (1 x n ) 在 x n1 ~ x n 平面上是一条抛物线， 抛物线高度由 值决定。
d 2x dt 2
k x
x3
0
由势能曲线知：
a. 在 k 时0仅有一个平衡点：
x 0
b.在 k时存0 在三个平衡点：
x0 x k
可见在参数 k = 0 处发生了一次从单解
转为三解的叉式分岔。
c.在这三个平衡点中， x ， 处k
在势能极小点，是稳定的； x 处0 在势
能极大点，是不稳定的平衡点。
4 霍夫型分岔
一、从横坐标 x0 处作竖直线 与抛物线相交，这点的纵坐标 高度即为 x1;
二、从此点作水平线与对角线 相交，此交点横坐标即为x1；
三、再由此点作竖直线，得到 与抛物线相交时的高度x2，再 将x2移植到对角线上，找到横 坐标x2。从这里作竖直线与抛 物线相交得x3，如此反复······
1.平方映射
对于稳定的不动点，应有 e n+1 e n , 即： m 1
3.1-分岔
4.霍夫型分岔
分岔分析
1/(2t C) 0
/(1 Ce2t ) 0
参数μ从负变到正，从焦点
产生出极限环,这种分岔称
霍夫分岔。分岔点位于μ=0。
4.霍夫型分岔
极限环
Hopf 分岔有超临界和亚临界的区别
二 平方映射与倍周期分岔
1. 平方映射 2. 平方映射的不动点及其稳定性 3. 平方映射的周期解
第 n 代数量: Nn 第 n+1 代数量: Nn1
A 如不考虑生存环境对种群生存的影响，第 n 代与第 n+1代有
如下关系：
Nn1 RN n
当 R >1，种群数量将线性地无限制增长。
1．平方映射
平方映射导出—生态平衡方程
B 种群受环境制约，数量有最大限额 N0 ，种群繁殖空间 N 0 N n ，第 n 代与第 n+1代关系
数学模型
dx
dt dy
y x[ (x 2 x y[ (x 2
y 2 )] y 2 )]
dt
引入极坐标
x2 y2
x cos
y
sin
求导
d
dt
(
2
)
d 1
dt
代入原方程
令正弦余弦系数 相等
dx
dt dy
d
dt
d
cos cos
sin cos
dt dt
Nn1 KNn (N0 Nn )
N n1
N n (1
Nn N0
)
K N0
xn1 xn (1 xn )
xn Nn / N0
1.平方映射
平方映射计算
方程展开 xn1 xn (1 xn ) xn xn2
fx2-1
y = f ( , x) 这里 为系统参数。设系统状态作等间隔 t，t+1，t+2，t+3，…变化，则时 间演化方程改写为： x (t + 1) = f ( , x (t ))
当时间间隔不取整数，各时刻写成 t1 = t0 + t , t2 = t0 + 2t ,L, tn = t0 + nt 相应的状态为：
d 2x dx + ε x2 1 +ω2x = 0 dt dt 2
(
)
第二节 平方映射与倍周期分岔
1. 平方映射 2. 平方映射的不动点及其稳定性 3. 平方映射的周期解及其稳定性 4. 倍周期分岔的功率谱
1．平方映射 ．
映射方程
物理学上一个动力学系统可以用连续变量表示，也可以用离散数表示。 一个以为连续变量的单参数的动力学系统：
x 0 = 0.1
x 1 = x 0 (1 x 0 ) = 2.4 0.1(1 0.1) = 0.216
x2 = x1 (1 x1 ) = 0.40642L
x3 = 0.578985L x4 = 0.5859465L
x5 = 0.58227L
x6 = 0.583755L
L
在此参数下，计算结果趋向一个终值： x∞ = 0.583335L
忽略高 阶量
dξ dx = = (ξ + x 0 )2 + dt dt dξ = 2ξx 0 dt
1. 切分岔
解的稳定性与相流
解 dξ ξ (t ) = ξ 0 exp( 2x 0t ) = 2ξx 0 dt ，当 t → ∞ 时， → 0 ，此解是稳定的，是稳定的结点 ξ 结点。 结点 ξ ，当 t → ∞ 时， → ∞ ，解是不稳定的，它是鞍点 鞍点。 鞍点
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
可见在参数 k = 0 处发生了一次从单解
转为三解的叉式分岔。
c.在这三个平衡点中， x ，k处 在势能极小点，是稳定的； x 0处在 势能极大点，是不稳定的平衡点。
4 霍夫型分岔
数学模型
dddxtyxyyx[[((xx22yy22))]]
dt
引入极坐标
x2 y2
x cos
• 设： x12 .4,x0 x(01 各 0 次x.0 1)计…算x值2 为x ：1(1x1)
x0 0.1 x 1 x 0 ( 1 x 0 ) 2 .4 0 .1 ( 1 0 .1 ) 0 .21 x26 x1(1x1)0.40 642
d2t
dt
第二节 平方映射与倍周期分岔
1. 平方映射 2. 平方映射的不动点及其稳定性 3. 平方映射的周期解及其稳定性 4. 倍周期分岔的功率谱
1．平方映射
映射方程
物理学上一个动力学系统可以用连续变量表示，也可以用离散数表示。 一个以为连续变量的单参数的动力学系统：
yf(,x) 这里 为系统参数。设系统状态作等间隔 t，t+1，t+2，t+3，…变化，则时
1．平方映射
映射方程计算
对一个映射 xn1f(,xn)
的计算采用的是迭代方法。即给定一个初值 x 0 将其代入映射计算得 x 1 ，
将 x 1 代入映射计算得 x 2 ，由 x 2 可算得 x 3 ，如此一直计算得：x n
例如： 一个简单映射 xn1 Axn 1 次迭代： x1 Ax0 2 次迭代： x2 A1 xA (A0)x A 2 x0 n 次迭代： xn A A (A0 )x A n x0
于是有：
x1,x2,x3,xn
如果将 x i 值看成为一条线上的一个点，则该组数值就构成一条轨道。
1．平方映射
映射与微分方程对应关系
动力学系统用连续变量表示为微分方程，离散数表示时为映射(map)，两者
对应关系为：
xn1 Axn
dx Ax dt
迭代计算 xn Anx0
解方程 x x0eAt
1．平方映射
平方映射导出—生态平衡方程
• 1838年，生物学家伏埃胡斯脱（Verhulst）在研究生物种群演化时提出一种 设想：一个世代交替的生物种群是在一个受制约的环境中生息繁衍的。
第 n 代有: N n 第 n+1 代有: N n 1 A 如不考虑生存环境对种群生存的影响，第 n 代与第 n+1代有如下关系：
间演化方程改写为：
x(t1 )f(,x(t))
当时间间隔不取整数，各时刻写成 t 1 t 0 t ,t 2 t 0 2 t , ,t n t 0 n t
相应的状态为：
x1,x2,xn
xn x(tn)
时间演化方程变成离散方程：
xn1f(,xn)
数学上称为映射的方程。在非线性发展史上第一个将映射方程用于研究系 统进入混沌状态的是美国科学家梅(May Robert)
y
sin
求导
d dt
Hale Waihona Puke (2)代入原方程
d 1 dt
令正弦余弦系数相等
dx ddyt
d cos
dt
d cos
sin cos
dt dt
4.霍夫型分岔
分岔分析
对方程
积分,可得：
d dt
(
2)
d
dt
1
1/(2t C) 0 /(1Ce2t ) 0 t t0
d d22xtx21d dxt2x0
对于平衡点 I1 邻域有:
I(t)I0exp t()
引进参数作用 量I 与角度量q
x 2I cosq
d d It 2 Ico 2qs1 2Isi2qn
相位求平均
/ C2
dI dt
I2 C
I
dI/dt0
平衡点： I1 0
I2 C 2
I0 是初始对 I1 的偏离小量。作用量I
二章分岔与奇怪吸引子
分岔与奇怪吸引子
第一节 简单数学分岔 第二节 平方映射与倍周期分岔 第三节 流体不稳定性与洛伦兹方程 第四节 李雅普诺夫指数与奇怪吸引子
3 叉式分岔
数学模型
利用方程： dxxx3
dt
由 dx/dt 得0平衡点
x0 0
x
0
分岔图形象一把叉子，故称岔式分岔。 解的稳定性：
Nn1 RNn 当 R >1，种群数量将线性地无限制增长。
B 种群受环境制约，数量有最大限额 N 0 ，种群繁殖空间 N0 Nn
第 n 代与第 n+1代关系
N n 1RnN (N 0N n)
Nn1 Nn(1N Nn0)
RN0
xn1 xn(1xn)
xn Nn/N0
1.平方映射
平方映射计算
方程展开
xn1xn(1xn)
C，t0 为积分常数。
参数μ从负变到正，从焦点产生出
极限环,这种分岔称霍夫分岔。分
岔点位于μ=0。
1.μ≤0，距离 随时间而缩短，当时间 t 时 。 说明0μ轴线上
各点是稳定的焦点。
2. μ＞0， 值随时间增长，不论初始 的大小；当 t 时
形成闭合圈即极限环
4.霍夫型分岔
范德玻耳方程分岔
值情况，即：
如果 为正值，相平面上坐标原点
是不稳定的焦点，而极限环是稳定的。 不论初始相点处于环内还是环外 ，t时总是趋向于极限环。
如果 为负值，情况刚好相反，坐
标原点变为稳定的焦点，为系统的不 动点，而极限环则是不稳定的。当
t 时，环内相点趋于不动点，环外 相点则远离环而去。
d2xx21dx2x0
随时间指数增长, I1是不稳定解, 为不 稳定焦点。
对于平衡点 I2 邻域有:
I(t) I0ex p t()
I为0 初始对I2 的偏离量。作用量
I 对的偏离量 随时间指数减小。当
t，I0 , I C, I2 是稳定
的解。
4.霍夫型分岔
范德玻耳方程分岔
结论 范德玻耳方程霍夫型分岔与参数
的 正负有关。上面讨论的是 为正
μ＜0时只有 x0= 0 的平衡点，经分析
方法可知它是稳定的。
相流图形
μ＞0有三个平衡点， x0= 0 是不稳定 的，解 x0 是稳定的。
3 叉式分岔
杜芬方程的叉式分岔
杜芬方程具有叉式分岔
d2x d t2
kx
x3
0
由势能曲线知：
a. 在 k 时0仅有一个平衡点： x 0 b.在 k时0存在三个平衡点：
x n 1x n (1 x n )x nx n 2
• xn+1 值与 xn 值是平方关系，称平方映射，文献中称洛吉斯蒂映射 (logistic
map), 该式是抛物线表示式，也称抛物映射。
• 由于亲、子两代种群数约化值，在0 ~ 1间，参数μ取值在[0，4]内。
• 离散映射采用迭代计算。即给定参数 值与初始值 x0 ，就有：