物理学的数学工具

物理学每一次重大突破，几乎都伴随着数学工具的革新。牛顿发明微积分来描述运动，麦克斯韦用偏微分方程统一电磁场，爱因斯坦借助黎曼几何构建广义相对论。数学不是物理的附属品，而是物理思想最精确的表达方式。掌握几类核心数学工具——偏微分方程、格林函数、变分法和对称性——能让人以更高的视角理解物理规律的共性结构。

这些工具看起来抽象，但每一个都有清晰的物理根源，并能解决实际问题。偏微分方程刻画波与热的传播，格林函数把复杂问题化为点源叠加，变分法揭示自然界“走最优路径”的内在原则，诺特定理则把守恒律的根源归结为对称性。

偏微分方程：描述连续系统的语言

普通微分方程（ODE）只有一个自变量，例如只依赖时间 $t$；偏微分方程（PDE）则同时处理多个自变量，例如既依赖空间坐标 $x$，又依赖时间 $t$ 的物理量 $u(x, t)$。物理中最重要的三类线性偏微分方程，各自对应一类本质不同的物理现象。

波动方程

弦振动是波动方程最直观的来源。一根均匀弦，线密度为 $\mu$（单位：$\text{kg/m}$），张力为 $T_0$（单位：$\text{N}$），在很小的横向位移 $u(x,t)$ 下，对一小段 $dx$ 做受力分析，可推导出一维波动方程：

$\frac{\partial^2 u}{\partial t^2} = v^2 \frac{\partial^2 u}{\partial x^2}, \quad v = \sqrt{\frac{T_0}{\mu}}$

其中 $v$ 是波的传播速度（单位：$\text{m/s}$）。方程左边是弦上质点的加速度（对时间的二阶导数），右边是弦的弯曲程度（对空间的二阶导数）。弦越弯，受到的净恢复力越大，加速度也越大——这正是波动的物理机制。

对于两端固定（$u(0,t)=0$，$u(L,t)=0$）、长度为 $L$ 的弦，固有的振动模式称为驻波，形式为：

$u_n(x,t) = A_n \sin\!\left(\frac{n\pi x}{L}\right)\cos(\omega_n t + \varphi_n), \quad \omega_n = \frac{n\pi v}{L}, \quad n = 1, 2, 3, \ldots$

各阶频率 $f_n = \omega_n/(2\pi) = nv/(2L)$，恰好是基频的整数倍。

热传导方程

将一根细棒局部加热后，热量会从高温处向低温处扩散。设棒的温度分布为 $T(x,t)$，热扩散率为 $\alpha$（单位：$\text{m}^2/\text{s}$），由傅里叶热传导定律推导出热传导方程（也叫扩散方程）：

$\frac{\partial T}{\partial t} = \alpha \frac{\partial^2 T}{\partial x^2}$

含义直观：一点温度的变化速率（左边）正比于该点处的温度曲率（右边）——曲率越大，说明该点与周围温差的变化越剧烈，热量流动越快。

热传导方程与波动方程的关键区别在于时间导数的阶数：波动方程是二阶（$\partial^2 u/\partial t^2$），热传导方程是一阶（$\partial T/\partial t$）。把 $t$ 换成 $-t$，热传导方程变为 ，与原方程不同——这表明热传导是过程，时间有方向；而理想弦的波动方程在时间反转下不变，是可逆的。

初始时刻在 $x=0$ 处集中加热后，温度分布随时间演化为高斯形：

$T(x,t) = \frac{Q}{\sqrt{4\pi\alpha t}}\,e^{-x^2/(4\alpha t)}$

分布的特征宽度 $\sigma = \sqrt{2\alpha t}$，即热量扩散距离与时间的平方根成正比，而非与 $t$ 成正比——这是扩散（而非匀速传播）的典型标志。

拉普拉斯方程

在无自由电荷的区域，电势 $\phi$ 满足拉普拉斯方程：

$\nabla^2 \phi = \frac{\partial^2 \phi}{\partial x^2} + \frac{\partial^2 \phi}{\partial y^2} + \frac{\partial^2 \phi}{\partial z^2} = 0$

拉普拉斯方程描述的是稳定平衡的场——没有波动，没有扩散，时间不出现。满足此方程的函数有一个重要性质：函数在任意一点的值，等于以该点为中心的任意球面上函数值的平均值。这意味着满足拉普拉斯方程的场没有孤立极值（除非在边界上），这正是导体内部无电场、静电屏蔽成立的数学根源。

格林函数：从点源到任意源

实际问题中，往往需要求解由非均匀“源”驱动的场方程，例如存在电荷分布 $\rho(\mathbf{r})$ 时的静电势。格林函数方法的核心思路是：先求单个点源产生的场（格林函数），再将任意源叠加。

点源的物理意义

静电学中，有电荷时电势满足泊松方程：

$\nabla^2 \phi = -\frac{\rho}{\varepsilon_0}$

在原点处放置点电荷 $q$，其电势为早已熟知的库仑势：

$\phi(r) = \frac{q}{4\pi\varepsilon_0 r}$

这个解就是泊松方程的格林函数（相差常数因子）：它告诉我们，当“源”是一个点时，场的形状是什么。

叠加求解任意分布

对于任意电荷密度分布 $\rho(\mathbf{r}')$，将其看作无数个点电荷的集合。位于 $\mathbf{r}'$ 处、电荷量为 $\rho(\mathbf{r}')\,dV'$ 的微元贡献电势：

$d\phi(\mathbf{r}) = \frac{\rho(\mathbf{r}')}{4\pi\varepsilon_0\,|\mathbf{r} - \mathbf{r}'|}\,dV'$

对全空间积分，得到任意电荷分布的电势：

$\phi(\mathbf{r}) = \frac{1}{4\pi\varepsilon_0}\int\frac{\rho(\mathbf{r}')}{|\mathbf{r} - \mathbf{r}'|}\,dV'$

以下是格林函数方法在不同物理场景中的对比：

变分法与最小作用量原理

费马原理：光走时间最短的路径

光从空气射入玻璃，路径发生偏折。费马原理指出，光在两点之间传播走的是所用时间最短的路径。在均匀介质中，最短时间对应直线；在界面处，最短时间路径给出折射定律：

$n_1\sin\theta_1 = n_2\sin\theta_2$

其中 $n_1$、$n_2$ 是两种介质的折射率，$\theta_1$、$\theta_2$ 是与法线的夹角。费马原理是变分法在光学中的最早应用——“在所有可能路径中找出最优的那条”正是变分法要解决的问题。

拉格朗日力学与最小作用量原理

在力学中，定义拉格朗日量：

$L = T - V$

其中 $T$ 是系统的动能，$V$ 是势能。对于从时刻 $t_1$ 运动到 $t_2$ 的力学系统，定义作用量：

$S = \int_{t_1}^{t_2} L\,dt$

最小作用量原理（哈密顿原理）指出：物理系统的真实运动轨迹，是使作用量 $S$ 取极值的那条路径，即对路径作任意微小变动时 $\delta S = 0$。

对 $\delta S = 0$ 进行推导，可得到欧拉-拉格朗日方程：

$\frac{d}{dt}\frac{\partial L}{\partial \dot{q}} - \frac{\partial L}{\partial q} = 0$

其中 $q$ 是广义坐标，$\dot{q}$ 是广义速度。这个方程等价于牛顿第二定律，但形式更普遍——无论坐标系如何选取，约束力自动消去，只需写出 $T$ 和 $V$ 即可得到运动方程。

用拉格朗日方法分析单摆

以单摆为例，摆长 $l$，摆角 $\theta$，忽略摆绳质量。

动能：$T = \dfrac{1}{2}ml^2\dot{\theta}^2$

势能（以最低点为零势能）：$V = mgl(1-\cos\theta)$

拉格朗日量：$L = \dfrac{1}{2}ml^2\dot{\theta}^2 - mgl(1-\cos\theta)$

以 $\theta$ 为广义坐标，代入欧拉-拉格朗日方程：

$\frac{d}{dt}(ml^2\dot{\theta}) + mgl\sin\theta = 0 \implies \ddot{\theta} + \frac{g}{l}\sin\theta = 0$

小角度近似下 $\sin\theta\approx\theta$，得简谐振动方程，角频率 $\omega = \sqrt{g/l}$，与牛顿力学结果完全一致。

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

定理的表述

1918年，德国数学家埃米·诺特（Emmy Noether）证明了一个深刻定理：若一个物理系统的作用量在某种连续变换下保持不变（即具有某种对称性），则必然存在一个对应的守恒量。

这三对最重要的对应关系将中学物理中分散学到的守恒律统一到了同一根源之下：

简单推导：平移对称性 $\Rightarrow$ 动量守恒

若 $L$ 不显含广义坐标 $q$（即 $\partial L/\partial q = 0$），欧拉-拉格朗日方程变为：

$\frac{d}{dt}\frac{\partial L}{\partial \dot{q}} = 0$

量 $p = \partial L/\partial \dot{q}$ 称为广义动量，上式说明它是常数（守恒量）。对于自由粒子 $L = \frac{1}{2}m\dot{x}^2$，广义动量 即线动量，不受外力时守恒——这正是牛顿第一定律的拉格朗日语言版本。

群论初步：对称操作的数学语言

把一个系统所有的对称操作（如绕某轴的旋转）收集在一起，这些操作满足一定的代数规则——任意两次操作的复合仍是对称操作，存在恒等操作，每个操作有逆操作——这样的集合在数学上称为群。

三维空间中所有旋转操作的集合记作 $SO(3)$ 群。量子力学中，氢原子的角动量量子数 $l = 0, 1, 2, 3, \ldots$ 恰好对应 $SO(3)$ 群的各个不可约表示，每个表示的维数为 $2l+1$，这给出了磁量子数 可取的个数：

这解释了为什么每个电子层中轨道的数目是 $1, 3, 5, 7, \ldots$ 这样的奇数序列，而不是任意的——这是旋转对称性在量子力学中的直接体现。

练习题

选择题

第 1 题 一根两端固定的弦，长度 $L = 1\,\text{m}$，波速 $v = 400\,\text{m/s}$，其第三泛音（$n = 4$）的频率为多少？

A. $100\,\text{Hz}$

B. $400\,\text{Hz}$

C. $800\,\text{Hz}$

D. $1600\,\text{Hz}$

第 2 题 关于热传导方程 $\partial T/\partial t = \alpha\,\partial^2 T/\partial x^2$，以下说法正确的是：

A. 热传导过程是时间可逆的，将 $t$ 替换为 $-t$ 方程形式不变

B. 初始时刻在 $x = 0$ 处集中加热，之后温度分布的特征宽度与时间 $t$ 成正比

C. 热传导方程与波动方程的主要区别在于时间导数的阶数不同

D. 热传导方程与拉普拉斯方程形式相同，物理意义也完全一样

第 3 题 对于弹簧振子（质量 $m$，劲度系数 $k$），拉格朗日量为 $L = \frac{1}{2}m\dot{x}^2 - \frac{1}{2}kx^2$，由欧拉-拉格朗日方程得到的运动方程是：

A. $m\ddot{x} - kx = 0$

B. $m\ddot{x} + kx = 0$

C. $\frac{1}{2}m\ddot{x} + kx = 0$

D. $m\dot{x} + kx = 0$

第 4 题 根据诺特定理，以下对应关系中错误的是：

A. 物理定律与时间起点无关（时间平移对称性）$\Rightarrow$ 能量守恒

B. 物理定律与空间位置无关（空间平移对称性）$\Rightarrow$ 动量守恒

C. 物理定律与空间方向无关（旋转对称性）$\Rightarrow$ 角动量守恒

D. 物理定律在匀速运动参考系中形式相同（伽利略不变性）$\Rightarrow$ 质量守恒

计算题

第 5 题 一根两端固定的弦，长度 $L = 0.8\,\text{m}$，线密度 $\mu = 5.0 \times 10^{-3}\,\text{kg/m}$，张力 $T_0 = 80\,\text{N}$。

（1）计算弦上波的传播速度 $v$（单位：$\text{m/s}$）。

（2）写出各阶驻波频率的通式 $f_n$，并计算基频 $f_1$ 和第二泛音的频率 $f_3$。

（3）若 $t = 0$ 时刻弦的位移为 $u(x,0) = 0.02\sin\!\left(\dfrac{\pi x}{L}\right)\,\text{m}$，速度为零，写出此后弦的运动表达式 。

第 6 题 对一维弹簧-滑块系统（弹簧劲度系数 $k = 50\,\text{N/m}$，滑块质量 $m = 0.2\,\text{kg}$，水平面光滑），用拉格朗日方法分析。

（1）以 $x$ 为广义坐标（弹簧原长位置为原点），写出拉格朗日量 $L$。

（2）写出欧拉-拉格朗日方程并化简，得到运动方程，指出其物理意义。

（3）计算系统的振动角频率 $\omega$ 和振动周期 $T_{\text{period}}$。

（4）若初始条件为 $x(0) = 0.1\,\text{m}$，$\dot{x}(0) = 0$，写出此后的运动方程 $x(t)$。

第三泛音对应 $n = 4$（基频 $n=1$，第一泛音 $n=2$，第二泛音 $n=3$，第三泛音 $n=4$）：

$f_4 = 200 \times 4 = 800\,\text{Hz}$

A 错：$n=1$ 时基频为 $200\,\text{Hz}$（不是 $100\,\text{Hz}$）；B 错：$n=2$，$f_2=400\,\text{Hz}$，是第一泛音；D 错：$n=8$，超出题目范围。

B 错：特征宽度 $\sigma = \sqrt{2\alpha t}$ 正比于 $\sqrt{t}$（时间平方根），而非与 $t$ 成正比，这是扩散过程与匀速传播的根本区别。

D 错：拉普拉斯方程不含时间项，描述稳态分布；热传导方程含时间项，描述动态扩散，两者物理意义截然不同。

$\frac{\partial L}{\partial \dot{x}} = m\dot{x} \implies \frac{d}{dt}(m\dot{x}) = m\ddot{x}$

$\frac{\partial L}{\partial x} = -kx$

代入得：$m\ddot{x} - (-kx) = 0$，即 $m\ddot{x} + kx = 0$。

这正是简谐运动方程，与牛顿第二定律 $ma = -kx$ 完全一致。A 错：符号有误；C 错：系数错误；D 错：方程中仅有一阶导，不描述振动。

（2）驻波频率：

两端固定弦的各阶驻波频率为：

$f_n = \frac{nv}{2L} = \frac{n \times 400}{2 \times 0.8} = 250n\,\text{Hz}, \quad n = 1, 2, 3, \ldots$

基频（$n=1$）：$f_1 = 250\,\text{Hz}$

第二泛音（$n=3$）：$f_3 = 3 \times 250 = 750\,\text{Hz}$

（3）运动表达式：

初始位移 $u(x,0) = 0.02\sin(\pi x/L)$ 恰好是 $n=1$ 驻波的空间分布，初速度为零，说明弦仅做基频振动。$n=1$ 的角频率为：

$\omega_1 = 2\pi f_1 = 2\pi \times 250 = 500\pi\,\text{rad/s}$

代入驻波表达式（初速度为零故余弦项初相为零）：

$u(x,t) = 0.02\sin\!\left(\frac{\pi x}{0.8}\right)\cos(500\pi\, t)\,\text{m}$

弦以振幅 $0.02\,\text{m}$、频率 $250\,\text{Hz}$ 做基频驻波振动，波腹在 $x = 0.4\,\text{m}$（弦中点）处，两端始终是波节。

$V = \frac{1}{2}kx^2 = \frac{1}{2}(50)x^2$

$L = T - V = \frac{1}{2}(0.2)\dot{x}^2 - \frac{1}{2}(50)x^2 = 0.1\dot{x}^2 - 25x^2$

（2）欧拉-拉格朗日方程：

$\frac{\partial L}{\partial \dot{x}} = 0.2\dot{x} \implies \frac{d}{dt}(0.2\dot{x}) = 0.2\ddot{x}$

$\frac{\partial L}{\partial x} = -50x$

代入欧拉-拉格朗日方程：

$0.2\ddot{x} - (-50x) = 0 \implies 0.2\ddot{x} + 50x = 0 \implies \ddot{x} + 250x = 0$

这是简谐振动方程（$\ddot{x} + \omega^2 x = 0$），描述弹簧振子在平衡位置附近的周期性振动，与牛顿第二定律 $m\ddot{x} = -kx$ 完全一致。

（3）角频率与周期：

$\omega = \sqrt{\frac{k}{m}} = \sqrt{\frac{50}{0.2}} = \sqrt{250} = 5\sqrt{10} \approx 15.8\,\text{rad/s}$

$T_{\text{period}} = \frac{2\pi}{\omega} = \frac{2\pi}{5\sqrt{10}} \approx 0.40\,\text{s}$

（4）运动方程：

简谐运动通解为 $x(t) = A\cos(\omega t + \varphi)$。代入初始条件 $x(0) = 0.1\,\text{m}$，$\dot{x}(0) = 0$：

$A = 0.1\,\text{m}, \quad \varphi = 0$

$x(t) = 0.1\cos(5\sqrt{10}\,t)\,\text{m} \approx 0.1\cos(15.8\,t)\,\text{m}$

滑块以振幅 $0.1\,\text{m}$、周期约 $0.40\,\text{s}$ 在平衡位置附近做简谐振动。