Matlab实现普朗克曲线绘制与峰值定位

1. 项目背景与目标

普朗克曲线是描述黑体辐射光谱能量分布的重要物理模型，在热力学、光学和天文观测等领域有广泛应用。这个项目看似简单，实则包含了物理建模、数值计算和可视化分析三个关键环节。通过Matlab实现普朗克曲线的绘制和峰值点定位，可以帮助我们：

• 直观理解黑体辐射的波长分布特性
• 验证维恩位移定律的数值关系
• 为后续的光谱分析建立基础工具链

我在实验室做红外测温研究时，经常需要快速验证不同温度下的辐射特性。传统查表法效率低下，而用Matlab动态生成曲线可以大幅提升工作效率。下面分享我的实现方案和踩坑经验。

2. 理论基础与公式解析

2.1 普朗克定律的数学表达

普朗克辐射定律的谱辐射度公式为：

$$
B_\lambda(\lambda, T) = \frac{2hc^2}{\lambda^5} \frac{1}{e^{hc/(\lambda k_B T)} - 1}
$$

其中各参数含义：

• $h = 6.626 \times 10^{-34}$ J·s（普朗克常数）
• $c = 2.998 \times 10^8$ m/s（光速）
• $k_B = 1.381 \times 10^{-23}$ J/K（玻尔兹曼常数）
• $\lambda$：波长（单位：米）
• $T$：绝对温度（单位：开尔文）

注意：公式中的指数项容易导致数值溢出，实际编程时需要特别注意处理量纲和计算顺序。

2.2 维恩位移定律

理论预测的峰值波长位置由维恩位移定律决定：

$$
\lambda_{max} = \frac{b}{T}
$$

其中维恩常数 $b ≈ 2.898 \times 10^{-3}$ m·K。我们后续的数值计算结果应当与该理论值吻合。

3. Matlab实现详解

3.1 基础绘图实现

matlab复制function plot_planck(T_K, lambda_range)
    % 常数定义
    h = 6.626e-34; c = 2.998e8; kB = 1.381e-23;
    
    % 波长数组（纳米转米）
    lambda_m = lambda_range * 1e-9;  
    
    % 计算谱辐射度
    exponent