matlab普朗克黑体辐射公式,MATLAB 黑体辐射规律的研究详解.doc

MATLAB 黑体辐射规律的研究详解

黑体辐射规律的研究

湖南大学 XX院系 XX 专业 XX年级

姓名学号

[问题] 黑体辐射的规律

在任何温度下对任意波长的电磁波只吸收不反射的物体称为绝对黑体，简称黑体。根据实验得出两个实验规律。黑体的总辐射本领(能力)为

E(T) = σT4，

这就是斯特潘-玻尔兹曼定律，其中，σ = 5.67×10-8W/(m2·K4)，σ称为斯特潘常数。黑体的单色辐射本领(能力)的峰值波长与温度的关系为

Tλm = b，

这就是维恩位移定律，其中，b = 2.897×10-3m·K，b称为维恩常数。

根据普朗克提出的黑体辐射公式，计算斯特潘常数和维恩常数。以温度为参数，单色辐射本领与波长的曲线有什么特点，峰值与波长的关系曲线有什么特点？

[数学模型]

黑体的单色辐射本领是在单位时间内从物体表面单位面积上所发射的波长在λ到λ + 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.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)

这就是维恩常数。理论值与实验值符号得很好。

[算法]用MATLAB的符号积分计算积分可直接计算I的值，从而计算斯特潘常数。设

， (14.1.5*)

用MATLAB的符号导数先求y的导数，求导数的零解而得xm之值，从而计算维恩常数。

取温度为参数向量，取波长为自变量向量，形成温度和波长矩阵，根据普朗克黑体辐射公式可直接计算单色辐射本领，画出以波长为自变量的曲线族。每条曲线的峰值坐标可用最大值指令获得，至于峰值曲线则根据维恩定律计算。

[程序]P.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=length(t); %向量长度

lambda=[0:0.01:5]*1e-6; %波长向量

lambda(1)=eps; %给零加一小量使分母不为零

[T,L]=meshgrid(t,lambda); %波长和温度矩阵

M=2*pi*h*c^2./(exp(h*c./(k*T.*L))-1)./L.^5;%单色辐射能力

figure

matlab普朗克黑体辐射公式,MATLAB 黑体辐射规律的研究详解.doc相关推荐

matlab普朗克黑体辐射公式,基于MATLAB的黑体辐射量计算
基于MATLAB的黑体辐射量计算介绍了黑体辐射理论,利用MATLAB求解了普朗克黑体辐射公式在有界域λ1-λ2内的积分,通过实例对梯形法和辛普森法两种数值积分方法进行了分析和比 (本文共3页) 阅读 ...
怎么用matlab画黑体辐射图像,matlab_黑体辐射规律的研究.doc
matlab_黑体辐射规律的研究.doc 1黑体辐射规律的研究湖南大学XX院系XX专业XX年级姓名学号问题黑体辐射的规律在任何温度下对任意波长的电磁波只吸收不反射的物体称为绝对黑体,简称黑体.根据实验 ...
普朗克黑体辐射公式使用指南
普朗克黑体辐射公式避坑(实用,告别疑惑): 最近想要使用普朗克公式反函数计算星上亮温.看了很多论文.教材和课件,发现一些常量的量纲甚至值都存在差异,找了半天也没个标准的,索性全试了一遍,结果就是一个都 ...
【张朝阳的物理课笔记】10. 由普朗克黑体辐射公式导出维恩定理，计算太阳表面的温度
10.1 太阳表面的温度上节课得到普朗克黑体辐射公式: I(ω)=ℏπ2c2ω3eℏω/KT−1I(\omega)=\frac\hbar{\pi^2c^2}\frac{\omega^3}{e^{\h ...
matlab车牌匹配时读取,基于Matlab的车牌识别(完整版)详解.doc
基于Matlab的车牌识别(完整版)详解.doc 基于Matlab的车牌识别摘要:车牌识别技术是智能交通系统的重要组成部分,在近年来得到了很大的发展.本文从预处理.边缘检测.车牌定位.字符分割.字符 ...
单相全控桥有源逆变电路matlab仿真,单相桥式全控整流及有源逆变电路的MATLAB仿真详解.doc...
单相桥式全控整流及有源逆变电路的MATLAB仿真详解.doc 学号天津城建大学控制系统仿真大作业单相整流-逆变电路仿真模型学生姓名班级成绩学院年月日一.仿真电路原理图:2 二.单相 ...
Matlab中画图以及plot函数及legend函数详解
Matlab中plot函数及legend函数详解 Matlab中plot函数全功能解析Matlab中plot函数及legend函数详解功能二维曲线绘图语法 plot(Y) plot(X1,Y1, ...
(Matlab源码)Matlab实现算术编码（Arithmetic coding）超级详解（每一段代码都可以看懂）
1.代码功能输入:一个字符串输出: codeword(码值) codeword所占的位数 2.代码框图 3.代码超详解统计字符串中的字符种类,调用函数,放入数组b中. b=unique(str1 ...
matlab自带SVM算法例子（附函数详解）
一.程序 SVM理论的学习,见文章:http://blog.csdn.net/ckzhb/article/details/68941695 load fisheriris ...
python gpu加速库比matlab快吗_Python之Unittest和Requests库详解
1.按类来执行 import unittest class f1(unittest.TestCase): def setUp(self): pass def tearDown(self): pass ...

matlab普朗克黑体辐射公式,MATLAB 黑体辐射规律的研究详解.doc

matlab普朗克黑体辐射公式,MATLAB 黑体辐射规律的研究详解.doc相关推荐

最新文章

热门文章