黃婷

摘 要：本文研究了一類三維非線性動力系統(tǒng)中心流形上流動隱函數(shù)形式級數(shù)計算方法，并給出了其代數(shù)公式，此公式是線性的，避免了復(fù)雜的積分運算，運用Mathematica系統(tǒng)軟件，基于代數(shù)遞推公式計算了三維微分系統(tǒng)中心流形上流的隱函數(shù)。該新算法為研究此類非線性動力系統(tǒng)的動力學(xué)問題，特別是穩(wěn)定性、Hopf分支問題提供了很大的便利。

關(guān)鍵詞：中心流形；穩(wěn)定性；Hopf分支；形式級數(shù)；計算機代數(shù)

中圖分類號：O322，O345

本文研究了一類三維非線性動力系統(tǒng)中心流形上流動隱函數(shù)形式級數(shù)計算方法，并給出了其代數(shù)公式，此公式是線性的，避免了復(fù)雜的積分運算，運用Mathematica系統(tǒng)軟件，基于代數(shù)遞推公式計算了三維微分系統(tǒng)中心流形上流的隱函數(shù)。該新算法為研究此類非線性動力系統(tǒng)的動力學(xué)問題，特別是穩(wěn)定性、Hopf分支問題提供了很大的便利。

1 三維系統(tǒng)中心流形的形式級數(shù)

文[6]中給出了方法的實質(zhì)，將二維分支系統(tǒng)形式級數(shù)的思想方法進行推廣，討論下列三維實解析系統(tǒng)：

其中，x，y，u，t，Akjl，Bkjl∈R（k，j，l∈N）。

同文[6]類似先利用復(fù)變換將系統(tǒng)（1）化為復(fù)自治微分系統(tǒng)：

（2）

其中，z，w，T，akjl ，bkjl，dkjl∈C（k，j，l∈N）。

并稱系統(tǒng)（1）與（2）互為伴隨系統(tǒng)。

定理1[ 6 ] 對系統(tǒng)（2），可逐項唯一確定形式級數(shù)：

確定；

確定；

2 算例與數(shù)學(xué)軟件的實現(xiàn)

考慮一類實三維系統(tǒng)中心流形的隱函數(shù)形式級數(shù)的計算問題。

其中a1，a2，b1，b2，d1，d2均為實變量。

同樣可以經(jīng)過復(fù)變換將系統(tǒng)（7）化為復(fù)自治微分系統(tǒng)：

根據(jù)定理1，有：

推論1

2，3，…，則可逐項確定形式級數(shù)（3），使得（4）成立，且當k≠j或者k=j，l≠0時，ckjl可由下列遞推公式：

確定，對任一正整數(shù)m，？滋m可由下列遞推公式：

應(yīng)用系統(tǒng)的強大符號運算功能，把推論1的公式編成運算程序，容易算出原點的前10個焦點量為：

實際上，系統(tǒng)（7）中心流形上流的方程在原點的所有焦點量為零，對應(yīng)原點為中心或者稱在其原點鄰域可積。

利用推論1的遞推公式逐步計算出系數(shù)ckjl并通過逆變換：

可得系統(tǒng)（7）中心流形上的流的隱函數(shù)形式為：

其中，

下面考慮另一個三維系統(tǒng)：

同樣可以經(jīng)過復(fù)變換將系統(tǒng)（12）化為復(fù)自治微分系統(tǒng)：

3 總結(jié)

本文研究了一類三維非線性動力系統(tǒng)中心流形上流動隱函數(shù)形式級數(shù)計算方法，并給出了其代數(shù)公式，此公式是線性的，避免了復(fù)雜的積分運算，運用Mathematica系統(tǒng)軟件，基于代數(shù)遞推公式計算了三維微分系統(tǒng)中心流形上流的隱函數(shù)。該新算法為研究此類非線性動力系統(tǒng)的動力學(xué)問題，特別是穩(wěn)定性、Hopf分支問題提供了很大的便利。

參考文獻：

[1] 張琪昌等.分岔與混沌理論及應(yīng)用[M].天津大學(xué)出版社， 2005，1～63.

[2] 于海，陳予恕.高維非線性動力學(xué)系統(tǒng)降維方法[J].力學(xué)進展，2009，39（2）：154-164.

[3] Carr J. Application of center Manifold Theory [M].Applied Mathematical Sciences 35，Springer verlag， NewYork，1981.

[4] 韓茂安.動力系統(tǒng)的周期解與分支理論[M].北京：科學(xué)出版社，2002.

[5] Yu P.Simplest normal forms of Hopf and generalized Hopf bifurcation [J]. International Journal of Bifurcation and Chaos， 1999， 10（9）：1917-1939.

[6] 劉一戎，李繼彬.論復(fù)自治微分系統(tǒng)的奇點量[J].中國科學(xué)（A輯），1989，19（3）：245-255.

[7] Q.wang.Y.Liu， H.Chen. Hopf bifurcation for a class of three-dimensional nonlinear dynamic systems[J].Bull.Sci.Math， 2010，134：786-798.

[8] Qinlong Wang，Wentao Huang， Bai-Lian Li. Limit cycles and singular point quantities for a3D Lotka-Valterra system[J].Applied Mathematics and Computation，2011 doi：10.1016.