黑体辐射出射度计算

根据普朗克定律计算给定温度黑体在某波长范围内的辐射出射度并绘图。

pro planck根据planck定律计算黑体辐射出射度并绘图

设定起始波长，终止波长，黑体温度等参数

ws=0.3

we=2.5

interval=0.01

t=6000

计算波长数目，创建波长数组

nums=ceil((we-ws)/interval)

wv=findgen(nums)*interval ws

计算出射度

M=Cal_emittance(wv,t)

绘制辐射出射度图

plot,wv,M,xtitle='Wavelength(um)',ytitle='M(J/s.cm^2.um)',$

title='Black body radiation',xrange=[min(wv),max(wv)],/nodata,$

xstyle=1,ystyle=16,color='000000'Xl,background='FFFFFF'xl

oplot,wv,M,color='0000EF'xl,thick=1.5

end

 Cal_emittance函数来计算辐射出射度

h=6.63e-34

c=2.998e8

k=1.38e-23

wv_m=wv*1e-6

普朗克定律计算公式，然后把结果的单位转化

M=2*!pi*double(h)*(c^2)/(wv_m^5*(exp(h*c/(wv_m*k*temp))-1))

return,M*1e-10

end

代码链接：https://pan.baidu.com/s/1lBhw-8XqaYWdYtGMVr0Xfg 

提取码：vsk4