python matlab 普朗克公式黑体光谱辐射出射度绘图

研志必有功

已于 2022-04-20 20:50:29 修改

阅读量1.4w

点赞数 23

分类专栏：物理文章标签： python 数学建模物理学

于 2020-11-11 00:21:00 首次发布

本文链接：https://blog.csdn.net/qq_44930937/article/details/109610045

版权

物理专栏收录该内容

1 篇文章

订阅专栏

公式

普兰克公式

python程序

import numpy as np 
import matplotlib.pyplot as plt 
#------------------------------------------------------------------------
#计算500K和800K黑体的光谱辐射出射度
def planck_formula(wavelength,      #波长
                   temperature,     #温度
                   c1=3.7414*10**8, #c1常量
                   c2=1.43879*10**4 #c2常量
                    ):
    return (c1/wavelength**5)*(1/(np.e**(c2/wavelength/temperature)-1))
#得到取样点以及输出取样点对应的值
wavelength_limit = np.linspace(0.001,100,100000)
out_500 = planck_formula(wavelength_limit,500)
out_800 = planck_formula(wavelength_limit,800)
plot_500 = plt.plot(wavelength_limit,out_500,label='500K')
plot_800 = plt.plot(wavelength_limit,out_800,label='800K')
#绘图
plt.xlim(0,20)
plt.ylim(0,8000)
plt.xlabel('wavelength: μm')
plt.ylabel('spectral radiant emission: w/(cm^2 * μm)')
plt.title('Blackbody spectral radiant emission curve:')
plt.legend()
plt.show()
#------------------------------------------------------------------------

结果

在这里插入图片描述

matlab程序

function IRwork1()
    function  emission = planck_formula(wavelength,temperature)
        c1=3.7414*10^8;
        c2=1.43879*10^4;
        emission = (c1./(wavelength().^5)).*(1./(exp(c2./wavelength./temperature)-1));
    end
wavelength = 0.1:0.01:25;
out_500 = planck_formula(wavelength,500);
out_800 = planck_formula(wavelength,800);
plot(wavelength,out_500);
hold on;
plot(wavelength,out_800);
title('Blackbody spectral radiant emission curve:');

end