平面二次系统极限环及其稳定与分岔的计算-中国力学学会.docVIP

平面二次系统极限环及其稳定与分岔的计算 黄赪彪 邬华东 （中山大学力学系，广州 510275 ） 摘要 将平面二次多项式微分系统极限环相图的x坐标假设为广义谐函数；用增量迭代法近似算出极限环的 y坐标、频率、周期、稳定性指标，以及关于参数分岔的表达式；给出绕奇点（0, 0）具有三个极限环的例 子。 关键词 平面二次多项式微分系统，极限环，频率，周期，稳定性，分岔 引言 一个多世纪前，D.Hilbert 提出的23个数学问题中，第16问题第二部分当时，微分系统可以写为 (1a) (1b) 其中，是实常量。关于系统（1）的极限环研究，一个余世纪以来，许多数学力学工作者不懈努力，或者从奇点类型的讨论、定性分析极限环的数目、分布及分岔，或者利用计算机进行数值搜索，取得了一系列重要成果。其中，文[3]、[4]各自独立地证明了。本文通过引进自变量的非线性变换，将极限环相图的坐标假设为广义谐函数，而坐标及极限环的频率、周期和稳定性判别指标，乃至极限环关于参数的分岔曲线，则由增量迭代法半解析半数值近似表示。文末给出绕奇点（0, 0）具有三个极限环的例子，计算结果与数值积分法吻合良好。 1 极限环的计算 改写式（1a）为 (2) 这里，（下同）。由式（2）得 (3) 将式（3）代入式（1b），得 (4) 其中， 而 国家自然科学基金资助项目。 引进自变量的非线性变换 满足 并设该系统（1）的极限环的x坐标为广义谐函数 (5) 则 (6) 其中，、b、、分别表征极限环的振幅、偏心距、相角和频率，。将式（5）、（6）代入式（4），得 (7a) 式（7a）两边对从0开始积分，得 (8a) 将式（6）代入式（1a），得 (9a) 其中 由的非奇异性，式（9a）中相关各量须满足一定条件。 分别令式（8a）中，可得确定, b的两个方程 , (10,11) 至此，系统（1）的极限环原则上可已由式（5）—（11）给出。不过难以解析表出。下面以迭代法近似求解之。为保证迭代法的收敛性，引进小参数，作用于式（7a）,（8a）,（9a），相应地，记为 (7b) (8b) (9b) 迭代格式如下 (12a) (12b) (12c) 式（12a）中， 由以下二式确定 (13) (14) 而则由 (15) (16) 二式给出。若数列收敛，则它们的极限即为所求 给小参数予增量，，记，，相应地可求出，直至为止。 具体地， 将分别展开为富氏级数，并设它们的前n项和可满足足够的精度，即 (17) (18) 将式（17）、（18）代入式（8a）、（9a），分别记为 , (19, 20) 其中，，，，，利用谐波平衡原理，可得 (21a)