2. Obliczanie niepewnoci pomiarów
2.3. Opracowanie wyników oraz niepewnoci pomiarów prostych
Dysponuj¹c wynikami x1, x2, ..., xn pomiarów wielkoci fizycznej X (jeli po-miary nie s¹ obarczone b³êdami systematycznymi), obliczamy redni¹ arytmetycz-n¹ x, która stanowi oszacowanie wartoci rzeczywistej µrz
. 1 n x x n i i
∑
= = (2.55)Jako niepewnoæ oceny x wielkoci µrz przyjmujemy odchylenie standardowe redniej arytmetycznej
( )
( )
1 , 1 2 − − = =∑
= n n x x n s s n i i x x (2.56)które przypominamy nosi nazwê niepewnoci standardowej. Wynik obliczeñ zapisujemy w postaci:
.
x
s x
x= ± (2.57)
Rys. 2.12. Histogram czêstoci Ok zapisanych w tabeli 2.9 oraz wykres gêstoci rozk³adu normalnego N(10,322, 10,322)
Ponownie przypomnijmy, ¿e wynik pomiaru jest wielkoci¹ mianowan¹, dlate-go obok wartoci liczbowej nale¿y podaæ jednostkê.
Wygodne jest wprowadzenie niepewnoci wzglêdnej
, x sx
w=
δ (2.58)
która jest wielkoci¹ bezwymiarow¹ podawan¹ najczêciej w procentach:
%. 100 %= ⋅ x sx δ (2.59)
2.3.1. Obliczanie niepewnoci w przypadku ma³ej liczby pomiarów
za pomoc¹ d³ugoci przedzia³ów ufnoci
Jeli seria z³o¿ona jest z niewielkiej liczby pomiarów, prawdopodobieñstwo, ¿e przedzia³y
(x−ks
x,x+ks
x)
, gdzie k =1, 2, 3, obejm¹ wartoæ rzeczywist¹ µrz jest mniejsze ni¿ dla du¿ej liczby pomiarów (patrz podrozdzia³ 2.2.9). Aby otrzymaæ ta-kie same prawdopodobieñstwa dla ma³ej liczby pomiarów, jak dla du¿ej, nale¿y sx pomno¿yæ przez pewien wspó³czynnik t(n,p), zale¿ny od liczby pomiarów n i ¿¹-danego prawdopodobieñstwa p. Tak skonstruowany przedzia³(x t(n, p)sx, x + t(n, p)sx),
nazywany przedzia³em ufnoci, obejmuje wartoæ rzeczywist¹ µrz z zadanym praw-dopodobieñstwem p. Prawdopodobieñstwo p z jakim przedzia³ ufnoci obejmuje µrz nosi nazwê poziomu ufnoci.
Kilka wartoci wspó³czynników t(n,p) dla ma³ych n oraz trzech prawdopodo-bieñstw p podano w tabeli 2.10. Obszerniejszy zestaw wspó³czynników t(n,p) po-dano w tabeli 1 (Dodatek). W æwiczeniach laboratoryjnych korzystamy najczêciej z dwóch poziomów ufnoci: p = 0,6827 oraz p = 0,9973.
Wspó³czynniki t(n,p) s¹ kwantylami rozk³adu Studenta (pseudonim W.S. Gos-seta (18761937)). Teoria przedzia³ów ufnoci zosta³a stworzona przez polskiego uczonego T. Sp³awê-Neymana (18941981) w latach trzydziestych dwudziestego wieku i nale¿y do wa¿nych osi¹gniêæ nauki wiatowej.
Tabela 2.10. Wartoci wspó³czynników Studenta dla trzech poziomów ufnoci p
n 2 3 4 5 6 7 8 9 10 20 41
p = 0,6827 1,837 1,321 1,197 1,142 1,111 1,091 1,077 1,067 1,059 1,024 1,013 p = 0,9545 13,968 4,527 3,307 2,869 2,649 2,517 2,429 2,366 2,320 2,126 2,064 p = 0,9973 235,777 19,206 9,219 6,620 5,507 4,904 4,530 4,277 4,094 3,400 3,199
Jak ju¿ wiemy, jako miarê niepewnoæ redniej arytmetycznej x mo¿na przyj¹æ jej odchylenie standardowe, czyli sx. Inn¹ miar¹ niepewnoci redniej arytmetycz-nej, któr¹ u¿ywa siê, gdy seria pomiarów jest niedu¿a, jest po³owa przedzia³u ufno-ci, czyli wielkoæ t(n,p)sx. Podaj¹c niepewnoæ redniej arytmetycznej nale¿y za-znaczyæ, w jaki sposób zosta³a ona obliczona, a w przypadku, gdy wyra¿a siê j¹ przez po³owê przedzia³u ufnoci, jaki poziom ufnoci p zosta³ przyjêty. Nastêpny przyk³ad pokazuje jak okreliæ niepewnoæ redniej arytmetycznej w przypadku krótkiej se-rii pomiarów oporu za pomoc¹ mostka.
Przyk³ad 9
Czterokrotnie powtórzono pomiar pewnego oporu R, uzyskuj¹c nastêpuj¹ce wy-niki (w Ω): x1 = 72,3, x2 = 71,9, x3 = 72,0, x4 = 71,8. Nale¿y oszacowaæ mierzony opór R za pomoc¹ redniej arytmetycznej x i podaæ niepewnoæ tej oceny wyra¿on¹ przez po³owê przedzia³u ufnoci t(n,p)sx przyjmuj¹c poziom ufnoci p = 0,6827. Korzystaj¹c ze wzorów (2.6) i (2.10) otrzymujemy x = 72Ω oraz sx = 0,108 Ω, a z tabeli 2.10 odczytujemy t(4, 0,6287) = 1,197.
W takim razie po³owa przedzia³u ufnoci wynosi t(4, 0,6187)sx = 1,197·0,108 Ω = 1,293 Ω.
Poszukiwan¹ niepewnoci¹ oceny x = 78 Ω mierzonego oporu R, odpowiadaj¹-c¹ poziomowi ufnoci p = 0,6827, jest wielkoæ 1,293 Ω.
2.3.2. Okrelanie niepewnoci na podstawie klasy przyrz¹dów
W przypadku pomiarów, których celem jest wyznaczenie zale¿noci miêdzy dwo-ma lub wiêksz¹ liczb¹ wielkoci fizycznych czêsto brak jest czasu na to, by dla ka¿-dego punktu pomiarowego wykonaæ odpowiednio du¿¹ liczbê pomiarów umo¿liwia-j¹c¹ obliczenie wartoci redniej oraz niepewnoci pomiaru poszczególnych wiel-koci. Podobna sytuacja ma miejsce w przypadku pomiarów wykonywanych w wa-runkach dynamicznych (zmieniaj¹cych siê w czasie prowadzenia dowiadczeñ), np. pomiary zale¿noci oporu metalu lub pó³przewodnika od temperatury. Czêsto wska-zania mierników powtarzaj¹ siê i wykonywanie serii pomiarowej jest niemo¿liwe. W takich przypadkach korzystamy z pojêcia klasy przyrz¹du pomiarowego oznacza-nej symbolem kl.Klasa przyrz¹du pomiarowegokl jest to wyra¿ona w procentach wartoæ bezwzglêdna ze stosunku maksymalnego dopuszczalnego na danym zakresie po-miarowym b³êdu pomiaru ∆xmax do zakresu pomiarowego Z:
%. 100 max Z x kl= ∆ (2.60)
Aby obliczyæ maksymaln¹ niepewnoæ (b³¹d) bezwzglêdn¹ pomiaru wykonane-go za pomoc¹ miernika o znanej klasie, nale¿y skorzystaæ z równania
. 100
max kl Z x =m
∆ (2.61)
Klasa miernika jest podawana obok jego skali (dotyczy wiêkszoci mierników analogowych) lub w instrukcji obs³ugi.
Na niepewnoæ pomiaru obliczon¹ na podstawie klasy miernika sk³adaj¹ siê za-równo niepewnoci systematyczne jak i niepewnoci przypadkowe. Nale¿y podkre-liæ, ¿e klasa miernika jest podana dla warunków pomiaru okrelonych w instrukcji obs³ugi, np. dla danego zakresu temperatur lub wilgotnoci, pionowego lub pozio-mego ustawienia miernika, okrelonego oporu doprowadzeñ itp.
Korzystaj¹c z mierników analogowych o zmiennych zakresach nale¿y dobie-raæ zakres pomiarowy tak, aby wskazówka znajdowa³a siê (o ile jest to mo¿li-we) powy¿ej 2/3 skali. Zakresy pomiarowe mierników analogowych s¹ tak dobie-rane, aby spe³nienie tego zalecenia by³o mo¿liwe.
Aby uzasadniæ przedstawione wy¿ej zalecenie, przypuæmy, ¿e mierzymy pr¹d o natê¿eniu I = 5 mA za pomoc¹ miliamperomierza klasy 0,5 na zakresie 100 mA. Maksymalna niepewnoæ pomiaru obliczona na podstawie klasy miernika
mA, 5 , 0 100 mA 100 5 , 0 max = ⋅ = ∆I
a maksymalna niepewnoæ wzglêdna δ%max = 10%. Je¿eli pomiar wykonamy na za-kresie 10 mA, to ∆Imax = 0,05 mA aδ%max = 1%.
Korzystaj¹c z mierników elektrycznych nale¿y oprócz wyników pomiarów, zanotowaæ zakresy, na których wykonano pomiary oraz klasê miernika. Bez tych danych nie bêdziemy w stanie oszacowaæ niepewnoci pomiarów!
2.3.3. Niepewnoci pomiarów mierników cyfrowych
Niepewnoæ pomiaru wykonanego miernikiem cyfrowym sk³ada siê z niepewnoci przetwarzania uk³adów analogowych miernika (np. nieliniowoci wzmacniacza) oraz niepewnoci dyskretyzacji. Aby wyjaniæ pojêcie niepewnoci dyskretyzacji, zwróæmy uwagê na to, ¿e korzystaj¹c z mierników analogowych (wska-zówkowych) jestemy w stanie odczytaæ wynik z wiêksz¹ dok³adnoci¹ ni¿ dok³ad-noæ dzia³ek elementarnych (mo¿emy stwierdziæ, ¿e wskazówka znajduje siê w po-³o¿eniu odpowiadaj¹cym np. 15,5 dzia³ek skali). Odczytuj¹c wynik z wywietlacza miernika cyfrowego nie mamy takiej mo¿liwoci, poniewa¿ istnieje pewien zakres sygna³ów wejciowych, dla których wskazania wywietlacza s¹ takie same. Niepew-noæ wynikaj¹c¹ z dyskretnoci miernika nazywamy niepewnoci¹ dyskretyzacji lub rozdzielczoci¹ miernika oznaczana dalej za pomoc¹ symbolu rozdz.
Niepewnoæ przetwarzania mierników cyfrowych podawana jest w procentach wartoci mierzonej x i nazywana czêsto klas¹ miernika. Klasê miernika cyfrowego bêdziemy oznaczaæ przez kld (w celu odró¿nienia od klasy miernika analogowego). Niepewnoæ pomiaru miernika cyfrowego obliczamy korzystaj¹c ze wzoru
. 100
max kl x rozdz
x = d +
∆ (2.62)
Nale¿y zaznaczyæ, ¿e rozdzielczoæ miernika cyfrowego zale¿y od zakresu po-miarowego.
Dodajmy, ¿e odchylenie standardowe u miernika cyfrowego, wynikaj¹ce z dyskre-tyzacji wynosi . 12 rozdz u= (2.63) Przyk³ad 10
Mierz¹c napiêcie za pomoc¹ woltomierza cyfrowego klasy 0,05 odczytano z wy-wietlacza wynik U = 225,3 V. Obliczyæ maksymaln¹ niepewnoæ ∆Umax pomiaru napiêcia. Rozdzielczoæ miernika rozdz = 0,1V. Korzystaj¹c z wzoru (2.62) otrzy-mujemy V. 22 , 0 ) 1 , 0 11265 , 0 ( V 1 , 0 V 3 , 225 100 05 , 0 max = + = + ≈ ∆U
Klasa mierników cyfrowych jest podawana w instrukcji obs³ugi (w laboratorium studenckim jest podana w instrukcji roboczej do danego stanowiska).