李寅杰 徐 慧

（南京師范大學(xué)附屬中學(xué) 江蘇 南京 210003）

1 引言

在物理教學(xué)中，在講解一些具有復(fù)雜表達(dá)式的物理問(wèn)題時(shí)，常受限制于學(xué)生當(dāng)前的知識(shí)水平，只能一帶而過(guò)，導(dǎo)致學(xué)生對(duì)知識(shí)死記硬背，理解淺顯.Mathematica是Wolfram公司開(kāi)發(fā)的科學(xué)計(jì)算軟件，具有強(qiáng)大的數(shù)值計(jì)算能力［1］.利用 Mathematica軟件可解決數(shù)學(xué)計(jì)算的困難，還可以得到直觀的圖像，讓很多復(fù)雜問(wèn)題清晰直觀，使學(xué)生輕松、深入地理解物理概念、公式.本文以黑體輻射中的普朗克公式為例，通過(guò)數(shù)值擬合和符號(hào)運(yùn)算，驗(yàn)證在低頻率范圍與瑞利-金斯公式的符合度問(wèn)題.

2 正文

2.1 概述

瑞利-金斯公式是由英國(guó)物理學(xué)家瑞利和金斯提出，用于描述黑體輻射問(wèn)題光譜能量分布，表達(dá)式為

其中ρ為熱輻射能量密度，ν為熱輻射頻率，T為熱力學(xué)溫度，c為真空中光速，h為普朗克常量，k為玻爾茲曼常數(shù).該公式在長(zhǎng)波區(qū)域與實(shí)驗(yàn)符合較好，但在短波區(qū)域，會(huì)導(dǎo)致能量無(wú)限大的荒謬結(jié)果，也就是當(dāng)時(shí)著名的“紫外災(zāi)難”.

普朗克公式是德國(guó)物理學(xué)家普朗克提出，用于描述黑體輻射問(wèn)題光譜能量分布的公式.表達(dá)式為

其中ρ為熱輻射能量密度，ν為熱輻射頻率，T為熱力學(xué)溫度，c為真空中光速，h為普朗克常量，k為玻爾茲曼常數(shù).這一公式彌補(bǔ)了之前瑞利-金斯公式和維恩位移定律的不足，與實(shí)驗(yàn)吻合度高.

我們利用Mathematica軟件進(jìn)行模擬時(shí)，為了使數(shù)據(jù)直觀且便于輸入，取能量的單位為10－20J，頻率ν的單位為1013Hz，熱力學(xué)溫度T的單位為100 K，普朗克常量h取6.626 069 572 9×10－34J·s［2］，光速c取299 792 458m·s－1［3］，玻爾茲曼常數(shù)k取1.380 648 813×10－23J·K－1［4］.代入數(shù)據(jù)后，普朗克公式的表達(dá)式可寫(xiě)為

2.2 數(shù)值模擬

將式（1）輸入，證明普朗克公式在低頻率區(qū)域轉(zhuǎn)化為瑞利-金斯公式，按照普朗克公式在低頻率區(qū)域取值后，可以擬合為瑞利-金斯公式.首先考慮溫度為100K的情況，取溫度為1，從ν＝0.000 1開(kāi)始，每次間隔0.000 1，一直到0.002，取點(diǎn).得到

（0.000 1，2.049 18×10－8）

（0.000 2，8.194 77×10－8）

（0.000 3，1.843 38×10－7）

（0.000 4，3.276 33×10－7）

（0.000 5，5.118 04×10－7）

（0.000 6，7.368 21×10－7）

（0.000 7，1.002 65×10－6）

（0.000 8，1.309 28×10－6）

（0.000 9，1.656 65×10－6）

（0.001 0，2.044 76×10－6）

（0.001 1，2.473 57×10－6）

（0.001 2，2.943 04×10－6）

（0.001 3，3.453 16×10－6）

（0.001 4，4.003 88×10－6）

（0.001 5，4.595 19×10－6）

（0.001 6，5.227 05×10－6）

（0.001 7，5.899 43×10－6）

（0.001 8，6.612 30×10－6）

（0.001 9，7.365 64×10－6）

（0.002 0，8.159 42×10－6）

由于瑞利-金斯公式的形式為ρ～Tν2，所以我們對(duì)這些數(shù)據(jù)用二次函數(shù)擬合.將散點(diǎn)和擬合出的函數(shù)繪制出，并顯示在同一坐標(biāo)系中，可以清晰地看出，完美符合二次曲線，如圖1所示.（在這里Mathematica為了使圖像美觀，將y軸右移，實(shí)質(zhì)上該圖像ν≥0.）

圖1

我們依次考慮溫度為300K，500K，900K和1 000K的情況，代碼與之前幾乎相同，擬合的結(jié)果為：

300K：R＝－2.610 7×10－9＋

1.332 26 ×10－5ν＋6.133 54ν2

500K：R＝－2.611 68×10－9＋

1.332 71 ×10－5ν＋10.232 9ν2

700K：R＝－2.612 1×10－9＋

1.332 9 ×10－5ν＋14.332 2ν2

900K：R＝－2.612 33×10－9＋

1.333 01 ×10－5ν＋18.431 6ν2

1 000 K：R＝－2.612 41×10－9＋

1.333 04 ×10－5ν＋20.481 3ν2

我們著手進(jìn)行兩公式的轉(zhuǎn)化，觀察不同溫度下R函數(shù)，一次項(xiàng)系數(shù)、常數(shù)項(xiàng)的數(shù)量級(jí)與二次項(xiàng)系數(shù)的數(shù)量級(jí)相差巨大，而且瑞利-金斯公式中只含二次項(xiàng)，故舍去常數(shù)項(xiàng)和一次項(xiàng).現(xiàn)在R函數(shù)與瑞利－金斯公式在形式上相一致，只需驗(yàn)證二次項(xiàng)系數(shù).由于二次項(xiàng)系數(shù)與T成正比，所以對(duì)R函數(shù)的系數(shù)與溫度（表1）進(jìn)行線性擬合.

表1 函數(shù)的系數(shù)與溫度值

由于一次項(xiàng)系數(shù)遠(yuǎn)大于常數(shù)項(xiàng)，故舍去常數(shù)項(xiàng)，Rho＝2.049 67t.這與現(xiàn)在單位制下的瑞利-金斯公式：ρ＝2.049 66Tν2，幾乎完全相同，由此可知，普朗克公式在長(zhǎng)波區(qū)域可轉(zhuǎn)化為瑞利-金斯公式.

而當(dāng)我們將散點(diǎn)范圍增大時(shí)，可見(jiàn)到在高頻區(qū)域，二者出現(xiàn)較大偏差.取溫度為1，ν＝0.000 1開(kāi)始，每次間隔0.01，一直到0.2，取點(diǎn).將瑞利-金斯公式的圖像和普朗克公式的散點(diǎn)繪出，并顯示在同一坐標(biāo)系中，得到圖2.清晰地看出，在短波區(qū)域，二者的差別不斷增大.

圖2

3 反思

觀察數(shù)值模擬的數(shù)據(jù)，在第6位有效數(shù)字處與公式不同，分析原因及改進(jìn)措施：

（1）Mathematica使用過(guò)程中分布計(jì)算，導(dǎo)致第一步保留的有效數(shù)字不足.應(yīng)當(dāng)一次性鍵入算式，只保留一次有效數(shù)字.

經(jīng)過(guò)驗(yàn)證，由于Mathematica本身可以儲(chǔ)存較多位數(shù)，所以對(duì)于此處的5位小數(shù)結(jié)果，是否一步計(jì)算并無(wú)差別.

（2）在二次擬合時(shí)選用了ax2＋bx＋c的形式，可以嘗試使用ax2，ax2＋bx，ax2＋c等形式，看是否可以減少誤差.在線性擬合時(shí)選用了kx＋b的形式，可以嘗試使用kx形式，看是否可以減少誤差.

經(jīng)過(guò)驗(yàn)證，當(dāng)選擇ax2和kx時(shí)，得到的結(jié)果是2.049 66T，確實(shí)可以減少誤差.我們從中可以受到啟發(fā)，既然是驗(yàn)證性質(zhì)的，就可以根據(jù)已知進(jìn)行盡量精確的擬合，這樣才能提高結(jié)果的準(zhǔn)確性.

4 總結(jié)

Mathematica具有清晰、直觀的特點(diǎn)，可較好地演示.同時(shí)其擁有強(qiáng)大的計(jì)算、繪圖功能，可以快速、準(zhǔn)確地完成計(jì)算和推導(dǎo).對(duì)于一些難以計(jì)算、推導(dǎo)過(guò)程的知識(shí)點(diǎn)，Mathematica的使用可能會(huì)帶來(lái)意想不到的效果.在科學(xué)飛速發(fā)展的今天，Mathematica為物理教學(xué)開(kāi)辟了一條新思路.

1 https：∥zh.m.wikipedia.org／wiki／Mathematica

2 http：∥ physics.nist.gov ［Thursday，02－Jun－2011 21：00：12EDT］. National Institute of Standards and Technology，Gaithersburg，MD 20899

3 Penrose，R （2004）.The Road to Reality：A Complete Guide to the Laws of the Universe.Vintage Books.pp.410－1.ISBN 978－0－679－77631－4

4 http：∥physics.nist.gov／constants ［Thursday，02－Jun－2011 21：00：12EDT］. National Institute of Standards and Technology，Gaithersburg，MD 20899