matlab代码绘制普朗克分布定律将绝对黑体的辐射能力(发射能量的数量)定义为温度和波长函数: 在c1 = 3.742х10 ^ w - 8μm ^ 4 / m ^ 2, c2 = 1.439x10 ^ 4µm - k, t—开尔文温度(k)或μm波长。 编写一个脚本文件，允许创建一个名为“Рlank's law”的图形窗口(与图中所描述的最大相似性)。 l需要使用对数刻度(例如，团队可以做到:set(gca,'xscale','log'))。 使用的字体是Courier黑体，字体大小为10pt，图形窗口的颜色为浅蓝色（使用RGB代码）。

下面是用 MATLAB 代码绘制普朗克分布定律的实现过程：

% 定义常数
c1 = 3.742e-16; % m^2
c2 = 1.439e4;   % um K

% 定义波长范围
lambda = logspace(-7, -3, 1000); % m

% 定义温度范围
T = [300, 600, 900, 1200, 1500, 1800]; % K

% 计算辐射能力
B = zeros(length(T), length(lambda));
for i = 1:length(T)
    B(i,:) = (2 * pi * c1) ./ (lambda .^ 5 .* (exp(c2 ./ (lambda * T(i))) - 1));
end

% 绘图
figure('Name', 'Plank''s law', 'Color', [0.7, 0.85, 1]);
set(gca, 'XScale', 'log', 'FontName', 'Courier', 'FontSize', 10);
hold on;
for i = 1:length(T)
    plot(lambda * 1e6, B(i,:) * 1e-6, 'LineWidth', 1.5);
end
xlabel('Wavelength (\mum)');
ylabel('Spectral radiance (W sr^{-1} m^{-3})');
legend(cellstr(num2str(T', '%-4d K')), 'Location', 'NorthEast');

解释一下代码：

首先定义了两个常数：c1和c2，分别对应普朗克分布定律中的常数。然后定义了波长范围和温度范围，使用logspace函数生成了一个从10^-7到10^-3的1000个等比数列，表示波长范围；温度范围则是一个包含6个元素的向量，表示6个不同的温度。接下来使用一个循环计算了不同温度下的辐射能力B，并将其存储在一个矩阵中。最后绘制图形，使用了set函数设置了对数刻度、字体和字体大小等属性，使用plot函数绘制了不同温度下的曲线，并使用legend函数添加了图例。