写一个matlab代码，计算某一温度下某个波段内的黑体光谱辐射亮度和光谱辐射强度

这里提供一个matlab代码示例，计算某一温度下某个波段内的黑体光谱辐射亮度和光谱辐射强度。代码如下：

%% 参数设置
T = 1000; % 温度，单位：K
lambda_start = 1e-9; % 波段起始波长，单位：m
lambda_end = 5e-6; % 波段结束波长，单位：m
lambda_num = 10000; % 波长个数
h = 6.626e-34; % 普朗克常数，单位：J·s
c = 2.998e8; % 光速，单位：m/s
k = 1.381e-23; % 玻尔兹曼常数，单位：J/K

%% 计算
lambda = linspace(lambda_start, lambda_end, lambda_num); % 生成波长序列
B_lambda = 2*h*c^2./lambda.^5./(exp(h*c./(lambda*k*T))-1); % 计算黑体光谱辐射亮度
I_lambda = B_lambda*pi; % 计算光谱辐射强度

%% 绘图
figure;
plot(lambda*1e9, B_lambda);
xlabel('波长（nm）');
ylabel('光谱辐射亮度（W/(m^2·nm·sr)）');
title(['T=', num2str(T), 'K']);

figure;
plot(lambda*1e9, I_lambda);
xlabel('波长（nm）');
ylabel('光谱辐射强度（W/(m^2·nm)）');
title(['T=', num2str(T), 'K']);

使用该代码可以得到某一温度下某个波段内的黑体光谱辐射亮度和光谱辐射强度的计算结果，并且可以绘制出光谱辐射亮度和光谱辐射强度随波长的变化趋势。