前言
这篇文章是交通大学物理学院开设的理论力学的课程PPT缩略版。哈密顿力学是泛函的一大体现之一,而泛函作为如今十大数学技能应该被学习。为了方便大家的阅读,特制了本篇网页版的董兵老师的PPT。
提示:本篇文章为私人笔记,若有侵权,请联系我删除
哈密顿原理
经典变分问题
最速落径问题
在垂直平面内,连接不在同一铅直线上的两点A、B间的曲线,使质点在自重作用下无初速以最短时间从A滑至B。
最短线程问题
曲面上给定两点间的长度最短的曲线。
等周问题
长度固定的平面封闭曲线所围面积最大的曲线形状。
数学基础
区别函数和泛函
函数:f ( x ) f(x) f ( x ) ,自变量x x x ,定义域为数。
泛函:F [ f , x ] = F [ f ( x ) , x ] F[f, x]=F[f(x), x] F [ f , x ] = F [ f ( x ) , x ] ,自变量f ( x ) , f(x), f ( x ) , 定义域为函数。
变分—自变量不变条件下函数自身的变化
微分 x → x + d x f ( x ) → f ( x + d x ) d f = f ( x + d x ) − f ( x ) = f ′ ( x ) d x \begin{aligned} &x \rightarrow x+d x \\ & f(x) \rightarrow f(x+d x) \\ d f=& f(x+d x)-f(x)=f^{\prime}(x) d x \end{aligned} d f = x → x + d x f ( x ) → f ( x + d x ) f ( x + d x ) − f ( x ) = f ′ ( x ) d x
变分 f ( x ) → f ˉ ( x ) f(x) \rightarrow \bar{f}(x) f ( x ) → f ˉ ( x )
δ f = f ˉ ( x ) − f ( x ) \delta f=\bar{f}(x)-f(x) δ f = f ˉ ( x ) − f ( x )
变分运算规则
{ δ ( y 1 ± y 2 ) = δ y 1 + δ y 2 δ ( y 1 ⋅ y 2 ) = y 1 δ y 2 + y 2 δ y 1 δ ( y 1 y 2 ) = y 2 δ y 1 − y 1 δ y 2 y 2 2 δ ( k y ) = k δ y δ ( d y ) = d ( δ y ) δ ( d y d x ) = d d x ( δ y ) ( δ x = 0 ) δ ∫ t 1 t 2 y d t = ∫ t 1 t 2 δ y d t \left\{\begin{array}{l}\delta\left(y_{1} \pm y_{2}\right)=\delta y_{1}+\delta y_{2} \\ \delta\left(y_{1} \cdot y_{2}\right)=y_{1} \delta y_{2}+y_{2} \delta y_{1} \\ \delta\left(\frac{y_{1}}{y_{2}}\right)=\frac{y_{2} \delta y_{1}-y_{1} \delta y_{2}}{y_{2}^{2}} \\ \delta(k y)=k \delta y \\ \delta(d y)=d(\delta y) \\ \delta\left(\frac{d y}{d x}\right)=\frac{d}{d x}(\delta y) \quad(\delta x=0) \\ \delta \int_{t_{1}}^{t_{2}} y d t=\int_{t_{1}}^{t_{2}} \delta y d t\end{array}\right. ⎩ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎧ δ ( y 1 ± y 2 ) = δ y 1 + δ y 2 δ ( y 1 ⋅ y 2 ) = y 1 δ y 2 + y 2 δ y 1 δ ( y 2 y 1 ) = y 2 2 y 2 δ y 1 − y 1 δ y 2 δ ( k y ) = k δ y δ ( d y ) = d ( δ y ) δ ( d x d y ) = d x d ( δ y ) ( δ x = 0 ) δ ∫ t 1 t 2 y d t = ∫ t 1 t 2 δ y d t
积分形式的泛函及其极值问题
例如泛函F [ y ] = ∫ x 1 x 2 L ( x , y , y ′ ) d x F[y]=\int_{x_{1}}^{x_{2}} L\left(x, y, y^{\prime}\right) d x F [ y ] = ∫ x 1 x 2 L ( x , y , y ′ ) d x 为两端点固定的积分型泛函,就像在我们研究y = f ( x ) y=f(x) y = f ( x ) 时我们希望其导数为0一样,我们也希望泛函满足条件δ F = 0 \delta F=0 δ F = 0 。
最速降线问题
为了更好地明白这个问题,我们使用最速降线来举例子
在铅直平面内在所有连接两个定点A和B的曲线中,找出一 条曲线来,使得初速度为零的质点,在重力作用下,自全点沿 它无摩擦地滑下时,以最短时间到达B点。
设曲线AB方程为y = f ( x ) y=f(x) y = f ( x ) ,我们可以得知质点沿着曲线运动的距离为
v = 2 g y = d s d t = ( d x ) 2 + ( d y ) 2 d t = ( 1 ± y 2 ) d t d x v=\sqrt{2 g y}=\frac{d s}{d t}=\frac{\sqrt{(d x)^{2}+(d y)^{2}}}{d t}=\frac{\left(\sqrt{1 \pm y^{2}}\right)}{d t} d x v = 2 g y = d t d s = d t ( d x ) 2 + ( d y ) 2 = d t ( 1 ± y 2 ) d x
我们可以写出质点自A沿曲线y(x)自由滑至B点所需的时间为:T [ y ] = ∫ x A x B d t = ∫ x A x B 1 + y ′ 2 2 g y T[y]=\int_{x_A}^{x_B}dt=\int_{x_A}^{x_B}\frac{\sqrt{1+y^{\prime 2}}}{\sqrt{2gy}} T [ y ] = ∫ x A x B d t = ∫ x A x B 2 g y 1 + y ′ 2 可知T的值与曲线的形状有关,T是y的泛函,则我们希望δ T = 0 \delta T=0 δ T = 0
变分问题的欧拉方程
对于F [ y ] = ∫ x 1 x 2 L ( x , y , y ′ ) d x F[y]=\int_{x_{1}}^{x_{2}} L\left(x, y, y^{\prime}\right) d x F [ y ] = ∫ x 1 x 2 L ( x , y , y ′ ) d x 我们要找到一个合适的δ T = 0 \delta T=0 δ T = 0 使得δ F = 0 \delta F=0 δ F = 0 。对于固定边界的问题有
δ y ∣ x = x 1 = δ y ∣ x = x 2 = 0 \left.\delta y\right|_{x=x_{1}}=\left.\delta y\right|_{x=x_{2}}=0 δ y ∣ x = x 1 = δ y ∣ x = x 2 = 0
我们由已知的变分运算法则得到
δ y = y ˉ ( x ) − y ( x ) , ( δ y ) ′ = δ y ′ = y ˉ ′ ( x ) − y ′ ( x ) \delta y=\bar{y}(x)-y(x),(\delta y)^{\prime}=\delta y^{\prime}=\bar{y}^{\prime}(x)-y^{\prime}(x) δ y = y ˉ ( x ) − y ( x ) , ( δ y ) ′ = δ y ′ = y ˉ ′ ( x ) − y ′ ( x )
δ F = F [ y ˉ ( x ) ] − F [ y ( x ) ] \delta F=F[\bar{y}(x)]-F[y(x)] δ F = F [ y ˉ ( x ) ] − F [ y ( x ) ]
于是有
F [ y ˉ ( x ) ] = ∫ x 1 x 2 L ( x , y ˉ , y ˉ ′ ) d x = ∫ x 1 x 2 L ( x , y + δ y , y ′ + δ y ′ ) d x F[\bar{y}(x)]=\int_{x_{1}}^{x_{2}} L\left(x, \bar{y}, \bar{y}^{\prime}\right) d x=\int_{x_{1}}^{x_{2}} L\left(x, y+\delta y, y^{\prime}+\delta y^{\prime}\right) d x F [ y ˉ ( x ) ] = ∫ x 1 x 2 L ( x , y ˉ , y ˉ ′ ) d x = ∫ x 1 x 2 L ( x , y + δ y , y ′ + δ y ′ ) d x
对上式做Taylor展开有
F [ y ˉ ( x ) ] = ∫ x 1 x 2 L ( x , y , y ′ ) d x + ∫ x 1 x 2 ( ∂ L ∂ y δ y + ∂ L ∂ y ′ δ y ′ ) d x + F[\bar{y}(x)]=\int_{x_{1}}^{x_{2}} L\left(x, y, y^{\prime}\right) d x+\int_{x_{1}}^{x_{2}}\left(\frac{\partial L}{\partial y} \delta y+\frac{\partial L}{\partial y^{\prime}} \delta y^{\prime}\right) d x+ F [ y ˉ ( x ) ] = ∫ x 1 x 2 L ( x , y , y ′ ) d x + ∫ x 1 x 2 ( ∂ y ∂ L δ y + ∂ y ′ ∂ L δ y ′ ) d x + higher terms
下面计算变分
δ F = ∫ x 1 x 2 L ( x , y ˉ , y ˉ ′ ) d x − ∫ x 1 x 2 L ( x , y , y ′ ) d x = ∫ x 1 x 2 ( ∂ L ∂ y δ y + ∂ L ∂ y ′ δ y ′ ) d x \delta F=\int_{x_{1}}^{x_{2}} L\left(x, \bar{y}, \bar{y}^{\prime}\right) d x-\int_{x_{1}}^{x_{2}} L\left(x, y, y^{\prime}\right) d x=\int_{x_{1}}^{x_{2}}\left(\frac{\partial L}{\partial y} \delta y+\frac{\partial L}{\partial y^{\prime}} \delta y^{\prime}\right) d x δ F = ∫ x 1 x 2 L ( x , y ˉ , y ˉ ′ ) d x − ∫ x 1 x 2 L ( x , y , y ′ ) d x = ∫ x 1 x 2 ( ∂ y ∂ L δ y + ∂ y ′ ∂ L δ y ′ ) d x
其中
∫ x 1 x 2 ∂ L ∂ y ′ δ y ′ d x = ∫ x 1 x 2 ∂ L ∂ y ′ δ ( d y d x ) d x = ∫ x 1 x 2 ∂ L ∂ y ′ d δ y d x d x \int_{x_{1}}^{x_{2}} \frac{\partial L}{\partial y^{\prime}} \delta y^{\prime} d x=\int_{x_{1}}^{x_{2}} \frac{\partial L}{\partial y^{\prime}} \delta\left(\frac{d y}{d x}\right) d x=\int_{x_{1}}^{x_{2}} \frac{\partial L}{\partial y^{\prime}} \frac{d \delta y}{d x} d x ∫ x 1 x 2 ∂ y ′ ∂ L δ y ′ d x = ∫ x 1 x 2 ∂ y ′ ∂ L δ ( d x d y ) d x = ∫ x 1 x 2 ∂ y ′ ∂ L d x d δ y d x
= ∫ x 1 x 2 d d x ( ∂ L ∂ y ′ δ y ) d x − ∫ x 1 x 2 ( d d x ∂ L ∂ y ′ ) δ y d x =\int_{x_{1}}^{x_{2}} \frac{d}{d x}\left(\frac{\partial L}{\partial y^{\prime}} \delta y\right) d x-\int_{x_{1}}^{x_{2}}\left(\frac{d}{d x} \frac{\partial L}{\partial y^{\prime}}\right) \delta y d x = ∫ x 1 x 2 d x d ( ∂ y ′ ∂ L δ y ) d x − ∫ x 1 x 2 ( d x d ∂ y ′ ∂ L ) δ y d x
= ∂ L ∂ y ′ δ y ∣ x 1 x 2 − ∫ x 1 x 2 ( d d x ∂ L ∂ y ′ ) δ y d x =\left.\frac{\partial L}{\partial y^{\prime}} \delta y\right|_{x_{1}} ^{x_{2}}-\int_{x_{1}}^{x_{2}}\left(\frac{d}{d x} \frac{\partial L}{\partial y^{\prime}}\right) \delta y d x = ∂ y ′ ∂ L δ y ∣ ∣ ∣ x 1 x 2 − ∫ x 1 x 2 ( d x d ∂ y ′ ∂ L ) δ y d x
于是整理得到:δ F = ∫ x 1 x 2 ( ∂ L ∂ y − d d x ∂ L ∂ y ′ ) δ y d x = 0 \delta F=\int_{x_{1}}^{x_{2}}\left(\frac{\partial L}{\partial y}-\frac{d}{d x} \frac{\partial L}{\partial y^{\prime}}\right) \delta y d x=0 δ F = ∫ x 1 x 2 ( ∂ y ∂ L − d x d ∂ y ′ ∂ L ) δ y d x = 0
即有d d x ∂ L ∂ y ′ − ∂ L ∂ y = 0 \frac{d}{d x} \frac{\partial L}{\partial y^{\prime}}-\frac{\partial L}{\partial y}=0 d x d ∂ y ′ ∂ L − ∂ y ∂ L = 0 ,我们称其为欧拉方程
上面的问题中只有一个未知量,对于多元泛函也是差不多的处理方案。对多元泛函δ F = ∫ x 1 x 2 ( ∂ L ∂ y − d d x ∂ L ∂ y ′ ) δ y d x = 0 \delta F=\int_{x_{1}}^{x_{2}}\left(\frac{\partial L}{\partial y}-\frac{d}{d x} \frac{\partial L}{\partial y^{\prime}}\right) \delta y d x=0 δ F = ∫ x 1 x 2 ( ∂ y ∂ L − d x d ∂ y ′ ∂ L ) δ y d x = 0 ,若其边界也是固定的δ y i ∣ x = x 1 = δ y i ∣ x = x 2 = 0 , i = 1 , 2 , ⋯ n \left.\delta y_{i}\right|_{x=x_{1}}=\left.\delta y_{i}\right|_{x=x_{2}}=0, i=1,2, \cdots n δ y i ∣ x = x 1 = δ y i ∣ x = x 2 = 0 , i = 1 , 2 , ⋯ n ,则它的极值问题由Euler方程组决定:
d d x ∂ L ∂ y i ′ − ∂ L ∂ y i = 0 , i = 1 , 2 , ⋯ n \frac{d}{d x} \frac{\partial L}{\partial y_{i}^{\prime}}-\frac{\partial L}{\partial y_{i}}=0, \quad i=1,2, \cdots n d x d ∂ y i ′ ∂ L − ∂ y i ∂ L = 0 , i = 1 , 2 , ⋯ n
哈密顿原理的优点
统一、简洁、完美,具有坐标变换的不变性;
具有很强的普适性可推广至无限自由度以及物理学其他领域;是积分形式的变分原理;
可用于创建新的理论。根据假设构造出拉格朗日函数,用哈密顿原理导出运动方程,由实践检验其正确性;
任何理论有一定的适用范围,这里的哈密顿原理的表述方式也并非对于任意力学系统成立,实际上对力学系统内外部的相互作用有一定的限制,要求相互作用可表示为一标量函数,一般来说,物理学通常关心的正是这种体系
哈密顿正则方程
回顾Lagrange力学:以广义坐标为独立变量,运动方程是构型空间中二阶微分方程组— Lagrange方程;在数学上为了处理问题的方便,往往将一个二阶微分方程化为两个一阶微分方程从而在相空间里讨论问题(降阶法)。对于 Lagrange方程,即将s个阶微分方程化为2s个一阶微分方程来处理,比如
q ˙ α = f α ( q , X , t ) , X ˙ α = g α ( q , X , t ) \dot{q}_{\alpha}=f_{\alpha}(q, X, t), \dot{X}_{\alpha}=g_{\alpha}(q, X, t) q ˙ α = f α ( q , X , t ) , X ˙ α = g α ( q , X , t )
Hamilton力学:以广义坐标和广义动量为独立变量,运动方程是相空间 中的一阶微分方程组——Hamilton正则方程:分析力学的第二个理论形式 。
广义动量与相空间
我们首先定义广义动量p α ≡ ∂ L ∂ q ˙ α , α = 1 , 2 , ⋯ , s p_{\alpha} \equiv \frac{\partial L}{\partial \dot{q}_{\alpha}}, \alpha=1,2, \cdots, s p α ≡ ∂ q ˙ α ∂ L , α = 1 , 2 , ⋯ , s 称 q α , p α q_{\alpha}, p_{\alpha} q α , p α 为相互共轩的正则变量。 q α q_{\alpha} q α 描写系统 的位置状态, p α p_{\alpha} p α 描写系统的运动状态。由相互独 立的 ( q α , p α ) \left(q_{\alpha}, p_{\alpha}\right) ( q α , p α ) 组成的2s维空间,才是系统状态 的完整描述,称为相空间(相即状态)。
一般Lagrangian L ( q , q ˙ , t ) L(q, \dot{q}, t) L ( q , q ˙ , t ) 中 q , q ˙ q, \dot{q} q , q ˙ 的地位不对 称,如果通过变换把 q ˙ α \dot{q}_{\alpha} q ˙ α 变为 p α = p α ( q , q ˙ , t ) p_{\alpha}=p_{\alpha}(q, \dot{q}, t) p α = p α ( q , q ˙ , t ) ,
即做反函数 q ˙ α = q ˙ α ( q , p , t ) \dot{q}_{\alpha}=\dot{q}_{\alpha}(q, p, t) q ˙ α = q ˙ α ( q , p , t ) 并代入L,则可以 把运动微分方程降价,这可行吗?
我们并不可以直接代换,我们需要更多的操作!
Legrand变换
法国数学家与天文学家 A.M. LeGendre于1787年间在研究最小曲面的启发下提出如下问题:
有双元函数 f ( x , y ) , f(x, y), f ( x , y ) , 其全微分为 d f ( x , y ) = u d x + v d y d f(x, y)=u d x+v d y d f ( x , y ) = u d x + v d y
我们今天希望能得到一个g ( u , y ) g(u, y) g ( u , y ) 函数,其中x = ∂ g ∂ u , v = − ∂ g ∂ y x=\frac{\partial g}{\partial u}, v=-\frac{\partial g}{\partial y} x = ∂ u ∂ g , v = − ∂ y ∂ g ,满足x = ∂ g ∂ u , v = − ∂ g ∂ y x=\frac{\partial g}{\partial u}, v=-\frac{\partial g}{\partial y} x = ∂ u ∂ g , v = − ∂ y ∂ g ,也就是说x , u x,u x , u 为一个共轭变量对。一个典型的勒让德变换是这样的
d g ( u , y ) = ∂ g ∂ u d u + ∂ g ∂ y d y d g(u, y)=\frac{\partial g}{\partial u} d u+\frac{\partial g}{\partial y} d y d g ( u , y ) = ∂ u ∂ g d u + ∂ y ∂ g d y
其中有g ( u , y ) = x u − f ( x , y ) = x ( u , y ) u − f ( x ( u , y ) , y ) g(u, y)=x u-f(x, y)=x(u, y) u-f(x(u, y), y) g ( u , y ) = x u − f ( x , y ) = x ( u , y ) u − f ( x ( u , y ) , y ) ,也就是说函数f , g f,g f , g 的微分关系构成互逆关系。
对于全变量的Legrand变换,独立变量( x , y ) → ( u , v ) (x, y) \rightarrow(u, v) ( x , y ) → ( u , v ) ,如下表示
h ( u , v ) = x u + y v − f ( x , y ) h(u, v)=x u+y v-f(x, y) h ( u , v ) = x u + y v − f ( x , y )
其中读者可以自行证明有h ( u , v ) = x u + y v − f ( x , y ) h(u, v)=x u+y v-f(x, y) h ( u , v ) = x u + y v − f ( x , y ) ,即有x = ∂ h ∂ u , y = ∂ h ∂ v x=\frac{\partial h}{\partial u}, y=\frac{\partial h}{\partial v} x = ∂ u ∂ h , y = ∂ v ∂ h
Legrand变换的核心是针对“对偶(共轭)变量对”换变量,建立两组变量及对应的函数之间的变换关系。这两组变量和函数同等地描述系统的性质;变换后新函数的极值点就是原来的极值点。
单变量 Legrand变换的几何意义:用斜率与截距来描述同一条曲线。
热力学的拉格朗日变换
我们由热力学第一定律d U ( S , V ) = T d S − P d V d U(S, V)=T d S-P d V d U ( S , V ) = T d S − P d V ,有
{ T = ( ∂ U ∂ S ) V P = − ( ∂ U ∂ V ) S \left\{\begin{array}{l}T=\left(\frac{\partial U}{\partial S}\right)_{V} \\ P=-\left(\frac{\partial U}{\partial V}\right)_{S}\end{array}\right. { T = ( ∂ S ∂ U ) V P = − ( ∂ V ∂ U ) S
于是可以得到
Helmholtz自由能 f ( T , V ) = U − T S \quad f(T, V)=U-T S f ( T , V ) = U − T S
焓 h ( S , P ) = U + P V \quad h(S, P)=U+P V h ( S , P ) = U + P V
Gibbs自由能 g ( T , P ) = U − T S + P V \quad g(T, P)=U-T S+P V g ( T , P ) = U − T S + P V
多元函数和Legrand变换关系
旧系统
新系统
变量
x i , α i x_{i}, \alpha_{i} x i , α i
G ( y i , α i ) G\left(y_{i}, \alpha_{i}\right) G ( y i , α i )
函数
F ( x i , α i ) F\left(x_{i}, \alpha_{i}\right) F ( x i , α i )
G ( y i , α i ) G\left(y_{i}, \alpha_{i}\right) G ( y i , α i )
变换
G ( y i , α i ) G\left(y_{i}, \alpha_{i}\right) G ( y i , α i )
y i = ∂ F ∂ x i , β i = ∂ F ∂ α i y_{i}=\frac{\partial F}{\partial x_{i}}, \beta_{i}=\frac{\partial F}{\partial \alpha_{i}} y i = ∂ x i ∂ F , β i = ∂ α i ∂ F
G = ∑ i x i y i − F G=\sum_{i} x_{i} y_{i}-F G = ∑ i x i y i − F
哈密顿正则方程
一般Lagrangian L ( q , q ˙ , t ) L(q, \dot{q}, t) L ( q , q ˙ , t ) 中,通过Legrand变换,把 q ˙ α \dot{q}_{\alpha} q ˙ α 变为 p α , p_{\alpha}, p α , 即引入新的函数
H ( p , q , t ) = ∑ α = 1 s q ˙ α p α − L ( q ˙ , q , t ) H(p, q, t)=\sum_{\alpha=1}^{s} \dot{q}_{\alpha} p_{\alpha}-L(\dot{q}, q, t) H ( p , q , t ) = ∑ α = 1 s q ˙ α p α − L ( q ˙ , q , t )
其中 q ˙ α = q ˙ α ( q , p , t ) , \dot{q}_{\alpha}=\dot{q}_{\alpha}(q, p, t) \quad, q ˙ α = q ˙ α ( q , p , t ) , 这样广义能量 H的自变量 变为 ( q , p , t ) (q, p, t) ( q , p , t ) , 称为哈密顿量(Hamiltonian),是描述系统动力学性质的新的特性函数。
原先,我们有p α ≡ ∂ L ∂ q ˙ α p_{\alpha} \equiv \frac{\partial L}{\partial \dot{q}_{\alpha}} p α ≡ ∂ q ˙ α ∂ L ,由 Legrand变换关系:
q ˙ α = ∂ H ∂ p α , ∂ H ∂ q α = − ∂ L ∂ q α , ∂ H ∂ t = − ∂ L ∂ t \dot{q}_{\alpha}=\frac{\partial H}{\partial p_{\alpha}}, \frac{\partial H}{\partial q_{\alpha}}=-\frac{\partial L}{\partial q_{\alpha}}, \frac{\partial H}{\partial t}=-\frac{\partial L}{\partial t} q ˙ α = ∂ p α ∂ H , ∂ q α ∂ H = − ∂ q α ∂ L , ∂ t ∂ H = − ∂ t ∂ L
即有p ˙ α = d d t ∂ L ∂ q ˙ α = ∂ L ∂ q α = − ∂ H ∂ q α \dot{p}_{\alpha}=\frac{d}{d t} \frac{\partial L}{\partial \dot{q}_{\alpha}}=\frac{\partial L}{\partial q_{\alpha}}=-\frac{\partial H}{\partial q_{\alpha}} p ˙ α = d t d ∂ q ˙ α ∂ L = ∂ q α ∂ L = − ∂ q α ∂ H
同时有∂ H ∂ t = − ∂ L ∂ t \frac{\partial H}{\partial t}=-\frac{\partial L}{\partial t} ∂ t ∂ H = − ∂ t ∂ L ,这说明了L不显含t,则H不显含t; L+H肯定不显含t。
注明
正则“ Canonica”一词由 Jacobi于1837年采用,英语中词头 Canon有“正经”的含义,该词对西方人有神圣感;
H量是2s+1个变量, ( q , p , t ) (q, p, t) \quad ( q , p , t ) 的函数 ; 最大优点是 在H中 ( q , p ) (q, p) \quad ( q , p ) 处于平等的地位; H的物理含义是广义能量, H = T 2 − T 0 + V , H=T_{2}-T_{0}+V, H = T 2 − T 0 + V , 在稳定约束下 H = T 2 + V H=T_{2}+V H = T 2 + V
3.“正则方程”更对称,更简洁,更优美,是由2s个阶微分方程组成的方程组,而“ Lagrange方程”是由s个二阶微分方程组成。在“正则方程”中尽管方程的数目成倍增加,但方程的阶数却以相同的比例减少.“正则方程”的更重要的优越性却是方程组之间的对称性(辛对称);这种对称性正是它在物理学其他领域例如量子力学和统计力学中得到广泛应用的理论基础,从这一意义上来说,“ Hamilton力学”不同于“ Lagrange力学”之处,就在于它更侧重于理论体系的构建。
正则方程的“辛”形式
引入2s维空间中的列矢量η ^ = ( q p ) = ( q 1 ⋮ q s p 1 ⋮ p s ) \begin{array}{c}\hat{\eta}=\left(\begin{array}{c}q \\ p\end{array}\right)= & \left(\begin{array}{l}q_{1} \\ \vdots \\ q_{s} \\ p_{1} \\ \vdots \\ p_{s}\end{array}\right)\end{array} η ^ = ( q p ) = ⎝ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎛ q 1 ⋮ q s p 1 ⋮ p s ⎠ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎞ ,则正则方程可以写成矩阵形式(辛形式)
η ^ ˙ = J ^ ∂ H ∂ η ^ , J ^ = ( 0 s × s 1 s × s − 1 s × s 0 s × s ) \dot{\hat{\eta}}=\hat{J} \frac{\partial H}{\partial \hat{\eta}}, \quad \hat{J}=\left(\begin{array}{cc} 0_{s \times s} & 1_{s \times s} \\ -1_{s \times s} & 0_{s \times s}\end{array}\right) η ^ ˙ = J ^ ∂ η ^ ∂ H , J ^ = ( 0 s × s − 1 s × s 1 s × s 0 s × s )
辛矩阵J性质:J ^ 2 = − 1 2 s × 2 s , J ^ T = J ^ − 1 = − J ^ , \hat{J}^{2}=-1_{2 s \times 2 s}, \hat{J}^{T}=\hat{J}^{-1}=-\hat{J}, J ^ 2 = − 1 2 s × 2 s , J ^ T = J ^ − 1 = − J ^ , det J ^ = 1 \hat{J}=1 J ^ = 1
正则方程的初积分
若H函数不显含t,即H ( p , q ) , ∂ H ∂ t = 0 H(p, q), \frac{\partial H}{\partial t}=0 H ( p , q ) , ∂ t ∂ H = 0 ,我们有
d H d t = ∂ H ∂ t + ∑ α = 1 s ( ∂ H ∂ q α q ˙ α + ∂ H ∂ p α p ˙ α ) = ∂ H ∂ t = 0 \frac{d H}{d t}=\frac{\partial H}{\partial t}+\sum_{\alpha=1}^{s}\left(\frac{\partial H}{\partial q_{\alpha}} \dot{q}_{\alpha}+\frac{\partial H}{\partial p_{\alpha}} \dot{p}_{\alpha}\right)=\frac{\partial H}{\partial t}=0 d t d H = ∂ t ∂ H + ∑ α = 1 s ( ∂ q α ∂ H q ˙ α + ∂ p α ∂ H p ˙ α ) = ∂ t ∂ H = 0
故有广义能量积分为H ( p , q ) = E = cons H(p, q)=E=\operatorname{cons} H ( p , q ) = E = c o n s
若H函数不显含某个广义坐标,如∂ H ∂ q α = 0 \frac{\partial H}{\partial q_{\alpha}}=0 ∂ q α ∂ H = 0 则由正则方程有p ˙ α = − ∂ H ∂ q α = 0 , p α = c o n s \dot{p}_{\alpha}=-\frac{\partial H}{\partial q_{\alpha}}=0, p_{\alpha}= cons p ˙ α = − ∂ q α ∂ H = 0 , p α = c o n s 此即广义动量守恒。
Routhian及Rouh方程—混合H-L方法
基本的拉朗德变换有时候是很复杂的,有时候我们不需要所有自由度的变换,我们需要方法的改进!
Routhian
Routh于1876年应用循环积分使 Lagrange方程降维!s个自由度的力学系统中,如只要求其中的m个广义速度通过定义关系变换为广义动量,而保留其余s-m个广义速度不变,此时部分变量Legendre变换为(Routhian)
R = ∑ k = 1 m q ˙ k p k − L ( q ˙ , q , t ) = R ( q α , q ˙ l , p k , t ) R=\sum_{k=1}^{m} \dot{q}_{k} p_{k}-L(\dot{q}, q, t)=R\left(q_{\alpha}, \dot{q}_{l}, p_{k}, t\right) R = ∑ k = 1 m q ˙ k p k − L ( q ˙ , q , t ) = R ( q α , q ˙ l , p k , t )
α = 1 , ⋯ , s ; k = 1 , ⋯ , m ; l = m + 1 , ⋯ , s \alpha=1, \cdots, s ; k=1, \cdots, m ; l=m+1, \cdots, s α = 1 , ⋯ , s ; k = 1 , ⋯ , m ; l = m + 1 , ⋯ , s
全部变量拉朗德变换为
H = ∑ l = m + 1 s q ˙ l p l + ∑ k = 1 m q ˙ k p k − L ( q ˙ , q , t ) = R + ∑ l = m + 1 s q ˙ l p L H=\sum_{l=m+1}^{s} \dot{q}_{l} p_{l}+\sum_{k=1}^{m} \dot{q}_{k} p_{k}-L(\dot{q}, q, t)=R+\sum_{l=m+1}^{s} \dot{q}_{l} p_{L} H = ∑ l = m + 1 s q ˙ l p l + ∑ k = 1 m q ˙ k p k − L ( q ˙ , q , t ) = R + ∑ l = m + 1 s q ˙ l p L
Routh方程
对以上变换我们有
{ ∂ R ∂ t = − ∂ L ∂ t = ∂ H ∂ t ∂ R ∂ q α = − ∂ L ∂ q α = ∂ H ∂ q α \left\{\begin{array}{l}\frac{\partial R}{\partial t}=-\frac{\partial L}{\partial t}=\frac{\partial H}{\partial t} \\ \frac{\partial R}{\partial q_{\alpha}}=-\frac{\partial L}{\partial q_{\alpha}}=\frac{\partial H}{\partial q_{\alpha}}\end{array}\right. { ∂ t ∂ R = − ∂ t ∂ L = ∂ t ∂ H ∂ q α ∂ R = − ∂ q α ∂ L = ∂ q α ∂ H { ∂ R ∂ q ˙ l = − ∂ L ∂ q ˙ l = − p l ∂ R ∂ p k = q ˙ k = ∂ H ∂ p k \left\{\begin{array}{l}\frac{\partial R}{\partial \dot{q}_{l}}=-\frac{\partial L}{\partial \dot{q}_{l}}=-p_{l} \\ \frac{\partial R}{\partial p_{k}}=\dot{q}_{k}=\frac{\partial H}{\partial p_{k}}\end{array}\right. { ∂ q ˙ l ∂ R = − ∂ q ˙ l ∂ L = − p l ∂ p k ∂ R = q ˙ k = ∂ p k ∂ H
则由 Lagrange方程有
d d t ∂ R ∂ q ˙ l − ∂ R ∂ q l = 0 , ( l = m + 1 , ⋯ , s ) \frac{d}{d t} \frac{\partial R}{\partial \dot{q}_{l}}-\frac{\partial R}{\partial q_{l}}= 0 ,(l=m+1, \cdots, s) d t d ∂ q ˙ l ∂ R − ∂ q l ∂ R = 0 , ( l = m + 1 , ⋯ , s )
由 Hamilton正则方程有
{ p ˙ k = − ∂ H ∂ q k = − ∂ R ∂ q k q ˙ k = ∂ R ∂ p k \left\{\begin{array}{l}\dot{p}_{k}=-\frac{\partial H}{\partial q_{k}}=-\frac{\partial R}{\partial q_{k}} \\ \dot{q}_{k}=\frac{\partial R}{\partial p_{k}}\end{array}\right. { p ˙ k = − ∂ q k ∂ H = − ∂ q k ∂ R q ˙ k = ∂ p k ∂ R
泊松括号与泊松定狸
Poisson括号的引入
有了正则方程,如何求解正则方程是物理上必须考虑的问题。如果选取的广义坐标使H函数中能多出现些循环坐标,就能在正则方程中多得出一些积分,对微分方程的求解就更有利,否则徼分方程的求解就变得十分困难。而可遗坐标代表着守恒量,如何系统地寻找系统的守恒量? Poisson括号及 Poisson定理就是系统地寻找守恒量的理论方法。与正则变换,H-J方程一起构成 Hamilton力学的主要内容。
寻找守恒量,就是要考察任意一个力学量的时间变化率。设F = F ( q , p , t ) F=F(q, p, t) F = F ( q , p , t ) 代表一力学量,计算其时间变化率,即计算其时间全导数:
d F d t = ∂ F ∂ t + ∑ α = 1 s ( ∂ F ∂ q α q ˙ α + ∂ F ∂ p α p ˙ α ) = ∂ F ∂ t + ∑ α = 1 s ( ∂ F ∂ q α ∂ H ∂ p α − ∂ F ∂ p α ∂ H ∂ q α ) = ∂ F ∂ t + [ F , H ] \begin{aligned} \frac{d F}{d t} &=\frac{\partial F}{\partial t}+\sum_{\alpha=1}^{s}\left(\frac{\partial F}{\partial q_{\alpha}} \dot{q}_{\alpha}+\frac{\partial F}{\partial p_{\alpha}} \dot{p}_{\alpha}\right) \\ &=\frac{\partial F}{\partial t}+\sum_{\alpha=1}^{s}\left(\frac{\partial F}{\partial q_{\alpha}} \frac{\partial H}{\partial p_{\alpha}}-\frac{\partial F}{\partial p_{\alpha}} \frac{\partial H}{\partial q_{\alpha}}\right) \\ &=\frac{\partial F}{\partial t}+[F, H] \end{aligned} d t d F = ∂ t ∂ F + α = 1 ∑ s ( ∂ q α ∂ F q ˙ α + ∂ p α ∂ F p ˙ α ) = ∂ t ∂ F + α = 1 ∑ s ( ∂ q α ∂ F ∂ p α ∂ H − ∂ p α ∂ F ∂ q α ∂ H ) = ∂ t ∂ F + [ F , H ]
定义缩写符号:[ F , H ] ≡ ∑ α = 1 s ( ∂ F ∂ q α ∂ H ∂ p α − ∂ F ∂ p α ∂ H ∂ q α ) [F, H] \equiv \sum_{\alpha=1}^{s}\left(\frac{\partial F}{\partial q_{\alpha}} \frac{\partial H}{\partial p_{\alpha}}-\frac{\partial F}{\partial p_{\alpha}} \frac{\partial H}{\partial q_{\alpha}}\right) [ F , H ] ≡ ∑ α = 1 s ( ∂ q α ∂ F ∂ p α ∂ H − ∂ p α ∂ F ∂ q α ∂ H ) ,我们称其为力学量F与H间的Poisson括号。
我们可以做推广:任意两力学量F F F 与G G G 均是( q , p , t ) (q, p, t) ( q , p , t ) 的函数,则可以定义它们间的Poisson括号为[ F , G ] ≡ ∑ α = 1 s ( ∂ F ∂ q α ∂ G ∂ p α − ∂ F ∂ p α ∂ G ∂ q α ) [F, G] \equiv \sum_{\alpha=1}^{s}\left(\frac{\partial F}{\partial q_{\alpha}} \frac{\partial G}{\partial p_{\alpha}}-\frac{\partial F}{\partial p_{\alpha}} \frac{\partial G}{\partial q_{\alpha}}\right) [ F , G ] ≡ ∑ α = 1 s ( ∂ q α ∂ F ∂ p α ∂ G − ∂ p α ∂ F ∂ q α ∂ G )
其辛形式为[ F , G ] ( q , p ) ≡ ( ∂ F ∂ η ^ ) T J ∂ G ∂ η ^ [F, G]_{(q, p)} \equiv\left(\frac{\partial F}{\partial \hat{\eta}}\right)^{T} J \frac{\partial G}{\partial \hat{\eta}} [ F , G ] ( q , p ) ≡ ( ∂ η ^ ∂ F ) T J ∂ η ^ ∂ G ,采用辛形式的好处可以简化后面一系列定理、性质等的证明。
注意 :这里要求( q , p ) (q, p) ( q , p ) 相互独立,即有∂ p β ∂ q α = 0 , ∂ p β ∂ p α = δ α β , ∂ q β ∂ q α = δ α β \frac{\partial p_{\beta}}{\partial q_{\alpha}}=0, \frac{\partial p_{\beta}}{\partial p_{\alpha}}=\delta_{\alpha \beta}, \frac{\partial q_{\beta}}{\partial q_{\alpha}}=\delta_{\alpha \beta} ∂ q α ∂ p β = 0 , ∂ p α ∂ p β = δ α β , ∂ q α ∂ q β = δ α β
Poisson括号的性质
基本性质
和常数c c c 的运算:[ f , c ] = 0 , [ f , c g ] = c [ f , g ] [f, c]=0, \quad[f, c g]=c[f, g] [ f , c ] = 0 , [ f , c g ] = c [ f , g ]
反对称性:[ f , f ] = 0 , [ f , g ] = − [ g , f ] [f, f]=0, \quad[f, g]=-[g, f] [ f , f ] = 0 , [ f , g ] = − [ g , f ]
分配律:[ f 1 + f 2 , g ] = [ f 1 , g ] + [ f 2 , g ] \left[f_{1}+f_{2}, g\right]=\left[f_{1}, g\right]+\left[f_{2}, g\right] [ f 1 + f 2 , g ] = [ f 1 , g ] + [ f 2 , g ]
[ f , g 1 + g 2 ] = [ f , g 1 ] + [ f , g 2 ] \left[f, g_{1}+g_{2}\right]=\left[f, g_{1}\right]+\left[f, g_{2}\right] [ f , g 1 + g 2 ] = [ f , g 1 ] + [ f , g 2 ]
结合律:[ f 1 f 2 , g ] = f 1 [ f 2 , g ] + [ f 1 , g ] f 2 \left[f_{1} f_{2}, g\right]=f_{1}\left[f_{2}, g\right]+\left[f_{1}, g\right] f_{2} [ f 1 f 2 , g ] = f 1 [ f 2 , g ] + [ f 1 , g ] f 2
[ f , g 1 g 2 ] = g 1 [ f , g 2 ] + [ f , g 1 ] g 2 \left[f, g_{1} g_{2}\right]=g_{1}\left[f, g_{2}\right]+\left[f, g_{1}\right] g_{2} [ f , g 1 g 2 ] = g 1 [ f , g 2 ] + [ f , g 1 ] g 2
基本Poisson括号关系(正则变量满足的关系)
[ p α , p β ] = 0 , [ q α , q β ] = 0 , [ q α , p β ] = δ α β \left[p_{\alpha}, p_{\beta}\right]=0, \quad\left[q_{\alpha}, q_{\beta}\right]=0, \quad\left[q_{\alpha}, p_{\beta}\right]=\delta_{\alpha \beta} [ p α , p β ] = 0 , [ q α , q β ] = 0 , [ q α , p β ] = δ α β
该性质的独特性在于这是唯一体现q和p变量在量的消长方面严格互补的定量制约关系的关系—量子力学不确定关系。由 Poisson括号性质及三个基本 Poisson括号关系可以计算任意两个力学量间的 Poisson括号。
其可导出微商性质:
∂ ∂ t [ f , g ] = [ ∂ f ∂ t , g ] + [ f , ∂ g ∂ t ] \frac{\partial}{\partial t}[f, g]=\left[\frac{\partial f}{\partial t}, g\right]+\left[f, \frac{\partial g}{\partial t}\right] ∂ t ∂ [ f , g ] = [ ∂ t ∂ f , g ] + [ f , ∂ t ∂ g ]
[ q α , f ] = ∂ f ∂ p α [ p α , f ] = − ∂ f ∂ q α \left[q_{\alpha}, f\right]=\frac{\partial f}{\partial p_{\alpha}} \quad\left[p_{\alpha}, f\right]=-\frac{\partial f}{\partial q_{\alpha}} [ q α , f ] = ∂ p α ∂ f [ p α , f ] = − ∂ q α ∂ f
同时我们有Jacobi恒等式:
[ f , [ g , h ] ] + [ g , [ h , f ] ] + [ h , [ f , g ] ] = 0 [f,[g, h]]+[g,[h, f]]+[h,[f, g]]=0 [ f , [ g , h ] ] + [ g , [ h , f ] ] + [ h , [ f , g ] ] = 0
Poisson括号的高级性质
我们对任意力学量F = F ( q , p , t ) F=F(q, p, t) F = F ( q , p , t ) 的时间全导数(EOM)有d F d t = ∂ F ∂ t + [ F , H ] \frac{d F}{d t}=\frac{\partial F}{\partial t}+[F, H] d t d F = ∂ t ∂ F + [ F , H ] ,这可以启发我们
Hamiltonian掌管了系统力学量随时间的演化,是描述系统时间演化的特性函数。
力学量F守恒的充要条件为∂ F ∂ t + [ F , H ] = 0 \frac{\partial F}{\partial t}+[F, H]=0 ∂ t ∂ F + [ F , H ] = 0
如果力学量F不显含时间,则其守恒的充要条件为[ F , H ] = 0 [F, H]=0 [ F , H ] = 0
同时,H也是力学量,故其时间演化为d H d t = ∂ H ∂ t + [ H , H ] = ∂ H ∂ t \frac{d H}{d t}=\frac{\partial H}{\partial t}+[H, H]=\frac{\partial H}{\partial t} d t d H = ∂ t ∂ H + [ H , H ] = ∂ t ∂ H ,也就是说H的时间全导数等于其对时间的偏导数。如果H不显含时间t,则其偏导数等于0,同时L也不显含时间,意味着其时间全导数等于0,即广义能量守恒。
另外,正则方程可用泊松括号表示:
q ˙ α = ∂ H ∂ p α ⟹ q ˙ α = d q α d t = [ q α , H ] \dot{q}_{\alpha}=\frac{\partial H }{\partial p _{\alpha}} \Longrightarrow \dot{ q }_{\alpha}=\frac{d q _{\alpha}}{d t}=\left[ q _{\alpha}, H \right] q ˙ α = ∂ p α ∂ H ⟹ q ˙ α = d t d q α = [ q α , H ]
p ˙ α = − ∂ H ∂ q α ⟹ p ˙ α = d p α d t = [ p α , H ] \dot{p}_{\alpha}=-\frac{\partial H }{\partial q _{\alpha}} \Longrightarrow \dot{p}_{\alpha}=\frac{d p_{\alpha}}{d t}=\left[p_{\alpha}, H \right] p ˙ α = − ∂ q α ∂ H ⟹ p ˙ α = d t d p α = [ p α , H ]
注意到,用泊松括号表示的正则方程完全对称!
我们还可以引出Poisson定理:如果f f f 和g g g 是系统守恒量,则[ f , g ] [f, g] [ f , g ] 也是系统守恒量,证明如下:
d d t [ f , g ] = ∂ ∂ t [ f , g ] + [ [ f , g ] , H ] = [ ∂ f ∂ t , g ] + [ f , ∂ g ∂ t ] − [ g , [ f , H ] ] + [ f , [ g , H ] ] = [ ∂ f ∂ t , g ] + [ [ f , H ] , g ] + [ f , ∂ g ∂ t ] + [ f , [ g , H ] ] = [ ∂ f ∂ t + [ f , H ] , g ] + [ f , ∂ g ∂ t + [ g , H ] ] = [ d f d t , g ] + [ f , d g d t ] = 0 \begin{aligned} \frac{d}{d t}[f, g] &=\frac{\partial}{\partial t}[f, g]+[[f, g], H] \\ &=\left[\frac{\partial f}{\partial t}, g\right]+\left[f, \frac{\partial g}{\partial t}\right]-[g,[f, H]]+[f,[g, H]] \\ &=\left[\frac{\partial f}{\partial t}, g\right]+[[f, H], g]+\left[f, \frac{\partial g}{\partial t}\right]+[f,[g, H]] \\ &=\left[\frac{\partial f}{\partial t}+[f, H], g\right]+\left[f, \frac{\partial g}{\partial t}+[g, H]\right]=\left[\frac{d f}{d t}, g\right]+\left[f, \frac{d g}{d t}\right]=0 \end{aligned} d t d [ f , g ] = ∂ t ∂ [ f , g ] + [ [ f , g ] , H ] = [ ∂ t ∂ f , g ] + [ f , ∂ t ∂ g ] − [ g , [ f , H ] ] + [ f , [ g , H ] ] = [ ∂ t ∂ f , g ] + [ [ f , H ] , g ] + [ f , ∂ t ∂ g ] + [ f , [ g , H ] ] = [ ∂ t ∂ f + [ f , H ] , g ] + [ f , ∂ t ∂ g + [ g , H ] ] = [ d t d f , g ] + [ f , d t d g ] = 0
用泊松括号可以从已知守恒量构造新的守恒量!当然需要注意的是构造的守恒量不一定是独立的全新的守恒量否则会导致理论的内在矛盾,存在任意多独立守恒量。
由上面的性质我们可以很轻易地使用Poisson括号证明角动量之间的关系[ L x , L y ] = L z \left[L_{x}, L_{y}\right]=L_{z} [ L x , L y ] = L z ,[ L y , L z ] = L x \left[L_{y}, L_{z}\right]=L_{x} [ L y , L z ] = L x ,[ L z , L x ] = L y \left[L_{z}, L_{x}\right]=L_{y} [ L z , L x ] = L y ,由 Poisson定理可知,只要有两个相互垂直的角动量分量是运动积分,则上式表明其第三个垂直分量必是运动积分。说明与线动量不同,角动量三个分量在运动学上不是完全独立的(是指在变化上的不独立),这符合如下事实:两个垂直轴上的转动可以合成绕第三个垂直轴上的转动第二,更为重要的是上式表明,与线动量不同,角动量三个分量不能同时取作为正则变量,因为这违反基本 Poisson括号公理。这种限制在量子力学上的反映:不存在L x , L y , L z L_{x}, L_{y}, L_{z} L x , L y , L z 同时有确切取值的状态(L = 0 L=0 L = 0 除外)。这说明:正则描述更深刻地描述了自然。
动力学问题的纯代数解法或者称为算子解法量子代数的前身是Taylor展开的 Poisson括号表示:若f ( q , p ) f(q, p) f ( q , p ) 不显含t t t ,则有d f d t = [ f , H ] , d 2 f d t 2 = [ [ f , H ] , H ] , ⋯ \frac{d f}{d t}=[f, H], \frac{d^{2} f}{d t^{2}}=[[f, H], H], \cdots d t d f = [ f , H ] , d t 2 d 2 f = [ [ f , H ] , H ] , ⋯ 。若定义算符D g f = [ f , g ] D_{g} f=[f, g] D g f = [ f , g ] ,则有D g f = [ f , g ] , D g 0 f = f , D g 2 f = [ [ f , g ] ] g ] \left.\left.D_{g} f=[f, g], \quad D_{g}^{0} f=f, D_{g}^{2} f=[[f,g]\right] g\right] D g f = [ f , g ] , D g 0 f = f , D g 2 f = [ [ f , g ] ] g ] ,于是有
f ( t ) = D H 0 f ∣ 0 + D H 1 f ∣ 0 t + D H 2 f ∣ 0 1 2 ! t 2 + ⋯ = ∑ n = 0 ∞ D H n f ∣ 0 1 n ! t n = e D H t f ∣ 0 f(t)=\left.D_{H}^{0} f\right|_{0}+\left.D_{H}^{1} f\right|_{0} t+\left.D_{H}^{2} f\right|_{0} \frac{1}{2 !} t^{2}+\cdots=\left.\sum_{n=0}^{\infty} D_{H}^{n} f\right|_{0} \frac{1}{n !} t^{n}=\left.e^{D_{H} t} f\right|_{0} f ( t ) = D H 0 f ∣ ∣ 0 + D H 1 f ∣ ∣ 0 t + D H 2 f ∣ ∣ 0 2 ! 1 t 2 + ⋯ = ∑ n = 0 ∞ D H n f ∣ 0 n ! 1 t n = e D H t f ∣ ∣ 0 注意到最后的表达是很简洁的。
正则变换
选取不同的广义坐标,所得的微分方程的形式不同,求解方程的难易程度不同。如果选取的广义坐标使H函数中能多出现一些循环坐标,就能在正则方程中多得出一些积分,对微分方程的求解就更有利,否则微分方程的求解就变得十分困难。但有无可遗坐标是与坐标系的选择有关的。因此,选取广义坐标是理论力学中最富技术性的环节。
正则变换即是一种坐标变换,使尽可能多的广义坐标及广义动量变为可遗坐标,同时要保证正则方程形式不变。
正则变换理论的目的就在于保证正则方程形式不变条件下寻找到尽可能多的广义坐标和广义动量为可遗坐标的理论。
原则上,对可积系统总可使所有广义坐标全部是循环坐标,从而使正则方程求解大大简化。
点变换与接触变换
L ( q , q ˙ , t ) L(q, \dot{q}, t) L ( q , q ˙ , t ) 中( q , q ˙ ) (q, \dot{q}) ( q , q ˙ ) 为动力学的独立变量,而从运动学角度看,两者又是有关联的,所以我们称它们为在“动力学”上的独立变量,而在“运动学”上是相关变量。若对 Lagrange方程中的广义坐标作坐标变换:
q ~ α = q ~ α ( q β ) ⇒ q ~ ˙ α = q ~ ˙ α ( q ˙ β , q β ) \tilde{q}_{\alpha}=\tilde{q}_{\alpha}\left(q_{\beta}\right) \Rightarrow \dot{\tilde{q}}_{\alpha}=\dot{\tilde{q}}_{\alpha}\left(\dot{q}_{\beta}, q_{\beta}\right) q ~ α = q ~ α ( q β ) ⇒ q ~ ˙ α = q ~ ˙ α ( q ˙ β , q β )
L ( q , q ˙ , t ) ⇒ L ~ ( q ~ , q ~ ˙ , t ) L(q, \dot{q}, t) \Rightarrow \tilde{L}(\tilde{q}, \dot{\tilde{q}}, t) L ( q , q ˙ , t ) ⇒ L ~ ( q ~ , q ~ ˙ , t )
d d t ∂ L ∂ q ˙ α − ∂ L ∂ q α = 0 ⇒ d d t ∂ L ~ ∂ q ˙ α − ∂ L ~ ∂ q ~ ˙ α = 0 \frac{d}{d t} \frac{\partial L}{\partial \dot{q}_{\alpha}}-\frac{\partial L}{\partial q_{\alpha}}=0 \Rightarrow \frac{d}{d t} \frac{\partial \tilde{L}}{\partial \dot{q}_{\alpha}}-\frac{\partial \tilde{L}}{\partial \dot{\tilde{q}}_{\alpha}}=0 d t d ∂ q ˙ α ∂ L − ∂ q α ∂ L = 0 ⇒ d t d ∂ q ˙ α ∂ L ~ − ∂ q ~ ˙ α ∂ L ~ = 0
即 Lagrange方程形式保持不变,此种变换称为Lagrange方程的“点变换”。
对一力学系统而言,描述它的 Lagrange方程是不会随此系统广义坐标的选择不同而改变其形式的。因此点变换不会影响L氏方程的结构。
然而,对 Hamiltonian H = H ( q , p , t ) H=H(q, p, t) H = H ( q , p , t ) ,描述系统的坐标变为 Lagrangian的两倍,为( q , p ) (q, p) ( q , p ) ,故相应的坐标变换为“接触变换”{ q ~ α = q ~ α ( q , p , t ) p ~ α = p ~ α ( q , p , t ) \left\{\begin{array}{l}\tilde{q}_{\alpha}=\tilde{q}_{\alpha}(q, p, t) \\ \tilde{p}_{\alpha}=\tilde{p}_{\alpha}(q, p, t)\end{array}\right. { q ~ α = q ~ α ( q , p , t ) p ~ α = p ~ α ( q , p , t )
H ( q , p , t ) ⇒ H ~ ( q ~ , p ~ , t ) H(q, p, t) \Rightarrow \tilde{H}(\tilde{q}, \tilde{p}, t) H ( q , p , t ) ⇒ H ~ ( q ~ , p ~ , t )
此变换相对“点变换”任意性大,可选择的范围更多,但变换不能保证正则方程形式不变,即
{ q ˙ α = ∂ H ∂ p α p ˙ α = − ∂ H ∂ q α \left\{\begin{array}{l}\dot{q}_{\alpha}=\frac{\partial H}{\partial p_{\alpha}} \\ \dot{p}_{\alpha}=-\frac{\partial H}{\partial q_{\alpha}}\end{array}\right. { q ˙ α = ∂ p α ∂ H p ˙ α = − ∂ q α ∂ H 不能保证{ q ˙ α = ∂ H ∂ p α p ˙ α = − ∂ H ∂ q α \left\{\begin{array}{l}\dot{q}_{\alpha}=\frac{\partial H}{\partial p_{\alpha}} \\ \dot{p}_{\alpha}=-\frac{\partial H}{\partial q_{\alpha}}\end{array}\right. { q ˙ α = ∂ p α ∂ H p ˙ α = − ∂ q α ∂ H
我们思考是否存在满足正则方程的“接触变换”呢?如何能自动得到满足正则方程的“接触变换”呢?
正则变换理论正是一种寻找能自动满足正则方程形式不变的接触变换!
正则变换
定义 :正则变换就是既要保留 Hamilton力学的已有成果,如 Hamiltonian与 Lagrangian之间的Legrand变换关系和 Hamilton正则方程,同时又能增加可遗坐标的一种特殊接触变换
H ( q , p , t ) ⇒ H ~ ( q ~ , p ~ , t ) H ( p , q , t ) = ∑ α = 1 s q ˙ α p α − L ( q ˙ , q , t ) ⇒ H ~ ( p ~ , q ~ , t ) = ∑ α = 1 s q ~ ˙ α p ~ α − L ~ ( q ~ ˙ , q ~ , t ) \begin{aligned} H(q, p, t) & \Rightarrow \tilde{H}(\tilde{q}, \tilde{p}, t) \\ H(p, q, t) &=\sum_{\alpha=1}^{s} \dot{q}_{\alpha} p_{\alpha}-L(\dot{q}, q, t) \\ \Rightarrow \tilde{H}(\tilde{p}, \tilde{q}, t) &=\sum_{\alpha=1}^{s} \dot{\tilde{q}}_{\alpha} \tilde{p}_{\alpha}-\tilde{L}(\dot{\tilde{q}}, \tilde{q}, t) \end{aligned} H ( q , p , t ) H ( p , q , t ) ⇒ H ~ ( p ~ , q ~ , t ) ⇒ H ~ ( q ~ , p ~ , t ) = α = 1 ∑ s q ˙ α p α − L ( q ˙ , q , t ) = α = 1 ∑ s q ~ ˙ α p ~ α − L ~ ( q ~ ˙ , q ~ , t )
我们使用生成函数来得到符合以上性质的变换,回到开始的开始——生成函数:要求变换前后均满足Hamilton作用量原理:
δ S = δ ∫ L d t = δ ∫ [ ∑ α = 1 s q ˙ α p α − H ( p , q , t ) ] d t = 0 δ S ~ = δ ∫ L ~ d t = δ ∫ [ ∑ α = 1 s q ~ ˙ α p ~ α − H ~ ( p ~ , q ~ , t ) ] d t = 0 \begin{aligned} \delta S &=\delta \int L d t=\delta \int\left[\sum_{\alpha=1}^{s} \dot{q}_{\alpha} p_{\alpha}-H(p, q, t)\right] d t=0 \\ \delta \tilde{S} &=\delta \int \tilde{L} d t=\delta \int\left[\sum_{\alpha=1}^{s} \dot{\tilde{q}}_{\alpha} \tilde{p}_{\alpha}-\tilde{H}(\tilde{p}, \tilde{q}, t)\right] d t=0 \end{aligned} δ S δ S ~ = δ ∫ L d t = δ ∫ [ α = 1 ∑ s q ˙ α p α − H ( p , q , t ) ] d t = 0 = δ ∫ L ~ d t = δ ∫ [ α = 1 ∑ s q ~ ˙ α p ~ α − H ~ ( p ~ , q ~ , t ) ] d t = 0
由 Hamilton原理可以得到 Lagrange方程;将Hamiltonian与 Lagrangian之间的 Legrand变换代入其中,又可以得到 Hamilton正则方程,必有:
∑ α = 1 s q ˙ α p α − H ( p , q , t ) = [ ∑ α = 1 s q ~ ˙ α p ~ α − H ~ ( p ~ , q ~ , t ) ] + d F d t \sum_{\alpha=1}^{s} \dot{q}_{\alpha} p_{\alpha}-H(p, q, t)=\left[\sum_{\alpha=1}^{s} \dot{\tilde{q}}_{\alpha} \tilde{p}_{\alpha}-\tilde{H}(\tilde{p}, \tilde{q}, t)\right]+\frac{d F}{d t} ∑ α = 1 s q ˙ α p α − H ( p , q , t ) = [ ∑ α = 1 s q ~ ˙ α p ~ α − H ~ ( p ~ , q ~ , t ) ] + d t d F
即有
∑ α = 1 s p α d q α − ∑ α = 1 s p ~ α d q ~ α + [ H ~ ( p ~ , q ~ , t ) − H ( p , q , t ) ] d t = d F \sum_{\alpha=1}^{s} p_{\alpha} d q_{\alpha}-\sum_{\alpha=1}^{s} \tilde{p}_{\alpha} d \tilde{q}_{\alpha}+[\tilde{H}(\tilde{p}, \tilde{q}, t)-H(p, q, t)] d t=d F ∑ α = 1 s p α d q α − ∑ α = 1 s p ~ α d q ~ α + [ H ~ ( p ~ , q ~ , t ) − H ( p , q , t ) ] d t = d F
我们称其为Pfaff型方程,是微分形式的正则变换条件(F为任意形式的函数)。
注意式中H ~ , F \tilde{H}, F H ~ , F 都未确定,所以此条件较易满足。可先设定函数F F F ,据此确定正则变换关系,最后定出H H H 。
函数F F F 完全决定了整个正则变换,故称为正则变换的生成函数(母函数)。它必须同时与新老变量有关。
对Pfaff型方程进行不同的 Legrand变换,可得到4类不同类型的生成函数和正则变换(根据独立变量的选择方式不同而不同):
第一类生成函数及生成方程:Pfaff方程中母函数F F F 自变量是( q , q ~ , t ) (q, \tilde{q}, t) ( q , q ~ , t ) ,记为F 1 ( q , q ~ , t ) F_{1}(q, \tilde{q}, t) F 1 ( q , q ~ , t ) ,于是我们有:
∑ α = 1 S p α d q α − ∑ α = 1 p ~ α d q ~ α + ( H ~ − H ) d t = d F 1 ( q , q ~ , t ) \sum_{\alpha=1}^{S} p_{\alpha} d q_{\alpha}-\sum_{\alpha=1} \tilde{p}_{\alpha} d \tilde{q}_{\alpha}+(\tilde{H}-H) d t=d F_{1}(q, \tilde{q}, t) ∑ α = 1 S p α d q α − ∑ α = 1 p ~ α d q ~ α + ( H ~ − H ) d t = d F 1 ( q , q ~ , t )
其生成方程为{ p α = ∂ F 1 ∂ q α p ~ α = − ∂ F 1 ∂ q ~ α H ~ = H + ∂ F 1 ∂ t \left\{\begin{array}{l}p_{\alpha}=\frac{\partial F_{1}}{\partial q_{\alpha}} \\ \tilde{p}_{\alpha}=-\frac{\partial F_{1}}{\partial \tilde{q}_{\alpha}} \\ \tilde{H}=H+\frac{\partial F_{1}}{\partial t}\end{array}\right. ⎩ ⎪ ⎨ ⎪ ⎧ p α = ∂ q α ∂ F 1 p ~ α = − ∂ q ~ α ∂ F 1 H ~ = H + ∂ t ∂ F 1
可先从生成方程前两式中解出( q , p ) (q, p) ( q , p ) 为( q , p ) (q, p) ( q , p ) 的函数,再代入第三式,即可得到H ~ ( q ~ , p ~ , t ) \tilde{H}(\tilde{q}, \tilde{p}, t) H ~ ( q ~ , p ~ , t ) ,从而得到新的正则方程。
第二类生成函数及生成方程:Pfa方程中母函数F自变量是( q , p ~ , t ) (q, \tilde{p}, t) ( q , p ~ , t ) ,记为F 2 ( q , p ~ , t ) F_{2}(q, \tilde{p}, t) F 2 ( q , p ~ , t ) 。可对第一类生成函数进行部分变量 Legrand变换
d ( F 1 + ∑ α q ~ α p ~ α ) = ∑ α = 1 s p α d q α + ∑ α = 1 s q ~ α d p ~ α + ( H ~ − H ) d t = d F 2 ( q , q ~ , t ) d\left(F_{1}+\sum_{\alpha} \tilde{q}_{\alpha} \tilde{p}_{\alpha}\right)=\sum_{\alpha=1}^{s} p_{\alpha} d q_{\alpha}+\sum_{\alpha=1}^{s} \tilde{q}_{\alpha} d \tilde{p}_{\alpha}+(\tilde{H}-H) d t=d F_{2}(q, \tilde{q}, t) d ( F 1 + ∑ α q ~ α p ~ α ) = ∑ α = 1 s p α d q α + ∑ α = 1 s q ~ α d p ~ α + ( H ~ − H ) d t = d F 2 ( q , q ~ , t )
有生成方程{ p α = ∂ F 2 ∂ q α q ~ α = ∂ F 2 ∂ p ~ α H ~ = H + ∂ F 2 ∂ t \left\{\begin{array}{l}p_{\alpha}=\frac{\partial F_{2}}{\partial q_{\alpha}} \\ \tilde{q}_{\alpha}=\frac{\partial F_{2}}{\partial \tilde{p}_{\alpha}} \\ \tilde{H}=H+\frac{\partial F_{2}}{\partial t}\end{array}\right. ⎩ ⎪ ⎨ ⎪ ⎧ p α = ∂ q α ∂ F 2 q ~ α = ∂ p ~ α ∂ F 2 H ~ = H + ∂ t ∂ F 2
第三类生成函数及生成方程:Pfa方程中母函数F自变量是( p , q ~ , t ) (p, \tilde{q}, t) ( p , q ~ , t ) ,记为F 3 ( p , q ~ , t ) F_{3}(p, \tilde{q}, t) F 3 ( p , q ~ , t ) 。对第一类生成函数进行部分变量 Legrand变换
d ( F 1 − ∑ α q α p α ) = − ∑ α = 1 s q α d p α − ∑ α = 1 s p ~ α d q ~ α + ( H ~ − H ) d t = d F 3 ( p , q ~ , t ) d\left(F_{1}-\sum_{\alpha} q_{\alpha} p_{\alpha}\right)=-\sum_{\alpha=1}^{s} q_{\alpha} d p_{\alpha}-\sum_{\alpha=1}^{s} \tilde{p}_{\alpha} d \tilde{q}_{\alpha}+(\tilde{H}-H) d t=d F_{3}(p, \tilde{q}, t) d ( F 1 − ∑ α q α p α ) = − ∑ α = 1 s q α d p α − ∑ α = 1 s p ~ α d q ~ α + ( H ~ − H ) d t = d F 3 ( p , q ~ , t )
有生成方程{ q α = − ∂ F 3 ∂ p α p ~ α = − ∂ F 3 ∂ q ~ α H ~ = H + ∂ F 3 ∂ t \left\{\begin{array}{l}q_{\alpha}=-\frac{\partial F_{3}}{\partial p_{\alpha}} \\ \tilde{p}_{\alpha}=-\frac{\partial F_{3}}{\partial \tilde{q}_{\alpha}} \\ \tilde{H}=H+\frac{\partial F_{3}}{\partial t}\end{array}\right. ⎩ ⎪ ⎨ ⎪ ⎧ q α = − ∂ p α ∂ F 3 p ~ α = − ∂ q ~ α ∂ F 3 H ~ = H + ∂ t ∂ F 3
第四类生成函数及生成方程:Pfaff-方程中母函数F自变量是( p , p ~ , t ) (p, \tilde{p}, t) ( p , p ~ , t ) ,记为F 4 ( p , q ~ , t ) F_{4}(p, \tilde{q}, t) F 4 ( p , q ~ , t ) 。对一第类生成函数进行全变量 Legrand变换
d ( F 1 − ∑ α q α p α + ∑ α p ~ α q ~ α ) = − ∑ α = 1 s q α d p α + ∑ α = 1 s q ~ α d p ~ α + ( H ~ − H ) d t = d F 4 ( p , p ~ , t ) d\left(F_{1}-\sum_{\alpha} q_{\alpha} p_{\alpha}+\sum_{\alpha} \tilde{p}_{\alpha} \tilde{q}_{\alpha}\right)=-\sum_{\alpha=1}^{s} q_{\alpha} d p_{\alpha}+\sum_{\alpha=1}^{s} \tilde{q}_{\alpha} d \tilde{p}_{\alpha}+(\tilde{H}-H) d t=d F_{4}(p, \tilde{p}, t) d ( F 1 − ∑ α q α p α + ∑ α p ~ α q ~ α ) = − ∑ α = 1 s q α d p α + ∑ α = 1 s q ~ α d p ~ α + ( H ~ − H ) d t = d F 4 ( p , p ~ , t )
有生成方程{ q α = − ∂ F 4 ∂ p α q ~ α = ∂ F 4 ∂ p ~ α H ~ = H + ∂ F 4 ∂ t \left\{\begin{array}{l}q_{\alpha}=-\frac{\partial F_{4}}{\partial p_{\alpha}} \\ \tilde{q}_{\alpha}=\frac{\partial F_{4}}{\partial \tilde{p}_{\alpha}} \\ \tilde{H}=H+\frac{\partial F_{4}}{\partial t}\end{array}\right. ⎩ ⎪ ⎨ ⎪ ⎧ q α = − ∂ p α ∂ F 4 q ~ α = ∂ p ~ α ∂ F 4 H ~ = H + ∂ t ∂ F 4
正则变换例子—几种常见的基本正则变换
我们通常有几种正则变换
恒等变换:q ~ α = q α , p ~ α = p α \tilde{q}_{\alpha}=q_{\alpha}, \tilde{p}_{\alpha}=p_{\alpha} q ~ α = q α , p ~ α = p α
取第二类生成函数为:F 2 ( q , p ~ , t ) = ∑ α q α p ~ α F_{2}(q, \tilde{p}, t)=\sum_{\alpha} q_{\alpha} \tilde{p}_{\alpha} F 2 ( q , p ~ , t ) = ∑ α q