理论力学7-非线性动力学与混沌-讲义...ppt

理论力学7-非线性动力学与混沌-讲义...ppt

《理论力学7-非线性动力学与混沌-讲义...ppt》由会员分享，可在线阅读，更多相关《理论力学7-非线性动力学与混沌-讲义...ppt（106页珍藏版）》请在冰豆网上搜索。

第七章.非线性动力学与混沌Chapter7.NonlinearDynamicsandChaos,宋若龙吉林大学物理学院,参考书,刘秉正，非线性动力学与混沌基础，东北师范大学出版社，1994林振山，非线性力学与大气科学，南京大学出版社，1993刘式达，刘式适，非线性动力学和复杂现象，气象出版社，1989,7.1引言,一.“非线性动力学”的表观含义,线性非线性非线性,定义：

运动微分方程含有坐标或速度的非线性项的系统，称为非线性动力学系统，反之称为线性动力学系统。

例：

二.决定性系统与不可预测性,存在且唯一，可预测性,1.力学决定论及其伟大成就,设想一位智者在某一瞬间得知激励大自然所有力及组成它的物体的相互位置，如果这位智者又能对众多的数据进行分析，把宇宙间最庞大的物体和最轻微的原子的运动凝聚在一个公式中，没有什么事物是不确定的，将来就像过去一样清晰地展现在眼前。

Laplace，法国数学家，（1749-1827）,1757年，哈雷慧星（Hallycomet）按预测回归。

1846年，海王星在预言的位置被发现。

日月蚀的准确预测，宇宙探测器的成功发射与回收。

广义相对论，量子力学也是决定论的。

0.2352680.235,2.力学决定论不断受到挑战,1883年，英国流体力学家雷诺（Reynolds），湍流实验。

（烟）1903年，法国数学家昂利庞伽莱（HenriPoincare），三体问题，不存在统一的第一积分，混沌。

1963，美国气象学家洛仑兹（Lorenz），天气预报，“蝴蝶效应”：

巴西热带雨林中一只蝴蝶扇一下翅膀，两个星期后，就可能在美国得克萨斯州引起一场龙卷风。

洛仑兹方程,初值敏感演示,Duffing方程：

（带阻尼弹性系统的强迫振动）,初值敏感性,不可预测性，混沌,不可预测性=客观世界的非决定论？

线性系统是特殊的、近似的非线性系统是普遍的、本质的（ex:

弹簧、单摆）,振动、流体力学、声学、光学、气象学、天文学化学、生命学、生态学经济学、金融学、社会学,混沌现象是矛盾的结合体,决定性与随机性稳定与不稳定有序和无序,三.常微分方程的一般形式,1.自治方程与非自治方程,不显含时间，自治的显含时间，非自治的,2.常微分方程一般形式,

（1）自治的,2阶，1维,1阶，2维,

（2）非自治的,n维非自治,n+1维自治,例1：

Duffing方程,一阶常微分方程组,数值计算系统的状态相空间,优点：

四.相空间（相图）,相空间，也就是状态空间，是由广义坐标和广义动量（速度）张成的空间，也称相宇。

相空间中运动状态的变化轨迹称为相图。

弹簧振子,相图,时空轨迹,阻尼弹簧振子,通解,代入方程,当阻尼为正阻尼且很小时,阻尼弹簧振子,时空轨迹,7.2运动稳定性分析,一.非线性方程解的各种形式,1.定态解,平衡点，奇点,2.发散解,之一或几个随时间无限地偏离初值,爆炸，散射,3.振荡解,既不趋于无穷大，也不终止于某一点，而是在一定区域内不断变化。

周期振荡混沌,相轨迹没有确定的形状周期、貌似随机的运动。

闭合曲线,非闭合曲线,准周期振荡,F=0.1,F=0.29,F=0.32,二.解的稳定性,Lyapunov稳定性定义：

（1）设t=t0时方程的解为，t时为，另一受扰动而偏离它的解t0时为，t时为。

如果对于任意小的数，总有一小数存在，使得当时，必有则称解是Lyapunov意义下稳定的，简称Lyapunov稳定的或稳定的。

两矢量间的距离,

（2）如果解是稳定的，且则称此解是渐进稳定的。

（3）不满足上述条件的解是不稳定的。

例2.,解：

是Lyapunov稳定的,例3.,解：

渐进稳定的,三.线性稳定性分析,1.线性稳定性定理,设为方程的一个解（参考解），则为研究该解的稳定性，令为此解附件另一解，称扰动解。

若线性化方程的原点是渐进稳定的，则原非线性方程的参考态是渐进稳定的；若线性化方程的原点是不稳定的，则原非线性方程的参考态是不稳定的。

Lyapunov间接法,非线性方程组在参考态附近的线性化方程组,2.线性化方程组的解及其稳定性,试探解：

系数矩阵的迹,系数行列式的值,特征根,

（1）两特征根实部都是负的,参考态也是渐进稳定的。

是渐进稳定的,

（2）两特征根中至少有一个实部为正,是不稳定的,参考态也是不稳定的。

（3）两特征根中至少有一个实部为零，另一个实部为负,是Lyapunov稳定的,参考态处于临界情况。

奇点（平衡点，定态）的分类（取非线性方程的奇点为参考态）,

（1）,两根都是实的，且符号相同，此时奇点称为结点。

不稳定的结点,稳定的结点,

（2）,两根都是复的，此时奇点称为焦点。

不稳定的焦点,稳定的焦点,（3）,两根都是纯虚数，解是等幅振荡，此时奇点称为中心。

鞍点,（4）,两根都是实数，一正一负，此时奇点称为鞍点。

例4：

分析阻尼单摆定态的稳定性,解：

令,求定态解,两奇点,1.在奇点（0,0）处线性化方程组为,奇点（0,0）为结点,（过阻尼）,稳定的结点,不稳定的结点,稳定的焦点,不稳定的焦点,奇点（0,0）为焦点,（欠阻尼）,奇点（0,0）为中心,（无阻尼）,2.在奇点处线性化方程组为,奇点为鞍点,线性稳定性定理只适用于分析非线性方程奇点及其附近的解的性质，离奇点越远，线性化误差越大。

7.3极限环渐进稳定的周期振荡,一.定义,相空间里孤立的闭曲线，称为极限环,与初始条件有关的周期振荡不是极限环,极限环,与初始条件无关,此轨道极小邻域内不出现其它闭轨道,例5：

VanderPol方程（电子管振荡）,阻尼力与速度同向，负阻尼，对系统供能，振幅逐渐增大，振幅终将大于1。

阻尼力与速度反向，正阻尼，消耗能量，振幅逐渐减小，振幅只能等于1。

F=0,a=0.2,x0=4,F=0,a=0.2,x0=0.5,二.极限环存在的判据,庞伽莱-班狄克生判据（Poincare-Bendixsontheorem）:

有一解的相轨迹总是局限于相平面中不包含任何奇点的有限区域D内，则此轨迹或者是一极限环，或者趋于一极限环。

如果方程,（二维自治系统）,三.极限环的稳定性,稳定环,不稳环,半稳环,如果从包含极限环L的环形域（L的内侧和外侧）出发的任何轨线在时都渐近地趋于该极限环，则称极限环L是稳定的，否则称为不稳定的。

如果从包含L的环域内L的某一侧出发的轨线在时都渐近地逼近L，而从另一侧出发的轨线都远离L，则称L是半稳定的。

半稳定的极限环是不稳定极限环的一种。

例6：

求非线性系统,的极限环性解及其稳定性，c为参数。

解：

令,微分得,代入方程得,联立，令等式两侧的系数分别相等，得极坐标下方程：

在极坐标中系统相轨迹以常角速度旋转，由可求平衡态为：

奇点,极限环（为实数时）,为复数，只有平衡态,为稳定的焦点。

有两个平衡态,为半稳环（不稳环）。

为稳定的焦点，,有三个平衡态,为稳定的焦点，为不稳定极限环，为稳定极限环。

（心脏）,有两个平衡态,为不稳定的焦点，为稳定极限环。

四.极限环的特点,非线性系统周期振荡独有的特征；极限环在相空间中是孤立的；由系统的固有性质（运动方程及其参数）决定，与初始状态无关；包围不稳定奇点的极限环一定是稳定的，而包围稳定奇点的极限环一定是不稳定的；极限环只能包围结点和焦点，而不能包围鞍点。

Homework:

1.用线性稳定性定理讨论中心力场中圆轨道的稳定性。

2.求解如下常微分方程组的定态解、极限环型解，分析其稳定性，若有分岔现象，说明其分岔的类型。

3.用摄动方法求至1级近似解,7.4含弱非线性作用的一维振动摄动方法,一.无阻尼、无强迫力的一维弱非线性振动,为弱非线性作用,无因次化,摄动方法，设解为：

零级解,一级解,二级解,代入方程,的同次项相等,零级解方程,一级解方程,二级解方程,零级解方程为简谐振动方程，其解为,由得各级解初始条件为,可得零级解为,将零级解代入一级解方程,伪共振,非线性项导致系统固有频率改变,小量，可正可负,均为小量，可令,代回到原运动微分方程,将代入得,零级解为,一级解满足的方程,为避免伪共振，必有,一级方程变为,设特解,一级方程齐次方程通解可写为,非齐次方程解为,齐次方程通解,非齐次方程特解,把代入二级解方程，可得二级解。

当仅求至一级解时，非线性方程的解为,弱非线性作用下非线性振动的特点：

固有振动的频率由变为，且改变量与振幅a有关；整个振动除基频外，还有谐频，当进一步顾及高级近似解时，还有出现等奇数倍高次谐频振动；可推当非线性作用力为时会出现等偶数倍谐频振动；系统本来不受强迫力，但一级解满足的方程出现了强迫力，并且是3倍频的，这是由于非线性振动引起的。

7.5非线性强迫振动振幅破裂,Duffing方程,假定为小量，设试探解为,将试探解代入方程，仅保留至的一次项,利用关系式,令方程两端线性无关项的系数分别相等，可等待定系数满足的方程：

振幅A（近似为系统的振幅）随驱动频率的变化,当时,当时（考虑非线性），用数值方法求解，画出振幅频率响应曲线：

A,A,A,A,在段，同一频率下，振幅出现多值现象，CD段表示不稳定振动。

驱动频率逐渐增大或减小时，出现振幅跳跃（振幅破裂）现象。

例：

洗衣机甩干过程机械（汽车、飞机）,分谐振、组合频率谐振,取特解,小量,代入方程并保留到的一阶量,令方程两端线性无关项的系数分别相等：

分谐振,非线性系统受两个不同频率的外力同时作用时，系统除了以主要的频率振动外，还包含有频率为等组合谐振成份，若非线性项不太弱需要考虑高阶项时，振动将包含各种频率为的成份。

即谐频，则称组合频。

例：

耳膜。

当时，B有实根，特解存在，出现频率为的分谐振，也称为分频，且为主要振动。

例：

石英钟。

非线性受迫振动的特点,驱动频率连续变化时出现振幅跳跃现象。

驱动频率在某值处的微小改变，系统振幅发生剧烈变化。

谐频振动：

基频（驱动频率）为时，当非线性项为x的奇次幂时，会出现等奇数倍谐频；当非线性项为x的偶次幂时，会出现等偶数倍谐频。

当驱动频率远大于系统固有频率时（），会出现分频，也称为倍周期。

x的奇次幂，x的偶次幂。

当强迫力为两不同频率时，有组合频出现。

如耳膜。

7.6亥姆霍兹木马（Helmholtzcarousel）,2013年第26届国际青年物理学家竞赛IYPT题目,图片资料来源于https:

/,参考文献：

R.R.Boullusa,etal.,ThereactionforceonaHelmholtzresonatordrivenathighsoundpressureamplitudes,Am.J.Phys.,60（8）,pp722-726,1992,空气柱输入阻抗,动态粘滞系数,开口辐射阻抗,内口辐射阻抗,颈部空气柱总阻抗，实部为阻（能量损失），虚部为抗（等效质量）,=,m,等效质量,阻尼系数,若振动频率非常低使得可看作绝热过程：

比热容,Taylor展开,密闭空腔产生的恢复力,为小量,设,可解得,一级解将有的谐频振动。

当时振幅最大,密闭空腔对颈部气柱的恢复力,净力（一个周期内的平均）,7.7分岔（Bifurcation）,一.分岔的概念,1.定义：

对常微分方程组,为参数。

如果参数在某一值附近的微小变化将引起解的性质（相轨线的拓扑结构）发生突变，则此现象称为分岔。

称为临界值或分岔值。

在坐标轴上其对应点为分岔点。

例：

极限环求解,2.解的结构稳定性,指在参数发生微小变化时解的轨线仍维持在原轨线某一邻域内。

因次非线性系统在常点的解具有结构稳定性，而分岔点附近的解是结构不稳定的。

二.分岔的类型,1.叉式分岔,系统参数发生微小变化时，一个稳定的定态两个稳定的定态,例1：

水平滑动摆，弹簧原长l，参数a变化,定态,1个奇点,3个奇点,

（1）奇点,（0,0）为中心，Lyapunov稳定的（0,0）为鞍点，不稳定的,线性稳定性定理：

两奇点均为中心，Lyapunov稳定的,

（2）奇点,2.霍普夫（Hopf）分岔,系统参数发生微小变化时，稳定的定态