分数阶非线性系统动力学分析中数值算法可靠性与其诱导的复杂现象.pdfVIP

分数阶非线性系统动力学分析中数值算法可靠性 及其诱导的复杂现象1) *，**，2) * * 刘 杰 董鹏真 尚 钢 * （武汉科技学院非线性研究中心 武汉 430073 ） ** （武汉科技学院理学院 武汉 430073 ） 摘 要 分数阶系统中的混沌识别、动力学分析、控制与同步等问题最近已成为非线性科学研究领域的重要课题。目前， 研究者常用线性时不变的传递函数（LIT ）来近似分数阶的传递函数 1 ,(α∈(0,1)) 的方法进行相关数值模拟，并由 α s 此进行电路仿真设计。本文对于分数阶混沌系统的动力学分析近似方法的可靠性开展了初步分析讨论。首先，通过对 新近提出的 Liu 混沌系统的个案分析，研究发现：被众多文献广为采用的 Charef 和 Oustloup 几何近似方法所做出的线 性时不变的传递函数极可能导致具体数值分析中的“伪混沌”或者抑制混沌的“伪周期”情形发生，即“数值算法诱 导的复杂现象”。这一结果表明：在一些特定情形下，该类传统方法会丧失对现实的分数阶系统内在动力学进行准确刻 画的可靠性。进一步，本文指出，由于采用改进型 Oustloup 近似方法可在较大程度上改善传统的这两类几何近似方法 在拟合区间端点附近对分数阶微分算子进行近似的精确度，因而改进的 Oustloup 线性时不变的传递函数能够确保更高 程度的近似，从而在一定程度上避免“数值算法诱导的复杂现象”的发生。改进的 Oustloup 近似方法方法更适合工程 技术人员在混沌控制、同步研究等应用领域采用。当然，为保证可靠性，上述频域近似方法可以结合对应的时域上更 可靠的预估-校正方法使用。通过具体算例分析和仿真发现：改进型高精度线性时不变传递函数近似方法确实可以在一 定程度上更精确地接近原系统。因此，本文认为，在远离“由有序到混沌”的边缘处，采用改进型线性时不变传递函 数近似方法较之传统近似方法更具可行性，准确会得到提高；而在采用ADM预估-校正算法进行数值分析时，为避免产 生“伪混沌”现象，建议尽量设定足够小的迭代步长。最后，相关数值算例说明了本文分析结果的正确性。 关键词 改进的 Oustaloup 近似方法，ADM 预估-校正算法，分数阶混沌系统，稳定性分析，伪混沌 1、引言 最近十几年来，混沌现象的理论研究与实际应用一直是非线性科学研究中的热点问题[1]-[2]。经典 的微分方程理论表明，自治系统产生混沌吸引子的最小阶次数至少为 3 。而在引入了分数阶导数这一概 念后，研究者发现，自治系统产生混沌的最小阶次数可以小于 3，甚至更小。研究还表明：利用分数阶 动力系统能够较相应整数阶动力系统更准确地描述一些物理过程的特性，同时，这些结论在研究者提出 的典型的电学、热学实例中得到充分证实，如分数阶 Chua 电路可以用分数阶系统准确表示，且在阶次 为 2.7 就可能产生混沌现象，这些有趣的现象吸引了不同领域众多研究者的广泛兴趣[3]-[20]。 目前，寻找分析分数阶动力系统动力学行为的可靠数值方法备受研究者关注。较常用的数值分析方 法有基于波特图的“频域几何近似方法”和时域的广义 Adams-Bashforth-Moulton“预估-校正”算法[3]-[4]。 前者是“频域”近似方法，其基本思想是采用线性时不变传递函数近似分数阶传递函数，将分数阶动力 系统转化为更高维数的整数阶系统，再利用一系列分析整数阶动力系统的传统方法，如斯梅尔马蹄理论、 梅林科夫混沌解析分析等方法，或者计算李雅普诺夫指数、分析庞加莱截面、功率谱、计算信号的互信 息量、熵等分析方法相应的转化整数阶系统，进而反映原分数阶系统的动力学行为。其中，常用的分数 阶传递函数的近似方法为 Charef 方法或者 Oustaloup 方法[5]-[9]。目前，许多文献进行分数阶动力系统 数值仿真时，采用的就是上述文章中给出的近似表达式。后者则是一类基于算法稳定性、收敛性分析的 数值迭代算法，能够给出准确的、依赖于仿真步长和阶次参数的误差范围的估计，本质上是“时域”上