数学基础 · 附录 A 微积分要点 · A9

齐次函数与欧拉定理

Homogeneous Functions and Euler's Theorem
已完成速查更新于 2026.10.08统计物理讲义 v1.0

正文用到本节的地方:§1.9(欧拉关系 U=TS−pV+μNU = TS - pV + \mu N 与吉布斯–杜安关系)、§P1.2(哈密顿量等于动能加势能)、§5.5(能量均分定理)、§19.7 与 §20.4(临界点附近的标度形式,"广义齐次函数")。

§A9.1定义与例子

若对一切 λ>0\lambda>0

f(λx1,…,λxn)=λkf(x1,…,xn)(A9.1)f(\lambda x_1,\dots,\lambda x_n) = \lambda^kf(x_1,\dots,x_n) \tag{A9.1}

则称 ff 是 kk 次齐次函数。例:x2+3xyx^2 + 3xy(k=2k = 2);x2+y2\sqrt{x^2+y^2}(k=1k = 1);x/yx/y、ex/ye^{x/y}(k=0k = 0);1/∣r∣1/\lvert\mathbf r\rvert(k=−1k = -1)。x+1x + 1 不是齐次函数。

物理中的齐次函数:

  • 广延量是其广延自变量的一次齐次函数:把系统复制 λ\lambda 份并合在一起,U(λS,λV,λN)=λU(S,V,N)U(\lambda S,\lambda V,\lambda N) = \lambda U(S,V,N)。强度量是零次齐次函数:T(λS,λV,λN)=T(S,V,N)T(\lambda S,\lambda V,\lambda N) = T(S,V,N)。(这里假设表面效应、长程力可以忽略。)
  • 动能 12∑miq˙i2\frac12\sum m_i\dot q_i^2(或任意关于速度的二次型)是速度的二次齐次函数。
  • 幂律势 V∝rnV\propto r^n 是坐标的 nn 次齐次函数:谐振势 n=2n = 2,引力与库仑势 n=−1n = -1。

反例与吉布斯佯谬。 若把理想气体的熵写成 S=NkBln⁡V+32NkBln⁡(U/N)+Ns1S = Nk_{\mathrm B}\ln V + \frac32Nk_{\mathrm B}\ln(U/N) + Ns_1(没有 1/N!1/N! 修正的结果),则 S(λU,λV,λN)=λS+λNkBln⁡λ≠λSS(\lambda U,\lambda V,\lambda N) = \lambda S + \lambda Nk_{\mathrm B}\ln\lambda\ne\lambda S:它不是一次齐次的,熵不是广延量。这正是吉布斯佯谬(§3.7)的数学表现;把 ln⁡V\ln V 换成 ln⁡(V/N)\ln(V/N) 后恢复了齐次性。§1.11 推导 (1.27) 时正是用齐次性定出了积分常数对 NN 的依赖。

§A9.2欧拉定理

定理:若 ff 是 kk 次齐次函数(且可微),则

∑i=1nxi∂f∂xi=kf(A9.2)\sum_{i=1}^nx_i\frac{\partial f}{\partial x_i} = kf \tag{A9.2}

证明:把 (A9.1) 两边对 λ\lambda 求导。左边是复合函数,第 ii 个自变量 λxi\lambda x_i 对 λ\lambda 的导数是 xix_i,由链式法则 (A1.3) 得 ∑ixi (∂if)(λx)\sum_ix_i\,(\partial_if)(\lambda\mathbf x);右边是 kλk−1f(x)k\lambda^{k-1}f(\mathbf x)。令 λ=1\lambda = 1 即得 (A9.2)。□\square

逆定理:若 (A9.2) 对一切 x\mathbf x 成立,则 ff 是 kk 次齐次的。证明:令 ϕ(λ)=f(λx)/λk\phi(\lambda) = f(\lambda\mathbf x)/\lambda^k。在点 λx\lambda\mathbf x 处用 (A9.2),得 ∑ixi (∂if)(λx)=kf(λx)/λ\sum_ix_i\,(\partial_if)(\lambda\mathbf x) = kf(\lambda\mathbf x)/\lambda,于是

ϕ′(λ)=∑ixi(∂if)(λx)λk−kf(λx)λk+1=0\phi'(\lambda) = \frac{\sum_ix_i(\partial_if)(\lambda\mathbf x)}{\lambda^k} - \frac{kf(\lambda\mathbf x)}{\lambda^{k+1}} = 0

ϕ\phi 是常数,等于 ϕ(1)=f(x)\phi(1) = f(\mathbf x),这就是 (A9.1)。

检验:f=x2+3xyf = x^2 + 3xy,xfx+yfy=x(2x+3y)+y⋅3x=2x2+6xy=2fxf_x + yf_y = x(2x + 3y) + y\cdot3x = 2x^2 + 6xy = 2f。

§A9.3齐次函数的导数

若 ff 是 kk 次齐次的,则每个偏导数 ∂f/∂xj\partial f/\partial x_j 是 k−1k-1 次齐次的:把 (A9.1) 两边对 xjx_j 求导,左边由链式法则多出一个因子 λ\lambda,

λ (∂jf)(λx)=λk ∂jf(x)⟹(∂jf)(λx)=λk−1∂jf(x)(A9.3)\lambda\,(\partial_jf)(\lambda\mathbf x) = \lambda^k\,\partial_jf(\mathbf x)\quad\Longrightarrow\quad(\partial_jf)(\lambda\mathbf x) = \lambda^{k-1}\partial_jf(\mathbf x) \tag{A9.3}

热力学推论:T=∂U/∂ST = \partial U/\partial S、−p=∂U/∂V-p = \partial U/\partial V、μ=∂U/∂N\mu = \partial U/\partial N 都是零次齐次的,即强度量。取 λ=1/N\lambda = 1/N:

T(S,V,N)=T(SN,VN,1)≡T(s,v)T(S,V,N) = T\left(\frac SN,\frac VN,1\right)\equiv T(s,v)

强度量只依赖于"密度"(每个粒子的熵 ss、体积 vv),与系统大小无关。这就是 §1.9 末尾写 dμ=−s dT+v dpd\mu = -s\,dT + v\,dp 的根据。

§A9.4热力学:欧拉关系与吉布斯–杜安关系

U(S,V,N)U(S,V,N) 一次齐次,由 (A9.2):U=S ∂SU+V ∂VU+N ∂NU=TS−pV+μNU = S\,\partial_SU + V\,\partial_VU + N\,\partial_NU = TS - pV + \mu N,即 (1.23)。

只对部分变量齐次。F(T,V,N)F(T,V,N) 中 TT 是强度量,不随系统放大而改变,所以 FF 只对 (V,N)(V,N) 一次齐次:F(T,λV,λN)=λF(T,V,N)F(T,\lambda V,\lambda N) = \lambda F(T,V,N)。把欧拉定理只用于被放大的变量:

V∂F∂V+N∂F∂N=F⟹F=−pV+μN(A9.4)V\frac{\partial F}{\partial V} + N\frac{\partial F}{\partial N} = F\quad\Longrightarrow\quad F = -pV + \mu N \tag{A9.4}

与 G=F+pV=μNG = F + pV = \mu N 一致。同理 G(T,p,N)=N ∂G/∂N=μNG(T,p,N) = N\,\partial G/\partial N = \mu N;多组分时 G=∑iμiNiG = \sum_i\mu_iN_i。

例:理想气体。F=−NkBT[ln⁡VNλ3+1]F = -Nk_{\mathrm B}T\left[\ln\frac{V}{N\lambda^3} + 1\right](λ\lambda 为热波长,只依赖于 TT)。∂F/∂V=−NkBT/V=−p\partial F/\partial V = -Nk_{\mathrm B}T/V = -p;∂F/∂N=−kBTln⁡VNλ3=μ\partial F/\partial N = -k_{\mathrm B}T\ln\frac{V}{N\lambda^3} = \mu(与 (3.17) 一致)。于是 −pV+μN=−NkBT−NkBTln⁡VNλ3=F-pV + \mu N = -Nk_{\mathrm B}T - Nk_{\mathrm B}T\ln\frac{V}{N\lambda^3} = F,(A9.4) 成立。

吉布斯–杜安关系:对 (1.23) 求全微分,dU=T dS+S dT−p dV−V dp+μ dN+N dμdU = T\,dS + S\,dT - p\,dV - V\,dp + \mu\,dN + N\,d\mu;减去 dU=T dS−p dV+μ dNdU = T\,dS - p\,dV + \mu\,dN 即得 (1.25)。它表达的正是 §A9.3 的结论:三个强度量只依赖于两个密度,因而彼此不独立。

§A9.5力学:哈密顿量、能量均分与位力定理

(a) 哈密顿量就是能量(§P1.2)。动能 TT 是 q˙\dot q 的二次齐次函数,pi=∂T/∂q˙ip_i = \partial T/\partial\dot q_i(势能不含速度),由 (A9.2) 取 k=2k = 2:∑ipiq˙i=2T\sum_ip_i\dot q_i = 2T。所以 H=∑ipiq˙i−L=2T−(T−V)=T+VH = \sum_ip_i\dot q_i - L = 2T - (T - V) = T + V。

(b) 能量均分(§5.5)。广义均分定理 (5.15) 说 ⟨xi ∂H/∂xj⟩=δijkBT\langle x_i\,\partial\mathcal H/\partial x_j\rangle = \delta_{ij}k_{\mathrm B}T。若 H\mathcal H 是 ff 个变量的二次齐次函数,对 i=ji = j 求和并用 (A9.2):2⟨H⟩=fkBT2\langle\mathcal H\rangle = fk_{\mathrm B}T。更一般地,若 H\mathcal H 中依赖于某组 f′f' 个变量的部分 H1\mathcal H_1 是这些变量的 ss 次齐次函数(其余部分不含这些变量),同样的论证给出 ⟨H1⟩=f′kBT/s\langle\mathcal H_1\rangle = f'k_{\mathrm B}T/s。例:极端相对论性粒子 ε=c∣p∣\varepsilon = c\lvert\mathbf p\rvert 是 p\mathbf p 的一次齐次函数,所以 ⟨ε⟩=3kBT\langle\varepsilon\rangle = 3k_{\mathrm B}T(而非 32kBT\frac32k_{\mathrm B}T)。

(c) 位力定理。设粒子在有界区域内运动,势能 VV 是坐标的 nn 次齐次函数。考虑 G=∑ipi⋅riG = \sum_i\mathbf p_i\cdot\mathbf r_i:

dGdt=∑ip˙i⋅ri+∑ipi⋅r˙i=−∑iri⋅∇iV+2T=−nV+2T\frac{dG}{dt} = \sum_i\dot{\mathbf p}_i\cdot\mathbf r_i + \sum_i\mathbf p_i\cdot\dot{\mathbf r}_i = -\sum_i\mathbf r_i\cdot\nabla_iV + 2T = -nV + 2T

(第一项用了牛顿方程 p˙i=−∇iV\dot{\mathbf p}_i = -\nabla_iV 与 (A9.2),第二项用了 p⋅r˙=mr˙2=2×\mathbf p\cdot\dot{\mathbf r} = m\dot r^2 = 2\times动能。)对很长的时间取平均,有界运动中 GG 有界,dG/dt‾=[G(τ)−G(0)]/τ→0\overline{dG/dt} = [G(\tau) - G(0)]/\tau\to0,所以

2T‾=nV‾(A9.5)2\overline T = n\overline V \tag{A9.5}

谐振子(n=2n = 2):T‾=V‾\overline T = \overline V。引力或库仑系统(n=−1n = -1):2T‾=−V‾2\overline T = -\overline V,总能量 E=T‾+V‾=−T‾E = \overline T + \overline V = -\overline T。于是自引力系统(例如恒星)失去能量(辐射)时 T‾\overline T 反而增加——它"越冷却越热",热容为负。这类系统的熵不是能量的凹函数,正则系综与微正则系综不等价(参见 §A5.4)。第17章"位力展开"的名称来自克劳修斯(1870)定义的位力 −12∑iri⋅Fi-\frac12\sum_i\mathbf r_i\cdot\mathbf F_i(他的定理说:动能的时间平均等于位力)。

§A9.6广义齐次函数与标度

若存在指数 a1,…,ana_1,\dots,a_n 使

f(λa1x1,…,λanxn)=λkf(x1,…,xn)(A9.6)f(\lambda^{a_1}x_1,\dots,\lambda^{a_n}x_n) = \lambda^kf(x_1,\dots,x_n) \tag{A9.6}

则称 ff 为广义齐次函数。与 §A9.2 同样的论证给出广义欧拉定理 ∑iaixi ∂if=kf\sum_ia_ix_i\,\partial_if = kf,并且 ∂f/∂xj\partial f/\partial x_j 满足同样形式的关系,"次数"变为 k−ajk - a_j。

临界标度(§19.7、§20.4)。重整化群给出的 (20.5) 可以写成 fs(λytt,λyhh)=λdfs(t,h)f_{\mathrm s}(\lambda^{y_t}t,\lambda^{y_h}h) = \lambda^df_{\mathrm s}(t,h)(令 λ=b\lambda = b),即 fsf_{\mathrm s} 是 (t,h)(t,h) 的广义齐次函数,at=yta_t = y_t,ah=yha_h = y_h,k=dk = d。广义齐次函数的关键用法是选择 λ\lambda 使某个自变量变成常数:取 λ=∣t∣−1/yt\lambda = \lvert t\rvert^{-1/y_t},则 λytt=±1\lambda^{y_t}t = \pm1,

fs(t,h)=∣t∣d/ytfs(±1, h∣t∣−yh/yt)f_{\mathrm s}(t,h) = \lvert t\rvert^{d/y_t}f_{\mathrm s}\big(\pm1,\ h\lvert t\rvert^{-y_h/y_t}\big)

这就是维多姆形式 (19.21):一个两变量的函数被化为"幂律 × 一个单变量函数",所有临界指数都来自"次数"的记账。例如 m=−∂fs/∂hm = -\partial f_{\mathrm s}/\partial h 的次数是 d−yhd - y_h,零场时 m∝∣t∣(d−yh)/ytm\propto\lvert t\rvert^{(d-y_h)/y_t},所以 β=(d−yh)/yt\beta = (d - y_h)/y_t。

量纲分析是同一思想的另一种形式:物理规律在改变单位时形式不变,所以物理量必须是各参数的(广义)齐次组合。例如热波长 λ=h/2πmkBT\lambda = h/\sqrt{2\pi mk_{\mathrm B}T} 是由 hh、mm、kBTk_{\mathrm B}T 唯一能组成的长度(除数值因子外;2π\sqrt{2\pi} 是约定)。

自测题

  1. 判断下列函数是否齐次,若是,求次数:(a) x3+y3+z3−3xyzx^3 + y^3 + z^3 - 3xyz;(b) ln⁡x−ln⁡y\ln x - \ln y;(c) (x2+y2)−1/2(x^2 + y^2)^{-1/2};(d) x2+yx^2 + y。[答:3;0;−1-1;不是。]
  2. 一维经典粒子在势 V=A∣x∣nV = A\lvert x\rvert^n 中,与温度 TT 的热源接触。证明 ⟨V⟩=kBT/n\langle V\rangle = k_{\mathrm B}T/n。[提示:⟨x V′(x)⟩=kBT\langle x\,V'(x)\rangle = k_{\mathrm B}T,而 xV′=nVxV' = nV;或直接计算 qx∝∫e−βA∣x∣ndx∝β−1/nq_x\propto\int e^{-\beta A\lvert x\rvert^n}dx\propto\beta^{-1/n}。]
  3. 证明极端相对论性理想气体的 CV=3NkBC_V = 3Nk_{\mathrm B}。
  4. 由 (20.5) 在 t=0t = 0 处取 λ=h−1/yh\lambda = h^{-1/y_h},证明 m∝h(d−yh)/yhm\propto h^{(d-y_h)/y_h},从而 δ=yh/(d−yh)\delta = y_h/(d - y_h)。
  5. 用位力定理估计氢原子基态:库仑势 n=−1n = -1,基态能量 E=−13.6E = -13.6 eV。求动能与势能的平均值。[答:T‾=13.6\overline T = 13.6 eV,V‾=−27.2\overline V = -27.2 eV。(量子力学中位力定理对定态同样成立。)]