实验8(综合性的设计性实验)--安铎.docVIP

综合性设计性实验 班级：应物101 姓名：安铎 学号：10110020101 实验8 平衡辐射场的内能密度随频率和波长的变化 实验目的： 1、弄清平衡辐射场内能密度随频率和波长变化规律 2、推导辐射场能量取极大时对应的波长、频率满足的方程 3、编程求解上述方程 4、编程画出平衡辐射场内能密度随频率和波长变化曲线，理解辐射场能量取极大时对应的波长和频率间的差异 实验背景： 黑体辐射是量子力学重要的实验基础之一，利用光子假说和玻色统计可以给出黑体辐射的内能密度与光子频率或波长之间的关系。 光子是玻色子，自旋为1；光子静止质量为0；光子的动量、能量、频率、波矢、速度之间的关系为 光子数不守恒，，化学势为0。这样光子气体的统计分布为 光子的自旋在传播方向的投影有两个值，对应于左右两个偏振光。 在体积V内，动量为的量子态数为 在内的光子数为 乘以每一个光子的能量，则辐射场的内能为 (1) 若取V=1，则给出了辐射场的内能密度。 注意到光的频率与波长间的关系，容易给出内能密度与波长的关系 (2) 若令，则使(1)式有极大值的x满足 (3) 若令，则使(2)式有极大值的x满足 (4) 由(3)、(4)两式决定的波长和频率并未按 的关系对应。可以用程序计算内能密度随波长、频率的变化关系，用作图软件（例如Origin）严格画出这些关系，看两者是否对应并理解其原因。 实验报告： 对原理中给出的内能密度随波长、频率变化关系公式进行变量代换： 内能密度与频率关系： 令得内能密度： 提出无关项系数令：=a=200 则得出简化式： 内能密度与波长关系： 令 得内能密度： 提出无关项系数令：=b=15 则得出化简式： 编程求解以上方程 （1）内能密度随频率、波长变化关系初步求解 编写程序如下： program anduo implicit none real(8) i,x,y,z real,parameter:: a=200,b=15 open(10,file=anduo1.txt) open(20,file=anduo2.txt) do i=1,100 x=i*0.1 y=a*(x**3.0)/(exp(x)-1) write(10,*) x,y z=b*(x**5.0)/(exp(x)-1) write(20,*) x,z enddo write(*,*) 应物101班，安铎，10110020101 close(10) close(20) End 用Origin画出图像如下： 在内能密度与频率间的关系中： 若令，则使(1)式有极大值的x满足 用牛顿法求x的值，程序如下： PROGRAM anduo dl=1.0E-06 a=4.0 b=5.0 dx=b-a x0=(a+b)/2.0 istep=0 DO WHILE (ABS(dx).GT.dl) x1=x0-f(x0)/df(x0) dx=x1-x0 x0=x1 istep=istep+1 ENDDO WRITE(*,999) istep,x0,dx 999 FORMAT(I4,2F16.8) END SUBROUTINE secant(dl,x0,dx,istep) istep=0 x1=x0+dx DO WHILE (ABS(dx).GT.dl) d=f(x1)-f(x0) x2=x1-f(x1)*(x1-x0)/d x0=x1 x1=x2 dx=x1-x0 istep=istep+1 ENDDO END FUNCTION f(x) f=3*exp(-x)+x-3 write(*,*)应物101班，安铎，10110020101 END FUNCTION df(x) df=1-3*exp(-x) END 求出结果： 在内能密度与波长的关系中： 若令，则使该式有极大值的x满足 用牛顿法求x的值，编程如下： PROGRAM anduo2 dl=1.0E-06 a=4.0 b=5.0 dx=b-a x0=(a+b)/2.0 istep