方成鴻
(景德鎮(zhèn)陶瓷大學 信息工程學院,江西 景德鎮(zhèn) 333403)
關于平面二次多項式微分系統(tǒng)的定性分析,研究結果較豐富,文獻[1]有論述,對于三次系統(tǒng),多為分析不同形式的三次系統(tǒng),獲得系統(tǒng)的結構相圖或者極限環(huán)存在的條件,以便了解三次系統(tǒng)的各種全局結構.考慮如下形式的系統(tǒng):
(1)
文獻[2]包含c=0時的結論,以下設c=1,否則可作適當?shù)淖儞Q達成.本文分析點(a,b)位于第一、三象限的情形,這時系統(tǒng)(1)以(虛)橢圓1+ax2+by2=0為垂直等傾線,通過對奇點、無窮遠奇點的討論,獲得了系統(tǒng)(1)的全局結構相圖.
方程(1)的右端所定義的平面向量場關于坐標x軸對稱,故中心型奇點都是中心.
引理1系統(tǒng)(1)的奇點O(0, 0)是中心.
為應用方便,對文獻[4]的引理1進行推廣.設原點是
(2)
的孤立奇點,其中X(x,y)、Y(x,y)在原點附近具有任意階偏導數(shù),在原點處的函數(shù)值及一階偏導數(shù)值均為零.
定理1如圖1,設θ=θ0是特殊方向,曲線l1:y=f1(x)是系統(tǒng)(2)的水平等傾線,角域I內無其他等傾線,如果存在正數(shù)c使得當0 上述角域也可位于其他象限,相應修改不等號方向即可.同理可證. 定理2如圖1,設θ=θ0是特殊方向,曲線l2:y=f2(x)是系統(tǒng)(2)的垂直等傾線,角域Ⅱ內無其他等傾線,如果存在正數(shù)c使得當0 作Poincare變換u=y/x,z=1/x及時間變換dτ=dt/z2,系統(tǒng)(1)化為 (3) 系統(tǒng)(3)的奇點附近的軌線拓撲結構如圖2所示. 作變換v=x/y,z=1/y以及時間變換dτ=dt/z2,系統(tǒng)(1)化為 (4) 定理3當a=0時,系統(tǒng)(1)的全局結構相圖如圖4所示;當b=0時,系統(tǒng)(1)的全局結構相圖如圖5所示. 引理5當a=-1且-1/2≤b<0時,鞍點B2位于左半平面的不穩(wěn)定流形以A為ω極限點. 引理6當-1 定理4當a<-1且b<0或者a=-1且b<-a2/2時,系統(tǒng)(1)的全局結構相圖如圖6(1)所示;當a=-1且-a2/2≤b<0時,系統(tǒng)(1)的相圖如圖6(2)所示;當-1 引理7當-1 引理8當-1 引理10當-1 由于div(P,Q)=-2axy,可知a≠0時在任何一個象限,系統(tǒng)(1)都不存在閉軌線.根據(jù)引理9、10、11,結合無窮遠奇點的性態(tài)可得. 定理5設-1 關于ab<0的情形,另文分析,這時系統(tǒng)(1)的一條垂直等傾線是雙曲線1+ax2+by2=0.3 無窮遠奇點
4 全局結構