梁 雄，賴國忠

(龍巖學(xué)院 福建龍巖 364000)

塞曼效應(yīng)為研究原子結(jié)構(gòu)提供了重要途徑[1]，橫向塞曼效應(yīng)被用于檢測大氣汞[2]和煙氣汞[3]，除此之外，塞曼效應(yīng)還可用于測量超強(qiáng)磁場[4]。近代物理實(shí)驗(yàn)課程常開設(shè)塞曼效應(yīng)實(shí)驗(yàn)，用于測量電子的荷質(zhì)比，測量過程有個主要任務(wù)就是測量干涉環(huán)的直徑，測量過程數(shù)據(jù)比較多，耗時長。人們?yōu)榱吮阌谟^測或減輕觀測任務(wù)，提出了各種觀測方法：用阿貝比長儀觀測分裂譜底片[5]，利用五點(diǎn)中心測量法[6]優(yōu)化了測量過程，基于CCD傳感器技術(shù)，結(jié)合數(shù)字圖像處理技術(shù)提出了半自動[7-8]或自動[9-10]測直徑，也有基于LabVIEW結(jié)合手動取三點(diǎn)測直徑[11]。 雖然出現(xiàn)了全自動測直徑的方法，但是在學(xué)生實(shí)驗(yàn)中采用手動測直徑，一方面可以增強(qiáng)學(xué)生的動手能力，另一方面可以加深對物理規(guī)律的理解。三點(diǎn)法還被應(yīng)用到19JPC微機(jī)型萬用工具顯微鏡測樣品上的圓斑或圓孔的直徑，其組成與工作原理可參見文獻(xiàn)[12]。在塞曼效應(yīng)實(shí)驗(yàn)的測量過程中，發(fā)現(xiàn)大多數(shù)學(xué)生使用FD-FZ-C法拉第效應(yīng)塞曼效應(yīng)綜合實(shí)驗(yàn)儀時，常常按圖1 (a)或(b)所示的方式來選取A、B和C三點(diǎn)，只有少數(shù)學(xué)生采用如圖1(c)那種方式來取點(diǎn)。

圖1 三點(diǎn)取法對比示意圖

考慮到選點(diǎn)時受人眼(或儀器)分辨率的影響導(dǎo)致的位置誤差，這三種選點(diǎn)方式是否具有一樣的直徑測量精確度？本文從誤差傳遞公式及幾何知識來分析所取三點(diǎn)的連線構(gòu)成什么樣的三角形時，直徑測量誤差最小以及這三點(diǎn)選取誤差各自對直徑測量精確度有怎樣的影響,最后應(yīng)用Matlab軟件模擬干涉環(huán)的測量來驗(yàn)證討論所得的結(jié)論。

圖2 直角頂點(diǎn)偏離示意圖

1 三點(diǎn)法測直徑原理及其優(yōu)化

1.1 三點(diǎn)法測直徑原理

如圖1所示，在待測圓上任取三點(diǎn)A、B和C，該圓是三角形ABC的外接圓，弦AB的長度設(shè)為l， ∠ACB=α(下文中α均指∠ACB)，容易由正弦定理將待測圓的直徑表示成：

d=l/sinα，

(1)

除了按式(1)計算直徑外，還可根據(jù)所取的三點(diǎn)要滿足圓的方程而列出三個方程，從該方程組可先求得圓心坐標(biāo)，再求出圓心到這三點(diǎn)中的任一點(diǎn)之間距離即可求得相應(yīng)環(huán)的半徑[11]。一旦從待測圓弧上選定三點(diǎn)，不同計算法求解的直徑結(jié)果是一樣的，即與真值之間的偏差也相同。估算誤差時，由于不同計算方法對應(yīng)的誤差傳遞公式不同，估算誤差也會不完全一致，但是它們估算的誤差水平應(yīng)是相當(dāng)?shù)?，所以為了便于討論選點(diǎn)位置的誤差對直徑結(jié)果的影響，本文采用式(1)的誤差傳遞公式來討論三點(diǎn)應(yīng)該在圓上的什么位置選取可以將直徑的測量誤差最小化，而后一種方法便于在Matlab編程求解直徑，因此模擬測量過程將采用后面這種方法驗(yàn)證誤差分析所得的優(yōu)化方法。

1.2 三點(diǎn)法測直徑誤差及三點(diǎn)最佳取法

由誤差傳遞公式容易求出(1)式的直徑測量誤差為：

(2)

其中uc(l)和uc(α)分別代表測量弦長l和α的綜合不確定度。結(jié)合圖1可知，弦長l越大，角α就越接近90°，當(dāng)弦長取最大值時，有l(wèi)=d和α=90°，在等精度測量的前提下，此時上式根號里頭兩項(xiàng)都達(dá)到最小。由此看來，三點(diǎn)法測直徑時A、B兩點(diǎn)應(yīng)盡可能選取在直徑端點(diǎn)位置，它們與另外一點(diǎn)C的連線構(gòu)成直角三角形，此時直徑的測量誤差最小。

為了便于討論直角頂點(diǎn)C的誤差對直徑的影響，假設(shè)AB準(zhǔn)確地落在半徑為R待測圓直徑的端點(diǎn)上(即AB=2R)，而選取頂點(diǎn)C受人眼分辨率或儀器分辨率的影響將會導(dǎo)致該點(diǎn)偏離待測圓，距離記作Δ(如圖2中線段CD)，因此選取C點(diǎn)誤差會導(dǎo)致角α與直角存在一定的偏差。

不妨假設(shè)圓上不同位置選取C點(diǎn)時不確定度Δ相當(dāng)，令∠BOC=θ，∠BCO=α1，OC記作r。 在三角形BOC中由正弦定理得：

可求得

r=Rsin(θ+α1)/sinα1，

兩邊取微分得(在此假設(shè)θ為常數(shù)),并令dr=Δ，α1≈(π-θ)/2，整理得：

(3)

在三角形AOC中,令∠ACO=α2，同理可求得：

(4)

考慮判斷誤差Δ通常較小，即R?Δ或Δ/R?1，此時角α與直角有一個較小的偏離ε1：

(5)

其中負(fù)號含義：當(dāng)Δ>0，即選取直角頂點(diǎn)C在圓外，則導(dǎo)致α略小于π/2；反之，導(dǎo)致α略大于π/2.由此可見：倘若在圓弧上A、B兩點(diǎn)鄰近區(qū)域選取直角頂點(diǎn)C時，此時θ接近0°或180°，在同樣半徑偏離Δ情況下，將導(dǎo)致較大的角度誤差ε1。而在圓弧上θ=90°的位置選擇直角頂點(diǎn)C，結(jié)合圖2幾何關(guān)系可知：綜合誤差公式(2)的分析結(jié)果，所選三點(diǎn)的連線構(gòu)成一個等腰直角三角形，此時由選點(diǎn)位置誤差造成α角的誤差最小：ε1=-Δ/R， 將其代入誤差傳遞公式(2)可得測量直徑的相對誤差：

(6)

式(6)推導(dǎo)中應(yīng)用了cotα=cot(π/2+ε1)≈-ε1。

圖3 直徑方向上點(diǎn)偏離圓對角度的影響

1.3 直徑兩端取點(diǎn)對直徑測量誤差的影響

在直徑端點(diǎn)處取兩點(diǎn)A、B，考慮這兩點(diǎn)在圓弧上位置具有對稱性，它們位置誤差相當(dāng)，在此僅以B點(diǎn)為例進(jìn)行說明。如圖3所示，假設(shè)A、C兩點(diǎn)準(zhǔn)確無誤地落在圓上，在待測圓上位置D選擇B時，由于判斷誤差將會沿著直徑方向偏離到B點(diǎn)(注：也有可能落在圓內(nèi))，偏離距離記作Δ，容易知道∠ACB比90°大，導(dǎo)致直角誤差ε2。

(7)

對式(7)取微分可得：

(8)

(9)

(10)

2 干涉環(huán)直徑測量模擬

由于實(shí)際的塞曼效應(yīng)實(shí)驗(yàn)中無法事先確切地知道法布里—珀羅腔(F-P腔)干涉圖樣中干涉環(huán)的直徑真值，這將造成無法計算每次測量誤差，以及難以判斷多次測量算術(shù)平均值哪種取點(diǎn)方法更接近真值。為了克服這些不便，利用Matlab模擬塞曼效應(yīng)的干涉圖樣，各參數(shù)按照文獻(xiàn)[13]。具體選取情況：F-P腔的反射率設(shè)為90%，折射率n=1.5，腔厚d=5.55 mm，成像透鏡焦距f=200 mm，波長λ=555 nm，模擬F-P干涉圖樣的結(jié)果如圖4所示。測量過程選取圖4中最接近中央亮斑的干涉環(huán)，其直徑可由理論計算得：d0=3.266 mm(將它視作真實(shí)值，即可求出測量誤差)，應(yīng)用軟件中g(shù)input命令按圖1中三種采點(diǎn)方式取點(diǎn)，按從左到右將這三種方式簡稱M1、M2和M3，每種方法取30組(其中每組數(shù)據(jù)由三個點(diǎn)坐標(biāo)組成)，即每種方法測量圓直徑30次。

圖4 F-P腔干涉模擬圖

將相應(yīng)的直徑測量值減去理論值即可求得每次測量誤差，以數(shù)據(jù)編號i為橫坐標(biāo)，將三種取點(diǎn)方式的直徑測量誤差繪制成曲線(見圖5(a))：M1測量誤差曲線對應(yīng)左側(cè)縱坐標(biāo)軸u(d)，另外兩種方式M1、M3測量誤差曲線對應(yīng)右側(cè)縱坐標(biāo)軸u1(d)。

從圖5(a)中可以看出，M1和M2方式測量結(jié)果誤差的波動幅度都比較大，而M3方式取點(diǎn)每次測得直徑與理論值之間的偏差都比較小。這三種取點(diǎn)方法測得直徑的平均值分別為3.371 mm、3.272 mm和3.270 mm，其中M3平均值最接近待測干涉環(huán)直徑的真值(3.266 mm)。此外用標(biāo)準(zhǔn)偏差來衡量波動幅度，它們分別為： 0.730 mm(M1)、0.007 mm(M2)和0.004 mm(M3)。這是因?yàn)椋篗1取點(diǎn)所測的弦長比直徑小得多，按照式(1)可知測得直徑必須按某個比1大得多的因子(1/sinα)放大后求得，放大過程同樣將弦長l的測量誤差進(jìn)行了放大，而M2與M3雖然都有較大的弦長(甚至接近直徑)，但由于方式2中直角頂點(diǎn)非?？拷睆蕉它c(diǎn)，此時由于直角頂點(diǎn)選取誤差導(dǎo)致α偏離90°較大，導(dǎo)致M2的直徑測量誤差要大于M3的誤差，這與前面誤差分析所得的結(jié)論完全一致。

三種方式檢測的圓心坐標(biāo)情況，用測量所得的圓心坐標(biāo)與實(shí)際圓心坐標(biāo)(0,0)之間的距離e衡量圓心的檢測誤差，如圖5(b)所示，對比后容易發(fā)現(xiàn)：采用等腰直角三角形的方式取點(diǎn)，M3的圓心檢測誤差e要比其他兩種方式取點(diǎn)所得的結(jié)果要小得多。

(a) 直徑測量誤差曲線 (b) 圓心測量誤差曲線

3 結(jié)論

應(yīng)用三點(diǎn)法測量待測圓(例如干涉環(huán))的直徑或圓心，采用等腰直角三角形的選點(diǎn)方式具有較高的精確度。此外，建議用FD-FZ-C法拉第效應(yīng)塞曼效應(yīng)綜合實(shí)驗(yàn)儀做塞曼效應(yīng)實(shí)驗(yàn)時，首先在干涉圖樣中選取較大的干涉環(huán)，按照優(yōu)化后的方法測定環(huán)心坐標(biāo)，再沿著垂直直徑AB方向在每個待測環(huán)上選取一個點(diǎn)，這樣可以比較精確地測定相應(yīng)待測環(huán)直徑，同時也有利于減小電子的荷質(zhì)比測量誤差。順便提一下，在估算選點(diǎn)的誤差對直徑測量影響時，按照優(yōu)化取點(diǎn)方式，直徑測量誤差主要來源于直徑端點(diǎn)處所取兩點(diǎn)的位置誤差，而直角頂點(diǎn)的位置誤差對直徑測量精確度影響比較小，基本可以忽略。