【张朝阳的物理课笔记】9. 瑞丽金斯公式的推导(下),普朗克修正,黑体辐射公式

本文深入探讨了瑞丽金斯公式及其在高频区遇到的问题,即紫外灾难。通过普朗克的能级量子化理论,介绍了修正后的黑体辐射公式,解释了如何解决紫外灾难,并详细阐述了普朗克黑体辐射公式的推导过程。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

9.1 瑞丽金斯公式,紫外灾难

空腔内充斥着各种频率的波,每个谐振子都受到所有波的作用。
谐振子在频率为ω\omegaω的电场中受到的电场强度: E=E0eiωtE=E_0e^{i\omega t}E=E0et
用第8课里的阻尼运动方程再加上电场:
d2xdt2=−ω02x−γdxdt+qE0eiωtm\frac {d^2x}{dt^2}=-\omega_0^2x-\gamma\frac{dx}{dt}+\frac{qE_0e^{i\omega t}}mdt2d2x=ω02xγdtdx+mqE0et
空腔内能量守恒,所以上式第二项辐射能量和第三项吸收能量是相等的。
猜测运动频率是ω\omegaω,令x=x^0eiωtx=\hat x_0e^{i\omega t}x=x^0et,代入方程得:
x^0(ω02−ω2+iγω)=qE0m\hat x_0(\omega_0^2-\omega^2+i\gamma\omega)=\frac{qE_0}mx^0(ω02ω2+iγω)=mqE0
x^0=qE0m1ω02−ω2+iγω\hat x_0=\frac{qE_0}m\frac 1{\omega_0^2-\omega^2+i\gamma\omega}x^0=mqE0ω02ω2+iγω1
x^0\hat x_0x^0是复数(形式是:x0eiΔx_0 e^{i\Delta}x0eiΔ),电场是振动的,电荷位移被它拖着振动,所以会晚一个相位
x^02=x^x^∗=q2E02m21(ω02−ω2)2+γ2ω2\hat x_0^2=\hat x\hat x^*=\frac{q^2E_0^2}{m^2}\frac 1{(\omega_0^2-\omega^2)^2+\gamma^2\omega^2}x^02=x^x^=m2q2E02(ω02ω2)2+γ2ω21

继续分析:
ax=d2xdt2=x0ω2a_x=\frac{d^2x}{dt^2}=x_0\omega^2ax=dt2d2x=x0ω2
⟨ax2⟩=12x02ω2\langle a_x^2\rangle=\frac{1}2x_0^2\omega^2ax2=21x02ω2
代入:⟨P⟩=q26πϵ0c3⟨ax2⟩\langle P\rangle= \frac{q^2}{6\pi\epsilon_0c^3}\langle a_x^2\rangleP=6πϵ0c3

在MATLAB中,绘制比较普朗克黑体辐射公式、瑞利-金斯修正公式以及维恩定律的图形,通常涉及以下几个步骤: 1. **导入所需函数库**:开始时,你需要加载`syms`函数库以处理符号数学,这对于公式计算至关重要。 ```matlab syms T wavelength J_lambda ``` 2. **定义公式**: - **普朗克公式**:\(J_{\lambda} = \frac{2hc^2}{\lambda^5}\frac{1}{e^\frac{hc}{\lambda kT}-1}\) - **瑞利-金斯修正公式**:考虑到实验观测到的在短波长处辐射强度略低于普朗克公式预测的现象,可以添加一个高阶项,如\(J_{\lambda}^{RK} = J_{\lambda} + A\frac{\lambda^4}{(\lambda^5+B)^2}\) (A, B为常数) - **维恩公式**:描述了低温下黑体辐射峰值位置,\(B_T = \frac{2897.8}{T^4}\) 3. **设定变量范围**:选择合适的温度范围 \(T\) 和波长范围 \(\lambda\) 进行绘制。 ```matlab T = linspace(500, 5000, 100); % 温度范围 wavelength = linspace(0.1, 10, 100); % 波长范围 ``` 4. **计算辐射强度**:分别计算每个公式辐射强度。 5. **绘制比较图**: - 使用`plot`函数绘制三个函数,并在图上标明各自代表的公式。 ```matlab plot(wavelength, planksRadiation(T, wavelength), 'r', 'LineWidth', 2, 'DisplayName', '普朗克公式') hold on; plot(wavelength, rydbergsCorrection(T, wavelength, A, B), 'b', 'LineWidth', 2, 'DisplayName', '瑞利-金斯修正') plot(wavelength, viennesLaw(B_T), 'g', 'LineWidth', 2, 'DisplayName', '维恩定律') ``` 6. **设置轴标签和图例**: ```matlab xlabel('波长 (\(\mu m\)') ylabel('辐射强度 (W/m^2/\(\mu m\)') legend boxoff title('比较普朗克、瑞利-金斯及维恩辐射公式') ``` 7. **显示图形**: ```matlab hold off ```
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值