馮明春 王玉杰

(滁州學院機械與電氣工程學院 安徽 滁州 239000)

測量不確定度是指測量結果變化的不肯定，是表征被測量的真值在某個量值范圍的一個估計，是測量結果含有的一個參數(shù)，用以表示被測量值的分散性[1,2].測量不確定度的定義表明，一個完整的測量結果應包含被測量值的估計與分散性參數(shù)兩部分.例如，被測量Y的測量結果為y±u，其中y是被測量值的估計，u為被測量的不確定度.文獻[3～5]著重分析了不確定度的一些理論；而本文在不確定度理論基礎上再結合具體實驗數(shù)據(jù)，分別對力學、熱學和光學實驗中所對應的彈性模量、液體表面張力系數(shù)和牛頓環(huán)實驗中平凸透鏡的曲率半徑這3個實驗測量結果的不確定度進行了分析討論，這有助于教師實驗教學和學生分析實驗數(shù)據(jù)，使得研究更有意義.

1 分析過程

根據(jù)測量不確定度的定義，在測量實踐中如何對測量不確定度進行合理的評定，這是必須解決的基本問題.對于一個實際測量過程，測量不確定度的來源有多個，這些不同來源的不確定度在計算方法上只有兩類，一類稱為A類分量，它是用統(tǒng)計學方法計算的分量，是隨機誤差性質的不確定度；另一類稱為B類分量，是用非統(tǒng)計方法評定的分量，是系統(tǒng)誤差性質的不確定度[6].計算不確定度，常用計算標準差去表示，稱為標準不確定度.

1.1 直接測量值的標準不確定度的A類分量uA(x)

取x的標準不確定度的A類分量為

(1)

1.2 直接測量值的標準不確定度的B類分量uB(x)

設x誤差的某一項的誤差限為Δ，其標準差

則標準不確定度的B類分量

(2)

式(2)中，k稱為包含因子，k值與x的分布有關，其常用分布及其k值如表1所示.

表1 常用分布及其k

1.3 合成標準不確定度uc(x)或uc(y)

對某一物理量測量之后，要計算測量值的不確定度，由于其測量值的不確定度來源不止一個，所以要合成其標準不確定度.

對于直接測量，設被測量X的標準不確定度的來源有m項，則合成標準不確定度uc(x)為

(3)

式(3)中的u(x)可以是A類評定或B類評定.

對于間接測量，設被測量Y由n個不相關的直接被測量x1，x2，…，xn算出，它們的關系為y=y(x1，x2，…，xn)，各xi的標準不確定度為u(xi)，則y的合成標準不確定度uc(y)為

(4)

1.4 不確定度的報告

Y=y±uc(y)

(5)

式(5)就是被測量的不確定度的表現(xiàn)形式.

2 實驗分析

2.1 彈性模量的標準不確定度

金屬絲彈性模量E的表達式為[6]

(6)

在式(6)中，F(xiàn)為外力，l為金屬絲的長度，d為金屬絲的直徑，這3個物理量都易測量；由于金屬絲的伸長量δ太小，不容易測量，所以在計算金屬絲的彈性模量時，關鍵是如何解決測量金屬絲伸長量δ的問題.那么，在實際實驗中，巧妙地利用了光杠桿裝置來測量伸長量δ[6]，實驗裝置如圖1所示.

圖1 彈性模量實驗圖

在圖1中，光杠桿G前足尖到兩后足尖連線的垂直距離為d1，光杠桿平面鏡到直尺S的距離為d2，加砝碼m前后望遠鏡中直尺的讀數(shù)為A0和Am，根據(jù)幾何關系得到金屬絲伸長量δ為

(7)

根據(jù)式(6)、式(7)以及F=mg，得到伸長法測金屬絲彈性模量E的公式為

(8)

(9)

根據(jù)式(9)，可以推導出彈性模量E的標準不確定度u(E)為

(10)

在此實驗中測量的實驗數(shù)據(jù)，如表2所示.

在表2中一個砝碼質量m=360 g，根據(jù)表2數(shù)據(jù)采用3種方法來計算K.

(1)由逐差法計算求得

2.66×10-3m·kg-1

表2 彈性模量實驗數(shù)據(jù)

(2)由圖像法擬合求得

根據(jù)表2實驗數(shù)據(jù)繪制圖2.

圖2 圖像擬合法

(3)由最小二乘法求得

令xi=mi，yi=Ai進行直線擬合(y=a+bx)，此b值就是K，其中

∑xi∑x2i∑yi∑y2i∑xiyi10.0818.1440.438480.024070.56689

則令

計算得到

1.2.1.1 造口護理指導 示范更換造口袋的護理程序，指導患者造口袋滲漏的判斷和更換時機，觀看造口袋的護理錄像，發(fā)放造口護理手冊;指導患者使用收集袋的排放方法以及到廁所排放的方法，學會不同的排放氣體的方法;指導患者正確清潔造口袋;根據(jù)造口形狀、患者經(jīng)濟條件推薦不同的造口輔助產品，以滿足不同的需求;讓患者了解造口的各種并發(fā)癥的癥狀表現(xiàn)，示范處理方法，做到早期預防，告知患者定期隨診的重要性;指導選擇合適的造口袋，造口袋使用的注意事項，告知患者購買途徑以及造口袋的儲存方法及注意事項。

可以看出，3種方法求得的K誤差不大，再結合其他實驗數(shù)據(jù)l為38.52 cm，d1為9.800 cm，d2為161.50 cm，d為0.642 mm，根據(jù)式(9)求得E為

由測量儀器的誤差限可知

采用的都是均勻分布的B類分量，根據(jù)式(10)，計算得到E的標準不確定度

u(E)=0.770×1010N·m-2

那么，金屬絲的彈性模量E的標準不確定度報告為

E±u(E)=(14.46±0.77) ×1010N·m-2

2.2 液體表面張力系數(shù)的標準不確定度

實驗裝置如圖3所示.

圖3 實驗示意圖

在實驗中金屬框P中間拉一金屬細絲ab，將框及細絲浸入水中后慢慢地將其拉出水面，在細絲下面將帶起一水膜，當水膜即將破裂時，則有

F=W+2γl+ldhρg

(11)

在式(11)中，F(xiàn)為向上的拉力，W為框所受重力和浮力之差，l為細金屬絲的長度，d為金屬絲的直徑即水膜的厚度，h為水膜被拉斷前的高度，ρ為水的密度，g為重力加速度.ldhρg為水膜的重量，因為張力大小與長度成正比，同時水膜有前后兩面，所以式(11)中張力為2γl.從式(11)可得水的表面張力系數(shù)表達式為

(12)

在此實驗中，約利秤的勁度系數(shù)k為0.97 N·m-1，測量液體表面張力系數(shù)的其他實驗數(shù)據(jù)如表3所示.

表3 液體表面張力系數(shù)實驗數(shù)據(jù)

在表3中，s1為金屬絲剛好到水面時的位置，s2為水膜剛好破裂時的位置，L0為水膜剛好破裂時B柱上的位置，L1為當轉動E使金屬框緩慢下降到G回到零點(G的橫線、橫線的像和鏡面標線重合位置)時B柱上的位置，可知水膜高度

h=|s2-s1|

F-W=kL=k|L1-L0|

根據(jù)表3數(shù)據(jù)，再代入式(12)可求得水的表面張力系數(shù)為

51.99×10-3N·m-1

由測量儀器的誤差限可知

采用的也都是均勻分布的B類分量，由式(12)可知，γ的標準不確定度為

0.66×10-3N·m-1

那么，水表面張力系數(shù)的標準不確定度報告為

γ±u(γ)=(51.99±0.66)×10-3N·m-1

2.3 曲率半徑的標準不確定度

牛頓環(huán)儀是由平凸透鏡L和磨光的平玻璃板P疊合安裝在金屬框架F中構成的，如圖4所示；測量的實驗數(shù)據(jù)[7]，如表4所示.

(a)實驗裝置

(b)牛頓環(huán)干涉圖樣

表4 牛頓環(huán)實驗測量數(shù)據(jù)

根據(jù)表4中數(shù)據(jù)，可以得出

其中，環(huán)數(shù)差m=10，波長λ=589.3 nm，根據(jù)曲率半徑公式

計算得到R為877.1 mm.

下面計算曲率半徑R的標準不確定度u(R)

Δ的A類分量為

Δ的B類分量為

其中Δins=0.01 mm為移測顯微鏡的最小分度值.

合成后

同時，由于視覺疲勞所造成的環(huán)數(shù)誤差u(m)，不妨假設u(m)=0.1；再結合u(λ)=0.30 nm(鈉雙黃線波長λ1=589.0 nm和λ2=589.6 nm導致)，可得合成標準不確定度為

u(R)=

測量結果為

R±u(R)=(877.1±9.9) mm

而廠家提供的凸透鏡曲率半徑標準值為

AR=855.1 mm

根據(jù)測量結果

R=877.1 mm

與其標準值

AR=855.1 mm

可知，兩者之差不超過其標準不確定度u(R)=9.9 mm的3倍[6]，即

|R-AR|≤3u(R)

則可以認為測量結果和標準值在測量誤差范圍內是一致的.

3 分析和討論

金屬絲彈性模量E的不確定度表達式

但用逐差法求解

時，分母體現(xiàn)出來的是16m，所以計算時

而不是

4 結論

在分析金屬絲彈性模量E的標準不確定度u(E)和液體表面張力系數(shù)γ的標準不確定度u(γ)時，僅僅分析計算了B類標準不確定度作為總的不確定度，略去A類標準不確定度，主要是B類分量起主要作用；在分析曲率半徑R的標準不確定度u(R)時，同時考慮了A類不確定度和B類標準不確定度；對于隨機誤差為主的測量情況，可以忽略B類標準不確定度，只分析計算A類標準不確定度作為總的不確定度.