课程进度
12 / 12
上一节狭义相对论的经典力学
自在学
分类课程智能体
文章
订阅
分类课程AI导师
文章
价格
物理理论力学核心连续系统与场的拉格朗日-哈密顿方法

连续系统与场的拉格朗日-哈密顿方法

前面各章所处理的系统,都有有限个自由度——无论是单个质点、刚体,还是由有限个粒子组成的力学系统。然而,弦的横向振动、声波在空气中的传播、电磁场在空间中的演化,这些物理现象本质上涉及无穷多个自由度。把拉格朗日力学和哈密顿力学推广到这类系统,既是经典力学体系的自然延伸,也是理解场论和现代物理的入口。从弦的振动出发,可以清晰地看到离散体系如何平滑过渡到连续场,以及对称性如何通过诺特定理与守恒律产生深层联系。


从离散系统到连续系统

以一根紧绷的弦为例。先把它离散化:将弦分成 NNN 段,每段长度 Δx\Delta xΔx,每段端点处有一个质量为 m=μΔxm = \mu\Delta xm=μΔx 的质点(μ\muμ 为线质量密度,单位 kg/m\mathrm{kg/m}kg/m),相邻质点之间以劲度系数 κ\kappaκ 的弹簧连接。第 iii 个质点的横向位移记为 yi(t)y_i(t)yi​(t),其运动方程为:

μΔx y¨i=κ ⁣(yi+1−yi)−κ ⁣(yi−yi−1)=κ ⁣(yi+1−2yi+yi−1)\mu\Delta x\,\ddot{y}_i = \kappa\!\left(y_{i+1} - y_i\right) - \kappa\!\left(y_i - y_{i-1}\right) = \kappa\!\left(y_{i+1} - 2y_i + y_{i-1}\right)μΔxy¨​i​=κ(yi+1​−yi​)−κ(yi​−yi−1​)=κ(yi+1​−2yi​+yi−1​)

令 Δx→0\Delta x \to 0Δx→0,同时保持弦的张力 T=κΔxT = \kappa\Delta xT=κΔx 不变,右边的有限差分趋近于二阶偏导数:

yi+1−2yi+yi−1(Δx)2  ⟶  ∂2y∂x2\frac{y_{i+1} - 2y_i + y_{i-1}}{(\Delta x)^2} \;\longrightarrow\; \frac{\partial^2 y}{\partial x^2}(Δx)2yi+1​−2yi​+yi−1​​⟶∂x2∂2y​

于是离散的运动方程平滑过渡为连续弦的波动方程:

μ∂2y∂t2=T∂2y∂x2\mu\frac{\partial^2 y}{\partial t^2} = T\frac{\partial^2 y}{\partial x^2}μ∂t2∂2y​=T∂x2∂2y​

波速 v=T/μv = \sqrt{T/\mu}v=T/μ​(单位:m/s\mathrm{m/s}m/s),这就是经典的一维波动方程。

与此同时,离散系统的拉格朗日量也经历同样的过渡。离散系统的拉格朗日量:

L=∑i[12μΔx y˙i2−12κ ⁣(yi+1−yi)2]L = \sum_i \left[\frac{1}{2}\mu\Delta x\,\dot{y}_i^2 - \frac{1}{2}\kappa\!\left(y_{i+1}-y_i\right)^2\right]L=i∑​[21​μΔxy˙​i2​−21​κ(yi+1​−yi​)2]

连续极限下,对 xxx 的求和变成积分,得到:

L=∫L ⁣(y,y˙,y′)dx,L=12μy˙2−12T ⁣(∂y∂x) ⁣2L = \int \mathcal{L}\!\left(y, \dot{y}, y'\right)dx, \quad \mathcal{L} = \frac{1}{2}\mu\dot{y}^2 - \frac{1}{2}T\!\left(\frac{\partial y}{\partial x}\right)^{\!2}L=∫L(y,y˙​,y′)dx,L=21​μy˙​2−21​T(∂x∂y​)2

函数 L\mathcal{L}L 称为拉格朗日密度(单位:J/m\mathrm{J/m}J/m),它是描述连续系统动力学的基本量。


拉格朗日密度与场方程

对于一个广义场 ϕ(x1,x2,x3,t)=ϕ(xμ)\phi(x^1, x^2, x^3, t) = \phi(x^\mu)ϕ(x1,x2,x3,t)=ϕ(xμ)(μ=0,1,2,3\mu = 0,1,2,3μ=0,1,2,3),其中 x0=tx^0 = tx0=t 为时间坐标,系统的拉格朗日量写为:

L=∫L ⁣(ϕ,  ∂μϕ)d3xL = \int \mathcal{L}\!\left(\phi,\;\partial_\mu\phi\right) d^3xL=∫L(ϕ,∂μ​ϕ)d3x

对应的作用量为:

S=∫L dt=∫L ⁣(ϕ,  ∂μϕ)d4xS = \int L\,dt = \int \mathcal{L}\!\left(\phi,\;\partial_\mu\phi\right) d^4xS=∫Ldt=∫L(ϕ,∂μ​ϕ)d4x

对作用量取极值 δS=0\delta S = 0δS=0,经过与离散情形完全平行的变分计算,得到欧拉-拉格朗日场方程:

∂L∂ϕ−∂μ∂L∂(∂μϕ)=0\frac{\partial \mathcal{L}}{\partial \phi} - \partial_\mu\frac{\partial \mathcal{L}}{\partial(\partial_\mu\phi)} = 0∂ϕ∂L​−∂μ​∂(∂μ​ϕ)∂L​=0

这是连续场的运动方程,地位与有限自由度系统中拉格朗日方程完全相同。

  • 例一:弦的波动方程验证

welearn-62599639.png

弦的拉格朗日密度为 L=12μy˙2−12T(y′)2\mathcal{L} = \frac{1}{2}\mu\dot{y}^2 - \frac{1}{2}T(y')^2L=21​μy˙​2−21​T(y′)2,其中 y˙=∂y/∂t\dot{y} = \partial y/\partial ty˙​=∂y/∂t,y′=∂y/∂xy' = \partial y/\partial xy′=∂y/∂x。

计算各偏导数:∂L/∂y=0\partial\mathcal{L}/\partial y = 0∂L/∂y=0,∂L/∂y˙=μy˙\partial\mathcal{L}/\partial\dot{y} = \mu\dot{y}∂L/∂y˙​=μy˙​,∂L/∂y′=−Ty′\partial\mathcal{L}/\partial y' = -Ty'∂L/∂y′=−Ty′。代入欧拉-拉格朗日场方程:

0−∂∂t(μy˙)−∂∂x(−Ty′)=0  ⟹  μy¨=Ty′′0 - \frac{\partial}{\partial t}(\mu\dot{y}) - \frac{\partial}{\partial x}(-Ty') = 0 \implies \mu\ddot{y} = Ty''0−∂t∂​(μy˙​)−∂x∂​(−Ty′)=0⟹μy¨​=Ty′′

正好就是前面得到的波动方程,验证了场方程的正确性。

  • 例二:标量场的Klein-Gordon方程

welearn-66051265.png

在相对论场论中,质量为 mmm 的标量场 ϕ(xμ)\phi(x^\mu)ϕ(xμ) 的拉格朗日密度为:

L=12ημν∂μϕ ∂νϕ−12 ⁣(mcℏ) ⁣2ϕ2\mathcal{L} = \frac{1}{2}\eta^{\mu\nu}\partial_\mu\phi\,\partial_\nu\phi - \frac{1}{2}\!\left(\frac{mc}{\hbar}\right)^{\!2}\phi^2L=21​ημν∂μ​ϕ∂ν​ϕ−21​(ℏmc​)2ϕ2

其中 ημν=diag(+1,−1,−1,−1)\eta^{\mu\nu} = \mathrm{diag}(+1,-1,-1,-1)ημν=diag(+1,−1,−1,−1) 为闵可夫斯基度规。代入场方程后得到:

□ϕ+(mcℏ) ⁣2ϕ=0,□≡∂2∂t2−c2∇2\Box\phi + \left(\frac{mc}{\hbar}\right)^{\!2}\phi = 0, \quad \Box \equiv \frac{\partial^2}{\partial t^2} - c^2\nabla^2□ϕ+(ℏmc​)2ϕ=0,□≡∂t2∂2​−c2∇2

这就是Klein-Gordon方程,描述相对论性自旋为零的粒子的传播,m=0m=0m=0 时退化为标准的波动方程。

欧拉-拉格朗日场方程 ∂L/∂ϕ−∂μ(∂L/∂(∂μϕ))=0\partial\mathcal{L}/\partial\phi - \partial_\mu(\partial\mathcal{L}/\partial(\partial_\mu\phi)) = 0∂L/∂ϕ−∂μ​(∂L/∂(∂μ​ϕ))=0 与离散系统的拉格朗日方程在形式上完全对应:只需将广义坐标 qiq_iqi​ 替换为场 ϕ\phiϕ,对时间的导数替换为对全部时空坐标的偏导数,求和替换为对空间的积分。


连续系统的哈密顿形式

类比离散系统中从拉格朗日量到哈密顿量的过渡,连续系统中先定义正则动量密度 π\piπ,其定义为:

π(x,t)=∂L∂ϕ˙\pi(x,t) = \frac{\partial \mathcal{L}}{\partial \dot{\phi}}π(x,t)=∂ϕ˙​∂L​

再通过对空间的勒让德变换得到哈密顿密度(单位:J/m3\mathrm{J/m^3}J/m3):

H=πϕ˙−L\mathcal{H} = \pi\dot{\phi} - \mathcal{L}H=πϕ˙​−L

系统总哈密顿量(即总能量)为:

H=∫H ⁣(ϕ,π,∇ϕ)d3xH = \int \mathcal{H}\!\left(\phi, \pi, \nabla\phi\right) d^3xH=∫H(ϕ,π,∇ϕ)d3x

连续系统的哈密顿方程为:

ϕ˙=δHδπ,π˙=−δHδϕ\dot{\phi} = \frac{\delta H}{\delta \pi}, \quad \dot{\pi} = -\frac{\delta H}{\delta \phi}ϕ˙​=δπδH​,π˙=−δϕδH​

其中 δH/δϕ\delta H/\delta\phiδH/δϕ 表示泛函导数(functional derivative),其定义与离散情形中的偏导数平行,只是需要同时考虑 H\mathcal{H}H 对 ϕ\phiϕ 的直接依赖和通过 ∇ϕ\nabla\phi∇ϕ 的间接依赖。

  • 例三:弦的哈密顿密度

弦的正则动量密度 π=∂L/∂y˙=μy˙\pi = \partial\mathcal{L}/\partial\dot{y} = \mu\dot{y}π=∂L/∂y˙​=μy˙​,单位为 kg/s\mathrm{kg/s}kg/s(即单位长度弦所携带的横向动量,[kg/m]×[m/s]=[kg/s][\mathrm{kg/m}]\times[\mathrm{m/s}]=[\mathrm{kg/s}][kg/m]×[m/s]=[kg/s])。哈密顿密度为:

H=πy˙−L=μy˙2−(12μy˙2−12Ty′2)=12μy˙2+12Ty′2\mathcal{H} = \pi\dot{y} - \mathcal{L} = \mu\dot{y}^2 - \left(\frac{1}{2}\mu\dot{y}^2 - \frac{1}{2}T y'^2\right) = \frac{1}{2}\mu\dot{y}^2 + \frac{1}{2}Ty'^2H=πy˙​−L=μy˙​2−(21​μy˙​2−21​Ty′2)=21​μy˙​2+21​Ty′2

用 π\piπ 表示:H=π2/(2μ)+Ty′2/2\mathcal{H} = \pi^2/(2\mu) + Ty'^2/2H=π2/(2μ)+Ty′2/2。前一项是动能密度,后一项是弹性势能密度,总哈密顿量 H=∫H dxH = \int\mathcal{H}\,dxH=∫Hdx 就是弦的总机械能,在无耗散情形下守恒。

welearn-16966597.png

应力-能量张量与守恒定理

对于时空均匀的系统,存在一个称为应力-能量张量(或能量-动量张量)的量 TμνT^{\mu\nu}Tμν,定义为:

Tμν=∂L∂(∂μϕ)∂νϕ−ημνLT^{\mu\nu} = \frac{\partial\mathcal{L}}{\partial(\partial_\mu\phi)}\partial^\nu\phi - \eta^{\mu\nu}\mathcal{L}Tμν=∂(∂μ​ϕ)∂L​∂νϕ−ημνL

它是一个 4×44\times44×4 张量,包含了系统的能量密度、动量密度以及应力(动量流密度)等全部信息。

当场方程成立时,可以验证 TμνT^{\mu\nu}Tμν 满足守恒方程:

∂μTμν=0\partial_\mu T^{\mu\nu} = 0∂μ​Tμν=0

这四个方程(ν=0,1,2,3\nu = 0, 1, 2, 3ν=0,1,2,3)分别对应能量守恒和三个方向的动量守恒。

welearn-61399999.png

  • 例四:弦的能量守恒

弦的能量密度 H=12μy˙2+12Ty′2\mathcal{H} = \frac{1}{2}\mu\dot{y}^2 + \frac{1}{2}Ty'^2H=21​μy˙​2+21​Ty′2。对时间求导:

∂H∂t=μy˙y¨+Ty′y˙′=y˙(μy¨−Ty′′)+∂∂x(Ty′y˙)\frac{\partial\mathcal{H}}{\partial t} = \mu\dot{y}\ddot{y} + Ty'\dot{y}' = \dot{y}(\mu\ddot{y} - Ty'') + \frac{\partial}{\partial x}(Ty'\dot{y})∂t∂H​=μy˙​y¨​+Ty′y˙​′=y˙​(μy¨​−Ty′′)+∂x∂​(Ty′y˙​)

利用波动方程 μy¨=Ty′′\mu\ddot{y} = Ty''μy¨​=Ty′′,括号内第一项为零,得到:

∂H∂t+∂∂x ⁣(−Ty′y˙)=0\frac{\partial\mathcal{H}}{\partial t} + \frac{\partial}{\partial x}\!\left(-Ty'\dot{y}\right) = 0∂t∂H​+∂x∂​(−Ty′y˙​)=0

这是能量守恒的局域形式:能量密度的时间变化率等于能量流 JE=−Ty′y˙J_E = -Ty'\dot{y}JE​=−Ty′y˙​(单位:W/m2\mathrm{W/m^2}W/m2)在空间上的散度的负值。

welearn-54911636.png

应力-能量张量守恒方程 ∂μTμν=0\partial_\mu T^{\mu\nu} = 0∂μ​Tμν=0 是时空均匀性的直接后果:时间平移不变性给出能量守恒,空间平移不变性给出动量守恒。这一联系正是下一节诺特定理的具体体现。


诺特定理:对称性与守恒律

1915年,德国数学家埃米·诺特证明了一条极为深刻的定理:对于每一种连续对称性,系统都对应一个守恒量。 这条定理将对称性与守恒律的关系从经验现象提升为数学必然。

设作用量在如下无穷小变换下不变:

xμ→xμ+ϵ δxμ,ϕ→ϕ+ϵ δϕx^\mu \to x^\mu + \epsilon\,\delta x^\mu, \quad \phi \to \phi + \epsilon\,\delta\phixμ→xμ+ϵδxμ,ϕ→ϕ+ϵδϕ

其中 ϵ\epsilonϵ 是无穷小参数,δxμ\delta x^\muδxμ 和 δϕ\delta\phiδϕ 描述变换的方向。诺特定理指出,必然存在一个守恒流 jμj^\mujμ,满足:

∂μjμ=0\partial_\mu j^\mu = 0∂μ​jμ=0

对应的守恒荷(charge)为:

Q=∫j0 d3x=常数Q = \int j^0\,d^3x = \text{常数}Q=∫j0d3x=常数
  • 例五:时间平移不变性给出能量守恒

若拉格朗日密度 L\mathcal{L}L 不显含时间 ttt,则作用量在时间平移 t→t+ϵt \to t + \epsilont→t+ϵ 下不变,对应 δx0=1\delta x^0 = 1δx0=1,δϕ=0\delta\phi = 0δϕ=0。诺特流的时间分量正好就是哈密顿密度:

j0=πϕ˙−L=Hj^0 = \pi\dot{\phi} - \mathcal{L} = \mathcal{H}j0=πϕ˙​−L=H

守恒条件 ∂μjμ=0\partial_\mu j^\mu = 0∂μ​jμ=0 正是前面推导的能量守恒方程,守恒荷 Q=∫H d3x=HQ = \int\mathcal{H}\,d^3x = HQ=∫Hd3x=H 就是总能量。

  • 例六:相位旋转不变性给出电荷守恒

复标量场 ϕ\phiϕ 的拉格朗日密度 L=∣∂μϕ∣2−m2∣ϕ∣2\mathcal{L} = |\partial_\mu\phi|^2 - m^2|\phi|^2L=∣∂μ​ϕ∣2−m2∣ϕ∣2 在整体相位变换 ϕ→eiαϕ\phi \to e^{i\alpha}\phiϕ→eiαϕ 下不变(α\alphaα 为实常数)。对应的诺特流为:

jμ=i ⁣(ϕ∗∂μϕ−ϕ∂μϕ∗)j^\mu = i\!\left(\phi^*\partial^\mu\phi - \phi\partial^\mu\phi^*\right)jμ=i(ϕ∗∂μϕ−ϕ∂μϕ∗)

守恒荷 Q=∫j0 d3xQ = \int j^0\,d^3xQ=∫j0d3x 在量子场论中对应粒子数之差(粒子数减去反粒子数),物理上就是电荷守恒定律的经典对应。

诺特定理的深刻之处在于:它不仅说明某个守恒量存在,还给出了守恒量的明确表达式,以及守恒的局域形式(守恒流满足连续性方程)。对称性越基本,对应的守恒律就越普遍,这一思想贯穿了整个现代理论物理。


相对论场论基础

将连续场的分析力学框架与狭义相对论结合,便进入相对论场论的领域。相对论场论的核心原则是:拉格朗日密度 L\mathcal{L}L 必须是洛伦兹标量,即在洛伦兹变换下保持不变,从而保证场方程在所有惯性参考系中具有相同的形式。

根据场的洛伦兹变换性质,场可以分类:

  • 例七:电磁场的拉格朗日密度

电磁场由四维矢量势 Aμ=(φ/c,A)A^\mu = (\varphi/c, \mathbf{A})Aμ=(φ/c,A) 描述(单位:V⋅s/m\mathrm{V\cdot s/m}V⋅s/m),场强张量定义为:

Fμν=∂μAν−∂νAμF_{\mu\nu} = \partial_\mu A_\nu - \partial_\nu A_\muFμν​=∂μ​Aν​−∂ν​Aμ​

FμνF_{\mu\nu}Fμν​ 的分量包含了电场 E\mathbf{E}E 和磁场 B\mathbf{B}B(单位:V/m\mathrm{V/m}V/m 和 T\mathrm{T}T)。电磁场的拉格朗日密度为:

L=−14μ0FμνFμν=ε02 ⁣(E2−c2B2)\mathcal{L} = -\frac{1}{4\mu_0}F_{\mu\nu}F^{\mu\nu} = \frac{\varepsilon_0}{2}\!\left(E^2 - c^2B^2\right)L=−4μ0​1​Fμν​Fμν=2ε0​​(E2−c2B2)

将这个 L\mathcal{L}L 代入欧拉-拉格朗日场方程,直接得到麦克斯韦方程组的前两组(有源方程):

∂μFμν=μ0jν\partial_\mu F^{\mu\nu} = \mu_0 j^\nu∂μ​Fμν=μ0​jν

而后两组(无源方程,即 ∇⋅B=0\nabla\cdot\mathbf{B}=0∇⋅B=0 和法拉第定律)则是 FμνF_{\mu\nu}Fμν​ 反对称性的自动结果,无需额外引入。

相对论场论的内容仅作概念性介绍,完整的量子场论涉及场的量子化(正则量子化或路径积分量子化),超出经典力学的范围。这里的目的是展示拉格朗日-哈密顿方法的普适性:无论是一根振动的弦,还是充斥整个空间的电磁场,都可以纳入同一套变分原理的框架。


练习题

  • 1. 关于从离散系统过渡到连续系统,下列说法正确的是:

A. 连续系统的拉格朗日量等于拉格朗日密度 L\mathcal{L}L 本身,不需要对空间积分

B. 对 NNN 个质点组成的弦离散模型取极限 Δx→0\Delta x \to 0Δx→0 时,差分方程变为偏微分波动方程

C. 拉格朗日密度 L\mathcal{L}L 的单位与拉格朗日量 LLL 的单位相同,均为 J\mathrm{J}J

D. 连续弦的波速 v=T/μv = \sqrt{T/\mu}v=T/μ​ 中,TTT 是绝对温度

答案:B

分析:连续系统的拉格朗日量是拉格朗日密度对空间的积分 L=∫L dxL = \int\mathcal{L}\,dxL=∫Ldx,A 错误。拉格朗日密度 L\mathcal{L}L 的单位是 J/m\mathrm{J/m}J/m(一维)或 J/m3\mathrm{J/m^3}J/m3(三维),而 LLL 的单位是 J\mathrm{J}J,C 错误。弦波速公式中 TTT 是弦的张力(单位 N\mathrm{N}N),不是温度,D 错误。B 正确:对离散弦模型取连续极限,有限差分 (yi+1−2yi+yi−1)/(Δx)2(y_{i+1}-2y_i+y_{i-1})/(\Delta x)^2(yi+1​−2yi​+yi−1​)/(Δx)2 趋近于 ∂2y/∂x2\partial^2 y/\partial x^2∂2y/∂x2,运动方程变为波动方程 μ∂2y/∂t2=T∂2y/∂x2\mu\partial^2 y/\partial t^2 = T\partial^2 y/\partial x^2μ∂2y/∂t2=T∂2y/∂x2。

  • 2. 将欧拉-拉格朗日场方程应用于拉格朗日密度 L=12(∂μϕ)(∂μϕ)\mathcal{L} = \frac{1}{2}(\partial_\mu\phi)(\partial^\mu\phi)L=21​(∂μ​ϕ)(∂μϕ) 时,得到的场方程是:

A. ϕ=0\phi = 0ϕ=0(ϕ\phiϕ 等于零)

B. ∇2ϕ=0\nabla^2\phi = 0∇2ϕ=0(泊松方程,仅含空间导数)

C. □ϕ=0\Box\phi = 0□ϕ=0,即 ∂2ϕ/∂t2=c2∇2ϕ\partial^2\phi/\partial t^2 = c^2\nabla^2\phi∂2ϕ/∂t2=c2∇2ϕ(波动方程)

D. ∂ϕ/∂t=0\partial\phi/\partial t = 0∂ϕ/∂t=0(ϕ\phiϕ 不随时间变化)

答案:C

分析:对 L=12(∂μϕ)(∂μϕ)\mathcal{L} = \frac{1}{2}(\partial_\mu\phi)(\partial^\mu\phi)L=21​(∂μ​ϕ)(∂μϕ),有 ∂L/∂ϕ=0\partial\mathcal{L}/\partial\phi = 0∂L/∂ϕ=0,∂L/∂(∂μϕ)=∂μϕ\partial\mathcal{L}/\partial(\partial_\mu\phi) = \partial^\mu\phi∂L/∂(∂μ​ϕ)=∂μϕ。代入场方程 ∂μ(∂μϕ)=0\partial_\mu(\partial^\mu\phi) = 0∂μ​(∂μϕ)=0,即 □ϕ=0\Box\phi = 0□ϕ=0,就是无质量标量场满足的波动方程,C 正确。这里 □=∂2/∂t2−c2∇2\Box = \partial^2/\partial t^2 - c^2\nabla^2□=∂2/∂t2−c2∇2 是达朗贝尔算符,包含时间和空间的二阶导数,B(只有空间导数)和 D(只有时间导数)均不完整。

  • 3. 关于诺特定理,下列说法正确的是:

A. 诺特定理仅适用于有限自由度系统,不适用于场论

B. 空间平移不变性对应能量守恒,时间平移不变性对应动量守恒

C. 每一种连续对称性对应一个守恒流 jμj^\mujμ,满足 ∂μjμ=0\partial_\mu j^\mu = 0∂μ​jμ=0,对应守恒荷 Q=∫j0d3xQ = \int j^0 d^3xQ=∫j0d3x

D. 规范对称性(相位变换不变性)对应角动量守恒

答案:C

分析:诺特定理对有限自由度系统和场论都成立,A 错误。能量守恒来自时间平移不变性,动量守恒来自空间平移不变性,B 的说法正好相反,B 错误。规范对称性(如复场相位变换)对应电荷守恒,不是角动量守恒,D 错误。C 正确地描述了诺特定理的内容:连续对称性 ⇒\Rightarrow⇒ 守恒流 ∂μjμ=0\partial_\mu j^\mu = 0∂μ​jμ=0 ⇒\Rightarrow⇒ 守恒荷 QQQ 为常数。

  • 4. 弦的正则动量密度 π=∂L/∂y˙=μy˙\pi = \partial\mathcal{L}/\partial\dot{y} = \mu\dot{y}π=∂L/∂y˙​=μy˙​,哈密顿密度 H=π2/(2μ)+T(y′)2/2\mathcal{H} = \pi^2/(2\mu) + T(y')^2/2H=π2/(2μ)+T(y′)2/2。下列关于弦的哈密顿量的说法,正确的是:

A. 总哈密顿量 H=∫H dxH = \int\mathcal{H}\,dxH=∫Hdx 表示弦的总势能,不包含动能

B. H\mathcal{H}H 中的 π2/(2μ)\pi^2/(2\mu)π2/(2μ) 项对应弦的横向动能密度,T(y′)2/2T(y')^2/2T(y′)2/2 对应弹性势能密度

C. 弦的哈密顿量 HHH 不守恒,因为弦是连续系统

D. 正则动量密度 π\piπ 的单位与广义动量(kg⋅m/s\mathrm{kg\cdot m/s}kg⋅m/s)相同

答案:B

分析:哈密顿密度 H=π2/(2μ)+T(y′)2/2\mathcal{H} = \pi^2/(2\mu) + T(y')^2/2H=π2/(2μ)+T(y′)2/2 中,第一项 =μy˙2/2= \mu\dot{y}^2/2=μy˙​2/2 是单位长度的动能(动能密度),第二项是单位长度的弹性势能(势能密度),总哈密顿量是两者之和,A 错误。在无外力和无耗散情形下,弦的总机械能守恒,C 错误。正则动量密度 π=μy˙\pi = \mu\dot{y}π=μy˙​ 的单位是 kg/(m⋅s)\mathrm{kg/(m\cdot s)}kg/(m⋅s)(即动量除以长度),而不是 kg⋅m/s\mathrm{kg\cdot m/s}kg⋅m/s,D 错误。B 正确。

  • 5. 一根两端固定、长度为 LLL(单位 m\mathrm{m}m)的弦,线密度 μ\muμ(kg/m\mathrm{kg/m}kg/m),张力 TTT(N\mathrm{N}N),自然振动的驻波解形如 yn(x,t)=Ansin⁡(nπx/L)cos⁡(ωnt+φn)y_n(x,t) = A_n\sin(n\pi x/L)\cos(\omega_n t + \varphi_n)yn​(x,t)=An​sin(nπx/L)cos(ωn​t+φn​),其中 n=1,2,3,…n = 1,2,3,\ldotsn=1,2,3,…。(一)将驻波解代入波动方程 μy¨=Ty′′\mu\ddot{y} = Ty''μy¨​=Ty′′,求各阶角频率 ωn\omega_nωn​;(二)写出第 nnn 阶模式的拉格朗日密度 Ln\mathcal{L}_nLn​,并计算整弦的拉格朗日量 Ln=∫0LLn dxL_n = \int_0^L \mathcal{L}_n\,dxLn​=∫0L​Ln​dx;(三)说明为何不同阶模式之间的拉格朗日量相互独立(提示:利用正弦函数的正交性)。

解:

(一)将 yn=Ansin⁡(nπx/L)cos⁡(ωnt+φn)y_n = A_n\sin(n\pi x/L)\cos(\omega_n t+\varphi_n)yn​=An​sin(nπx/L)cos(ωn​t+φn​) 代入 μy¨=Ty′′\mu\ddot{y} = Ty''μy¨​=Ty′′:

−μωn2Ansin⁡nπxLcos⁡(ωnt+φn)=−T ⁣(nπL) ⁣2Ansin⁡nπxLcos⁡(ωnt+φn)-\mu\omega_n^2 A_n\sin\frac{n\pi x}{L}\cos(\omega_n t+\varphi_n) = -T\!\left(\frac{n\pi}{L}\right)^{\!2}A_n\sin\frac{n\pi x}{L}\cos(\omega_n t+\varphi_n)−μωn2​An​sinLnπx​cos(ωn​t+φn​)=−T(Lnπ​)2An​sinLnπx​cos(ωn​t+φn​)

两边约去公因子,得:

ωn2=Tμ ⁣(nπL) ⁣2  ⟹  ωn=nπLTμ=nπvL\omega_n^2 = \frac{T}{\mu}\!\left(\frac{n\pi}{L}\right)^{\!2} \implies \omega_n = \frac{n\pi}{L}\sqrt{\frac{T}{\mu}} = \frac{n\pi v}{L}ωn2​=μT​(Lnπ​)2⟹ωn​=Lnπ​μT​​=Lnπv​

其中 v=T/μv = \sqrt{T/\mu}v=T/μ​(m/s\mathrm{m/s}m/s)为波速,ωn\omega_nωn​ 的单位为 rad/s\mathrm{rad/s}rad/s。基频(n=1n=1n=1)为 ω1=πv/L\omega_1 = \pi v/Lω1​=πv/L,各阶频率均为基频的整数倍。

(二)第 nnn 阶模式的横向速度和斜率分别为:

y˙n=−Anωnsin⁡nπxLsin⁡(ωnt+φn),yn′=AnnπLcos⁡nπxLcos⁡(ωnt+φn)\dot{y}_n = -A_n\omega_n\sin\frac{n\pi x}{L}\sin(\omega_n t+\varphi_n), \quad y_n' = A_n\frac{n\pi}{L}\cos\frac{n\pi x}{L}\cos(\omega_n t+\varphi_n)y˙​n​=−An​ωn​sinLnπx​sin(ωn​t+φn​),yn′​=An​Lnπ​cosLnπx​cos(ωn​t+φn​)

拉格朗日密度:

Ln=12μy˙n2−12T(yn′)2\mathcal{L}_n = \frac{1}{2}\mu\dot{y}_n^2 - \frac{1}{2}T(y_n')^2Ln​=21​μy˙​n2​−21​T(yn′​)2=12μAn2ωn2sin⁡2 ⁣nπxLsin⁡2(ωnt+φn)−12TAn2 ⁣(nπL) ⁣2cos⁡2 ⁣nπxLcos⁡2(ωnt+φn)= \frac{1}{2}\mu A_n^2\omega_n^2\sin^2\!\frac{n\pi x}{L}\sin^2(\omega_n t+\varphi_n) - \frac{1}{2}T A_n^2\!\left(\frac{n\pi}{L}\right)^{\!2}\cos^2\!\frac{n\pi x}{L}\cos^2(\omega_n t+\varphi_n)=21​μAn2​ωn2​sin2Lnπx​sin2(ωn​t+φn​)−21​TAn2​(Lnπ​)2cos2Lnπx​cos2(ωn​t+φn​)

对 xxx 从 000 到 LLL 积分,利用 ∫0Lsin⁡2(nπx/L) dx=∫0Lcos⁡2(nπx/L) dx=L/2\int_0^L\sin^2(n\pi x/L)\,dx = \int_0^L\cos^2(n\pi x/L)\,dx = L/2∫0L​sin2(nπx/L)dx=∫0L​cos2(nπx/L)dx=L/2:

Ln=14μAn2ωn2Lsin⁡2(ωnt+φn)−14TAn2 ⁣(nπL) ⁣2Lcos⁡2(ωnt+φn)L_n = \frac{1}{4}\mu A_n^2\omega_n^2 L\sin^2(\omega_n t+\varphi_n) - \frac{1}{4}T A_n^2\!\left(\frac{n\pi}{L}\right)^{\!2} L\cos^2(\omega_n t+\varphi_n)Ln​=41​μAn2​ωn2​Lsin2(ωn​t+φn​)−41​TAn2​(Lnπ​)2Lcos2(ωn​t+φn​)

由 ωn2=T(nπ/L)2/μ\omega_n^2 = T(n\pi/L)^2/\muωn2​=T(nπ/L)2/μ,两系数的大小相等,故:

Ln=μAn2ωn2L4 ⁣[sin⁡2(ωnt+φn)−cos⁡2(ωnt+φn)]=−μAn2ωn2L4cos⁡(2ωnt+2φn)L_n = \frac{\mu A_n^2\omega_n^2 L}{4}\!\left[\sin^2(\omega_n t+\varphi_n) - \cos^2(\omega_n t+\varphi_n)\right] = -\frac{\mu A_n^2\omega_n^2 L}{4}\cos(2\omega_n t+2\varphi_n)Ln​=4μAn2​ωn2​L​[sin2(ωn​t+φn​)−cos2(ωn​t+φn​)]=−4μAn2​ωn2​L​cos(2ωn​t+2φn​)

(三)对于两个不同阶数 m≠nm \neq nm=n,混合项中含有 ∫0Lsin⁡(mπx/L)sin⁡(nπx/L) dx=0\int_0^L\sin(m\pi x/L)\sin(n\pi x/L)\,dx = 0∫0L​sin(mπx/L)sin(nπx/L)dx=0 和 ∫0Lcos⁡(mπx/L)cos⁡(nπx/L) dx=0\int_0^L\cos(m\pi x/L)\cos(n\pi x/L)\,dx = 0∫0L​cos(mπx/L)cos(nπx/L)dx=0(正弦和余弦函数族在 [0,L][0,L][0,L] 上满足正交性)。因此,总拉格朗日量 L=∫0LL dxL = \int_0^L\mathcal{L}\,dxL=∫0L​Ldx 在展开时,不同模式之间的交叉项全部消失,各阶模式的拉格朗日量直接相加,彼此独立——连续弦等效于无穷多个相互独立的谐振子。

  • 6. 对于拉格朗日密度 L=12μy˙2−12T(y′)2\mathcal{L} = \frac{1}{2}\mu\dot{y}^2 - \frac{1}{2}T(y')^2L=21​μy˙​2−21​T(y′)2,用诺特定理验证空间平移不变性对应动量守恒。(一)写出对应空间平移 x→x+ϵx \to x + \epsilonx→x+ϵ 的诺特流;(二)写出对应的守恒方程;(三)计算弦的总广义动量 P=∫j0 dxP = \int j^0\,dxP=∫j0dx 并说明其物理意义。

解:

(一)空间平移 x→x+ϵx \to x + \epsilonx→x+ϵ 对应 δx=1\delta x = 1δx=1,δt=0\delta t = 0δt=0,场的变化 δy=−y′\delta y = -y'δy=−y′(场随坐标平移的变化)。

对于空间平移,诺特流(动量流密度)的时间分量为:

j0=∂L∂y˙δy=π⋅(−y′)=−μy˙⋅y′j^0 = \frac{\partial\mathcal{L}}{\partial\dot{y}}\delta y = \pi\cdot(-y') = -\mu\dot{y}\cdot y'j0=∂y˙​∂L​δy=π⋅(−y′)=−μy˙​⋅y′

空间分量(即动量流密度,即应力)为:

j1=∂L∂y′δy−L⋅δx=(−Ty′)(−y′)−L=T(y′)2−12μy˙2+12T(y′)2=32T(y′)2−12μy˙2j^1 = \frac{\partial\mathcal{L}}{\partial y'}\delta y - \mathcal{L}\cdot\delta x = (-Ty')(-y') - \mathcal{L} = T(y')^2 - \frac{1}{2}\mu\dot{y}^2 + \frac{1}{2}T(y')^2 = \frac{3}{2}T(y')^2 - \frac{1}{2}\mu\dot{y}^2j1=∂y′∂L​δy−L⋅δx=(−Ty′)(−y′)−L=T(y′)2−21​μy˙​2+21​T(y′)2=23​T(y′)2−21​μy˙​2

更标准的推导给出应力-能量张量分量 T10=−Ty′y˙T^{10} = -Ty'\dot{y}T10=−Ty′y˙​(动量流)和 T00=HT^{00} = \mathcal{H}T00=H(能量密度),这里动量守恒的诺特流为:

j0=−μy˙ y′,j1=T11=12μy˙2−12T(y′)2=−Lj^0 = -\mu\dot{y}\,y', \quad j^1 = T^{11} = \frac{1}{2}\mu\dot{y}^2 - \frac{1}{2}T(y')^2 = -\mathcal{L}j0=−μy˙​y′,j1=T11=21​μy˙​2−21​T(y′)2=−L

(二)守恒方程为 ∂tj0+∂xj1=0\partial_t j^0 + \partial_x j^1 = 0∂t​j0+∂x​j1=0,即:

∂∂t(−μy˙ y′)+∂∂x ⁣(−L)=0\frac{\partial}{\partial t}(-\mu\dot{y}\,y') + \frac{\partial}{\partial x}\!\left(-\mathcal{L}\right) = 0∂t∂​(−μy˙​y′)+∂x∂​(−L)=0

验证:展开左边第一项:

∂∂t(−μy˙ y′)=−μy¨ y′−μy˙ y˙′\frac{\partial}{\partial t}(-\mu\dot{y}\,y') = -\mu\ddot{y}\,y' - \mu\dot{y}\,\dot{y}' ∂t∂​(−μy˙​y′)=−μy¨​y′−μy˙​y˙​′

第二项:

∂∂x(−L)=−∂L∂x∣显含−∂L∂yy′−∂L∂y˙y˙′−∂L∂y′y′′=μy˙y˙′+Ty′y′′\frac{\partial}{\partial x}(-\mathcal{L}) = -\frac{\partial\mathcal{L}}{\partial x}\bigg|_{\text{显含}} - \frac{\partial\mathcal{L}}{\partial y}y' - \frac{\partial\mathcal{L}}{\partial\dot{y}}\dot{y}' - \frac{\partial\mathcal{L}}{\partial y'}y'' = \mu\dot{y}\dot{y}' + Ty'y''∂x∂​(−L)=−∂x∂L​​显含​−∂y∂L​y′−∂y˙​∂L​y˙​′−∂y′∂L​y′′=μy˙​y˙​′+Ty′y′′

两项之和为 −μy¨y′+Ty′y′′=y′(Ty′′−μy¨)=0-\mu\ddot{y}y' + Ty'y'' = y'(T y'' - \mu\ddot{y}) = 0−μy¨​y′+Ty′y′′=y′(Ty′′−μy¨​)=0(利用波动方程),守恒方程成立。

(三)总广义动量:

P=∫0Lj0 dx=−∫0Lμy˙ y′ dx=−μ∫0L∂y∂t∂y∂x dxP = \int_0^L j^0\,dx = -\int_0^L \mu\dot{y}\,y'\,dx = -\mu\int_0^L\frac{\partial y}{\partial t}\frac{\partial y}{\partial x}\,dxP=∫0L​j0dx=−∫0L​μy˙​y′dx=−μ∫0L​∂t∂y​∂x∂y​dx

PPP 的单位为 kg⋅m/s\mathrm{kg\cdot m/s}kg⋅m/s,称为弦的“准动量”(pseudo-momentum)或“场动量”,它描述弦上振动波型在空间中整体平移所对应的守恒量。对于两端固定的弦(边界条件 y(0,t)=y(L,t)=0y(0,t) = y(L,t) = 0y(0,t)=y(L,t)=0),驻波情形下左行波和右行波的动量大小相等、方向相反,总场动量为零,与空间平移对称性在有限边界下受到破坏的情况相吻合。

  • 从离散系统到连续系统
  • 拉格朗日密度与场方程
  • 连续系统的哈密顿形式
  • 应力-能量张量与守恒定理
  • 诺特定理:对称性与守恒律
  • 相对论场论基础
  • 练习题

目录

  • 从离散系统到连续系统
  • 拉格朗日密度与场方程
  • 连续系统的哈密顿形式
  • 应力-能量张量与守恒定理
  • 诺特定理:对称性与守恒律
  • 相对论场论基础
  • 练习题
自在学

© 2025 - 2026 自在学,保留所有权利。

公网安备湘公网安备43020302000292号 | 湘ICP备2025148919号-1

关于我们隐私政策使用条款

© 2025 自在学,保留所有权利。

公网安备湘公网安备43020302000292号湘ICP备2025148919号-1