方淑婷，魯曉東

(東?？茖W技術(shù)學院，浙江 舟山　316000)

?

*通訊聯(lián)系人

用Mathematica計算油滴電荷的基本量

方淑婷，魯曉東*

(東海科學技術(shù)學院，浙江 舟山316000)

摘 要：密立根油滴實驗的主要目的是要從大量的實驗數(shù)據(jù)中證明電荷的不連續(xù)性，為提高計算效率，本文利用了Mathematica軟件對實驗數(shù)據(jù)進行了最小二乘處理，優(yōu)化了誤差的分布，得到了準確的電荷基本量值。

關(guān)鍵詞：密立根油滴;電荷量;基本帶電量;Mathematica;數(shù)據(jù)處理

密立根油滴實驗，是物理學史上的經(jīng)典實驗。該實驗精確地測定了基本電荷的數(shù)值[1]，證明了所有電荷都是基本電荷的整數(shù)倍，明確了電荷的不連續(xù)性。但在實際的實驗教學過程中，往往把該公認值當作已知值代入結(jié)果表達式,僅對實驗數(shù)據(jù)進行了驗證，雖然節(jié)省了實驗時間，但不符合物理實驗數(shù)據(jù)處理的一般規(guī)則，因為事先“電荷的不連續(xù)性”及“基本電荷量公認值”是未知的，而這兩點正是本實驗需要驗證的內(nèi)容。較科學的方法應該是：使用大量的實驗數(shù)據(jù)表達出數(shù)據(jù)的分布規(guī)律，然后用最大公約數(shù)法[2]、多次逐差法[3]，統(tǒng)計法[4]等對數(shù)據(jù)進行處理，找出表征電荷不連續(xù)的特征值；如果用手工計算或借助計算機高級語言C，Baisc進行編程處理[5]，這對于低年級的學生是有困難的，影響教學效果。所以在此提出使用Mathematica軟件，該軟件的特點是表達上貼近于原始的數(shù)學表達，有利于不懂計算機語言的人使用，另外該軟件還提供了大量的數(shù)據(jù)計算的工具包，包括優(yōu)化，作圖等，適合對實驗數(shù)據(jù)的處理和可視化表達。

1油滴電荷測量原理：

密立根油滴實驗通過控制靜電場中帶電油滴在空氣中的運動，檢測相應的參量并計算得到油滴的帶電量。不妨設一個質(zhì)量為m、帶電量為q的油滴處在兩塊平行極板之間，在平行極板未加電壓時，由于空氣阻力的作用，加速一段距離后，油滴將作勻速運動，速度為vg，這時重力與阻力平衡(空氣浮力忽略不計)，根據(jù)托克斯定律，粘滯阻力為

fr=6παηvg

(1)

式中：η為空氣的粘滯系數(shù)；α為油滴的半徑。通過將油滴放在平行板上加壓后，油滴將處于場強為E的靜電場中，并加速上升，但最終在空氣阻力影響下做勻速運動即重力電場力空空氣阻力三者平衡

fr=qE-mg

(2)

根據(jù)(1)(2)以及實驗參數(shù)值[1]，當測出其勻速上升時的速度即可解出q：

(3)

同時也可以測得同一油滴所帶電荷的改變量△q(可用紫外線或放射源照射油滴，使它所帶電荷改變)，得到另一組數(shù)據(jù)。但要證實電荷的不連續(xù)性，則要取多個電量值，計算它們最大公約數(shù)，得到基本電荷電量的數(shù)值，所以油滴所帶電量的計算和分析是得到實驗結(jié)果的關(guān)鍵。從上述原理中也可以看出，在實驗過程中將會出現(xiàn)一系列影響實驗結(jié)果的誤差，計時的精度、電壓的波動以及油滴本身的揮發(fā)等，而且其方程未知數(shù)較多使得計算量較大。

2數(shù)據(jù)處理方法分析：

2.1倒過來驗證法

當我們在實驗室中測得U、t后計算得出q，以基本電荷的公認值為最大公約數(shù)去除測得油滴的電荷數(shù)n，接下來便是通過e=q/n中找出這個e0，最終得出e0約為1.60×10-19c從而驗證電荷的不連續(xù)[6]。這種以已知值去驗證已知值的方法顛倒了過程的因果關(guān)系，所以不利于培養(yǎng)學生的探索精神。

2.2油滴電量平均值逐項相減法

該方法先將帶相同電荷數(shù)的電量取平均，然后再逐項相減則得平均電量差Q 值，最后對在誤差范圍內(nèi)近似相等的幾個最小Q 值取平均即可獲得基本電荷實驗值[3]。顯然,在平均值逐項相減法中體現(xiàn)了對隨機誤差抵償性的考慮,充分利用了實驗數(shù)據(jù),從而克服了在實驗數(shù)據(jù)中帶相同電荷數(shù)油滴電量測量值不止一組情況下逐項相減法所必定產(chǎn)生的基本電荷估計值偏小的弊病，由電量逐項相減法的結(jié)果可以看出，油滴所帶的電量是不連續(xù)的。從教學效果上看,它可以使學生充分利用有限的實驗數(shù)據(jù),并有效地克服實驗儀器精度不高的負面影響，解決了教學中的一個難題。但在此基礎上還可以在數(shù)據(jù)處理過程中繼續(xù)優(yōu)化，比如在計算平均值的時候，根據(jù)該電量出現(xiàn)不同次數(shù)加以權(quán)重，對估計值的控制時選用最小二乘的原則等，這些用手工計算是很大的工程，而借助于Mathematica軟件來處理則是僅需要你的想法與步驟。

2.3基于Mathematica的方法

方法的特點就是借助軟件的功能去處理一些繁瑣的重復的計算。當實驗測得大量的數(shù)據(jù)后，先對其進行統(tǒng)計分組，即按電荷量數(shù)據(jù)分成若干等間隔區(qū)域，并按間隔進行頻數(shù)的統(tǒng)計，并且用圖的方式表達。觀測其分布，當存在基本電荷量時，數(shù)據(jù)是一簇一簇的，有明顯的不連續(xù)分布特征，接下來便可以計算簇與簇之間的距離，這就是這些數(shù)據(jù)最大公約數(shù)的粗略估計，這時各量子數(shù)也可以確定。計算的過程中我們可采取最小二乘方法對電荷基本量估計作進一步的優(yōu)化，減少誤差的影響。這樣的過程是符合了實驗探索思路，盡管計算量很大但這些都可以借助軟件來完成。以我們13級學生的數(shù)據(jù)為例，對各班近150組的實驗數(shù)據(jù)進行如下的處理：(數(shù)據(jù)存放由于exdata.xlsx中)

(1)根據(jù)平衡電壓和運動時間計算油滴帶電量

ρ=981;g=9.793;η=1.8310-5;l=2.10-3;b=6.1710-6;p=76;d=5.10-3;(*實驗條件*)

mydata=Import["E:\exdata.xlsx",{"Data",1}];(*打開excel數(shù)據(jù)文件*)

myq=Apply[q,mydata,1];(*把實驗數(shù)據(jù)代入公式計算q值,存到變量myq*)

(2)根據(jù)電荷量進行統(tǒng)計分組

由實驗數(shù)據(jù)的范圍組數(shù)，以及儀器(MOD-5型)的精度和實驗過程中的不確定度的估計，取間隔為0.05(×10-19c)進行統(tǒng)計如圖所示，圖1顯示了各數(shù)據(jù)簇的不連續(xù)性，

TheArea=Table[0.05k,{k,0,200}];(*統(tǒng)計電荷量間隔的分布存到TheArea*)

Ncount=BinCounts[myq,{TheArea}];(*統(tǒng)計電荷量myq在各區(qū)間內(nèi)的頻數(shù)存到Ncount*)

listdata=Table[{TheArea[[n]],Ncount[[n]]},{n,200}];

ListPlot[listdata,Filling→Axis,PlotRange→{{0,10},{0,10}}}](*統(tǒng)計圖形顯示*)

(3)電荷基本量的估計

(4)加權(quán)最小二乘法

w=∑pi(qi-nie0)2→min

在軟件中可以直接調(diào)用數(shù)值計算工具[7]minimize[]就可以了。其使用格式為：Minimize[f,{x,y,...}]，f為要優(yōu)化的表達式即∑pi(qi-nie0)2，{x,y,…}為要估計的參數(shù)即e0，函數(shù)結(jié)果返回表達式的最小值。

p={0.087,0.196,0.348,0.239,0.13};(*出現(xiàn)的概率，見表1*)

q={1.57,3.125,4.94,6.26,8.11};(*各數(shù)據(jù)簇的期望值，見表1*)

n=Round[q/del];(*電荷量的量子數(shù)*)

結(jié)果為{0.01,{e0→1.604}}表示電荷的基本量為1.604×10-19c，擬合的殘差方和為0.01，與公認值比較百分差為0.25%。若按已知的理想值的公認值去估計本次實驗的結(jié)果則為1.59,百分差為0.6%，屬于同一級的精度，所以借助Mathematica計算的結(jié)果是比較理想的，并且其處理過程符合科學實驗的要求[8-9]。

3結(jié)束語

密立根油滴實驗是通過對宏觀量的觀測來揭示微觀量本身所具有的屬性，所以必定要通過對大量觀測數(shù)據(jù)的計算和統(tǒng)計。用計算機軟件輔助實驗數(shù)據(jù)處理，不僅僅是提高計算的效率，更重要的是可以挖掘數(shù)據(jù)背后大量的信息，揭示其中的物理規(guī)律。Mathematica 正是這樣的軟件，它提供了大量直接可調(diào)用的數(shù)學工具包，為實驗的進一步探索創(chuàng)造了條件。

參考文獻：

[1]竺江峰,蘆立娟,魯曉東.大學物理實驗[M].北京:中國科學技術(shù)出版社,2005.

[2]王廣濤,陳 健,魏建宇等.密立根油滴實驗數(shù)據(jù)的處理方法[J].物理實驗,2004,24(12):22-24.

[3]劉才明.密立根油滴實驗數(shù)據(jù)處理方法分析[J].浙江大學學報：自然科學版,1996,30(6):736-741.

[4]溫猛,洪朱旭,馮運軍,等.密立根油滴實驗的概率統(tǒng)計分析[J].中山大學學報：自然科學版,2004(S1):33-35.

[5]段蘋,王吉有,秦麗,等.計算機編程在密立根油滴實驗中的應用[J].大學物理實驗,2014,27(1):70-73.

[6]師文慶,呂楠,陳勁民.驗證法處理密立根油滴實驗探討[J].廣東海洋大學學報,2009,29(4):67-69.

[7]張韻華，王新茂.Mathematica 7 使用教程[M].安徽：中國科學技術(shù)大學出版社，2012.

[8]張鑫，郭勝，李鈺.C語言在密立根實驗數(shù)據(jù)處理中的應用[J].大學物理實驗，2015(1)：63-65.

[9]林曼虹，黃燒燦，林豐，等.利用精密秤測量微小電量[J].大學物理實驗，2015(1)：29-31.

Calculating the Base Charge of Oil Drop by Mathematica

FANG Shu-ting,LU Xiao-dong

(Donghai Science and Technology College,Zhejiang Zhoushan 316000)

Key words：Milliken oil drop；quantity of charge；basic charge；Mathematica；data processing

Abstract：Millikan oil drop experiment is aimed at proving the uncontinuous change carried by each oil drop.In order to improve calculation efficiency for these lots of experiment data,a method based on Mathematica is provided.By its way of MLS processing and error optimizing,the accurately value of base charge is calculated.

收稿日期：2016-02-24

基金項目：浙江省2015年度高等教育課堂教學改革項目(kg2015577)

文章編號：1007-2934(2016)03-0107-04

中圖分類號：O 4-39

文獻標志碼：A

DOI：10.14139/j.cnki.cn22-1228.2016.003.028