在研究平面系统时,首先要求解出奇点,然后分析奇点类型及其稳定性态,故而奇点性态分析在平面二次系统中占有重要地位。本文运用三种不同的方法对一类二次系统(II)类方程的奇点性态进行了细致分析。 In order to study the planar system, it is first required to find out the singular points, then analyze the singularity type and its stability. Therefore, the singularity behavior analysis plays an im-portant role in the planar quadratic system. In this paper, three different methods are used to an-alyze the singularities of a class of quadratic system (II) equations.
李丽君1,2,林建青3
1山东师范大学,数学与统计学院,山东 济南
2临沂大学,数学与统计学院,山东 临沂
3朔州师范高等专科学校,数计系,山西 朔州
收稿日期:2018年6月19日;录用日期:2018年7月11日;发布日期:2018年7月18日
在研究平面系统时,首先要求解出奇点,然后分析奇点类型及其稳定性态,故而奇点性态分析在平面二次系统中占有重要地位。本文运用三种不同的方法对一类二次系统(II)类方程的奇点性态进行了细致分析。
关键词 :二次系统,(II)类方程,奇点,定性分析
Copyright © 2018 by authors and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
常微分方程是数学中的一门重要学科。对于平面系统,我们在分析它的结构时首先要求解出奇点,对于奇点的类型大致分为以下几类:结点、焦点、中心、鞍点、高阶奇点等等。众所周知,奇点的线性化系数矩阵的特征值在奇点的初等分类中起着关键的作用。在有零特征值、纯虚数特征值和全零特征值的情况下有许多判断奇点的类型及稳定性的方法,本文对一类二次系统(II)类方程的奇点性态进行了分析。关于平面二次系统,有下述叶彦谦分类:
(I)类方程: x ˙ = − y + d x + l x 2 + m x y + n y 2 , y ˙ = x 。
(II)类方程: x ˙ = − y + d x + l x 2 + m x y + n y 2 , y ˙ = x ( 1 + a x ) , a ≠ 0 。
(III)类方程: x ˙ = − y + d x + l x 2 + m x y + n y 2 , y ˙ = x ( 1 + a x + b y ) , b ≠ 0 。
对于上述三类方程,当 d ≠ 0 时,以 O ( 0 , 0 ) 为粗焦点, d < 0 时O为稳定, d > 0 时O为不稳定。 d = 0 时,O为细焦点。
考虑平面二次系统
x ˙ = P ( x , y ) , y ˙ = Q ( x , y ) (2.1)
首先给出关于周期函数积分的一个引理。
引理1 ( [
F ( θ ) = ∫ 0 θ f ( s ) d s = g θ + φ ( θ ) (2.2)
其中 φ ( θ ) 仍以l为周期, g = 1 l ∫ 0 l f ( θ ) d θ 。
引理2 ( [
1) F正定: F ( 0 , 0 ) = 0 ; F ( x , y ) > 0 ,当 ( x , y ) ≠ ( 0 , 0 ) 。
2) d F d t | ( 2.1 ) < 0 ( > 0 ) ,当 ( x , y ) ≠ ( 0 , 0 ) 。
则系统(2.1)的平衡点 ( 0 , 0 ) 渐近稳定(不稳定)。
引理3 ( [
W 1 = m ( l + n ) − a ( b + 2 l ) .
W 2 = m a ( 5 a − m ) [ ( l + n ) 2 ( n + b ) − a 2 ( b + 2 l + n ) ] .
W 3 = m a 2 [ 2 a 2 + n ( l + 2 n ) ] [ ( l + n ) 2 ( n + b ) − a 2 ( b + 2 l + n ) ] .
有如下结论:
1) 当 d = 0 , W 1 ≠ 0 ,则点 O ( 0 , 0 ) 为一阶细焦点, W 1 < 0 时点 O ( 0 , 0 ) 为稳定, W 1 > 0 时点 O ( 0 , 0 ) 为不稳定;
2) 当 d = W 1 = 0 , W 2 ≠ 0 ,则点 O ( 0 , 0 ) 为二阶细焦点, W 2 < 0 时点 O ( 0 , 0 ) 为稳定, W 2 > 0 时点 O ( 0 , 0 ) 为不稳定;
3) 当 d = W 1 = W 2 = 0 , W 3 ≠ 0 ,则点 O ( 0 , 0 ) 为三阶细焦点, W 3 < 0 时点 O ( 0 , 0 ) 为稳定, W 3 > 0 时点 O ( 0 , 0 ) 为不稳定;
4) 当 d = W 1 = W 2 = W 3 时,点 O ( 0 , 0 ) 为中心。
本文考虑了一类二次系统(II)类方程在奇点的性态。
当 d = l = 0 , m = − a , n = − 1 时,(II)类方程化为 x ˙ = − y − a x y − y 2 , y ˙ = x + a x 2 。
解:令 − y − a x y − y 2 = 0 , x + a x 2 = 0 。得到下列4个奇点,分别为
O ( 0 , 0 ) , A ( 0 , − 1 ) , B ( − 1 a , 0 ) , C ( − 1 a , − 2 ) .
1) 由于 O ( 0 , 0 ) 是对应线性系统的中心,对其非线性系统在原点
方法一(后继函数法):令
消去dt,利用泰勒公式展开,得到
对充分小的c,求
其中
解得
故而可知
方法二(形式级数判别法):假设原系统具有下列级数形式的解
令
则
令三次项为0,则
取极坐标,令
上述化为
消去
因为
故
所以对应的三次齐次函数为
四次项显然为0,则
从而有
所以
方法三(焦点量判别法):利用焦点量公式,由于
2) 接下来我们对
令
代入原系统,得到在
特征方程为
3) 同理,对于
代入原系统,得到在
特征方程为
是鞍点,此解是不稳定的。
4) 对于
代入原系统,得到在
由于特征根为0,故而
接下来判断其在原点的性态。
解:由
代入比较系数,易得:
即
所以对于
山东省自然科学基金(ZR2018MA016)和国家自然科学基金(11601212)资助。
李丽君,林建青. 一类二次系统(II)类方程的奇点性态分析Singularity Analysis of a Class of Quadratic System (II) Equations[J]. 应用数学进展, 2018, 07(07): 788-793. https://doi.org/10.12677/AAM.2018.77095