文章

為什麼幾乎所有能量公式都是二次型?從穩定平衡到高斯分布

彈簧位能 ½kx²、動能 ½mv²、電場能量密度 ½ε₀E²……為什麼物理的能量公式幾乎都長成平方的樣子?本文從穩定平衡點的泰勒展開出發,證明這不是巧合,並串起正定矩陣、簡正模與高斯分布背後同一個結構。

為什麼幾乎所有能量公式都是二次型?從穩定平衡到高斯分布
\[\require{physics}\]

一個奇怪的共同點

彈簧位能:

\[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{高斯漲落}}\]

這四件事,其實是同一個數學結構的不同面貌。


本文章以 CC BY 4.0 授權