基于轨迹中点图的弹球系统混沌行为分析.docxVIP

基于轨迹中点图的弹球系统混沌行为分析 0 系统的轨迹中点图 布局系统反映了在封闭区域内传播辐射和其他物理过程的动态行为。作为一种通用的保守动力学系统，弹球系统的研究对力学、光学、声音等领域有着重要的意义。长期以来, 人们研究了各种形状的弹球系统, 如三角形[8,9,10,11,12,13,14]、多边形、正方形、矩形、圆形、半圆、1/4圆、椭圆、体育场形、柠檬形等等, 考察了系统轨迹的周期性、稳定性、遍历性、混沌性等多个方面的非线性动力学特性, 得到了许多有价值的结论。现有的研究中运用了理论分析、图解描述、解析模型数值计算等多种考察手法。本论文的研究将采用图解法。 典型的弹球系统是这样一个系统:一条射线在封闭区域内传播, 在到达边界之前, 射线在空间内沿直线传播, 遇到边界时, 遵照反射定律发生反射, 接着继续沿直线传播。射线在弹球系统中的整个传播轨迹为一折线, 在相邻2个反射点之间, 轨迹为一段直线。在对各种形状的弹球模型进行综合分析的过程中发现, 截取这些直线段的中点, 由这些中点组成的图形, 能够非常有效地反映弹球系统的动力学行为特性。这里将这种图形称为轨迹中点图。轨迹中点图的原理简单, 形象直观, 实现起来也很方便。只要利用几何学的方法, 借助于计算机程序, 得到系统的具体轨迹, 然后直接截取轨迹中每一直线段的中点, 当反射次数足够多时, 这些中点叠加在一起, 即构成完整的轨迹中点图。由轨迹中点图的形态, 就可以判断弹球系统的动力学特性。 1 设“左两端”边界 弹球系统中射线的一条轨迹由前后衔接的一系列直线段组成, 每一直线段可以用方位角和出射位置2个参量来描述。方位角φ为射线的前进方向与水平向右方向 (可以理解为x轴正向) 之间的夹角;出射位置可以用s来表示,s为从出射点所在的一段边界的左端点出发沿着该边界到达出射点的路程长度与该条边界总长度的比值 (对于圆形弹球系统, 由于边界不存在端点, 可以根据需要人为规定圆弧上某一点为其“左端点”) 。 基于上述规定, 弹球系统中射线的运动可以用以下形式的离散型动力学方程来描述: 式中sn为轨迹中第n个直线段的出射位置,φn为第n个直线段的方位角,n=0, 1, 2, …;整个轨迹的初始方位角用φ0表示, 初始出射位置用s0表示 (初始出射点若不在边界上, 则另行说明) 。函数Fs(·) 、Fφ(·) 的具体形式因弹球系统的具体模型而异, 给出每一种弹球模型的Fs(·) 、F?(·) 的具体形式均需要较大的篇幅, 所以本文将不给出Fs(·) 、Fφ(·) 的具体形式, 具体轨迹的计算将由计算机程序完成。本课题组此前撰写的论文中推导出了三角形、体育场模型的离散动力学方程的具体形式。 2 各种形式的弹簧系统轨道中心图 2.1 准周期的轨迹中点图 在考察一般的平行四边形之前, 可以先研究比较简单的矩形。图1 (a) 是矩形弹球系统中某一条射线的轨迹, 这里矩形长为1.5, 宽为1, 射线初始方位角φ0=62°, 初始出射点位于下底边且s0=0.2/1.5。对于矩形中这样一条初始方位角取非特殊值的射线, 射线传播途中其方位角φ只能取62°、118°、242°、298°等4个值, 轨迹由2组平行直线段构成 (φ=62°与φ=242°的线段平行, φ=118°与φ=298°的线段平行) , 前后各直线段不发生重叠, 并且随着传播无限进行下去, 轨迹将逐渐填满整个矩形, 这是一条准周期的轨迹。图1 (b) 是这一条轨迹的中点图, 它由若干条直线段组成。从这里可以看到, 直线型的轨迹中点图与弹球系统的准周期行为相对应。 研究表明, 当tanφ0=p/q (其中p、q为整数) 时, 正方形弹球系统呈现周期轨迹。如图2 (a) 所示, (图中正方形边长1,tanφ0=2/5, 初始出射点位于下底边且s0=0.1) , 可以看到这是一条周期轨迹, 整个轨迹由14段独立的直线段组成 (所谓“独立”这儿是指互不重叠) 。由于轨迹中点图是由轨迹直线段的中点组成, 所以它必然是由若干个离散的点构成 (点的个数通常等于独立直线段的数目, 某些情况下也可能小于线段数目, 因为不排除各个中点相互重叠的可能性) 。图2 (b) 给出了对应的轨迹中点图, 可以看出, 它由11个离散的点组成。这里的结论非常容易理解:离散点组成的轨迹中点图意味着弹球系统的周期行为。 一般平行四边形比矩形的情况稍显复杂一些, 通过观察发现, 平行四边形中的射线在传播途中其方位角只能取有限个不同的值, 也就是说整个轨迹由有限组平行直线段构成, 各组内的平行线段互不重叠, 随着反射的进行, 轨迹逐渐充满整个平行四边形。图3 (a) 所示的平行四边形底边长1.5, 高1, 锐内角72°, 射线初始方位角φ0=58°, 初始出射点位于下底边并