二维不稳定流形的计算-计算物理.PDFVIP

。．．。。。。 计门 算异 物仞 挫理 2005年翥碧臀期1 CHINESEJOURNALOFCOMPUTATIONALPHYSICS [文章编号】1001．246X(2005)06．0549—06 二维不稳定流形的计算 李清都1， 杨晓松2 (1．重庆邮电学院非线性系统研究所，重庆400065；2．华中科技大学数学系，湖北武汉430074) [摘要]提出了动力系统中稳定流形和不稳定流形的一种实用的快速算法，可以求得稳定流形和不稳定流形 的直观图像，从而从几何角度研究动力系统的动态行为和稳定性区域的边界特征．算法由两步构成：①在不稳 定流形上求得一些分布均匀的点，以精确反映流形的每个细节；②借助三角形剖分或二维单纯形剖分利用①的 算法将这些点画出直观流形图像． [关键词】 稳定流形；不稳定流形；Lorenz吸引子；动力系统 [中图分类号】0193 [文献标识码] A O 引言 稳定流形和不稳定流形的计算有着深刻的理论背景和广泛的实用意义．由于自然界中所有随时间而变 化的现象都可以用适当的动力系统模型描述(如微分方程、差分方程、随机过程等)，所以动力系统的一般理 论得到了广泛而深入的研究．在动力系统的研究中，把握稳定和不稳定流形的几何结构是动力系统理论的中 心问题之一．由于稳定和不稳定流形很难通过解析表达式求解，对其做近似计算就成为重要手段，这是当代 动力系统理论急需解决的主要问题之一．在动力系统理论中，不稳定流形往往预示着系统的未来走向．近期 人们在研究一类具有实际背景的系统，比如，基于天体动力学稳定和不稳定流形的航天器运行轨道的最优控 制设计…、电力系统中心1稳定区域的边界等，常可以通过计算某类轨道的稳定流形来实现．因此研究稳定和 不稳定流形对从整体上把握系统的性能有重要意义．如果能够直观地刻画它们，就会对上述研究更有利．然 而，目前的稳定和不稳定流形算法却无法快速有效的达到上述目的． 本文就三维常微分方程描述的自治系统的情形具体讨论稳定和不稳定流形的算法． 设微分动力系统 露=f(z)， (1) 其中髫∈R3，映射，：R3卜R3充分光滑，鞍点Xo处的稳定和不稳定流形定义为 W4(石o)={戈∈R3lira声。(z)=石o}， *—’+∞ IF(菇o)：{石∈R3lira声；(戈)：菇o}． 这里庐。代表方程(1)的流，即该方程的积分曲线或解．对于鞍形周期轨道工1的稳定和不稳定流形也可以用 同样的方法来定义．注意到稳定流形是不稳定流形的时间反向，通过求解不稳定流形就能得到稳定流形，因 此只讨论不稳定流形的情形就足够了．简单起见，我们主要介绍对于鞍点的不稳定流形算法，而对于周期轨 道的不稳定流形算法，只需稍做修改即可．稳定和不稳定流形的理论表明，在戈。的邻域存在的局部稳定和不 E“(戈。)相切．由于一维的不稳定流形仅由单个轨线构成，直接求解即可，比较简单，所以下面主要对二维的 稳定和不稳定流形进行详细讨论． [收稿日期】2004—08—04；【修回日期】2005—01—27 【基金项目】重庆邮电学院微电子学与固体电子学学科基金以及新世纪优秀人才支持计划资助项目 [作者简介】李清都(1980一)，男，湖北襄樊，硕士，主要从事混沌、电路与系统、神经网络等方面研究，重庆邮电学院非线性系统研究所， 400065． 550 计 算 物 理 第22卷 1过去的算法 目前，关于不稳定流形矿(石。)的计算方法比较多，但基本思想却很相似，大都是从雅可比矩阵of(菇。) 矿(省。)上向外求出满足一定条件的一个新圆，然后以新圆为初始圆，继续求解，逐步求出矿(菇。)． 根据数值求解方法的不同，对于目前连续系统的二维流形算法大致可分为两类． 人b1的做法是，在初始圆上选出一些均匀分布点来代替初始圆，对每个点求解初值问题，当达到某个事先约