龍梓軒，張 毅

(1.蘇州科技學(xué)院數(shù)理學(xué)院，江蘇蘇州 215009;2.蘇州科技學(xué)院土木工程學(xué)院，江蘇蘇州 215009)

分?jǐn)?shù)階微積分的歷史可以追溯到1695年Leibniz和L'Hospital之間的討論［1］，但是直到最近的40年間分?jǐn)?shù)階微積分在科學(xué)和工程各領(lǐng)域才得到廣泛的應(yīng)用［2－3］，包括:經(jīng)典和量子力學(xué)、場(chǎng)論和最優(yōu)控制等。1996年，Riewe［4］首先開(kāi)展了分?jǐn)?shù)階變分問(wèn)題的研究，提出了將分?jǐn)?shù)階微積分應(yīng)用于非保守系統(tǒng)動(dòng)力學(xué)和耗散系統(tǒng)的建模，初步形成了分?jǐn)?shù)階Euler-Lagrange方程和分?jǐn)?shù)階Hamilton方程。之后，分?jǐn)?shù)階變分問(wèn)題得到了應(yīng)用數(shù)學(xué)、物理學(xué)、動(dòng)力學(xué)與控制等研究領(lǐng)域的眾多學(xué)者的高度關(guān)注，并取得了一系列重要成果［5－9］。

為了建立非保守系統(tǒng)動(dòng)力學(xué)模型，El-Nabulsi于2005年在分?jǐn)?shù)階微積分的框架下基于Riemann-Liouville分?jǐn)?shù)階積分定義提出了一種新的建模方法［8］，稱之為類分?jǐn)?shù)階變分方法。該方法的特點(diǎn)在于分?jǐn)?shù)階時(shí)間積分僅引進(jìn)一個(gè)實(shí)參數(shù)α，所得到的Euler-Lagrange方程形式簡(jiǎn)單且類似于經(jīng)典的方程。該Euler-Lagrange方程的新穎之處在于存在一個(gè)作用在系統(tǒng)上的廣義分?jǐn)?shù)階外力，尤其是方程中不出現(xiàn)分?jǐn)?shù)階導(dǎo)數(shù)，而僅僅依賴于分?jǐn)?shù)階積分的階α。最近，El-Nabulsi將這種動(dòng)力學(xué)建模思想進(jìn)一步加以推廣，提出了基于按周期律拓展的分?jǐn)?shù)階積分的類分?jǐn)?shù)階模型［9］。

對(duì)稱性和守恒量的概念在物理學(xué)和數(shù)學(xué)中扮演著重要角色。對(duì)稱性在數(shù)學(xué)上是用變換群來(lái)描述的，即系統(tǒng)在變換群作用下保持的某種不變性，它對(duì)系統(tǒng)的動(dòng)力學(xué)行為和定性性質(zhì)具有深刻影響。守恒量在變分學(xué)和最優(yōu)控制方面的一個(gè)典型應(yīng)用是約化系統(tǒng)的自由度，從而實(shí)現(xiàn)維數(shù)的降低和微分方程積分的簡(jiǎn)化。1918年德國(guó)女?dāng)?shù)學(xué)家Noether提出了一個(gè)定理，揭示了Hamilton作用量在群的無(wú)限小變換下的不變性與守恒量之間的潛在關(guān)系。關(guān)于動(dòng)力學(xué)系統(tǒng)的對(duì)稱性與守恒量的研究是數(shù)學(xué)、物理學(xué)、分析力學(xué)的一個(gè)重要發(fā)展方向［10－16］。最近，基于El-Nabulsi提出的類分?jǐn)?shù)階模型，我們開(kāi)展了變分問(wèn)題的Noether對(duì)稱性研究［17－19］。本文進(jìn)一步研究基于按正弦周期律拓展的分?jǐn)?shù)階積分的類分?jǐn)?shù)階模型下的Noether理論，建立了類分?jǐn)?shù)階Noether定理，給出了類分?jǐn)?shù)階Noether對(duì)稱變換和Noether準(zhǔn)對(duì)稱變換的定義和判據(jù)，建立了類分?jǐn)?shù)階Noether對(duì)稱性與守恒量之間的關(guān)系。本文結(jié)果具有普遍意義，不僅經(jīng)典的Noether定理是其特例，而且可以進(jìn)一步推廣到各類約束力學(xué)系統(tǒng)，例如非完整系統(tǒng)、Birkhoff系統(tǒng)等。

1 基于按正弦周期律拓展的分?jǐn)?shù)階積分的類分?jǐn)?shù)階變分問(wèn)題

定義1［9］設(shè)f(t)是在區(qū)間［a，b］上的連續(xù)函數(shù)。對(duì)于t∈［a，b］，按正弦周期律拓展的階為α(α＞0)的左右分?jǐn)?shù)階積分定義為

假設(shè)力學(xué)系統(tǒng)的位形由n個(gè)廣義坐標(biāo)qk(k=1，…，n)來(lái)確定，系統(tǒng)的Lagrange函數(shù)為L(zhǎng)=L( τ ，q，q˙)。按照El-Nabulsi提出的非保守力學(xué)系統(tǒng)建模的類分?jǐn)?shù)階變分方法［9］，類分?jǐn)?shù)階變分問(wèn)題定義如下:

求積分泛函

在固定邊界條件

上述變分問(wèn)題稱為基于按正弦周期律拓展的分?jǐn)?shù)階積分的類分?jǐn)?shù)階變分問(wèn)題，泛函 (3)又可稱為基于按正弦周期律拓展的分?jǐn)?shù)階積分的類分?jǐn)?shù)階Hamilton作用量，當(dāng)α=1時(shí)，上述問(wèn)題成為經(jīng)典的力學(xué)系統(tǒng)的變分問(wèn)題。

對(duì)上式第二部分分部積分，并利用邊界條件 (4)，我們有

將式 (6)代入式 (5)，有

由積分區(qū)間［a，b］的任意性，我們得到

式 (8)可稱為基于按正弦周期律拓展的分?jǐn)?shù)階積分的類分?jǐn)?shù)階d'Alembert-Lagrange原理。原理 (8)既適用于完整約束系統(tǒng)，也適用于非完整約束系統(tǒng)。

對(duì)于完整約束系統(tǒng)，δqk(k=1，…，n)是相互獨(dú)立的，因此由式 (8)得出

方程 (9)是完整約束系統(tǒng)的類分?jǐn)?shù)階Euler-Lagrange 方程［9］。

2 類分?jǐn)?shù)階Hamilton作用量的變分

引進(jìn)無(wú)限小群變換

其展開(kāi)式為

其中εσ(σ=1，…，r)為無(wú)限小參數(shù)，，為無(wú)限小群變換的生成元或生成函數(shù)。

類分?jǐn)?shù)階作用量 (3)在變換前后的差為

對(duì)于任意函數(shù)F，其非等時(shí)變分Δ與等時(shí)變分δ之間有關(guān)系［10］

于是有

由式 (15)，式 (13)可表為

由式 (11)和 (15)，上式可進(jìn)一步表為

公式 (13)和 (17)是類分?jǐn)?shù)階Hamilton作用量變分的基本公式。

3 類分?jǐn)?shù)階Noether對(duì)稱變換

下面我們來(lái)建立基于按正弦周期律拓展的分?jǐn)?shù)階積分的類分?jǐn)?shù)階Noether對(duì)稱變換的定義和判據(jù)。

定義2 如果類分?jǐn)?shù)階Hamilton作用量 (3)是無(wú)限小群變換 (10)的不變量，即對(duì)每一個(gè)無(wú)限小變換，始終成立

則稱無(wú)限小變換是類分?jǐn)?shù)階Noether對(duì)稱變換。

由定義2和變分公式 (13)，我們有如下判據(jù)。

判據(jù)1 對(duì)于無(wú)限小群變換 (10)，如果滿足條件

則變換是系統(tǒng)的類分?jǐn)?shù)階Noether對(duì)稱變換。

由定義2和變分公式 (17)，我們有如下判據(jù)。

判據(jù)2 對(duì)于無(wú)限小群變換 (11)，如果滿足條件

則變換是系統(tǒng)的類分?jǐn)?shù)階Noether對(duì)稱變換。

利用關(guān)系

并考慮到無(wú)限小參數(shù)εσ的獨(dú)立性，則式 (19)可表為r個(gè)方程

式 (22)亦可作為系統(tǒng)對(duì)稱變換的判據(jù)。當(dāng)取r=1時(shí)，式 (22)可稱為類分?jǐn)?shù)階Noether等式

其次，我們來(lái)建立類分?jǐn)?shù)階Noether準(zhǔn)對(duì)稱變換的定義和判據(jù)。

設(shè)L′是另外的Lagrange函數(shù)，如果變換 (10)精確到一階小量滿足條件

則稱類分?jǐn)?shù)階Hamilton作用量 (3)是無(wú)限小群變換 (10)下的準(zhǔn)不變量。由此確定的L′和L具有同樣的運(yùn)動(dòng)微分方程，則變換稱為類分?jǐn)?shù)階Noether準(zhǔn)對(duì)稱變換，有

將式 (25)代入式 (24)，得到

式 (26)的左端在變換 (10)下是一階小量，因此，其右端應(yīng)為同階小量，可用ΔG代替G，而

因此有

定義3 如果類分?jǐn)?shù)階Hamilton作用量 (3)是無(wú)限小群變換 (10)的準(zhǔn)不變量，即對(duì)每一個(gè)無(wú)限小變換，始終成立

其中G=G(τ，q，)，則稱無(wú)限小變換是類分?jǐn)?shù)階Noether準(zhǔn)對(duì)稱變換。

由定義3和公式 (13)，我們得到如下判據(jù)。

判據(jù)3 對(duì)于無(wú)限小群變換 (10)，如果滿足條件

則變換 (29)是類分?jǐn)?shù)階Noether準(zhǔn)對(duì)稱變換。

由定義3和公式 (19)，我們得到如下判據(jù)。

判據(jù)4 對(duì)于無(wú)限小群變換 (11)，如果滿足條件

則變換是類分?jǐn)?shù)階Noether準(zhǔn)對(duì)稱變換。其中

利用關(guān)系 (21)， (31)，并考慮到εσ的獨(dú)立性，則式 (29)可表為r個(gè)方程

式 (32)亦可作為系統(tǒng)準(zhǔn)對(duì)稱變換的判據(jù)，函數(shù)Gσ=Gσ(τ，q)稱為規(guī)范函數(shù)。當(dāng)取r=1時(shí)，式(32)給出類分?jǐn)?shù)階Noether等式

利用判據(jù)1和判據(jù)2可以判斷所論系統(tǒng)的類分?jǐn)?shù)階Noether對(duì)稱性，利用判據(jù)3和判據(jù)4可以判斷所論系統(tǒng)的類分?jǐn)?shù)階Noether準(zhǔn)對(duì)稱性。

4 完整系統(tǒng)的類分?jǐn)?shù)階Noether定理

首先，基于上述類分?jǐn)?shù)階動(dòng)力學(xué)模型，給出完整系統(tǒng)的守恒量定義。

定義4 函數(shù)Iτ，q，˙( )q稱為完整系統(tǒng) (9)的守恒量，當(dāng)且僅當(dāng)沿著完整系統(tǒng)的類分?jǐn)?shù)階Euler-Lagrange方程 (9)的解曲線恒成立

對(duì)于完整系統(tǒng)，如果能找到其類分?jǐn)?shù)階Noether對(duì)稱變換或準(zhǔn)對(duì)稱變換，便可找到與之相應(yīng)的守恒量。有如下定理:

定理1 對(duì)于完整系統(tǒng) (9)，如果無(wú)限小群變換 (11)是定義2意義下的類分?jǐn)?shù)階Noether對(duì)稱變換，則系統(tǒng)存在r個(gè)線性獨(dú)立的守恒量，形如

證明 因?yàn)闊o(wú)限小群變換是系統(tǒng)的類分?jǐn)?shù)階Noether對(duì)稱變換，由定義2，有，即

將方程 (9)代入上式，由積分區(qū)間的任意性和參數(shù)εσ的獨(dú)立性，得到

積分之，便得式 (35)。證畢。

定理2 對(duì)于完整系統(tǒng) (9)，如果無(wú)限小群變換 (11)是定義3意義下的類分?jǐn)?shù)階Noether準(zhǔn)對(duì)稱變換，則系統(tǒng)存在r個(gè)線性獨(dú)立的守恒量，形如

定理1和定理2可稱為完整系統(tǒng)的類分?jǐn)?shù)階Noether定理。由Noether定理可知，對(duì)于完整約束系統(tǒng)，如能找到系統(tǒng)的一個(gè)類分?jǐn)?shù)階Noether對(duì)稱變換或準(zhǔn)對(duì)稱變換，便有可能得到系統(tǒng)的一個(gè)守恒量。

例 平面Kepler問(wèn)題的Lagrange函數(shù)為

試研究其基于按正弦周期律拓展的分?jǐn)?shù)階積分的類分?jǐn)?shù)階Noether對(duì)稱性與守恒量。

首先，尋找類分?jǐn)?shù)階Noether準(zhǔn)對(duì)稱變換。類分?jǐn)?shù)階廣義Noether等式 (32)給出

方程 (40)有解

生成元 (41)對(duì)應(yīng)于系統(tǒng)的類分?jǐn)?shù)階Noether對(duì)稱變換，生成元 (42)對(duì)應(yīng)于系統(tǒng)的類分?jǐn)?shù)階Noether準(zhǔn)對(duì)稱變換。

由定理1，相應(yīng)于生成元 (41)，系統(tǒng)存在如下守恒量

式 (43)是我們得到的系統(tǒng)的類分?jǐn)?shù)階Noether對(duì)稱性 (35)導(dǎo)致的守恒量。

由定理2，相應(yīng)于生成元 (42)，守恒量式(35)給出

因此，對(duì)應(yīng)生成元 (42)的無(wú)限小變換是平庸的。

［1］OLDHAM K B，SPANIER J.The fractional calculus［M］.San Diego:Academic Press，1974.

［2］PODLUBNY I.Fractional differential equations［M］.San Diego:Academic Press，1999.

［3］KILBAS A A，SRIVASTAVA H M，TRUJILLO J J.Theory and applications of fractional differential equations［M］.Amsterdam:Elsevier B V，2006.

［4］RIEWE F.Nonconservative lagrangian and Hamiltonian mechanics［J］.Phys Rev E，1996，53:1890－1899.

［5］AGRAWAL O P.Formulation of Euler-Lagrange equations for fractional variational problems［J］.J Math A-nal Appl，2002，272(1):368 －379.

［6］JUMARIE G.Fractional hamilton-jacobi equation for the optimal control of nonrandom fractional dynamics with fractional cost functions［J］.J Appl Math ＆ Computing，2007，23(1/2):215－228.

［7］ATANACKOVIC＇T M，KONJIK S，PILIPOVIC＇S，et al.Variational problems with fractional derivatives:Invariance conditions and Noether's theorem ［J］.Nonlinear Anal，2009，71:1504 －1517.

［8］EL-NABULSI A R.A fractional approach to nonconservative Lagrangian dynamical systems［J］.Fizika A，2005，14(4):289－298.

［9］EI-NABULSI A R.A periodic functional approach to the calculus of variations and the problem of time-dependent damped harmonic oscillators ［J］.Appl Math Lett，2011，24:1647－1653.

［10］梅鳳翔.李群和李代數(shù)對(duì)約束力學(xué)系統(tǒng)的應(yīng)用［M］.北京:科學(xué)出版社，1999.

［11］梅鳳翔.約束力學(xué)系統(tǒng)的對(duì)稱性與守恒量［M］.北京:北京理工大學(xué)出版社，2004.

［12］LUO S K，CAI J L，JIA L Q.Noether symmetry can lead to non-Noether conserved quantity of holonomic nonconservative systems in general Lie transformations［J］.Commun Theor Phys，2005，43(2):193－196.

［13］張毅.非完整力學(xué)系統(tǒng)的Hamilton對(duì)稱性［J］.中國(guó)科學(xué):物理學(xué) 力學(xué) 天文學(xué)，2010，40:1130－1137.

［14］賈利群，鄭世旺，張耀宇.事件空間中非Chetaev型非完整系統(tǒng)的Mei對(duì)稱性與Mei守恒量［J］.物理學(xué)報(bào)，2007，56(10):5575－5579.

［15］WU H B，MEI F X.Symmetry of lagrangians of holonomic systems in terms of quasi-coordinates［J］.Chin Phys B，2009，18(8):3145－3149.

［16］龍梓軒，張毅.Birkhoff系統(tǒng)Lie對(duì)稱性逆問(wèn)題的兩種提法和解法［J］.華東師范大學(xué)學(xué)報(bào):自然科學(xué)版，2012(3):49－55.

［17］ZHANG Y，ZHOU Y.Symmetries and conserved quantities for fractional action-like Pfaffian variational problems［J］.Nonlinear Dynamics，2013，73(1/2):783－793.

［18］LONG Z X，ZHANG Y.Noether's theorem for fractional variational problem from El-Nabulsi extended exponentially fractional integral in phase space［C］∥Acta Mechanica，2013，DOI 10.1007/s00707 －013 －0956 －5.

［19］張毅.相空間中類分?jǐn)?shù)階變分問(wèn)題的Noether對(duì)稱性與守恒量［J］.中山大學(xué)學(xué)報(bào):自然科學(xué)版，2013，52(4):45－50.