潘尹凡 瞿 霞 馮 杰 郭長江

(上海師范大學(xué)數(shù)理學(xué)院 上海 200234)

邁克耳孫-莫雷實驗的預(yù)期與實際模擬

潘尹凡 瞿 霞 馮 杰 郭長江*

(上海師范大學(xué)數(shù)理學(xué)院 上海 200234)

通過對邁克耳孫-莫雷實驗的預(yù)期與實際計算機模擬對比，凸顯經(jīng)典時空觀的問題.同時補充了普遍僅對中心位置分析計算的遺憾，給出了完整的模擬圖樣，以幫助理解該實驗的思想原理.

邁克耳孫-莫雷實驗 預(yù)期模擬 Matlab GUI

1900年4月27日，開爾文男爵在他的演講中道出了物理學(xué)歷史中經(jīng)久不衰的比喻——19世紀(jì)籠罩在熱和光動力理論上空的烏云.那兩朵烏云如同是日后物理學(xué)書籍中必不可少的一種定式化修辭，一次次重現(xiàn)著狂風(fēng)暴雨前夜物理學(xué)界自滿又不安的空氣[1].其中的一朵烏云是1887年由邁克耳孫與莫雷教授重做的“以太漂移”實驗，即著名的邁克耳孫-莫雷實驗，該實驗動搖了絕對靜止參考系的時空觀，對日后相對論的產(chǎn)生有著不可忽視的作用.

1 對條紋變化分析的不足

目前部分物理書籍在論述邁克耳孫-莫雷實驗時，為得出0.37個條紋移動的預(yù)期結(jié)果，多采用如下方式得出結(jié)論.

如圖1(a)所示,光經(jīng)分光板后，透射光先“追光”，后經(jīng)反射“迎光”，整個過程耗時

(1)

其中v為地球相對“以太系”的運行速度，c為光速，L1為透射光到反射鏡M1的距離.對于反射光，如圖1(b)所示，當(dāng)時認(rèn)為在以太系中光速為c，故在豎直方向上有分量，應(yīng)有

(2)

圖1 實驗光路圖

考慮到補償板已插入，兩路光的光程差即為Δδ=c(τ1-τ2).裝置旋轉(zhuǎn)90°后光路情況交換，故光程差變化為2Δδ，實驗中L1與L2取為11 m，若v大小取30 000 m/s，代入算式，旋轉(zhuǎn)前后光程變化對應(yīng)約0.37個波長，故應(yīng)有0.37個條紋移動[2].

上述方法僅選擇了中心位置進行分析，并未對其他近軸點進行一般分析，如此就可能產(chǎn)生疑問：“也許只有這中心點會有亮度改變，而旋轉(zhuǎn)90°對中心附近點的亮度分布幾乎沒有影響，導(dǎo)致我們?nèi)庋塾^察不到中心點的變化，所以才得出了零結(jié)果的結(jié)論？”或是：“是否因為旋轉(zhuǎn)前裝置與以太風(fēng)恰成45°角，導(dǎo)致旋轉(zhuǎn)90°后圖樣確實不會有變化？”誠然，從簡單定量角度來看，該方法可最快給出結(jié)論，然而在未嚴(yán)格給出總體變化規(guī)律前，將零維特殊點的變化直接上升到二維難免有失嚴(yán)謹(jǐn).

2 基于當(dāng)時時空觀的計算

我們可通過鏡像法將整個光路化為兩個點光源在觀察屏上方成干涉圖樣的問題，如圖2所示.

圖2 邁克耳孫-莫雷實驗的細化光路圖

需要注意的是，由于經(jīng)過了鏡像處理，部分光路等效的v方向會發(fā)生改變，同時在裝置旋轉(zhuǎn)90°時，其變化方向也不相同，具體規(guī)律如表1所示(為不至于混亂，暫以兩點光源重合于S″分析).

表1 等效光路理論變化情況

現(xiàn)不妨設(shè)L1與v所成夾角為α，此時由圖3可知：光在“以太系”中的走向為K→M1→K2，只需使用余弦定理即可求出KM1的長度，即對應(yīng)的實際光程.

由圖3得

則其實際光程

(3)

圖3 實際光程理論示意圖

由式(3)，只要知道某等效光線在“地球系”走過的長度L，及L與相對運動速度v所成的夾角α，便可得出在“以太系”下傳播的實際光程.

繼續(xù)計算各段的L與α.先來分析各段等效光路在“地球系”的傳播長度L.對觀察屏建立如圖4所示的三維直角坐標(biāo)系：定E點為坐標(biāo)原點，x軸正方向與v同向，K→M2方向為z軸正方向.為方便計算，設(shè)

h=KEH=S″E=SK+2KM2+KE

(4)

圖4 地球系下近軸光線的理論長度

(5)

則各段長為

(6)

則

(7)

式中h*=K′E,至此各段長度均可表示.繼續(xù)分析L與相對運動速度v所成夾角α的情況.如圖5所示，假設(shè)等效v沿HC方向(可以證明v恒在xOz平面內(nèi))，HJ為地球系下觀察到的光傳播方向，故α=∠CHJ，應(yīng)用兩次余弦定理便可求出.先對△COJ分析

CJ2=H2tan2θ+a2+b2+2aHtanθ

(8)

在△CHJ中

(9)

圖5 夾角計算

同理，v指向x正半軸時

3 基于Matlab的模擬

理論計算完備后，便可進行模擬，本文使用Matlab2011b版本運算，并為保證運算速度，采樣點選為28+1即257個.當(dāng)實驗光波長選擇為590 nm,臂長取為11 m與11.001 m,觀察屏半寬度0.8 m,分光板距觀察屏15 m,相對速度v取地球公轉(zhuǎn)速度30 000 m/s時，預(yù)期的不同角度下模擬圖樣應(yīng)如圖6所示.

圖6 在上述條件下，轉(zhuǎn)過不同角度時的預(yù)期圖樣

可以看出，干涉圖樣應(yīng)有明顯的變化，變化表現(xiàn)在條紋不斷被“吞入”.外圍的條紋逐漸向中心靠攏，越外圈的條紋移動程度越不明顯.有了該變化規(guī)律，才能說用中心點來計算預(yù)期變化是可行的，同時上述一系列的疑問自然也就解開了.

4 結(jié)論補充

上文雖已給出了邁克耳孫-莫雷實驗的預(yù)期模擬，但在實際中還存在著幾個值得注意的問題，接下來從3個問題切入，并結(jié)合模擬來逐一分析.

(1)表述問題.在上文中已得出了理論上變化的規(guī)律，應(yīng)是條紋“吞吐”的變化，而目前部分文章中表述為條紋的“移動”，移動是一個很模糊的描述，是需要避免的誤區(qū).

(2)臂長問題.在通過公式得出條紋變化條數(shù)時，我們常將L1與L2直接取為11 m,并得出0.37個條紋的結(jié)論.該步驟從物理意義角度來看并不妥當(dāng)，假若兩臂均取11 m,代入模擬并不能形成干涉條紋，因為兩個等效光源重合為了一點.因此,本實驗中臂長不等是產(chǎn)生現(xiàn)象的關(guān)鍵前提，當(dāng)L1為11 m時,L2必須為11 m±δ，δ是不可消除誤差，若不考慮該誤差，列出的計算式?jīng)]有對應(yīng)的實際意義.從這一角度來看，實驗中旋轉(zhuǎn)的操作也不是為了“消除”臂長誤差，主要是為交換兩光路的傳播方式.

圖7 v取750 m/s時，轉(zhuǎn)過角度后也幾乎沒有變化

(3)“零結(jié)果”問題.邁克耳孫-莫雷實驗的結(jié)果在當(dāng)時并沒有完全否定以太的存在，邁克耳孫本是為了測量地球相對以太的運動速度，但最后不成想只觀測到了極小的條紋變化，或才有了零結(jié)果這一說法.若我們在模擬中保持其他參數(shù)不變，將v修改

5 GUI模擬界面

利用Matlab軟件自帶的GUI功能生成一個模擬窗口(圖8)，可快速調(diào)整各個參數(shù)的數(shù)值并觀察其對圖樣的影響.本界面將旋轉(zhuǎn)角度、相對以太速度、兩臂長、屏寬這5個主要參數(shù)以滑塊方式集合.

圖8 邁克耳孫-莫雷實驗的預(yù)期模擬GUI界面

1 曹天元.量子物理史話 上帝擲骰子嗎.沈陽：遼寧教育出版社，2008.28～32

2 馮杰.大學(xué)物理專題硏究.北京：北京大學(xué)出版社，2011.291～297

3 漆安慎，杜嬋英.力學(xué)(第2版).北京：高等教育出版社，2005.422～425

4 姚啟鈞.光學(xué)教程(第4版).北京：高等教育出版社，2008.42～44

2017-02-25)

潘尹凡(1994- )，男，在讀碩士研究生，研究方向為物理學(xué)科教學(xué).

指導(dǎo)教師：馮杰(1961- )，男，教授，碩士生導(dǎo)師,研究方向為理論物理、非線性光學(xué)、物理課程與教學(xué)論.