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