2. 维恩位移定律的数学推导：从普朗克公式到峰值波长与温度的关系
说实话，很多教材一上来就扔出维恩位移定律的公式，然后告诉你「记住就行」。我不喜欢这种教法。你想想看，一个公式背后一定有它的物理逻辑。今天我们就从普朗克黑体辐射公式出发，一步步推导出维恩位移定律。这个过程不复杂，但能让你真正理解——为什么温度越高，峰值波长越短。
2.1 普朗克公式：一切的基础
1900年，普朗克提出了一个划时代的公式。它描述了黑体在温度 \( T \) 下，辐射能量密度随波长 \( \lambda \) 的分布：
\[
u(\lambda, T) = \frac{8\pi h c}{\lambda^5} \cdot \frac{1}{e^{hc/(\lambda k_B T)} - 1}
\]其中：
- \( h \)：普朗克常数（6.626×10⁻³⁴ J·s）
- \( c \)：光速（2.998×10⁸ m/s）
- \( k_B \)：玻尔兹曼常数（1.381×10⁻²³ J/K）
- \( T \)：绝对温度（K）
- \( \lambda \)：波长（m）
这个公式看着有点吓人，但核心就两个部分：前面的 \( 1/\lambda^5 \) 是经典物理的瑞利-金斯尾巴，后面的指数项是普朗克的量子修正。我当年第一次看到这个公式时，最困惑的就是那个「-1」——为什么不是简单的指数衰减？后来做项目才明白，没有这个「-1」，高温下的辐射能量就会发散，物理上根本说不通。
核心要点：普朗克公式是维恩位移定律的母公式。我们接下来要做的，就是找到这个函数在哪个波长处取最大值。
2.2 寻找峰值：求导的艺术
要找峰值，说白了就是求极值。我们对 \( u(\lambda, T) \) 关于 \( \lambda \) 求导，令导数为零：
\[
\frac{\partial u(\lambda, T)}{\partial \lambda} = 0
\]这个过程有点繁琐，我直接给你看关键步骤。先引入一个无量纲变量 \( x = \frac{hc}{\lambda k_B T} \)，这样公式会清爽很多：
\[
u(x) \propto \frac{x^5}{e^x - 1}
\]求导后得到：
\[
\frac{d}{dx} \left( \frac{x^5}{e^x - 1} \right) = 0
\]展开后：
\[
\frac{5x^4(e^x - 1) - x^5 e^x}{(e^x - 1)^2} = 0
\]分子为零：
\[
5(e^x - 1) - x e^x = 0
\]整理一下：
\[
(5 - x) e^x = 5
\]嗯，到这里就卡住了。这是一个超越方程，没有解析解。怎么办？数值求解。
我的经验：遇到超越方程别慌。我在做红外探测器项目时，经常需要解这类方程。用牛顿迭代法，几行代码就搞定了。关键是理解物理意义，而不是死磕解析解。
2.3 数值求解：得到关键常数
方程 \( (5 - x) e^x = 5 \) 的解是：
\[
x \approx 4.9651
\]回代到 \( x = \frac{hc}{\lambda_{\text{max}} k_B T} \)，得到：
\[
\lambda_{\text{max}} T = \frac{hc}{4.9651 k_B}
\]代入常数：
\[
\lambda_{\text{max}} T = \frac{(6.626 \times 10^{-34})(2.998 \times 10^8)}{4.9651 \times 1.381 \times 10^{-23}} \approx 2.898 \times 10^{-3} \, \text{m·K}
\]这就是维恩位移定律：
\[
\lambda_{\text{max}} T = b, \quad b \approx 2.898 \times 10^{-3} \, \text{m·K}
\]其中 \( b \) 就是维恩位移常数。
| 参数 | 符号 | 数值 | 单位 |
|---|---|---|---|
| 普朗克常数 | \( h \) | 6.626×10⁻³⁴ | J·s |
| 光速 | \( c \) | 2.998×10⁸ | m/s |
| 玻尔兹曼常数 | \( k_B \) | 1.381×10⁻²³ | J/K |
| 维恩位移常数 | \( b \) | 2.898×10⁻³ | m·K |
注意：这里的 \( \lambda_{\text{max}} \) 是辐射能量密度最大的波长。如果你用的是辐射通量密度（每单位波长），常数会略有不同。我曾经在项目里吃过这个亏——用了错误的常数，导致计算结果偏差了约 1.5%。
2.4 推导过程的可视化
为了让你更直观地理解整个推导逻辑，我画了一张流程图：
2.5 推导中的常见陷阱
我在实际教学中发现，很多人在推导时会犯几个典型错误：
- 求导对象搞错：有人对 \( \lambda \) 求导，有人对 \( \nu \) 求导，结果常数不一样。记住，我们用的是波长形式。
- 忽略常数因子：求导时前面的系数 \( 8\pi hc \) 是常数，可以直接扔掉，不影响极值点位置。
- 数值求解精度不够：\( x = 4.9651 \) 这个值，如果你只取 4.97，常数偏差约 0.1%。在工程上可能够用，但在科研中不行。
避坑指南：我曾经在给一个红外测温系统做校准时，用了 2.9×10⁻³ 这个近似值。结果在 1000K 时，峰值波长偏差了 0.7μm，直接导致测温误差超过 5°C。后来我改用 2.898×10⁻³，问题才解决。所以，常数能精确就精确。
2.6 从推导到工程直觉
推导完了，我们来建立一点直觉。维恩位移定律告诉我们：
- 温度翻倍，峰值波长减半
- 太阳表面（约 5778K）的峰值在 500nm 左右——正好是可见光波段
- 人体（约 310K）的峰值在 9.3μm——远红外波段
你想想看，为什么人类的眼睛进化到对可见光敏感？因为太阳的辐射峰值就在这个波段。这不是巧合，这是维恩位移定律在生物进化中的体现。
工程启示：设计红外探测器时，先算目标温度对应的峰值波长，再选探测器材料。比如探测 300K 的物体，峰值在 10μm 左右，那就选 HgCdTe 或微测辐射热计。选错了材料，灵敏度直接打对折。
好了，推导到这里就结束了。从普朗克公式到维恩位移定律，核心就是三步：引入变量、求导、数值求解。下次你看到 \( \lambda_{\text{max}} T = b \) 这个公式时，应该能想起它背后的物理逻辑——而不是死记硬背一个常数。