Matlab軟件在牛頓環(huán)實(shí)驗(yàn)數(shù)據(jù)處理方法上的創(chuàng)新

徐少剛，夏雪琴

(浙江海洋學(xué)院，浙江舟山 316022)

牛頓環(huán)實(shí)驗(yàn)是一個(gè)古老而又非常重要的光學(xué)實(shí)驗(yàn)，一般用它來(lái)測(cè)量平凸透鏡的曲率半徑十分方便，這個(gè)實(shí)驗(yàn)在我校已開設(shè)許多年了。目前數(shù)據(jù)處理方法主要有:逐差法和最小二乘法，用人工計(jì)算既繁瑣又篇幅大，占用了我們學(xué)生較多的學(xué)習(xí)時(shí)間，并且往往運(yùn)算容易出差錯(cuò)，得不到十分理想的實(shí)驗(yàn)結(jié)果，這樣必然會(huì)影響學(xué)生的學(xué)習(xí)積極性和教學(xué)效果。隨著計(jì)算機(jī)軟件的迅速發(fā)展，如用Matlab軟件的功能函數(shù)進(jìn)行最小二乘法曲線擬合，得到平凸透鏡的曲率半徑，而且還能進(jìn)行誤差理論分析，一方面簡(jiǎn)化了數(shù)據(jù)處理過程，提高了實(shí)驗(yàn)結(jié)果的準(zhǔn)確度;另一方面提高了大學(xué)生應(yīng)用軟件的能力，并且對(duì)實(shí)驗(yàn)數(shù)據(jù)處理方法上的一種創(chuàng)新，激發(fā)了大學(xué)生對(duì)學(xué)習(xí)大學(xué)物理實(shí)驗(yàn)的興趣，提高了教學(xué)質(zhì)量。

1 實(shí)驗(yàn)原理

在一塊平面玻璃上安放上一焦距很大的平凸透鏡，使其凸面與平面相接觸，在接觸點(diǎn)附近就形成一層空氣膜。當(dāng)用一平行的準(zhǔn)單色光垂直照射時(shí)，在空氣膜上表面反射的光束和下表面反射的光束在膜上表面相遇相干，形成以接觸點(diǎn)為圓心的明暗相間的環(huán)狀干涉圖樣，稱為牛頓環(huán)，其如圖1。

圖1 牛頓環(huán)光路示意圖

可改為 D2K=4KRλ(K=0，1，2，…)

用Matlab軟件作圖并處理實(shí)驗(yàn)數(shù)據(jù)，作出D2K～K的關(guān)系圖。

令y=D2K;x=K;k=4Rλ，可得到擬合直線方程y=kx+b，并從其斜率k中求出平凸透鏡曲率半徑R的值。

2 實(shí)驗(yàn)數(shù)據(jù)記錄

已知鈉光的波長(zhǎng) λ =5.893×10－7m，儀器誤差5 × 10－6m，R標(biāo)=0.880 0 m，見表1。

表1 透鏡曲率半徑測(cè)量數(shù)據(jù)

3 Matlab軟件處理實(shí)驗(yàn)數(shù)據(jù)

3.1 得到y(tǒng)與x的數(shù)據(jù)表

令DK=XL－XR，x=K，y=DK2，用Matlab軟件編寫程序計(jì)算:

＞ ＞ x= ［55 50 45 40 35 30 25 20 15 10］

＞ ＞ c1= ［0.032394 0.032112 0.031902 0.031638 0.031338 0.031036 0.030702 0.030318 0.029918 0.02941］

%XL的值

＞ ＞ c2= ［0.021538 0.021786 0.022062 0.022314 0.0226 0.022918 0.023338 0.023602 0.024026 0.024506］

%XR的值

＞＞c3=c1－c2

%DK=(XL－XR)的值

＞ ＞ y=c3.^2%

y-DK2的值

數(shù)據(jù)見表 2［2］。

表2 由D2K與K之間關(guān)系得到y(tǒng)與x的數(shù)據(jù)

3.2 用polyfit功能函數(shù)作y與x的擬合圖

3.2.1 用Matlab軟件最小二乘法公式求出斜率k、相關(guān)系數(shù)γ

y=kx+b用最小二乘法求出曲線的斜率k、相關(guān)系數(shù)γ，

Matlab程序:

＞ ＞ a=polyfit(x，y，1)

結(jié)果得:k=0.2079e－ 5，b=0.3307e－ 5

3.2.2 用plot函數(shù)繪y與x的擬合圖2

Matlab程序:

＞ ＞x1=10:5:55;k=2，

＞＞y1=a(1)*x1+a(2);

＞ ＞ plot(x，y，’* ’，x1，y1，’k’)，xlabel(‘P’)，ylabel(‘I’)，grid on，

Matlab擬合圖，見圖2。

圖2 Matlab繪y與x的擬合圖

所以擬合直線方程為:y=0.207 9×10－5x+0.330 7 × 10－5

4 Matlb軟件最小二乘法誤差分析計(jì)算

4.1 Matlab 程序

＞ ＞ yn=a(1)*x+a(2)，xn=(y－a(2))/a(1)，

＞ ＞ sigmayn=sqrt((yn－y)*(yn－y)’/(k*(k－1)))，

%DK的絕對(duì)誤差ΔY

＞ ＞ sigmaxn=sqrt((xn－x)*(xn－x)’/(k*(k－1)))，

%K的絕對(duì)誤差ΔX

＞ ＞ xp=mean(x)，yp=mean(y)，

Lxx=(x－xp)*(x－xp)’，

%

＞ ＞ Lxy=(x－xp)*(y－yp)’，

Lyy=(y－yp)*(y－yp)’，

＞ ＞ r=Lxy/sqrt(Lxx*Lyy)，

線性相關(guān)系數(shù)r

＞ ＞ R=a(1)/(4*5.893*0.0000001)，

%曲率半徑ˉR

＞ ＞ R0=0.88，deth=(R － R0)/R0，

%百分差ER

＞ ＞ y2=0.000107，x2=50，b=a(2)，

＞＞sigmah=R*sqrt((sigmayn/(y2－b))^2+(sigmaxn/x2)^2)，

%ΔR

＞＞deth1=sigmah/R，

%相對(duì)不確定度EˉR

4.2 誤差結(jié)果得

sigmayn=9.9912e － 007

sigmaxn=0.4807

sigmah=0.012

r=0.9999

R=0.8818

deth=0.0020

death1=0.0136

線性相關(guān)系數(shù)r=0.9999接近于1，說明k、b的值準(zhǔn)確性高且可用，擬合直線方程為:y=0.2079 ×10－5x+0.3307 ×10－5線性擬合度高。

曲率半徑的最佳值:ˉR0.881 8 m，

D2K的絕對(duì)誤差 ΔY=sigmayn=9.991 2e－007，

K 的絕對(duì)誤差 ΔX=sigmaxn=0.480 7，R 的絕對(duì)誤差 ΔR=sigmah=0.012。

5 實(shí)驗(yàn)結(jié)論

通過上面用Matlab軟件的功能函數(shù)對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行處理，很快得到如下實(shí)驗(yàn)結(jié)果:

平凸透鏡的曲率半徑:

R=ˉR±ΔR=(8.82±0.12)×10－1m

相對(duì)不確定度:

百分差:

根據(jù)上面的數(shù)據(jù)處理和誤差分析可知，Matlab軟件的功能函數(shù)最小二乘法對(duì)牛頓環(huán)實(shí)驗(yàn)進(jìn)行數(shù)據(jù)處理精度高，而且程序運(yùn)行簡(jiǎn)便，計(jì)算快捷、高效;因此，大學(xué)生在物理實(shí)驗(yàn)中用Matlab軟件對(duì)牛頓環(huán)實(shí)驗(yàn)進(jìn)行數(shù)據(jù)處理，是對(duì)傳統(tǒng)的數(shù)據(jù)處理方法上的一種創(chuàng)新，提高了大學(xué)生學(xué)習(xí)積極性，同時(shí)也提高了大學(xué)物理實(shí)驗(yàn)的教學(xué)質(zhì)量。

［1］王曉雄.利用Origin擬合方法處理牛頓環(huán)數(shù)據(jù)［J］.大學(xué)物理實(shí)驗(yàn)，2011，24(2).

［2］劉浩.新型基于彩色熒光燈的牛頓環(huán)實(shí)驗(yàn)［J］.大學(xué)物理實(shí)驗(yàn)，2012，25(5).