使用matlab进行普朗克公式计算

对进行红外仿真必经环节，通过它可以算出黑体相应温度的辐射出射度

clear all;
close all;
c1=3.741832;%第一辐射常数
c2=14387.86;%第二辐射常数
for T=270:10:370 % 设置辐射温度（K）
l=0.0001:0.02:15; % 设置波长范围及计算步长
M=1e+4.*（c1./（l.^5）./（exp（c2./（l.*T））-1））; % 计算指定温度光谱辐出度
plot（lM‘-b‘‘LineWidth‘1.4） % 绘制光谱辐出度曲线
maxM = max（M） % 找出指定温度最大光谱辐出度
i=find（maxM==M）; % 找峰值波长点
text（l（i+20）M（i+20）[num2str（T）‘K‘]‘VerticalAlignment‘...
‘baseline‘‘HorizontalAlignment‘‘left‘‘fontsize‘7）
hold on % 在指定位置按给定方式标记对应温度
end
set（gca‘XTick‘[0123456789101112131415]） %设置横坐标点
xlabel（‘\lambda / \mum‘） % 横坐标名称及单位
ylabel（‘M_{b\lambda} / W\cdotcm^{-2}\cdot\mum^{-1}‘） % 纵坐标名称及单位
hold off

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        722  2010-11-02 08:26  普朗克公式计算代码\pulangke.m

     文件        551  2009-05-13 17:23  普朗克公式计算代码\quxiannihe.m

     文件       1342  2009-05-14 09:03  普朗克公式计算代码\yingshehanshu.m

     文件        722  2010-11-02 08:24  普朗克公式计算代码\普朗克公式.m

     文件        530  2009-05-13 16:48  普朗克公式计算代码\辐射度计算.m

     目录          0  2011-04-17 10:19  普朗克公式计算代码

----------- ---------  ---------- -----  ----

                 3867                    6