四次平面多项式系统的Lyapunov量的复算法.doc

四次平面多项式系统的Lyapunov量的复算法 陈莹1 李静1 张伟2 1北京工业大学 应用数理学院, 北京100022, leejing@ 2北京工业大学 机械与电子工程学院, 北京100022, sandyzhang0@ 摘要： (1) 的Lyapunov量的复算法。这一结果可用于判定系统在原点的极限环的个数，对平面多项式系统的多极限环分岔的研究具有重要的理论指导意义。 二、Lyapunov量及其复算法 在平面多项式复系统 (2) 中，如果存在形式级数 和实数列使得，且有 (3) 则称(3)中的是(2)在原点处的第个Lyapunov量，级数为复系统(2)的Lyapunov级数。虽然Lyapunov级数不唯一，但第一个非零的Lyapunov量是唯一的[6]。 性质1 如果所有的Lyapunov量都为零，则原点是中心；否则，是细焦点；如果是第一个非零的Lyapunov量，则原点被称为是阶细焦点。 性质2 如果，则原点是稳定的；如果，则原点是不稳定的[7] 。 为了计算(2)在原点的Lyapunov量，只需找到(3)式中关于(2)的一个Lyapunov级数，等价地只需找出的系数向量。令是维向量，是级数中 项的系数。 定义1　定义阶矩阵，， (4) 在(4)中 是(2)中中项的系数。 定义2　定义的矩阵， 　　 (5) 引理1 如果是(2)的一个Lyapunov级数，则存在一系列实数，使得的系数向量满足下述的递推公式： (6) ， (7a) , (7b) 这里为维标准单位向量，当是偶数时，它的第个分量为1；且满足.反过来，任给实数列，由递归关系(6)-(7b)定义了向量序列，则以为系数向量的形式级数是(2)式的一个Lyapunov级数。 定理1 方程(2)在原点的第k个Lyapunov量是 (8) 其中由(6),(7a),(7b)式定义。 注 使用公式(8)计算Lyapunov量时，总可以取. 三、计算流程 计算一般平面多项式复系统的Lyapunov量的算法流程可按以下三个步骤进行： 第一步：参照在不同条件下不同方程的，按照定义1，依次计算及其相关矩阵. 第二步：按照定义2，计算相关矩阵，这时应注意到在程序中矩阵名称的设定会有所改变；结合第一步所得结果依引理1，计算. 第三步：结合第一、二步计算结果，由定理1所给计算公式，计算方程在原点的Lyapunov量 . 利用Maple计算Lyapunov量的算法流程图如下： 图1 ， (9) 这里 其中 均为实常数。 借助于Maple符号程序，获得复方程（9）的第一个的Lyapunov量： . 令，获得方程的第二个Lyapunov量： 令，获得方程的第三个Lyapunov量： . 令，获得方程的第四个Lyapunov量： 在中, . 令，获得方程的第五个Lyapunov量： 在中, 是一具有89项的表达式，限于篇幅这里忽略它的表出。 令，获得方程的第六个Lyapunov量： 令，获得方程的第七个Lyapunov量： . 令 获得方程的第八个Lyapunov量： . 令 则有 . 令取 ，是的最小的根，获得方程的第九个Lyapunov量： 在中, ， . 在以上的计算过程中，由于相关的计算复杂，本例计算到了. 三、结论 本文研究了一类缺二次项的四次平面多项式系统的Lyapunov量的复计算方法和Maple符号计算程序。给出了四次复系统的前九个Lyapunov量，这一结果说明通过扰动系数可以判定方程在原点至少存在九个小极限环。 参考文献 D. Wang. A Recursive Formula and Its Appli