基于Foliation条件的离散动力系统二维流形计算.pdfVIP

第 31卷 第 4期 计 算 物 理 CHINESEJOURNALOF COMPUTATIONAL PHYSICS Vo1．31．NO．4 2014年7月 Ju1．，2014 文章编号：1001-246X(2014)04-0495．10 基于Foliation条件的离散动力系统二维流形计算 贾 蒙 (新乡学院机电工程学院，河南 新乡 453003) 摘 要：研究离散动力系统双 曲不动点的二维流形计算，利用不变流形轨道上 Jacobian矩阵能够传递导数这一特 殊性质 ，提出一种新的一维流形计算方法，通过预测 一校正两个步骤迅速确定流形上新 网格点，避免重复计算，并 简化精度控制条件．在此基础上，将基于流形面Foliation条件进行推广 ，推广后的Foliation条件能够控制二维流形 上的一维子流形的增长速度，从而实现二维流形在各个方向上的均匀增长．此外，算法可以同时用于二维稳定和不 稳定流形的计算．以超混沌三维 H6non映射和具有蝶形吸引子的Lorenz系统为例验证了算法的有效性． 关键词：离散动力系统；稳定流形；不稳定流形；导数传递；三维H6non映射；Lorenz系统；混沌吸引子 中图分类号：059 文献标志码 ：A 0 引言 稳定流形和不稳定流形在分析系统的动力学特性中起着非常重要的作用，它们充当着不同吸引子吸引 域的边界 ，将全空间划分为多个具有不同动力特性的不变子空间；而 当稳定与不稳定流形相交时，就会 引起 同宿、异宿以及混沌等复杂动力学行为的出现． 流形计算作为动力系统分析的一种几何方法，一直以来都受到研究人员的重视．研究的对象主要有两类 系统 ：连续动力系统(向量场)和离散动力系统．连续动力系统一般以微分方程 的形式来表示 ，其一维流形只 需通过积分求解初值问题就能得到，比较简单，二维流形计算就变得困难多了：其 中一个最主要的问题是当 系统的两个稳定(或不稳定)特征值的比值较大时，流形轨道在各个方 向上的增长速度不均 ，而且轨道朝其 中一个特征值方向汇集，从而使得计算的结果变差．针对这个问题，国内外的学者给出了许多不同解决办法， 详见文献 [1—7]． 离散动力系统一般 以映射的形式出现，所 以有时也称作映射动力系统．离散动力系统 的流形计算比较困 难 ，与连续动力系统的流形计算相 比，它的一个重要特点是轨道上新网格点对轨道上 已有的网格点的依赖 性，也就是说，新网格点必须是通过已有网格点或其插值经过映射后才能得到；而连续动力系统的流形计算 是通过求解初值问题来计算的，所以只有给定一个出发点，就能计算之后的轨道，而无需考虑出发点之前轨 道的状态．这种依赖性的存在 ，使得离散动力系统的流形计算更为困难． You ，Sire6 ，Parker和 Chua…，Hobson 给 出了基域迭代 的改进算法．BerndKrauskopf和 Hinke Osinga¨ ”通过每步向流形 中加入一个离散点来增长流形 ，结合 曲率控制技术 ，使得流形的增长速度 由流 形的局部 曲率来决定，该算法能够有效展示流形的细节．与此思想类似 ，文献 [14]中提出的算法能够在逆映 射不能显式表达的情况下计算系统的一维稳定流形，而之前提到的大部分算法都是把系统的稳定流形当作 其逆映射的不稳定流形来计算 的．Dellnitz和 Hohmann提出的细分法 ¨也能有效地计算一维流形 ；Danny Fundinger¨的算法与此思想类似．以上算法在计算的快速性上都有待提高．一维流形计算的问题没有有效 地解决，所以二维流形的计算就显得更为困难 了，除了要面对计算 向量场二维流形时的困难 ，还要牵扯到插 值点原像的搜索问题 ，并且插值点之后的轨道计算对于插值点之前的轨道存在依赖性，所 以就不得不对插值 点之前原本不存在的轨道也进行插值，从而使计算变得尤为复杂．目前对于离散系统二维流形的算法很少， 文献[17]中提出了一种应用Foliation条件 的算法，该算法 比较好地解决了流形各 向增长速度不均的问题， 收稿B期 ：2013—07—22；修 回日期 ：2013—11—26