1062 字
5 分钟
哈密顿方程

一、勒让德变换#

NOTE

将原有的独立变量的一部分或者全部,换为相应的共轭变量,从而得到一个新函数,这种变换称为勒让德变换。

1、一元函数勒让德变换#

函数f(x)f(x)是变量xx的函数,且u(x)=ddxf(x)u(x)=\frac{\mathrm{d}}{\mathrm{d}x}f(x),可以构造一个新的以uu为自变量的函数:

g=xufg=xu-f

gg的全微分:

dg=xdu+udxdf=xdu\mathrm{d}g=x\mathrm{d}u+u\mathrm{d}x-\mathrm{d}f=x\mathrm{d}u

也就是:

x=dgdux=\frac{\mathrm{d}g}{\mathrm{d}u}

于是,gg就是变量uu的函数:

dg(u)=dgdudu=x(u)du\mathrm{d}g(u)=\frac{\mathrm{d}g}{\mathrm{d}u}\mathrm{d}u=x(u)\mathrm{d}u

我们称g=xufg=xu-f称为函数ff的勒让德变换,它存在的条件,就是xx可以表示为uu的函数,也就是说uu存在反函数,根据反函数定理,要求:

dudx0\frac{\mathrm{d}u}{\mathrm{d}x}\neq0

而在真实的物理系统中,总是存在勒让德变换。

2、二元函数勒让德变换#

(1) 变换一#

函数f(x,y)f(x,y)是变量xxyy的函数,且u=dfdxv=dfdyu=\frac{\mathrm{d}f}{\mathrm{d}x},v=\frac{\mathrm{d}f}{\mathrm{d}y},故有:

df=udx+vdy\mathrm{d}f=u\mathrm{d}x+v\mathrm{d}y

为了得到一个以u,vu,v为自变量的新函数,可以构造:

g=xu+vyfg=xu+vy-f

gg的全微分:

dg=xdu+udx+vdy+ydvdf=xdu+vdy\mathrm{d}g=x\mathrm{d}u+u\mathrm{d}x+v\mathrm{d}y+y\mathrm{d}v-\mathrm{d}f=x\mathrm{d}u+v\mathrm{d}y

根据全微分定义,得到:

x=gu,y=gyx=\frac{\partial g}{\partial u},y=\frac{\partial g}{\partial y}

于是,gg就是变量uuyy的函数。我们称g=xu+vyfg=xu+vy-f称为函数ff的勒让德变换,它存在的条件,就是x,yx,y可以表示为u,vu,v的函数,也就是说u,vu,v存在反函数,根据反函数定理,要求:

uxuyvxvy0\left| \begin{matrix} \frac{\partial u}{\partial x} & \frac{\partial u}{\partial y} \\ \frac{\partial v}{\partial x} & \frac{\partial v}{\partial y} \end{matrix} \right| \neq 0

(2) 变换二#

函数f(x,y)f(x,y)是变量xxyy的函数,且u=dfdxv=dfdyu=\frac{\mathrm{d}f}{\mathrm{d}x},v=\frac{\mathrm{d}f}{\mathrm{d}y},为了得到一个以u,yu,y为自变量的新函数,可以构造:

g=xufg=xu-f

gg的全微分:

dg=xdu+udxdf=xduvdy\mathrm{d}g=x\mathrm{d}u+u\mathrm{d}x-\mathrm{d}f=x\mathrm{d}u-v\mathrm{d}y

根据全微分定义,得到:

x=gu,y=gyx=\frac{\partial g}{\partial u},y=-\frac{\partial g}{\partial y}

于是,gg就是变量uuyy的函数。我们称g=xufg=xu-f称为函数ff的勒让德变换,它存在的条件,就是xx可以表示为uu的函数,也就是说uu存在反函数,根据反函数定理,要求:

dudx0\frac{\mathrm{d}u}{\mathrm{d}x}\neq0

3、勒让德变换的一般方法#

函数f(x1,x2,,xn;a1,a2,,am)f(x_1,x_2,\cdots,x_n;a_1,a_2,\cdots,a_m)是变量x1,x2,,xn;a1,a2,,amx_1,x_2,\cdots,x_n;a_1,a_2,\cdots,a_m的函数,有:

yi=fxi(i=1,2,,n)y_i=\frac{\partial f}{\partial x_i}\quad(i=1,2,\cdots,n)

将变量x1,x2,,xn;a1,a2,,amx_1,x_2,\cdots,x_n;a_1,a_2,\cdots,a_m变换为y1,y2,,yn;a1,a2,,amy_1,y_2,\cdots,y_n;a_1,a_2,\cdots,a_m,相应的函数ff变为g(y1,y2,,yn;a1,a2,,am)g(y_1,y_2,\cdots,y_n;a_1,a_2,\cdots,a_m),可以构造:

g=i=1nxiyifg=\sum_{i=1}^n x_iy_i-f

勒让德变换存在的条件为x,yx,y的雅可比行列式不为00,在真实的物理系统中都存在勒让德变换。

勒让德变换的逆变换:

f=i=1nxiyigf=\sum_{i=1}^n x_iy_i-g

同样也适合勒让德逆变换

二、哈密顿方程#

1、保守系统的哈密顿方程#

设质点系由nn个质点组成,受到ss个完整约束,具有k=3nsk=3n-s个自由度。理想约束情况下,保守系统的拉格朗日方程为:

ddt(Lq˙i)=Lqi\frac{\mathrm{d}}{\mathrm{d}t}\left(\frac{\partial L}{\partial \dot q_i}\right)=\frac{\partial L}{\partial q_i}

引入广义动量:

pi=Lq˙i(i=1,2,,k)L=L(q˙1,q˙2,,q˙k;t)p_i=\frac{\partial L}{\partial \dot q_i}\quad(i=1,2,\cdots,k)\quad L=L(\dot{q}_1,\dot q_2,\cdots,\dot q_k;t)

以广义动量p1,p2,,pkp_1,p_2,\cdots,p_k和广义坐标q1,q2,,qkq_1,q_2,\cdots,q_k作为描述系统状态变量,建立系统的运动微分方程组:

H=H(p1,p2,,pk;q1,q2,,qk;t)=(i=1kpiq˙iL)q˙ipiH=H(p_1,p_2,\cdots,p_k;q_1,q_2,\cdots,q_k;t)=\left(\sum_{i=1}^k p_i\dot q_i-L\right)_{\dot q_i\to p_i}

其中的函数HH称为哈密顿函数。对HH取全微分,得到:

q˙i=Hpi\dot q_i=\frac{\partial H}{\partial p_i}Hqi=Lqi\frac{\partial H}{\partial q_i}=-\frac{\partial L}{\partial q_i}

利用拉格朗日方程:

p˙i=ddt(Lq˙i)=Lqi=Hqi\dot p_i=\frac{\mathrm{d}}{\mathrm{d}t}\left(\frac{\partial L}{\partial \dot q_i}\right)=\frac{\partial L}{\partial q_i}=-\frac{\partial H}{\partial q_i}

综合即为哈密顿方程,是关于系统状态变量q1,q2,,qkq_1,q_2,\cdot,q_kp1,p2,,pkp_1,p_2,\cdots,p_k2k2k个一阶微分方程组:

q˙i=Hpi\dot q_i=\frac{\partial H}{\partial p_i}p˙i=Hqi\dot p_i=-\frac{\partial H}{\partial q_i}

哈密顿方程的数量是系统自由度的2倍。

2、非保守系统的哈密顿方程#

对于非保守系统,主动力可以分为有势力和非有势力两类,系统的拉格朗日方程为:

ddt(Lq˙i)Lqi=Q~\frac{\mathrm{d}}{\mathrm{d}t}\left(\frac{\partial L}{\partial \dot q_i}\right)-\frac{\partial L}{\partial q_i}=\tilde{Q}

那么广义动量关于时间的全导数为:

p˙i=ddt(Lq˙i)=Lqi+Q~=Hqi+Q~\dot p_i=\frac{\mathrm{d}}{\mathrm{d}t}\left(\frac{\partial L}{\partial \dot q_i}\right)=\frac{\partial L}{\partial q_i}+\tilde{Q}=-\frac{\partial H}{\partial q_i}+\tilde{Q}

于是非保守系统的哈密顿方程为:

q˙i=Hpi\dot q_i=\frac{\partial H}{\partial p_i}p˙i=Hqi+Q~\dot p_i=-\frac{\partial H}{\partial q_i}+\tilde{Q}

3、哈密顿函数#

广义速度表示的动能为:

T=12i=1kj=1kaijq˙iq˙j+i=1kbiq˙i+c=T2+T1+T0\begin{aligned} T &= \frac{1}{2} \sum_{i=1}^{k} \sum_{j=1}^{k} a_{ij} \dot{q}_i \dot{q}_j + \sum_{i=1}^{k} b_i \dot{q}_i + c \\ &= T_2 + T_1 + T_0 \end{aligned}

哈密顿函数HH是通过广义动量和广义坐标表示的广义能量:

H=(i=1kpiq˙iL)q˙ipi=(i=1kLq˙iq˙iL)q˙ipi=(i=1k(j=1kaijq˙j+bi)q˙iL)q˙ipi=[2T2+T1(T2+T1+T0V)]q˙ipi=(T2+VT0)q˙ipi\begin{aligned} H &= \left( \sum_{i=1}^{k} p_i \dot{q}_i - L \right)_{\dot{q}_i \to p_i} = \left( \sum_{i=1}^{k} \frac{\partial L}{\partial \dot{q}_i} \dot{q}_i - L \right)_{\dot{q}_i \to p_i} \\ &= \left( \sum_{i=1}^{k} \left( \sum_{j=1}^{k} a_{ij} \dot{q}_j + b_i \right) \dot{q}_i - L \right)_{\dot{q}_i \to p_i} = \left[ 2T_2 + T_1 - (T_2 + T_1 + T_0 - V) \right]_{\dot{q}_i \to p_i} \\ &= (T_2 + V - T_0)_{\dot{q}_i \to p_i} \end{aligned}

那么对于定常约束系统,满足T1=T0=0,T=T2T_1=T_0=0,T=T_2,有:

H=(T+V)q˙ipiH=(T+V)_{\dot q_i\to p_i}

即通过广义动量和广义坐标表示的机械能。

分享

如果这篇文章对你有帮助,欢迎分享给更多人!

哈密顿方程
https://www.laoguantx.cn/posts/hamiltonsequations/
作者
老官童鞋gogo
发布于
2025-12-25
许可协议
CC BY-NC-SA 4.0

部分信息可能已经过时

封面
Sample Song
Sample Artist
封面
Sample Song
Sample Artist
0:00 / 0:00