Conserved Quantity and Adiabatic Invariant for Hamiltonian System with Variable Order

Abstract

Hamiltonian mechanics plays an important role in the development of nonlinear science. This paper aims for a fractional Hamiltonian system of variable order. Several issues are discussed, including differential equation of motion, Noether symmetry, and perturbation to Noether symmetry. As a result, fractional Hamiltonian mechanics of variable order are established, and conserved quantity and adiabatic invariant are presented. Two applications, fractional isotropic harmonic oscillator model of variable order and fractional Lotka biochemical oscillator model of variable order are given to illustrate the Methods and Results.

1. Introduction

Fractional calculus was firstly proposed because the order of a function’s integral or derivative is a non-integer constant. In fact, fractional calculus started more than 300 years ago when L’Hopital and Leibniz were discussing the meaning of ${\mathrm{d}}^{1/2}y/\mathrm{d}{x}^{1/2}$. The development of fractional calculus is much slower than that of integer calculus though many famous mathematicians, such as Fourier, Euler, Riemann, Liouville, Letnikov, Grunwald, etc., have contributed to it (see [1,2,3] for the history of fractional calculus). In recent decades, because of fractional calculus’s many applications in various fields of science, engineering, biomechanics, economics, and so forth, it has gained more attention (see [4,5,6,7] for a review). For example, fractional calculus has been applied to physics [8], quantum mechanics [9,10,11], field theory [12,13], etc. Among several definitions of the fractional derivatives, which are generally nonlocal operators and are historically applied to study time-dependent or nonlocal processes, Riemann–Liouville fractional derivative and Caputo fractional derivative are the most famous ones.

Furthermore, Samko and Ross [14] introduced a generalization of fractional calculus. They considered the order of a function’s integral or derivative as $\alpha \left(·,·\right)$, where $\alpha \left(·,·\right)$ is a function rather than a constant. Because of the nonlocality and memory of fractional order calculus, it is reasonable that the order of a function’s integral or derivative may vary. Afterwards, several works were dedicated to the properties and applications of the fractional operators with variable order [15,16,17,18,19,20]. For instance, Coimbra [17] examined dynamical behaviors with frictional force using the variable order model. Diaz and Coimbra [18] utilized the differential equation with variable order to reveal some dynamical control behaviors of the nonlinear oscillator. Sun et al. [20] investigated a class of fractional models with variable order. Nowadays, more and more fractional models with variable order have been applied to mechanics, mathematics, and other related subjects (see [21,22,23,24] and references therein).

In 1996, Riewe started the theory of the fractional calculus of variations by considering that the first-order derivative terms should come from the half-order derivative terms [25,26]. From then on, scholars began to use fractional derivatives to describe the dissipative forces such as frictional forces for the nonconservative systems. Then, fractional calculus of variations was studied by several scholars such as Klimek [27], Agrawal [28,29], Baleanu et al. [30,31], Torres [32], etc. Moreover, the development of the fractional calculus of variations with variable order also made great progress. In [33,34,35], several results for the fractional calculus of variations with variable order were obtained. Particularly, motivated by references [36,37,38], a linear combination of the fractional derivative of variable order was introduced. Based on a combined Caputo fractional derivative of variable order (CCVO), the necessary optimality conditions for variational problems were established [39,40], the fractional variational problem of Herglotz type with variable order was studied [41], and two fractional isoperimetric problems and a new variational problem subject to a holonomic constraint were presented [42].

After differential equations of motion are established through the calculus of variations, the next task is to find the solutions to them in dynamics. In fact, the solution to the equations can be given if one can find all of the integrals of the equations. An integral is a conserved quantity; therefore, people try their best to find all of the conserved quantities of a mechanical system. By means of the analysis of forces, Newtonian mechanics gives three conservation laws, i.e., the conservation of momentum, the conservation of mechanical energy, and the conservation of moment of momentum. By means of the analysis of the form of Lagrangian, Lagrangian mechanics gives two conservation laws, i.e., the conservation of generalized energy and the conservation of generalized momentum. The conservation of generalized momentum may be a conservation of momentum, a conservation of moment of momentum, or neither. The physical meaning of the conservation law in Lagrangian mechanics is less clear than that in Newtonian mechanics, whose three conservation laws have very clear physical meaning. However, the conserved quantities deduced by Lagrangian mechanics are more than those deduced by Newtonian mechanics. Since German mathematician Emmy Noether published her famous paper [43], the Noether symmetry method has become a modern method for seeking the conservation law of mechanical systems (see [44] for a review). Similarly, although the physical meaning of the conservation law obtained from the Noether symmetry method is less than that obtained from Lagrangian mechanics, the numbers of the conserved quantities deduced by the Noether symmetry method are more than those deduced by Lagrangian mechanics.

Recently, fractional Noether symmetry and fractional conserved quantity are also under strong research. There are two different definitions of fractional conserved quantity. One was introduced by Frederico and Torres [45] by means of a bilinear fractional operator $D (D(C)=0)$, and the other was introduced by Atanacković [46] by means of the classical definition $(\mathrm{d}C/\mathrm{d}t=0)$, where C means a fractional conserved quantity. Fractional Noether symmetry and fractional conserved quantity have been investigated through both definitions. For example, works [47,48] were done based on the former definition, and results [49,50,51,52,53] were obtained on the basis of the latter one. Furthermore, the study of the fractional Noether symmetry and fractional conserved quantity with variable order has also begun. For instance, Odzijewicz, Malinowska, and Torres [54,55] investigated Noether theorem and the second Noether theorem for the fractional Lagrangian system with variable order. Yan and Zhang [56,57] studied Noether symmetry and conserved quantity for the fractional Birkhoffian system in terms of Caputo fractional derivative of variable order (CVO), Riemann–Liouville fractional derivative of variable order (RLVO), Riesz–Caputo fractional derivative of variable order (RCVO), and Riesz–Riemann–Liouville fractional derivative of variable order (RVO).

When a dynamical system is disturbed by small forces, the conserved quantity, which refers to the integrability of the system, may also change. Therefore, the research on the perturbation to Noether symmetry and adiabatic invariants is also of great significance for a dynamical system [58,59]. A classical adiabatic invariant means a certain physical quantity that changes more slowly than the parameter that varies very slowly. That is, when a dynamical system is disturbed by small forces and varies very slowly, a slower physical quantity, i.e., adiabatic invariant, will be found. Some important results on the study of perturbation to Noether symmetry and adiabatic invariants for constrained mechanical systems have been obtained [60,61]. Very recently, we studied the perturbation to fractional Noether symmetry for the Hamiltonian system [53] as well as the Birkhoffian system [62]. Besides, the results on perturbation to Noether symmetry with RLVO for the fractional generalized Birkhoffian system [63] were also presented.

The purpose of this paper is to generalize the problem of the calculus of variations, Noether theory, and perturbation to Noether symmetry to the fractional Hamiltonian systems in terms of combined Riemann–Liouville fractional derivatives of variable order (CRLVO) and combined Caputo fractional derivatives of variable order (CCVO), respectively. Here, the CRLVO and the CCVO are the first time to be introduced, and many results obtained before are the special cases of this paper. As the main results, Hamilton equations with the CRLVO and the CCVO are established. Then, Noether symmetry and conserved quantities for the fractional Hamiltonian systems with the CRLVO and the CCVO are presented. Perturbation to Noether symmetry and adiabatic invariants for the disturbed fractional Hamiltonian systems with the CRLVO and the CCVO are also investigated, and two applications are given to illustrate the Methods and Results of this paper.

In a word, there are mainly three factors for presenting this paper. Firstly, analytical mechanics. Lagrange was the first man to publish the famous book Analytical Mechanics. Lagrange felt proud when he successfully described mechanics problems mathematically without any diagrams because within the framework of Newtonian mechanics, mechanical phenomena are treated mostly schematically. In fact, the representation by means of analytical mechanics is more general as well as convenient. Secondly, Noether theory. Noether theory not only has mathematical importance but can also be helpful in understanding some inherent physical properties of a dynamical system. Generally, conserved quantity helps reduce the degrees of freedom and make it easier to find the solutions to the differential equations of motion. Thirdly, fractional calculus. Fractional calculus is a useful tool to deal with the dissipative systems, and, as early as 1993, Miller and Ross [64] once said that fractional calculus has been involved in almost every field of science and engineering.

This paper is constructed as follows: Combined fractional derivatives of variable order and their properties are given in Section 2. Differential equations of motion for the Hamiltonian system of variable order are established in Section 3. In Section 4 and Section 5, Noether symmetry and conserved quantity, and perturbation to Noether symmetry and adiabatic invariants are investigated. The Results and Methods are illustrated in Section 6 with two applications and we finish with Section 7 for Conclusions.

2. Combined Fractional Derivative of Variable Order

Following the references [54,57], we list the left and right Riemann–Liouville fractional integrals of variable order, the left and right RLVO, and the left and right CVO by

$${}_{t_1}I_t^{\alpha \left(·,·\right)}f\left(t\right)={\displaystyle {\int }_{t_1}^t\frac{1}{\mathsf{\Gamma }\left(\alpha \left(t,\tau \right)\right)}{\left(t-\tau \right)}^{\alpha \left(t,\tau \right)-1}f\left(\tau \right)\mathrm{d}\tau }$$

(1)

$${}_{t}I_{t_2}^{\alpha \left(·,·\right)}f\left(t\right)={\displaystyle {\int }_t^{t_2}\frac{1}{\mathsf{\Gamma }\left(\alpha \left(\tau ,t\right)\right)}{\left(\tau -t\right)}^{\alpha \left(\tau ,t\right)-1}f\left(\tau \right)\mathrm{d}\tau }$$

(2)

$${}_{t_1}{}^{RL}D_t^{\alpha \left(·,·\right)}f\left(t\right)=\frac{\mathrm{d}}{\mathrm{d}t}{}_{t_1}I_t^{1-\alpha \left(·,·\right)}f\left(t\right)=\frac{\mathrm{d}}{\mathrm{d}t}{\displaystyle {\int }_{t_1}^t\frac{1}{\mathsf{\Gamma }\left(1-\alpha \left(t,\tau \right)\right)}{\left(t-\tau \right)}^{-\alpha \left(t,\tau \right)}f\left(\tau \right)\mathrm{d}\tau }$$

(3)

$${}_t{}^{RL}D_{t_2}^{\alpha \left(·,·\right)}f\left(t\right)=\left(-\frac{\mathrm{d}}{\mathrm{d}t}\right){}_{t}I_{t_2}^{1-\alpha \left(·,·\right)}f\left(t\right)=\left(-\frac{\mathrm{d}}{\mathrm{d}t}\right){\displaystyle {\int }_t^{t_2}\frac{1}{\mathsf{\Gamma }\left(1-\alpha \left(\tau ,t\right)\right)}{\left(\tau -t\right)}^{-\alpha \left(\tau ,t\right)}f\left(\tau \right)\mathrm{d}\tau }$$

(4)

$${}_{t_1}{}^{C}D_t^{\alpha \left(·,·\right)}f\left(t\right)={\displaystyle {\int }_{t_1}^t\frac{1}{\mathsf{\Gamma }\left(1-\alpha \left(t,\tau \right)\right)}{\left(t-\tau \right)}^{-\alpha \left(t,\tau \right)}\frac{\mathrm{d}}{\mathrm{d}\tau }f\left(\tau \right)\mathrm{d}\tau }$$

(5)

$${}_t{}^{C}D_{t_2}^{\alpha \left(·,·\right)}f\left(t\right)={\displaystyle {\int }_t^{t_2}\frac{1}{\mathsf{\Gamma }\left(1-\alpha \left(\tau ,t\right)\right)}{\left(\tau -t\right)}^{-\alpha \left(\tau ,t\right)}\left(-\frac{\mathrm{d}}{\mathrm{d}\tau }\right)f\left(\tau \right)\mathrm{d}\tau }$$

(6)

where $0<\alpha \left(·,·\right)<1$, $\mathsf{\Gamma }(·)$ is the Gamma function, $f\left(t\right)\in L_1\left[t_1,t_2\right]$, ${}_{t_1}I_t^{1-\alpha \left(·,·\right)}f\left(t\right)\in AC\left[t_1,t_2\right]$, ${}_{t}I_{t_2}^{1-\alpha \left(·,·\right)}f\left(t\right)\in AC\left[t_1,t_2\right]$. In this paper, we assume that $\alpha \left(t,\tau \right)=\alpha \left(t-\tau \right)$, $\alpha \left(\tau ,t\right)=\alpha \left(\tau -t\right)$.

Based on the definitions of the RLVO and the CVO, the RVO and the RCVO are defined in [57] as

$${}_{t_1}{}^{R}D_{t_2}^{\alpha \left(·,·\right)}f\left(t\right)=\frac{1}{2}\left[{}_{t_1}{}^{RL}D_t^{\alpha \left(·,·\right)}f\left(t\right)-{}_t{}^{RL}D_{t_2}^{\alpha \left(·,·\right)}f\left(t\right)\right]$$

(7)

$${}_{t_1}{}^{RC}D_{t_2}^{\alpha \left(·,·\right)}f\left(t\right)=\frac{1}{2}\left[{}_{t_1}{}^{C}D_t^{\alpha \left(·,·\right)}f\left(t\right)-{}_t{}^{C}D_{t_2}^{\alpha \left(·,·\right)}f\left(t\right)\right]$$

(8)

In this paper, CRLVO and CCVO are constructed as

$${}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}f\left(t\right)=\gamma {}_{t_1}{}^{RL}D_t^{\alpha \left(·,·\right)}f\left(t\right)-\left(1-\gamma \right){}_t{}^{RL}D_{t_2}^{\beta \left(·,·\right)}f\left(t\right)$$

(9)

$${}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}f\left(t\right)=\gamma {}_{t_1}{}^{C}D_t^{\alpha \left(·,·\right)}f\left(t\right)-\left(1-\gamma \right){}_t{}^{C}D_{t_2}^{\beta \left(·,·\right)}f\left(t\right)$$

(10)

where $0<\alpha \left(·,·\right)<1$, $0<\beta \left(·,·\right)<1$, ${}_{t_1}D_t^{\alpha \left(·,·\right)}$ and ${}_{t}D_{t_2}^{\beta \left(·,·\right)}$ show the arrow of time, and $\gamma$ is a parameter, $0\le \gamma \le 1$, which determines the different quantities of information from the past and the future.

The CRLVO and the CCVO are the most general operators with variable order to some extent. From Equations (9) and (10), the RVO, the RCVO, the left and right Riemann–Liouville fractional derivatives (RL), the left and right Caputo fractional derivatives (C), the Riesz–Riemann–Liouville fractional derivative (R), the Riesz–Caputo fractional derivative (RC), the combined Riemann–Liouville fractional derivative (CRL), the combined Caputo fractional derivative (CC), and the classical integer derivative can be deduced as special cases. We list them as follows:

Let $\beta =\alpha$, $\gamma =\frac{1}{2}$ in Equations (9) and (10), so the RVO and the RCVO can be achieved.

When $\alpha \left(·,·\right)\to \alpha$, $\beta \left(·,·\right)\to \beta$, with two constants $0<\alpha ,\beta <1$, the left and right RL, the left and right C, the R, the RC, the CRL, and the CC can be obtained as

$${}_{t_1}{}^{RL}D_t^{\alpha }f\left(t\right)=\frac{1}{\mathsf{\Gamma }\left(1-\alpha \right)}\frac{\mathrm{d}}{\mathrm{d}t}{\displaystyle {\int }_{t_1}^t{\left(t-\tau \right)}^{-\alpha }f\left(\tau \right)\mathrm{d}\tau }$$

(11)

$${}_t{}^{RL}D_{t_2}^{\alpha }f\left(t\right)=\frac{1}{\mathsf{\Gamma }\left(1-\alpha \right)}\left(-\frac{\mathrm{d}}{\mathrm{d}t}\right){\displaystyle {\int }_t^{t_2}{\left(\tau -t\right)}^{-\alpha }f\left(\tau \right)\mathrm{d}\tau }$$

(12)

$${}_{t_1}{}^{C}D_t^{\alpha }f\left(t\right)=\frac{1}{\mathsf{\Gamma }\left(1-\alpha \right)}{\displaystyle {\int }_{t_1}^t{\left(t-\tau \right)}^{-\alpha }\frac{\mathrm{d}}{\mathrm{d}\tau }f\left(\tau \right)\mathrm{d}\tau }$$

(13)

$${}_t{}^{C}D_{t_2}^{\alpha }f\left(t\right)=\frac{1}{\mathsf{\Gamma }\left(1-\alpha \right)}{\displaystyle {\int }_t^{t_2}{\left(\tau -t\right)}^{-\alpha }\left(-\frac{\mathrm{d}}{\mathrm{d}\tau }\right)f\left(\tau \right)\mathrm{d}\tau }$$

(14)

$${}_{t_1}{}^{R}D_{t_2}^{\alpha }f\left(t\right)=\frac{1}{2}\left[{}_{t_1}{}^{RL}D_t^{\alpha }f\left(t\right)-{}_t{}^{RL}D_{t_2}^{\alpha }f\left(t\right)\right]$$

(15)

$${}_{t_1}{}^{RC}D_{t_2}^{\alpha }f\left(t\right)=\frac{1}{2}\left[{}_{t_1}{}^{C}D_t^{\alpha }f\left(t\right)-{}_t{}^{C}D{}_{t_2}^{\alpha }f\left(t\right)\right]$$

(16)

$${}^{RL}D_{\gamma }^{\alpha ,\beta }f\left(t\right)=\gamma {}_{t_1}{}^{RL}D_t^{\alpha }f\left(t\right)-\left(1-\gamma \right){}_t{}^{RL}D_{t_2}^{\beta }f\left(t\right)$$

(17)

$${}^{C}D_{\gamma }^{\alpha ,\beta }f\left(t\right)=\gamma {}_{t_1}{}^{C}D_t^{\alpha }f\left(t\right)-\left(1-\gamma \right){}_t{}^{C}D_{t_2}^{\beta }f\left(t\right)$$

(18)

Equations (11) to (18) can be found in [65].

When $\alpha \left(·,·\right)\to 1$, we have

$${}_{t_1}{}^{RL}D_t^{\alpha \left(·,·\right)}={}_{t_1}{}^{C}D_t^{\alpha \left(·,·\right)}=\mathrm{d}/\mathrm{d}t,{}_t{}^{RL}D_{t_2}^{\alpha \left(·,·\right)}={}_t{}^{C}D_{t_2}^{\alpha \left(·,·\right)}=-\mathrm{d}/\mathrm{d}t$$

(19)

The formulae of integration by parts for the CVO are

$${\displaystyle {\int }_{t_1}^{t_2}g\left(t\right)\left({}_{t_1}{}^{C}D_t^{\alpha \left(·,·\right)}f\left(t\right)\right)}\mathrm{d}t={\displaystyle {\int }_{t_1}^{t_2}f\left(t\right)\left({}_t{}^{RL}D_{t_2}^{\alpha \left(·,·\right)}g\left(t\right)\right)}\mathrm{d}t+\left(f\left(t\right)\left({}_{t}I_{t_2}^{1-\alpha \left(·,·\right)}g\left(t\right)\right)\right)|{}_{t_1}^{t_2}$$

(20)

$${\displaystyle {\int }_{t_1}^{t_2}g\left(t\right)\left({}_t{}^{C}D_{t_2}^{\alpha \left(·,·\right)}f\left(t\right)\right)}\mathrm{d}t={\displaystyle {\int }_{t_1}^{t_2}f\left(t\right)\left({}_{t_1}{}^{RL}D_t^{\alpha \left(·,·\right)}g\left(t\right)\right)}\mathrm{d}t-\left(f\left(t\right)\left({}_{t_1}I_t^{1-\alpha \left(·,·\right)}g\left(t\right)\right)\right)|{}_{t_1}^{t_2}$$

(21)

The formulae of integration by parts for the RLVO, RVO, and RCVO are

$${\displaystyle {\int }_{t_1}^{t_2}g\left(t\right)\left({}_{t_1}{}^{RL}D_t^{\alpha \left(·,·\right)}f\left(t\right)\right)}\mathrm{d}t={\displaystyle {\int }_{t_1}^{t_2}f\left(t\right)\left({}_t{}^{RL}D_{t_2}^{\alpha \left(·,·\right)}g\left(t\right)\right)}\mathrm{d}t-\left[\left({}_{t_1}I_t^{1-\alpha \left(·,·\right)}f\left(t\right)\right)g\left(t\right)\right]|{}_{t_1}$$

(22)

$${\displaystyle {\int }_{t_1}^{t_2}g\left(t\right)\left({}_t{}^{RL}D_{t_2}^{\alpha \left(·,·\right)}f\left(t\right)\right)}\mathrm{d}t={\displaystyle {\int }_{t_1}^{t_2}f\left(t\right)\left({}_{t_1}{}^{RL}D_t^{\alpha \left(·,·\right)}g\left(t\right)\right)}\mathrm{d}t-\left[\left({}_{t}I_{t_2}^{1-\alpha \left(·,·\right)}f\left(t\right)\right)g\left(t\right)\right]|{}_{t_2}$$

(23)

$$\begin{array}{ll}{\displaystyle {\int }_{t_1}^{t_2}g\left(t\right)\left({}_{t_1}{}^{R}D_{t_2}^{\alpha \left(·,·\right)}f\left(t\right)\right)}\mathrm{d}t=& -{\displaystyle {\int }_{t_1}^{t_2}f\left(t\right)\left({}_{t_1}{}^{R}D_{t_2}^{\alpha \left(·,·\right)}g\left(t\right)\right)}\mathrm{d}t+\frac{1}{2}\left[\left({}_{t}I_{t_2}^{1-\alpha \left(·,·\right)}f\left(t\right)\right)g\left(t\right)\right]|{}_{t_2}\\ & -\frac{1}{2}\left[\left({}_{t_1}I_t^{1-\alpha \left(·,·\right)}f\left(t\right)\right)g\left(t\right)\right]|{}_{t_1}\end{array}$$

(24)

$${\displaystyle {\int }_{t_1}^{t_2}g\left(t\right)\left({}_{t_1}{}^{RC}D_{t_2}^{\alpha \left(·,·\right)}f\left(t\right)\right)}\mathrm{d}t=-{\displaystyle {\int }_{t_1}^{t_2}f\left(t\right)\left({}_{t_1}{}^{R}D_{t_2}^{\alpha \left(·,·\right)}g\left(t\right)\right)}\mathrm{d}t+\left[f\left(t\right)\left({}_{t_1}{}^{R}I_{t_2}^{1-\alpha \left(·,·\right)}g\left(t\right)\right)\right]|{}_{t_1}^{t_2}$$

(25)

For more details about the above Equations (19)–(25), please refer to [54,57,65].

3. Hamilton Equation of Variable Order

A unified Hamilton action and Hamilton principle of variable order will be given first, and then the Hamilton action and Hamilton principle in terms of the CRLVO and the CCVO will be discussed, respectively.

Assuming the Lagrangian of variable order is $L_U\left(t,q_j,D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_j\right)$, then the generalized momentum of variable order and the Hamiltonian of variable order can be listed as

$$\begin{array}{c}{p}_{Ui}=\frac{\partial L_U\left(t,q_j,D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_j\right)}{\partial D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i},\\ H_U=p_{Ui}·D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i-L_U\left(t,q_j,D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_j\right), i,j=1,2,\cdots ,n\end{array}$$

(26)

And the functional

$$S_{HU}={\displaystyle {\int }_{t_1}^{t_2}\left[p_{Ui}·D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i-H_U\left(t,q_j,p_{Uj}\right)\right]}\mathrm{d}t$$

(27)

Represents a unified Hamilton action of variable order, where $q_j$, $p_{Uj}\in {\mathrm{C}}^1$, $H_U$ is a ${\mathrm{C}}^1$ function with respect to all its arguments.

The formula

$$\delta S_{HU}=\delta {\displaystyle {\int }_{t_1}^{t_2}\left[p_{Ui}·D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i-H_U\left(t,q_j,p_{Uj}\right)\right]}\mathrm{d}t=0$$

(28)

with [57]

$$\delta {}_{t_1}D_t^{\alpha \left(·,·\right)}{q}_i={}_{t_1}D_t^{\alpha \left(·,·\right)}\delta q_i, \delta {}_{t}D_{t_2}^{\beta \left(·,·\right)}{q}_i={}_{t}D_{t_2}^{\beta \left(·,·\right)}\delta q_i, \delta q_i|{}_{t=t_1}=\delta q_i|{}_{t=t_2}=0$$

(29)

is called a unified Hamilton principle of variable order.

3.1. Hamilton Equation with the CRLVO

When we consider the CRLVO, the Lagrangian of variable order can be denoted by $L_{RL}\left(t,q_j,{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_j\right)$. The generalized momentum of variable order and the Hamiltonian of variable order are

$$\begin{array}{c}{p}_{RLi}=\frac{\partial L_{RL}\left(t,q_j,{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_j\right)}{\partial {}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i},\\ H_{RL}=p_{RLi}·{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i-L_{RL}\left(t,q_j,{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_j\right)\end{array}$$

(30)

Under this condition, Hamilton action of variable order is

$$S_{HRL}={\displaystyle {\int }_{t_1}^{t_2}\left[p_{RLi}·{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i-H_{RL}\left(t,q_j,p_{RLj}\right)\right]}\mathrm{d}t$$

(31)

Hamilton principle of variable order is

$$\delta S_{HRL}=\delta {\displaystyle {\int }_{t_1}^{t_2}\left[p_{RLi}·{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i-H_{RL}\left(t,q_j,p_{RLj}\right)\right]}\mathrm{d}t=0$$

(32)

with [57]

$$\delta {}_{t_1}{}^{RL}D_t^{\alpha \left(·,·\right)}{q}_i={}_{t_1}{}^{RL}D_t^{\alpha \left(·,·\right)}\delta q_i, \delta {}_t{}^{RL}D_{t_2}^{\beta \left(·,·\right)}{q}_i={}_t{}^{RL}D_{t_2}^{\beta \left(·,·\right)}\delta q_i, \delta q_i|{}_{t=t_1}=\delta q_i|{}_{t=t_2}=0$$

(33)

It follows from Equation (32) that

$$\begin{array}{l}\delta S_{HRL}=\delta {\displaystyle {\int }_{t_1}^{t_2}\left[p_{RLi}·{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i-H_{RL}\left(t,q_j,p_{RLj}\right)\right]}\mathrm{d}t\\ ={\displaystyle {\int }_{t_1}^{t_2}\left[\delta p_{RLi}·{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i+p_{RLi}·\delta {}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i-\delta H_{RL}\left(t,q_j,p_{RLj}\right)\right]\mathrm{d}t}\\ ={\displaystyle {\int }_{t_1}^{t_2}\left(\delta p_{RLi}·{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i-\delta q_i·{}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}{p}_{RLi}-\frac{\partial H_{RL}}{\partial q_i}\delta q_i-\frac{\partial H_{RL}}{\partial p_{RLi}}\delta p_{RLi}\right)\mathrm{d}t}\\ ={\displaystyle {\int }_{t_1}^{t_2}\left[\delta p_{RLi}·\left({}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i-\frac{\partial H_{RL}}{\partial p_{RLi}}\right)-\delta q_i·\left({}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}{p}_{RLi}+\frac{\partial H_{RL}}{\partial q_i}\right)\right]\mathrm{d}t}=0\end{array}$$

(34)

where

$$\begin{array}{l}{\displaystyle {\int }_{t_1}^{t_2}{p}_{RLi}·\delta {}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i\mathrm{d}t}={\displaystyle {\int }_{t_1}^{t_2}{p}_{RLi}·\delta \left[\gamma {}_{t_1}{}^{RL}D_t^{\alpha \left(·,·\right)}{q}_i-\left(1-\gamma \right){}_t{}^{RL}D_{t_2}^{\beta \left(·,·\right)}{q}_i\right]\mathrm{d}t}\\ ={\displaystyle {\int }_{t_1}^{t_2}\left[\gamma p_{RLi}·{}_{t_1}{}^{RL}D_t^{\alpha \left(·,·\right)}\delta q_i-\left(1-\gamma \right)p_{RLi}·{}_t{}^{RL}D_{t_2}^{\beta \left(·,·\right)}\delta q_i\right]\mathrm{d}t}\\ ={\displaystyle {\int }_{t_1}^{t_2}\left[\delta q_i·\gamma {}_t{}^{RL}D_{t_2}^{\alpha \left(·,·\right)}{p}_{RLi}-\delta q_i·\left(1-\gamma \right){}_{t_1}{}^{RL}D_t^{\beta \left(·,·\right)}{p}_{RLi}\right]\mathrm{d}t}\\ =-{\displaystyle {\int }_{t_1}^{t_2}\left[\delta q_i·{}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}{p}_{RLi}\right]\mathrm{d}t}\end{array}$$

(35)

From Equations (30) and (34), we have

$${}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i=\frac{\partial H_{RL}}{\partial p_{RLi}}, {}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}{p}_{RLi}=-\frac{\partial H_{RL}}{\partial q_i}$$

(36)

Equation (36) is called the Hamilton equation with the CRLVO.

3.2. Hamilton Equation with the CCVO

When we consider the CCVO, the Lagrangian of variable order can be denoted by $L_C\left(t,q_j,{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_j\right)$. The generalized momentum of variable order and the Hamiltonian of variable order are

$$p_{Ci}=\frac{\partial L_C\left(t,q_j,{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_j\right)}{\partial {}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i}, H_C=p_{Ci}·{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i-L_C\left(t,q_j,{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_j\right)$$

(37)

Under this condition, the Hamilton action of variable order is

$$S_{HC}={\displaystyle {\int }_{t_1}^{t_2}\left[p_{Ci}·{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i-H_C\left(t,q_j,p_{Cj}\right)\right]}\mathrm{d}t$$

(38)

The Hamilton principle of variable order is

$$\delta S_{HC}=\delta {\displaystyle {\int }_{t_1}^{t_2}\left[p_{Ci}·{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i-H_C\left(t,q_j,p_{Cj}\right)\right]}\mathrm{d}t=0$$

(39)

with [57]

$$\delta {}_{t_1}{}^{C}D_t^{\alpha \left(·,·\right)}{q}_i={}_{t_1}{}^{C}D_t^{\alpha \left(·,·\right)}\delta q_i, \delta {}_t{}^{C}D_{t_2}^{\beta \left(·,·\right)}{q}_i={}_t{}^{C}D_{t_2}^{\beta \left(·,·\right)}\delta q_i, \delta q_i|{}_{t=t_1}=\delta q_i|{}_{t=t_2}=0$$

(40)

Similarly, from Equations (37) and (39), we can get

$${}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i=\frac{\partial H_C}{\partial p_{Ci}}, {}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}{p}_{Ci}=-\frac{\partial H_C}{\partial q_i}$$

(41)

Equation (41) is called the Hamilton equation with the CCVO.

Remark 1.

When$\beta \left(·,·\right)=\alpha \left(·,·\right)$,$\gamma =\frac{1}{2}$, the Hamilton equation in terms of the RVO can be obtained from Equation (36) as

$${}_{t_1}{}^{R}D_{t_2}^{\alpha \left(·,·\right)}{q}_i=\frac{\partial H_R}{\partial p_{Ri}}, {}_{t_1}{}^{R}D_{t_2}^{\alpha \left(·,·\right)}{p}_{Ri}=-\frac{\partial H_R}{\partial q_i}$$

(42)

where$p_{Ri}=\frac{\partial L_R\left(t,q_j,{}_{t_1}{}^{R}D_{t_2}^{\alpha \left(·,·\right)}{q}_j\right)}{\partial {}_{t_1}{}^{R}D_{t_2}^{\alpha \left(·,·\right)}{q}_i}$,$H_R=p_{Ri}·{}_{t_1}{}^{R}D_{t_2}^{\alpha \left(·,·\right)}{q}_i-L_R\left(t,q_j,{}_{t_1}{}^{R}D_{t_2}^{\alpha \left(·,·\right)}{q}_j\right)$.

Similarly, the Hamilton equation in terms of the RCVO can be obtained from Equation (41) as

$${}_{t_1}{}^{RC}D_{t_2}^{\alpha \left(·,·\right)}{q}_i=\frac{\partial H_{RC}}{\partial p_{RCi}}, {}_{t_1}{}^{R}D_{t_2}^{\alpha \left(·,·\right)}{p}_{RCi}=-\frac{\partial H_{RC}}{\partial q_i}$$

(43)

where$p_{RCi}=\frac{\partial L_{RC}\left(t,q_j,{}_{t_1}{}^{RC}D_{t_2}^{\alpha \left(·,·\right)}{q}_j\right)}{\partial {}_{t_1}{}^{RC}D_{t_2}^{\alpha \left(·,·\right)}{q}_i}$,$H_{RC}=p_{RCi}·{}_{t_1}{}^{RC}D_{t_2}^{\alpha \left(·,·\right)}{q}_i-L_{RC}\left(t,q_j,{}_{t_1}{}^{RC}D_{t_2}^{\alpha \left(·,·\right)}{q}_j\right)$.

Moreover, Equations (42) and (43) can also be obtained through their respective Hamilton principles.

Remark 2.

From Equations (30) and (36), the Lagrange equation with the CRLVO can be obtained as

$$\frac{\partial L_{RL}\left(t,q_j,{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_j\right)}{\partial q_i}-{}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}\frac{\partial L_{RL}\left(t,q_j,{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_j\right)}{\partial {}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i}=0$$

(44)

It is noted that Equation (44) can be obtained through

$$\begin{array}{c}\delta S_{LRL}=\delta {\displaystyle {\int }_{t_1}^{t_2}{L}_{RL}\left(t,q_j,{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_j\right)}\mathrm{d}t=0,\\ \delta {}_{t_1}{}^{RL}D_t^{\alpha \left(·,·\right)}{q}_i={}_{t_1}{}^{RL}D_t^{\alpha \left(·,·\right)}\delta q_i, \delta {}_t{}^{RL}D_{t_2}^{\beta \left(·,·\right)}{q}_i={}_t{}^{RL}D_{t_2}^{\beta \left(·,·\right)}\delta q_i, \delta q_i|{}_{t=t_1}=\delta q_i|{}_{t=t_2}=0\end{array}$$

(45)

Similarly, from Equations (37) and (41), the Lagrange equation with the CCVO can be obtained as

$$\frac{\partial L_C}{\partial q_i}-{}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}\frac{\partial L_C\left(t,q_j,{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_j\right)}{\partial {}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i}=0$$

(46)

And Equation (46) can also be obtained through

$$\begin{array}{c}\delta S_{LC}=\delta {\displaystyle {\int }_{t_1}^{t_2}{L}_C\left(t,q_j,{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_j\right)}\mathrm{d}t=0, \delta {}_{t_1}{}^{C}D_t^{\alpha \left(·,·\right)}{q}_i={}_{t_1}{}^{C}D_t^{\alpha \left(·,·\right)}\delta q_i,\\ \delta {}_t{}^{C}D_{t_2}^{\beta \left(·,·\right)}{q}_i={}_t{}^{C}D_{t_2}^{\beta \left(·,·\right)}\delta q_i, \delta q_i|{}_{t=t_1}=\delta q_i|{}_{t=t_2}=0\end{array}$$

(47)

Remark 3.

When$\alpha \left(·,·\right)\to \alpha$,$\beta \left(·,·\right)\to \beta$, Hamilton equations of variable order, i.e., Equations (36) and (41) reduce to

$${}^{RL}D_{\gamma }^{\alpha ,\beta }{q}_i=\frac{\partial H_{RLF}}{\partial p_{RLiF}}, {}^{RL}D_{1-\gamma }^{\beta ,\alpha }{p}_{RLiF}=-\frac{\partial H_{RLF}}{\partial q_i}$$

(48)

and

$${}^{C}D_{\gamma }^{\alpha ,\beta }{q}_i=\frac{\partial H_{CF}}{\partial p_{CiF}}, {}^{RL}D_{1-\gamma }^{\beta ,\alpha }{p}_{CiF}=-\frac{\partial H_{CF}}{\partial q_i}$$

(49)

Equations (48) and (49) are the Hamilton equations with the CRL and the CC, respectively.

Lagrange equations of variable order, i.e., Equations (44) and (46) reduce to

$$\frac{\partial L_{RLF}\left(t,q_j,{}^{RL}D_{\gamma }^{\alpha ,\beta }{q}_j\right)}{\partial q_i}-{}^{RL}D_{1-\gamma }^{\beta ,\alpha }\frac{\partial L_{RLF}\left(t,q_j,{}^{RL}D_{\gamma }^{\alpha ,\beta }{q}_j\right)}{\partial {}^{RL}D_{\gamma }^{\alpha ,\beta }{q}_i}=0$$

(50)

and

$$\frac{\partial L_{CF}}{\partial q_i}-{}^{RL}D_{1-\gamma }^{\beta ,\alpha }\frac{\partial L_{CF}\left(t,q_j,{}^{C}D_{\gamma }^{\alpha ,\beta }{q}_j\right)}{\partial {}^{C}D_{\gamma }^{\alpha ,\beta }{q}_i}=0$$

(51)

Equations (50) and (51) are the Lagrange equations with the CRL and the CC, respectively. It is noted that Equations (48) to (51) are consistent with the results in [53].

Remark 4.

When$\alpha \left(·,·\right)$,$\beta \left(·,·\right)\to 1$, Hamilton equations of variable order, i.e., Equations (36) and (41) reduce to

$${\dot{q}}_i=\frac{\partial H\left(t,q_j,p_j\right)}{\partial p_i}, {\dot{p}}_i=-\frac{\partial H\left(t,q_j,p_j\right)}{\partial q_i}$$

(52)

Lagrange equations of variable order, i.e., Equations (44) and (46) reduce to

$$\frac{\partial L\left(t,q_j,{\dot{q}}_j\right)}{\partial q_i}-\frac{\mathrm{d}}{\mathrm{d}t}\frac{\partial L\left(t,q_j,{\dot{q}}_j\right)}{\partial {\dot{q}}_i}=0$$

(53)

Equations (52) and (53) are the classical Hamilton equation and the classical Lagrange equation, respectively, which is consistent with the results in [66].

4. Noether Symmetry and Conserved Quantity

The classical Noether symmetry means the invariance of the classical Hamilton action under the infinitesimal transformations of time and coordinates. A Noether symmetry can lead to a conserved quantity. Here, we generalize the classical definition of Noether symmetry to fractional form with variable order. Then Noether symmetry and conserved quantity will be discussed in detail with the CRLVO and the CCVO, respectively.

Under the infinitesimal transformations

$$\overline{t}=t+\Delta t, {\overline{q}}_i\left(\overline{t}\right)=q_i\left(t\right)+\Delta q_i, {\overline{p}}_{Ui}\left(\overline{t}\right)=p_{Ui}\left(t\right)+\Delta p_{Ui}$$

(54)

Whose linear parts are

$$\begin{array}{c}\overline{t}=t+\theta {\xi }_{U0}^0\left(t,q_j,p_{Uj}\right), {\overline{q}}_i\left(\overline{t}\right)=q_i\left(t\right)+\theta {\xi }_{Ui}^0\left(t,q_j,p_{Uj}\right),\\ {\overline{p}}_{Ui}\left(\overline{t}\right)=p_{Ui}\left(t\right)+\theta {\eta }_{Ui}^0\left(t,q_j,p_{Uj}\right)\end{array}$$

(55)

The unified Hamilton action of variable order, i.e., Equation (27) changes to be

$${\overline{S}}_{HU}={\displaystyle {\int }_{{\overline{t}}_1}^{{\overline{t}}_2}\left[{\overline{p}}_{Ui}·D_{\gamma }^{\overline{\alpha }\left(·,·\right),\overline{\beta }\left(·,·\right)}{\overline{q}}_i-H_U\left(\overline{t},{\overline{q}}_j,{\overline{p}}_{Uj}\right)\right]}\mathrm{d}\overline{t}$$

(56)

where $\theta$ is an infinitesimal parameter, and ${\xi }_{U0}^0$, ${\xi }_{Ui}^0$, and ${\eta }_{Ui}^0$ are infinitesimal generators.

Let $\Delta S_{HU}$ denote the linear part of ${\overline{S}}_{HU}-S_{HU}$, if there exists infinitesimal generators satisfying $\Delta S_{HU}=-{\displaystyle {\int }_{t_1}^{t_2}\frac{\mathrm{d}}{\mathrm{d}t}\left(\Delta G_U^0\right)}\mathrm{d}t$, then fractional Noether symmetry with variable order can be determined by the corresponding infinitesimal transformations, where $\Delta G_U^0=\theta G_U^0\left(t,q_j,p_{Uj}\right)$, and $G_U^0$ is called a gauge function.

For a certain mechanical system, conserved quantity can be deduced from Noether symmetry. A conserved quantity can be seen as a solution to the differential equation of motion of the mechanical system, and help reduce the freedom of the mechanical system.

A quantity $I_U$ is called a conserved quantity if and only if $\frac{\mathrm{d}}{\mathrm{d}t}{I}_U=0$.

4.1. Noether Symmetry and Conserved Quantity with the CRLVO

When we consider the CRLVO, the infinitesimal transformations can be expressed as

$$\begin{array}{c}\overline{t}=t+\Delta t, {\overline{q}}_i\left(\overline{t}\right)=q_i\left(t\right)+\Delta q_i, {\overline{p}}_{RLi}\left(\overline{t}\right)=p_{RLi}\left(t\right)+\Delta p_{RLi}\\ (\overline{t}=t+{\theta }_{RL}{\xi }_{RL0}^0\left(t,q_j,p_{RLj}\right), {\overline{q}}_i\left(\overline{t}\right)=q_i\left(t\right)+{\theta }_{RL}{\xi }_{RLi}^0\left(t,q_j,p_{RLj}\right),\\ {\overline{p}}_{RLi}\left(\overline{t}\right)=p_{RLi}\left(t\right)+{\theta }_{RL}{\eta }_{RLi}^0\left(t,q_j,p_{RLj}\right)).\end{array}$$

(57)

Under the action of Equation (57), the Hamilton action of variable order, i.e., Equation (31) becomes

$${\overline{S}}_{HRL}={\displaystyle {\int }_{{\overline{t}}_1}^{{\overline{t}}_2}\left[{\overline{p}}_{RLi}·{}^{RL}D_{\gamma }^{\overline{\alpha }\left(·,·\right),\overline{\beta }\left(·,·\right)}{\overline{q}}_i-H_{RL}\left(\overline{t},{\overline{q}}_j,{\overline{p}}_{RLj}\right)\right]}\mathrm{d}\overline{t}$$

(58)

Then we get

$$\begin{array}{c}\Delta S_{HRL}={\overline{S}}_{HRL}-S_{HRL}\\ ={\displaystyle {\int }_{t_1}^{t_2}\{-H_{RL}\left(t+\Delta t,q_j+\Delta q_j,p_{RLj}+\Delta p_{RLj}\right)+[{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i+\Delta t\frac{\mathrm{d}}{\mathrm{d}t}{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i}\\ -\gamma q_i\left(t_1\right)\Delta t_1\frac{\mathrm{d}}{\mathrm{d}t}\frac{{\left(t-t_1\right)}^{-\alpha \left(t,t_1\right)}}{\mathsf{\Gamma }\left(1-\alpha \left(t,t_1\right)\right)}+{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}\delta q_i+\left(1-\gamma \right)q_i\left(t_2\right)\Delta t_2\frac{\mathrm{d}}{\mathrm{d}t}\frac{{\left(t_2-t\right)}^{-\beta \left(t_2,t\right)}}{\mathsf{\Gamma }\left(1-\beta \left(t_2,t\right)\right)}]\\ ·\left(p_{RLi}+\Delta p_{RLi}\right)\}·\left(1+\frac{\mathrm{d}}{\mathrm{d}t}\Delta t\right)\mathrm{d}t-{\displaystyle {\int }_{t_1}^{t_2}\left[p_{RLi}·{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i-H_{RL}\left(t,q_j,p_{RLj}\right)\right]}\mathrm{d}t\\ ={\displaystyle {\int }_{t_1}^{t_2}[p_{RLi}·{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}\delta q_i+\left(p_{RLi}·\frac{\mathrm{d}}{\mathrm{d}t}{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i-\frac{\partial H_{RL}}{\partial t}\right)\Delta t}\\ -\gamma p_{RLi}{q}_i\left(t_1\right)\Delta t_1\frac{\mathrm{d}}{\mathrm{d}t}\frac{{\left(t-t_1\right)}^{-\alpha \left(t,t_1\right)}}{\mathsf{\Gamma }\left(1-\alpha \left(t,t_1\right)\right)}+\left(p_{RLi}·{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i-H_{RL}\right)·\frac{\mathrm{d}}{\mathrm{d}t}\Delta t\\ +\left(1-\gamma \right)p_{RLi}{q}_i\left(t_2\right)\Delta t_2\frac{\mathrm{d}}{\mathrm{d}t}\frac{{\left(t_2-t\right)}^{-\beta \left(t_2,t\right)}}{\mathsf{\Gamma }\left(1-\beta \left(t_2,t\right)\right)}-\frac{\partial H_{RL}}{\partial q_i}\Delta q_i]\mathrm{d}t\end{array}$$

(59)

where

$$\begin{array}{c}{}^{RL}D_{\gamma }^{\overline{\alpha }\left(·,·\right),\overline{\beta }\left(·,·\right)}{\overline{q}}_i={}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i+{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}\delta q_i+\Delta t\frac{\mathrm{d}}{\mathrm{d}t}{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i\\ -\gamma q_i\left(t_1\right)\Delta t_1\frac{\mathrm{d}}{\mathrm{d}t}\frac{{\left(t-t_1\right)}^{-\alpha \left(t,t_1\right)}}{\mathsf{\Gamma }\left(1-\alpha \left(t,t_1\right)\right)}+\left(1-\gamma \right)q_i\left(t_2\right)\Delta t_2\frac{\mathrm{d}}{\mathrm{d}t}\frac{{\left(t_2-t\right)}^{-\beta \left(t_2,t\right)}}{\mathsf{\Gamma }\left(1-\beta \left(t_2,t\right)\right)}\end{array}$$

(60)

Let $\Delta S_{HRL}=-{\displaystyle {\int }_{t_1}^{t_2}\frac{\mathrm{d}}{\mathrm{d}t}\left(\Delta G_{RL}^0\right)}\mathrm{d}t$, and we obtain

$$\begin{array}{c}{p}_{RLi}·{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}\left({\xi }_{RLi}^0-{\dot{q}}_i{\xi }_{RL0}^0\right)+\left(p_{RLi}·\frac{\mathrm{d}}{\mathrm{d}t}{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i-\frac{\partial H_{RL}}{\partial t}\right){\xi }_{RL0}^0\\ +\left(p_{RLi}·{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i-H_{RL}\right)·{\dot{\xi }}_{RL0}^0-\frac{\partial H_{RL}}{\partial q_i}{\xi }_{RLi}^0+{\dot{G}}_{RL}^0\left(t,q_j,p_{RLj}\right)\\ -\gamma p_{RLi}{q}_i\left(t_1\right){\xi }_{RL0}^0\left(t_1,q_j\left(t_1\right),p_{RLj}\left(t_1\right)\right)\frac{\mathrm{d}}{\mathrm{d}t}\frac{{\left(t-t_1\right)}^{-\alpha \left(t,t_1\right)}}{\mathsf{\Gamma }\left(1-\alpha \left(t,t_1\right)\right)}\\ +\left(1-\gamma \right)p_{RLi}{q}_i\left(t_2\right){\xi }_{RL0}^0\left(t_2,q_j\left(t_2\right),p_{RLj}\left(t_2\right)\right)\frac{\mathrm{d}}{\mathrm{d}t}\frac{{\left(t_2-t\right)}^{-\beta \left(t_2,t\right)}}{\mathsf{\Gamma }\left(1-\beta \left(t_2,t\right)\right)}=0\end{array}$$

(61)

Equation (61) is called Noether identity with the CRLVO.

Theorem 1.

For the Hamiltonian system with the CRLVO (Equation (36)), if there exists a gauge function$G_{RL}^0$such that the infinitesimal generators${\xi }_{RL0}^0$and${\xi }_{RLi}^0$satisfy Equation (61), then there exists a conserved quantity.

$$\begin{array}{c}{I}_{HRL0}={\displaystyle {\int }_{t_1}^t\left[p_{RLi}·{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}\left({\xi }_{RLi}^0-{\dot{q}}_i{\xi }_{RL0}^0\right)+\left({\xi }_{RLi}^0-{\dot{q}}_i{\xi }_{RL0}^0\right)·{}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}{p}_{RLi}\right]}\mathrm{d}\tau +(p_{RLi}·\\ {}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i-H_{RL})·{\xi }_{RL0}^0-{\displaystyle {\int }_{t_1}^t{p}_{RLi}[\gamma q_i\left(t_1\right){\xi }_{RL0}^0\left(t_1,q_j\left(t_1\right),p_{RLj}\left(t_1\right)\right)\frac{\mathrm{d}}{\mathrm{d}\tau }\frac{{\left(\tau -t_1\right)}^{-\alpha \left(\tau ,t_1\right)}}{\mathsf{\Gamma }\left(1-\alpha \left(\tau ,t_1\right)\right)}}\\ -\left(1-\gamma \right)q_i\left(t_2\right){\xi }_{RL0}^0\left(t_2,q_j\left(t_2\right),p_{RLj}\left(t_2\right)\right)\frac{\mathrm{d}}{\mathrm{d}\tau }\frac{{\left(t_2-\tau \right)}^{-\beta \left(t_2,\tau \right)}}{\mathsf{\Gamma }\left(1-\beta \left(t_2,\tau \right)\right)}]\mathrm{d}\tau +G_{RL}^0\end{array}$$

(62)

Proof. 

Using Equations (36) and (61), we have

$$\frac{\mathrm{d}}{\mathrm{d}t}{I}_{HRL0}=\left({\dot{p}}_{RLi}·{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i+p_{RLi}·\frac{\mathrm{d}}{\mathrm{d}t}{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i-\frac{\partial H_{RL}}{\partial t}-\frac{\partial H_{RL}}{\partial q_i}{\dot{q}}_i-\frac{\partial H_{RL}}{\partial p_{RLi}}{\dot{p}}_{RLi}\right)$$

$$·{\xi }_{RL0}^0+\left(p_{RLi}·{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i-H_{RL}\right)·{\dot{\xi }}_{RL0}^0+p_{RLi}·{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}\left({\xi }_{RLi}^0-{\dot{q}}_i{\xi }_{RL0}^0\right)+({\xi }_{RLi}^0$$

$$-{\dot{q}}_i{\xi }_{RL0}^0)·{}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}{p}_{RLi}-\gamma p_{RLi}{q}_i\left(t_1\right){\xi }_{RL0}^0\left(t_1,q_j\left(t_1\right),p_{RLj}\left(t_1\right)\right)\frac{\mathrm{d}}{\mathrm{d}t}\frac{{\left(t-t_1\right)}^{-\alpha \left(t,t_1\right)}}{\mathsf{\Gamma }\left(1-\alpha \left(t,t_1\right)\right)}$$

$$+\left(1-\gamma \right)p_{RLi}{q}_i\left(t_2\right){\xi }_{RL0}^0\left(t_2,q_j\left(t_2\right),p_{RLj}\left(t_2\right)\right)\frac{\mathrm{d}}{\mathrm{d}t}\frac{{\left(t_2-t\right)}^{-\beta \left(t_2,t\right)}}{\mathsf{\Gamma }\left(1-\beta \left(t_2,t\right)\right)}+{\dot{G}}_{RL}^0$$

$$=\left({\dot{p}}_{RLi}·{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i-\frac{\partial H_{RL}}{\partial q_i}{\dot{q}}_i-\frac{\partial H_{RL}}{\partial p_{RLi}}{\dot{p}}_{RLi}\right)·{\xi }_{RL0}^0+\left({\xi }_{RLi}^0-{\dot{q}}_i{\xi }_{RL0}^0\right)·{}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}{p}_{RLi}+\frac{\partial H_{RL}}{\partial q_i}{\xi }_{RLi}^0$$

$$=\left({\xi }_{RLi}^0-{\dot{q}}_i{\xi }_{RL0}^0\right)·\left({}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}{p}_{RLi}+\frac{\partial H_{RL}}{\partial q_i}\right)=0$$

 □

4.2. Noether Symmetry and Conserved Quantity with the CCVO

When we consider the CCVO, the infinitesimal transformations can be expressed as

$$\begin{array}{l}\overline{t}=t+\Delta t, {\overline{q}}_i\left(\overline{t}\right)=q_i\left(t\right)+\Delta q_i, {\overline{p}}_{Ci}\left(\overline{t}\right)=p_{Ci}\left(t\right)+\Delta p_{Ci},\\ (\overline{t}=t+{\theta }_C{\xi }_{C0}^0\left(t,q_j,p_{Cj}\right), {\overline{q}}_i\left(\overline{t}\right)=q_i\left(t\right)+{\theta }_C{\xi }_{Ci}^0\left(t,q_j,p_{Cj}\right),\\ {\overline{p}}_{Ci}\left(\overline{t}\right)=p_{Ci}\left(t\right)+{\theta }_C{\eta }_{Ci}^0\left(t,q_j,p_{Cj}\right).\end{array}$$

(63)

Under the action of Equation (63), the Hamilton action of variable order, i.e., Equation (38) becomes

$${\overline{S}}_{HC}={\displaystyle {\int }_{{\overline{t}}_1}^{{\overline{t}}_2}\left[{\overline{p}}_{Ci}·{}^{C}D_{\gamma }^{\overline{\alpha }\left(·,·\right),\overline{\beta }\left(·,·\right)}{\overline{q}}_i-H_C\left(\overline{t},{\overline{q}}_j,{\overline{p}}_{Cj}\right)\right]}\mathrm{d}\overline{t}$$

(64)

Then let $\Delta S_{HC}={\overline{S}}_{HC}-S_{HC}=-{\displaystyle {\int }_{t_1}^{t_2}\frac{\mathrm{d}}{\mathrm{d}t}\left(\Delta G_C^0\right)}\mathrm{d}t$, and we obtain

$$\begin{array}{c}{p}_{Ci}·{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}\left({\xi }_{Ci}^0-{\dot{q}}_i{\xi }_{C0}^0\right)+\left(p_{Ci}·\frac{\mathrm{d}}{\mathrm{d}t}{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i-\frac{\partial H_C}{\partial t}\right){\xi }_{C0}^0\\ -\gamma p_{Ci}{\xi }_{C0}^0\left(t_1,q_j\left(t_1\right),p_{Cj}\left(t_1\right)\right)\frac{{\left(t-t_1\right)}^{-\alpha \left(t,t_1\right)}}{\mathsf{\Gamma }\left(1-\alpha \left(t,t_1\right)\right)}{\dot{q}}_i\left(t_1\right)+\left(p_{Ci}·{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i-H_C\right)·{\dot{\xi }}_{C0}^0\\ +\left(1-\gamma \right)p_{Ci}{\xi }_{C0}^0\left(t_2,q_j\left(t_2\right),p_{Cj}\left(t_2\right)\right)\frac{{\left(t_2-t\right)}^{-\beta \left(t_2,t\right)}}{\mathsf{\Gamma }\left(1-\beta \left(t_2,t\right)\right)}{\dot{q}}_i\left(t_2\right)-\frac{\partial H_C}{\partial q_i}{\xi }_{Ci}^0+{\dot{G}}_C^0\left(t,q_j,p_{Cj}\right)=0\end{array}$$

(65)

Equation (65) is called Noether identity with the CCVO.

Theorem 2.

For the Hamiltonian system with the CCVO (Equation (41)), if there exists a gauge function$G_C^0$such that the infinitesimal generators${\xi }_{C0}^0$and${\xi }_{Ci}^0$satisfy Equation (65), then there exists a conserved quantity.

$$\begin{array}{c}{I}_{HC0}={\displaystyle {\int }_{t_1}^t\left[p_{Ci}·{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}\left({\xi }_{Ci}^0-{\dot{q}}_i{\xi }_{C0}^0\right)+\left({\xi }_{Ci}^0-{\dot{q}}_i{\xi }_{C0}^0\right)·{}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}{p}_{Ci}\right]}\mathrm{d}\tau +(p_{Ci}\\ ·{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i-H_C)·{\xi }_{C0}^0-{\displaystyle {\int }_{t_1}^t{p}_{Ci}[\gamma {\xi }_{C0}^0\left(t_1,q_j\left(t_1\right),p_{Cj}\left(t_1\right)\right){\dot{q}}_i\left(t_1\right)\frac{{\left(\tau -t_1\right)}^{-\alpha \left(\tau ,t_1\right)}}{\mathsf{\Gamma }\left(1-\alpha \left(\tau ,t_1\right)\right)}}\\ -\left(1-\gamma \right){\dot{q}}_i\left(t_2\right){\xi }_{C0}^0\left(t_2,q_j\left(t_2\right),p_{Cj}\left(t_2\right)\right)\frac{{\left(t_2-\tau \right)}^{-\beta \left(t_2,\tau \right)}}{\mathsf{\Gamma }\left(1-\beta \left(t_2,\tau \right)\right)}]\mathrm{d}\tau +G_C^0\end{array}$$

(66)

Proof. 

The proof is similar to that in Theorem 1 and thus is omitted. □

4.3. Some Explanations

We will present some special cases in this section on the basis of Theorems 1 and 2, which are the main results of fractional Noether symmetry and fractional conserved quantity with variable order in this paper.

Remark 5.

Let$\beta \left(·,·\right)=\alpha \left(·,·\right)$,$\gamma =\frac{1}{2}$, and conserved quantities for the Hamiltonian systems with variable order (Equations (42) and (43)) can be achieved as follows:

Theorem 3.

For the Hamiltonian system with the RVO (Equation (42)), if there exists a gauge function$G_R^0\left(t,q_j,p_{Rj}\right)$such that the infinitesimal generators${\xi }_{R0}^0\left(t,q_j,p_{Rj}\right)$and${\xi }_{Ri}^0\left(t,q_j,p_{Rj}\right)$satisfy

$$\begin{array}{c}{p}_{Ri}·{}_{t_1}{}^{R}D_{t_2}^{\alpha \left(·,·\right)}\left({\xi }_{Ri}^0-{\dot{q}}_i{\xi }_{R0}^0\right)+\left(p_{Ri}·\frac{\mathrm{d}}{\mathrm{d}t}{}_{t_1}{}^{R}D_{t_2}^{\alpha \left(·,·\right)}{q}_i-\frac{\partial H_R}{\partial t}\right){\xi }_{R0}^0+\left(p_{Ri}·{}_{t_1}{}^{R}D_{t_2}^{\alpha \left(·,·\right)}{q}_i-H_R\right)·{\dot{\xi }}_{R0}^0\\ -\gamma p_{Ri}{q}_i\left(t_1\right){\xi }_{R0}^0\left(t_1,q_j\left(t_1\right),p_{Rj}\left(t_1\right)\right)\frac{\mathrm{d}}{\mathrm{d}t}\frac{{\left(t-t_1\right)}^{-\alpha \left(t,t_1\right)}}{\mathsf{\Gamma }\left(1-\alpha \left(t,t_1\right)\right)}-\frac{\partial H_R}{\partial q_i}{\xi }_{Ri}^0+{\dot{G}}_R^0\\ +\left(1-\gamma \right)p_{Ri}{q}_i\left(t_2\right){\xi }_{R0}^0\left(t_2,q_j\left(t_2\right),p_{Rj}\left(t_2\right)\right)\frac{\mathrm{d}}{\mathrm{d}t}\frac{{\left(t_2-t\right)}^{-\beta \left(t_2,t\right)}}{\mathsf{\Gamma }\left(1-\beta \left(t_2,t\right)\right)}=0\end{array}$$

(67)

then there exists a conserved quantity.

$$\begin{array}{c}{I}_{HR0}={\displaystyle {\int }_{t_1}^t\left[p_{Ri}·{}_{t_1}{}^{R}D_{t_2}^{\alpha \left(·,·\right)}\left({\xi }_{Ri}^0-{\dot{q}}_i{\xi }_{R0}^0\right)+\left({\xi }_{Ri}^0-{\dot{q}}_i{\xi }_{R0}^0\right)·{}_{t_1}{}^{R}D_{t_2}^{\alpha \left(·,·\right)}{p}_{Ri}\right]}\mathrm{d}\tau +(p_{Ri}·{}_{t_1}{}^{R}D_{t_2}^{\alpha \left(·,·\right)}{q}_i\\ -H_R)·{\xi }_{R0}^0-{\displaystyle {\int }_{t_1}^t{p}_{Ri}[\gamma q_i\left(t_1\right){\xi }_{R0}^0\left(t_1,q_j\left(t_1\right),p_{Rj}\left(t_1\right)\right)\frac{\mathrm{d}}{\mathrm{d}\tau }\frac{{\left(\tau -t_1\right)}^{-\alpha \left(\tau ,t_1\right)}}{\mathsf{\Gamma }\left(1-\alpha \left(\tau ,t_1\right)\right)}}\\ -\left(1-\gamma \right)q_i\left(t_2\right){\xi }_{R0}^0\left(t_2,q_j\left(t_2\right),p_{Rj}\left(t_2\right)\right)\frac{\mathrm{d}}{\mathrm{d}\tau }\frac{{\left(t_2-\tau \right)}^{-\beta \left(t_2,\tau \right)}}{\mathsf{\Gamma }\left(1-\beta \left(t_2,\tau \right)\right)}]\mathrm{d}\tau +G_R^0\end{array}$$

(68)

Theorem 4.

For the Hamiltonian system with the RCVO (Equation (43)), if there exists a gauge function$G_{RC}^0\left(t,q_j,p_{RCj}\right)$such that the infinitesimal generators${\xi }_{RC0}^0\left(t,q_j,p_{RCj}\right)$and${\xi }_{RCi}^0\left(t,q_j,p_{RCj}\right)$satisfy

$$\begin{array}{c}{p}_{RCi}·{}_{t_1}{}^{RC}D_{t_2}^{\alpha \left(·,·\right)}\left({\xi }_{RCi}^0-{\dot{q}}_i{\xi }_{RC0}^0\right)+\left(p_{RCi}·\frac{\mathrm{d}}{\mathrm{d}t}{}_{t_1}{}^{RC}D_{t_2}^{\alpha \left(·,·\right)}{q}_i-\frac{\partial H_{RC}}{\partial t}\right){\xi }_{RC0}^0\\ -\gamma p_{RCi}{\xi }_{RC0}^0\left(t_1,q_j\left(t_1\right),p_{RCj}\left(t_1\right)\right)\frac{{\left(t-t_1\right)}^{-\alpha \left(t,t_1\right)}}{\mathsf{\Gamma }\left(1-\alpha \left(t,t_1\right)\right)}{\dot{q}}_i\left(t_1\right)+\left(p_{RCi}·{}_{t_1}{}^{RC}D_{t_2}^{\alpha \left(·,·\right)}{q}_i-H_{RC}\right)·{\dot{\xi }}_{RC0}^0\\ +\left(1-\gamma \right)p_{RCi}{\xi }_{RC0}^0\left(t_2,q_j\left(t_2\right),p_{RCj}\left(t_2\right)\right)\frac{{\left(t_2-t\right)}^{-\beta \left(t_2,t\right)}}{\mathsf{\Gamma }\left(1-\beta \left(t_2,t\right)\right)}{\dot{q}}_i\left(t_2\right)-\frac{\partial H_{RC}}{\partial q_i}{\xi }_{RCi}^0+{\dot{G}}_{RC}^0=0\end{array}$$

(69)

then there exists a conserved quantity.

$$\begin{array}{c}{I}_{HRC0}={\displaystyle {\int }_{t_1}^t\left[p_{RCi}·{}_{t_1}{}^{RC}D_{t_2}^{\alpha \left(·,·\right)}\left({\xi }_{RCi}^0-{\dot{q}}_i{\xi }_{RC0}^0\right)+\left({\xi }_{RCi}^0-{\dot{q}}_i{\xi }_{RC0}^0\right)·{}_{t_1}{}^{R}D_{t_2}^{\alpha \left(·,·\right)}{p}_{RCi}\right]}\mathrm{d}\tau +(p_{RCi}\\ ·{}_{t_1}{}^{RC}D_{t_2}^{\alpha \left(·,·\right)}{q}_i-H_{RC})·{\xi }_{RC0}^0-{\displaystyle {\int }_{t_1}^t{p}_{RCi}[\gamma {\xi }_{RC0}^0\left(t_1,q_j\left(t_1\right),p_{RCj}\left(t_1\right)\right){\dot{q}}_i\left(t_1\right)\frac{{\left(\tau -t_1\right)}^{-\alpha \left(\tau ,t_1\right)}}{\mathsf{\Gamma }\left(1-\alpha \left(\tau ,t_1\right)\right)}}\\ -\left(1-\gamma \right){\dot{q}}_i\left(t_2\right){\xi }_{RC0}^0\left(t_2,q_j\left(t_2\right),p_{RCj}\left(t_2\right)\right)\frac{{\left(t_2-\tau \right)}^{-\beta \left(t_2,\tau \right)}}{\mathsf{\Gamma }\left(1-\beta \left(t_2,\tau \right)\right)}]\mathrm{d}\tau +G_{RC}^0\end{array}$$

(70)

Remark 6.

Conserved quantities for the Lagrangian systems with variable order (i.e. Equations (44) and (46)) can be achieved as follows:

Theorem 5.

For the Lagrangian system with the CRLVO (Equation (44)), if there exists a gauge function$G_{LRL}^0\left(t,q_j,{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_j\right)$such that the infinitesimal generators${\xi }_{LRL0}^0\left(t,q_j\right)$and${\xi }_{LRLi}^0\left(t,q_j\right)$satisfy

$$\begin{array}{c}\frac{\partial L_{RL}\left(t,q_j,{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_j\right)}{\partial {}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i}·{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}\left({\xi }_{LRLi}^0-{\dot{q}}_i{\xi }_{LRL0}^0\right)+L_{RL}·{\dot{\xi }}_{RL0}^0+{\dot{G}}_{LRL}^0+\frac{\partial L_{RL}}{\partial q_i}{\xi }_{LRLi}^0+(\frac{\partial L_{RL}}{\partial t}\\ +\frac{\partial L_{RL}}{\partial {}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i}·\frac{\mathrm{d}}{\mathrm{d}t}{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i){\xi }_{LRL0}^0-\frac{\gamma q_i\left(t_1\right)\partial L_{RL}}{\partial {}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i}{\xi }_{LRL0}^0\left(t_1,q_j\left(t_1\right)\right)\frac{\mathrm{d}}{\mathrm{d}t}\frac{{\left(t-t_1\right)}^{-\alpha \left(t,t_1\right)}}{\mathsf{\Gamma }\left(1-\alpha \left(t,t_1\right)\right)}\\ +\frac{\left(1-\gamma \right)q_i\left(t_2\right)\partial L_{RL}}{\partial {}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i}{\xi }_{LRL0}^0\left(t_2,q_j\left(t_2\right)\right)\frac{\mathrm{d}}{\mathrm{d}t}\frac{{\left(t_2-t\right)}^{-\beta \left(t_2,t\right)}}{\mathsf{\Gamma }\left(1-\beta \left(t_2,t\right)\right)}=0\end{array}$$

(71)

then there exists a conserved quantity.

$$\begin{array}{c}{I}_{LRL0}={\displaystyle {\int }_{t_1}^t[\left({\xi }_{LRLi}^0-{\dot{q}}_i{\xi }_{LRL0}^0\right){}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}\frac{\partial L_{RL}}{\partial {}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i}+\frac{\partial L_{RL}}{\partial {}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i}}\\ ·{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}\left({\xi }_{LRLi}^0-{\dot{q}}_i{\xi }_{LRL0}^0\right)]\mathrm{d}\tau +L_{RL}·{\xi }_{LRL0}^0-{\displaystyle {\int }_{t_1}^t\frac{\partial L_{RL}}{\partial {}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i}[\gamma q_i\left(t_1\right){\xi }_{LRL0}^0\left(t_1,q_j\left(t_1\right)\right)}\\ ·\frac{\mathrm{d}}{\mathrm{d}\tau }\frac{{\left(\tau -t_1\right)}^{-\alpha \left(\tau ,t_1\right)}}{\mathsf{\Gamma }\left(1-\alpha \left(\tau ,t_1\right)\right)}-\left(1-\gamma \right)q_i\left(t_2\right){\xi }_{LRL0}^0\left(t_2,q_j\left(t_2\right)\right)\frac{\mathrm{d}}{\mathrm{d}\tau }\frac{{\left(t_2-\tau \right)}^{-\beta \left(t_2,\tau \right)}}{\mathsf{\Gamma }\left(1-\beta \left(t_2,\tau \right)\right)}]\mathrm{d}\tau +G_{LRL}^0\end{array}$$

(72)

Theorem 6.

For the Lagrangian system with the CCVO (Equation (46)), if there exists a gauge function$G_{LC}^0\left(t,q_j,{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_j\right)$such that the infinitesimal generators${\xi }_{LC0}^0\left(t,q_j\right)$and${\xi }_{LCi}^0\left(t,q_j\right)$satisfy

$$\begin{array}{c}\frac{\partial L_C\left(t,q_j,{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_j\right)}{\partial {}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i}·{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}\left({\xi }_{LCi}^0-{\dot{q}}_i{\xi }_{LC0}^0\right)+(\frac{\partial L_C}{\partial {}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i}·\frac{\mathrm{d}}{\mathrm{d}t}{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i\\ +\frac{\partial L_C}{\partial t}){\xi }_{LC0}^0+\left(1-\gamma \right)\frac{\partial L_C}{\partial {}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i}{\xi }_{LC0}^0\left(t_2,q_j\left(t_2\right)\right)\frac{{\left(t_2-t\right)}^{-\beta \left(t_2,t\right)}}{\mathsf{\Gamma }\left(1-\beta \left(t_2,t\right)\right)}{\dot{q}}_i\left(t_2\right)+\frac{\partial L_C}{\partial q_i}{\xi }_{LCi}^0\\ -\gamma \frac{\partial L_C}{\partial {}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i}{\xi }_{LC0}^0\left(t_1,q_j\left(t_1\right)\right)\frac{{\left(t-t_1\right)}^{-\alpha \left(t,t_1\right)}}{\mathsf{\Gamma }\left(1-\alpha \left(t,t_1\right)\right)}{\dot{q}}_i\left(t_1\right)+L_C·{\dot{\xi }}_{LC0}^0+{\dot{G}}_{LC}^0=0\end{array}$$

(73)

then there exists a conserved quantity.

$$\begin{array}{c}{I}_{LC0}={\displaystyle {\int }_{t_1}^t[\frac{\partial L_C}{\partial {}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i}·{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}\left({\xi }_{LCi}^0-{\dot{q}}_i{\xi }_{LC0}^0\right)+\left({\xi }_{LCi}^0-{\dot{q}}_i{\xi }_{LC0}^0\right)}\\ ·{}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}\frac{\partial L_C}{\partial {}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i}]\mathrm{d}\tau -{\displaystyle {\int }_{t_1}^t\frac{\partial L_C}{\partial {}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i}[\gamma {\xi }_{LC0}^0\left(t_1,q_j\left(t_1\right)\right){\dot{q}}_i\left(t_1\right)\frac{{\left(\tau -t_1\right)}^{-\alpha \left(\tau ,t_1\right)}}{\mathsf{\Gamma }\left(1-\alpha \left(\tau ,t_1\right)\right)}}\\ -\left(1-\gamma \right){\dot{q}}_i\left(t_2\right){\xi }_{LC0}^0\left(t_2,q_j\left(t_2\right)\right)\frac{{\left(t_2-\tau \right)}^{-\beta \left(t_2,\tau \right)}}{\mathsf{\Gamma }\left(1-\beta \left(t_2,\tau \right)\right)}]\mathrm{d}\tau +L_C·{\xi }_{LC0}^0+G_{LC}^0\end{array}$$

(74)

Remark 7.

When$\alpha \left(·,·\right)\to \alpha$,$\beta \left(·,·\right)\to \beta$, conserved quantities for the Hamiltonian systems with the CRL and the CC (Equations (48) and (49)) and the Lagrangian systems with the CRL and the CC (Equations (50) and (51)) can be deduced as follows:

Theorem 7.

For the Hamiltonian system with the CRL (Equation (48)), if there exists a gauge function$G_{RLF}^0\left(t,q_j,p_{RLjF}\right)$such that the infinitesimal generators${\xi }_{RL0F}^0\left(t,q_j,p_{RLjF}\right)$and${\xi }_{RLiF}^0\left(t,q_j,p_{RLjF}\right)$satisfy

$$\begin{array}{c}{p}_{RLiF}·{}^{RL}D_{\gamma }^{\alpha ,\beta }\left({\xi }_{RLiF}^0-{\dot{q}}_i{\xi }_{RL0F}^0\right)+\left(p_{RLiF}·\frac{\mathrm{d}}{\mathrm{d}t}{}^{RL}D_{\gamma }^{\alpha ,\beta }{q}_i-\frac{\partial H_{RLF}}{\partial t}\right){\xi }_{RL0F}^0+(p_{RLiF}·{}^{RL}D_{\gamma }^{\alpha ,\beta }{q}_i\\ -H_{RLF})·{\dot{\xi }}_{RL0F}^0-\frac{\partial H_{RLF}}{\partial q_i}{\xi }_{RLiF}^0+{\dot{G}}_{RLF}^0-\frac{\gamma p_{RLiF}}{\mathsf{\Gamma }\left(1-\alpha \right)}{q}_i\left(t_1\right){\xi }_{RL0F}^0\left(t_1,q_j\left(t_1\right),p_{RLjF}\left(t_1\right)\right)\frac{\mathrm{d}}{\mathrm{d}t}{\left(t-t_1\right)}^{-\alpha }\\ +\frac{\left(1-\gamma \right)p_{RLiF}}{\mathsf{\Gamma }\left(1-\beta \right)}{q}_i\left(t_2\right){\xi }_{RL0F}^0\left(t_2,q_j\left(t_2\right),p_{RLjF}\left(t_2\right)\right)\frac{\mathrm{d}}{\mathrm{d}t}{\left(t_2-t\right)}^{-\beta }=0\end{array}$$

(75)

then there exists a conserved quantity.

$$\begin{array}{c}{I}_{HRL0F}={\displaystyle {\int }_{t_1}^t\left[p_{RLiF}·{}^{RL}D_{\gamma }^{\alpha ,\beta }\left({\xi }_{RLiF}^0-{\dot{q}}_i{\xi }_{RL0F}^0\right)+\left({\xi }_{RLiF}^0-{\dot{q}}_i{\xi }_{RL0F}^0\right)·{}^{RL}D_{1-\gamma }^{\beta ,\alpha }{p}_{RLiF}\right]}\mathrm{d}\tau +(p_{RLiF}\\ ·{}^{RL}D_{\gamma }^{\alpha ,\beta }{q}_i-H_{RLF})·{\xi }_{RL0F}^0-{\displaystyle {\int }_{t_1}^t{p}_{RLiF}[\frac{\gamma q_i\left(t_1\right)}{\mathsf{\Gamma }\left(1-\alpha \right)}{\xi }_{RL0F}^0\left(t_1,q_j\left(t_1\right),p_{RLjF}\left(t_1\right)\right)\frac{\mathrm{d}}{\mathrm{d}\tau }{\left(\tau -t_1\right)}^{-\alpha }}\\ -\frac{\left(1-\gamma \right)q_i\left(t_2\right)}{\mathsf{\Gamma }\left(1-\beta \right)}{\xi }_{RL0F}^0\left(t_2,q_j\left(t_2\right),p_{RLjF}\left(t_2\right)\right)\frac{\mathrm{d}}{\mathrm{d}\tau }{\left(t_2-\tau \right)}^{-\beta }]\mathrm{d}\tau +G_{RLF}^0\end{array}$$

(76)

Theorem 8.

For the Hamiltonian system with the CC (Equation (49)), if there exists a gauge function$G_{CF}^0\left(t,q_j,p_{CjF}\right)$such that the infinitesimal generators${\xi }_{C0F}^0\left(t,q_j,p_{CjF}\right)$and${\xi }_{CiF}^0\left(t,q_j,p_{CjF}\right)$satisfy

$$\begin{array}{c}{p}_{CiF}·{}^{C}D_{\gamma }^{\alpha ,\beta }\left({\xi }_{CiF}^0-{\dot{q}}_i{\xi }_{C0F}^0\right)+\left(p_{CiF}·\frac{\mathrm{d}}{\mathrm{d}t}{}^{C}D_{\gamma }^{\alpha ,\beta }{q}_i-\frac{\partial H_{CF}}{\partial t}\right){\xi }_{C0F}^0\\ -\gamma p_{CiF}{\xi }_{C0F}^0\left(t_1,q_j\left(t_1\right),p_{CjF}\left(t_1\right)\right)\frac{{\left(t-t_1\right)}^{-\alpha }}{\mathsf{\Gamma }\left(1-\alpha \right)}{\dot{q}}_i\left(t_1\right)+\left(p_{CiF}·{}^{C}D_{\gamma }^{\alpha ,\beta }{q}_i-H_{CF}\right)·{\dot{\xi }}_{C0F}^0\\ +\left(1-\gamma \right)p_{CiF}{\xi }_{C0F}^0\left(t_2,q_j\left(t_2\right),p_{CjF}\left(t_2\right)\right)\frac{{\left(t_2-t\right)}^{-\beta }}{\mathsf{\Gamma }\left(1-\beta \right)}{\dot{q}}_i\left(t_2\right)-\frac{\partial H_{CF}}{\partial q_i}{\xi }_{CiF}^0+{\dot{G}}_{CF}^0=0\end{array}$$

(77)

then there exists a conserved quantity.

$$\begin{array}{c}{I}_{HC0F}={\displaystyle {\int }_{t_1}^t\left[p_{CiF}·{}^{C}D_{\gamma }^{\alpha ,\beta }\left({\xi }_{CiF}^0-{\dot{q}}_i{\xi }_{C0F}^0\right)+\left({\xi }_{CiF}^0-{\dot{q}}_i{\xi }_{C0F}^0\right)·{}^{RL}D_{1-\gamma }^{\beta ,\alpha }{p}_{CiF}\right]}\mathrm{d}\tau +(p_{CiF}\\ ·{}^{C}D_{\gamma }^{\alpha ,\beta }{q}_i-H_{CF})·{\xi }_{C0F}^0-{\displaystyle {\int }_{t_1}^t{p}_{CiF}[\gamma {\xi }_{C0F}^0\left(t_1,q_j\left(t_1\right),p_{CjF}\left(t_1\right)\right){\dot{q}}_i\left(t_1\right)\frac{{\left(\tau -t_1\right)}^{-\alpha }}{\mathsf{\Gamma }\left(1-\alpha \right)}}\\ -\left(1-\gamma \right){\dot{q}}_i\left(t_2\right){\xi }_{C0F}^0\left(t_2,q_j\left(t_2\right),p_{CjF}\left(t_2\right)\right)\frac{{\left(t_2-\tau \right)}^{-\beta }}{\mathsf{\Gamma }\left(1-\beta \right)}]\mathrm{d}\tau +G_{CF}^0\end{array}$$

(78)

Theorem 9.

For the Lagrangian system with the CRL (Equation (50)), if there exists a gauge function$G_{LRLF}^0\left(t,q_j,{}^{RL}D{}_{\gamma }^{\alpha ,\beta }{q}_j\right)$such that the infinitesimal generators${\xi }_{LRL0F}^0\left(t,q_j\right)$and${\xi }_{LRLiF}^0\left(t,q_j\right)$satisfy

$$\begin{array}{c}\frac{\partial L_{RLF}\left(t,q_j,{}^{RL}D_{\gamma }^{\alpha ,\beta }{q}_j\right)}{\partial {}^{RL}D_{\gamma }^{\alpha ,\beta }{q}_i}·{}^{RL}D_{\gamma }^{\alpha ,\beta }\left({\xi }_{LRLiF}^0-{\dot{q}}_i{\xi }_{LRL0F}^0\right)+L_{RLF}·{\dot{\xi }}_{RL0F}^0+\frac{\partial L_{RLF}}{\partial q_i}{\xi }_{LRLiF}^0+(\frac{\partial L_{RLF}}{\partial t}\\ +\frac{\partial L_{RLF}}{\partial {}^{RL}D_{\gamma }^{\alpha ,\beta }{q}_i}·\frac{\mathrm{d}}{\mathrm{d}t}{}^{RL}D_{\gamma }^{\alpha ,\beta }{q}_i){\xi }_{LRL0F}^0-\frac{\gamma q_i\left(t_1\right)}{\mathsf{\Gamma }\left(1-\alpha \right)}\frac{\partial L_{RLF}}{\partial {}^{RL}D_{\gamma }^{\alpha ,\beta }{q}_i}{\xi }_{LRL0F}^0\left(t_1,q_j\left(t_1\right)\right)\frac{\mathrm{d}}{\mathrm{d}t}{\left(t-t_1\right)}^{-\alpha }\\ +\frac{\left(1-\gamma \right)q_i\left(t_2\right)}{\mathsf{\Gamma }\left(1-\beta \right)}\frac{\partial L_{RLF}}{\partial {}^{RL}D_{\gamma }^{\alpha ,\beta }{q}_i}{\xi }_{LRL0F}^0\left(t_2,q_j\left(t_2\right)\right)\frac{\mathrm{d}}{\mathrm{d}t}{\left(t_2-t\right)}^{-\beta }+{\dot{G}}_{LRLF}^0=0\end{array}$$

(79)

then there exists a conserved quantity.

$$\begin{array}{c}{I}_{LRL0F}={\displaystyle {\int }_{t_1}^t[\left({\xi }_{LRLiF}^0-{\dot{q}}_i{\xi }_{LRL0F}^0\right){}^{RL}D_{1-\gamma }^{\beta ,\alpha }\frac{\partial L_{RLF}}{\partial {}^{RL}D_{\gamma }^{\alpha ,\beta }{q}_i}+\frac{\partial L_{RLF}}{\partial {}^{RL}D_{\gamma }^{\alpha ,\beta }{q}_i}}\\ ·{}^{RL}D_{\gamma }^{\alpha ,\beta }\left({\xi }_{LRLiF}^0-{\dot{q}}_i{\xi }_{LRL0F}^0\right)]\mathrm{d}\tau +L_{RLF}·{\xi }_{LRL0F}^0-{\displaystyle {\int }_{t_1}^t\frac{\partial L_{RLF}}{\partial {}^{RL}D_{\gamma }^{\alpha ,\beta }{q}_i}[\gamma q_i\left(t_1\right){\xi }_{LRL0F}^0\left(t_1,q_j\left(t_1\right)\right)}\\ ·\frac{\mathrm{d}}{\mathrm{d}\tau }\frac{{\left(\tau -t_1\right)}^{-\alpha }}{\mathsf{\Gamma }\left(1-\alpha \right)}-\left(1-\gamma \right)q_i\left(t_2\right){\xi }_{LRL0F}^0\left(t_2,q_j\left(t_2\right)\right)\frac{\mathrm{d}}{\mathrm{d}\tau }\frac{{\left(t_2-\tau \right)}^{-\beta }}{\mathsf{\Gamma }\left(1-\beta \right)}]\mathrm{d}\tau +G_{LRLF}^0\end{array}$$

(80)

Theorem 10.

For the Lagrangian system with the CC (Equation (51)), if there exists a gauge function$G_{LCF}^0\left(t,q_j,{}^{C}D_{\gamma }^{\alpha ,\beta }{q}_j\right)$such that the infinitesimal generators${\xi }_{LC0F}^0\left(t,q_j\right)$and${\xi }_{LCiF}^0\left(t,q_j\right)$satisfy

$$\begin{array}{c}\frac{\partial L_{CF}\left(t,q_j,{}^{C}D_{\gamma }^{\alpha ,\beta }{q}_j\right)}{\partial {}^{C}D_{\gamma }^{\alpha ,\beta }{q}_i}·{}^{C}D_{\gamma }^{\alpha ,\beta }\left({\xi }_{LCiF}^0-{\dot{q}}_i{\xi }_{LC0F}^0\right)+(\frac{\partial L_{CF}}{\partial {}^{C}D_{\gamma }^{\alpha ,\beta }{q}_i}·\frac{\mathrm{d}}{\mathrm{d}t}{}^{C}D_{\gamma }^{\alpha ,\beta }{q}_i\\ +\frac{\partial L_{CF}}{\partial t}){\xi }_{LC0F}^0+\left(1-\gamma \right)\frac{\partial L_{CF}}{\partial {}^{C}D_{\gamma }^{\alpha ,\beta }{q}_i}{\xi }_{LC0F}^0\left(t_2,q_j\left(t_2\right)\right)\frac{{\left(t_2-t\right)}^{-\beta }}{\mathsf{\Gamma }\left(1-\beta \right)}{\dot{q}}_i\left(t_2\right)+\frac{\partial L_{CF}}{\partial q_i}{\xi }_{LCiF}^0\\ -\gamma \frac{\partial L_{CF}}{\partial {}^{C}D_{\gamma }^{\alpha ,\beta }{q}_i}{\xi }_{LC0F}^0\left(t_1,q_j\left(t_1\right)\right)\frac{{\left(t-t_1\right)}^{-\alpha }}{\mathsf{\Gamma }\left(1-\alpha \right)}{\dot{q}}_i\left(t_1\right)+L_{CF}·{\dot{\xi }}_{LC0F}^0+{\dot{G}}_{LCF}^0=0\end{array}$$

(81)

then there exists a conserved quantity.

$$\begin{array}{c}{I}_{LC0F}={\displaystyle {\int }_{t_1}^t[\frac{\partial L_{CF}}{\partial {}^{C}D_{\gamma }^{\alpha ,\beta }{q}_i}·{}^{C}D_{\gamma }^{\alpha ,\beta }\left({\xi }_{LCiF}^0-{\dot{q}}_i{\xi }_{LC0F}^0\right)+\left({\xi }_{LCiF}^0-{\dot{q}}_i{\xi }_{LC0F}^0\right)}\\ ·{}^{RL}D_{1-\gamma }^{\beta ,\alpha }\frac{\partial L_{CF}}{\partial {}^{C}D_{\gamma }^{\alpha ,\beta }{q}_i}]\mathrm{d}\tau -{\displaystyle {\int }_{t_1}^t\frac{\partial L_{CF}}{\partial {}^{C}D_{\gamma }^{\alpha ,\beta }{q}_i}[\gamma {\xi }_{LC0F}^0\left(t_1,q_j\left(t_1\right)\right){\dot{q}}_i\left(t_1\right)\frac{{\left(\tau -t_1\right)}^{-\alpha }}{\mathsf{\Gamma }\left(1-\alpha \right)}}\\ -\left(1-\gamma \right){\dot{q}}_i\left(t_2\right){\xi }_{LC0F}^0\left(t_2,q_j\left(t_2\right)\right)\frac{{\left(t_2-\tau \right)}^{-\beta }}{\mathsf{\Gamma }\left(1-\beta \right)}]\mathrm{d}\tau +L_{CF}·{\xi }_{LC0F}^0+G_{LCF}^0\end{array}$$

(82)

Remark 8.

The results of Theorems 7–10 are consistent with those in [53].

Remark 9.

When$\alpha \left(·,·\right),\beta \left(·,·\right)\to 1$, conserved quantities for the Hamiltonian system (52) and the Lagrangian system (53) can be obtained as follows:

Theorem 11.

For the classical Hamiltonian system (Equation (52)), if there exists a gauge function$G^0\left(t,q_j,p_j\right)$such that the infinitesimal generators${\xi }_0^0\left(t,q_j,p_j\right)$and${\xi }_i^0\left(t,q_j,p_j\right)$satisfy

$$p_i{\dot{\xi }}_i^0-\frac{\partial H}{\partial t}{\xi }_0^0-H{\dot{\xi }}_0^0-\frac{\partial H}{\partial q_i}{\xi }_i^0+{\dot{G}}^0=0$$

(83)

then there exists a conserved quantity.

$$I_0=p_i{\xi }_i^0-H{\xi }_0^0+G^0$$

(84)

Theorem 12.

For the classical Lagrangian system (Equation (53)), if there exists a gauge function$G_L^0\left(t,q_j,{\dot{q}}_j\right)$such that the infinitesimal generators${\xi }_{L0}^0\left(t,q_j\right)$and${\xi }_{Li}^0\left(t,q_j\right)$satisfy

$$\frac{\partial L\left(t,q_j,{\dot{q}}_j\right)}{\partial {\dot{q}}_i}·\left({\dot{\xi }}_{Li}^0-{\dot{q}}_i{\dot{\xi }}_{L0}^0\right)+L·{\dot{\xi }}_{L0}^0+{\dot{G}}_L^0+\frac{\partial L}{\partial q_i}{\xi }_{Li}^0+\frac{\partial L}{\partial t}{\xi }_{L0}^0=0$$

(85)

then there exists a conserved quantity.

$$I_{L0}=\frac{\partial L}{\partial {\dot{q}}_i}\left({\xi }_{Li}^0-{\dot{q}}_i{\xi }_{L0}^0\right)+L·{\xi }_{L0}^0+G_L^0$$

(86)

Remark 10.

The results of Theorems 11 and 12 are consistent with those in [66].

5. Perturbation to Noether Symmetry and Adiabatic Invariants

Noether symmetry is one of the methods of finding the solutions to the differential equations of motion. However, if the mechanical system is disturbed by small forces, the original conserved quantity may also change. Perturbation to Noether symmetry and adiabatic invariants will be studied in this section because they also have a close relationship with the integration of the mechanical system.

Supposing the mechanical system is disturbed by $\epsilon W_{Ui}$, the way of disturbing the infinitesimal generators and the gauge function is set as

$$\begin{array}{c}{G}_U=G_U^0+\epsilon G_U^1+{\epsilon }^2{G}_U^2+\cdots , {\xi }_{U0}={\xi }_{U0}^0+\epsilon {\xi }_{U0}^1+{\epsilon }^2{\xi }_{U0}^2+\cdots \\ {\xi }_{Ui}={\xi }_{Ui}^0+\epsilon {\xi }_{Ui}^1+{\epsilon }^2{\xi }_{Ui}^2+\cdots , {\eta }_{Ui}={\eta }_{Ui}^0+\epsilon {\eta }_{Ui}^1+{\epsilon }^2{\eta }_{Ui}^2+\cdots \end{array}$$

(87)

where $G_U$ is the gauge function of the disturbed system, and ${\xi }_{U0}$, ${\xi }_{Ui}$, and ${\eta }_{Ui}$ are the infinitesimal generators of the disturbed system.

If a quantity $I_z$, which has a parameter $\epsilon$ and the highest power of $\epsilon$ as $z$, satisfies the rule that $\mathrm{d}{I}_z/\mathrm{d}t$ is in direct proportion to ${\epsilon }^{z+1}$, then the quantity $I_z$ is called an adiabatic invariant.

Then we have:

Theorem 13.

For the disturbed Hamiltonian system in terms of the CRLVO

$${}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i=\frac{\partial H_{RL}}{\partial p_{RLi}}, {}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}{p}_{RLi}=-\frac{\partial H_{RL}}{\partial q_i}-\epsilon W_{RLi}\left(t,q_j,p_{RLj}\right)$$

(88)

If there exists a gauge function$G_{RL}^m\left(t,q_j,p_{RLj}\right)$such that the infinitesimal generators${\xi }_{RL0}^m\left(t,q_j,p_{RLj}\right)$and${\xi }_{RLi}^m\left(t,q_j,p_{RLj}\right)$satisfy

$$\begin{array}{c}{p}_{RLi}·{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}\left({\xi }_{RLi}^m-{\dot{q}}_i{\xi }_{RL0}^m\right)+\left(p_{RLi}·\frac{\mathrm{d}}{\mathrm{d}t}{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i-\frac{\partial H_{RL}}{\partial t}\right){\xi }_{RL0}^m+(p_{RLi}·{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i\\ -H_{RL})·{\dot{\xi }}_{RL0}^m-\frac{\partial H_{RL}}{\partial q_i}{\xi }_{RLi}^m+{\dot{G}}_{RL}^m-\gamma p_{RLi}{q}_i\left(t_1\right){\xi }_{RL0}^m\left(t_1,q_j\left(t_1\right),p_{RLj}\left(t_1\right)\right)\frac{\mathrm{d}}{\mathrm{d}t}\frac{{\left(t-t_1\right)}^{-\alpha \left(t,t_1\right)}}{\mathsf{\Gamma }\left(1-\alpha \left(t,t_1\right)\right)}\\ +\left(1-\gamma \right)p_{RLi}{q}_i\left(t_2\right){\xi }_{RL0}^m\left(t_2,q_j\left(t_2\right),p_{RLj}\left(t_2\right)\right)\frac{\mathrm{d}}{\mathrm{d}t}\frac{{\left(t_2-t\right)}^{-\beta \left(t_2,t\right)}}{\mathsf{\Gamma }\left(1-\beta \left(t_2,t\right)\right)}-W_{RLi}\left({\xi }_{RLi}^{m-1}-{\dot{q}}_i{\xi }_{RL0}^{m-1}\right)=0\end{array}$$

(89)

then there exists an adiabatic invariant.

$$\begin{array}{c}{I}_{HRLz}={\displaystyle \sum _{m=0}^z{\epsilon }^m\{{\displaystyle {\int }_{t_1}^t\left[p_{RLi}·{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}\left({\xi }_{RLi}^m-{\dot{q}}_i{\xi }_{RL0}^m\right)+\left({\xi }_{RLi}^m-{\dot{q}}_i{\xi }_{RL0}^m\right)·{}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}{p}_{RLi}\right]}\mathrm{d}\tau +(p_{RLi}}\\ ·{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i-H_{RL})·{\xi }_{RL0}^m-{\displaystyle {\int }_{t_1}^t{p}_{RLi}[\gamma q_i\left(t_1\right){\xi }_{RL0}^m\left(t_1,q_j\left(t_1\right),p_{RLj}\left(t_1\right)\right)\frac{\mathrm{d}}{\mathrm{d}\tau }\frac{{\left(\tau -t_1\right)}^{-\alpha \left(\tau ,t_1\right)}}{\mathsf{\Gamma }\left(1-\alpha \left(\tau ,t_1\right)\right)}}\\ -\left(1-\gamma \right)q_i\left(t_2\right){\xi }_{RL0}^m\left(t_2,q_j\left(t_2\right),p_{RLj}\left(t_2\right)\right)\frac{\mathrm{d}}{\mathrm{d}\tau }\frac{{\left(t_2-\tau \right)}^{-\beta \left(t_2,\tau \right)}}{\mathsf{\Gamma }\left(1-\beta \left(t_2,\tau \right)\right)}]\mathrm{d}\tau +G_{RL}^m\}\end{array}$$

(90)

where${\xi }_{RLi}^{m-1}={\xi }_{RL0}^{m-1}=0$when$m=0$.

Proof. 

Using Equations (88) and (89), we have

$$\frac{\mathrm{d}}{\mathrm{d}t}{I}_{HRLz}={\displaystyle \sum _{m=0}^z{\epsilon }^m[({\dot{p}}_{RLi}·{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i+p_{RLi}·\frac{\mathrm{d}}{\mathrm{d}t}{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i-\frac{\partial H_{RL}}{\partial q_i}{\dot{q}}_i-\frac{\partial H_{RL}}{\partial p_{RLi}}{\dot{p}}_{RLi}}$$

$$-\frac{\partial H_{RL}}{\partial t})·{\xi }_{RL0}^m+\left(p_{RLi}·{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i-H_{RL}\right)·{\dot{\xi }}_{RL0}^m+p_{RLi}·{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}\left({\xi }_{RLi}^m-{\dot{q}}_i{\xi }_{RL0}^m\right)$$

$$+\left({\xi }_{RLi}^m-{\dot{q}}_i{\xi }_{RL0}^m\right)·{}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}{p}_{RLi}-\gamma p_{RLi}{q}_i\left(t_1\right){\xi }_{RL0}^m\left(t_1,q_j\left(t_1\right),p_{RLj}\left(t_1\right)\right)\frac{\mathrm{d}}{\mathrm{d}t}\frac{{\left(t-t_1\right)}^{-\alpha \left(t,t_1\right)}}{\mathsf{\Gamma }\left(1-\alpha \left(t,t_1\right)\right)}$$

$$+\left(1-\gamma \right)p_{RLi}{q}_i\left(t_2\right){\xi }_{RL0}^m\left(t_2,q_j\left(t_2\right),p_{RLj}\left(t_2\right)\right)\frac{\mathrm{d}}{\mathrm{d}t}\frac{{\left(t_2-t\right)}^{-\beta \left(t_2,t\right)}}{\mathsf{\Gamma }\left(1-\beta \left(t_2,t\right)\right)}+{\dot{G}}_{RL}^m]$$

$$={\displaystyle \sum _{m=0}^z{\epsilon }^m\left[\left({\xi }_{RLi}^m-{\dot{q}}_i{\xi }_{RL0}^m\right)·\left({}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}{p}_{RLi}+\frac{\partial H_{RL}}{\partial q_i}\right)+W_{RLi}\left({\xi }_{RLi}^{m-1}-{\dot{q}}_i{\xi }_{RL0}^{m-1}\right)\right]}$$

$$={\displaystyle \sum _{m=0}^z{\epsilon }^m\left[-\epsilon W_{RLi}\left({\xi }_{RLi}^m-{\dot{q}}_i{\xi }_{RL0}^m\right)+W_{RLi}\left({\xi }_{RLi}^{m-1}-{\dot{q}}_i{\xi }_{RL0}^{m-1}\right)\right]}=-{\epsilon }^{z+1}{W}_{RLi}\left({\xi }_{RLi}^z-{\dot{q}}_i{\xi }_{RL0}^z\right)$$

 □

Theorem 14.

For the disturbed Hamiltonian system in terms of the CCVO

$${}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i=\frac{\partial H_C}{\partial p_{Ci}}, {}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}{p}_{Ci}=-\frac{\partial H_C}{\partial q_i}-\epsilon W_{Ci}\left(t,q_j,p_{Cj}\right)$$

(91)

If there exists a gauge function$G_C^m\left(t,q_j,p_{Cj}\right)$such that the infinitesimal generators${\xi }_{C0}^m\left(t,q_j,p_{Cj}\right)$and${\xi }_{Ci}^m\left(t,q_j,p_{Cj}\right)$satisfy

$$\begin{array}{c}{p}_{Ci}·{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}\left({\xi }_{Ci}^m-{\dot{q}}_i{\xi }_{C0}^m\right)+\left(p_{Ci}·\frac{\mathrm{d}}{\mathrm{d}t}{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i-\frac{\partial H_C}{\partial t}\right){\xi }_{C0}^m-W_{Ci}\left({\xi }_{Ci}^{m-1}-{\dot{q}}_i{\xi }_{C0}^{m-1}\right)\\ -\gamma p_{Ci}{\xi }_{C0}^m\left(t_1,q_j\left(t_1\right),p_{Cj}\left(t_1\right)\right)\frac{{\left(t-t_1\right)}^{-\alpha \left(t,t_1\right)}}{\mathsf{\Gamma }\left(1-\alpha \left(t,t_1\right)\right)}{\dot{q}}_i\left(t_1\right)+\left(p_{Ci}·{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i-H_C\right)·{\dot{\xi }}_{C0}^m\\ +\left(1-\gamma \right)p_{Ci}{\xi }_{C0}^m\left(t_2,q_j\left(t_2\right),p_{Cj}\left(t_2\right)\right)\frac{{\left(t_2-t\right)}^{-\beta \left(t_2,t\right)}}{\mathsf{\Gamma }\left(1-\beta \left(t_2,t\right)\right)}{\dot{q}}_i\left(t_2\right)-\frac{\partial H_C}{\partial q_i}{\xi }_{Ci}^m+{\dot{G}}_C^m=0\end{array}$$

(92)

Then there exists an adiabatic invariant.

$$\begin{array}{c}{I}_{HCz}={\displaystyle \sum _{m=0}^z{\epsilon }^m\{{\displaystyle {\int }_{t_1}^t\left[p_{Ci}·{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}\left({\xi }_{Ci}^m-{\dot{q}}_i{\xi }_{C0}^m\right)+\left({\xi }_{Ci}^m-{\dot{q}}_i{\xi }_{C0}^m\right)·{}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}{p}_{Ci}\right]}\mathrm{d}\tau +(p_{Ci}}\\ ·{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i-H_C)·{\xi }_{C0}^m-{\displaystyle {\int }_{t_1}^t{p}_{Ci}[\gamma {\xi }_{C0}^m\left(t_1,q_j\left(t_1\right),p_{Cj}\left(t_1\right)\right){\dot{q}}_i\left(t_1\right)\frac{{\left(\tau -t_1\right)}^{-\alpha \left(\tau ,t_1\right)}}{\mathsf{\Gamma }\left(1-\alpha \left(\tau ,t_1\right)\right)}}\\ -\left(1-\gamma \right){\dot{q}}_i\left(t_2\right){\xi }_{C0}^m\left(t_2,q_j\left(t_2\right),p_{Cj}\left(t_2\right)\right)\frac{{\left(t_2-\tau \right)}^{-\beta \left(t_2,\tau \right)}}{\mathsf{\Gamma }\left(1-\beta \left(t_2,\tau \right)\right)}]\mathrm{d}\tau +G_C^m\}\end{array}$$

(93)

where${\xi }_{Ci}^{m-1}={\xi }_{C0}^{m-1}=0$when$m=0$.

Theorems 13 and 14 are the main results of this paper for the disturbed Hamiltonian systems with variable order. Based on the two theorems, adiabatic invariants for the disturbed Lagrangian systems with variable order can also been obtained.

Theorem 15.

For the disturbed Lagrangian system in terms of the CRLVO

$$\frac{\partial L_{RL}\left(t,q_j,{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_j\right)}{\partial q_i}-{}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}\frac{\partial L_{RL}}{\partial {}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i}=\epsilon W_{LRLi}\left(t,q_j,{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_j\right)$$

(94)

If there exists a gauge function$G_{LRL}^m\left(t,q_j,{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_j\right)$such that the infinitesimal generators${\xi }_{LRL0}^m\left(t,q_j\right)$and${\xi }_{LRLi}^m\left(t,q_j\right)$satisfy

$$\begin{array}{c}\frac{\partial L_{RL}}{\partial {}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i}·{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}\left({\xi }_{LRLi}^m-{\dot{q}}_i{\xi }_{LRL0}^m\right)+L_{RL}·{\dot{\xi }}_{RL0}^m+{\dot{G}}_{LRL}^m+\frac{\partial L_{RL}}{\partial q_i}{\xi }_{LRLi}^m+(\frac{\partial L_{RL}}{\partial t}\\ +\frac{\partial L_{RL}}{\partial {}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i}·\frac{\mathrm{d}}{\mathrm{d}t}{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i){\xi }_{LRL0}^m-\frac{\gamma q_i\left(t_1\right)\partial L_{RL}}{\partial {}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i}{\xi }_{LRL0}^m\left(t_1,q_j\left(t_1\right)\right)\frac{\mathrm{d}}{\mathrm{d}t}\frac{{\left(t-t_1\right)}^{-\alpha \left(t,t_1\right)}}{\mathsf{\Gamma }\left(1-\alpha \left(t,t_1\right)\right)}\\ +\frac{\left(1-\gamma \right)q_i\left(t_2\right)\partial L_{RL}}{\partial {}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i}{\xi }_{LRL0}^m\left(t_2,q_j\left(t_2\right)\right)\frac{\mathrm{d}}{\mathrm{d}t}\frac{{\left(t_2-t\right)}^{-\beta \left(t_2,t\right)}}{\mathsf{\Gamma }\left(1-\beta \left(t_2,t\right)\right)}-W_{LRLi}\left({\xi }_{LRLi}^{m-1}-{\dot{q}}_i{\xi }_{LRL0}^{m-1}\right)=0\end{array}$$

(95)

Then there exists an adiabatic invariant.

$$\begin{array}{c}{I}_{LRLz}={\displaystyle \sum _{m=0}^z{\epsilon }^m\{{\displaystyle {\int }_{t_1}^t[\left({\xi }_{LRLi}^m-{\dot{q}}_i{\xi }_{LRL0}^m\right){}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}\frac{\partial L_{RL}}{\partial {}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i}+\frac{\partial L_{RL}}{\partial {}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i}}}\\ ·{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}\left({\xi }_{LRLi}^m-{\dot{q}}_i{\xi }_{LRL0}^m\right)]\mathrm{d}\tau +L_{RL}·{\xi }_{LRL0}^m-{\displaystyle {\int }_{t_1}^t\frac{\partial L_{RL}}{\partial {}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i}[\gamma q_i\left(t_1\right){\xi }_{LRL0}^m\left(t_1,q_j\left(t_1\right)\right)}\\ ·\frac{\mathrm{d}}{\mathrm{d}\tau }\frac{{\left(\tau -t_1\right)}^{-\alpha \left(\tau ,t_1\right)}}{\mathsf{\Gamma }\left(1-\alpha \left(\tau ,t_1\right)\right)}-\left(1-\gamma \right)q_i\left(t_2\right){\xi }_{LRL0}^m\left(t_2,q_j\left(t_2\right)\right)\frac{\mathrm{d}}{\mathrm{d}\tau }\frac{{\left(t_2-\tau \right)}^{-\beta \left(t_2,\tau \right)}}{\mathsf{\Gamma }\left(1-\beta \left(t_2,\tau \right)\right)}]\mathrm{d}\tau +G_{LRL}^m\}\end{array}$$

(96)

where${\xi }_{LRLi}^{m-1}={\xi }_{LRL0}^{m-1}=0$when$m=0$.

Theorem 16.

For the disturbed Lagrangian system in terms of the CCVO

$$\frac{\partial L_C}{\partial q_i}-{}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}\frac{\partial L_C\left(t,q_j,{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_j\right)}{\partial {}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i}=\epsilon W_{LCi}\left(t,q_j,{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_j\right)$$

(97)

If there exists a gauge function$G_{LC}^m\left(t,q_j,{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_j\right)$such that the infinitesimal generators${\xi }_{LC0}^m\left(t,q_j\right)$and${\xi }_{LCi}^m\left(t,q_j\right)$satisfy

$$\begin{array}{c}\frac{\partial L_C}{\partial {}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i}·{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}\left({\xi }_{LCi}^m-{\dot{q}}_i{\xi }_{LC0}^m\right)+\left(\frac{\partial L_C}{\partial {}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i}·\frac{\mathrm{d}}{\mathrm{d}t}{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i+\frac{\partial L_C}{\partial t}\right){\xi }_{LC0}^m\\ +\left(1-\gamma \right)\frac{\partial L_C}{\partial {}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i}{\xi }_{LC0}^m\left(t_2,q_j\left(t_2\right)\right)\frac{{\left(t_2-t\right)}^{-\beta \left(t_2,t\right)}}{\mathsf{\Gamma }\left(1-\beta \left(t_2,t\right)\right)}{\dot{q}}_i\left(t_2\right)+\frac{\partial L_C}{\partial q_i}{\xi }_{LCi}^m+L_C·{\dot{\xi }}_{LC0}^m+{\dot{G}}_{LC}^m\\ -\gamma \frac{\partial L_C}{\partial {}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i}{\xi }_{LC0}^m\left(t_1,q_j\left(t_1\right)\right)\frac{{\left(t-t_1\right)}^{-\alpha \left(t,t_1\right)}}{\mathsf{\Gamma }\left(1-\alpha \left(t,t_1\right)\right)}{\dot{q}}_i\left(t_1\right)-W_{LCi}\left({\xi }_{LCi}^{m-1}-{\dot{q}}_i{\xi }_{LC0}^{m-1}\right)=0\end{array}$$

(98)

Then there exists an adiabatic invariant.

$$\begin{array}{c}{I}_{LCz}={\displaystyle \sum _{m=0}^z{\epsilon }^m\{{\displaystyle {\int }_{t_1}^t[\frac{\partial L_C}{\partial {}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i}·{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}\left({\xi }_{LCi}^m-{\dot{q}}_i{\xi }_{LC0}^m\right)+\left({\xi }_{LCi}^m-{\dot{q}}_i{\xi }_{LC0}^m\right)}}\\ ·{}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}\frac{\partial L_C}{\partial {}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i}]\mathrm{d}\tau -{\displaystyle {\int }_{t_1}^t\frac{\partial L_C}{\partial {}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_i}[\gamma {\xi }_{LC0}^m\left(t_1,q_j\left(t_1\right)\right){\dot{q}}_i\left(t_1\right)\frac{{\left(\tau -t_1\right)}^{-\alpha \left(\tau ,t_1\right)}}{\mathsf{\Gamma }\left(1-\alpha \left(\tau ,t_1\right)\right)}}\\ -\left(1-\gamma \right){\dot{q}}_i\left(t_2\right){\xi }_{LC0}^m\left(t_2,q_j\left(t_2\right)\right)\frac{{\left(t_2-\tau \right)}^{-\beta \left(t_2,\tau \right)}}{\mathsf{\Gamma }\left(1-\beta \left(t_2,\tau \right)\right)}]\mathrm{d}\tau +L_C·{\xi }_{LC0}^m+G_{LC}^m\}\end{array}$$

(99)

where${\xi }_{LCi}^{m-1}={\xi }_{LC0}^{m-1}=0$when$m=0$.

Moreover, adiabatic invariants for the disturbed Hamiltonian systems with the RVO and the RCVO, the disturbed Hamiltonian systems with the CRL and the CC, the disturbed Lagrangian systems with the CRL and the CC, the disturbed classical Hamiltonian system, and the disturbed classical Lagrangian system can also be deduced. We only present them in the forms of remarks rather than describe them in detail.

Remark 11.

Let$\beta \left(·,·\right)=\alpha \left(·,·\right)$,$\gamma =\frac{1}{2}$, so adiabatic invariants for the disturbed Hamiltonian systems in terms of the RVO and the RCVO can be achieved.

Remark 12.

Let$\alpha \left(·,·\right)\to \alpha$,$\beta \left(·,·\right)\to \beta$, so adiabatic invariants for the disturbed Hamiltonian systems with the CRL and the CC as well as the disturbed Lagrangian systems with the CRL and the CC can be deduced.

Remark 13.

Let$\alpha \left(·,·\right),\beta \left(·,·\right)\to 1$, so adiabatic invariants for the disturbed classical Hamiltonian system and the disturbed classical Lagrangian system can be deduced.

6. Applications

In this section, we present two applications to illustrate the Methods and Results.

Application 1.

The differential equation of the Lotka biochemical oscillator [67] of variable order is

$${}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{x}_1={\alpha }_1+{\beta }_1\exp x_2, {}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}{x}_2={\alpha }_2+{\beta }_2\exp x_1$$

(100)

where${\alpha }_1$,${\alpha }_2$,${\beta }_1$, and${\beta }_2$ are constants.

Let $x_1=q$, $x_2=p_{RL}$, so the Hamiltonian can be expressed as

$$H_{RL}={\alpha }_1{p}_{RL}-{\alpha }_{2}q+{\beta }_1\exp p_{RL}-{\beta }_2\exp q$$

(101)

From Equation (61), we have

$$\begin{array}{c}{p}_{RL}·{\xi }_{RL0}^0·\frac{\mathrm{d}}{\mathrm{d}t}{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}q+\left(p_{RL}·{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}q-{\alpha }_1{p}_{RL}+{\alpha }_{2}q-{\beta }_1\exp p_{RL}+{\beta }_2\exp q\right)·{\dot{\xi }}_{RL0}^0\\ -p_{RL}·\gamma q\left(t_1\right){\xi }_{RL0}^0\left(t_1,q_j\left(t_1\right),p_{RLj}\left(t_1\right)\right)\frac{\mathrm{d}}{\mathrm{d}t}\frac{{\left(t-t_1\right)}^{-\alpha \left(t,t_1\right)}}{\mathsf{\Gamma }\left(1-\alpha \left(t,t_1\right)\right)}+p_{RL}·{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}\left({\xi }_{RL}^0-\dot{q}{\xi }_{RL0}^0\right)+{\dot{G}}_{RL}^0\\ +\left({\alpha }_2+{\beta }_2\exp q\right){\xi }_{RL}^0+p_{RL}·\left(1-\gamma \right)q\left(t_2\right){\xi }_{RL0}^0\left(t_2,q_j\left(t_2\right),p_{RL}{}_j\left(t_2\right)\right)\frac{\mathrm{d}}{\mathrm{d}t}\frac{{\left(t_2-t\right)}^{-\beta \left(t_2,t\right)}}{\mathsf{\Gamma }\left(1-\beta \left(t_2,t\right)\right)}=0\end{array}$$

(102)

It is obvious that

$${\xi }_{RL0}^0=-1, {\xi }_{RL}^0=0, G_{RL}^0=0$$

(103)

satisfies Equation (102), where [57]

$$\frac{\mathrm{d}}{\mathrm{d}t}{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}q={}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}\dot{q}+\gamma q\left(t_1\right)\frac{\mathrm{d}}{\mathrm{d}t}\frac{{\left(t-t_1\right)}^{-\alpha \left(t,t_1\right)}}{\mathsf{\Gamma }\left(1-\alpha \left(t,t_1\right)\right)}-\left(1-\gamma \right)q\left(t_2\right)\frac{\mathrm{d}}{\mathrm{d}t}\frac{{\left(t_2-t\right)}^{-\beta \left(t_2,t\right)}}{\mathsf{\Gamma }\left(1-\beta \left(t_2,t\right)\right)}$$

(104)

Then, a conserved quantity can be obtained from Theorem 1 as

$$I_{HRL0}={\displaystyle {\int }_{t_1}^t\left(p_{RL}·\frac{\mathrm{d}}{\mathrm{d}\tau }{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}q+\dot{q}{}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}{p}_{RL}\right)}\mathrm{d}\tau -\left(p_{RL}·{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}q-H_{RL}\right)$$

(105)

When the Lotka biochemical oscillator model is disturbed as

$${}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}q={\alpha }_1+{\beta }_1\exp p_{RL}, {}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}{p}_{RL}={\alpha }_2+{\beta }_2\exp q-\epsilon \left(2q+1\right)$$

(106)

We can get the following solution to Equation (89)

$${\xi }_{RL0}^1=1, {\xi }_{RL}^1=0, G_{RL}^1=q^2+q$$

(107)

Then, the first order adiabatic invariant can be obtained from Theorem 13 as

$$\begin{array}{c}{I}_{HRL1}={\displaystyle {\int }_{t_1}^t\left(p_{RL}·\frac{\mathrm{d}}{\mathrm{d}\tau }{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}q+\dot{q}{}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}{p}_{RL}\right)}\mathrm{d}\tau -(p_{RL}·{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}q-{\alpha }_1{p}_{RL}+{\alpha }_{2}q\\ -{\beta }_1\exp p_{RL}+{\beta }_2\exp q)+\epsilon [q^2+q+\left(p_{RL}·{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}q-{\alpha }_1{p}_{RL}+{\alpha }_{2}q-{\beta }_1\exp p_{RL}+{\beta }_2\exp q\right)\\ -{\displaystyle {\int }_{t_1}^t\left(p_{RL}·\frac{\mathrm{d}}{\mathrm{d}\tau }{}^{RL}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}q+\dot{q}{}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}{p}_{RL}\right)}\mathrm{d}\tau ]\end{array}$$

(108)

When $\alpha$, $\beta \to 1$, the classical differential equation of the Lotka biochemical oscillator, the classical Hamiltonian, the classical conserved quantity, and the classical adiabatic invariant can be obtained, which are consistent with the results in [66].

Application 2.

The Lagrangian of the two-dimensional isotropic harmonic oscillator of variable order is

$$L_C=\frac{1}{2}m\left[{\left({}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_1\right)}^2+{\left({}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_2\right)}^2\right]-\frac{1}{2}k\left[{\left(q_1\right)}^2+{\left(q_2\right)}^2\right]$$

(109)

where$k$and$m$are constants.

Firstly, the generalized momentum and the Hamiltonian of variable order can be obtained from Equation (37):

$$p_{C1}=\frac{\partial L_C\left(t,q_j,{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_j\right)}{\partial {}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_1}=m{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_1$$

$$p_{C2}=\frac{\partial L_C\left(t,q_j,{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_j\right)}{\partial {}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_2}=m{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_2$$

$$\begin{array}{c}{H}_C=p_{C1}·{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_1+p_{C2}·{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_2-L_C\left(t,q_j,{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_j\right)\\ =p_{C1}·\frac{p_{C1}}{m}+p_{C2}·\frac{p_{C2}}{m}-\frac{1}{2}m\left[{\left(\frac{p_{C1}}{m}\right)}^2+{\left(\frac{p_{C2}}{m}\right)}^2\right]+\frac{1}{2}k\left[{\left(q_1\right)}^2+{\left(q_2\right)}^2\right]\\ =\frac{{\left(p_{C1}\right)}^2}{2m}+\frac{{\left(p_{C2}\right)}^2}{2m}+\frac{1}{2}k\left[{\left(q_1\right)}^2+{\left(q_2\right)}^2\right]\end{array}$$

(110)

From Equation (41), we have

$$\begin{array}{c}{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_1=\frac{p_{C1}}{m}, {}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_2=\frac{p_{C2}}{m}\\ {}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}{p}_{C1}=-k{q}_1, {}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}{p}_{C2}=-k{q}_2\end{array}$$

(111)

And from Equation (65), we find that

$${\xi }_{C0}^0=-1, {\xi }_{C1}^0={\xi }_{C2}^0=0, G_C^0=0$$

(112)

is a proper solution. A conserved quantity can be obtained from Theorem 2 as

$$\begin{array}{c}{I}_{HC0}=-\left[p_{C1}·{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_1+p_{C2}·{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_2-\frac{{\left(p_{C1}\right)}^2}{2m}-\frac{{\left(p_{C2}\right)}^2}{2m}-\frac{1}{2}k\left({\left(q_1\right)}^2+{\left(q_2\right)}^2\right)\right]\\ +{\displaystyle {\int }_{t_1}^t(p_{C1}·\frac{\mathrm{d}}{\mathrm{d}\tau }{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_1+p_{C2}·\frac{\mathrm{d}}{\mathrm{d}\tau }{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_2+{\dot{q}}_1·{}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}{p}_{C1}}\\ +{\dot{q}}_2·{}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}{p}_{C2})\mathrm{d}\tau \end{array}$$

(113)

When the isotropic harmonic oscillator model of variable order is disturbed as

$$\begin{array}{c}{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_1=\frac{p_{C1}}{m}, {}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_2=\frac{p_{C2}}{m}, {}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}{p}_{C1}=-k{q}_1-\epsilon q_2\\ {}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}{p}_{C2}=-k{q}_2-\epsilon q_1\end{array}$$

(114)

we obtain

$${\xi }_{C0}^1=1, {\xi }_{C1}^1={\xi }_{C2}^1=0, G_C^1=q_1{q}_2$$

(115)

which satisfies Equation (92). Then, the first order adiabatic invariant can be obtained from Theorem 14 as

$$\begin{array}{c}{I}_{HC1}=-\left[p_{C1}·{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_1+p_{C2}·{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_2-\frac{{\left(p_{C1}\right)}^2}{2m}-\frac{{\left(p_{C2}\right)}^2}{2m}-\frac{1}{2}k\left({\left(q_1\right)}^2+{\left(q_2\right)}^2\right)\right]\\ +{\displaystyle {\int }_{t_1}^t(p_{C1}·\frac{\mathrm{d}}{\mathrm{d}\tau }{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_1+p_{C2}·\frac{\mathrm{d}}{\mathrm{d}\tau }{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_2+{\dot{q}}_1·{}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}{p}_{C1}+{\dot{q}}_2}\\ ·{}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}{p}_{C2})\mathrm{d}\tau -\epsilon \{{\displaystyle {\int }_{t_1}^t(p_{C1}·\frac{\mathrm{d}}{\mathrm{d}\tau }{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_1+p_{C2}·\frac{\mathrm{d}}{\mathrm{d}\tau }{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_2+{\dot{q}}_1}\\ ·{}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}{p}_{C1}+{\dot{q}}_2·{}^{RL}D_{1-\gamma }^{\beta \left(·,·\right),\alpha \left(·,·\right)}{p}_{C2})\mathrm{d}\tau -[p_{C1}·{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_1+p_{C2}·{}^{C}D_{\gamma }^{\alpha \left(·,·\right),\beta \left(·,·\right)}{q}_2\\ -\frac{{\left(p_{C1}\right)}^2}{2m}-\frac{{\left(p_{C2}\right)}^2}{2m}-\frac{1}{2}k\left({\left(q_1\right)}^2+{\left(q_2\right)}^2\right)]-q_1{q}_2\}\end{array}$$

(116)

When $\alpha$, $\beta \to 1$, the classical Lagrangian, the classical Hamiltonian of the two-dimensional isotropic harmonic oscillator, the classical conserved quantity, and the classical adiabatic invariant can be obtained, which is consistent with the results in [66].

7. Conclusions

Equations of motion, Noether symmetry and conserved quantities, and perturbation to Noether symmetry and adiabatic invariants are investigated here. Hamilton equations with the CRLVO (Equation (36)), CCVO (Equation (41)), RVO (Equation (42)), RCVO (Equation (43)), CRL (Equation (48)), and CC (Equation (49)) are established. Then, Noether symmetry and conserved quantities are studied for those six Hamiltonian systems (Theorem 1, Theorem 2, Theorem 3, Theorem 4, Theorem 7, and Theorem 8). Meanwhile, Lagrange equations with the CRLVO (Equation (44)), CCVO (Equation (46)), CRL (Equation (50)), and CC (Equation (51)) are presented. Then Noether symmetry and conserved quantities are investigated for those four Lagrangian systems (Theorem 5, Theorem 6, Theorem 9, and Theorem 10). As for perturbation to Noether symmetry and adiabatic invariants, only adiabatic invariants with the CRLVO and CCVO for the Hamiltonian systems and the Lagrangian systems are studied and described in detail (Theorem 13, Theorem 14, Theorem 15, and Theorem 16).

Among the results obtained in this paper, Equations (36), (41)–(44), and (46), Theorems 1–6, and Theorems 13–16 are new. Equations (48)–(53) and Theorems 7–12, which are deduced from the main results of this paper, are consistent with the existing results.

It is generally known that there are three kinds of symmetry in analytical mechanics, i.e., Noether symmetry, Lie symmetry, and Mei symmetry. Noether symmetry plays an important role in finding the solutions to the differential equations of motion for mechanical systems because conserved quantities can be deduced from it. In fact, Lie symmetry and Mei symmetry are also able to deduce conserved quantities under certain conditions. The Lie symmetry means the invariance of the differential equations of motion under the infinitesimal transformations of time and coordinates. The Mei symmetry means the invariance under which the transformed dynamical functions still satisfy the original differential equations of motion. The conserved quantities deduced directly by the Lie symmetry and the Mei symmetry are called the Hojman conserved quantity and the Mei conserved quantity, respectively. However, only Noether symmetry is considered here, so Lie symmetry, Mei symmetry, perturbation to Lie symmetry, and perturbation to Mei symmetry are important and valuable aspects waiting to be studied.