為什麼幾乎所有能量公式都是二次型?從穩定平衡到高斯分布
彈簧位能 ½kx²、動能 ½mv²、電場能量密度 ½ε₀E²……為什麼物理的能量公式幾乎都長成平方的樣子?本文從穩定平衡點的泰勒展開出發,證明這不是巧合,並串起正定矩陣、簡正模與高斯分布背後同一個結構。
一個奇怪的共同點
彈簧位能:
\[U=\frac12kx^2\]平移動能:
\[T=\frac12mv^2\]電磁場的能量密度:
\[u=\frac{\epsilon_0}{2}E^2+\frac{1}{2\mu_0}B^2\]統計力學中的高斯分布:
\[p(x)\propto\exp\left(-\frac{x^2}{2\sigma^2}\right)\]這些式子來自完全不同的領域——力學、電磁學、統計力學——卻都長成同一個樣子:某個量的平方。這是巧合,還是背後有一個共同的數學原因?
答案是:這不是巧合,而是穩定平衡這個條件必然導出的結果。
直覺:平衡點附近,一次項必須消失
想像一個系統的位能是某個座標 $q$ 的函數 $U(q)$。系統受到的力是
\[F=-\dv{U}{q}\]在平衡點 $q_0$,系統不受力,所以
\[\left.\dv{U}{q}\right|_{q_0}=0\]現在把 $U$ 在 $q_0$ 附近做泰勒展開:
\[U(q)=U(q_0)+\underbrace{\left.\dv{U}{q}\right|_{q_0}}_{=0}(q-q_0)+\frac12\left.\dv[2]{U}{q}\right|_{q_0}(q-q_0)^2+\cdots\]一次項因為平衡條件自動消失。所以只要 $U$ 夠平滑,展開式裡第一個不是常數的項,幾乎必然是二次項:
\[U-U_0\approx\frac12k(q-q_0)^2,\qquad k=\left.\dv[2]{U}{q}\right|_{q_0}\]如果這個平衡點是穩定的(系統偏離後會被推回來,而不是愈跑愈遠),則 $k>0$。這正是簡諧振子的位能形式——彈簧公式 $U=\frac12kx^2$ 不是 Hooke 定律的巧合,而是任何系統在穩定平衡附近,泰勒展開的必然結果。
任何足夠光滑的系統,在穩定平衡點附近,看起來都像簡諧振子。
正式推導:從一維到多自由度的二次型
多自由度:Hessian 矩陣登場
如果系統有 $n$ 個座標
\[\vb{q}=\begin{pmatrix}q_1\\q_2\\\vdots\\q_n\end{pmatrix}\]同樣在平衡點附近展開,一次項(梯度)為零,留下的二次項是
\[U(\vb{q})\approx U_0+\frac12\sum_{i,j}H_{ij}q_iq_j\]其中
\[H_{ij}=\pdv[2]{U}{q_i}{q_j}\]是 Hessian 矩陣——由所有二階偏導數組成。寫成矩陣形式:
\[\boxed{U=\frac12\vb{q}^{\mathsf T}K\vb{q}}\]這正是線性代數裡的二次型(quadratic form)。例如兩個自由度耦合:
\[U=\frac12\left(k_1x^2+2k_{12}xy+k_2y^2\right)\]交叉項 $xy$ 代表兩個自由度互相耦合。而穩定平衡這個物理條件,翻譯成線性代數語言就是:矩陣 $K$ 必須是正定矩陣——這正好呼應正定矩陣一文的核心主題:正定矩陣就是「所有方向都往上凸」的矩陣,對應的正是「所有方向都穩定」的平衡點。
把 $K$ 對角化:
\[K=PDP^{\mathsf T}\]換成新座標(也就是簡正模,normal modes)之後,能量會拆成互不干擾的獨立項:
\[U=\frac12\sum_i\lambda_iQ_i^2\]一個複雜的耦合系統,在小振動近似下,就這樣拆解成一堆彼此獨立的簡諧振子——這是分子振動、晶格振動、聲學工程裡反覆出現的核心結構。
動能為什麼也是二次型?
平移動能 $T=\frac12mv^2$,多維寫成
\[T=\frac12m\vb{v}^{\mathsf T}\vb{v}\]更一般地,用廣義座標表示:
\[\boxed{T=\frac12\dot{\vb{q}}^{\mathsf T}M(\vb{q})\dot{\vb{q}}}\]這背後有一個幾何原因:速度是一個向量,而一個純量能量,若要不依賴座標軸方向、又要在 $\vb{v}\to-\vb{v}$ 時保持不變,最簡單的非平凡構造就是向量與自己的內積:
\[\vb{v}\vdot\vb{v}=v^2\]也就是二次的。換句話說:
長度平方,本來就是向量空間裡最自然的純量。
在一般曲線座標(微分幾何語言)中,這件事寫成
\[T=\frac12g_{ij}\dot q^i\dot q^j\]動能本質上就是「組態空間上的度量所定義的速度平方」。
場的能量也是「場或梯度的平方」
同樣的邏輯延伸到連續系統。電磁場:
\[u=\frac12\epsilon_0E^2+\frac{1}{2\mu_0}B^2\]弦的振動能量:
\[E=\int\left[\frac12\mu\left(\pdv{y}{t}\right)^2+\frac12T\left(\pdv{y}{x}\right)^2\right]\dd{x}\]自由純量場的拉格朗日密度:
\[\mathcal L=\frac12\partial_\mu\phi\,\partial^\mu\phi-\frac12m^2\phi^2\]背後的道理完全一致:在某個背景場 $\phi_0$ 附近做小擾動 $\phi=\phi_0+\delta\phi$,展開能量泛函
\[E[\phi]=E[\phi_0]+\delta E+\frac12\delta^2E+\cdots\]若 $\phi_0$ 是平衡解,一次變分 $\delta E=0$,所以最低階非平凡項又是
\[\delta E\sim(\delta\phi)^2\]線性化理論與二次能量幾乎永遠一起出現,這就是原因。
例題一:單擺的小角度近似
單擺的位能
\[U(\theta)=mgl(1-\cos\theta)\]不是二次型——它是三角函數。但在 $\theta=0$ 附近展開:
\[1-\cos\theta=\frac12\theta^2-\frac{1}{24}\theta^4+\cdots\]保留最低階非零項:
\[U\approx\frac12mgl\theta^2\]於是回復力矩
\[\tau=-\dv{U}{\theta}\approx-mgl\theta\]運動方程式線性化為
\[mL^2\ddot\theta+mgl\theta=0\]這就是簡諧運動方程式。單擺的完整理論是非線性的(大角度時週期會隨振幅改變),但在平衡點附近,二次近似已經足夠精準——這正是「二次型是平衡點附近的普遍語言」這句話的具體示範。
例題二:兩個耦合彈簧的簡正模
考慮兩個質量 $m$ 的物體,中間與兩側各接一根彈簧(力常數皆為 $k$),位移分別為 $x_1,x_2$。總位能是
\[U=\frac12kx_1^2+\frac12k(x_2-x_1)^2+\frac12kx_2^2\]寫成矩陣形式:
\[U=\frac12\begin{pmatrix}x_1&x_2\end{pmatrix}\begin{pmatrix}2k&-k\\-k&2k\end{pmatrix}\begin{pmatrix}x_1\\x_2\end{pmatrix}\]這個 $2\times2$ 矩陣的特徵值是 $k$ 和 $3k$,對應的簡正模分別是「兩物體同向擺動」($x_1=x_2$)和「兩物體反向擺動」($x_1=-x_2$)。在這兩組新座標下,能量各自變成互不干擾的單一二次項——這就是對角化的物理意義:把一個看起來複雜的耦合系統,拆成幾個獨立的簡諧振子。
二次能量 ⟺ 線性運動方程
還有一個值得記住的對應關係:能量是二次的,力就是一次的,運動方程式就是線性的。
單自由度:
\[U=\frac12kx^2\quad\Longrightarrow\quad F=-kx\quad\Longrightarrow\quad m\ddot x+kx=0\]多自由度:
\[U=\frac12\vb{q}^{\mathsf T}K\vb{q}\quad\Longrightarrow\quad M\ddot{\vb{q}}+K\vb{q}=0\]線性系統帶來一整套好用的工具:疊加原理、簡正模態、傅立葉分析、頻譜分解——這是為什麼物理學總是優先從二次型開始:它是非平凡、卻仍然可解的最低階理論。
延伸連結:二次能量與高斯分布
統計力學中,系統處於某狀態的機率正比於玻茲曼因子:
\[p\propto e^{-\beta E}\]如果平衡點附近 $E-E_0=\frac12\vb{x}^{\mathsf T}K\vb{x}$,機率分布就自動變成高斯分布:
\[p(\vb{x})\propto\exp\left(-\frac{\beta}{2}\vb{x}^{\mathsf T}K\vb{x}\right)\]這也是為什麼卡方統計量
\[\chi^2=\sum_i\frac{(x_i-\mu_i)^2}{\sigma_i^2}\]和最小平方法
\[\min\lVert A\vb{x}-\vb{b}\rVert^2\]都長著同一張「平方」的臉——它們全部來自同一套「平衡點附近二次近似」的邏輯。
但不是所有能量本質上都是二次的
萬有引力位能
\[U(r)=-\frac{GMm}{r}\]和 Coulomb 位能
\[U(r)=\frac{kq_1q_2}{r}\]都不是二次型;它們是完整的非線性理論。二次型出現的地方,永遠是「在某個平衡點附近」,而不是理論本身處處都是二次的:
\[\boxed{\text{完整理論通常是非線性的}\ \longrightarrow\ \text{平衡點附近的最低階近似通常是二次的}}\]總結
| 領域 | 二次型的具體樣貌 | 平衡點/展開點 |
|---|---|---|
| 力學位能 | $U=\frac12kx^2$ | 力平衡點(一次項消失) |
| 多自由度位能 | $U=\frac12\vb{q}^{\mathsf T}K\vb{q}$,$K$ 為 Hessian | 各方向皆穩定 → $K$ 正定 |
| 動能 | $T=\frac12m\vb{v}^{\mathsf T}\vb{v}$ | 向量長度平方,與座標軸無關 |
| 場論 | $\mathcal L\sim(\partial\phi)^2-m^2\phi^2$ | 背景場的小擾動 |
| 統計力學 | $p\propto e^{-\beta\cdot\text{二次型}}$ | 玻茲曼分布 + 二次能量 |
一句話濃縮整篇文章:能量在穩定平衡點是極小值,所以一次變分消失,泰勒展開的第一個非平凡項通常就是二次項;而二次型又正好描述向量的「平方長度」,因此動能、位能、場能與微小擾動的能量都普遍呈現二次型。
\[\boxed{\text{二次型}=\text{局部幾何}=\text{線性化物理}=\text{高斯漲落}}\]這四件事,其實是同一個數學結構的不同面貌。
→ 歡迎預約家教課,直接針對你的問題討論