MATLAB黑体辐射规律的研究.docxVIP

PAGE PAGE # 4 4 0 3 4 4 0 3 PAGE PAGE # 黑体辐射规律的研究 湖南大学 XX院系XX专业XX年级 姓名学号 姓名 ［问题］黑体辐射的规律 在任何温度下对任意波长的电磁波只吸收不反射的物体称为绝对黑体， 简称黑体。根据 实验得出两个实验规律。黑体的总辐射本领 (能力)为 4 E(T)=打, 这就是斯特潘-玻尔兹曼定律，其中， 卢5.67 W8W/(m2 K4), 称为斯特潘常数。黑体的单 色辐射本领(能力)的峰值波长与温度的关系为 T Am = b, 这就是维恩位移定律，其中， b = 2.897 W3m-K , b称为维恩常数。 根据普朗克提出的黑体辐射公式， 计算斯特潘常数和维恩常数。 以温度为参数，单色辐 射本领与波长的曲线有什么特点，峰值与波长的关系曲线有什么特点？ ［数学模型］ 黑体的单色辐射本领是在单位时间内从物体表面单位面积上所发射的波长在 A到入+ d入 范围内的辐射能量 dE(入T)与波长间隔d入之比 dE (■ ,T) M「，T) , (14.1.1) d扎 M(入T)表示在单位时间内从物体表面单位面积发射的波长在附近单位波长间隔内的辐射本 领，是波长和温度的函数，其单位是 W/m3。普朗克提出的黑体单色辐射本领的公式为 (14.1.2)2 nic2 (14.1.2) ■5[exp(乎)-1]kT扎 其中，k为玻尔兹曼常数，h为普朗克常数，c为真空中的光速。 对波长从0到无穷大积分就得总辐射本领，即：黑体单位面积辐射能量的功率 E(T)二 M C ,T)d ?, E(T)二 M C ,T)d ?, 0 设 (14.1.2)式可化为 he kT - (14.1.3) (14.1.4) M(x,T) 3 2 x h c (e -1) (14.1.5) he 由(14.1.4)式得dx 2d,所以(14.1.3)式可化为 kT九 (14.1.6)2 n4其中，-疋为常数，］为积分3x—dx ° °e -1(14.1.7)手工计算积分 (14.1.6) 2 n4 其中，-疋为常数，］为积分 3 x —dx ° °e -1 (14.1.7) 手工计算积分I比较麻烦，其步骤如下 3 . x X e xdx 01 - e QO 二 X 0 oO oO 八e』xdx八 x n =0 n 三1 0 oO e n卫 设y = nx，可得 DO . °O八 A .y3n 3 n 0 e_y dy 二 =:1 ~4 n ：1 n 3! 4 n —。 15 由此可得 (14.1.8)这就是斯特潘常数。当波长趋于零时,CI = 5.6688 1 氐8 (14.1.8) 这就是斯特潘常数。 当波长趋于零时, 单色辐射本领趋于零；当波长趋于无穷大时，单色辐射本领也趋于零。 因此单色辐射本领随波长的变化有极值。令 dM (x,T)/dx = 0 ,可得方程： _x x = 5(1 - e ), (14.1.9) 一般用迭代算法计算上式之值， 除了零解之外，可得Xm的值为4.965,由(14.1.4)式可得维恩 常数 hckxm= hc kxm =0.0029。 (14.1.10) 这就是维恩常数。理论值与实验值符号得很好。 ［算法］用MATLAB的符号积分计算积分可直接计算 I的值，从而计算斯特潘常数。设 5 (14.1.5*)x (14.1.5*) y =^r~ e -1 用MATLAB的符号导数先求y的导数，求导数的零解而得 Xm之值，从而计算维恩常数。 取温度为参数向量，取波长为自变量向量，形成温度和波长矩阵，根据普朗克黑体辐 射公式可直接计算单色辐射本领， 画出以波长为自变量的曲线族。 每条曲线的峰值坐标可用 最大值指令获得，至于峰值曲线则根据维恩定律计算。 ［程序］Plank1.m如下。 %温度不同的普朗克黑体单色辐射能力与波长的曲线 clear %清除变量 k=1.38054e_23; %玻尔兹曼常数 h=6.626e_34; %普朗克常数 C=2.997925e8; %光速 sigma=5.6688e_008; %斯特潘常数 b=0.0029;t=1400:100:2000;n=len gth(t);lambda=[0:0.01:5]*1e-6;lambda(1)=eps;[T,L]=meshgrid(t,lambda);M=2*pi*h*cA2./(exp(h*c./(k*T.*L))-1)./L.A5;figureplot(lambda*1e6,M)%维恩常数%热力学温度向量%向量长度%波长向量%给零加一小量使分母不为零% b=0.0029; t=1400:100:2000; n=len gth(t); lambda=[0:0.01:5]*1e-6; lam