北京大学学报自然科学版   2016, Vol. 52 Issue(4): 681-686

文章信息

楼智美, 王元斌, 谢志堃
LOU Zhimei, WANG Yuanbin, XIE Zhikun
典型微扰力学系统的近似Lie对称性、近似Noether对称性和近似Mei对称性
Approximate Lie Symmetries, Approximate Noether Symmetries and Approximate Mei Symmetries of Typical Perturbed Mechanical System
北京大学学报(自然科学版), 2016, 52(4): 681-686
Acta Scientiarum Naturalium Universitatis Pekinensis, 2016, 52(4): 681-686

文章历史

收稿日期: 2015-10-09
修回日期: 2016-02-13
网络出版日期: 2016-07-14
典型微扰力学系统的近似Lie对称性、近似Noether对称性和近似Mei对称性
楼智美1, 王元斌2, 谢志堃1     
1. 绍兴文理学院物理系, 绍兴 312000;
2. 绍兴文理学院数学系, 绍兴 312000
摘要: 利用3种近似对称性方法(近似Lie对称性法、近似Noether对称性法和近似Mei对称性法)研究典型微扰力学系统的一阶近似对称性和近似守恒量。结果表明, 利用近似Lie对称性法找到的6个一阶近似对称性和近似守恒量与利用近似Noether对称性法找到的相同, 而利用近似Mei对称性法只能找到其中5个一阶近似对称性和近似守恒量。
关键词: 微扰力学系统     近似Lie对称性     近似Noether对称性     近似Mei对称性     近似守恒量    
Approximate Lie Symmetries, Approximate Noether Symmetries and Approximate Mei Symmetries of Typical Perturbed Mechanical System
LOU Zhimei1, WANG Yuanbin2, XIE Zhikun1     
1. Department of Physics, Shaoxing University, Shaoxing 312000;
2. Department of mathematics, Shaoxing University, Shaoxing 312000
Corresponding author: LOU Zhimei, louzhimei@usx.edu.cn
Abstract: Three methods, which are approximate Lie symmetry method, approximate Noether symmetry method and approximate Mei symmetry method, are adopted to study the first order approximate symmetries and approximate conserved quantities of a typical perturbed mechanical system. Six identical first order approximate symmetries and approximate conserved quantities of the typical perturbed mechanical system are obtained by approximate Lie symmetry method and approximate Noether symmetry method, but only five of them can be obtained by approximate Mei symmetry method.
Key words: perturbed mechanical system     approximate Lie symmetry     approximate Noether symmetry     approximate Mei symmetry     approximate conserved quantity    

分析力学中研究力学系统的对称性与守恒量有3种对称性方法[1-2]: Lie对称性法、Noether对称性法和Mei对称性法。引进群无限小变换, 微分方程在此变换下保持不变为Lie对称性, 哈密顿作用量在此变换下保持不变为Noether对称性, 力学系统的动力学函数在此变换下仍然满足运动方程为Mei对称性。事实上, 许多实际力学系统的某些参数常常会随着位移、速度和时间的变化发生微小的变化, 即力学系统受到微扰作用, 这样的力学系统称为微扰力学系统。此类系统的近似对称性和近似守恒量研究对于研究力学系统的特性至关重要。目前, 研究微扰力学系统近似对称性与近似守恒量有两种近似对称性法:近似Lie对称性法[3]和近似Noether对称性法[4]。引进近似的群无限小变换, 微分方程在此变换下近似保持不变为近似Lie对称性; 哈密顿作用量在此变换下近似保持不变则为近似Noether对称性, 所得的守恒量为近似守恒量。

近年来, 关于常微分方程、偏微分方程近似对称性和近似守恒量的研究已取得很多成果[3-16], 文献[3-14]采用的方法都是近似Lie对称性法或近似Noether对称性法, 其中文献[3-11]侧重近似对称性理论的研究, 文献[12-14]侧重近似对称性理论的实际应用, 文献[15]提出用直接积分法求近似守恒量的方法, 文献[16]将直接积分法应用于求二阶近似守恒量。但是, 到目前为止, 还没有建立用近似Mei对称性研究近似守恒量的理论, 关于3个近似对称性理论相互关系方面也没有深入的研究。事实上, 近似Lie对称性法和近似Noether对称性法是在Lie对称性法和Noether对称性法的基础上发展起来的, 那么在Mei对称性法的基础上也可以建立相应的近似Mei对称性法, 即引进近似的群无限小变换, 力学系统的动力学函数在此变换下近似满足运动方程为近似Mei对称性。

本文分析近似Lie对称性、近似Noether对称性和Mei对称性理论, 讨论3种对称性间的关系, 并以频率比为2:1的弱非线性耦合谐振子为例, 研究系统的一阶近似对称性和近似守恒量。

1 近似Lie对称性、近似Noether对称性和近似Mei对称性理论 1.1 近似Lie对称性理论[3]

具有n个自由度的微扰Lagrange力学系统, 其Lagrange函数可以表示为$L=L (t, {q_s}, {\dot q_s}, \varepsilon)$, 其中qs为广义坐标, ${\dot q_s}$为广义速度, $0 < \varepsilon \ll 1$为微扰系数。微扰Lagrange力学系统的运动微分方程可以表示为

$ \frac{\text{d}}{\text{d}t}\frac{{\partial L}}{{\partial {{\dot q}_s}}}-\frac{{\partial L}}{{\partial {q_s}}} = 0\;\;{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} (s = 1, \;2, \;..., \;n), $ (1)

方程(1)可简写为

$ {E_s}(L) = 0, $ (2)

其中${E_s}=\frac{\text{d}}{\text{d}t}\frac{\partial }{{\partial {{\dot q}_s}}}-\frac{\partial }{{\partial {q_s}}}$为Euler算子。假设系统(式(1))非奇异, 即设

$ \det \left( {\frac{{{\partial ^2}L}}{{\partial {{\dot q}_s}\partial {{\dot q}_k}}}} \right) \ne 0, $ (3)

则可求出所有的广义加速度:

$ {\ddot q_s} = {a_s}(t, {q_s}, {\dot q_s}, \varepsilon )。 $ (4)

引进近似的群无限小变换

$ \left\{ \begin{array}{l} {t^ * } = t + \delta \tau (t, {q_s}, {{\dot q}_s}, \varepsilon ), \\ q_s^ * ({t^ * }) = {q_s}(t) + \delta {\xi _s}(t, {q_s}, {{\dot q}_s}, \varepsilon ), \end{array} \right. $ (5)

其中s=1, 2, …, n, δ为无限小参数, τξs为无限小变换生成元。式(5)的无限小生成元向量为

$ {X^{(0)}} = \tau \frac{\partial }{{\partial t}} + \sum\limits_{s = 1}^n {{\xi _s}\frac{\partial }{{\partial {q_s}}}}, $ (6)

式(6)的一次扩展为

$ {X^{(1)}} = {X^{(0)}} + \sum\limits_{s = 1}^n {({{\dot \xi }_s}-{{\dot q}_s}\dot \tau )\frac{\partial }{{\partial {{\dot q}_s}}}}, $ (7)

二次扩展为

$ {X^{(2)}} = {X^{(1)}} + \sum\limits_{s = 1}^n {({{\ddot \xi }_s}-{{\dot q}_s}\ddot \tau-2{{\ddot q}_s}\dot \tau )\frac{\partial }{{\partial {{\ddot q}_s}}}} 。 $ (8)

式(5)~(8)中,

$ \tau = \sum\limits_{i = 0}^k {{\varepsilon ^i}} {\tau _i}, {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \;{\xi _s} = \sum\limits_{i = 0}^k {{\varepsilon ^i}{\xi _{si}}}, $ (9)

其中s=1, 2, …, n, k为微扰项的阶数。

运动微分方程(1)的k阶近似Lie对称性是指式(4)在近似的群无限小变换(式(5))下近似保持不变[3], 即

$ {X^{(2)}}({\ddot q_s}-{a_s}) = O({\varepsilon ^{k + 1}})。 $ (10)

若存在规范函数

$ \begin{array}{l} G = G({q_s}, {{\dot q}_s}, \varepsilon ) = \sum\limits_{i = 0}^k {{\varepsilon ^i}{G_i}} \\ {\kern 1pt} (s = 1, \;2, \;...\;, \;n) \end{array} $ (11)

满足

$ \begin{array}{l} \;\;\;\frac{{\partial L}}{{\partial t}}\tau + \sum\limits_{s = 1}^n {\frac{{\partial L}}{{\partial {q_s}}}{\xi _s} + \sum\limits_{s = 1}^n {\frac{{\partial L}}{{\partial {{\dot q}_s}}}{{\dot \xi }_s}} + } \left( {L-\sum\limits_{s = 1}^n {\frac{{\partial L}}{{\partial {{\dot q}_s}}}} {{\dot q}_s}} \right)\dot \tau \\ =-\dot G, \end{array} $ (12)

则系统存在k阶近似守恒量

$ I = \sum\limits_{i = 0}^k {{\varepsilon ^i}} {I_i}, $ (13)

$ I = L\tau + \sum\limits_{s = 1}^n {\frac{{\partial L}}{{\partial {{\dot q}_s}}}({\xi _s}-{{\dot q}_s}\tau ) + G}, $ (14)

满足

$ \frac{{{\rm{d}}I}}{{{\rm{d}}t}} = O({\varepsilon ^{k + 1}})。 $ (15)

将式(4)代入式(10)并展开, 令${\varepsilon ^0}, {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\varepsilon ^1}, \; ..., \; {\varepsilon ^k}$的系数为0, 可求得无限小生成元τξs。将所得生成元代入式(12), 并比较等式两边${\varepsilon ^0}, {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\varepsilon ^1}, \; ..., \; {\varepsilon ^k}$的系数可求得规范函数G。将Lagrange函数L及求得的无限小生成元τξs以及规范函数G代入式(14), 可求得系统的近似守恒量I

在式(13)中, 若${I_0} \ne 0, \, \, \, \, {I_i}\; (i=1, \; 2, \; ..., \; k)$均为0, 则称I为未受微扰力学系统的精确守恒量, 相应的对称性为未受微扰力学系统的精确对称性。若${I_0} \ne 0, \, \, $且至少有一个${I_i}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} (i=1, \; 2, \; ..., \; k)$不为0, 则称I为微扰力学系统稳定的近似守恒量, 相应的对称性为稳定的近似对称性。若${I_0}=0, $且至少有一个${I_i}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} (i=1, \; 2, \; ..., \; k)$不为0, 则称I为微扰力学系统平凡的近似守恒量, 相应的对称性为平凡的近似对称性。

1.2 近似Noether对称性理论

近似Noether对称性[4]指在近似的群无限小变换(式(5))下, 哈密顿作用量在此变换下近似保持不变, 即

$ \begin{array}{l} \;\;\;\int_{{t_1}}^{{t_2}} {L(t, \;{q_s}, } \;{{\dot q}_s}, \;\varepsilon ){\rm{d}}t\\ = \int_{t_1^*}^{t_2^*} {L(t, \;q_s^*, } \;\dot q_s^*, \;\varepsilon ){\rm{d}}{t^*} + O({\varepsilon ^{k + 1}})\;, \end{array} $ (16)

则称无限小变换(式(5))为近似Noether对称变换。

若无限小生成元(式(9))和规范函数(式(11))满足确定方程(12), 则存在相应的近似守恒量(式(14))。

将式(9)和(11)代入式(12), 并比较式(12)两边${\varepsilon ^0}, {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\varepsilon ^1}, \; ..., \; {\varepsilon ^k}$的系数, 可求得无限小生成元τξs以及规范函数G。将Lagrange函数L及求得的无限小生成元τξs以及规范函数G代入式(14), 可求得系统的近似守恒量I

1.3 近似Mei对称性理论

假设经无限小变换(5)后, 系统的Lagrange函数$L (t, \; q, \; \dot q, \; \varepsilon)$变成L*:

$ \begin{array}{l} {L^*} = L\left( {{t^*}, {q^*}, \frac{{{\rm{d}}{q^*}}}{{{\rm{d}}{t^*}}}, \;\varepsilon } \right)\\ \;\;\;\; = L(t, q, \dot q, \varepsilon ) + \delta {X^{(1)}}(L) + O({\delta ^2})。 \end{array} $ (17)

如果用变换后的Lagrange函数L*代替变换前的L时, 方程(2)的形式近似保持不变, 即

$ {E_s}({L^*}) = O({\varepsilon ^{k + 1}}), $ (18)

那么称这种不变性为系统的k阶近似Mei对称性。

将式(17)代入式(18), 忽略δ2及更高阶小量, 并利用式(2), 得到

$ {E_s}({X^{(1)}}(L)) = O({\varepsilon ^{k + 1}}) $ (19)

对于微扰Lagrange力学系统(式(2)), 如果无限小生成元$\tau, {\xi _s}$满足方程(19), 则相应的近似不变性为系统的k阶近似Mei对称性。方程(19)称为近似Mei对称性的判据方程。

若无限小生成元(式(9))和规范函数(11)满足确定方程(12), 则存在相应的近似守恒量(式(14))。

将Lagrange函数L代入式(19)并展开, 令ε0, ε1, …, εk的系数为0, 可求得无限小生成元τξs。将所得生成元代入式(12), 并比较等式两边ε0, ε1, …, εk的系数可求得规范函数G。将Lagrange函数L及求得的无限小生成元τξs以及规范函数G代入式(14), 可求得系统的近似守恒量I

1.4 3种对称性的关系

若无限小生成元(式(9))满足式(10), 说明系统具有近似Lie对称性。若无限小生成元(式(9))满足式(10)且满足式(12), 并能找到相应的规范函数, 则说明系统既具有近似Lie对称性, 又具有近似Noether对称性, 找到的近似守恒量既是近似Lie对称性守恒量, 又是近似Noether对称性守恒量。

若无限小生成元(式(9))满足式(19), 则说明系统具有近似Mei对称性。若无限小生成元(式(9))同时满足式(10)和(19), 则说明系统同时具有近似Lie对称性和近似Mei对称性。

若无限小生成元(式(9))同时满足式(10), (12)和(19), 并能找到相应的规范函数, 则说明系统同时具有近似Lie对称性、近似Mei对称性和近似Noether对称性, 找到的近似守恒量既是近似Lie对称性守恒量, 又是近似Mei对称性守恒量, 也是近似Noether对称性守恒量。

2 典型微扰力学系统的一阶近似对称性和一阶近似守恒量

频率比为2:1的弱非线性耦合谐振子的Lagrange函数[13]

$ L = \frac{1}{2}\dot x_1^2 + \frac{1}{2}\dot x_2^2-2x_1^2-\frac{1}{2}x_2^2 + \varepsilon {x_1}x_2^2, $ (20)

运动微分方程为

$ \left\{ \begin{array}{l} {{\ddot x}_1} =-4{x_1} + \varepsilon x_2^2, \\ {{\ddot x}_2} =-{x_2} + 2\varepsilon {x_1}{x_2}。 \end{array} \right. $ (21)
2.1 一阶近似Lie对称性与近似守恒量

研究系统的一阶近似Lie对称性与近似守恒量。将式(21)代入式(10)并展开, 令${\varepsilon ^0}, {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\varepsilon ^1}$的系数为0, 可求得如下6组生成元[13]:

$ \left\{ \begin{array}{l} {\tau _0} =-1, \;\;{\tau _1} = {\xi _{11}} = {\xi _{21}} = 0, \\ {\xi _{10}} = {\xi _{20}} = 0, \end{array} \right. $ (22a)
$ \left\{ \begin{array}{l} {\tau _0} = {\xi _{10}} = {\xi _{20}} = 0, \;{\tau _1} =-1, \;\\ {\xi _{11}} = {\xi _{21}} = 0, \end{array} \right. $ (22b)
$ \left\{ \begin{array}{l} {\tau _0} = {\xi _{10}} = {\xi _{20}} = 0, \;{\tau _1} =-1, \;\\ {\xi _{11}} = {\xi _{21}} = 0, \end{array} \right. $ (22c)
$ \left\{ \begin{array}{l} {\tau _0} = {\xi _{10}} = {\xi _{20}} = 0, \;{\tau _1} =-1, \;\\ {\xi _{11}} =-{{\dot x}_1}, \;{\xi _{21}} = 0, \end{array} \right. $ (22d)
$ \left\{ \begin{array}{l} {\tau _0} = {\xi _{10}} = {\xi _{20}} = 0, \;{\tau _1} = 0, \\ {\xi _{11}} = {x_2}{{\dot x}_2}, \;{\xi _{21}} = {x_2}{{\dot x}_1}-2{x_1}{{\dot x}_2}, \end{array} \right. $ (22e)
$ \left\{ \begin{array}{l} {\tau _0} = 0, \;{\tau _1} = 0, \\ {\xi _{10}} = {x_2}{{\dot x}_2}, \;{\xi _{20}} = {x_2}{{\dot x}_1}-2{x_1}{{\dot x}_2}, \;\\ {\xi _{11}} = \frac{{-8{x_1}{x_2}{{\dot x}_2} + 5x_2^2{{\dot x}_1}-3{{\dot x}_1}\dot x_2^2}}{8}, \\ {\xi _{21}} = \frac{1}{8}( - 4x_1^2{{\dot x}_2} - 8{x_1}{x_2}{{\dot x}_1} + \\ \;\;\;\;\;\;3x_2^2{{\dot x}_2} - 3\dot x_1^2{{\dot x}_2} + 3\dot x_2^3)。 \end{array} \right. $ (22f)

说明频率比为2:1的弱非线性耦合谐振子系统具有6个一阶近似Lie对称性。将式(20)和(22)代入式(12), 并比较等式两边${\varepsilon ^0}, {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\varepsilon ^1}$的系数, 可求得与上述6组生成元相应的规范函数[13]:

$ {G_0} = 0, \;\;{G_1} = 0, $ (23a)
$ {G_0} = 0, \;\;{G_1} = 0, $ (23b)
$ {G_0} = 0, \;\;\;{G_1} = \frac{1}{2}\dot x_2^2-\frac{1}{2}x_2^2, $ (23c)
$ {G_0} = 0, \;\;\;{G_1} = \frac{1}{2}\dot x_1^2-2x_1^2, $ (23d)
$ {G_0} = 0, \;{G_1} = {x_1}x_2^2 + {x_1}\dot x_2^2-{x_2}{\dot x_1}{\dot x_2}, $ (23e)
$ \left\{ \begin{array}{l} {G_0} = {x_1}x_2^2 + {x_1}\dot x_2^2-{x_2}{{\dot x}_1}{{\dot x}_2}, \\ {G_1} = \frac{1}{{32}}(8x_1^2\dot x_2^2-8x_1^2x_2^2 + \\ \;\;\;\;\;\;32{x_1}{x_2}{{\dot x}_1}{{\dot x}_2}-5x_2^4 - 10x_2^2\dot x_1^2 - \\ \;\;\;\;\;\;6x_2^2\dot x_2^2 + 18\dot x_1^2\dot x_2^2 - 9\dot x_2^4)。 \end{array} \right. $ (23f)

将式(20), (22)和(23)代入式(14), 得到6个一阶近似守恒量[13]:

$ {I^1} = \frac{1}{2}\dot x_1^2 + \frac{1}{2}\dot x_2^2 + 2x_1^2 + x_2^2-\varepsilon {x_1}x_2^2, $ (24a)
$ {I^2} = \varepsilon \left( {\frac{1}{2}\dot x_1^2 + \frac{1}{2}\dot x_2^2 + 2x_1^2 + x_2^2} \right) = \varepsilon I_0^1, $ (24b)
$ {I^3} = \varepsilon \left( {\frac{1}{2}\dot x_1^2 + 2x_1^2} \right) = \varepsilon I_0^2, $ (24c)
$ {I^4} = \varepsilon \left( {\frac{1}{2}\dot x_2^2 + \frac{1}{2}x_2^2} \right) = \varepsilon I_0^3, $ (24d)
$ {I^5} = \varepsilon ({x_1}x_2^2-{x_1}\dot x_2^2 + {x_2}{\dot x_1}{\dot x_2}) = \varepsilon I_0^4, $ (24e)
$ \begin{array}{l} {I^6} = {x_1}x_2^2-{x_1}\dot x_2^2 + {x_2}{{\dot x}_1}{{\dot x}_2} + \\ \;\;\;\;\;\;\frac{\varepsilon }{{32}}(10x_2^2\dot x_1^2 + 6x_2^2\dot x_2^2 + 3\dot x_2^4-\\ \;\;\;\;\;\;8x_1^2x_2^2-8x_1^2\dot x_2^2 - 32{x_1}{x_2}{{\dot x}_1}{{\dot x}_2} - \\ \;\;\;\;\;\;5x_2^4 - 6\dot x_1^2\dot x_2^2), \end{array} $ (24f)

其中,

$ I_0^1 = \frac{1}{2}\dot x_1^2 + \frac{1}{2}\dot x_2^2 + 2x_1^2 + x_2^2, $ (25a)
$ I_0^2 = \frac{1}{2}\dot x_1^2 + 2x_1^2, $ (25b)
$ I_0^3 = \frac{1}{2}\dot x_2^2 + \frac{1}{2}x_2^2, $ (25c)
$ I_0^4 = {x_1}x_2^2-{x_1}\dot x_2^2 + {x_2}{\dot x_1}{\dot x_2}。 $ (25d)

I01, I02, I03I04为不受非线性耦合项作用的二维各向异性谐振子的4个守恒量, 均满足$\frac{{{\rm{d}}I_0^i}}{{{\rm{d}}t}}=0{\kern 1pt} $(i=1, 2, 3, 4), 其中I01为系统的总能量, I02I03为两个分振子的能量, 这4个守恒量不是相互独立的, 存在$I_0^1=I_0^2 + I_0^3$的关系。I1为弱非线性耦合谐振子的哈密顿函数, I2, I3, I4I5为平凡的一阶近似守恒量, I6为稳定的一阶近似守恒量。

2.2 一阶近似Noether对称性与近似守恒量

研究系统的一阶近似Noether对称性与一阶近似守恒量。式(22)表示的6组生成元和式(23)表示的6个规范函数均符合Noether恒等式(12), 因此, 系统同时具有式(22)表示的6个一阶近似Noether对称性和式(24)表示的6个一阶近似Noether守恒量。

2.3 一阶近似Mei对称性与近似守恒量

研究系统的一阶近似Mei对称性与一阶近似守恒量。将式(20)代入式(19), 只能求得如下5组生成元:

$ \left\{ \begin{array}{l} {\tau _0} =-1, \;\;{\tau _1} = {\xi _{11}} = {\xi _{21}} = 0, \\ {\xi _{10}} = {\xi _{20}} = 0, \end{array} \right. $ (26a)
$ \left\{ \begin{array}{l} {\tau _0} = {\xi _{10}} = {\xi _{20}} = 0, \;{\tau _1} =-1, \;\\ {\xi _{11}} = {\xi _{21}} = 0, \end{array} \right. $ (26b)
$ \left\{ \begin{array}{l} {\tau _0} = {\xi _{10}} = {\xi _{20}} = 0, \;{\tau _1} =-1, \;\\ {\xi _{11}} = 0, \;{\xi _{21}} =-{{\dot x}_2}, \end{array} \right. $ (26c)
$ \left\{ \begin{array}{l} {\tau _0} = {\xi _{10}} = {\xi _{20}} = 0, \;{\tau _1} =-1, \;\\ {\xi _{11}} =-{{\dot x}_1}, \;{\xi _{21}} = 0, \end{array} \right. $ (26d)
$ \left\{ \begin{array}{l} {\tau _0} = {\xi _{10}} = {\xi _{20}} = 0, \;{\tau _1} = 0, \\ {\xi _{11}} = {x_2}{{\dot x}_2}, \;{\xi _{21}} = {x_2}{{\dot x}_1}-2{x_1}{{\dot x}_2}。 \end{array} \right. $ (26e)

这5组生成元与式(22)中的前5组相同。将式(20)和(26)代入式(12), 并比较等式两边${\varepsilon ^0}, {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\varepsilon ^1}$的系数, 可求得与上述5组生成元相应的规范函数:

$ {G_0} = 0, \;\;{G_1} = 0, $ (27a)
$ {G_0} = 0, \;\;{G_1} = 0, $ (27b)
$ {G_0} = 0, \;\;{G_1} = \frac{1}{2}\dot x_2^2-\frac{1}{2}x_2^2, $ (27c)
$ {G_0} = 0, \;\;{G_1} = \frac{1}{2}\dot x_1^2-2x_1^2, $ (27d)
$ {G_0} = 0, \;\;{G_1} = {x_1}x_2^2 + {x_1}\dot x_2^2-{x_2}{\dot x_1}{\dot x_2}。 $ (27e)

这5个规范函数与式(23)中前5个规范函数相同。

将式(20), (26)和(27)代入式(14), 得到5个一阶近似守恒量:

$ {I^1} = \frac{1}{2}\dot x_1^2 + \frac{1}{2}\dot x_2^2 + 2x_1^2 + x_2^2-\varepsilon {x_1}x_2^2, $ (28a)
$ {I^2} = \varepsilon \left( {\frac{1}{2}\dot x_1^2 + \frac{1}{2}\dot x_2^2 + 2x_1^2 + x_2^2} \right) = \varepsilon I_0^1, $ (28b)
$ {I^3} = \varepsilon \left( {\frac{1}{2}\dot x_1^2 + 2x_1^2} \right) = \varepsilon I_0^2, $ (28c)
$ {I^4} = \varepsilon \left( {\frac{1}{2}\dot x_2^2 + \frac{1}{2}x_2^2} \right) = \varepsilon I_0^3, $ (28d)
$ {I^5} = \varepsilon ({x_1}x_2^2-{x_1}\dot x_2^2 + {x_2}{\dot x_1}{\dot x_2}) = \varepsilon I_0^4。 $ (28e)

这5个守恒量与(24)式中的前5个守恒量相同。由近似Mei对称性只能求得微扰力学系统的哈密顿函数和4个平凡的一阶近似守恒量, 不能求得稳定的一阶近似守恒量。

3 结论

本文阐述了3种近似对称性理论及其相互关系, 并用3种近似对称性理论研究了频率比为2:1的弱非线性耦合谐振子的近似对称性与近似守恒量, 结果表明, 利用近似Lie对称性法和近似Noether对称性法能找到6个相同的一阶近似对称性和近似守恒量。6个近似守恒量中, 1个是系统的哈密顿函数; 4个是平凡的一阶近似守恒量, 1个是稳定的一阶近似守恒量。用近似Mei对称性法只能找到5个一阶近似对称性和近似守恒量, 且5个近似守恒量与用近似Lie对称性法和近似Noether对称性法找到的其中5个相同。用近似Mei对称性法只找到哈密顿函数和4个平凡的一阶近似守恒量, 不能找到稳定的一阶近似守恒量。结果说明, 系统要具有近似Mei对称性的条件更严。

参考文献
[1] 梅凤翔. 李群和李代数对约束力学系统的应用. 北京: 科学出版社, 1999 .
[2] 梅凤翔. 约束力学系统的对称性与守恒量. 北京: 北京理工大学出版社, 2004 .
[3] Leach P G L, Moyo S, Cotsakis S, et al. Symmetry, singularities and integrability in complex dynamics Ⅲ: approximate symmetries and invariants. Journal of Nonlinear Mathematical Physics , 2001, 8 (1) : 139–156 DOI:10.2991/jnmp.2001.8.1.11 .
[4] Govinder K S, Heil T G, Uzer T. Approximate Noether symmetries. Physics Letters A , 1998, 240 (3) : 127–131 DOI:10.1016/S0375-9601(98)00067-X .
[5] Naeem I, Mahomed F M. Approximate first integrals for a system of two coupled van der Pol oscillators with linear diffusive coupling. Mathematical and Computational Applications , 2010, 15 (4) : 720–731 .
[6] Unal G. Approximate generalized symmetries, normal forms and approximate first integrals. Physics Letters A , 2000, 266 (2) : 106–122 .
[7] Dolapci I T, Pakdemirli M. Approximate symmetries of creeping flow equations of a second grade fluid. International Journal of Non-linear Mechanics , 2004, 39 (10) : 1603–1619 DOI:10.1016/j.ijnonlinmec.2004.01.002 .
[8] Kara A H, Mahomed F M, Qadir A. Approximate symmetries and conservation laws of the geodesic equations for the Schwarzschild metric. Nonlinear Dynamics , 2008, 51 (1/2) : 183–188 .
[9] Grebenev V N, Oberlack M. Approximate Lie symmetries of the Navier-Stokes equations. Journal of Non-linear Mathematical Physics , 2007, 14 (2) : 157–163 DOI:10.2991/jnmp.2007.14.2.1 .
[10] Johnpillai A G, Kara A H, Mahomed F M. Approximate Noether-type symmetries and conser-vation laws via partial Lagrangians for PDEs with a small parameter. Journal of Computational and Applied Mathematics , 2009, 223 (1) : 508–518 DOI:10.1016/j.cam.2008.01.020 .
[11] Zhang Z Y, Yong X L, Chen Y F. A new method to obtain approximate symmetry of nonlinear evolution equation form perturbations. Chinese Physics B , 2009, 18 (7) : 2629–2633 DOI:10.1088/1674-1056/18/7/002 .
[12] 楼智美. 两自由度弱非线性耦合系统的一阶近似Lie对称性与近似守恒量. 物理学报 , 2013, 62 (22) : 220202.
[13] 楼智美, 梅凤翔, 陈子栋. 弱非线性耦合二维各向异性谐振子的一阶近似Lie对称性与近似守恒量. 物理学报 , 2012, 61 (11) : 110204.
[14] 楼智美. 微扰Kepler系统轨道微分方程的近似Lie对称性与近似不变量. 物理学报 , 2010, 59 (10) : 6764–6769.
[15] 楼智美. 含非线性微扰项的二阶动力学系统的一阶近似守恒量的一种新求法. 物理学报 , 2014, 63 (6) : 060202.
[16] 楼智美. 两自由度微扰力学系统的二阶近似守恒量. 动力学与控制学报 , 2015, 13 (3) : 165–169.