Full Text PDF

Vol.
122 (2012)
ACTA PHYSICA POLONICA A
No. 1
Perturbation to Noether Symmetry and Noether
Adiabatic Invariants of Discrete Dierence
Variational Hamilton Systems
M.-J. Zhanga,b,∗ and J.-H. Fangc
Purple Mountain Observatory, Chinese Academy of Sciences, Nanjing 210008, China
Key Laboratory of Space Object and Debris Observation, PMO, CAS, Nanjing 210008, China
a
b
c
College of Science, China University of Petroleum, Qingdao 266555, China
(Received April 15, 2011; in nal form March 2, 2012)
The perturbation to the Noether symmetry and the Noether adiabatic invariants of discrete dierence
variational Hamilton systems are investigated. The discrete the Noether exact invariant induced directly by the
the Noether symmetry of the system without perturbation is given. The concept of discrete high-order adiabatic
invariant is presented, the criterion of the perturbation to the Noether symmetry is established, and the discrete
the Noether adiabatic invariant induced directly by the perturbation to the Noether symmetry is obtained. Lastly,
an example is discussed to illustrate the application of the results.
PACS: 02.20.Sv, 11.30.−j, 45.05.+x
1. Introduction
adiabatic invariants play a very important role in the research on the quasi-integrability of mechanical systems.
A classical adiabatic invariant is a certain physical quantity that changes more slowly than some slowly-varying
parameter of the system [17]. In fact, the parameter
varying very slowly is equivalent to the action of a small
disturbance. At present, studies in this eld have become
very active, and many important results have been obtained [1823]. But all these studies have focused on perturbation to symmetries and adiabatic invariants of the
continuous mechanical systems. Recently, Zhang et al.
[24] presented the concept of discrete high-order adiabatic invariant, and studied the perturbation to the the
Noether symmetry and the Noether adiabatic invariant
of the general discrete holonomic system. Wang and Zhu
[25] further discussed perturbation to symmetry and adiabatic invariants of general discrete holonomic dynamical systems on a uniform lattice. However, perturbation
to symmetries and adiabatic invariants of the discrete
Hamilton systems has never been studied so far.
In this paper, based on the concept of discrete
high-order adiabatic invariant, the perturbation to the
Noether symmetry and the Noether adiabatic invariants
of discrete dierence variational Hamilton systems are
studied. The discrete the Noether exact invariant induced directly by the the Noether symmetry of the system without perturbation is given. The criterion of the
perturbation to the Noether symmetry is established, and
the discrete the Noether adiabatic invariant induced directly by the perturbation to the Noether symmetry is
obtained. Meanwhile, an example is discussed to illustrate the application of the results.
The research on symmetries and conserved quantities of a mechanical system possesses important mathematical and physical signicance, and is also an important development direction in analytical mechanics. In
Refs. [1, 2] symmetries and conserved quantities of constrained mechanical systems have been thoroughly investigated. Moreover, symmetries and conserved quantities of more general systems (optimal control problems)
are also studied in depth [35]. Recently, researches
on symmetries and conserved quantities have been extended to a general time scale (including, as particular
cases, both discrete and continuous settings). Levi et al.
[68] rst extended the Lie symmetry to discrete systems; Dorodnitsyn [9] adapted the Noether's theory to
discrete Lagrangian system; Torres et al. [1012] studied
the Noether symmetry theorems for an arbitrary time
scale; Shi et al. [13, 14] investigated Lie symmetry and the
Noether conserved quantities (or exact invariants) of discrete non-conservative mechanical systems and discrete
dierence variational Hamilton systems without perturbations; based upon the property of the discrete models entirely inheriting the symmetry of the continuous
systems, Fu et al. [15] presented the Hojman conserved
quantities and Lie symmetries of discrete mechanico-electrical coupling systems by the Lie groups of transformations of continuous systems.
As we know, even a tiny disturbance acting on the mechanical systems, that we can call a perturbation, may
inuence the original symmetries and conserved quantities of mechanical systems. Pioneered in this area, Burgers [16] proposed adiabatic invariants for a special kind of
the Hamilton systems. Perturbation to symmetries and
∗ corresponding author; e-mail:
2. the Noether symmetry and discrete the
Noether exact invariant
For brevity of notation, we consider a one-dimensional discrete dierence variational Hamil-
[email protected]
(25)
M.-J. Zhang, J.-H. Fang
26
ton system.
The Hamiltonian of the system is
Hd (tk , tk+1 , qk , qk+1 , pk , pk+1 ) (k = 0, 1, . . . , N − 1), tk ,
qk , and pk are discrete time, discrete generalized coordinates, and discrete generalized momenta, respectively.
The discrete equations of motion of the system can be
written as [14]:
1
(pk−1 − pk+1 ) − D3 Hd (φk )(tk+1 − tk )
2
− D4 Hd (φk−1 )(tk − tk−1 ) = 0,
1
2
(1)
(qk+1 − qk−1 ) − D5 Hd (φk )(tk+1 − tk )
− D6 Hd (φk−1 )(tk − tk−1 ) = 0,
(2)
Hd (φk ) − Hd (φk−1 ) − D1 Hd (φk )(tk+1 − tk )
− D2 Hd (φk−1 )(tk − tk−1 ) = 0,
(3)
where Dj is the partial derivative of the discrete function with respect to the argument j and
φk = (tk , tk+1 , qk , qk+1 , pk , pk+1 ) represents the discrete
sequence.
We introduce the innitesimal transformations with
respect to discrete time tk , discrete generalized coordinates qk , and discrete generalized momenta pk as
t∗k = tk + ∆tk = tk + ετk0 (tk , qk , pk ),
qk∗ = qk + ∆qk = qk + εξk0 (tk , qk , pk ),
p∗k = pk + ∆pk = qk + εηk0 (tk , qk , pk ),
(4)
τk0 ,
where ε is an innitesimal group parameter,
and
ηk0 are the discrete innitesimal generators. The vector
eld of generators is
∂
∂
∂
(0)
X0,d = τk0
+ ξk0
+ ηk0
,
(5)
∂tk
∂qk
∂pk
which can be prolonged to the two-point scheme
∂
∂
(1)
(0)
0
0
X0,d = X0,d + τk+1
+ ξk+1
∂tk+1
∂qk+1
∂
0
+ ηk+1
.
(6)
∂pk+1
ξk0
The recursive and derivative operators of discrete systems for any discrete variable or function are represented
as
R± f (zk ) = f (zk±1 ),
(7)
R+ f (zk ) − f (zk )
.
(8)
tk+1 − tk
The the Noether symmetry is an invariance of the
Hamiltonian action functional under the innitesimal
transformations.
Then the requirement of the the
Noether symmetry of the discrete dierence variational
Hamilton system gives
Dd f (zk ) =
(1)
Hd (φk )Dd (τk0 ) + X0,d [Hd (φk )] − 21 (pk+1 + pk )Dd (ξk0 )
+ 12 (qk+1 + qk )Dd (ηk0 ) + Dd (G0N k ) = 0,
where
(9)
is a discrete gauge function.
Criterion 1. If a discrete gauge function G0N k exists
such that the innitesimal generators τk0 , ξk0 and ηk0 satisfy
G0N k
=
G0N k (tk , qk , pk )
Eq. (9), the invariance is the Noether symmetry of the
discrete dierence variational Hamilton system (1), (2)
and (3).
Equation (9) is called the discrete the Noether identity
of the discrete dierence variational Hamilton system (1),
(2) and (3).
Theorem 1 [14]. For the discrete dierence varia-
tional Hamilton system (1), (2) and (3), if the innitesimal generators τk0 , ξk0 and ηk0 of the the Noether symmetry
and discrete gauge function G0N k satisfy the discrete the
Noether identity (9), the the Noether symmetry of the
system can induce the discrete the Noether exact invariant
IN0,d = τk0 (tk − tk−1 )D2 [R− Hd (φk )]
+ ξk0 (tk
(10)
− tk−1 )D4 [R− Hd (φk )]
+ ηk0 (tk − tk−1 )D6 [R− Hd (φk )] + τk0 R− Hd (φk )
− 21 (pk−1 + pk )ξk0 + 12 (qk−1 + qk )ηk0 + G0N k = const.
3. Perturbation to the Noether symmetry and
discrete the Noether adiabatic invariant
Small force acting on the discrete dierence variational
Hamilton system may result in a small change in its symmetries, i.e. perturbation to symmetries. The conserved
quantity associated with the symmetries, under a corresponding change, is an adiabatic invariant. In analytical mechanics, we study perturbation to symmetries and
adiabatic invariants of mechanical systems based on the
concept of high-order adiabatic invariant. According to
the concept of adiabatic invariant [17], for the discrete
dierence variational Hamilton system, we have
Denition. If Lz,d (tk , tk+1 , qk , qk+1 , pk , pk+1 , ε) is
a physical quantity including ε in which the highest
power is z in a discrete dierence variational Hamil-
ton system, and its derivative with respect to discrete
time tk is directly proportional to εz+1 , Iz,d is called
a z -th order adiabatic invariant of the discrete dierence variational Hamilton system. Denition. If
Lz,d (tk , tk+1 , qk , qk+1 , pk , pk+1 , ε) is a physical quantity
including ε in which the highest power is z in a discrete
dierence variational Hamilton system, and its derivative
with respect to discrete time tk is directly proportional to
εz+1 , Iz,d is called a z -th order adiabatic invariant of the
discrete dierence variational Hamilton system.
Suppose the discrete dierence variational Hamilton
system (1), (2) and (3) is perturbed by small quantity εWd (tk , tk+1 , qk , qk+1 , pk , pk+1 ) = εWd (φk ), then we
have the following total variation of corresponding discrete Hamiltonian action functional:
N
−1
∑
∆Sd +
εWd (φk )(tk+1 − tk )∆qk
−
k=0
N
−1
∑
εWd (φk )(qk+1 − qk )∆tk = 0,
k=0
where ∆ is a total variation symbol, and
(11)
Perturbation to Noether Symmetry and Noether Adiabatic Invariants . . .
Sd =
N
−1
∑
1
2
(pk+1 + pk )(qk+1 − qk )
(12)
k=0
−
N
−1
∑
Hd (tk , tk+1 , qk , qk+1 , pk , pk+1 )(tk+1 − tk ).
From (11), using the similar deduction in [14], we obtain the discrete equations of motion of the discrete difference variational Hamilton system perturbed by small
quantity εWd (φk ):
1
2
(pk−1 − pk+1 ) − D3 Hd (φk )(tk+1 − tk )
− D4 Hd (φk−1 )(tk − tk−1 )
+ εWd (φk )(tk+1 − tk ) = 0,
1
2
(14)
− D2 Hd (φk−1 )(tk − tk−1 )
− εWd (φk )(qk+1 − qk ) = 0.
(15)
Because of the action of εWd (φk ), the original symmetries and invariants of the system may vary. Assume that
the variation is a small perturbation based on the symmetrical transformations of the system without perturbation, and τk (tk , qk , pk ), ξk (tk , qk , pk ) and ηk (tk , qk , pk )
express the generators of innitesimal transformations after being perturbed, then
τk = τk0 + ετk1 + ε2 τk2 + . . . ,
ξk =
ηk =
ηk0
+
εξk1
+
εηk1
+
ε2 ξk2
+ ...,
+
ε2 ηk2
+ ...
(16)
The innitesimal transformations become
t∗k = tk + ∆tk = tk + ετk (tk , qk , pk ),
qk∗ = qk + ∆qk = qk + εξk (tk , qk , pk ),
p∗k = pk + ∆pk = pk + εηk (tk , qk , pk ).
(1)
Xd = εm Xm,d (m = 0, 1, . . . , z),
(21)
where
∂
∂
∂
∂
(1)
m
+ ξkm
+ ηkm
+ τk+1
Xm,d = τkm
∂tk
∂qk
∂pk
∂tk+1
m
+ ξk+1
∂
∂
m
+ ηk+1
.
∂qk+1
∂pk+1
(22)
Substituting (16) into (18), noticing (19)(22), and
making the coecients of εm equal, we obtain
(23)
Hd (φk )Dd (τkm ) + Xm,d [Hd (φk )]
− (pk+1 +
+ (qk+1 +
[
]
− Wd (φk ) ξkm−1 − Dd (qk )τkm−1 + Dd (Gm
N k ) = 0,
1
2
pk )Dd (ξkm )
1
2
qk )Dd (ηkm )
when m = 0, we note that τkm−1 = ξkm−1 = ηkm−1 = 0
holds, then Eq. (23) turns into Eq. (9).
Hd (φk ) − Hd (φk−1 ) − D1 Hd (φk )(tk+1 − tk )
ξk0
Substituting (16) into (19), we have
(1)
(13)
(qk+1 − qk−1 ) − D5 Hd (φk )(tk+1 − tk )
− D6 Hd (φk−1 )(tk − tk−1 ) = 0,
(20)
GN k = G0N k + εG1N k + ε2 G2N k + . . .
(1)
k=0
27
(17)
According to the the Noether symmetry theory, the
the Noether identity of the discrete dierence variational
Hamilton system perturbed by small quantity εWd (φk )
becomes
Criterion 2. For the discrete dierence variational
Hamilton system (1), (2) and (3) perturbed by small
quantity εWd (φk ), if the discrete gauge function Gm
N k exists such that the innitesimal generators τkm , ξkm and ηkm
satisfy (23), the corresponding variety of the the Noether
symmetry is called the perturbation to the the Noether
symmetry.
Equation (23) is the determining equation of the perturbation to the the Noether symmetry of the discrete
dierence variational Hamilton system (1), (2) and (3).
Theorem 2. For the discrete dierence variational Hamilton system (1), (2) and (3) perturbed by
small quantity εWd (φk ), if the innitesimal generators
τkm , ξkm and ηkm , and discrete gauge function Gm
N k satisfy
the determining of Eq. (23), the perturbation to the the
Noether symmetry of the system can induce a z -th order
discrete adiabatic invariant
{
INz,d = εm τkm (tk − tk−1 )D2 [R− Hd (φk )]
+ ξkm (tk − tk−1 )D4 [R− Hd (φk )]
+ ηkm (tk − tk−1 )D6 [R− Hd (φk )] + τkm R− Hd (φk )
}
m
m
m
1
1
− 2 (pk−1 + pk )ξk + 2 (qk−1 + qk )ηk + GN k
(1)
Hd (φk )Dd (τk ) + Xd [Hd (φk )] − 21 (pk+1 + pk )Dd (ξk )
+ 21 (qk+1 + qk )Dd (ηk ) − εWd (φk )[ξk − Dd (qk )τk ]
+ Dd (GN k ) = 0,
where
∂
∂
∂
∂
(1)
+ ξk
+ ηk
+ τk+1
Xd = τk
∂tk
∂qk
∂pk
∂tk+1
+ ξk+1
∂
∂
+ ηk+1
,
∂qk+1
∂pk+1
(18)
(24)
(m = 0, 1, . . . , z).
Proof: Using representations (7), (8) and (22), the
expansion of (23) gives
(1)
Hd (φk )Dd (τkm ) + Xm,d [Hd (φk )]
− 12 (pk+1 + pk )Dd (ξkm ) + 12 (qk+1 + qk )Dd (ηkm )
(19)
and GN k = GN k (tk , qk , pk ) is a gauge function. After
being perturbed, the gauge function comes into
[
]
− Wd (φk ) ξkm−1 − Dd (qk )τkm−1 + Dd (Gm
Nk)
= Hd (φk )
m
τk+1
− τkm
∂Hd (φk )
m ∂Hd (φk )
+ τkm
+ τk+1
tk+1 − tk
∂tk
∂tk+1
M.-J. Zhang, J.-H. Fang
28
+ ξkm
∂Hd (φk )
∂Hd (φk )
m ∂Hd (φk )
+ ξk+1
+ ηkm
∂qk
∂qk+1
∂pk
m
+ ηk+1
ξ m − ξkm
∂Hd (φk ) 1
− 2 (pk+1 + pk ) k+1
∂pk+1
tk+1 − tk
+ 12 (qk+1 + qk )
m
ηk+1
− ηkm
tk+1 − tk
]
[
− Wd (φk ) ξkm−1 − Dd (qk )τkm−1 + Dd (Gm
Nk)
[
∂Hd (φk ) tk − tk−1 ∂Hd (φk−1 )
+
∂pk
tk+1 − tk
∂pk
]
qk+1
qk−1
− 21
+ 21
tk+1 − tk
tk+1 − tk
[
∂Hd (φk−1 )
+ Dd τkm (tk − tk−1 )
∂tk
+ηkm
+ ξkm (tk − tk−1 )
∂Hd (φk−1 )
∂qk
∂Hd (φk−1 )
∂pk
= τkm
∂Hd (φk )
tk − tk−1 ∂Hd (φk−1 )
+ τkm
∂tk
tk+1 − tk
∂tk
+ ηkm (tk − tk−1 )
+ τkm
Hd (φk−1 )
Hd (φk )
− τkm
tk+1 − tk
tk+1 − tk
+ τkm Hd (φk−1 ) − 21 (pk−1 + pk )ξkm
+ ξkm
∂Hd (φk )
tk − tk−1 ∂Hd (φk−1 )
+ ξkm
∂qk
tk+1 − tk
∂qk
+ 21 (qk−1 + qk )ηkm + Gm
Nk
− 12
pk−1 ξkm
pk+1 ξkm
+ 12
tk+1 − tk
tk+1 − tk
+ ηkm
− 21
[
]
− Wd (φk ) ξkm−1 − Dd (qk )τkm−1 = 0.
∂Hd (φk )
tk − tk−1 ∂Hd (φk−1 )
+ ηkm
∂pk
tk+1 − tk
∂pk
qk+1 ηkm
qk−1 ηkm
+ 12
tk+1 − tk
tk+1 − tk
m
+ τk+1
∂Hd (φk )
tk − tk−1 ∂Hd (φk−1 )
− τkm
∂tk+1
tk+1 − tk
∂tk
m
+ ξk+1
∂Hd (φk )
tk − tk−1 ∂Hd (φk−1 )
− ξkm
∂qk+1
tk+1 − tk
∂qk
m
+ ηk+1
∂Hd (φk )
tk − tk−1 ∂Hd (φk−1 )
− ηkm
∂pk+1
tk+1 − tk
∂pk
m
+ τk+1
Hd (φk )
Hd (φk−1 )
− τkm
tk+1 − tk
tk+1 − tk
pk−1 ξkm
1
+2
+ 12
tk+1 − tk
pk ξkm
tk+1 − tk
−
−
p
ξm
1 k+1 k+1
2
tk+1 − tk
qk−1 ηkm
1
2
tk+1 − tk
−
1
2
+
1
2
m
pk ξk+1
tk+1 − tk
m
qk+1 ηk+1
tk+1 − tk
m
qk ηk+1
qk ηkm
− 12
+ Dd (Gm
Nk)
tk+1 − tk
tk+1 − tk
[
]
− Wd (φk ) ξkm−1 − Dd (qk )τkm−1
[
∂Hd (φk ) tk − tk−1 ∂Hd (φk−1 )
+
∂tk
tk+1 − tk
∂tk
]
Hd (φk )
Hd (φk−1 )
−
+
tk+1 − tk
tk+1 − tk
[
tk − tk−1 ∂Hd (φk−1 )
m ∂Hd (φk )
+
+ ξk
∂qk
tk+1 − tk
∂qk
]
pk−1
pk+1
+ 12
− 12
tk+1 − tk
tk+1 − tk
=
(25)
Using (13)(15), from (25), we obtain
{ [
∂Hd (φk−1 )
Dd εm τkm (tk − tk−1 )
∂tk
∂H
(φ
)
d
k−1
+ ξkm (tk − tk−1 )
∂qk
∂Hd (φk−1 )
+ ηkm (tk − tk−1 )
+ τkm Hd (φk−1 )
∂pk
]}
− 21 (pk−1 + pk )ξkm + 12 (qk−1 + qk )ηkm + Gm
Nk
{ [
= Dd εm τkm (tk − tk−1 )D2 (R− Hd (φk ))
+ ξkm (tk − tk−1 )D4 (R− Hd (φk ))
+ 12
τkm
]
+ ηkm (tk − tk−1 )D6 (R− Hd (φk )) + τkm R− Hd (φk )
]}
m
m
m
1
1
− 2 (pk−1 + pk )ξk + 2 (qk−1 + qk )ηk + GN k
= Dd INz,d
{
qk+1 − qk m
= εm εWd (φk )
τ − εWd (φk )ξkm
tk+1 − tk k
}
[ m−1
]
m−1
+ Wd (φk ) ξk
− Dd (qk )τk
{
= εm − εWd (φk )[ξkm − Dd (qk )τkm ]
[
]}
+ Wd (φk ) ξkm−1 − Dd (qk )τkm−1
= −εz+1 Wd (φk )[ξkz − Dd (qk )τkz ].
(26)
It shows that Dd INz,d is directly proportional to εz+1 ,
so INz,d is a z -th order adiabatic invariant of the discrete dierence variational Hamilton system, which can
be called the discrete the Noether adiabatic invariant.
When z = 0, which mean that εWd (φk ) vanish, the system becomes the one without perturbation. Then the
Perturbation to Noether Symmetry and Noether Adiabatic Invariants . . .
discrete the Noether adiabatic invariant (24) turns into
the discrete the Noether exact invariant (10).
τk0 (tk , qk , pk ) = 2tk ,
Take the Emden equation to illustrate the application
of the above results. The discrete Hamiltonian of the
system and its recursive operator are
R− Hd (φk ) = Hd (φk−1 ) =
(27)
1 p2k + p2k−1
2 t2k + t2k−1
1
6
+ (t2k qk6 + t2k−1 qk−1
),
(28)
12
and the system is perturbed by the small quantity
qk+1 pk+1 − qk pk
εWd (φk ) = −ε
.
(29)
3qk tk − qk tk+1 − 2qk+1 tk
Let us study the perturbation to the the Noether symmetry and the the Noether adiabatic invariant of the discrete
dierence variational Hamilton system.
We calculate the Dj Hd (φk ) (j = 1, 2, 3, 4, 5, 6) as
tk (p2 + p2 )
D1 Hd (φk ) = − 2 k+1 2 k2 + 16 tk qk6 ,
(tk+1 + tk )
Firstly, we seek the discrete the Noether exact invariant. Suppose that the generators τk0 , ξk0 and ηk0 are linear,
i.e.
τk0 (tk , qk , pk ) = C1 tk + C2 qk + C3 pk + C4 ,
ξk0 (tk , qk , pk ) = C5 tk + C6 qk + C7 pk + C8 ,
(31)
where C1 C12 are constants.
Substituting (30) and (31) into the the Noether identity (9) of the system without perturbation, we have
[ 2
]
1 pk+1 + p2k
1 2
6
2 6
0
(t
q
+
t
q
)
+
k k Dd (τk )
2 t2k+1 + t2k
12 k+1 k+1
[
]
tk (p2 + p2 )
+ 16 tk qk6 − 2 k+1 2 k2 τk0
(tk+1 + tk )
]
[
tk+1 (p2k+1 + p2k ) 0
6
τk+1
−
+ 16 tk+1 qk+1
(t2k+1 + t2k )2
5
0
+ 12 t2k qk5 ξk0 + 12 t2k+1 qk+1
ξk+1
+
+
0
pk+1 ηk+1
t2k+1 + t2k
pk ηk0
2
tk+1 + t2k
− 12 (pk+1 + pk )Dd (ξk0 )
(32)
ξk0 (tk , qk , pk ) = −qk ,
(33)
ηk0 (tk , qk , pk ) = pk ,
the function
G0N k = −qk pk
(34)
satises the the Noether identity (9). According to Theorem 1, the discrete dierence variational Hamilton system
has the following discrete the Noether exact invariant:
IN0,d =
+
tk (p2k + p2k−1 ) 2t2k (tk − tk−1 )(p2k + p2k−1 )
−
t2k + t2k−1
(t2k + t2k−1 )2
(tk − tk−1 )p2k
6
+ 61 tk tk−1 (tk qk6 + tk−1 qk−1
)
t2k + t2k−1
(35)
+ 21 (pk−1 qk + pk qk−1 ) = const.
Secondly, we seek the rst order discrete the Noether
adiabatic invariant. We also suppose that the generators
τk1 , ξk1 and ηk1 are linear, i.e.
τk1 (tk , qk , pk ) = C1′ tk + C2′ qk + C3′ pk + C4′ ,
ξk1 (tk , qk , pk ) = C5′ tk + C6′ qk + C7′ pk + C8′ ,
tk+1 (p2k+1 + p2k ) 1
6
,
D2 Hd (φk ) = −
+ 6 tk+1 qk+1
(t2k+1 + t2k )2
5
,
D3 Hd (φk ) = 12 t2k qk5 , D4 Hd (φk ) = 12 t2k+1 qk+1
pk
pk+1
D5 Hd (φk ) = 2
, D6 Hd (φk ) = 2
.(30)
tk+1 + t2k
tk+1 + t2k
ηk0 (tk , qk , pk ) = C9 tk + C10 qk + C11 pk + C12 ,
+ 12 (qk+1 + qk )Dd (ηk0 ) + Dd (G0N k ) = 0.
It follows that when
4. An example
1 p2k+1 + p2k
1
6
Hd (φk ) =
+ (t2k+1 qk+1
+ t2k qk6 ),
2 t2k+1 + t2k
12
29
′
′
′
ηk1 (tk , qk , pk ) = C9′ tk + C10
qk + C11
pk + C12
,
(36)
′
C1′ C12
are constants.
where
Substituting (30), (31) and (36) into the determining
Eq. (23) of the perturbation to the the Noether symmetry of the system, we have
[ 2
]
1 pk+1 + p2k
1 2
6
2 6
1
+
(t
q
+
t
q
)
k k Dd (τk )
2 t2k+1 + t2k
12 k+1 k+1
[
]
tk (p2 + p2 )
+ 16 tk qk6 − 2 k+1 2 k2 τk1
(tk+1 + tk )
[
]
tk+1 (p2k+1 + p2k ) 1
6
−
+ 16 tk+1 qk+1
τk+1 + 21 t2k qk5 ξk1
(t2k+1 + t2k )2
5
1
+ 12 t2k+1 qk+1
ξk+1
+
pk ηk1
2
tk+1 + t2k
+
1
pk+1 ηk+1
2
tk+1 + t2k
− 21 (pk+1 + pk )Dd (ξk1 ) + 12 (qk+1 + qk )Dd (ηk1 )
+
[ 0
]
qk+1 pk+1 − qk pk
ξk − Dd (qk )τk0
3qk tk − qk tk+1 − 2qk+1 tk
(37)
+ Dd (G1N k ) = 0.
It follows that when
τk1 (tk , qk , pk ) = 2tk ,
ηk1 (tk , qk , pk ) = pk ,
ξk1 (tk , qk , pk ) = −qk ,
(38)
the function
G1N k = −2qk pk
(39)
satises the determining Eq. (23) of the perturbation to
M.-J. Zhang, J.-H. Fang
30
the Noether symmetry. According to Theorem 2, the
discrete dierence variational Hamilton system has the
following rst order discrete the Noether adiabatic invariant:
[
tk (p2k + p2k−1 ) 2t2k (tk − tk−1 )(p2k + p2k−1 )
IN1,d =
−
t2k + t2k−1
(t2k + t2k−1 )2
(
)
(tk − tk−1 )p2k
6
+ 16 tk tk−1 tk qk6 + tk−1 qk−1
2
2
tk + tk−1
]
+ 12 (pk−1 qk + pk qk−1 ) (1 + ε) − εqk pk .
+
(40)
Further we can obtain more high-order discrete adiabatic
invariants. In addition, it is worth noting that there are
many calculations above, which can be also done in an
automatic way using a computer algebra system stated
in Refs. [26, 27].
[3] D.F.M. Torres, in:
[4]
[5]
[6]
[7]
[8]
[9]
[10]
[11]
[12]
[13]
The perturbation to the the Noether symmetry and
the the Noether adiabatic invariants of discrete dierence
variational systems are studied in this paper. We obtain the discrete the Noether adiabatic invariant induced
directly by the perturbation to the Noether symmetry
of the system. When z = 0, this means that εWd (φk )
vanish, the system becomes the one without perturbations, and the discrete the Noether adiabatic invariant
will turn into the discrete the Noether exact invariant
(i.e. the so-called discrete the Noether conserved quantity in Ref. [14]) naturally. The results of this paper
have important signicance for further study on the discrete mechanical systems and symmetrical perturbation
theory.
[14]
Acknowledgments
[22]
This research was supported by the National Natural
Science Foundation of China (grant No. 11103086) and
Fundamental Research Funds for the Central Universities
of China (grant No. 09CX04018A).
[23]
[1] F.X. Mei,
Symmetries and Conserved Quantities of
, Beijing Institute of
Technology Press, Beijing 2004 (in Chinese).
[2] S.K. Luo, Y.F. Zhang, et al., Advances in the Study of
Dynamics of Constrained Mechanics Systems, Science
Press, Beijing 2008 (in Chinese).
Constrained Mechanical Systems
, Eds. F. Colonius, L. Grüne, Vol. 273,
Springer, Berlin 2002, p. 287.
D.F.M. Torres, Eur. J. Control. 8, 56 (2002).
D.F.M. Torres, Commun. Pure Appl. Anal. 3, 491
(2004).
D. Levi, R. Yamilov, J. Math. Phys. 38, 6648 (1997).
D. Levi, S. Tremblay, P. Winternitz, J. Phys. A,
Math. Gen. 33, 8507 (2000).
D. Levi, S. Tremblay, P. Winternitz, J. Phys. A,
Math. Gen. 34, 9507 (2001).
V. Dorodnitsyn, Appl. Numer. Math. 39, 307 (2001).
Z. Bartosiewicz, D.F.M. Torres, J. Math. Anal. Appl.
342, 1220 (2008).
N. Martins, D.F.M. Torres, Appl. Math. Lett. 23,
1432 (2010).
Z. Bartosiewicz, N. Martins, D.F.M. Torres, Eur. J.
Control. 17, 9 (2011).
S.Y. Shi, J.L. Fu, L.Q. Chen, Chin. Phys. B 17, 385
(2008).
S.Y. Shi, X.H. Huang, X.B. Zhang, L. Jin, Acta Phys.
Sin. 58, 3625 (2009) (in Chinese).
J.L. Fu, H. Fu, R.W. Liu, Phys. Lett. A 374, 1812
(2010).
J.M. Burgers, Ann. Phys. 357, 195 (1917).
Y.Y. Zhao, F.X. Mei, Symmetries and Invariants of
Mechanical Systems, Science Press, Beijing 1999 (in
Chinese).
J.L. Fu, L.Q. Chen, Phys. Lett. A 324, 95 (2004).
X.W. Chen, Y.M. Li, Y.H. Zhao, Phys. Lett. A 337,
274 (2005).
M.J. Zhang, J.H. Fang, X.N. Zhang, K. Lu, Chin.
Phys. B 17, 1957 (2008).
N. Ding, X.F. Chen, J.H. Fang, C.Z. Liu, Phys.
Lett. A 373, 3005 (2009).
M.J. Zhang, J.H. Fang, K. Lu, T. Pang, P. Lin, Commun. Theor. Phys. 51, 961 (2009).
M.J. Zhang, J.H. Fang, K. Lu, Int. J. Theor. Phys.
49, 427 (2010).
M.J. Zhang, J.H. Fang, K. Lu, K.J. Zhang, Y. Li,
Chin. Phys. Lett. 26, 120201 (2009).
P. Wang, H.J. Zhu, Acta Phys. Pol. A 119, 298
(2011).
P.D.F. Gouveia, D.F.M. Torres, E.A.M. Rocha, Control Cybernet. 35, 831 (2006).
P.D.F. Gouveia, D.F.M. Torres, Nonlinear Anal. 71,
e138 (2009).
Inform. Sci.
5. Conclusion
References
Dynamics, Bifurcations, and Con-
trol, Kloster Irsee, 2001, Lecture Notes in Control and
[15]
[16]
[17]
[18]
[19]
[20]
[21]
[24]
[25]
[26]
[27]