Symmetries of Fractional Constrained Hamiltonian System Described by the Singular Lagrangian

Cai WANG; Chuanjing SONG

doi:10.1051/wujns/2024296547

Open Access

Issue		Wuhan Univ. J. Nat. Sci. Volume 29, Number 6, December 2024


Page(s)		547 - 557
DOI		https://doi.org/10.1051/wujns/2024296547
Published online		07 January 2025

Wuhan University Journal of Natural Sciences, 2024, Vol.29 No.6, 547-557

Mathematics

CLC number: O316

Symmetries of Fractional Constrained Hamiltonian System Described by the Singular Lagrangian

由奇异Lagrange函数描述的约束Hamilton系统的对称性

Cai WANG (王采) and Chuanjing SONG (宋传静)^†

School of Mathematical Sciences, Suzhou University of Science and Technology, Suzhou 215009, Jiangsu, China

^† Corresponding author. E-mail: This email address is being protected from spambots. You need JavaScript enabled to view it.

Received: 28 September 2023

Abstract

Singular systems within combined fractional derivatives are established. Firstly, the fractional Lagrange equation is analyzed. Secondly, the fractional primary constraint is given. Thirdly, the Noether and Lie symmetry methods of fractional constrained Hamiltonian system are studied. Finally, the obtained results are illustrated with an example.

摘要

基于联合分数阶导数建立了奇异系统。首先，分析了分数阶Lagrange方程；其次，给出了分数阶初级约束；第三，研究了分数阶约束Hamilton系统的Noether对称性和Lie对称性方法；最后，举例对所得结果进行阐述。

Key words: combined fractional derivative / constrained Hamilton equation / conserved quantity / Noether symmetry / Lie symmetry

关键字 : 联合分数阶导数 / 约束Hamilton方程 / 守恒量 / Noether对称性 / Lie对称性

Cite this article: WANG Cai, SONG Chuanjing. Symmetries of Fractional Constrained Hamiltonian System Described by the Singular Lagrangian[J]. Wuhan Univ J of Nat Sci, 2024, 29(6): 547-557.

Biography: WANG Cai, male, Master candidate, research direction: analytical mechanics. E-mail: This email address is being protected from spambots. You need JavaScript enabled to view it.

Foundation item: Supported by the National Natural Science Foundation of China (12172241, 12272248) and the Qing Lan Project of Colleges and Universities in Jiangsu Province

© Wuhan University 2024

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

0 Introduction

With the continuous development of mathematics and physics, fractional calculus has been widely discussed in recent years. Fractional-order model enables us to describe complex dynamics and physical behaviors better. Fractional calculus has important applications in fluid mechanics, nuclear magnetic resonance imaging, mechanics of complex viscoelastic materials, and so on^[1-5]. Riewe^[6,7] introduced fractional derivatives to deal with dissipative forces. After that, research on fractional derivatives has mainly focused on the left and right Riemann-Liouville, the left and right Caputo, the Riesz-Riemann-Liouville and the Riesz-Caputo fractional derivatives. Recently, a more general definition of fractional derivative, which is called combined fractional derivative, was introduced. For example, Malinowska and Torres^[8] studied fractional calculus of variations based on a combined Caputo fractional derivative. Zhang^[9] established the fractional differential equations of motion using a combined Riemann-Liouville fractional derivative. In 2015, Luo et al^[10] presented the dynamics of a Birkhoffian system with both combined Caputo fractional derivatives and combined Riemann-Liouville fractional derivatives. Specifically, the four derivatives mentioned above are all special cases of the combined fractional derivative. Therefore, the combined fractional derivative is general.

The singular system refers to the procedure described by the singular Lagrangian. The singular system plays a vital role in modern quantum field theory, such as gravity theory, string (membrane) field theory, Yang-Mills theory, supersymmetry, and supergravity^[11]. When a singular system is described by the canonical variables, there are several inherent constraints. In this case, the system is called the constrained Hamiltonian system, which is essential in condensed matter theory, gauge field theory and many other aspects^[12-14].

Noether symmetry and Lie symmetry are two commonly used methods to study constrained mechanics systems. Noether symmetry was put forward by Noether^[15]. Lie symmetry was introduced by Lutzky^[16]. Then Noether and Lie symmetry were studied further^[17,18]. Particularly, for singular systems, Lie symmetry and conserved quantity based on Riemann-Liouville fractional derivative^[19], Noether symmetry and conserved quantity based on mixed order derivative and Caputo fractional derivative^[20] are studied. Besides, based on the combined fractional derivative, Noether symmetry and conserved quantity for the Birkhoffian system are studied^[21]. In this paper, we intend to study Noether symmetry and Lie symmetry for the singular system based on the combined fractional derivative.

Section 1 gives the preliminaries. In Section 2, the fractional Lagrange equation is analyzed using the combined fractional derivative. In Section 3 and Section 4, the fractional primary constraint is established and the Hamilton equation of the fractional constraint is given. In Section 5 and Section 6, the fractional Noether symmetry, Lie symmetry, and conserved quantities are studied, respectively. An example is given in Section 7 and the conclusion is given in Section 8.

1 Preliminaries

The following preliminaries are about the combined fractional derivatives, which can be seen in Ref. [8], Ref. [22] and Ref. [23]. The combined fractional derivative includes the combined Riemann-Liouville fractional derivative and the combined Caputo fractional derivative.

The combined Riemann-Liouville fractional derivative and combined Caputo fractional derivative are

$^{R L} D_{γ}^{α, β} f (t) = γ_{t_{1}}^{R L} D_{t}^{α} f (t) + {(- 1)}^{n} {(1 - γ)}_{t}^{R L} D_{t_{2}}^{β} f (t)$ Mathematical equation (1)

$^{C} D_{γ}^{α, β} f (t) = γ_{t_{1}}^{C} D_{t}^{α} f (t) + {(- 1)}^{n} {(1 - γ)}_{t}^{C} D_{t_{2}}^{β} f (t)$ Mathematical equation (2)

where $_{t_{1}}^{R L} D_{t}^{α} f (t)$ Mathematical equation , $_{t}^{R L} D_{t_{2}}^{β} f (t)$ , $_{t_{1}}^{C} D_{t}^{α} f (t)$ and $_{t}^{C} D_{t_{2}}^{β} f (t)$ are the left and right Riemann-Liouville and Caputo fractional derivatives, $t \in [t_{1}, t_{2}]$ , $0 \leq γ \leq 1$ , $n - 1 \leq α, β < n$ .

Under the condition $0 < α, β < 1$ Mathematical equation , we have

$_{t_{1}}^{R L} D_{t}^{α} f (t) = {}_{t_{1}}^{C}D_{t}^{α} f (t) - \frac{1}{Γ (1 - α)} \frac{f (t_{1})}{{(t - t_{1})}^{α}}$ Mathematical equation (3)

$_{t}^{R L} D_{t_{2}}^{β} f (t) = {}_{t}^{C}D_{t_{2}}^{β} f (t) + \frac{1}{Γ (1 - β)} \frac{f (t_{2})}{{(t_{2} - t)}^{β}}$ Mathematical equation (4)

Additionally, the fractional partial integral formulae are

$\int_{t_{1}}^{t_{2}} [*] {}_{t_{1}}^{R L}D_{t}^{α} η d t = \int_{t_{1}}^{t_{2}} η {}_{t}^{C}D D_{t_{2}}^{α} [*] d t - {\sum_{j = 0}^{n - 1} {(- 1)}^{n + j} {}_{t_{1}}^{R L}D D_{t}^{α + j - n} η (t) D^{n - 1 - j} [*] |}_{t_{1}}^{t_{2}}$ Mathematical equation (5)

$\int_{t_{1}}^{t_{2}} [*] {}_{t}^{R L}D_{t_{2}}^{β} η d t = \int_{t_{1}}^{t_{2}} η {}_{t_{1}}^{C}D D_{t}^{β} [*] d t - {\sum_{j = 0}^{n - 1} {}_{t}^{R L} D_{t_{2}}^{β + j - n} η (t) D^{n - 1 - j} [*] |}_{t_{1}}^{t_{2}}$ Mathematical equation (6)

$\int_{t_{1}}^{t_{2}} [*] {}_{t_{1}}^{C}D_{t}^{α} η d t = \int_{t_{1}}^{t_{2}} η {}_{t}^{R L}D D_{t_{2}}^{α} [*] d t + {\sum_{j = 0}^{n - 1} {}_{t}^{R L} D_{t_{2}}^{α + j - n} [*] D^{n - 1 - j} η (t) |}_{t_{1}}^{t_{2}}$ Mathematical equation (7)

$\int_{t_{1}}^{t_{2}} [*] {}_{t_{1}}^{C}D_{t}^{β} η d t = \int_{t_{1}}^{t_{2}} η {}_{t_{1}}^{R L}D D_{t}^{β} [*] d t + {\sum_{j = 0}^{n - 1} {(- 1)}^{n + j} {}_{t_{1}}^{R L}D D_{t}^{β + j - n} [*] D^{n - 1 - j} η (t) |}_{t_{1}}^{t_{2}}$ Mathematical equation (8)

$\int_{t_{1}}^{t_{2}} [*] {}_{t_{1}}^{R}D_{t_{2}}^{α} η d t = {(- 1)}^{n} \int_{t_{1}}^{t_{2}} η {}_{t_{1}}^{R C}D D_{t_{2}}^{α} [*] d t - {\sum_{j = 0}^{n - 1} {(- 1)}^{n + j} {}_{t_{1}}^{R}D D_{t_{2}}^{α + j - n} η (t) D^{n - 1 - j} [*] |}_{t_{1}}^{t_{2}}$ Mathematical equation (9)

$\int_{t_{1}}^{t_{2}} [*] {}_{t_{1}}^{R C}D_{t_{2}}^{α} η d t = {(- 1)}^{n} \int_{t_{1}}^{t_{2}} η {}_{t_{1}}^{R}D D_{t_{2}}^{α} [*] d t + {\sum_{j = 0}^{n - 1} {(- 1)}^{j} {}_{t_{1}}^{R}D D_{t_{2}}^{α + j - n} [*] D^{n - 1 - j} η (t) |}_{t_{1}}^{t_{2}}$ Mathematical equation (10)

where $D = d / d t$ Mathematical equation represents the integer order derivative.

It should be noted that in this paper, we set $0 < α, β < 1$ Mathematical equation .

2 Fractional Lagrange Equation

In this section, the fractional variational problems under combined fractional derivatives are studied.

We give a function

$I_{R L} [q_{R L} (\cdot)] = \int_{t_{1}}^{t_{2}} L_{R L} (t, q_{R L}, {}^{R L}D D_{γ}^{α, β} q_{R L}) d t, q_{R L} (t_{1}) = q_{R L 1}, q_{R L} (t_{2}) = q_{R L 2},$ Mathematical equation (11)

where $q_{R L} = (q_{R L 1}, q_{R L 2}, \dots, q_{R L n})$ Mathematical equation , ${}^{R L}_{γ}^{α, β} q_{R L} = ({}^{R L} D_{γ}^{α, β} q_{R L 1}, {}^{R L}D D_{γ}^{α, β} q_{R L 2}, \dots,$ ${}^{R L} D_{γ}^{α, β} q_{R L n})$ , $q_{R L 1} = (q_{R L 11}, q_{R L 12}, \dots, q_{R L 1 n})$ , $q_{R L 2} = (q_{R L 21}, q_{R L 22}, \dots, q_{R L 2 n})$ Mathematical equation , and the Lagrangian $L_{R L}$ is assumed to be a $C^{2}$ function. If $q_{R L} (\cdot)$ is an extremal of Eq. (11), we obtain

${\frac{d}{d ε_{R L}} I_{R L} [q_{R L} + ε_{R L} h_{R L}] |}_{ε_{R L} = 0} = 0,$ Mathematical equation (12)

where $ε_{R L}$ Mathematical equation is a small parameter, $h_{R L} = (h_{R L 1}, h_{R L 2}, \dots, h_{R L n})$ , $h_{R L} (t_{1}) = h_{R L} (t_{2}) = 0$ . Further, we get

$\int_{t_{1}}^{t_{2}} (\frac{\partial L_{R L}}{\partial q_{R L k}} \cdot h_{R L k} + \frac{\partial L_{R L}}{\partial {}^{R L}D D_{γ}^{α, β} q_{R L k}} \cdot {}^{R L}D D_{γ}^{α, β} h_{R L k}) d t = 0, k = 1,2, \dots, n .$ Mathematical equation (13)

Using fractional partial integral formulae (Eqs. (5) and (6)) in the second terms of Eq. (13), we obtain

$\begin{array}{l} \int_{t_{1}}^{t_{2}} \frac{\partial L_{R L}}{\partial {}^{R L}D D_{γ}^{α, β} q_{R L k}} \cdot {}^{R L}D D_{γ}^{α, β} h_{R L k} d t = - \int_{t_{1}}^{t_{2}} (h_{R L k} \cdot {}^{C}D D_{1 - γ}^{β, α} \frac{\partial L_{R L}}{\partial^{R L} D_{γ}^{α, β} q_{R L k}}) d t + \frac{γ}{Γ (1 - α)} \int_{t_{1}}^{t_{2}} {(t_{2} - t)}^{- α} h_{R L k} d t \cdot \frac{\partial L_{R L} (t_{2})}{\partial^{R L} D_{γ}^{α, β} q_{R L k}} \\ - \frac{1 - γ}{Γ (1 - β)} \int_{t_{1}}^{t_{2}} {(t - t_{1})}^{- β} h_{R L k} d t \cdot \frac{\partial L_{R L} (t_{1})}{\partial^{R L} D_{γ}^{α, β} q_{R L k}} \end{array}$ Mathematical equation (14)

where $\frac{\partial L_{R L} (t_{1})}{\partial^{R L} D_{γ}^{α, β} q_{R L k}} = \frac{\partial L_{R L} (t_{1}, q_{R L} (t_{1}),^{R L} D_{γ}^{α, β} q_{R L} (t_{1}))}{\partial^{R L} D_{γ}^{α, β} q_{R L k}}$ Mathematical equation , $\frac{\partial L_{R L} (t_{2})}{\partial^{R L} D_{γ}^{α, β} q_{R L k}} = \frac{\partial L_{R L}}{\partial^{R L} D_{γ}^{α, β} q_{R L k}}$ $(t_{2}, q_{R L} (t_{2}),^{R L} D_{γ}^{α, β} q_{R L} (t_{2}))$ . Substituting Eq. (14) into Eq. (13), we obtain

$\int_{t_{1}}^{t_{2}} [\frac{\partial L_{R L}}{\partial q_{R L k}} - {}^{C}D D_{1 - γ}^{β, α} \frac{\partial L_{R L}}{\partial^{R L} D_{γ}^{α, β} q_{R L k}} + \frac{γ {(t_{2} - t)}^{- α}}{Γ (1 - α)} \frac{\partial L_{R L} (t_{2})}{\partial^{R L} D_{γ}^{α, β} q_{R L k}} - \frac{(1 - γ) {(t - t_{1})}^{- β}}{Γ (1 - β)} \frac{\partial L_{R L} (t_{1})}{\partial^{R L} D_{γ}^{α, β} q_{R L k}}] \cdot h_{R L k} d t = 0 .$ Mathematical equation (15)

According to the fundamental lemma of variational calculation^[24], we can deduce that

$\frac{\partial L_{R L}}{\partial q_{R L k}} - {}^{C}D_{1 - γ}^{β, α} \frac{\partial L_{R L}}{\partial^{R L} D_{γ}^{α, β} q_{R L k}} + \frac{γ {(t_{2} - t)}^{- α}}{Γ (1 - α)} \frac{\partial L_{R L} (t_{2})}{\partial^{R L} D_{γ}^{α, β} q_{R L k}} - \frac{(1 - γ) {(t - t_{1})}^{- β}}{Γ (1 - β)} \frac{\partial L_{R L} (t_{1})}{\partial^{R L} D_{γ}^{α, β} q_{R L k}} = 0, k = 1,2, \dots, n .$ Mathematical equation (16)

Eq. (16) is called the fractional Lagrange equation within the combined Riemann-Liouville fractional derivative.

The fractional Lagrange equation has another form, which can be obtained from Eqs. (3), (4), (7), (8) and (13) as

$\frac{\partial L_{R L}}{\partial q_{R L k}} - {}^{R L}D_{1 - γ}^{β, α} \frac{\partial L_{R L}}{\partial {}^{R L}D D_{γ}^{α, β} q_{R L k}} = 0, k = 1,2, \dots \cdot n .$ Mathematical equation (17)

Similarly, the fractional Lagrange equation within the combined Caputo fractional derivative can also be obtained as

$\frac{\partial L_{C}}{\partial q_{C k}} - {}^{R L}D_{1 - γ}^{β, α} \frac{\partial L_{C}}{\partial^{C} D_{γ}^{α, β} q_{C k}} = 0, k = 1,2, \dots, n,$ Mathematical equation (18)

where the Lagrangian $L_{C} = L_{C} (t, q_{C}, {}^{C}D D_{γ}^{α, β} q_{C})$ Mathematical equation , $L_{C}$ is assumed to be a $C^{2}$ function, $q_{C} = (q_{C 1}, q_{C 2}, \dots, q_{C n})$ , ${}^{C}_{γ}^{α, β} q_{C} = ({}^{C} D_{γ}^{α, β} q_{C 1}, {}^{C}D D_{γ}^{α, β} q_{C 2}, \dots, {}^{C}D D_{γ}^{α, β} q_{C n})$ .

Remark 1 Eqs. (16) and (18) are general and universal, so we can select different values of $γ$ Mathematical equation to get different results.

3 Fractional Primary Constraint

This section presents fractional primary constraints within the combined Riemann-Liouville fractional derivative and combined Caputo fractional derivative.

Define fractional generalized momentum and Hamiltonian as

$p_{R L k} = \frac{\partial L_{R L} (t, q_{R L}, {}^{R L}D D_{γ}^{α, β} q_{R L})}{\partial {}^{R L}D D_{γ}^{α, β} q_{R L k}},$ Mathematical equation (19)

$H_{R L} = p_{R L k} \cdot {}^{R L}D_{γ}^{α, β} q_{R L k} - L_{R L} (t, q_{R L}, {}^{R L}D D_{γ}^{α, β} q_{R L}), k = 1,2, \dots, n .$ Mathematical equation (20)

Here, we consider that the $L_{R L} (t, q_{R L}, {}^{R L}D D_{γ}^{α, β} q_{R L})$ Mathematical equation is singular, meaning that only a portion of the terms ${}^{R L}_{γ}^{α, β} q_{R L k}$ , $k = 1,2, \dots, n$ can be solved. In this case, we assume the number is $R$ , $0 \leq R < n$ .

Firstly, when $1 \leq R < n$ Mathematical equation , we have

$^{R L} D_{γ}^{α, β} q_{R L σ} = f_{R L}^{σ} (t, q_{R L}, p_{R L E}, {}^{R L}D D_{γ}^{α, β} q_{R L ρ}), σ, E = 1,2, \dots, R, ρ = R + 1, \dots, n,$ Mathematical equation (21)

where $p_{R L E} = (p_{R L 1}, p_{R L 2}, \dots, p_{R L R})$ Mathematical equation , $^{R L} D_{γ}^{α, β} q_{R L ρ} = (^{R L} D_{γ}^{α, β} q_{R L (R + 1)},^{R L} D_{γ}^{α, β} q_{R L (R + 2)}, \dots,$ $^{R L} D_{γ}^{α, β} q_{R L n})$ , $1 \leq R < n$ . Substituting Eq. (21) into Eq. (19), we have

$p_{R L k} = g_{R L k} (t, q_{R L}, f_{R L}^{σ} (t, q_{R L}, p_{R L E}, {}^{R L}D D_{γ}^{α, β} q_{R L ρ}), {}^{R L}D D_{γ}^{α, β} q_{R L ρ}) = g_{R L k} (t, q_{R L}, p_{R L E}, {}^{R L}D D_{γ}^{α, β} q_{R L ρ}), k = 1,2, \dots, n .$ Mathematical equation (22)

If $k = 1,2, \dots, R$ Mathematical equation , Eq. (22) obviously holds. If $k = R + 1, \dots, n$ , we have

$p_{R L ρ} = g_{R L ρ} (t, q_{R L}, p_{R L E}) o r p_{R L F} = g_{R L F} (t, q_{R L}, p_{R L E}), 1 \leq R < n,$ Mathematical equation (23)

where $p_{R L F} = (p_{R L (R + 1)}, \dots, p_{R L n})$ Mathematical equation , $g_{R L F} = (g_{R L (R + 1)}, \dots, g_{R L n})$ , $ρ = R + 1, \dots, n$ . Another form of Eq. (23) is

$ϕ_{R L} (t, q_{R L}, p_{R L}) = p_{R L F} - g_{R L F} (t, q_{R L}, p_{R L E}) = 0,$ Mathematical equation (24)

where $ϕ_{R L} = (ϕ_{R L 1}, ϕ_{R L 2}, \dots, ϕ_{R L (n - R)})$ Mathematical equation , $p_{R L} = (p_{R L 1}, p_{R L 2}, \dots, p_{R L n})$ , $1 \leq R < n$ .

Secondly, when $R = 0$ Mathematical equation , ${}^{R L}_{γ}^{α, β} q_{R L k}$ cannot be solved. Therefore, from Eq. (19), we have

$p_{R L k} = g_{R L k} (t, q_{R L}), k = 1,2, \dots, n .$ Mathematical equation (25)

$ϕ_{R L a} (t, q_{R L}, p_{R L}) = p_{R L a} - g_{R L a} (t, q_{R L}) = 0, a = 1,2, \dots, n .$ Mathematical equation (26)

From Eq. (24) and Eq. (26), we can obtain

$ϕ_{R L a} (t, q_{R L}, p_{R L}) = 0, 0 \leq R < n, a = 1,2, \dots, n - R .$ Mathematical equation (27)

Eq. (27) is called fractional primary constraint within a combined Riemann-Liouville fractional derivative.

Similarly, we can define

$p_{C k} = \frac{\partial L_{C} (t, q_{C}, {}^{C}D D_{γ}^{α, β} q_{C})}{\partial {}^{C}D D_{γ}^{α, β} q_{C k}}, k = 1,2, \dots, n,$ Mathematical equation (28)

$H_{C} = p_{C k} {}^{C}D_{γ}^{α, β} q_{C k} - L_{C} (t, q_{C}, {}^{C}D D_{γ}^{α, β} q_{C}), k = 1,2, \dots, n .$ Mathematical equation (29)

The fractional primary constraint within the combined Caputo fractional derivative can be obtained as

$ϕ_{C a} (t, q_{C}, p_{C}) = 0, 0 \leq R < n, a = 1,2, \dots, n - R, p_{C} = (p_{C 1}, p_{C 2}, \dots, p_{C n}) .$ Mathematical equation (30)

Remark 2 According to the definition of the fractional generalized momenta (Eqs. (19) and (28)) rather than the fractional Euler-Lagrange equations (Eqs. (16) and (18)), the fractional primary constraints (Eqs. (27) and (30)) are obtained.

4 Fractional Constrained Hamilton Equation

Firstly, the fractional constrained Hamilton equation within the combined Riemann-Liouville fractional derivative is studied.

From Eq. (20), we obtain

$δ H_{R L} = δ p_{R L k} \cdot {}^{R L}D_{γ}^{α, β} q_{R L k} - \frac{\partial L_{R L}}{\partial q_{R L k}} δ q_{R L k}, k = 1,2, \dots, n .$ Mathematical equation (31)

Besides, because of the Hamiltonian $H_{R L} = H_{R L} (t, q_{R L}, p_{R L})$ Mathematical equation , we have

$δ H_{R L} = \frac{\partial H_{R L}}{\partial q_{R L k}} \cdot δ q_{R L k} + \frac{\partial H_{R L}}{\partial p_{R L k}} \cdot δ p_{R L k}, k = 1,2, \dots, n .$ Mathematical equation (32)

It follows from Eqs. (31) and (32) that

$({}^{R L} D_{γ}^{α, β} q_{R L k} - \frac{\partial H_{R L}}{\partial p_{R L k}}) δ p_{R L k} - (\frac{\partial L_{R L}}{\partial q_{R L k}} + \frac{\partial H_{R L}}{\partial q_{R L k}}) δ q_{R L k} = 0 .$ Mathematical equation (33)

Making use of Eqs. (16) and (19), the term $\partial L_{R L} / \partial q_{R L k}$ Mathematical equation in Eq. (33) can be replaced by ${}^{C}_{1 - γ}^{β, α} p_{R L k} - \frac{γ {(t_{2} - t)}^{- α}}{Γ (1 - α)} p_{R L k} (t_{2}) + \frac{(1 - γ) {(t - t_{1})}^{- β}}{Γ (1 - β)} p_{R L k} (t_{1})$ , thus we have

$(^{R L} D_{γ}^{α, β} q_{R L k} - \frac{\partial H_{R L}}{\partial p_{R L k}}) \cdot δ p_{R L k} - [{}^{C} D_{1 - γ}^{β, α} p_{R L k} - \frac{γ {(t_{2} - t)}^{- α}}{Γ (1 - α)} p_{R L k} (t_{2}) + \frac{(1 - γ) {(t - t_{1})}^{- β}}{Γ (1 - β)} p_{R L k} (t_{1}) + \frac{\partial H_{R L}}{\partial q_{R L k}}] δ q_{R L k} = 0, k = 1,2, \dots, n,$ Mathematical equation (34)

where $p_{R L k} (t_{1}) = \frac{\partial L_{R L} (t_{1})}{\partial^{R L} D_{γ}^{α, β} q_{R L k}}$ Mathematical equation , $p_{R L k} (t_{2}) = \frac{\partial L_{R L} (t_{2})}{\partial^{R L} D_{γ}^{α, β} q_{R L k}}$ . When the system (Eq. (16)) is singular, we have

$λ_{R L a} \frac{\partial ϕ_{R L a}}{\partial q_{R L k}} \cdot δ q_{R L k} + λ_{R L a} \frac{\partial ϕ_{R L a}}{\partial p_{R L k}} \cdot δ p_{R L k} = 0 .$ Mathematical equation (35)

It follows from Eqs. (34) and (35) that

${}^{C}_{1 - γ}^{β, a} p_{R L k} = - \frac{\partial H_{R L}}{\partial q_{R L k}} + \frac{γ {(t_{2} - t)}^{- α}}{Γ (1 - α)} p_{R L k} (t_{2}) - \frac{(1 - γ) {(t - t_{1})}^{- β}}{Γ (1 - β)} p_{R L k} (t_{1}) - λ_{R L a} \frac{\partial ϕ_{R L a}}{\partial q_{R L k}},$ Mathematical equation

${}^{R L}_{γ}^{α, β} q_{R L k} = \frac{\partial H_{R L}}{\partial p_{R L k}} + λ_{R L a} \frac{\partial ϕ_{R L a}}{\partial p_{R L k}}, a = 1,2, \dots, n - R, 0 \leq R < n, k = 1,2, \dots, n .$ Mathematical equation (36)

Eq. (36) is called fractional constrained Hamilton equation within a combined Riemann-Liouville fractional derivative.

Similarly, we can also get the Hamilton equation with fractional constraints within the combined Caputo fractional derivative

${}^{R L}_{1 - γ}^{α, β} p_{C k} = - \frac{\partial H_{C}}{\partial q_{C k}} - λ_{C a} \frac{\partial ϕ_{C a}}{\partial q_{C k}}, {}^{C}D_{γ}^{α, β} q_{C k} = \frac{\partial H_{C}}{\partial p_{C k}} + λ_{C a} \frac{\partial ϕ_{C a}}{\partial p_{C k}}, k = 1,2, \dots, n, a = 1,2, \dots, n - R, 0 \leq R < n .$ Mathematical equation (37)

Remark 3 Because of the various values of $γ$ Mathematical equation , different results can be obtained from Eqs. (36) and (37).

According to the above method of establishing the Hamilton equation with fractional constraints, it is very important to solve the Lagrange multiplier. Next, we will calculate the Lagrange multiplier by fractional Poisson bracket.

Let $F = F (t, q, p)$ Mathematical equation , $G = G (t, q, p)$ , $q = (q_{1}, q_{2}, \dots, q_{n})$ , $p = (p_{1}, p_{2}, \dots, p_{n})$ , then we define

${F, G} = \frac{\partial F}{\partial q_{k}} \frac{\partial G}{\partial p_{k}} - \frac{\partial F}{\partial p_{k}} \frac{\partial G}{\partial q_{k}}, k = 1,2, \dots, n .$ Mathematical equation (38)

Then we have

$\begin{matrix} ({ϕ_{R L a}, H_{R L}} + λ_{R L b} {ϕ_{R L a}, ϕ_{R L b}}) \cdot \frac{\partial ϕ_{R L a}}{\partial q_{R L j}} {\dot{q}}_{R L j} + (\frac{\partial ϕ_{R L a}}{\partial t} + \frac{\partial ϕ_{R L a}}{\partial p_{R L j}} {\dot{p}}_{R L j}) \times \frac{\partial ϕ_{R L a}}{\partial q_{R L k}} {}^{R L}D_{γ}^{α, β} q_{R L k} - (\frac{\partial ϕ_{R L a}}{\partial p_{R L k}} \cdot \frac{\partial ϕ_{R L a}}{\partial q_{R L j}} {\dot{q}}_{R L j}) [{}^{C} D_{1 - γ}^{β, α} p_{R L k} - \frac{γ {(t_{2} - t)}^{- α}}{Γ (1 - α)} p_{R L k} (t_{2}) + \frac{(1 - γ) {(t - t_{1})}^{- β}}{Γ (1 - β)} p_{R L k} (t_{1})] = 0, a, b = 1,2, \dots, n - R, k, j = 1,2, \dots, n, 0 \leq R < n, \end{matrix}$ Mathematical equation (39)

and

$({ϕ_{C a}, H_{C}} + λ_{C b} {ϕ_{C a}, ϕ_{C b}}) \frac{\partial ϕ_{C a}}{\partial q_{C j}} {\dot{q}}_{C j} + (\frac{\partial ϕ_{C a}}{\partial t} + \frac{\partial ϕ_{C a}}{\partial p_{C j}} \cdot {\dot{p}}_{C j}) \times \frac{\partial ϕ_{C a}}{\partial q_{C k}} {}^{C}D_{γ}^{α, β} q_{C k} - \frac{\partial ϕ_{C a}}{\partial p_{C k}} \cdot \frac{\partial ϕ_{C a}}{\partial q_{C j}} {\dot{q}}_{C j} \cdot {}^{R L}D_{1 - γ}^{β, α} p_{C k} = 0, a, b = 1,2, \dots, n - R, k, j = 1,2, \dots, n, 0 \leq R < n .$ Mathematical equation (40)

From Eqs. (39) and (40), the Lagrange multipliers can be calculated.

5 Noether Theorem

Definition 1 A quantity $C$ Mathematical equation is called a conserved quantity if and only if $d C / d t = 0$ .

5.1 Noether Theorem within the Combined Riemann-Liouville Fractional Derivative

The Hamilton action within the combined Riemann-Liouville fractional derivative is

$I_{R L} = \int_{t_{1}}^{t_{2}} [p_{R L k} {}^{R L}D D_{γ}^{α, β} q_{R L k} - H_{R L} (t, q_{R L}, p_{R L})] d t .$ Mathematical equation (41)

The infinitesimal transformations are

$\bar{t} = t + Δ t, {\bar{q}}_{R L k} (\bar{t}) = q_{R L k} (t) + Δ q_{R L k}, {\bar{p}}_{R L k} (\bar{t}) = p_{R L k} (t) + Δ p_{R L k}, k = 1,2, \dots, n .$ Mathematical equation (42)

Namely,

$\begin{array}{l} \bar{t} = t + θ_{R L} ξ_{R L 0} (t, q_{R L}, p_{R L}) + o (θ_{R L}), {\bar{q}}_{R L k} (\bar{t}) = q_{R L k} (t) + θ_{R L} ξ_{R L k} (t, q_{R L}, p_{R L}) + o (θ_{R L}), \\ {\bar{p}}_{R L k} (\bar{t}) = p_{R L k} (t) + θ_{R L} η_{R L k} (t, q_{R L}, p_{R L}) + o (θ_{R L}), \end{array}$ Mathematical equation (43)

where $θ_{R L}$ Mathematical equation is a small parameter, $ξ_{R L 0}$ , $ξ_{R L k}$ , and $η_{R L k}$ are infinitesimal generators, and $o (θ_{R L})$ means the higher order of $θ_{R L}$ .

Letting $Δ I_{R L}$ Mathematical equation be the linear part of ${\bar{I}}_{R L} - I_{R L}$ , and without considering the higher order of $θ_{R L}$ , we obtain

$\begin{array}{l} Δ I_{R L} = θ_{R L} \int_{t_{1}}^{t_{2}} [p_{R L k} {}^{R L}D D_{γ}^{α, β} (ξ_{R L k} - {\dot{q}}_{R L k} ξ_{R L 0}) + (p_{R L k} \cdot {}^{R L}D D_{γ}^{α, β} q_{R L k} - H_{R L}) {\dot{ξ}}_{R L 0} + (p_{R L k} \frac{d}{d t} {}^{R L}D D_{γ}^{α, β} q_{R L k} - \frac{\partial H_{R L}}{\partial t}) ξ_{R L 0} \\ - \frac{\partial H_{R L}}{\partial q_{R L k}} ξ_{R L k} + λ_{R L a} \frac{\partial ϕ_{R L a}}{\partial p_{R L k}} \cdot η_{R L k} + q_{R L k} (t_{2}) ξ_{R L 0} (t_{2}) \times \frac{(1 - γ) p_{R L k}}{Γ (1 - β)} \frac{d}{d t} {(t_{2} - t)}^{- β} - q_{R L k} (t_{1}) ξ_{R L 0} (t_{1}) \frac{γ p_{R L k}}{Γ (1 - α)} \frac{d}{d t} {(t - t_{1})}^{- α}] d t \end{array}$ Mathematical equation (44)

where $ξ_{R L 0} (t_{1}) = ξ_{R L 0} (t_{1}, q_{R L} (t_{1}), p_{R L} (t_{1}))$ Mathematical equation and $ξ_{R L 0} (t_{2}) = ξ_{R L 0} (t_{2}, q_{R L} (t_{2}), p_{R L} (t_{2}))$ . Let $Δ I_{R L} = 0$ , then Eq. (44) gives

$\begin{array}{l} {p_{R L k}}^{R L} D_{γ}^{α, β} (ξ_{R L k} - {\dot{q}}_{R L k} ξ_{R L 0}) + (p_{R L k} {}^{R L}D D_{γ}^{α, β} q_{R L k} - H_{R L}) {\dot{ξ}}_{R L 0} + λ_{R L a} \frac{\partial ϕ_{R L a}}{\partial p_{R L k}} \cdot η_{R L k} + (p_{R L k} \frac{d}{d t} {}^{R L}D D_{γ}^{α, β} q_{R L k} - \frac{\partial H_{R L}}{\partial t}) ξ_{R L 0} \\ - \frac{\partial H_{R L}}{\partial q_{R L k}} ξ_{R L k} + q_{R L k} (t_{2}) ξ_{R L 0} (t_{2}) \times \frac{(1 - γ) p_{R L k}}{Γ (1 - β)} \frac{d}{d t} {(t_{2} - t)}^{- β} - q_{R L k} (t_{1}) ξ_{R L 0} (t_{1}) \frac{γ p_{R L k}}{Γ (1 - α)} \frac{d}{d t} {(t - t_{1})}^{- α} = 0 \end{array}$ Mathematical equation (45)

Equation (45) is called the fractional Noether identity within the combined Riemann-Liouville fractional derivative.

Theorem 1 If $ξ_{R L 0}$ Mathematical equation , $ξ_{R L k}$ , and $η_{R L k}$ satisfy Eq. (45), then there exists a conserved quantity

$\begin{array}{l} C_{R L} = (p_{R L k} {}^{R L}D D_{γ}^{α, β} q_{R L k} - H_{R L}) ξ_{R L 0} + \int_{t_{1}}^{t} {p_{R L k} {}^{R L}D D_{γ}^{α, β} (ξ_{R L k} - {\dot{q}}_{R L k} ξ_{R L 0}) \\ + (ξ_{R L k} - {\dot{q}}_{R L k} ξ_{R L 0}) [{}^{C} D_{1 - γ}^{β, α} p_{R L k} - \frac{γ {(t_{2} - τ)}^{- α}}{Γ (1 - α)} p_{R L k} (t_{2}) + \frac{(1 - γ) {(τ - t_{1})}^{- β}}{Γ (1 - β)} p_{R L k} (t_{1})]} d τ \\ + q_{R L k} (t_{2}) ξ_{R L 0} (t_{2}) \frac{1 - γ}{Γ (1 - β)} \int_{t_{1}}^{t} p_{R L k} \frac{d}{d τ} {(t_{2} - τ)}^{- β} d τ - q_{R L k} (t_{1}) ξ_{R L 0} (t_{1}) \frac{γ}{Γ (1 - α)} \int_{t_{1}}^{t} p_{R L k} \frac{d}{d τ} {(τ - t_{1})}^{- α} d τ \end{array}$ Mathematical equation (46)

for the system within the combined Riemann-Liouville fractional derivative.

Proof Using Eqs. (27), (36), and (45), we have

$\frac{d C_{R L}}{d t} = (p_{R L k} {}^{R L}D D_{γ}^{α, β} q_{R L k} - H_{R L}) {\dot{ξ}}_{R L 0} + ξ_{R L 0} ({\dot{p}}_{R L k} {}^{R L}D D_{γ}^{α, β} q_{R L k} + p_{R L k} \frac{d}{d t} {}^{R L}D D_{γ}^{α, β} q_{R L k} - \frac{\partial H_{R L}}{\partial t} - \frac{\partial H_{R L}}{\partial q_{R L k}} \cdot {\dot{q}}_{R L k} - \frac{\partial H_{R L}}{\partial p_{R L k}} \cdot {\dot{p}}_{R L k})$ Mathematical equation

$+ {p_{R L k}}^{R L} D_{γ}^{α, β} (ξ_{R L k} - {\dot{q}}_{R L k} ξ_{R L 0}) + (ξ_{R L k} - {\dot{q}}_{R L k} ξ_{R L 0}) [{}^{C} D_{1 - γ}^{β, α} p_{R L k} - \frac{γ {(t_{2} - t)}^{- α}}{Γ (1 - α)} p_{R L k} (t_{2}) + \frac{(1 - γ) {(t - t_{1})}^{- β}}{Γ (1 - β)} p_{R L k} (t_{1})]$ Mathematical equation

$+ q_{R L k} (t_{2}) ξ_{R L 0} (t_{2}) \frac{1 - γ}{Γ (1 - β)} p_{R L k} \frac{d}{d t} {(t_{2} - t)}^{- β} - q_{R L k} (t_{1}) ξ_{R L 0} (t_{1}) \frac{γ}{Γ (1 - α)} p_{R L k} \frac{d}{d t} {(t - t_{1})}^{- α}$ Mathematical equation

$\begin{array}{l} = \frac{\partial H_{R L}}{\partial q_{R L k}} ξ_{R L k} - λ_{R L a} \frac{\partial ϕ_{R L a}}{\partial p_{R L k}} η_{R L k} + ξ_{R L 0} ({\dot{p}}_{R L k} \cdot λ_{R L a} \frac{\partial ϕ_{R L a}}{\partial p_{R L k}} - \frac{\partial H_{R L}}{\partial q_{R L k}} {\dot{q}}_{R L k}) + (ξ_{R L k} - {\dot{q}}_{R L k} ξ_{R L 0}) (- \frac{\partial H_{R L}}{\partial q_{R L k}} - λ_{R L a} \frac{\partial ϕ_{R L a}}{\partial q_{R L k}}) \\ = - λ_{R L a} \frac{\partial ϕ_{R L a}}{\partial p_{R L k}} (η_{R L k} - {\dot{p}}_{R L k} ξ_{R L 0}) - λ_{R L a} \frac{\partial ϕ_{R L a}}{\partial q_{R L k}} (ξ_{R L k} - {\dot{q}}_{R L k} ξ_{R L 0}) = - λ_{R L a} \cdot δ ϕ_{R L a} = 0 . \end{array}$ Mathematical equation

The proof is completed.

5.2 Noether Theorem within the Combined Caputo Fractional Derivative

The Hamilton action within the combined Caputo fractional derivative is

$I_{C} = \int_{t_{1}}^{t_{2}} [p_{C k} \cdot {}^{C}D D_{γ}^{α, β} q_{C k} - H_{C} (t, q_{C}, p_{C})] d t .$ Mathematical equation (47)

The infinitesimal transformations are

$\bar{t} = t + Δ t, {\bar{q}}_{C k} (\bar{t}) = q_{C k} (t) + Δ q_{C k}, {\bar{p}}_{C k} (\bar{t}) = p_{C k} (t) + Δ p_{C k}, k = 1,2, \dots, n .$ Mathematical equation (48)

Namely,

$\bar{t} = t + θ_{C} ξ_{C 0} (t, q_{C}, p_{C}) + o (θ_{C}), {\bar{q}}_{C k} (\bar{t}) = q_{C k} (t) + θ_{C} ξ_{C k} (t, q_{C}, p_{C}) + o (θ_{C}), {\bar{p}}_{C k} (\bar{t}) = p_{C k} (t) + θ_{C} η_{C k} (t, q_{C}, p_{C}) + o (θ_{C})$ Mathematical equation (49)

where $θ_{C}$ Mathematical equation is a small parameter, $ξ_{C 0}$ , $ξ_{C k}$ , and $η_{C k}$ are infinitesimal generators. Letting $Δ I_{C} = {\bar{I}}_{C} - I_{C} = 0$ , we have

$\begin{matrix} p_{C k} {}^{C}D_{γ}^{α, β} (ξ_{C k} - {\dot{q}}_{C k} ξ_{C 0}) + (p_{C k} {}^{C}D D_{γ}^{α, β} q_{C k} - H_{C}) {\dot{ξ}}_{C 0} + (p_{C k} \frac{d}{d t} {}^{C}D D_{γ}^{α, β} q_{C k} - \frac{\partial H_{C}}{\partial t}) ξ_{C 0} - \frac{\partial H_{C}}{\partial q_{C k}} ξ_{C k} + λ_{C a} \frac{\partial ϕ_{C a}}{\partial p_{C k}} η_{C k} \\ + {\dot{q}}_{C k} (t_{2}) ξ_{C 0} (t_{2}) \cdot \frac{(1 - γ) p_{C k}}{Γ (1 - β)} {(t_{2} - t)}^{- β} - {\dot{q}}_{C k} (t_{1}) ξ_{C 0} (t_{1}) \cdot \frac{γ \cdot p_{C k}}{Γ (1 - α)} {(t - t_{1})}^{- α} = 0 \end{matrix}$ Mathematical equation (50)

Equation (50) is called fractional Noether identity within a combined Caputo fractional derivative.

Theorem 2 If $ξ_{C 0}$ Mathematical equation , $ξ_{C k}$ , and $η_{C k}$ satisfy Eq. (50), then there exists a conserved quantity

$\begin{array}{l} C_{C} = (p_{C k} {}^{C}D D_{γ}^{α, β} q_{C k} - H_{C}) ξ_{C 0} + \int_{t_{1}}^{t} [p_{C k} {}^{C}D D_{γ}^{α, β} (ξ_{C k} - q_{C k} ξ_{C 0}) + (ξ_{C k} - {\dot{q}}_{C k} ξ_{C 0}) \cdot {}^{R L}D D_{1 - γ}^{β, α} p_{C k}] d τ \\ + {\dot{q}}_{C k} (t_{2}) ξ_{C 0} (t_{2}) \frac{1 - γ}{Γ (1 - β)} \int_{t_{1}}^{t} p_{C k} {(t_{2} - τ)}^{- β} d τ - {\dot{q}}_{C k} (t_{1}) ξ_{C 0} (t_{1}) \frac{γ}{Γ (1 - α)} \int_{t_{1}}^{t} p_{C k} {(τ - t_{1})}^{- α} d τ \end{array}$ Mathematical equation (51)

for the system within a combined Caputo fractional derivative.

Proof Using Eqs. (30), (37), and (50), it is easy to obtain

$\begin{array}{l} \frac{d C_{C}}{d t} = (p_{C k} {}^{C}D D_{γ}^{α, β} q_{C k} - H_{C}) {\dot{ξ}}_{C 0} + ξ_{C 0} ({\dot{p}}_{C k} {}^{C}D D_{γ}^{α, β} q_{C k} + p_{C k} \frac{d}{d t} {}^{C}D D_{γ}^{α, β} q_{C k} - \frac{\partial H_{C}}{\partial t} - \frac{\partial H_{C}}{\partial q_{C k}} \cdot {\dot{q}}_{C k} - \frac{\partial H_{C}}{\partial p_{C k}} \cdot {\dot{p}}_{C k}) \\ + p_{C k} \cdot {}^{C}D_{γ}^{α, β} (ξ_{C k} - {\dot{q}}_{C k} ξ_{C 0}) + (ξ_{C k} - {\dot{q}}_{C k} ξ_{C 0}) {}^{R L}D_{1 - γ}^{β, α} p_{C k} + {\dot{q}}_{C k} (t_{2}) ξ_{C 0} (t_{2}) \frac{(1 - γ) p_{C k}}{Γ (1 - β)} {(t_{2} - t)}^{- β} \\ - {\dot{q}}_{C k} (t_{1}) ξ_{C 0} (t_{1}) \frac{γ \cdot p_{C k}}{Γ (1 - α)} {(t - t_{1})}^{- α} + {\dot{q}}_{C k} (t_{2}) ξ_{C 0} (t_{2}) \frac{(1 - γ) p_{C k}}{Γ (1 - β)} {(t_{2} - t)}^{- β} - {\dot{q}}_{C k} (t_{1}) ξ_{C 0} (t_{1}) \frac{γ \cdot p_{C k}}{Γ (1 - α)} {(t - t_{1})}^{- α} \\ = \frac{\partial H_{C}}{\partial q_{C k}} ξ_{C k} - λ_{C a} \frac{\partial ϕ_{C a}}{\partial p_{C k}} η_{C k} + ξ_{C 0} ({\dot{p}}_{C k} \cdot λ_{C a} \frac{\partial ϕ_{C a}}{\partial p_{C k}} - \frac{\partial H_{C}}{\partial q_{C k}} {\dot{q}}_{C k}) - (ξ_{C k} - q_{C k} ξ_{C 0}) (\frac{\partial H_{C}}{\partial q_{C k}} + λ_{C a} \frac{\partial ϕ_{C a}}{\partial q_{C k}}) \\ = - λ_{C a} \frac{\partial ϕ_{C a}}{\partial p_{C k}} (η_{C k} - ξ_{C 0} {\dot{p}}_{C k}) - λ_{C a} \frac{\partial ϕ_{C a}}{\partial q_{C k}} (ξ_{C k} - {\dot{q}}_{C k} ξ_{C 0}) = - λ_{C a} \cdot δ ϕ_{C a} = 0 \end{array}$ Mathematical equation

The proof is completed.

6 Lie Symmetry Conserved Quantity

6.1 Lie Symmetry Conserved Quantity within the Combined Riemann-Liouville Fractional Derivative

Eq. (36) can be expressed as

$^{R L} D_{γ}^{α, β} q_{R L k} = h_{R L k} (t, q_{R L}, p_{R L}), k = 1,2, \dots, n,$ Mathematical equation (52)

${}^{C}_{1 - γ}^{β, α} p_{R L k} = f_{R L k} (t, q_{R L}, p_{R L}), k = 1,2, \dots, n .$ Mathematical equation (53)

We then study Eqs. (52) and (53) under the infinitesimal transformations (Eq. (43)). For Eq. (52), we obtain

$\begin{array}{l} ^{R L} D_{\bar{γ}}^{α, β} {\bar{q}}_{R L k} - h_{R L k} (\bar{t}, {\bar{q}}_{R L k}, {\bar{p}}_{R L k}) = {}^{R L}D_{γ}^{α, β} q_{R L k} - h_{R L k} (t, q_{R L}, p_{R L}) + θ_{R L} [^{R L} D_{γ}^{α, β} (ξ_{R L k} - {\dot{q}}_{R L k} ξ_{R L 0}) + ξ_{R L 0} \frac{d}{d t} {}^{R L}D D_{γ}^{α, β} q_{R L k} - X_{R L}^{(0)} (h_{R L k}) \\ - q_{R L k} (t_{1}) \frac{γ ξ_{R L 0} (t_{1})}{Γ (1 - α)} \frac{d}{d t} {(t - t_{1})}^{- α} + q_{R L k} (t_{2}) ξ_{R L 0} (t_{2}) \frac{1 - γ}{Γ (1 - β)} \frac{d}{d t} {(t_{2} - t)}^{- β}] \end{array}$ Mathematical equation (54)

where $X_{R L}^{(0)} = ξ_{R L 0} \frac{\partial}{\partial t} + ξ_{R L i} \frac{\partial}{\partial q_{R L i}} + η_{R L i} \frac{\partial}{\partial p_{R L i}}, i = 1,2, \dots, n . F o r E q . (53), w e h a v e$ Mathematical equation

$\begin{array}{l} {}^{C}_{1 - γ}^{β, α} {\bar{p}}_{R L k} = f_{R L k} (\bar{t}, {\bar{q}}_{R L}, {\bar{p}}_{R L}) = {}^{C}D_{1 - γ}^{β, α} p_{R L k} - f_{R L k} (t, q_{R L}, p_{R L}) + θ_{R L} [^{C} D_{1 - γ}^{β, α} (η_{R L k} - {\dot{p}}_{R L k} ξ_{R L 0}) + ξ_{R L 0} \frac{d}{d t} {}^{C}D D_{1 - γ}^{β, α} p_{R L k} - X_{R L}^{(0)} (f_{R L k}) \\ - \frac{1 - γ}{Γ (1 - β)} {(t - t_{1})}^{- β} {\dot{p}}_{R L k} (t_{1}) ξ_{R L 0} (t_{1}) + \frac{γ}{Γ (1 - α)} {(t_{2} - t)}^{- α} {\dot{p}}_{R L k} (t_{2}) ξ_{R L 0} (t_{2})] \end{array}$ Mathematical equation (55)

The fractional primary constraint (Eq. (27)) gives

$ϕ_{R L a} (t, {\bar{q}}_{R L}, {\bar{p}}_{R L}) = ϕ_{R L a} (t, q_{R L}, p_{R L}) + θ_{R L} X_{R L}^{(0)} (ϕ_{R L a}) .$ Mathematical equation (56)

Therefore, we have

$\begin{array}{l} ^{R L} D_{γ}^{α, β} (ξ_{R L k} - {\dot{q}}_{R L k} ξ_{R L 0}) + ξ_{R L 0} \frac{d}{d t} {}^{R L}D_{γ}^{α, β} q_{R L k} - X_{R L}^{(0)} (h_{R L k}) \\ - q_{R L k} (t_{1}) \frac{γ ξ_{R L 0} (t_{1})}{Γ (1 - α)} \frac{d}{d t} {(t - t_{1})}^{- α} + q_{R L k} (t_{2}) ξ_{R L 0} (t_{2}) \frac{1 - γ}{Γ (1 - β)} \frac{d}{d t} {(t_{2} - t)}^{- β} = 0, \end{array}$ Mathematical equation (57)

$\begin{array}{l} {}^{C}_{1 - γ}^{β, α} (η_{R L k} - {\dot{p}}_{R L k} ξ_{R L 0}) + ξ_{R L 0} \frac{d}{d t} {}^{C}D_{1 - γ}^{β, α} p_{R L k} - X_{R L}^{(0)} (f_{R L k}) \\ - \frac{1 - γ}{Γ (1 - β)} {(t - t_{1})}^{- β} {\dot{p}}_{R L k} (t_{1}) ξ_{R L 0} (t_{1}) + \frac{γ}{Γ (1 - α)} {(t_{2} - t)}^{- α} {\dot{p}}_{R L k} (t_{2}) ξ_{R L 0} (t_{2}) = 0, \end{array}$ Mathematical equation (58)

and

$X_{R L}^{(0)} (ϕ_{R L a}) = 0 .$ Mathematical equation (59)

Equations (57) and (58) are called the determined equations within the combined Riemann-Liouville fractional derivative, and Eq. (59) is called the limited equation within the combined Riemann-Liouville fractional derivative.

Besides, if we consider the deduction process of the system (Eq. (36)), an extra additional limited equation

$\frac{\partial ϕ_{R L a}}{\partial q_{R L i}} (ξ_{R L i} - {\dot{q}}_{R L i} ξ_{R L 0}) + \frac{\partial ϕ_{R L a}}{\partial p_{R L i}} (η_{R L i} - {\dot{p}}_{R L i} ξ_{R L 0}) = 0$ Mathematical equation (60)

can be obtained. Then we have

Theorem 3 For the system within the combined Riemann-Liouville fractional derivative (Eq. (36)), when Eqs. (45), (57) and (58) are satisfied, Eq. (46) gives a Lie symmetry conserved quantity; when Eqs. (45), (57)-(59) are satisfied, Eq. (46) gives a weak Lie symmetry conserved quantity; when Eqs. (45), (57)-(60) are satisfied, Eq. (46) gives a strong Lie symmetry conserved quantity.

6.2 Lie Symmetry Conserved Quantity within the Combined Caputo Fractional Derivative

Eq. (37) can be expressed as

${}^{C}_{γ}^{α, β} q_{C k} = h_{C k} (t, q_{C}, p_{C}), k = 1,2, \dots, n,$ Mathematical equation (61)

${}^{R L}_{1 - γ}^{β, α} p_{C k} = f_{C k} (t, q_{C}, p_{C}), k = 1,2, \dots, n .$ Mathematical equation (62)

Using the similar way, we can obtain the determined equations within the combined Caputo fractional derivative

$^{C} D_{γ}^{α, β} (ξ_{C k} - {\dot{q}}_{C k} ξ_{C 0}) + ξ_{C 0} \frac{d}{d t} {}^{C}D_{γ}^{α, β} q_{C k} - X_{C}^{(0)} (h_{C k}) - {\dot{q}}_{C k} (t_{1}) \frac{γ ξ_{C 0} (t_{1})}{Γ (1 - α)} {(t - t_{1})}^{- α} + {\dot{q}}_{C k} (t_{2}) ξ_{C 0} (t_{2}) \frac{1 - γ}{Γ (1 - β)} {(t_{2} - t)}^{- β} = 0,$ Mathematical equation (63)

$\begin{array}{l} {}^{R L}_{1 - γ}^{β, α} (η_{C k} - {\dot{p}}_{C k} ξ_{C 0}) + ξ_{C 0} \frac{d}{d t} {}^{R L}D_{1 - γ}^{β, α} p_{C k} - X_{C}^{(0)} - \frac{1 - γ}{Γ (1 - β)} p_{C k} (t_{1}) ξ_{C 0} (t_{1}) \frac{d}{d t} {(t - t_{1})}^{- β} \\ + \frac{γ}{Γ (1 - α)} p_{C k} (t_{2}) ξ_{C 0} (t_{2}) \frac{d}{d t} {(t_{2} - t)}^{- α} = 0, \end{array}$ Mathematical equation (64)

the limited equation within the combined Caputo fractional derivative

$X_{C}^{(0)} (ϕ_{C a}) = 0,$ Mathematical equation (65)

and the additional limited equation within the combined Caputo fractional derivative

$\frac{\partial ϕ_{C a}}{\partial q_{C i}} (ξ_{C i} - {\dot{q}}_{C i} ξ_{C 0}) + \frac{\partial ϕ_{C a}}{\partial p_{C i}} (η_{C i} - {\dot{p}}_{C i} ξ_{C 0}) = 0,$ Mathematical equation (66)

where $X_{C}^{(0)} = ξ_{C 0} \frac{\partial}{\partial t} + ξ_{C i} \frac{\partial}{\partial q_{C i}} + η_{C i} \frac{\partial}{\partial p_{C i}}, i = 1,2, \dots, n$ Mathematical equation . Then we have

Theorem 4 For the system within the combined Caputo fractional derivative (Eq. (37)), when Eqs. (50), (63) and (64) are satisfied, Eq. (51) gives a Lie symmetry conserved quantity; when Eqs. (50), (63)-(65) are satisfied, Eq. (51) gives a weak Lie symmetry conserved quantity; when Eqs. (50), (63)-(66) are satisfied, Eq. (51) gives a strong Lie symmetry conserved quantity.

7 An Example

The singular Lagrangian is

$L_{R L} = - \frac{c}{2} {(q_{R L 1})}^{2} + \frac{b}{2} {(q_{R L 2})}^{2} + \frac{1}{2} ({q_{R L 2}}^{R L} D_{γ}^{α, β} q_{R L 1} - {q_{R L 1}}^{R L} D_{γ}^{α, β} q_{R L 2}) .$ Mathematical equation (67)

Try to study its Noether symmetry and conserved quantity.

From Eqs. (19) and (20), we can obtain

$p_{R L 1} = \frac{\partial L_{R L}}{\partial^{R L} D_{γ}^{α, β} q_{R L 1}} = \frac{1}{2} q_{R L 2}, p_{R L 2} = \frac{\partial L_{R L}}{\partial^{R L} D_{γ}^{α, β} q_{R L 2}} = - \frac{1}{2} q_{R L 1},$ Mathematical equation (68)

$H_{R L} = {p_{R L 1}}^{R L} D_{γ}^{α, β} q_{R L 1} + {p_{R L 2}}^{R L} D_{γ}^{α, β} q_{R L 2} - L_{R L} = \frac{c}{2} {(q_{R L 1})}^{2} - \frac{b}{2} {(q_{R L 2})}^{2} .$ Mathematical equation (69)

Then Eq. (27) gives two fractional primary constraints

$ϕ_{R L 1} = p_{R L 1} - \frac{1}{2} q_{R L 2} = 0, ϕ_{R L 2} = p_{R L 2} + \frac{1}{2} q_{R L 1} = 0 .$ Mathematical equation (70)

From Eq. (39), we obtain

$b {\dot{q}}_{R L 1} q_{R L 2} + λ_{R L 1} {\dot{q}}_{R L 1} + {\dot{p}}_{R L 2}^{R L} D_{γ}^{α, β} q_{R L 1} - {\dot{q}}_{R L 1} [{}^{C} D_{1 - γ}^{β, α} p_{R L 2} - \frac{{(t_{2} - t)}^{- α}}{Γ (1 - α)} p_{R L 2} (t_{2}) + \frac{(1 - γ) {(t - t_{1})}^{- β}}{Γ (1 - β)} p_{R L 2} (t_{1})] = 0,$ Mathematical equation (71)

$c q_{R L 1} {\dot{q}}_{R L 2} + λ_{R L 2} {\dot{q}}_{R L 2} - {\dot{p}}_{R L 1}^{R L} D_{γ}^{α, β} q_{R L 2} + {\dot{q}}_{R L 2} [{}^{C} D_{1 - γ}^{β, α} p_{R L 1} - \frac{γ {(t_{2} - t)}^{- α}}{Γ (1 - α)} p_{R L 1} (t_{2}) + \frac{(1 - γ) {(t - t_{1})}^{- β}}{Γ (1 - β)} p_{R L 1} (t_{1})] = 0 .$ Mathematical equation (72)

Then, using Eq. (36), the fractional constrained Hamilton equation within the combined Riemann-Liouville fractional derivative is

$^{R L} D_{γ}^{α, β} q_{R L 1} = - 2 b q_{R L 2} + 2 [{}^{C} D_{1 - γ}^{β, α} p_{R L 2} - \frac{γ {(t_{2} - t)}^{- α}}{Γ (1 - α)} p_{R L 2} (t_{2}) + \frac{(1 - γ) {(t - t_{1})}^{- β}}{Γ (1 - β)} p_{R L 2} (t_{1})],$ Mathematical equation

$^{R L} D_{γ}^{α, β} q_{R L 2} = - 2 c q_{R L 1} - 2 [{}^{C} D_{1 - γ}^{β, α} p_{R L 1} - \frac{γ {(t_{2} - t)}^{- α}}{Γ (1 - α)} p_{R L 1} (t_{2}) + \frac{(1 - γ) {(t - t_{1})}^{- β}}{Γ (1 - β)} p_{R L 1} (t_{1})] .$ Mathematical equation (73)

From Noether identity (Eq. (45)), we have

$\begin{array}{l} p_{R L 1} {}^{R L}D_{γ}^{α, β} (ξ_{R L 1} - {\dot{q}}_{R L 1} ξ_{R L 0}) + {p_{R L 2}}^{R L} D_{γ}^{α, β} (ξ_{R L 2} - {\dot{q}}_{R L 2} ξ_{R L 0}) - c q_{R L 1} ξ_{R L 1} + b q_{R L 2} ξ_{R L 2} \\ + (p_{R L 1} \frac{d}{d t} {}^{R L}D D_{γ}^{α, β} q_{R L 1} + p_{R L 2} \frac{d}{d t} {}^{R L}D D_{γ}^{α, β} q_{R L 2}) ξ_{R L 0} + λ_{R L 1} η_{R L 1} + λ_{R L 2} η_{R L 2} + ({p_{R L 1}}^{R L} D_{γ}^{α, β} q_{R L 1} + {p_{R L 2}}^{R L} D_{γ}^{α, β} q_{R L 2} - H_{R L}) {\dot{ξ}}_{R L 0} \\ + q_{R L 1} (t_{2}) ξ_{R L 0} (t_{2}) \times \frac{(1 - γ) p_{R L 1}}{Γ (1 - β)} \frac{d}{d t} {(t_{2} - t)}^{- β} - q_{R L 1} (t_{1}) ξ_{R L 0} (t_{1}) \frac{γ p_{R L 1}}{Γ (1 - α)} \frac{d}{d t} {(t - t_{1})}^{- α} \\ + q_{R L 2} (t_{2}) ξ_{R L 0} (t_{2}) \times \frac{(1 - γ) p_{R L 2}}{Γ (1 - β)} \frac{d}{d t} {(t_{2} - t)}^{- β} - q_{R L 2} (t_{1}) ξ_{R L 0} (t_{1}) \frac{γ p_{R L 2}}{Γ (1 - α)} \frac{d}{d t} {(t - t_{1})}^{- α} = 0 \end{array}$ Mathematical equation (74)

Then we can verify that

$ξ_{R L 0} = - 1, ξ_{R L 1} = ξ_{R L 2} = 0, η_{R L 1} = η_{R L 2} = 0$ Mathematical equation (75)

satisfy Eq. (74). Therefore, we have

$\begin{array}{l} C_{R L} = - ({p_{R L 1}}^{R L} D_{γ}^{α, β} q_{R L 1} + {p_{R L 2}}^{R L} D_{γ}^{α, β} q_{R L 2} - \frac{c}{2} q_{R L 1}^{2} + \frac{b}{2} q_{R L 2}^{2}) + \int_{t_{1}}^{t} {p_{R L 1} \frac{d}{d τ} {}^{R L}D D_{γ}^{α, β} q_{R L 1} + p_{R L 2} \frac{d}{d τ} {}^{R L}D D_{γ}^{α, β} q_{R L 2} \\ + {\dot{q}}_{R L 1} [{}^{C} D_{1 - γ}^{β, α} p_{R L 1} - \frac{γ {(t_{2} - τ)}^{- α}}{Γ (1 - α)} p_{R L 1} (t_{2}) + \frac{(1 - γ) {(τ - t_{1})}^{- β}}{Γ (1 - β)} p_{R L 1} (t_{1})] \\ + {\dot{q}}_{R L 2} [{}^{C} D_{1 - γ}^{β, α} p_{R L 2} - \frac{γ {(t_{2} - τ)}^{- α}}{Γ (1 - α)} p_{R L 2} (t_{2}) + \frac{(1 - γ) {(τ - t_{1})}^{- β}}{Γ (1 - β)} p_{R L 2} (t_{1})]} d τ \end{array}$ Mathematical equation (76)

Specially, we degenerate the above fractional order to the integer order, namely $α, β \to 1$ Mathematical equation , so we can get

$C_{R L 0} = - \frac{c}{2} q_{R L 1}^{2} + \frac{b}{2} q_{R L 2}^{2} .$ Mathematical equation (77)

If $t \in [1,10]$ Mathematical equation , and let $b = c = 1, q_{R L 1} (1) = 0, q_{R L 2} (1) = 1,$ the trajectory of $C_{R L 0}$ can be drawn as follows.

It follows from Fig. 1 that $C_{R L 0}$ Mathematical equation is a constant.

Thumbnail: Fig. 1 Refer to the following caption and surrounding text.

Fig. 1 The trajectory of

C_{R L 0}

8 Conclusion

Fractional calculus has important applications in various fields, and it is also a research hotspot. In this paper, based on the combined Riemann-Liouville fractional derivative and combined Caputo fractional derivative, we establish Noether theorem and Lie theorem for the fractional constrained Hamiltonian system. Besides, Noether identity, determined equation, limited equation and additional limited equation are presented.

In this paper, Noether-type conserved quantity is obtained from Lie symmetry, whether the Lie symmetry of the system can lead to another conserved quantity, called Hojman conserved quantity, remains to be studied. At the same time, Mei symmetry is also an effective method to solve the differential equation. Therefore, Mei symmetry is also a future research direction.

References

Oldham K B, Spanier J. The Fractional Calculus [M]. New York: Academic Press, 1974. [MathSciNet] [Google Scholar]
Luo S K. Analytical mechanics method of fractional dynamics and its applications [J]. Journal of Dynamics and Control, 2019, 17(5): 432-438. [Google Scholar]
Yin Y D, Zhao D M, Liu J L, et al. Study on the rate dependency of acrylic elastomer based fractional viscoelastic model [J]. Chinese Journal of Theoretical and Applied Mechanics, 2022, 54(1): 155-163. [Google Scholar]
Si H, Zheng Y A. Adaptive predictive synchronization of fractional order chaotic systems [J]. Journal of Dynamics and Control, 2021, 19(5): 8-12(Ch). [Google Scholar]
Song C J. Generalized fractional forced Birkhoff equations [J]. Journal of Dynamics and Control, 2019, 17(5): 446-452. [Google Scholar]
Riewe F. Nonconservative Lagrangian and Hamiltonian mechanics [J]. Phys Rev E, 1996, 53(2): 1890-1899. [CrossRef] [MathSciNet] [PubMed] [Google Scholar]
Riewe F. Mechanics with fractional derivatives [J]. Phys Rev E, 1997, 55(3): 3581-3592. [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
Malinowska A B, Torres D F M. Fractional calculus of variations for a combined Caputo derivative [J]. Fract Calc Appl Anal, 2011, 14(4): 523-537. [Google Scholar]
Zhang Y. Fractional differential equations of motion in terms of combined Riemann-Liouville derivatives [J]. Chin Phys B, 2012, 21(8): 84502. [Google Scholar]
Luo S K, Xu Y L. Fractional Birkhoffian mechanics [J]. Acta Mech, 2015, 226(3): 829-844. [CrossRef] [MathSciNet] [Google Scholar]
Li Z P. Classical and Quantal Dynamics of Constrained Systems and Their Symmetrical Properties [M]. Beijing: Beijing Polytechnic University Press, 1993(Ch). [Google Scholar]
Dirac P A M. Lecture on Quantum Mechanics [M]. New York: Yeshiva University, 1964. [Google Scholar]
Li Z P. Contrained Hamiltonian Systems and Their Symmetrical Properties [M]. Beijing: Beijing Polytechnic University Press, 1999. [Google Scholar]
Li Z P, Jiang J H. Symmetries in Constrained Canonical Systems [M]. Beijing: Science Press, 2002. [Google Scholar]
Noether A E. Invariante variationsprobleme [J]. Nachrichten von der Gesellschaft der Gesellschaft der Wissenschaften zu Göttingen, Mathematisch-Physikalische Klasse, 1918, KI(2): 235-257. [Google Scholar]
Lutzky M. Dynamical symmetries and conserved quantities [J]. Journal of Physics A: Mathematical and General, 1979, 12(7): 973-981. [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
Mei F X, Wu H B. Dynamics of Constrained Mechanical Systems [M]. Beijing: Beijing Institute of Technology Press, 2009. [Google Scholar]
Mei F X. Analytical Mechanics (II) [M]. Beijing: Beijing Institute of Technology Press, 2013. [Google Scholar]
Song C J, Wang J H. Conserved quantity for fractional constrained Hamiltonian system [J]. Wuhan University Journal of Natural Sciences, 2022, 27(3): 201-210. [CrossRef] [EDP Sciences] [Google Scholar]
Song C J, Zhai X H. Noether theorem for fractional singular systems [J]. Wuhan University Journal of Natural Sciences, 2023, 28(3): 207-216. [CrossRef] [EDP Sciences] [Google Scholar]
Song C J, Zhang Y. Noether symmetry and conserved quantity for fractional Birkhoffian mechanics and its applications [J]. Frac Cal Appl Anal, 2018, 21(2): 509-526. [CrossRef] [MathSciNet] [Google Scholar]
Podlubny I. Fractional Differential Equations [M]. San Diego: Academic Press, 1999. [Google Scholar]
Wu Q, Huang J H. Fractional Order Calculus [M]. Beijing: Tsinghua University Press, 2016(Ch). [Google Scholar]
Gelfand I M, Fomin S V. Calculus of Variations [M]. London: Prentice-Hall, 1963. [Google Scholar]

All Figures

	Fig. 1 The trajectory of $C_{R L 0}$
In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.

[R1] Oldham K B, Spanier J. The Fractional Calculus [M]. New York: Academic Press, 1974. [MathSciNet] [Google Scholar]

[R2] Luo S K. Analytical mechanics method of fractional dynamics and its applications [J]. Journal of Dynamics and Control, 2019, 17(5): 432-438. [Google Scholar]

[R3] Yin Y D, Zhao D M, Liu J L, et al. Study on the rate dependency of acrylic elastomer based fractional viscoelastic model [J]. Chinese Journal of Theoretical and Applied Mechanics, 2022, 54(1): 155-163. [Google Scholar]

[R4] Si H, Zheng Y A. Adaptive predictive synchronization of fractional order chaotic systems [J]. Journal of Dynamics and Control, 2021, 19(5): 8-12(Ch). [Google Scholar]

[R5] Song C J. Generalized fractional forced Birkhoff equations [J]. Journal of Dynamics and Control, 2019, 17(5): 446-452. [Google Scholar]

[R6] Riewe F. Nonconservative Lagrangian and Hamiltonian mechanics [J]. Phys Rev E, 1996, 53(2): 1890-1899. [CrossRef] [MathSciNet] [PubMed] [Google Scholar]

[R7] Riewe F. Mechanics with fractional derivatives [J]. Phys Rev E, 1997, 55(3): 3581-3592. [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]

[R8] Malinowska A B, Torres D F M. Fractional calculus of variations for a combined Caputo derivative [J]. Fract Calc Appl Anal, 2011, 14(4): 523-537. [Google Scholar]

[R9] Zhang Y. Fractional differential equations of motion in terms of combined Riemann-Liouville derivatives [J]. Chin Phys B, 2012, 21(8): 84502. [Google Scholar]

[R10] Luo S K, Xu Y L. Fractional Birkhoffian mechanics [J]. Acta Mech, 2015, 226(3): 829-844. [CrossRef] [MathSciNet] [Google Scholar]

[R11] Li Z P. Classical and Quantal Dynamics of Constrained Systems and Their Symmetrical Properties [M]. Beijing: Beijing Polytechnic University Press, 1993(Ch). [Google Scholar]

[R12] Dirac P A M. Lecture on Quantum Mechanics [M]. New York: Yeshiva University, 1964. [Google Scholar]

[R13] Li Z P. Contrained Hamiltonian Systems and Their Symmetrical Properties [M]. Beijing: Beijing Polytechnic University Press, 1999. [Google Scholar]

[R14] Li Z P, Jiang J H. Symmetries in Constrained Canonical Systems [M]. Beijing: Science Press, 2002. [Google Scholar]

[R15] Noether A E. Invariante variationsprobleme [J]. Nachrichten von der Gesellschaft der Gesellschaft der Wissenschaften zu Göttingen, Mathematisch-Physikalische Klasse, 1918, KI(2): 235-257. [Google Scholar]

[R16] Lutzky M. Dynamical symmetries and conserved quantities [J]. Journal of Physics A: Mathematical and General, 1979, 12(7): 973-981. [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]

[R17] Mei F X, Wu H B. Dynamics of Constrained Mechanical Systems [M]. Beijing: Beijing Institute of Technology Press, 2009. [Google Scholar]

[R18] Mei F X. Analytical Mechanics (II) [M]. Beijing: Beijing Institute of Technology Press, 2013. [Google Scholar]

[R19] Song C J, Wang J H. Conserved quantity for fractional constrained Hamiltonian system [J]. Wuhan University Journal of Natural Sciences, 2022, 27(3): 201-210. [CrossRef] [EDP Sciences] [Google Scholar]

[R20] Song C J, Zhai X H. Noether theorem for fractional singular systems [J]. Wuhan University Journal of Natural Sciences, 2023, 28(3): 207-216. [CrossRef] [EDP Sciences] [Google Scholar]

[R21] Song C J, Zhang Y. Noether symmetry and conserved quantity for fractional Birkhoffian mechanics and its applications [J]. Frac Cal Appl Anal, 2018, 21(2): 509-526. [CrossRef] [MathSciNet] [Google Scholar]

[R22] Podlubny I. Fractional Differential Equations [M]. San Diego: Academic Press, 1999. [Google Scholar]

[R23] Wu Q, Huang J H. Fractional Order Calculus [M]. Beijing: Tsinghua University Press, 2016(Ch). [Google Scholar]

[R24] Gelfand I M, Fomin S V. Calculus of Variations [M]. London: Prentice-Hall, 1963. [Google Scholar]