连续时间系统同宿轨的搜索算法及其应用-物理学报.pdf

物理学报 Acta Phys. Sin. Vol. 62, No. 10 (2013) 100501 连续时间系统同宿轨的搜索算法及其应用* 12† 2 1 杨芳艳 胡明 姚尚平 1) ( 重庆邮电大学, 工业物联网与网络化控制教育部重点实验室, 重庆 400065 ) 2) ( 重庆邮电大学非线性电路与系统研究所, 重庆 400065 ) ( 2012年8月31 日收到; 2013年1月10 日收到修改稿) 同宿轨的求解是非线性系统领域的核心问题之一, 特别是对动力系统分岔与混沌的研究有重要意义. 根据同 宿轨的几何特点, 采用轨线逼近的方式, 通过定义逼近轨线与鞍点的距离, 将同宿轨的求解转化为求距离最小值的 无约束非线性优化问题. 为了提高优化结果的完整性, 还提出了基于区间细分的搜索算法和实现方法, 并找出了 Lorenz 系统, Shimizu-Morioka 系统和超混沌Lorenz 系统等的多个同宿轨道和对应参数, 验证了本文方法的有效性. 关键词: 混沌, 同宿轨, 非线性系统, 数值计算 PACS: 05.10.−a, 05.45.−a, 05.45.Pq DOI: 10.7498/aps.62.100501 系统, 因难以获得合适的解析表达式而失效. 关于 1 引言 同宿轨的数值计算, 由于同宿轨可看作是过鞍点且 有无穷周期的周期轨, 因此人们最早采用周期轨进 同宿轨是从一个鞍点到其本身的轨线, 是研究 21 行逼近 . 但由于不适用于鞍点特征值实部较小 系统状态结构和稳定性的基础, 与非线性动力系统 22 的情形, 例如Lorenz 系统 , 从而引发人们开始 1 的分岔和混沌有着密切关系, 在Hamilton 力学 、 从不同角度进行尝试. 随着当今计算机科学和数值 2−4 56 非线性震荡 、混沌信号产生器设计 中有重 计算的发展, 国内外涌现出不少方法, 从核心思想 要应用. 同宿轨常常是混沌运动的组织中心, 同宿 ´ 上, 可分为以下几类: Rodrıguez-Luis 和Freire 的打 分岔和同宿相交又往往是动力系统产生混沌的源 靶法 2324 , Beyn, Friedman, Sandstede, Bao, Witte 7 25−29 泉. 借助于同宿轨, 人们可以利用Silnikov 定理 、 等的边值法 , Lenci, Tigan 等基于流形的方 89 10−12 30−33 Smale-Birkhoff 定理 、Melnikov 方法 、拓 法 . 这些算法各具特色, 但也有不足. 打靶法 扑马蹄理论13−15 来判定混沌行为, 分析混沌产生 简单直观, 但是容易遇到截面选取的困难. 因为轨 的内在机制. 因此, 同宿