金世欣，張 毅

(1．南京理工大學(xué)理學(xué)院，江蘇 南京 210094；2．蘇州科技學(xué)院土木工程學(xué)院，江蘇 蘇州 215011)

?

基于Caputo分?jǐn)?shù)階導(dǎo)數(shù)的含時(shí)滯的非保守系統(tǒng)動(dòng)力學(xué)的Noether對(duì)稱性*

金世欣1，張 毅2

(1．南京理工大學(xué)理學(xué)院，江蘇 南京 210094；2．蘇州科技學(xué)院土木工程學(xué)院，江蘇 蘇州 215011)

提出并研究基于Caputo分?jǐn)?shù)階導(dǎo)數(shù)的含時(shí)滯的力學(xué)系統(tǒng)的Noether對(duì)稱性與守恒量。建立了含時(shí)滯的非保守系統(tǒng)的分?jǐn)?shù)階運(yùn)動(dòng)微分方程；根據(jù)系統(tǒng)的含時(shí)滯的分?jǐn)?shù)階Hamilton作用量在無(wú)限小群變換下的泛函不變性，給出了含時(shí)滯的分?jǐn)?shù)階Noether對(duì)稱變換，Noether準(zhǔn)對(duì)稱變換以及Noether廣義準(zhǔn)對(duì)稱變換的定義判據(jù)；研究了含時(shí)滯的分?jǐn)?shù)階Noether對(duì)稱性與守恒量之間的聯(lián)系，并舉例說(shuō)明結(jié)果的應(yīng)用。

非保守系統(tǒng)；時(shí)滯；Caputo分?jǐn)?shù)階導(dǎo)數(shù)；Noether對(duì)稱性；守恒量

近年來(lái)，Baleanu等[17-19]研究了分?jǐn)?shù)階模型下的含時(shí)滯的變分和最優(yōu)化控制問(wèn)題。2012年，F(xiàn)rederico和Torres[20]首次討論了含時(shí)滯的變分和最優(yōu)化控制問(wèn)題的Noether定理；隨后，張毅等[21-25]研究建立了含時(shí)滯的約束力學(xué)系統(tǒng)的Noether對(duì)稱性與守恒量理論。盡管含時(shí)滯的分?jǐn)?shù)階變分和最優(yōu)控制問(wèn)題的研究已經(jīng)取得了一些重要成果，但是研究狀態(tài)變量或控制變量具有時(shí)滯的分?jǐn)?shù)階變分和控制系統(tǒng)的對(duì)稱性與守恒量問(wèn)題還是一個(gè)開(kāi)放的課題，特別是在不同的分?jǐn)?shù)階模型下的對(duì)稱性與守恒量的問(wèn)題。本文將進(jìn)一步研究基于Caputo分?jǐn)?shù)階導(dǎo)數(shù)定義下的含時(shí)滯的非保守系統(tǒng)動(dòng)力學(xué)的Noether對(duì)稱性與守恒量。建立相應(yīng)力學(xué)系統(tǒng)的分?jǐn)?shù)階Noether對(duì)稱性的定義和判據(jù)，并導(dǎo)出含時(shí)滯的分?jǐn)?shù)階Noether定理。

1 分?jǐn)?shù)階導(dǎo)數(shù)的定義及其若干性質(zhì)

在這一部分我們簡(jiǎn)單地回憶一下將要用到的一些Riemann-Lionville以及Caputo分?jǐn)?shù)階導(dǎo)數(shù)的定義和性質(zhì)。詳細(xì)的討論和證明，可參考文獻(xiàn)[1-3]。

左Riemann-Lionville分?jǐn)?shù)階導(dǎo)數(shù)的定義如下

(1)

右Remiann-Lionville分?jǐn)?shù)階導(dǎo)數(shù)的定義如下

(2)

相應(yīng)的Caputo分?jǐn)?shù)階導(dǎo)數(shù)的定義如下：左Caputo分?jǐn)?shù)階導(dǎo)數(shù)

(3)

右Caputo分?jǐn)?shù)階導(dǎo)數(shù)

(4)

其中Γ(*)為Gamma函數(shù)，滿足k-1≤α

(5)

(6)

以及

(7)

(8)

(9)

2 含時(shí)滯的分?jǐn)?shù)階運(yùn)動(dòng)微分方程

假設(shè)力學(xué)系統(tǒng)由n個(gè)廣義坐標(biāo)qs(s=1,2,…,n)來(lái)確定。考慮非保守力學(xué)系統(tǒng)的Hamilton原理

(10)

其中，Lagrange函數(shù)以及非勢(shì)廣義力為

(11)

(12)

且時(shí)滯常量τ

(13)

(14)

其中Ωs(t)為[t1-τ,t1]上的已知分段光滑函數(shù)。則原理(10)可寫(xiě)為

(15)

將式(15)的第三項(xiàng)，第五項(xiàng)進(jìn)行變量替換t=θ+τ，并考慮初始條件(13)，得到

(16)

考慮分部積分公式(6)和(7)以及(8)和(9)，并考慮條件(13)和(14)，則有

(17)

且有

(18)

以及

(19)

當(dāng)滿足條件

(20)

將(17)，(18)以及(19)式代入(16)式，并考慮到積分區(qū)間的任意性以及δqs的獨(dú)立性，得到

(21)

滿足式(20)，則方程(21)可稱為Caputo導(dǎo)數(shù)下的含時(shí)滯的非保守系統(tǒng)的分?jǐn)?shù)階運(yùn)動(dòng)微分方程。如果廣義非勢(shì)力Q″s=0，則方程(21)就成為

(22)

滿足式(20)，則方程(22)為Caputo導(dǎo)數(shù)下的含時(shí)滯的分?jǐn)?shù)階Euler-Lagrange方程。

3 含時(shí)滯的分?jǐn)?shù)階Hamilton作用量變分

含時(shí)滯的分?jǐn)?shù)階Hamilton作用量為

(23)

引入r參數(shù)的有限變換群的無(wú)限小變換

(24)

其展開(kāi)式為

(25)

(26)

(27)

(28)

對(duì)(28)式的第五項(xiàng)，第六項(xiàng)進(jìn)行變量替換t=θ+τ，并考慮邊界條件(13)，得到

(29)

注意到關(guān)系式

(30)

并考慮到分部積分公式(8)和(9)以及(6)和(7)，則(29)式可變?yōu)?/p>

(31)

滿足條件(20)。其中

(32)

式(29)和(31)是Caputo導(dǎo)數(shù)下的含時(shí)滯的分?jǐn)?shù)階Hamilton作用量變分的兩個(gè)基本公式。

4 含時(shí)滯的分?jǐn)?shù)階Noether對(duì)稱性

首先，給出Caputo導(dǎo)數(shù)下的含時(shí)滯的分?jǐn)?shù)階Noether對(duì)稱變換的定義和判據(jù)。

定義1 如果含時(shí)滯的分?jǐn)?shù)階Hamilton作用量(23)，在無(wú)限小群變換(24)作用下，滿足條件

(33)

則稱無(wú)限小變換為含時(shí)滯的分?jǐn)?shù)階Noether對(duì)稱變換。

由定義1和公式(29)和(31)，可得判據(jù)：

判據(jù)1 對(duì)于無(wú)限小變換(24)，當(dāng)t1≤t≤t2-τ時(shí)，滿足條件

(34)

當(dāng)t2-τ

(35)

則變換(24)是含時(shí)滯的力學(xué)系統(tǒng)的分?jǐn)?shù)階Noether對(duì)稱變換．

式(34)和(35)可表為：當(dāng)t1≤t≤t2-τ時(shí)，有

(36)

當(dāng)t2-τ

(37)

當(dāng)r=1時(shí)，式(36)和(37)稱為含時(shí)滯的力學(xué)系統(tǒng)的分?jǐn)?shù)階Noether等式。

其次，研究Caputo導(dǎo)數(shù)下的含時(shí)滯的力學(xué)系統(tǒng)的分?jǐn)?shù)階Noether準(zhǔn)對(duì)稱性。

設(shè)L1是某個(gè)另外的Lagrange函數(shù)，如果變換(24)精確到一階小量滿足如下關(guān)系

(38)

那么稱這種不變性為含時(shí)滯的分?jǐn)?shù)階Hamilton作用量(23)在無(wú)限小變換(24)下的準(zhǔn)不變性。由此確定的L1與L具有相同的運(yùn)動(dòng)微分方程，因而變換(24)可稱為含時(shí)滯的力學(xué)系統(tǒng)的分?jǐn)?shù)階準(zhǔn)對(duì)稱變換。于是有

定義2 如果含時(shí)滯的分?jǐn)?shù)階Hamilton作用量(23)，在無(wú)限小群變換(24)作用下，滿足條件

(39)

由定義2和公式(29)和(31)，可得判據(jù)：

判據(jù)2 對(duì)于無(wú)限小變換(24)，當(dāng)t1≤t≤t2-τ時(shí)，滿足條件

(40)

當(dāng)t2-τ

(41)

則變換(24)是含時(shí)滯的力學(xué)系統(tǒng)的分?jǐn)?shù)階Noether準(zhǔn)對(duì)稱變換．

式(40)和(41)可表為：當(dāng)t1≤t≤t2-τ時(shí),

(42)

當(dāng)t2-τ

(43)

其中ΔG=εσGσ。當(dāng)r=1時(shí)，式(42)和(43)稱為含時(shí)滯的力學(xué)系統(tǒng)的分?jǐn)?shù)階Noether等式。

最后，討論Caputo導(dǎo)數(shù)下的含時(shí)滯的非保守系統(tǒng)的分?jǐn)?shù)階廣義Noether準(zhǔn)對(duì)稱性。

假設(shè)Caputo導(dǎo)數(shù)下的含時(shí)滯的非保守力學(xué)系統(tǒng)受到廣義非勢(shì)力Q″s的作用，如果精確到一階小量滿足如下條件

(44)

則相應(yīng)不變性稱為含時(shí)滯的分?jǐn)?shù)階Hamilton作用量(23)在無(wú)限小變換(24)下的廣義準(zhǔn)不變性，而變換(24)稱為力學(xué)系統(tǒng)的含時(shí)滯的分?jǐn)?shù)階廣義準(zhǔn)對(duì)稱變換。于是有

定義3 如果含時(shí)滯的分?jǐn)?shù)階Hamilton作用量(23)，在無(wú)限小群變換(24)作用下，滿足條件

(45)

則稱無(wú)限小變換為含時(shí)滯的力學(xué)系統(tǒng)的分?jǐn)?shù)階Noether廣義準(zhǔn)對(duì)稱變換。

由定義3和公式(29)和(31)，可得判據(jù)：

判據(jù)3 對(duì)于無(wú)限小變換(24)，當(dāng)t1≤t≤t2-τ時(shí)，滿足條件

(46)

當(dāng)t2-τ

(47)

則變換(24)是含時(shí)滯的力學(xué)系統(tǒng)的分?jǐn)?shù)階廣義Noether準(zhǔn)對(duì)稱變換。

式(46)和(47)可表為：當(dāng)t1≤t≤t2-τ時(shí),

(48)

當(dāng)t2-τ

(49)

當(dāng)r=1時(shí)，式(48)和(49)稱為含時(shí)滯的力學(xué)系統(tǒng)的分?jǐn)?shù)階Noether等式。

利用判據(jù)1-判據(jù)3或含時(shí)滯的分?jǐn)?shù)階Noether等式(36)和(37)，(42)和(43)，(48)和(49)可以判斷含時(shí)滯的力學(xué)系統(tǒng)的分?jǐn)?shù)階Noether對(duì)稱性。

5 含時(shí)滯的分?jǐn)?shù)階Noether定理

本節(jié)我們研究Caputo導(dǎo)數(shù)下的含時(shí)滯的力學(xué)系統(tǒng)的分?jǐn)?shù)階守恒量。首先給出所論含時(shí)滯的力學(xué)系統(tǒng)的分?jǐn)?shù)階守恒量的定義。

(50)

對(duì)于含時(shí)滯的Lagrange系統(tǒng)(22)，如果能找到系統(tǒng)的分?jǐn)?shù)階Noether對(duì)稱變換或分?jǐn)?shù)階Noether準(zhǔn)對(duì)稱變換，便可求得相應(yīng)的分?jǐn)?shù)階守恒量。于是，有如下定理。

定理1 對(duì)于含時(shí)滯的Lagrange系統(tǒng)(22)，如果無(wú)限小變換(24)是定義1下的分?jǐn)?shù)階Noether對(duì)稱變換，則系統(tǒng)存在r個(gè)線性獨(dú)立的分?jǐn)?shù)階守恒量，當(dāng)t1≤t≤t2-τ時(shí)，形如

(51)

當(dāng)t2-τ

(52)

證明 由于無(wú)限小變換(24)是系統(tǒng)的分?jǐn)?shù)階Noether對(duì)稱變換，由定義1，以及式(31)，并將方程(31)代入式(33)，由積分區(qū)間的任意性和εσ的獨(dú)立性，并利用含時(shí)滯的Lagrange方程(22)，得到，當(dāng)t1≤t≤t2-τ時(shí)，有

(53)

當(dāng)t2-τ

(54)

對(duì)(53)和(54)式進(jìn)行積分，便得到結(jié)果。

定理2 對(duì)于含時(shí)滯的Lagrange系統(tǒng)(22)，如果無(wú)限小變換(24)是定義2下的分?jǐn)?shù)階Noether準(zhǔn)對(duì)稱變換，則系統(tǒng)存在r個(gè)線性獨(dú)立的分?jǐn)?shù)階守恒量，當(dāng)t1≤t≤t2-τ時(shí)，形如

Gσ=const.

(55)

當(dāng)t2-τ

(56)

證明 由于無(wú)限小變換(22)是系統(tǒng)的Noether準(zhǔn)對(duì)稱變換，由定義2，以及式(39)，并將方程(31)代入式(39)，由積分區(qū)間的任意性和εσ的獨(dú)立性，并利用含時(shí)滯的Lagrange方程(22)，易知定理成立。證畢。

下面，我們進(jìn)一步討論含時(shí)滯的非保守力學(xué)系統(tǒng)的分?jǐn)?shù)階Noether定理。

定理3 對(duì)于含時(shí)滯的非保守系統(tǒng)(21)，如果無(wú)限小變換(24)是定義3下的分?jǐn)?shù)階廣義Noether準(zhǔn)對(duì)稱變換，則系統(tǒng)存在r個(gè)線性獨(dú)立的分?jǐn)?shù)階守恒量(53)和(54)。

證明 由于無(wú)限小變換(24)是系統(tǒng)的Noether廣義準(zhǔn)對(duì)稱變換，由定義3，以及(45)式，并將方程(31)代入式(45)，由積分區(qū)間的任意性和εσ的獨(dú)立性，并利用含時(shí)滯的Lagrange方程(21)，易知定理成立。

定理1-3稱為含時(shí)滯的非保守力學(xué)系統(tǒng)的分?jǐn)?shù)階Noether定理。由Noether定理可知，對(duì)于所論分?jǐn)?shù)階模型下的含時(shí)滯的非保守系統(tǒng)，如果能找到系統(tǒng)的一個(gè)含時(shí)滯的分?jǐn)?shù)階Noether對(duì)稱變換，便有可能得到系統(tǒng)的一個(gè)分?jǐn)?shù)階守恒量。

6 算 例

例 已知力學(xué)系統(tǒng)的Lagrange函數(shù)和非勢(shì)廣義力為

(57)

其中質(zhì)量m及阻尼系數(shù)c均為常數(shù)。且t∈[t1,t2]，時(shí)滯常量τ

系統(tǒng)的運(yùn)動(dòng)微分方程給出

(58)

由含時(shí)滯的分?jǐn)?shù)階Noether等式(48)和(49)給出

(59)

方程(59)有解

(60)

(61)

生成元(60)和(61)都相應(yīng)于系統(tǒng)的Noether廣義準(zhǔn)對(duì)稱變換．由定理3，得到

I1=0

(62)

(63)

因此，生成元(62)相應(yīng)的分?jǐn)?shù)階守恒量是平庸的．若分?jǐn)?shù)階導(dǎo)數(shù)不存在時(shí)，方程(63)就成為含時(shí)滯的運(yùn)動(dòng)微分方程

(64)

式(63)就成為

(65)

式(65)是含時(shí)滯的非保守系統(tǒng)的廣義Noether準(zhǔn)對(duì)稱性相應(yīng)的守恒量。若時(shí)滯常量τ=0時(shí)，式(65)就成為經(jīng)典力學(xué)系統(tǒng)的相應(yīng)的守恒量

(66)

7 結(jié) 論

提出并研究了Caputo導(dǎo)數(shù)下含時(shí)滯的非保守動(dòng)力學(xué)系統(tǒng)的分?jǐn)?shù)階Noether對(duì)稱性與守恒量。建立了含時(shí)滯的分?jǐn)?shù)階Hamilton原理(15)和(16)，并由此進(jìn)一步導(dǎo)出了含時(shí)滯的分?jǐn)?shù)階Lagrange方程(21)。給出了含時(shí)滯的分?jǐn)?shù)階Hamilton作用量變分的兩個(gè)基本公式，建立了含時(shí)滯的分?jǐn)?shù)階Noether對(duì)稱性的定義和判據(jù)，并得到了相應(yīng)的分?jǐn)?shù)階Noether定理。文章的方法和結(jié)果具有普遍性，可進(jìn)一步應(yīng)用于Caputo導(dǎo)數(shù)下含時(shí)滯的分?jǐn)?shù)階非完整力學(xué)系統(tǒng)，分?jǐn)?shù)階Hamilton系統(tǒng)以及分?jǐn)?shù)階Birkhoff系統(tǒng)等。

[1] OLDHAM K B, SPANIER J. The fractional calculus[M]. San Diego: Academic Press, 1974.

[2] SAMKO S G, KILBAS A A, MARICHER O I. Fractional integrals and derivatives: theory and applications[M]. Breach, 1993.

[3] PODLUBNY I. Fractional differential equations[M]. San Diego: Academic Press，1999.

[4] RIEWE F. Nonconservation lagrangian and Hamiltonian mechanics[J]. Phys Rev E, 1996,53(2):1890-1899.

[5] RIEWE F. Mechanics with fractional derivatives[J]. Phys Rev E, 1997，55(3): 3581-3592.

[6] AGRAWALl O P. Formulation of Euler-Lagrange equations for fractional variational problems [J]. J Math Anal Appl, 2002, 272(1): 368-379.

[7] AGRAWALl O P. Fractional variational calculus in terms of Riesz fractional derivatives[J]. J Phys A: Math Theor, 2007, 40(24): 6283-6303.

[8] FREDERICO G S F, TORRES D F M. Fractional conservation laws in optimal control theory[J]. Nonlinear Dyn, 2008, 53(3):215-222.

[9] FREDERICO G S F, TORRES D F M. Fractional optimal control in the sense of caputo and the fractional Noether’s theorem[J]. Int Math Forum, 2008, 3 (10): 479 - 493..

[11] ZHANG Y, ZHAI X H. Noether symmetries and conserved quantities for fractional Birkhoffian systems[J]. Nonlinear Dyn, 2015, 81(1/2):469-480.

[12] BALEANU D. Fractional Hamilton formalism within Caputo’s derivative[J]. Czech J Phys, 2006, 56 (10/11); 1087-1092.

[13] MALINOWSKA A B, TORRES D F M. Fractional variational calculus in terms of a combined Caputo derivative[J]. Fract Calcu Appl Anal, 2011, 14(4): 523-537.

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

[15] 龍梓軒,張毅.基于正弦周期律拓展的分?jǐn)?shù)階積分的變分問(wèn)題的Noether定理[J].中山大學(xué)學(xué)報(bào):自然科學(xué)版, 2013, 52(5): 51-56.

[16] 何勝鑫,朱建清.基于分?jǐn)?shù)階模型的相空間中非保守力學(xué)系統(tǒng)的Noether準(zhǔn)對(duì)稱性[J].中山大學(xué)學(xué)報(bào):自然科學(xué)版, 2015, 54(4): 5-11.

[17] BALEANU D, MARAABA T, JARAD F. Fractional variational principles with delay[J]. J Phys A:Math Theor, 2008, 41(31): 315403(1-8).

[18] JARAD F, ABDELJAWAD T, BALEANU D. Fractional variational principles with delay within Caputo derivatives[J]. Math Phys , 2010, 1(65): 17-28.

[19] JARAD F, ABDELJAWAD T, BALEANU D. Fractional variational optimal control problems with delayed argument[J]. Nonlinear Dyn, 2010, 62(3): 609-614.

[20] FREDERICO G S F, TORRES D F M. Noether’s symmetry theorem for variational and optimal control problem with time delay[J]. Num Algebra, Contr Optim, 2012, 2(3): 619-630.

[21] 張毅,金世欣.含時(shí)滯的非保守系統(tǒng)動(dòng)力學(xué)的Noether對(duì)稱性[J].物理學(xué)報(bào),2013, 62(23): 214502.

[22] JIN S X, ZHANG Y. Noether symmetry and conserved quantity for a Hamilton system with time delay[J]. Chin Phys B, 2014, 23(5): 054501.

[23] 金世欣,張毅.相空間中含時(shí)滯的非保守力學(xué)系統(tǒng)的Noether定理[J].中山大學(xué)學(xué)報(bào):自然科學(xué)版, 2014, 53(4): 18-23.

[24] ZHAI X H, ZHANG Y. Noether symmetries and conserved quantities for Birkhoffian systems with time delay[J]. Nonlinear Dyn, 2014, 77(1/2): 73-86.

[25] JIN S X, ZHANG Y. Noether theorem for non-conservative Lagrange systems with time delay based on fractional model[J]. Nonlinear Dyn, 2015, 79(2): 1169-1183.

Noether Symmetries for Non-Conservative Lagrange Systems with Time Delay Based on Caputo Fractional Derivative

JINShixin1，ZHANGYi2

(1.College of Physics, Nanjing University of Science and Technology, Nanjing 210094, China；2.College of Civil Engineering, Suzhou University of Science and Technology, Suzhou 215011, China)

The Noether symmetries and the conserved quantities of a mechanical system with time delay based on Caputo fractional derivatives are proposed and studied. Firstly, the fractional Lagrange equations with time delay are established. Secondly, based upon the invariance of the fractional Hamilton action with time delay under the group of infinitesimal transformations, the fractional Noether symmetric transformations, the definitions and criteria of the Noether quasi-symmetric transformations and generalized Noether quasi-symmetric transformations with time delay are given. Finally, the relationship between the fractional symmetries and the fractional conserved quantities with time delay are studied. At the end, an example is given to illustrate the application of the results.

nonconservative system; time delay; Caputo fractional derivative; Noether symmetry; conserved quantity

10.13471/j.cnki.acta.snus.2015.05.011

2015-02-03

國(guó)家自然科學(xué)基金資助項(xiàng)目(10972151，11272227)

金世欣(1987年生)，男；研究方向：一般力學(xué)與力學(xué)基礎(chǔ)；通訊作者：張毅;E-mail:weidiezh@gmail.com

O316

A

0529-6579(2015)05-0049-08