277. PTB-Seminar Berechnung der Messunsicherheit – Empfehlungen für die Praxis,
11. und 12. März 2014, PTB Berlin
Version 1.0 © 2014 by ZMK GmbH
1
Die richtige Wahl von Verteilungen
N. Schiering, ZMK GmbH Sachsen-Anhalt
277. PTB-Seminar Berechnung der Messunsicherheit – Empfehlungen für die Praxis,
11. und 12. März 2014, PTB Berlin
Version 1.0 © 2014 by ZMK GmbH
Agenda
� Einleitung� Standardmessunsicherheiten Typ A und Typ B� Normalverteilung� Rechteckverteilung� Dreieckverteilung� Trapezverteilung � U-Verteilung� Poisson-Verteilung� Zusammenfassung
2
277. PTB-Seminar Berechnung der Messunsicherheit – Empfehlungen für die Praxis,
11. und 12. März 2014, PTB Berlin
Version 1.0 © 2014 by ZMK GmbH
Einleitung (1)
� Exakte Kontrolle des Messprozesses und Angabe der Messbedingungen in der messtechnischen Praxis nicht möglich� Messgröße kann daher nicht nur ein einziger Wert zugeordnet werden
� Beschreibung der nicht vollständigen Kenntnisse durch Verteilungen von Werten und Einschätzung deren Variation
� Die Wahrscheinlichkeitsdichteverteilung bestimmt das Gewicht, das einem Wert der Größe X aufgrund der vorhandenen Kenntnisse beigemessen werden kann
� Grundlagen der Wahrscheinlichkeitsrechnung bzw. Statistik
notwendig
3
277. PTB-Seminar Berechnung der Messunsicherheit – Empfehlungen für die Praxis,
11. und 12. März 2014, PTB Berlin
Version 1.0 © 2014 by ZMK GmbH
Einleitung (2)
� Abschätzung der Verteilung aufgrund von speziellen Kenntnissen bzw. grundsätzlichen Überlegungen.
� Der GUM geht von Wahrscheinlichkeitsverteilungen der verträglichen Werte aus. Kann für die Größe Xi auf theoretischer oder empirischer Grundlage eine Wahrscheinlichkeitsverteilung angenommen werden, so sind der Erwartungswert und die Quadratwurzel der Varianz dieser Verteilung als Schätzwert xi
bzw. die ihm beigeordnete Standardmessunsicherheit u(xi ) zu verwenden.
4
( ) ( )
( ) ( ) ( )XVARxudxxgxXVar
XExdxxgxXE
i
i
=−=
===
∫
∫∞
∞−
∞
∞−
)(und)(
und)(
2µ
µ
277. PTB-Seminar Berechnung der Messunsicherheit – Empfehlungen für die Praxis,
11. und 12. März 2014, PTB Berlin
Version 1.0 © 2014 by ZMK GmbH
StandardmessunsicherheitTyp A
� Typ A: Anwendung von Analysemethoden der Statistik für Messreihen unter Wiederholbedingungen zur Ermittlung des Wertes und der ihm beigeordneten Standardmessunsicherheit.
� Für die Normalverteilung gilt: � bester Schätzwert = arithmetischer Mittelwert � beizuordnende Unsicherheit = empirische Standardabweichung
(Wurzel aus Varianz dividiert durch Quadratwurzel aus Anzahl der Messungen):
und
5
∑=
=n
iix
nx
1
1 ( ) ( )∑=
−−
=n
ii xx
nnxu
1
2
1
11
277. PTB-Seminar Berechnung der Messunsicherheit – Empfehlungen für die Praxis,
11. und 12. März 2014, PTB Berlin
Version 1.0 © 2014 by ZMK GmbH
StandardmessunsicherheitTyp B
� Typ B: Ermittlung des Wertes und der ihm beigeordneten Standardmessunsicherheit beruht auf anderen Erkenntnissen und kann aus folgenden Informationen eingeschätzt werden:
� Daten und Ergebnisse aus vorangegangenen Messungen � allgemeine Kenntnisse und Erfahrungen über die
Eigenschaften und das Verhalten von Messinstrumenten und Materialien
� Herstellerangaben � Kalibrierscheine oder anderen Zertifikate � Referenzdaten / Literaturwerte
6
277. PTB-Seminar Berechnung der Messunsicherheit – Empfehlungen für die Praxis,
11. und 12. März 2014, PTB Berlin
Version 1.0 © 2014 by ZMK GmbH
Typ A versus Typ B
� Verwendung nach Methode B nur möglich, wenn ausreichende Erfahrungen und allgemeine Kenntnisse vorhanden sind. Diese Fertigkeit wird in der messtechnischen Praxis erlernt.
� Eine gut fundierte Ermittlung der Standardmessunsicherheit nach Methode B ist genau so verlässlich wie eine Ermittlung nach Methode A. Insbesondere, wenn die Ermittlung nach Methode A nur auf einer verhältnismäßig geringen Anzahl von statistisch unabhängigen Beobachtungen beruht.
� Typ A: Angabe der Standardunsicherheit bereits als Standardabweichung des Mittelwertes.
� Typ B: Berechnung der Standardunsicherheit aller Eingangs-größen aus deren angenommenen Wahrscheinlichkeitsverteilungen
7
277. PTB-Seminar Berechnung der Messunsicherheit – Empfehlungen für die Praxis,
11. und 12. März 2014, PTB Berlin
Version 1.0 © 2014 by ZMK GmbH
Normalverteilung
� Verteilungsform ist eine glockenförmige Normalverteilung(Gauß’sche Glockenkurve)
� Verwendung: bester Schätzwert und beigeordnete Unsicherheit, z. B. aus früheren eigenen oder anderen Messungen oder aus Tabellen bekannt. Diese werden unverändert übernommen und man nimmt für die möglichen Werte eine Normalverteilung an.
– , wobei k der Erweiterungsfaktor ist –
8
k
Uxu
xn
xn
ii
=
= ∑=
)(
1
1
277. PTB-Seminar Berechnung der Messunsicherheit – Empfehlungen für die Praxis,
11. und 12. März 2014, PTB Berlin
Version 1.0 © 2014 by ZMK GmbH
9
NormalverteilungGrafische Darstellung
Quelle: QUAM:2000.P1
277. PTB-Seminar Berechnung der Messunsicherheit – Empfehlungen für die Praxis,
11. und 12. März 2014, PTB Berlin
Version 1.0 © 2014 by ZMK GmbH
NormalverteilungBeispiele
� Unsicherheitsangabe in einem Kalibrierschein
� Unsicherheitsangabe in einem Prüfschein
� Unsicherheitsangabe in einem Eichschein
10
277. PTB-Seminar Berechnung der Messunsicherheit – Empfehlungen für die Praxis,
11. und 12. März 2014, PTB Berlin
Version 1.0 © 2014 by ZMK GmbH
Rechteckverteilung
� In vielen Fällen lässt sich für den Wert einer Größe nur die Ober-und Untergrenze a+ und a- angeben, wobei alle Werte innerhalb der Grenzen als gleich wahrscheinlich angesehen werden können.
� Beschreibung dieses Sachverhaltes durch eine rechteckförmigeWahrscheinlichkeitsverteilung und es gilt:
–– , wobei a Halbweite der Verteilung
11
( )
( )332
1)(
2
1
aaaxu
aax
=−=
+=
−+
+−
277. PTB-Seminar Berechnung der Messunsicherheit – Empfehlungen für die Praxis,
11. und 12. März 2014, PTB Berlin
Version 1.0 © 2014 by ZMK GmbH
12
RechteckverteilungGrafische Darstellung
Quelle: QUAM:2000.P1
277. PTB-Seminar Berechnung der Messunsicherheit – Empfehlungen für die Praxis,
11. und 12. März 2014, PTB Berlin
Version 1.0 © 2014 by ZMK GmbH
RechteckverteilungBeispiele (1)
� Die Rechteckverteilung wird beispielsweise angenommen:
� wenn in Kalibrier-, Prüf- und Eichscheinen oder Herstellerdokumentationen Fehlergrenzen (MPE) angegeben sind
� wenn Temperaturgrenzen gegeben sind� zur Beschreibung der Auflösung einer Skala oder
digitalen Anzeige� oder um eine mögliche zeitliche Drift zu
berücksichtigen
13
277. PTB-Seminar Berechnung der Messunsicherheit – Empfehlungen für die Praxis,
11. und 12. März 2014, PTB Berlin
Version 1.0 © 2014 by ZMK GmbH
RechteckverteilungBeispiele (2)
� Herstellerangaben in Bedienanleitungen:
� Spezifikation � Auflösung� Drift
14
Quelle: www.fluke.com
277. PTB-Seminar Berechnung der Messunsicherheit – Empfehlungen für die Praxis,
11. und 12. März 2014, PTB Berlin
Version 1.0 © 2014 by ZMK GmbH
Dreieckverteilung
� Liegen die Werte mit größerer Wahrscheinlichkeit in der Mitte des Bereiches, dann ist die Annahme einer dreieckförmigen Verteilung sinnvoll:
wobei a Halbweite der Verteilung
� Mathematisch ergibt sich die Dreieckverteilung aus der Summe oder Differenz zweier Rechteckverteilungen mit gleicher Halbweite.
15
( )
( )662
1)(
2
1
aaaxu
aax
=−=
+=
−+
+−
277. PTB-Seminar Berechnung der Messunsicherheit – Empfehlungen für die Praxis,
11. und 12. März 2014, PTB Berlin
Version 1.0 © 2014 by ZMK GmbH
16
DreieckverteilungGrafische Darstellung
T
Quelle: QUAM:2000.P1
277. PTB-Seminar Berechnung der Messunsicherheit – Empfehlungen für die Praxis,
11. und 12. März 2014, PTB Berlin
Version 1.0 © 2014 by ZMK GmbH
DreieckverteilungBeispiel
� Schwankung des Luftdruckes bei der Kalibrierung von Kolbenhubpipettten mit Luftpolster
� Interpolationsabweichung bei der statischen Kalibrierung von anzeigenden Drehmomentschlüsseln
17
Quelle: DAkkS-DKD-R 3-7
277. PTB-Seminar Berechnung der Messunsicherheit – Empfehlungen für die Praxis,
11. und 12. März 2014, PTB Berlin
Version 1.0 © 2014 by ZMK GmbH
Trapezverteilung
� Die trapezförmige Verteilung ist ein weiterer Sonderfall, die sich aus der Summe oder der Differenz zweier Rechteckverteilungen mit unterschiedlicher Halbweite ergibt. Sie kommt in der Praxis selten vor.
wobei a Halbweite der Verteilung und 0 ≤ b ≤ 1
18
( )
216
)(
2
1
β+⋅=
+= +−
axu
aax
Quelle: QZ-onlinde.de
277. PTB-Seminar Berechnung der Messunsicherheit – Empfehlungen für die Praxis,
11. und 12. März 2014, PTB Berlin
Version 1.0 © 2014 by ZMK GmbH
U-Verteilung
� Wenn bei harmonischen Schwingungen nichts über die Phase bekannt ist, dann ist für die Phase eine Rechteckverteilung anzunehmen. Mathematisch folgt daraus für die Amplitude eine U-förmige Verteilung.
wobei a hier gleich der Amplitude A0 ist, z.B. A(x) = A0 sin ωx
19
( )
2)(
2
1
axu
aax
=
+= +−
Quelle: Meßunsicherheit, Wolfgang Kessel,
277. PTB-Seminar Berechnung der Messunsicherheit – Empfehlungen für die Praxis,
11. und 12. März 2014, PTB Berlin
Version 1.0 © 2014 by ZMK GmbH
U-VerteilungBeispiele
� Bei Messungen im Bereich der elektromagnetischen Verträglichkeit (EMV)
� Bei Messungen im Bereich der Hochfrequenztechnik
� Bei der Kalibrierung von Beschleunigungsaufnehmern
� Einfluss der Hebelarmlänge bei der statischen Kalibrierung von Kalibriereinrichtungen für Drehmomentschraubwerkzeuge
20
Quelle: DAkkS-DKD-R 3-8
277. PTB-Seminar Berechnung der Messunsicherheit – Empfehlungen für die Praxis,
11. und 12. März 2014, PTB Berlin
Version 1.0 © 2014 by ZMK GmbH
Poissonverteilung
� Als letztes Beispiel sei die Poissonverteilung genannt. Für diese gilt:
wobei N .. Anzahl der Ereignisse.
� Es ist sinnvoll die relative Unsicherheit zu betrachten. Sie ist gegeben durch 1/√N.
� Die Poissonverteilung ist nicht symmetrisch, aber für große Werte von N geht sie in eine Normalverteilung über.
21
Nxu =)(
277. PTB-Seminar Berechnung der Messunsicherheit – Empfehlungen für die Praxis,
11. und 12. März 2014, PTB Berlin
Version 1.0 © 2014 by ZMK GmbH
PoissonverteilungBeispiele
� Eine Poissonverteilung wird z.B. angenommen für die Anzahl
� von Zerfällen radioaktiver Substanzen in einem Zeitintervall
� von Ereignissen in einem Teilchendetektor
22
277. PTB-Seminar Berechnung der Messunsicherheit – Empfehlungen für die Praxis,
11. und 12. März 2014, PTB Berlin
Version 1.0 © 2014 by ZMK GmbH
Zusammenfassung (1)
� Erfahrungen und allgemeine Kenntnisse aus der messtechnischen Praxis müssen vorhanden sein, um anhand der verfügbaren Informationen der Eingangsgröße die richtige Verteilungsfunktion auswählen zu können.
� Manchmal wird mehr als eine Angabe berücksichtigt. Dann ist es nicht eindeutig, welche Wahrscheinlichkeitsverteilung man wählen soll.
� Im GUM Supplement 1, Tabelle 1 werden weitere Beispiele gegeben, wie aus vorhandenen Informationen einer Typ B –Eingangsgröße auf deren Wahrscheinlichkeitsdichtefunktion (PDF) geschlossen werden kann.
23
277. PTB-Seminar Berechnung der Messunsicherheit – Empfehlungen für die Praxis,
11. und 12. März 2014, PTB Berlin
Version 1.0 © 2014 by ZMK GmbH
24
Zusammenfassung (2)
z.B. Spezifikationsangaben, Auflösung
z.B. Messgerät mit Kalibrierschein Quelle: JCGM 102:2011
277. PTB-Seminar Berechnung der Messunsicherheit – Empfehlungen für die Praxis,
11. und 12. März 2014, PTB Berlin
Version 1.0 © 2014 by ZMK GmbH
Zusammenfassung- Übersicht -
25
Verteilungsform Standard-messunsicherheit
Normal oder
Rechteck
Dreieck
Trapez
U-förmig
Poisson
6)(
axu =
3)(
axu =
2)(
axu =
Nxu =)(
216
)( β+⋅= axu
σ=)(xuk
Uxu =)(
277. PTB-Seminar Berechnung der Messunsicherheit – Empfehlungen für die Praxis,
11. und 12. März 2014, PTB Berlin
Version 1.0 © 2014 by ZMK GmbH
2626
Danke für Ihre Aufmerksamkeit !Zentrum für Messen und Kalibrieren GmbH
Sachsen-Anhalt
www.zmk-wolfen.de