第十九章 拉格朗日、哈密顿:旧力学为什么能再生新形式
一、同一只摆,为何值得换一种语言
上一章把公共物理的成立放回钟、尺、天平及其校准关系之中。现在,测量程序已经比较稳定,研究者能够反复记录一只摆的位置,问题似乎只剩把数据代入力学方程。然而,当摆连着另一只摆,支点又可以移动,或者摆球被限制在一条曲线上时,最熟悉的空间坐标未必还是最合适的组织方式。测得准确,并不保证问题已经被有效地提出。
设一颗小球由不可伸长的轻杆连接到固定支点,只能在竖直平面内运动。用横、纵坐标描述它,必须时刻保证两坐标满足同一个圆的条件,还要引入杆对小球的约束作用。若改用杆与竖直方向的夹角,允许的位置则由一个数给出。少掉的不是某块真实空间,而是一个不必独立指定的变量。物理约束被吸收到变量的选择之中。
这样的改变有时只是计算方便,有时却会改写研究者看得见的问题。角度让振动的往复结构更清楚,两个角度可以把双摆组织为一个整体状态,进一步的坐标变换又可能显出共同振动与相对振动。新语言的价值不由符号是否新颖决定,而由它是否使原先被纠缠在一起的任务获得可操作的分解与重新连接来决定。
本章因此既要说明分析力学的创造,也要限制一种过快的本体论宣布:把牛顿方程改写成拉格朗日或哈密顿形式,并不自动产生新的经验预测,更不自动证明一个新世界已经从无到有。我们需要追踪,形式改变如何参与新的对象组织、计算操作和检验安排;只有这些环节实际发生变化,形式再生才超出对既有结果的重新排版。
二、广义坐标怎样让约束进入对象
广义坐标并不神秘。它们是一组足以在研究范围内指定系统构形的独立参数,可以是距离,也可以是角度或其他适当量。对于平面单摆,一个角度即可;对于两颗可独立运动的自由质点,通常需要更多参数。坐标数量所表达的是模型保留的自由度,不是器物表面看上去包含多少零件。
必须补上“在研究范围内”这个条件。一套坐标可能在某些位置失效,或者不能用一个数值唯一识别所有构形。一个角度增加一整圈,杆的空间朝向可以相同;若模型还要记录累计转过多少圈,就需说明是在圆上取位置,还是在展开的角度轴上记录运动历史。相同符号的数值差,未必代表相同类型的物理差异。
约束的性质也很重要。轻杆长度固定,可以写成位置之间的关系,适合通过减少坐标处理;某些滚动约束涉及速度,未必能积分为单纯的位置限制。若把所有约束都当作“少写一个变量”,可能改变允许的运动集合。分析力学的整齐依赖对约束的准确分类,而不是把不方便的作用藏进一个新坐标后便不再负责。
以摆为例,沿允许位移方向,理想杆的约束作用不做虚功,这使我们能够暂时不显式求出杆的拉力而得到角运动方程。但如果研究问题恰好是杆是否断裂,张力就必须重新求出。把约束作用从某组运动方程中消去,不等于作用在现实中消失,也不等于模型已经回答关于装置的一切问题。
这种对象重组对SDE具有直接意义:变量不是任意标签,它决定哪些差异能够被独立改变,哪些差异通过约束一起变化。把一组笛卡尔坐标换成广义坐标,可能重新组织D中的可实施变化以及E中的约束关系。若只把坐标变换看作纸面换名,就容易遗漏其中被重新封装的操作条件;若把每次换名都称为存在革命,又会把真正的创造门槛降得过低。
三、从受力列表到一个函数
对本章讨论的保守系统,可把拉格朗日函数取为动能减势能,即L=K−V。用广义坐标q和广义速度q̇表达,运动方程写作:
d/dt(∂L/∂q̇) − ∂L/∂q = 0。
多自由度时,对每个独立坐标分别写一个方程。现代教材用这套形式系统处理坐标和约束,明确说明它与满足相应条件的牛顿运动方程相容;这里采用现代记号,不把教材排式冒称拉格朗日原稿的逐字复制。[19-1]
为什么一个函数能够替代许多分量关系?因为它不是一个随便挑选的曲线,而是已经把系统的动能、势能与坐标结构收进同一个表达。求导操作再从整体表达中取出各方向的变化关系。如果漏掉某一部分动能,或者把势能的方向写反,统一形式不会替我们修复错误。函数的紧凑恰恰要求对它的物理构成更认真。
对于单摆,取摆球质量m、杆长l、重力加速度g,以最低点势能为零,则:
L = ½ml²θ̇² − mgl(1−cosθ)。
对θ̇求偏导再对时间求导,得到ml²θ̈;对θ求偏导,得到−mgl sinθ。代回运动方程,有θ̈+(g/l)sinθ=0。这个结果没有先求出每时每刻的杆张力,却保留了沿允许运动方向的动力学。推导中每个系数都来自一个明确的构成条件:质量集中于摆球,杆质量忽略,长度不变,支点固定,空气阻力和摩擦暂时忽略。
若支点开始运动,或者杆可以伸长,不能只把旧方程中的某个常数随意换成时间函数便认为处理完成。动能的表达可能出现新的交叉项,自由度也可能增加。实际装置改变时,首先应回到函数的构成检查,再求运动方程。分析力学的优势,是提供一种可以重建的组织程序,而不是一条永远不用重新理解装置的公式。
四、一个能够复算的摆动检验
为了把上述推导与数字联系起来,取杆长一米,g取每二次方秒九点八米。小角度时用sinθ近似θ,得到θ̈+九点八θ=0,角频率约为每秒三点一三,周期约为二点零零七秒。设初始角为零点一弧度、初始角速度为零,则在线性近似中:
θ(t) = 0.1 cos(√9.8 t)。
从释放到第一次经过最低点的时间约为零点五零二秒。此时角速度大小约为每秒零点三一三弧度,摆球速率也约为每秒零点三一三米。用小角度势能近似计算,初始每单位质量能量为二分之一乘九点八再乘零点一的平方,即零点零四九焦耳每千克;最低点的动能给出同一结果。
还可以用未线性化的势能作一次独立检查。释放角仍为零点一弧度,由能量关系得到最低点速率√[2gl(1−cos0.1)],约为每秒零点三一二九二米。它与小角度结果接近,却不完全相同。这个微小差别不是两套力学互相矛盾,而是近似所删去的高阶影响开始显露。
如果把释放角增大,却仍要求周期严格保持二点零零七秒,问题就不是实验者不够熟练,而是线性近似被过度外推。一个良好的检验应当改变振幅,检查偏差是否随之系统增长,同时控制阻力和支点运动。若偏差只与计时方式有关,便不支持振幅解释;若在独立计时与不同装置中都跟随振幅变化,解释才获得更强约束。
这个算例同时划定创造的边界。牛顿方法、拉格朗日方法和能量方法在各自适用条件内给出相容结果,这不是三次独立实验支持。它们共享同一个模型假设。不过,三种表达能互相揭示计算错误,帮助我们辨别某个失败来自代数、近似还是物理构成。形式之间的相容属于理论自检,实验之间的相容属于经验检验,二者都重要,但证据地位不同。
五、驻作用量并不是自然预先知道终点
拉格朗日方程还可以从一个整体路径条件得到。为避免与SDE的显露S混淆,本章把作用量记为A:
A[q] = ∫L(q,q̇,t)dt。
在起止时刻及端点构形固定的变分问题中,真实运动路径满足一阶变分δA=0。这通常称为驻作用量条件。“最小作用量”是广泛使用的名称,但一般不能把驻值直接等同于所有候选路径中的绝对最小值。边界条件、允许的变分范围和函数本身的性质,都参与这一判断。[19-1]
一个常见误解由此产生:物体似乎需要先知道终点,再比较所有路线,选择最好的一条。实际上,变分形式是对运动关系的另一种数学组织,不要求自然物体实施一次有意识的规划。在适当正则条件下,它导出局部运动方程;给定初始位置与速度,又可以按局部方程求后续运动。端点变分问题与初值预测问题有联系,却不是把物体的行动解释为预知未来。
可以用最简单的自由质点直接验证。令质量为m,在时长τ内从位置零走到a。候选路径取x(t)=at/τ+b sin(πt/τ),其中第二项在起止时刻都为零。将速度平方积分,交叉项的积分为零,得到:
A = ma²/(2τ) + mb²π²/(4τ)。
在这组候选中,b=0给出较小作用量,即匀速直线路径。这里没有神秘选择,只有一个可复算的比较:加入额外的弯折速度,会增加该模型中的积分值。这个特例确有最小值,但不能从它推出所有系统的驻作用量都必定是全局最小。
如果路径中含有摩擦耗散,单纯使用上述L=K−V及封闭系统条件就未必够用,需要加入适当广义力、扩展系统或采用其他正确处理。驻作用量的威力不在于把任何过程硬塞进一个看似美丽的积分,而在于用明确条件连接局部方程与整体路径,允许我们发现哪些结构能够共存。
六、从构形走向相空间
只知道摆球在最低点,还不能知道下一刻往哪边走。它可能从左向右经过,也可能从右向左经过,甚至正停在那里。构形给出位置安排,尚未完整给出经典初值问题所需要的状态。把广义坐标与相应动量同时纳入,就得到相空间中的状态描述。
定义共轭动量p=∂L/∂q̇。若速度到动量的关系能够局部反解,便可以进行勒让德变换,写出哈密顿函数:
H(q,p,t) = pq̇ − L(q,q̇,t)。
右侧的速度最终须用q、p和t表示。多自由度时,对各pq̇项求和。哈密顿方程为q̇=∂H/∂p,ṗ=−∂H/∂q。这里的变换需要正则性;遇到退化拉格朗日函数或约束系统,不能不加处理地照用反解步骤。[19-2]
动量也不总是我们最初熟悉的“质量乘直线速度”。在单摆角坐标下,pθ=ml²θ̇,它具有角动量的量纲。坐标变化会改变共轭变量的形式,存在速度相关相互作用时更须分清正则动量与机械动量。若见到字母p就自动代入mv,反而会被新语言误导。
哈密顿的研究并非只是把微分方程换行排列。他在一八三四年《动力学的一般方法》的引言中,把复杂运动问题转向特征函数的研究,并明确承认困难可能只是由一类积分转移到另一类问题。原文还指出其方法与此前光学研究的联系。[19-3] 这种自我限制值得保留:中心函数的出现,没有使所有多体问题自动得到简单答案。
相空间的意义在于把位置与继续运动的条件放进同一个状态点,把运动组织为状态空间中的流。这个空间不是日常三维空间的另一名称。一只一维振子的相空间已经有位置和动量两个方向;许多自由度组成更高维的状态组织。维数增加不等于宇宙多出若干肉眼可见方向,而是模型需要同时指定更多独立差异。
七、同一振子在两种空间里的不同显露
取一个一维弹簧振子,质量为一千克,劲度系数为每米四牛顿。忽略阻尼,平衡点取为位置零,释放时位移为零点一米、速度为零。哈密顿函数为H=p²/2+2q²,以国际单位制数值计算。哈密顿方程给出q̇=p,ṗ=−4q,因而:
q(t) = 0.1 cos(2t),p(t) = −0.2 sin(2t)。
总能量为零点零二焦耳。物理空间里的运动是在两端点之间往复;相空间里的轨迹则是一条椭圆,满足(q/0.1)²+(p/0.2)²=1。经过四分之一周期,即π/4秒,位置为零、动量为负零点二千克米每秒。过最低势能点并不表示状态没有差异,动量方向记录了它将怎样继续。
这条椭圆把不同信息同时呈现出来。它让转折点、最大速率与能量约束互相定位,也使不同初始条件对应的轨道族容易比较。但不能因此把椭圆当成振子在房间里真正绕出的路线。视觉图上的几何是状态关系的几何,其物理意义由坐标轴和测量程序共同规定。图画得更直观,不等于图中的所有方向都是空间方向。
给模型加入阻尼后,能量会转移到未显式追踪的自由度,简单振子的相轨迹可以向平衡点收缩。若仍用无阻尼的封闭二变量哈密顿描述,却把收缩说成标准哈密顿流的性质,就混淆了模型层级。为了在更完整层面讨论能量守恒,需要把相关环境自由度纳入,或者明确承认当前使用的是有效开放系统描述。
这个比较带来一项可败条件:如果宣称装置属于无阻尼模型,就应在适当误差范围内看到能量及相轨迹的稳定;若连续记录显示系统性收缩,必须检查阻尼、测量反馈与参数漂移。不能一边保持原模型的全部前提,一边把所有不相容结果都归入“复杂性”。相空间的价值,是把这种不相容变得更容易辨认。
八、守恒结构为什么比一条轨迹更有迁移力
对于一个没有显式时间依赖的正则哈密顿系统,沿运动有dH/dt=0。直接代入两条哈密顿方程,就能看到坐标变化与动量变化的贡献相消。如果H显式依赖时间,结果则变为dH/dt=∂H/∂t,不能再因它名叫哈密顿函数就宣布它永远守恒。[19-2]
同样,如果某个广义坐标没有显式出现在拉格朗日函数中,相应共轭动量可以守恒。这个结构使我们不必每次先求出完整轨迹,再从结果中猜测有没有不变量。我们可以在求解之前检查函数怎样依赖变量,由此发现降维和比较的可能。这是形式再生的重要收益:变化从哪些地方进入,哪些量能够在变化中保持,获得了系统的提问方式。
但是“不出现”需要相对于一套明确表示来理解。一个坐标变换可以把原先显然的对称藏起来,也可以把原先不明显的对称显露出来。若变换显式依赖时间,还要处理附加项和能量解释。数学形式的简洁并不是对象本身突然变得简单,而是研究者找到了与某种结构更相容的组织方式。
在自拟振子算例中,守恒椭圆可以帮助检查数值计算。假如使用一种粗糙算法,计算轨迹逐圈向外漂移,研究者应先检查离散化是否人为增加能量,而不是立即宣布物理装置能够自行增能。相反,算法恰好保存某个量也不保证所有轨迹细节都准确;相位误差仍可能积累。守恒检查是必要的结构检验之一,不是取代全部误差分析的通行证。
由一条具体轨迹转向轨道族、守恒量和变换结构,意味着研究对象的组织层级提高了。它让不同机械系统之间出现新的可比关系,为后来跨领域的理论迁移提供了条件。但这种可能性需要后续真实研究来兑现,不能从形式优美直接跳到所有未来成果都已包含其中的历史叙事。
九、数学等价为什么不是创造的自动证明
设有两套形式,在相同对象定义、相同初始条件和同样适用范围内,能够一一映射并给出相同可观测结果。就这一范围而言,它们没有通过互相改写增加一条新的经验事实。如果只在新形式中重新算出旧轨道,再宣布物理世界已经改变,论证是不够的。
然而,经验等价也不意味着认知与研究操作毫无变化。一套形式可能使约束系统更易处理,使守恒关系先于解显露,使模型能够与另一领域发生结构性接续。这些变化改变了人—仪器—数学复合体能够稳定完成的工作。发生论应在这里提出具体证据:哪些问题以前难以组织,现在可以系统提出;哪些预测或装置设计因此能够实现;哪些失败可以更早被识别。
一个竞争解释会把这些都称为工具改进,认为对象及世界完全不变。本书不靠贬低工具来反驳它,而是追问工具改进何时进入对象规定本身。若相空间只是画出已知轨迹的辅助纸面,它的本体论分量有限;若它参与规定统计状态、可达区域或后续理论中的基本关系,研究对象就可能不再由原先的孤立位置列表充分表达。这里需要逐案论证,而非一次性宣布所有形式都是同一个层级的创造。
判断也须接受反向可能。如果新形式只带来更多符号,既不改善推导,也不增加可检验结构,所谓再生可能只是修辞。如果形式变换不满足可逆条件,却被错误宣称与旧模型完全等价,则它可能悄悄增删了物理条件。创新评价不能只看新语言的气势,要看它改变了哪些实际发生能力,并且是否如实承担了改变后的检验责任。
十、两条路径怎样把形式改写变成研究
可以从约束摆的问题具体实行D—E—S的倒逼路径。原有逐分量计算D不断被张力和几何限制缠住,研究者于是重组E,把允许构形、动能和势能组织为广义坐标中的关系;再由新结构产生可与原实验比较的角运动S。如果推导无法回到同一装置的记录,程序便没有完成。少写了张力并不算成功,正确保留它对允许运动的限制才算。
进一步实行E—S—D的锚定路径:从已建立的哈密顿结构E先推定能量轨道的某种显露S,再安排同时记录位置与速度的D,检查轨迹是闭合、漂移还是收缩。这与事后在数据上画一条椭圆不同。新的测量安排需要处理速度估计、采样误差和模型范围,形式预测才真正进入可失败的操作。
三方程在这两条路径中共同运转。S=F(D,E)要求轨道显露由实际序列与结构支撑;D=G(S,E)要求用偏差修改测量或计算;E=H(S,D)则允许稳定的不相容结果反过来改变原先的模型。不是每一次偏差都要抛弃理论,也不是每一条理论都可以靠调整数据永远保留。
三原理进一步要求,两维间的张力推动第三维改变并返回检查。S与D不相容,可推动审查E中的无阻尼假设;S与E不相容,可推动改变D的采样与初始准备;D与E不能共同支持旧状态描述,可推动S由单位置曲线扩充为相轨迹。每次扩充都须说明新增信息来自什么操作,而不是让抽象名词代替观察。
十一、旧力学的再生留下了什么边界
拉格朗日与哈密顿形式的重要性,不在于它们让牛顿力学从此没有困难,而在于它们改变了困难能够怎样被组织。约束不再只能以额外作用逐个处理,状态不再只能由空间位置图表达,守恒关系不再只能在解出轨迹以后偶然发现。新的数学组织为物理世界的继续发生提供了可以复用的条件。
与此同时,本章一直保留几条边界:驻值不等于全局最小,广义坐标不等于任意删去自由度,共轭动量不总等于质量乘速度,哈密顿函数不在所有表示和条件下都等于守恒机械能,形式等价不等于经验新增。这些边界不会削弱理论的美,反而使美能够与可靠的推导和失败检验共同成立。
还剩一个更深的问题。无论使用哪种形式,经典模型仍需规定参考系、惯性响应与允许的运动结构。我们可以优雅地描述一个系统怎样运动,但惯性为什么具有这种地位?它只是物体独自带着的性质,还是需要与其他物体以及整体世界的组织相联系?分析形式的成功没有自动消除这个问题;下一章将看到,马赫正是在经典力学最稳定的根基处重新提出了发生性的追问。
本章参考资料
[19-1] David Tong,Classical Dynamics,Chapter 2,The Lagrangian Formalism,§§2.1—2.3。https://www.damtp.cam.ac.uk/user/tong/dynamics/dynhtml/S2.html 。本章仅简述一般形式与条件;单摆、自由质点变分与数值算例为自拟推导。 [19-2] David Tong,同上,Chapter 4,The Hamiltonian Formalism,§4.1。https://www.damtp.cam.ac.uk/user/tong/dynamics/dynhtml/S4.html 。 [19-3] William Rowan Hamilton,On a General Method in Dynamics,1834,Philosophical Transactions,247—308;David R. Wilkins校订转录,2000,重点引言与开头。https://www.maths.tcd.ie/pub/HistMath/People/Hamilton/Dynamics/GenMeth.pdf 。未声称逐页精读全文。