劉兆龍 胡海云

(北京理工大學(xué)物理學(xué)院，北京 100081)

在大學(xué)物理課程中嵌入數(shù)值計(jì)算的捷徑

劉兆龍 胡海云

(北京理工大學(xué)物理學(xué)院，北京 100081)

數(shù)值計(jì)算已經(jīng)成為研究復(fù)雜系統(tǒng)物理規(guī)律時(shí)不可缺少的手段，但基于各個(gè)方面的困難，它很難在大學(xué)物理課程中有所體現(xiàn)。本文介紹了借助Excel軟件將數(shù)值計(jì)算嵌入大學(xué)物理課程的方法，給出了利用Excel軟件實(shí)現(xiàn)可視化、數(shù)值積分和求微分方程數(shù)值解的方法。文中所述的方法簡捷易行，可操作性強(qiáng)，且不必耗費(fèi)大量學(xué)時(shí)，非常適用于面向低年級(jí)理工科學(xué)生的大學(xué)物理課程。采用此方法，可以使低年級(jí)學(xué)生較早地熟悉數(shù)值計(jì)算，更好地掌握物理知識(shí)和應(yīng)用物理知識(shí)解決問題，是大學(xué)物理課程教學(xué)改革的一種有益嘗試。

大學(xué)物理；數(shù)值計(jì)算；教學(xué)改革；Excel

1 背景

伴隨著計(jì)算機(jī)技術(shù)的不斷進(jìn)步，數(shù)值計(jì)算已經(jīng)被公認(rèn)為是研究復(fù)雜系統(tǒng)物理規(guī)律時(shí)不可缺少的手段，也成為了高校物理教學(xué)的重要內(nèi)容之一。從20世紀(jì)80年代起，國內(nèi)外大學(xué)的物理系陸續(xù)開設(shè)了計(jì)算物理學(xué)課程，教授物理專業(yè)的學(xué)生利用數(shù)值計(jì)算方法去解決問題。說到物理教育與教學(xué)，受眾最廣的群體是非物理專業(yè)的學(xué)生，特別是理工科學(xué)生，其群體數(shù)量遠(yuǎn)遠(yuǎn)大于物理專業(yè)的學(xué)生。在面向該群體的大學(xué)物理課程中，數(shù)值計(jì)算一直沒有像近代物理那樣得到大規(guī)模的普及。如何將計(jì)算物理學(xué)普物化，即將計(jì)算物理學(xué)的研究方法引入面向非物理專業(yè)的理工科大學(xué)物理課程甚至是文科物理課程中，將之與原有教學(xué)內(nèi)容有機(jī)結(jié)合，加強(qiáng)大學(xué)物理課程的應(yīng)用性，使大學(xué)物理課程更現(xiàn)代化，是一個(gè)值得研究的課題，也是大學(xué)物理課程教學(xué)的改革方向之一。

現(xiàn)今，隨著各校教學(xué)設(shè)施的改善和手提電腦的日益普及，實(shí)施計(jì)算物理普物化的時(shí)機(jī)日益成熟，其主要困難在于計(jì)算軟件的使用。大學(xué)物理課程面向低年級(jí)學(xué)生，學(xué)生們的編程能力和計(jì)算機(jī)語言基礎(chǔ)普遍有限，而教師不能在物理課上花費(fèi)大量課時(shí)講解編程，這一矛盾是計(jì)算物理普物化的困難之一。此外，對(duì)大部分學(xué)生來說，大學(xué)物理學(xué)課程本身難度相對(duì)較大，如果再花費(fèi)很多時(shí)間用于編程，會(huì)增加學(xué)生的學(xué)習(xí)負(fù)擔(dān)或減少用于物理知識(shí)本身的學(xué)習(xí)時(shí)間，有喧賓奪主之憂。

2 Excel及其在大學(xué)物理課程數(shù)值計(jì)算中的應(yīng)用

很顯然要在大學(xué)物理中引入數(shù)值計(jì)算，需要借助相對(duì)易學(xué)的通用軟件，在現(xiàn)有的教學(xué)框架中有機(jī)地融入計(jì)算物理學(xué)的思想與方法，凸顯數(shù)值計(jì)算的魅力，使低年級(jí)學(xué)生較快地解決一些“相對(duì)復(fù)雜”的問題。國外大學(xué)物理課程中已經(jīng)在此方向進(jìn)行了有益的嘗試[1]。目前可用于科學(xué)計(jì)算的軟件較多，例如Matlab, Mathematic, Maple, Python計(jì)算平臺(tái)等，但是使用這些軟件需要學(xué)生們學(xué)習(xí)程序設(shè)計(jì)語言，并且學(xué)生們?cè)诟咧袝r(shí)很少接觸它們。

Microsoft Excel是微軟公司Microsoft office的組件之一，是一款功能強(qiáng)大的電子表格處理軟件，包含大量公式和函數(shù)，可以執(zhí)行計(jì)算。許多教師對(duì)于此軟件并不陌生，常用之統(tǒng)計(jì)學(xué)生成績或在大學(xué)物理實(shí)驗(yàn)課上處理實(shí)驗(yàn)數(shù)據(jù)，而且有些學(xué)生在高中的信息技術(shù)課上接觸過Excel。將此軟件用于低年級(jí)本科生的大學(xué)物理課程之中，可以起到事半功倍的作用，使得我們可以基于傳統(tǒng)大學(xué)物理課程的知識(shí)框架，回避對(duì)學(xué)生進(jìn)行系統(tǒng)地編程訓(xùn)練，幫助學(xué)生了解建立模型、設(shè)計(jì)計(jì)算方法、編程和結(jié)果分析這一計(jì)算物理學(xué)研究問題的過程，并進(jìn)行簡單應(yīng)用。它還有助于學(xué)生對(duì)物理現(xiàn)象和規(guī)律的理解。

下面通過大學(xué)物理課程中幾個(gè)常見的實(shí)例，看如何使用Excel實(shí)現(xiàn)數(shù)值計(jì)算。

2.1 在大學(xué)物理課程中實(shí)現(xiàn)“可視化”

使用Excel軟件，可以輕松地解決這個(gè)問題，即使對(duì)于沒有學(xué)過編程或?qū)ｉT的計(jì)算機(jī)語言的學(xué)生，也能快速掌握曲線的作圖方法，獨(dú)立繪出漂亮的函數(shù)圖線，并利用它探究其中的物理知識(shí)。

例如要繪制氧氣在0℃時(shí)的麥克斯韋速率分布曲線，完成步驟為

圖1 利用Excel完成的麥克斯韋速率分布曲線(截取了局部工作表)。繪制其他氣體曲線時(shí)，只要在C2單元格中輸入相應(yīng)的質(zhì)量。在C3單元中輸入不同的數(shù)值，可以得到不同溫度下的分布曲線。為了便于比較，圖中給出了溫度為T和T/2的兩條曲線

? 新建Excel工作表，在第一行中敲入注釋，解釋此工作的內(nèi)容，如：麥克斯韋速率分布函數(shù)，如圖1。在A2單元格中輸入文字“輸入質(zhì)量：”，C2單元中敲入氧氣分子的質(zhì)量5.32E-26，并在名稱框中將該單元的名稱定義為m，隨后在D2單元格中輸入其單位kg。依此法在第3、4和5 行中分別輸入溫度、溫度的步長值和玻爾茲曼常數(shù)，名稱分別定義為T、DELTA和k。若要研究其他氣體，只要直接在C2單元格中敲入其分子或原子質(zhì)量的數(shù)值即可。要改變溫度，也只需在C3單元格中直接敲入所需數(shù)值。同理，改變C4的數(shù)值就可以改變溫度的步長

? 建立表頭，如速率、速率分布函數(shù)等，見圖1中的第7行

圖2 利用Excel完成的拍頻繪圖(截取了局部工作表)

? 在A列中輸入速率。例如若速率范圍是0～1200m/s，步長為25m/s，則在A8單元格中輸入0，A9單元格中輸入=A8+DELTA，將鼠標(biāo)放在A9的右下角，當(dāng)鼠標(biāo)成黑十字星時(shí)，按住鼠標(biāo)的左鍵豎直向下拖拽(本文中將此操作簡稱為豎直向下拖拽鼠標(biāo))，直至速率到達(dá)1200m/s，即設(shè)置好了速率的取值

? 在B8單元中輸入麥克斯韋速率分布函數(shù)=4*3.14*(m/(2*3.14*k*T))^(3/2)*A8^2*EXP(-m*A8^2/(2*k*T))

? 從B8單元格的右下角豎直向下拖拽鼠標(biāo)，完成第B列各單元格的計(jì)算

? 選中A、B列中所有數(shù)據(jù)，點(diǎn)擊工具欄中的繪圖按鈕，選擇繪制XY散點(diǎn)圖，得到曲線

? 編輯曲線，加上說明文字，完成曲線繪制，見圖1

利用這個(gè)小小的程序，可以進(jìn)行許多教學(xué)討論。例如：溫度、質(zhì)量、速率區(qū)間、步長等量的變化對(duì)曲線形狀的影響；以及數(shù)值計(jì)算中如何處理無限大等。還可以計(jì)算各個(gè)速率區(qū)間內(nèi)的分子數(shù)占總分子數(shù)的百分比，并繪出相應(yīng)的分布，見圖1中最右側(cè)的圖。

采用相同的方法可以很容易地繪出許多曲線，例如拍頻(圖2)、李薩如圖形(圖3)、單縫衍射光強(qiáng)曲線(圖4)、普朗克黑體輻射曲線(圖5)等許多曲線，在此不做更多的贅述。對(duì)于學(xué)習(xí)大學(xué)物理課程的低年級(jí)學(xué)生來說，用Excel實(shí)現(xiàn)的初級(jí)可視化易學(xué)、好用。

圖3 利用Excel完成的李薩如圖形(截取了局部工作表)

圖4 利用Excel完成的單縫衍射光強(qiáng)曲線(截取了局部工作表)

圖5 普朗克黑體輻射曲線(截取了局部工作表)

2.2 利用Excel做積分

積分是大學(xué)物理課程中常用的計(jì)算，有些積分難以使用現(xiàn)成的積分公式，需要做數(shù)值計(jì)算。這一點(diǎn)限制了教材中的例題和習(xí)題的設(shè)置。使用Excel可以解決此困難。例如：某容器內(nèi)裝有溫度為T=300K的氦氣，系統(tǒng)處于平衡態(tài)。已知其中的原子個(gè)數(shù)為N=106，原子的質(zhì)量為m=6.65×10-27kg，求： (1)速率區(qū)間vp～vp+40m/s內(nèi)的原子個(gè)數(shù)，(2)速率區(qū)間10vp～10vp+40m/s內(nèi)的原子個(gè)數(shù)，vp為最概然速率。

可以利用Excel，練習(xí)數(shù)值計(jì)算中的梯形法積分法。很容易算出vp=1120m/s，要計(jì)算所求速率區(qū)間內(nèi)的原子數(shù)DN，需計(jì)算積分

圖6 利用Excel做數(shù)值積分

? 新建Excel文件，設(shè)置原子質(zhì)量MH、溫度TH、原子總數(shù)N和玻爾茲曼常數(shù)k，并定義它們的名稱，如圖6(a)。根據(jù)積分范圍，在A和C列中分別輸入所需，即速率

? 在B8和D8單元格中輸入被積函數(shù)=4*3.14*N*(MH/(2*3.14*k*TH))^(3/2)*A8^2*EXP(-MH*A8^2/(2*k*TH))

? 將鼠標(biāo)放置在B8單元格的右下角向下拖拽，計(jì)算出各個(gè)節(jié)點(diǎn)處的y值。在D列重復(fù)此操作

? 在B和D列已有數(shù)值下方的單元格中輸入計(jì)算ΔN的公式，完成積分，如圖6(b)

由B49單元中的值可知第1問的結(jié)果為30000個(gè)，由D49單元中的值可以得到第2問的結(jié)果是：原子基本上不出現(xiàn)在這個(gè)區(qū)間內(nèi)。

在此基礎(chǔ)上，可改變質(zhì)量、溫度、原子數(shù)、速率區(qū)間等量，得到所需的積分值。還可改變步長，討論它對(duì)計(jì)算結(jié)果的影響?？梢钥闯?，將解析解和數(shù)值解混搭，擴(kuò)大了解決問題的范圍。

2.3 利用Excel解微分方程

Excel可用于解一些簡單的微分方程，使學(xué)生們了解有關(guān)數(shù)值計(jì)算的簡單知識(shí)。以阻尼振動(dòng)為例。方程為：

(1)

設(shè)ω0=20s-1，對(duì)于欠阻尼，取阻尼系數(shù)β1=3s-1；對(duì)過阻尼，取β2=40s-1。初始條件為：t=0,x0=A0=0.3×10-2m，初速度v0=0。對(duì)應(yīng)該初始條件的解析解為

(2)

(3)

(4)

(5)

對(duì)于欠阻尼，用Excel實(shí)現(xiàn)計(jì)算的步驟為：

? 新建Excel文件，設(shè)置時(shí)間步長(名稱定義為DELTA)、初始條件、阻尼系數(shù)b、ω0等量的名稱和數(shù)值等，在第D6中敲入=SQRT(OMEGA0^2-BETA1^2)，也就是式(3)，計(jì)算ω的值，并定義其名稱

? 在第10行中輸入各個(gè)變量的名稱，如T、x、y等，如圖7

? 在A11、B11、C11、D11單元格中分別輸入各個(gè)量的初值

? A12單元格中輸入=A11+DELTA，并豎直拖拽鼠標(biāo)，在A列中完成時(shí)間的數(shù)據(jù)輸入

? B12單元中輸入式(2)，C12單元內(nèi)輸入式(4)，在D12單元內(nèi)輸入式(5)

圖7 利用Excel解阻尼振動(dòng)微分方程圖中x0為欠阻尼的解析解；x、 x1 和x2分別為欠阻尼、臨界阻尼和過阻尼的數(shù)值解。圖中欠阻尼的數(shù)值解幅度大于其解析解

? 同時(shí)選中B12、C12和D12單元，用鼠標(biāo)向下拖拽，拖拽的同時(shí)觀察x的值，當(dāng)它趨于零時(shí)，停止拖拽，完成計(jì)算

? 同時(shí)選中B、C、D列，點(diǎn)擊工具欄的繪圖按鈕，選擇XY散點(diǎn)圖，得圖7

? 添加文字，編輯所得圖線

通過上述步驟，采用迭代法，利用Excel解出了這個(gè)微分方程，并繪制出振動(dòng)曲線，見圖7右側(cè)。將數(shù)值解與解析解進(jìn)行對(duì)照，可驗(yàn)證數(shù)值解的正確性。還可對(duì)步長、頻率等量進(jìn)行討論，看曲線如何變化，如何影響數(shù)值解的精確度等。Excel解微分方程的步驟和迭代法的應(yīng)用很直觀。利用這種方法，學(xué)生還可以研究有空氣阻力時(shí)拋體的運(yùn)動(dòng)和行星軌道等等一些問題。

可以看出，將Excel軟件移植進(jìn)大學(xué)物理課程，進(jìn)行數(shù)值計(jì)算并非難事，除了上面提到的用途外，還可以用來解n次冪方程、非線性方程和線性方程組等。

3 對(duì)教學(xué)的促進(jìn)

利用Excel將數(shù)值計(jì)算移植到大學(xué)物理課程教學(xué)的優(yōu)勢(shì)是簡單易行并且行之有效的。對(duì)于低年級(jí)學(xué)生來說，它可提供一種解決問題的方便工具，擴(kuò)大了解決問題的范圍，并且可以借助它更好地理解物理。對(duì)于教師來說，它豐富了講授方法，使得教學(xué)更加直觀生動(dòng)。此外，它還給大學(xué)物理的教材編寫和習(xí)題配置帶來新的生機(jī)。例如；在國外一些大學(xué)物理教材中，數(shù)值計(jì)算已經(jīng)進(jìn)入課后習(xí)題，成為常規(guī)的作業(yè)題目[1]。

說到在大學(xué)物理課程中引入數(shù)值計(jì)算，我國的物理教育工作者已經(jīng)在此方向進(jìn)行了有益的探索，例如采用C語言或計(jì)算平臺(tái)[2]、Matlab[3]來實(shí)現(xiàn)等。相比于這些軟件，Excel功能差一些，如不能實(shí)現(xiàn)三維繪圖以及動(dòng)畫等。但是，它秉承了“所見即所得”的微軟Windows系統(tǒng)理念，簡單、易學(xué)，可以很輕松地融入現(xiàn)有的大學(xué)物理教學(xué)體系之中，適用于低年級(jí)理工科學(xué)生，甚至是文科學(xué)生，是在大學(xué)物理課程中普及數(shù)值計(jì)算的得力工具。

4 結(jié)語

借助Excel在大學(xué)物理課程中實(shí)現(xiàn)數(shù)值計(jì)算對(duì)于教師和學(xué)生來說都是簡單易行的，是在大學(xué)物理課程中嵌入數(shù)值計(jì)算的一條捷徑。盡管Excel的功能有限，但是對(duì)于面向低年級(jí)學(xué)生的大學(xué)物理課程來說，它能夠處理課程中所需的常用計(jì)算，并且不耗費(fèi)很多課時(shí)。打個(gè)比方來說，如果我們要去遠(yuǎn)方旅行，自駕、高鐵、飛機(jī)是我們常常要考慮的出行工具，但是如果我們要在附近活動(dòng)時(shí)，選擇騎自行車往往是明智之舉。

[1] Paul A.Tipler. Physics for scientists and engineers. Six Edition[M]. New York: W.H.Freeman and Company, 2008.

[2] 李元杰，陸果. 大學(xué)物理學(xué)[M].2版.北京：高等教育出版社, 2008.

[3] 周群益,侯兆陽,劉讓蘇. Matlab可視化大學(xué)物理學(xué)[M].北京：清華大學(xué)出版社，2011.

■

A SIMPLE WAY TO EMBED NUMERICAL CALCULATIONS INTO INTRODUCTORY PHYSICS COURSES

Liu Zhaolong Hu Haiyun

(School of Physics, Beijing Institute of Technology, Beijing 100081)

Numerical calculation is an indispensable means to study complex physics systems, but it is difficult to be reflected in the course of introductory physics due to various reasons. This paper introduces the method of embedding numerical calculation into introductory physics curriculum by using Excel software. It shows how to do visualization, numerical integration and solve numerical differential equations by Excel. The method described in this paper is simple, easy to operate, and does not need to spend a great deal of teaching hours. It is very suitable for introductory physics courses designed for junior college students. This method is aimed at helping the lower grade students know the numerical calculation earlier and learn physics more effectively, it also helps students relate what they are learning to real-world problems. It is a useful attempt on the teaching reform of the introductory physics courses.

introductory physics course; physics education; teaching reform; Excel

2017-06-14

劉兆龍，女，教授，主要從事大學(xué)物理教學(xué)和材料物理方面的研究,zlliu@bit.edu.cn。

劉兆龍，胡海云. 在大學(xué)物理課程中嵌入數(shù)值計(jì)算的捷徑[J]. 物理與工程，2017，27(5)：71-76.