一類非自治廣義Birkhoff系統(tǒng)的穩(wěn)定性和分岔

曹秋鵬　陳向煒

1. 蘇州科技大學(xué)數(shù)理學(xué)院, 蘇州 215009; 2. 商丘師范學(xué)院物理與電氣信息學(xué)院, 商丘 476000;? 通信作者, E-mail: hnchenxw@163.com

?

一類非自治廣義Birkhoff系統(tǒng)的穩(wěn)定性和分岔

曹秋鵬1陳向煒2，?

1. 蘇州科技大學(xué)數(shù)理學(xué)院, 蘇州 215009; 2. 商丘師范學(xué)院物理與電氣信息學(xué)院, 商丘 476000;? 通信作者, E-mail: hnchenxw@163.com

研究一類非自治廣義Birkhoff系統(tǒng)的分岔。將該系統(tǒng)轉(zhuǎn)化為梯度系統(tǒng), 利用梯度系統(tǒng)的性質(zhì)研究這一類系統(tǒng)平衡點的穩(wěn)定性。研究表明, 當系統(tǒng)含有某個參數(shù)時, 系統(tǒng)平衡點的數(shù)目和穩(wěn)定性將會隨參數(shù)的變化而發(fā)生改變, 從而產(chǎn)生分岔現(xiàn)象。

廣義Birkhoff系統(tǒng); 梯度系統(tǒng); 分岔

北京大學(xué)學(xué)報（自然科學(xué)版）第52卷第4期2016年7月

Acta Scientiarum Naturalium Universitatis Pekinensis, Vol. 52, No. 4 (July 2016)

1927 年, Birkhoff[1]在《Dynamical systems》(動力系統(tǒng))一書中給出一類新型的積分變分原理和運動微分方程, Santilli[2]稱之為Pfaff-Birkhoff原理和Birkhoff方程。Birkhoff系統(tǒng)動力學(xué)是Hamilton力學(xué)的自然推廣。近年來, 對Birkhoff系統(tǒng)動力學(xué)的研究非?；钴S, 取得一些重要進展, 主要集中在Birkhoff系統(tǒng)的積分理論[3]、逆問題[4]、穩(wěn)定性[5]和對稱性[6]等方面。1993 年, 梅鳳翔[7]研究了 Birkhoff 方程增加一個附加項的情形, 稱為廣義 Birkhoff方程。廣義Birkhoff系統(tǒng)動力學(xué)也取得豐富的研究成果, 同樣集中在廣義Birkhoff系統(tǒng)的逆問題[8]、積分理論[9]、對稱性[10]、穩(wěn)定性[11]等方面,但很少涉及系統(tǒng)分岔的問題。開展廣義Birkhoff系統(tǒng)的分岔研究, 可將非線性動力學(xué)相關(guān)理論推廣應(yīng)用到廣義Birkhoff系統(tǒng), 探討該系統(tǒng)的動力學(xué)行為,進一步完善 Birkhoff 系統(tǒng)的理論體系。因此, 對廣義Birkhoff系統(tǒng)分岔問題的研究有重要意義。

非線性動力學(xué)的分岔是動力系統(tǒng)的重要性質(zhì),是流體力學(xué)、電力學(xué)、非線性振動理論、控制理論、生態(tài)學(xué)等領(lǐng)域研究的重點內(nèi)容[12-14]。2000年,陳向煒等[15-16]首次研究Birkhoff系統(tǒng)的分岔問題,分析了二階自治Birkhoff系統(tǒng)的極限點分岔、跨臨界分岔以及叉形分岔。梅鳳翔[17]研究了二階自治廣義Birkhoff系統(tǒng)平衡點分岔的相關(guān)問題, 也討論了系統(tǒng)的極限點分岔、跨臨界分岔以及叉形分岔。但是, 這些研究僅針對自治 Birkhoff 系統(tǒng), 未涉及非自治情形。

梯度系統(tǒng)是一類重要的動力學(xué)系統(tǒng)[18], 在研究運動穩(wěn)定性問題時具有很好的應(yīng)用價值。一旦動力學(xué)系統(tǒng)能夠轉(zhuǎn)化成梯度系統(tǒng), 便可以借助梯度系統(tǒng)的性質(zhì)研究動力學(xué)系統(tǒng)的穩(wěn)定性。文獻[19-23]研究了各類力學(xué)系統(tǒng)的梯度表示, 利用梯度系統(tǒng)的性質(zhì)分析了這些系統(tǒng)的平衡穩(wěn)定性。Mei等[24]利用梯度系統(tǒng)的性質(zhì), 研究了自治廣義 Birkhoff 系統(tǒng)的分岔。

在上述研究的基礎(chǔ)上, 本文利用梯度系統(tǒng)的性質(zhì), 進一步研究含參數(shù)非自治廣義Birkhoff系統(tǒng)的穩(wěn)定性和分岔。研究表明, 當系統(tǒng)所含的某個參數(shù)發(fā)生變化時, 系統(tǒng)平衡點的數(shù)目和穩(wěn)定性將會隨參數(shù)的變化而發(fā)生改變, 從而產(chǎn)生分岔現(xiàn)象。

1　廣義Birkhoff系統(tǒng)

廣義Birkhoff 系統(tǒng)的運動微分方程[17]為

可以寫成如下形式:

其中,

為 Birkhoff 協(xié)變張量,μνΩ為 Birkhoff 逆變張量,它們之間有關(guān)系:

2　梯度系統(tǒng)

梯度系統(tǒng)的微分方程有以下形式[17]:

其中, V=V(x)稱為勢函數(shù)。

梯度系統(tǒng)有以下兩個重要性質(zhì)。

性質(zhì) 1對于系統(tǒng)(4)所有的x, 都有0V≤˙, 當且僅當x為系統(tǒng)的平衡點時, 0V=˙。

性質(zhì) 2對于系統(tǒng)(4)的線性化系統(tǒng), 在任意平衡點處其特征方程只有實根。

由Lyapunov 一次近似理論可得如下命題。

命題 1如果梯度系統(tǒng)(4)的一次近似特征方程的根皆為負數(shù), 則平衡位置是漸近穩(wěn)定的; 如果有正根, 則平衡位置是不穩(wěn)定的; 如果存在單根 0 且無正根, 則平衡位置是穩(wěn)定的, 但不是漸近穩(wěn)定的;如果存在重根0, 則平衡位置是不穩(wěn)定的。

3　系統(tǒng)的梯度表示

一般情況下, 非自治廣義Birkhoff系統(tǒng)(1)并不是梯度系統(tǒng)。如果系統(tǒng)(1)滿足

同時有

此時, 方程(1)可以找到勢函數(shù)V=V(a), 使得

這樣, 非自治廣義 Birkhoff 系統(tǒng)(1)便成為一個梯度系統(tǒng)。于是, 我們利用梯度系統(tǒng)的性質(zhì)研究該系統(tǒng)的分岔。

4　系統(tǒng)平衡點的靜態(tài)分岔

假設(shè)非自治廣義Birkhoff系統(tǒng)(1)的Birkhoff函數(shù) B, Birkhoff函數(shù)組Rv或附加項νΛ含有某一個常參數(shù)μ, 則可以將非自治廣義Birkhoff系統(tǒng)(1)寫成如下形式:

其中,

且等式右端滿足式(6)。

對于研究內(nèi)容 1, 若非自治廣義Birkhoff系統(tǒng)(8)能夠成為一個梯度系統(tǒng), 那么平衡點的穩(wěn)定性由命題1判定, 隨著參數(shù)μ的變化, 梯度系統(tǒng)的一次近似特征方程的根的正負性可能發(fā)生變化。

對于研究內(nèi)容 2, 當 F(a, μ)=0 的解是一些曲線時, 分岔點是某兩條解曲線的交點。下面給出系統(tǒng)(8)發(fā)生研究內(nèi)容 2 情況的一個必要條件。

定理1對于系統(tǒng)(8), 設(shè)點(a0, μ0)處有 F(a0, μ0)=0, 在(a0, μ0)的鄰域內(nèi) F(a, μ)對 a 可微, 且F(a, μ)和DF(a, μ)對a, μ 連續(xù)。假如(a0, μ0)是F(a, μ)=0的一個分岔點, 則|DF(a0, μ0)|=0。其中DF(a, μ)表示F(a, μ)關(guān)于a的Jacobian矩陣。

證明: 假設(shè)|DF(a0, μ0)|≠0, 那么由隱函數(shù)定理,可以得到當|μ-μ0|<<1時, F(a, μ)=0唯一地確定了一條解曲線 a =a(μ), 使得a0=a(μ0)。此時, (a, μ0) 不能成為分岔點與(a0, μ0) 是分岔點矛盾。定理 1 得證。

下面利用兩個算例, 分別從系統(tǒng)(8)的平衡點穩(wěn)定性的變化和平衡點個數(shù)的變化, 說明系統(tǒng)(8)的靜態(tài)分岔。

5　算例

例1非自治廣義Birkhoff系統(tǒng):

其中μ是參數(shù)。

下面, 試寫出該系統(tǒng)的運動微分方程, 并分析μ對系統(tǒng)平衡位置穩(wěn)定性的影響。

由方程(8), 得到系統(tǒng)的微分方程:

顯然這是一個梯度系統(tǒng)。此時系統(tǒng)的一次近似特征方程為

當1μ＞時, 方程(11)有兩個負實根, 此時平衡位置是漸近穩(wěn)定的; 當1μ=時, 方程(11)有一個根是 0,一個根是-2, 此時平衡位置是穩(wěn)定的; 當1μ＜時,方程(11)兩個根一正一負, 此時平衡位置是不穩(wěn)定的。因此, 當1μ=時, 系統(tǒng)發(fā)生分岔, 1μ=是系統(tǒng)的分岔值。

例2非自治廣義Birkhoff系統(tǒng):

其中μ是參數(shù)。

下面, 試寫出該系統(tǒng)的運動微分方程, 并分析μ 對系統(tǒng)平衡位置穩(wěn)定性的影響。

由方程(8)得到系統(tǒng)的微分方程:

顯然, 方程(12)滿足式(5)和(6), 是一個梯度系統(tǒng)。系統(tǒng)(12)的平衡點個數(shù)與參數(shù)μ的取值有關(guān)。

對于任意的μ, (0, μ)總是方程(12)的解, 且

由定理 1 可知, 點(,0)0可能成為系統(tǒng)(12)的分岔點。下面具體分析分岔情況。

當0μ=時, 方程(12)有一個平衡點(0, 0), 方程(12)在(0, 0)點處的一次近似特征方程為

此方程有根0λ=和1λ=-, 由命題1可知此時平衡點是穩(wěn)定的, 但不是漸近穩(wěn)定的。

1) 平衡點（）0,0處的一次近似特征方程為

此方程有根λ=-μ2和λ=-1, 由命題1可知此平衡點是漸近穩(wěn)定的。

2) 我們對平衡點(μ, 0)做如下處理: 令

則方程(12)變成

同樣, 此平衡點是漸近穩(wěn)定的。

那么方程(12)變成

處的一次近似特征方程為

由命題1可知此平衡點是不穩(wěn)定的。

我們發(fā)現(xiàn), 含有參數(shù)μ的系統(tǒng)(12), 當0μ=時,平衡點只有1個, 當0μ≠時, 平衡點有3個, 平衡點的個數(shù)發(fā)生了改變, 并且對于同一平衡點(0, 0),兩種情況下其穩(wěn)定性并不相同。因此, (0, 0)成為系統(tǒng)(12)的分岔點, 0μ=是該系統(tǒng)的分岔值。

6　結(jié)論

本文把利用梯度系統(tǒng)研究穩(wěn)定性的方法推廣應(yīng)用到一類非自治廣義 Birkhoff 系統(tǒng), 該系統(tǒng)在滿足條件(5)和(6)情況下就可轉(zhuǎn)化為一個梯度系統(tǒng), 于是可以利用梯度系統(tǒng)的性質(zhì)研究這類非自治廣義Birkhoff 系統(tǒng)的穩(wěn)定性和分岔。本文的例子說明,隨著參數(shù)的變化, 系統(tǒng)平衡位置的穩(wěn)定性以及平衡點的個數(shù)會隨之變化, 系統(tǒng)發(fā)生分岔。當然, 這里的分岔指系統(tǒng)平衡點處的靜態(tài)分岔。

[1] Birkhoff G D. Dynamical systems. Providence: AMS College Publisher, 1927

[2] Santilli R M. Foundations of theoretical mechanicsⅠ. New York: Springer, 1978

[3] Zhang Hongbin, Chen Liqun, Gu Shulong, et al. The discrete variational principle and the first integrals of Birkhoff systems. Chinese Physics B, 2007, 16(3): 582-587

[4] 張永發(fā), 梅鳳翔. Birkhoff系統(tǒng)動力學(xué)逆問題的兩種提法和解法. 北京理工大學(xué)學(xué)報, 1996, 16(4): 352-356

[5] 傅景禮, 陳立群, 薛紜, 等. 相對論 Birkhoff 系統(tǒng)的平衡穩(wěn)定性. 物理學(xué)報, 2002, 51(12): 2683-2689

[6] 張毅. 相對論性力學(xué)系統(tǒng)的 Birkhoff 對稱性與守恒量. 物理學(xué)報, 2012, 61(21): 214501

[7] 梅鳳翔. Birkhoff系統(tǒng)的Noether理論. 中國科學(xué): A輯, 1993, 23(7): 709-717

[8] 梅鳳翔, 解加芳, 冮鐵強. 廣義 Birkhoff 系統(tǒng)動力學(xué)的一類逆問題. 物理學(xué)報, 2008, 57(8): 4649-4651

[9] 葛偉寬, 梅鳳翔. 廣義Birkhoff系統(tǒng)的時間積分定理. 物理學(xué)報, 2009, 58(2): 699-702

[10] Li Yanmin. Lie symmetries, perturbation to symmetries and adiabatic invariants of a generalized Birkhoff system. Chinese Physics Letters, 2010, 27(1): 010202

[11] 張毅. 自治廣義 Birkhoff 系統(tǒng)的平衡穩(wěn)定性. 物理學(xué)報, 2010, 59(1): 20-24

[12] Wen Dengzhe, Chen Yushu. Bifurcation analysis of fan casing under rotating air flow excitation. Applied Mathematics and Mechanics, 2014, 35(9): 1099-1114 [13] Ma Meiling, Min Fuhong. Bifurcation behavior and coexisting motions in a time-delayed power system. Chinese Physics B, 2015, 24(3): 030501

[14] Stiefs D, Gross T, Steuer R, et al. Computation and visualization of bifurcation surfaces. International Journal of Bifurcation and Chaos, 2008, 18(8): 2191-2206

[15] Chen Xiangwei, Mei Fengxiang. Existence of periodic solutions for higher order autonomous Birkhoff systems. Journal of Beijing Institute of Technology, 2000, 9(2): 125-130

[16] 陳向煒, 羅紹凱, 梅鳳翔. 二階自治 Birkhoff 系統(tǒng)的平衡點分岔. 固體力學(xué)學(xué)報, 2000, 21(3): 251-255

[17] 梅鳳翔. 廣義 Birkhoff 系統(tǒng)動力學(xué). 北京: 北京理工大學(xué)出版社, 2013

[18] Hirsch M W, Smale S, Devaney R L. Differential equations, dynamical systems and an introduction to chaos. Singapore: Elsevier, 2008

[19] 梅鳳翔. 關(guān)于梯度系統(tǒng). 力學(xué)與實踐, 2012, 34(1): 89-90

[20] 梅鳳翔, 吳惠彬. 一階 Lagrange 系統(tǒng)的梯度表示.物理學(xué)報, 2013, 62(21): 214501

[21] 陳向煒, 李彥敏, 梅鳳翔. 雙參數(shù)對廣義 Hamilton系統(tǒng)穩(wěn)定性的影響. 應(yīng)用數(shù)學(xué)和力學(xué), 2014, 35(12): 1392-1397

[22] 曹秋鵬, 張毅, 陳向煒. 約束自治廣義 Birkhoff 系統(tǒng)平衡穩(wěn)定性的梯度系統(tǒng)方法. 云南大學(xué)學(xué)報: 自然科學(xué)版, 2015, 37(2): 228-232

[23] 梅鳳翔, 吳惠彬. 廣義 Birkhoff 系統(tǒng)的梯度表示.動力學(xué)與控制學(xué)報, 2012, 10(4): 289-292

[24] Mei Fengxiang, Wu Huibin. Bifurcation for the generalized Birkhoffian system. Chinese Physics B, 2015, 24(5): 054501

Stability and Bifurcation for a Type of Non-autonomous Generalized Birkhoffian Systems

CAO Qiupeng1， CHEN Xiangwei2，?

1. School of Mathematics and Physics, Suzhou University of Science and Technology, Suzhou 215009; 2. Department of Physics and Information Engineering, Shangqiu Normal University, Shangqiu 476000; ? Corresponding author, E-mail: hnchenxw@163.com

Bifurcation for a type of non-autonomous generalized Birkhoffian systems is studied. Gradient representations for this type of non-autonomous generalized Birkhoffian systems are given. The stability of equilibrium point of these systems is discussed by the characteristic of the gradient system. Further the systems which contain some parameter are studied. The stability and the number of equilibrium point will change along with the change of the parameter to produce the bifurcation phenomenon.

generalized Birkhoffian system; gradient system; bifurcation

O316

10.13209/j.0479-8023.2016.067

國家自然科學(xué)基金(11372169)和蘇州科技大學(xué)研究生科研創(chuàng)新計劃(SKCX14_056)資助

2015-10-08;

2016-02-05; 網(wǎng)絡(luò)出版日期: 2016-07-14