ANALIZA NIEPEWNOŚCI

POMIAROWYCH

1. Pomiary wielkości fizycznych

2. Błędy i niepewności pomiarowe

3. Metody określania niepewności pomiarowych

4. Zapis wyników pomiaru

5. Przykład opracowania wyników doświadczenia

6. Dodatek: Zestawienie najwaŜniejszych elementów Międzynarodowej Normy Oceny Niepewności Pomiarowej

1

1. Pomiary wielkości fizycznych

Pomiar wielkości fizycznej polega na porównaniu jej z wielkością tego samego rodzaju przyjętą za jednostkę. Zatem liczba otrzymana jako wynik pomiaru zaleŜy od wyboru jednostki (przykład: pomiar długości w cm, m, ft, in itp.). Wynik pomiaru musi więc zawsze składać się z dwóch części: wartości liczbowej oraz jednostki.

Pomiary wielkości fizycznych dzielimy na bezpośrednie i pośrednie. Pomiary

bezpośrednie są najprostsze – polegają wprost na porównaniu danej wielkości z odpowiednią miarą wzorcową np. pomiar wymiarów ciała za pomocą linijki, suwmiarki, śruby mikrometrycznej itp., pomiar czasu trwania jakiegoś procesu przy uŜyciu stopera, pomiar natęŜenia prądu amperomierzem. W przypadku pomiarów pośrednich wartość badanej wielkości wyznaczana jest na podstawie pomiarów bezpośrednich innych wielkości fizycznych, które są z nią związane znanym nam prawem fizycznym. Na przykład – chcemy

wyznaczyc wartość przyspieszenia ziemskiego na podstawie okresu drgań wahadła

matematycznego. Jak wiadomo okres drgań wahadła opisuje wzór: T = 2π l / g , stąd 2

4π l

g =

. W celu wyznaczenia wartości g musimy zatem dokonać pomiarów

2

T

(bezpośrednich) okresu drgań wahadła ( T) oraz długości nici ( L). Innym przykładem jest wyznaczanie natęŜenia prądu elektrycznego na podstawie pomiarów spadku napięcia na oporniku wzorcowym oraz prawa Ohma I = U / R . Widzimy, Ŝe w zaleŜności od wyboru metody pomiarowej, wartości niektórych wielkości fizycznych mogą być wyznaczane

zarówno drogą pomiarów bezpośrednich, jak i pośrednich.

2. Błędy i niepewności pomiarowe

NiezaleŜnie od metody pomiarów nie moŜemy nigdy bezwzględnie dokładnie wyznaczyć rzeczywistej wartości wielkości fizycznej. RóŜnicę pomiędzy wynikiem pomiaru, a

rzeczywistą wartością mierzonej wielkości nazywamy błędem pomiaru. Błędy pomiarów tradycyjnie dzielimy na grube (omyłki), przypadkowe oraz systematyczne.

Błę dy grube powstają zwykle na skutek nieuwagi lub niestaranności obserwatora przy odczytywaniu lub zapisywaniu wyników lub w wyniku nagłej zmiany warunków pomiaru (np. wstrząsy). Jeśli mamy serię pomiarów wyniki obarczone błędem grubym są łatwe do wykrycia i usunięcia.

Błę dy systematyczne wynikają z niedoskonałości przyrządów i metod pomiarowych. MoŜna je redukować stosując bardziej doskonałe i precyzyjne metody i przyrządy, jednak całkowite 2

wyeliminowanie błędów systematycznych jest niemoŜliwe. Rozpoznane błędy systematyczne naleŜy uwzględniać poprzez wprowadzenie odpowiednich poprawek do wyniku, np. kiedy waŜymy na wadze, której wskazanie bez obciąŜenia wynosi m 0 zamiast 0 to m 0 jest błędem systematycznym, który naleŜy odjąć od wyniku waŜenia, innym typowym przykładem jest poprawka na opór wewnętrzny woltomierza przy pomiarze napięcia .

Z błę dami przypadkowymi mamy do czynienia zawsze. Wynikają one z róŜnych przypadkowych i nie dających się uwzględnić czynników (np. wahania temperatury, lub ruch powietrza w pobliŜu przyrządu pomiarowego). Inną przyczyną moŜe być niezgodność przyjętego modelu z obiektem mierzonym – np. gdy mamy zmierzyć średnicę pręta, zakładamy milcząco, Ŝe jest on idealnym walcem, co nie jest prawdą. O istnieniu błędów przypadkowych świadczy niepowtarzalność wyników pomiaru jednej i tej samej wielkości.

Błędy przypadkowe redukuje się poprzez wielokrotne powtarzanie pomiaru – zachodzi wówczas częściowa kompensacja przypadkowych zawyŜających i zaniŜających odchyłek

wyniku.

PoniewaŜ nigdy nie znamy rzeczywistej wartości wielkości mierzonej, więc

posługiwanie się w praktyce pojęciem błędu pomiaru nie jest wygodne. Obecnie przy opracowywaniu wyników pomiarów naleŜy stosować się do zaleceń Międzynarodowej

Normy Oceny Niepewności Pomiaru. Norma ta uzgodniona w 1995 r. i przyjęta ustawowo w

Polsce w 1999 r. znajduje zastosowanie w róŜnych dziedzinach nauki i techniki.

Międzynarodowa Norma zaleca posługiwanie się terminem niepewność pomiarowa zdefiniowanym jako parametr charakteryzujący wątpliwości dotyczące wartości wyniku pomiarowego. Miarą niepewności pomiarowej jest niepewność standardowa, która moŜe być szacowana na 2 sposoby: typu A wykorzystujący analizę statystyczną serii pomiarów oraz typu B oparty na naukowym osądzie obserwatora. Symbolem niepewności standardowej jest u (od ang. uncertainty), który moŜna zapisywać na 3 róŜne sposoby, np. u, u( x) lub u(stęŜenie NaCl). Zaletą tego zapisu jest to, Ŝe informacja o wielkości mierzonej moŜe być wyraŜona słownie, co ułatwia tworzenie dokumentacji pomiaru. NaleŜy jednak pamiętać, Ŝe u nie jest funkcją tylko liczbą!

3. Metody określania niepewności pomiarowych

3.1. Niepewność standardowa pomiarów bezpośrednich

Przypuśćmy, Ŝe wykonaliśmy serię n pomiarów bezpośrednich wielkości fizycznej X

otrzymując wyniki X 1, X 2 ... X n. Jeśli wyniki pomiarów nie są takie same, wówczas za 3

najbardziej zbliŜoną do wartości prawdziwej przyjmujemy średnią arytmetyczną ze wszystkich wyników pomiarów:

n

1

X ≈ X =

∑ X

(1)

i

n i=1

Stwierdzenie to jest tym bardziej słuszne im większa jest liczba przeprowadzonych pomiarów (dla n → ∞ , X → X ). W celu określenia niepewności standardowej posługujemy się w tym wypadku sposobem typu A, czyli korzystamy ze wzoru na odchylenie standardowe średniej

∑ n( Xi − X )2

u( X )

2

i 1

= s =

=

(2)

X

n( n − )

1

Jeśli natomiast wyniki pomiarów nie wykazują rozrzutu, czyli X

, lub teŜ gdy

1 = X

= ...

2

= X n

istnieje tylko jeden wynik pomiaru, wówczas niepewność standardową szacujemy sposobem

typu B. MoŜna np. wykorzystać informację o niepewności maksymalnej ∆ X określonej przez producenta przyrządu pomiarowego, jeśli nie mamy innych dodatkowych informacji,

wówczas niepewność standardową obliczamy ze wzoru

X

∆

u( X ) =

(3)

3

Dla prostych przyrządów (tj. linijka, śruba mikrometryczna czy termometr) jako ∆ X moŜna przyjąć działkę elementarną przyrządu. W elektronicznych przyrządach cyfrowych

niepewność maksymalna podawana jest przez producenta w instrukcji obsługi i jest zwykle kilkakrotnie większa od działki elementarnej. Najczęściej zaleŜy ona od wielkości mierzonej X i zakresu na którym mierzymy Z:

X

∆ = c X + c Z

1

2

Gdy występują oba typy niepewności (tzn. zarówno rozrzut wyników jak i

niepewność wzorcowania) i Ŝadna z nich nie moŜe być zaniedbana (tzn. obie są tego samego rzędu), wówczas niepewność standardową (całkowitą) obliczamy ze wzoru

∆ X

2

( )2

u( X ) = s +

.

(4)

X

3

3.2. Niepewność standardowa pomiarów pośrednich – niepewność złoŜona ( uc) W przypadku pomiarów pośrednich wielkość mierzoną Y obliczamy korzystając ze związku funkcyjnego, który moŜna zapisać w ogólnej postaci: Y = f ( X , X ,..., X ) , gdzie symbolami 1

2

k

X , X ,..., X oznaczamy k wielkości fizycznych mierzonych bezpośrednio. Zakładamy, Ŝe 1

2

k

4

znane są wyniki pomiarów tych wielkości X 1, X 2 ,..., X k oraz ich niepewności standardowe u( X ), u( X ),..., u( X ) . Wynik (końcowy) pomiaru oblicza się wówczas ze wzoru: 1

2

k

Y ≈ Y = f ( X 1, X 2 ,..., X k )

W przypadku pomiarów pośrednich nieskorelowanych (tzn. gdy kaŜdą z wielkości

X , X ,..., X mierzy się niezaleŜnie) niepewność złoŜoną wielkości Y szacujemy przy 1

2

k

pomocy przybliŜonego wzoru:

2

k





∂ f

u Y

( )

X 1 , X 2 ,..., X

u 2 X

(5)

c

= ∑ 

(

k )

( j)

1

X

j = 

∂ j



3.3. Niepewność rozszerzona

Niepewność standardowa całkowicie i jednoznacznie określa wartość wyniku, jednak do wnioskowania o zgodności wyniku pomiaru z innymi rezultatami (np. z wartością

tabelaryczną) oraz dla celów komercyjnych i do ustalania norm przemysłowych, zdrowia, bezpieczeństwa itp. Międzynarodowa Norma wprowadza pojęcie niepewności rozszerzonej oznaczanej symbolem U (dla pomiarów bezpośrednich), lub Uc (dla pomiarów pośrednich).

Wartość niepewności rozszerzonej oblicza się ze wzoru

U ( X ) = ku( X ) lub U ( X ) = ku ( X ) (6)

c

c

Liczba k, zwana współczynnikiem rozszerzenia, jest umownie przyjętą liczbą wybraną tak, aby w przedziale X ± U ( X ) znalazła się większość wyników pomiaru potrzebna dla danych zastosowań. Wartość współczynnika rozszerzenia mieści się najczęściej w przedziale 2-3. W

większości zastosowań zaleca się przyjmowanie umownej wartości k = 2 .

4. Zapis wyników pomiaru

Wyniki pomiaru zapisujemy zawsze łącznie z niepewnością i jednostką. Niepewność podajemy zawsze z dokładnością do dwu cyfr, zaś liczbę cyfr znaczących wyniku dobieramy tak, aby ostatnia cyfra rezultatu i niepewności naleŜały do tego samego rzędu. Dla niepewności standardowych zalecany jest zapis z uŜyciem nawiasów, zaś dla niepewności rozszerzonej stosowany jest zapis z uŜyciem symbolu ± .

Przykłady zapisu

Dobrze:

Niepewność standardowa:

m = 100,0214 g, u( m) = 3,5 mg

5

m = 100,0214(35) g

m = 100,0214(0,0035) g

Niepewność rozszerzona:

m = 100,0214 g, U ( m) = 0,0070 g

m = (100,0214 ± ,

0 0070 ) g

Źle:

m = 100,0214 g – nie podano niepewności,

m = 100,021(0,0035) g – ostatnie cyfry rezultatu i niepewności nie są tego samego rzędu, m = 100,021 g, u( m) = 3 mg – przy zapisie niepewności podano zbyt mało cyfr, m = 100,02147(0,00352) g - przy zapisie niepewności podano zbyt duŜo cyfr.

5. Przykład opracowania wyników doświadczenia

Celem wyznaczenia przyspieszenia ziemskiego przeprowadzono pomiary czasu spadku ciała

z pewnej wysokości. Wysokość spadku h zmierzono 3-krotnie taśmą mierniczą z podziałką milimetrową uzyskując za kaŜdym razem wynik 1270 mm. Czas spadku t zmierzono 5 razy otrzymując następujące wyniki (w s) t = 5

,

0 09 , t = 5

,

0 12 , t = 5

,

0 10 , t = 5

,

0 04 ,

1

2

3

4

t = 5

,

0 01. Dokładność czasomierza wynosiła 0,001 s, zaś niepewność systematyczną

5

związaną z wyborem chwili włączenia i wyłączenia oszacowano na 0,01 s. Obliczyć z tych danych przyspieszenie ziemskie i jego niepewność.

2 h

Przyspieszenie ziemskie będziemy obliczać ze wzoru g =

. Wartość g otrzymamy

2

t

wstawiając do powyŜszego równania średnie arytmetyczne wysokości spadku ( h ) oraz czasu spadku ( t ) (wzór (1)). Dla danych z tego przykładu mamy:

1

h = 1270mm = 1,27 m, t = ( 5

,

0 09 + 5

,

0 12 + 5

,

0 10 + 5

,

0 04 + 5

,

0 0 )

1 s = 0,5072 s,

5

2 ⋅ ,

1 27 m

m

stąd

g =

≈ 8

,

9 7

4

2

2

(0,507

2 s

)

s

Aby obliczyć niepewność złoŜoną pomiaru pośredniego g musimy najpierw określić niepewności standardowe pomiaru czasu i wysokości.

6

Oszacowanie niepewności standardowej (bezpośredniego) pomiaru czasu u( t): Ocena typu A: Korzystając ze wzoru (2) oraz z tabeli obliczamy odchylenie standardowe średniej t :

Nr pomiaru t

2

i [s]

t − t [ms]

i

t − t

[ms2]

i

1

0,509

1,8

3,24

2

0,512

4,8

23,04

3

0,510

2,8

7,84

4

0,504

3,2

10,24

5

0,501

6,2

38,44

Suma: 82,80

82 8

, 0 m

s2

2

s =

=

1

,

4 4 m

s

= 2,0 ms

t

5 ⋅ 4

Ocena typu B: MoŜemy przyjąć, Ŝe niepewność maksymalna związana z pomiarem czasu wynika przede wszystkim z niepewności chwili włączenia i wyłączenia, a zatem wynosi

∆ t = ,

0 01 s = 10 ms (zaniedbujemy przy tym 10-krotnie mniejszą niepewność związaną z

∆ t

dokładnością czasomierza). Niepewność standardowa typu B wynosi zatem

= 5,8 ms

3

(wzór (3)). Jak widać w tym wypadku naleŜy uwzględnić oba typy niepewności

standardowych (poniewaŜ są one tego samego rzędu). Ostatecznie więc całkowita niepewność standardowa pomiaru czasu wynosi (wzór (4)):

u( t) = ( 0

,

2 2 + 8

,

5 2 m

)

s 2 ≈ 1

,

6 ms = 0,0061 s.

Końcowy wynik pomiaru czasu moŜna zapisać w postaci: t = 0,5072(0,0061) s.

Oszacowanie niepewności standardowej (bezpośredniego) pomiaru wysokości u( h): PoniewaŜ w tym wypadku nie wystąpił rozrzut wyników więc poprzestaniemy na określeniu

niepewności standardowej typu B. Najmniejsza działka przyrządu pomiarowego wynosi w tym wypadku 1 mm. PoniewaŜ jednak pewien wpływ na wynik pomiaru moŜe mieć równieŜ

sposób ustawienia miarki oraz sposób odczytu, rozsądnie będzie przyjąć, Ŝe niepewność maksymalna tego pomiaru jest większa od działki elementarnej np. dwukrotnie: ∆ h = 2 mm.

Zgodnie ze wzorem (3), niepewność standardowa pomiaru wysokości wynosi zatem:

u( h) = ∆ h / 3 = 1,2 mm = 0,0012 m, a więc h = 1270,0(1,2) mm.

7

Oszacowanie niepewności złoŜonej pomiaru pośredniego uc( g): Korzystamy ze wzoru (5). Obliczmy najpierw pochodne cząstkowe:

g

∂

2

g

∂

4 h

( t, h) =

,

( t, h) =

. Aby niepewność u ( g) wyraŜona była w m/s2, przy

2

c

h

∂

3

t

t

∂

t

podstawianiu danych do wzoru (5) musimy pamiętać o uzgodnieniu jednostek ( t i u( t) naleŜy wyrazić w s, zaś h i u( h) naleŜy wyrazić w m).

2

2

2

2

 2 

 4 



h

2



 4 ⋅ ,12700 m



2

2

u ( g)





c

=

u ( h)





+

u ( t) =



⋅ 0

,

0 012 m

 + 

⋅ 0

,

0 061s



 2 

 3 

 t 

 t 

( 5

,

0 07 s

)2



 ( 5

,

0 07 s

)3



u ( g)

c

≈ (

−

m

m

m

5

7

,

8 ⋅10

+ 0

,

0 57) 2

2

≈

0

,

0 57

≈ ,

0 24

4

4

2

s

s

s

Jak widać, przyczynek do niepewności złoŜonej uc( g) związany z niepewnością pomiaru wysokości okazał się zaniedbywalnie mały.

Obliczenie niepewności rozszerzonej Uc( g):

m

m

Podstawiając dane do wzoru (6) otrzymujemy: U ( g)

.

c

= 2 u ( g)

c

= 2 ⋅ ,

0 2

4

= ,

0 48

2

2

s

s

Ostatecznie końcowy rezultat pomiaru przyspieszenia ziemskiego, który moŜemy

porównywać z wielkością tablicową, wygląda następująco:

g =( 9,87 ± 0,48) m/s2

Literatura

1. A. Zięba, 2001 : Natura rachunku niepewności pomiarowych a jego nowa kodyfikacja.

Postępy fizyki 52, nr 5, s. 238-247

2. H. Szydłowski, 2000: Międzynarodowe normy oceny niepewności pomiarowych.

Postępy fizyki 51, nr 2, s. 92-97

3. Guide to Expression of Uncertainty in Measurement, ISO 1995, Switzerland.

Tłumaczenie: WyraŜanie niepewności pomiaru. Przewodnik (Główny Urząd Miar

Warszawa 1999)

8

Dodatek: Zestawienie najwaŜniejszych elementów Międzynarodowej Normy Oceny Niepewności Pomiarowej

Wielkość

Symbol i sposób obliczania oraz nr wzoru w tekście

Niepewność standardowa:

Podstawa: statystyczna analiza serii pomiarów.

ocena typu A

Dla serii n równowaŜnych pomiarów (wzory (2) i (1)):

(pomiary bezpośrednie)

∑ n ( X i − X )2

n

1

u( X )

2

i 1

= s =

=

, gdzie X ≈ X =

∑ X

X

n( n − )

1

i

n i=1

Niepewność standardowa:

Podstawa: naukowy osąd eksperymentatora.

ocena typu B

X

∆

u( X ) =

(3)

(pomiary bezpośrednie)

3

(gdy znana jest niepewność maksymalna ∆ X)

Niepewność standardowa całkowita

∆ X

2

( )2

u( X ) = s +

(4)

ocena typu A oraz typu B

X

3

(pomiary bezpośrednie)

(gdy niepewności typu A i typu B są tego samego rzędu)

Niepewność złoŜona

Dla wielkości Y = f ( X , X ,..., X ) :

1

2

k

(pomiary pośrednie)

2

k





∂ f

u Y

( )

X 1 , X 2 ,..., X

u 2 X

(5)

c

= ∑ 

(

k )

( j )

1

X

j = 

∂ j



(gdy wszystkie wielkości Xi są nieskorelowane)

Współczynnik rozszerzenia

k ≥ 2

Niepewność rozszerzona

U ( X ) = ku( X ) lub U ( X ) = ku ( X ) (6) c

c

Zalecany

zapis

niepewności standardowa: g = ,

9 781 m/s2, u ( g)

m/s2

c

= 0

,

0 76

(przykład)

g = ,

9 78 (

1 76) m/s2

g = ,

9 78 (

1 ,

0 076) m/s2

rozszerzona: g = 7

,

9 8 m/s2, U ( g)

m/s2

c

= 1,

0 5

g =

,

9

( 78 ± 1

,

0

)

5 m/s2

(obowiązuje zasada podawania 2 cyfr znaczących

niepewności)

9