黑体辐射规律的研究

1、黑体辐射规律的研究湖南大学 XX院系 XX 专业 XX年级姓名 学号问题 黑体辐射的规律在任何温度下对任意波长的电磁波只吸收不反射的物体称为绝对黑体，简称黑体。根据实验得出两个实验规律。黑体的总辐射本领(能力)为E(T) = T4，这就是斯特潘-玻尔兹曼定律，其中， = 5.67×10-8W/(m2·K4)，称为斯特潘常数。黑体的单色辐射本领(能力)的峰值波长与温度的关系为Tm = b，这就是维恩位移定律，其中，b = 2.897×10-3m·K，b称为维恩常数。根据普朗克提出的黑体辐射公式，计算斯特潘常数和维恩常数。以温度为参数，单色辐射本领与波长的曲

2、线有什么特点，峰值与波长的关系曲线有什么特点？数学模型黑体的单色辐射本领是在单位时间内从物体表面单位面积上所发射的波长在到 + d范围内的辐射能量dE(,T)与波长间隔d之比， (14.1.1)M(,T)表示在单位时间内从物体表面单位面积发射的波长在附近单位波长间隔内的辐射本领，是波长和温度的函数，其单位是W/m3。普朗克提出的黑体单色辐射本领的公式为， (14.1.2)其中，k为玻尔兹曼常数，h为普朗克常数，c为真空中的光速。对波长从0到无穷大积分就得总辐射本领，即：黑体单位面积辐射能量的功率， (14.1.3)设 ， (14.1.4)(14.1.2)式可化为， (14.1.5)由(14.1

3、.4)式得，所以(14.1.3)式可化为， (14.1.6)其中，为常数，I为积分。 (14.1.7)手工计算积分I比较麻烦，其步骤如下，设y = nx，可得。由此可得CI = 5.6688×10-8。 (14.1.8)这就是斯特潘常数。当波长趋于零时，单色辐射本领趋于零；当波长趋于无穷大时，单色辐射本领也趋于零。因此单色辐射本领随波长的变化有极值。令dM (x,T)/dx = 0，可得方程：x = 5(1 - e-x)， (14.1.9)一般用迭代算法计算上式之值，除了零解之外，可得xm的值为4.965，由(14.1.4)式可得维恩常数= 0.0029。 (14.1.10)这就是维

4、恩常数。理论值与实验值符号得很好。算法用MATLAB的符号积分计算积分可直接计算I的值，从而计算斯特潘常数。设， (14.1.5*)用MATLAB的符号导数先求y的导数，求导数的零解而得xm之值，从而计算维恩常数。取温度为参数向量，取波长为自变量向量，形成温度和波长矩阵，根据普朗克黑体辐射公式可直接计算单色辐射本领，画出以波长为自变量的曲线族。每条曲线的峰值坐标可用最大值指令获得，至于峰值曲线则根据维恩定律计算。程序Plank1.m如下。%温度不同的普朗克黑体单色辐射能力与波长的曲线clear %清除变量k=1.38054e-23; %玻尔兹曼常数h=6.626e-34; %普朗克常数c=2.

5、997925e8; %光速sigma=5.6688e-008; %斯特潘常数b=0.0029; %维恩常数t=1400:100:2000; %热力学温度向量n=length(t); %向量长度lambda=0:0.01:5*1e-6; %波长向量lambda(1)=eps; %给零加一小量使分母不为零T,L=meshgrid(t,lambda); %波长和温度矩阵M=2*pi*h*c2./(exp(h*c./(k*T.*L)-1)./L.5;%单色辐射能力figure %创建图形窗口plot(lambda*1e6,M) %画曲线hl=legend(repmat('itTrm='

6、,n,1),num2str(t'),repmat('K',n,1);%标记图例fs=16; %字体大小set(hl,'fontsize',fs) %设置图例大小grid on %加网格title('普朗克黑体单色辐射能力与波长的关系','fontsize',fs)%标题xlabel('波长itlambdarm/mum','fontsize',fs)%横坐标yl='单色辐射能力itMrm(itlambdarm,itTrm)/(Wcdotm-3)'%纵坐标字符串ylabel(yl

7、,'fontsize',fs) %纵坐标txt='itbrm=' num2str(b) 'mcdotK' %维恩常数文本txt=txt ',itsigmarm=' num2str(sigma) 'W/(m2cdotK4)'%斯特潘常数文本text(0,max(M(:)/10,txt,'fontsize',fs) %显示常数hold on %保持图像mx,ix=max(M); %找最大值和下标stem(lambda(ix)*1e6,mx,'-','filled') %画

8、直杆图text(lambda(ix)*1e6,mx,num2str(lambda(ix)'*1e6),'fontsize',fs)%显示峰值波长t=1300:2020; %较密的温度向量lambda=b./t; %波长向量m=2*pi*h*c2./(exp(h*c./(k*t.*lambda)-1)./lambda.5;%单色辐射能力向量plot(lambda*1e6,m) %画峰值曲线图示(1)取温度为参数，黑体的单色辐射本领与波长的关系如P图所示。不论温度是多少，单色辐射本领随波长的增加先增加再减小，峰值波长与温度的关系遵守维恩位移律：峰值波长与温度成反比。温度升高

9、时，峰值波长变短，峰变高。(2)曲线下的面积表示总辐射本领，温度越高，曲线下的面积越大，总辐射本领越强。P图 程序Plank2.m如下。%普朗克黑体单色辐射斯特潘常数和维恩常数求法clear %清除变量k=1.38054e-23; %玻尔兹曼常数h=6.626e-34; %普朗克常数c=2.997925e8; %光速syms x %定义符号变量y=x3/(exp(x)-1); %被积函数i=int(y,0,inf) %求积分sigma=eval(2*pi*k4/h3/c2*i) %求斯特潘常数y=x5/(exp(x)-1); %普朗克约化公式d=diff(y) %求符号导数s=solve(d) %求符号零解e=eval(s) %求零解的数值b=h*c/e/k %求维恩常数结果i =1/15*pi4sigma = 5.6688e-008d =5*x4/(exp(x)-1)-x5/(