计算周期解的数值方法.pdfVIP

第 16 卷 第 4 期 　 计 　　算 　　物 　　理 　　 Vol . 16 ,No . 4 1999 年 7 月 　 CH IN ESE J OU RNAL OF COM PU TA T IONAL PH YSICS 　　 　J ul . , 1999 　 计算周期解的数值方法 张锁春 ( 中国科学院应用数学研究所 ,北京 　100080) 摘 　要 　比较系统地叙述计算常微分方程 自治系统周期解的数值方法 ,重点介绍方法或算法的原 始构想和过程 。 关键词 　常微分方程 　自治系统 　周期解 　数值方法 中图分类号 　O24 1 0 　引言 周期振荡现象在 自然界普遍存在着 ,在物理 、化学 、生物 、电子学和工程技术等问题中经常 遇到 ,在数学 中求解微分方程 的周期解是一个古老而又 困难 的问题 。著名 的希尔伯特 ( Hilbert ) 第 16 问题就是有关微分方程周期解存在性的判定问题 ,至今进展甚少 ,可见难度之 大 。 由于计算机的普及和计算技术的发展 ,人们的研究思路和方法已发生改变 ,往往是先计算 后分析 。从计算所绘的图象中观察到有周期解的存在 ,然后再从理论上给予严格的证明。一 ( ) [ 1 ] 些来自实际问题 的研究 , 如生物化学 中的布鲁塞尔振子 Brusselator 、著名的 Belousov ( ) [2 ] Zhabotinskii 化学反应中产生的俄勒冈振子 Oregonator 等 ,人们的兴趣不在于周期解存在 性的理论证明 ,更多地是关心周期解的位置 、形状及周期的大小 。近十多年来发展迅速的非线 性现象的研究中 ,与耗散结构 、Hopf 分歧和混沌等有关的问题 ,常常需要数值地计算周期轨 线 。这就意味着在理论分析先天不足的条件下 ,还要在计算机上数值地计算出周期解 。于是 计算方法的研究极为重要 。 1 　难点 数学上求微分方程周期解和周期的问题可归结为求解下列的数学问题 : ( ) 寻找{ u t , T} , 使得 ( ) d u t ( ( ) ) = A u t d t ( ) 1 ( ) ( ) u t = u t + T ( ) n ( ) 成立 。其中 u t 是 R 中的实值向量函数, t 是时间变量, A ·是一个非线性算子, T 是最小