单摆混沌性质的计算机方法的研究.pdfVIP

‘国际物理教育通讯》 NTPE专辑 单摆混沌性质的计算机方法研究 吴宏途吴中光 福建师范大学物理系福州350007 摘要：本文介绍用计算机数值模拟方法研究单摆的运动状态由确定论走向混沌的全过程及由此得到的 单摆混沌运动的性质。 关键词：计算机数值模拟相平面庞加莱截匾混沌单摆 一．引言 混沌是非线性系统在远离平衡状态下，由于运动轨迹高度不稳定经多次分岔达到的一种运 动状态，是～种有统计规律的序。哈肯曾指出：“混沌性来源于决定论性方程的无规则运动。” 那么确定论系统又是如何产生混沌的呢?由自然界中广泛存在的混沌现象，我们可以知道系统 产生混沌归咎于非线性条件。单摆是物理学中的一个简单模型，是确定论的典型例子，但在有 阻尼的情况下，受周期性驱动力作用之后，其运动规律就将发生变化，出现周期的分岔，在一 定大小的驱动力条件下出现混沌行为。为了将这一由确定论走向混沌的全过程完整地演示出 来，并研究其中的一些性质，我们可利用具有高速运算能力的计算机作为实现工具，借助微机 的数值模拟来实现。本文介绍了单摆运动周期变化过程的计算机模拟程序的设计思想及其运行 结果。 二．理论依据 单摆在有阻尼情况下，受周期性驱动力作用的运动规律可由以下方程组给出 02∞ (_)=一2B∞-sin0+fCOS巾 (1) r 单摆混沌性质的计算机方法研究 由=Q 式中0表示单摆的摆角，o表示单摆的瞬时角速度，B是无量纲的阻尼因子，f．Q都是无 量纲的驱动力参量；物理量上的圆点表示对时间的导数。 为了研究单摆的运动规律，我们引入相平面分析法，根据方程组(1)在0一。坐标平面上 作相图。我们发现单摆在有阻尼的情况下，受周期性驱动力作用之后，其相轨道首先有—暂态 过程，而后才进入我们所要研究的终态，并且对于不同的驱动力参量f，其相轨迹表现出不同 的终态：当f1。0时，相轨迹只有一条，但相图己由原先对称的交得不对称了；随着f的 继续增大，当f1．066时，相轨迹由一条演化为两条，即周期开始分岔，由一倍周期变为 二倍周期；当fL079时，周期再次分岔，．由二倍周期分为四倍周期；当f1．082时，周 期又发生了分岔。随着f值的不断增大，周期越分越多，最终出现混沌行为。单摆就是这样由 确定论进入混沌状态的。但随着f的继续增大，相轨迹又发生了变化，当f3．00之后，相 轨迹又变成了一条，只是此时的单摆不再是左右摆动，而是卓B着～个方向运动。当然，这一倍 周期也不稳定，随着f的增大，周期发生分岔进入二倍周期．四倍周期…．．．最终回到混沌状 态。 由以上的分析，我们不仅可以了解确定论是如何进入混沌状态的，而且还会发现在远离平 衡的情况下，从分岔．失稳到形成自组织与经多次分岔而进入混沌状态二者是可以相互转化的， 并且经常互相嵌套。 在研究单摆的复杂运动时，相图固然可以使运动的描绘简化，但在周期倍数很高时，其简化程 度就远远不够了。此时轨迹是银乱的，较难与非周期的运动区分开来。因而我们要引入～种新 的分析方法一一～庞加莱截面法。用此法可以较简单地描述运动的周弱性。庞自n莱截面法是在二 维相平面的基础上增加一个变量由，将坐标扩展蔓l三维，并且在巾方向上每隔2Ⅱ周期取一 个截面，这样我们馒可以从落在截面上的点的数目来判断周期数。当f；1．025时，截面上只 会出现一个点，而当f=1．15时，截面上出现～系列点，这正与单摆的运动状态相吻合。利 用这种截面法可以将连续的动力学系统简化为离散的动力学系统，使研究变得简单，周期数一 目了然。并且由实验可以发现在相同参量的条件下，不 同的