单摆混沌现象研究 ——寻找周期解与混沌解.ppt

1.1492往后出现的有界相图变成形状奇特的曲线。这些就是奇异吸引子。 跳出混沌以后，相图会在角位移上无界。仔细研究这些无界的图形会发现一些有趣的事情 三张相图都是单方向旋转，但是相图的振荡有不同的振幅数。 第一张A=1.35是一个； 第二张A=1.45有两个； 第三张A=1.47有四个。 相图虽然没有了极限环，但是行为却有了周期性。 往后的相图又出现了收敛的极限环，以及更复杂的混沌解。 四、往后的情况 再画出A2.0的情况来探讨其规律性 单摆相图变化的规律大概就是： 单摆混沌现象研究 ——寻找周期解与混沌解 07～08 学年力学课电子作业 06级 4班 0610736 曾昱 单摆是在悬挂的细线的另一端连接着一个小球（如图所示）。单摆又称数学摆，是物理学中最简单的模型之一。 一、单摆的物理模型 可以认为，细线的质量可以忽略，且是刚性的。系统质量集中在可视为质点的小球上。设摆长为l，小球质量为m，相对于平衡的下垂位置的角度为θ，重力加速度为g。在无阻尼无驱动情形其运动方程为 有阻尼无驱动情形 有阻尼有驱动情形 对下面的非线性微分方程进行一些化减 Q称作质量因子 初值条件为： 并取： 当?从0.5到2.0变化时，单摆的相图也会相应地发生一系列变化。 开始的一段暂态过程过去后，周期运动的相轨会趋于闭合曲线，即极限环。中间将会出现一些分岔周期解。 后来随着?的增大，会出现混沌运动的相轨。 下面来探讨出现周期分岔和混沌的?的取值。 二、数值方法解微分方程及相图的绘制 下面的通过matlab6.5解微分方程并用其自带的绘图程序绘制的相图 （1）无阻尼无驱动情形： （2）有阻尼无驱动情形： （3）有阻尼有驱动情形 随着?的取值变化下面画出相图，每个相图相差0.02. (1)A=[0.5,1.0] Amg 是的相图都收敛在一个规则的极限环内。 （2）A=[1.00,1.48] [1.00,1.14]内似乎有较规则的极限环。 [1.16,1.26]内图形有界，但是“极限环”形状已经非常不规则，出现奇异吸引子。 [1.28,1.48]内相图呈振荡装，随时间变化角位移是无界的。即单向旋转 （3）A=[1.50,1.98]: [1.50,1.66]内图像再度收敛， [1.56,1.66]内出现规则极限环。 [1.68,1.98]区间图像再度出现混沌解。 三、寻找周期解 可以发现，单摆的相图随驱动力大小的变化是非常复杂的。我们在上述几个区间内寻找出 用matlab的循环语句可以较方便的找出一个区间的相图。 此图为A=1.025是的相图。这张图的t值已经取得很大，但可以清楚的看到轨迹收敛到一个很狭窄的区域。 下面为了寻找周期分叉，我以0.001为单位绘图。 中间对每个图形检查，在A=1.07处发现二倍分岔，下图为A=1.072时的解，它的图形更为清晰。 继续往后画 … … 这些图中出现了一个奇怪的图形。观察会发现，以前的相图有类似的相图，只是这回乡相图收敛到这个图形上。其实我们可以把它看作时两个类似椭圆的叠加。这个图只有一个环，但是若看作两个椭圆状的叠加，是否还是可以说时二倍周期呢？ 就着样一直画下去，逐个图检查直到图形非常“混乱”，也就是混沌了也没有找见二倍以上的周期解。