NABL 141 - Guidelines For Estimation and Expression of Uncertainty in Measurement
NABL 141 - Guidelines For Estimation and Expression of Uncertainty in Measurement
NABL 141 - Guidelines For Estimation and Expression of Uncertainty in Measurement
NABL
NATIONAL ACCREDITATION
BOARD FOR TESTING AND
CALIBRATION LABORATORIES
EXPRESSION OF
UNCERTAINTY IN
MEASUREMENT
ISSUE NO : 02
ISSUE DATE: 02.04.2000
AMENDMENT NO : 03
AMENDMENT DATE: 18.08.2000
GOVT. OF INDIA
MINISTRY OF SCIENCE & TECHNOLOGY
Technology Bhavan, New Mehrauli Road,
New Delhi - 110016
FOREWARD
The expression of Uncertainty in Measurements is an integral component of the
accreditation certificate being issued to the calibration laboratories. Globalization of trade and
technology implies the need for interchangeability of components, which must be produced
with a high degree of exactness in measurement system. This concept is equally true for all
other fundamental units of measurement. The International Bureau of Weights and Measures
(BIPM), in consultation with various international bodies, have arrived at a new ISO standard
on Expression of Uncertainty in Measurements, in 1995.
I am glad to dedicate the document of NABL on Guidelines for Estimation and Expression of
Uncertainty in Measurement to the cause of calibration laboratories in the country. I take this
opportunity to congratulate the scientists who have made handsome contributions in bringing
out this document based on the latest ISO standard.
New Delhi
2nd April, 2000
V. S. Ramamurthy,
Chairman, NABL
and Secretary, DST
AMENDMENT SHEET
Sl
no
1
Page
No.
28
Clause
Date of
Amendment made
Reasons
No.
Amendment
Appen- 18.08.00 Not deleted from Printing
dix B
Type B evaluation Mistake/
in Note: /
APLAC
evaluation
Signature
QM
Sd/-
Signature
Director
Sd/-
29
Appen- 18.08.00
dix B
Sd/-
Sd/-
Interpretation
to For Better clarity
effective degrees of on selection of
freedom is added
effective degrees
of freedom /
APLAC
evaluation
Sd/-
Sd/-
10
Page No: ii
Contents
Sl.
Section
Page No.
Introduction
10
16
18
19
20
21
10
22
11
27
12
32
Page No: ii
1.
Introduction
1.1
Purpose
The purpose of the document is to harmonize procedures for evaluating uncertainty in
measurements and for stating the same in calibration certificates as are being followed
by the NABL with the contemporary international approach. The document will apprise
calibration laboratories of the current requirements for evaluating and reporting
uncertainty and will assist accreditation bodies with a coherent assignment of test
measurement capability to calibration laboratories accredited by them. The document
will also provide broad guidelines to all those who are concerned with measurements
about uncertainty in measurement, estimation and apportionment of uncertainty and
interpretation of uncertainty. In fact, the purpose is to provide guidelines to users about
contemporary requirements for global acceptance of various kinds of measurements.
Attempts have been made to make the provisions of this document easy to understand
and ready for implementation. The present document will replace NABLs document
141 (1992).
1.2
Scope
Provisions of this document apply to measurements of all sorts as are carried out in
calibration laboratories. For specialized measurements, these may have to be
supplemented by more specific details and, in some cases, appropriately modified
forms of the concerned formulae. Measurements which can be treated as outputs of
several correlated inputs have been excluded from the scope of this document. The
document covers the following topics:
-
Page No: 1/ 70
1.3
Normative References :
This document is based primarily on the Guide to the expression of uncertainty in
measurement (1993) jointly prepared by BIPM, IEC, ISO and OIML for definition of
various terms and phrases. One should refer to ISO 3534-I (1993) part I probability
and general statistical terms.
1. Guidelines for estimation and statement of overall uncertainty in measurement
results, NABL 141, Department of Science and Technology, New Delhi (India),
(1992).
2. Guide to the expression of uncertainty in measurement, International Bureau of
Weights and Measures (BIPM), International Organization for Standardization (ISO)
et. al., Switzerland, 1995.
3. International vocabulary of basic and general terms in metrology, International
Bureau of Weights and Measures (BIPM), International Organization for
Standardization (ISO) et. al. ., Switzerland , 1993.
4. Expression of the uncertainty of measurement in calibration, European Cooperation
for Accreditation of laboratories (EAL R-2), 1997
5. Guidelines on the evaluation and expression of the measurement uncertainty,
Singapore Institute of Standards and Industrial Research, Singapore 1995.
6. International standard ISO 3534 I, statistics vocabulary and symbols Part I.
Probability and general statistical terms, first edition, International Organization for
Standardization (ISO) , Switzerland ,1993.
Page No: 2/ 70
2.
2.1
Concept
2.1.1
Quality of measurements has assumed great significance in view of the fact that
measurements (in a broad sense) provide the very basis of all control actions.
Incidentally, the word measurement should be understood to mean both a process and
the output of that process.
2.1.2
It is widely recognized that the true value of a measurand (or a duly specified quantity to
be measured) is indeterminate, except when known in terms of theory. What we obtain
from the concerned measurement process is at best an estimate of or approximation to
the true value. Even when appropriate corrections for known or suspected components
of error have been applied, there still remains an uncertainty, that is, a doubt about how
well the result of measurement represents the true value of the quantity being
measured.
2.1.3
2.1.4
2.2
Source
2.2.1
2.2.2
Random errors presumably arise from unpredictable and spatial variations of influence
quantities, for example:
-
Page No: 3/ 70
2.2.3
Various other kinds of errors, recognized as systematic, are also observed. Some
common type of these errors are due to:
-
2.2.4
It should be pointed out that errors, which can be recognized as systematic and can be
isolated in one case, may simply pass of as random in another case.
2.3
Measures
2.3.1
(2.1)
The model function f represents the procedure of the measurement and the method of
evaluation. It describes how values of the output quantity Y are obtained from values of
the input quantities Xi.
2.3.2
(2.2)
It is understood that the input values are best estimates that have been corrected for all
effects significant for the model. If not, necessary corrections have been introduced as
separate input quantities.
2.3.3
2.3.4
The standard uncertainty of the result of a measurement, when that result is obtained
from the values of a number of other quantities is termed combined standard
uncertainty.
2.3.5
Page No: 4/ 70
3.
The guide explains explicitly a large number of metrological terms which are used in practice.
A few terms of general interest have been taken from the International Vocabulary of Basic
and General terms in Metrology and EAL document [3-4]. To facilitate the reader, various
terms and phrases are arranged in alphabetical order
accepted reference value
a value that serves as an agreed upon reference for comparison.
accuracy of measurement
the closeness of agreement between a test result and the accepted reference value
arithmetic mean
The sum of values divided by the number of values
combined standard uncertainty (uc)
standard uncertainty of the result of a measurement when that result is obtained from the
values of a number of other quantities, equal to the positive square root of a sum of terms, the
terms being the variances or covariances of these other quantities weighted according to how
the measurement result varies with changes in these quantities
conventional true value (of a quantity)
a value of a quantity which for a given purpose, may be substituted for the true value
correction
value added algebraically to the uncorrected result of a measurement to compensate for
systematic error
correction factor
numerical factor by which the uncorrected result of a measurement is multiplied to
compensate for a systematic error
correlation
the relationship between two or several random variables within a distribution of two or
more random variables
correlation coefficient
the ratio of the covariance of two random variables to the product of their standard deviations.
Page No: 5/ 70
covariance
The sum of the products of the deviations of xi and yi from their respective averages divided by
one less than the number of observed pairs:
)(
1 N
s xy =
x i x yi y
n 1 i =1
(3.1)
Page No: 6/ 70
expectation
the expectation of a function g(z) over a probability density function p(z) of the random
variables z is defined by
E[g(z)] = g(z)p(z)dz
(3.2)
the expectation of the random variable z, denoted by z and which is also termed as the
expected value or the mean of z. It is estimated statistically by z , the arithmetic mean or
average of n independent observations zi of the random variable z, the probability density
function of which is p(z)
z=
1 n
zi
n i =1
(3.3)
s (q j ) =
(q
n
j =1
n 1
(3.4)
qj being the result of the jth measurement and q being the arithmetic mean of the n results
considered.
measurand
a quantity subject to measurement
probability distribution
a function giving the probability that a random variable takes any given value or belongs to a
given set of values
probability density function
the derivative (when it exits) of the distribution function :
f(x) = dF(x) /dx
(3.5)
(3.6)
probability function
a function giving for every value x, the probability that the random variable X takes value x :
F(x) = Pr(X = x)
(3.7)
Page No: 7/ 70
random error
result of a measurement minus the mean that would result from an infinite number of
measurements of the same measurand carried out under repeatability conditions
Notes :
1.
2.
Page No: 8/ 70
standard uncertainty
uncertainty of the result of a measurement expressed as a standard deviation
systematic error
mean that would result from an infinite number of measurements of the same measurand
carried out under repeatability conditions minus acceptance reference value of the measurand.
Note:
Systematic error is equal to error minus random error
true value (of a quantity)
the value which characterized a quantity perfectly defined in the conditions which exist when
that quantity is considered
Note:
The true value is a theoretical concept, and, in general, can not be known exactly
Type A evaluation (of uncertainty)
Method of evaluation of uncertainty by the statistical analysis of series of observations
Type B evaluation (of uncertainty)
Method of evaluation of uncertainty by means other than the statistical analysis of series of
observations.
uncertainty (in measurement)
parameter, associated with the result of a measurement that characterizes the dispersion of the
values that could reasonably be attributed to the measurand.
variance
A measure of dispersion, which is the sum of the squared deviations of observations from their
average divided by one less than the number of observations.
Page No: 9/ 70
4.
4.1
General considerations
4.1.1
Examples:
Case I : Digital multimeter (DMM)
Let us consider, an experiment in which a high accuracy reference standard e.g. a 6
digit stable meter calibrator is used to calibrate a device of much lower accuracy
like 4 digit DMM . The readings of the test DMM may remain unchanged or
undergo flicker 1 count due to its digitizing process. In this case, the Type A
evaluation of the uncertainty may be taken to be negligible, and the uncertainty on
account of repeatable observations can be treated as Type B on the basis of the
resolution error of the test DMM.
4.2.1
4.2.2
Let us denote by Q the repeatedly measured input quantity Xi. With n statistically
independent observations (n > 1), the estimate of Q is q, the arithmetic mean of the
individual observed values qj (j = 1, 2,.n).
q=
1 n
qj
n j =1
(4.1)
s 2 (q ) =
1 n
qj q
n 1 j =1
(4.2)
The positive square root of s2 (q) is termed experimental standard deviation. The best
estimate of the variance of the arithmetic mean q is given by
s 2 (q )
s q =
n
2
()
(4.3)
Table 4.1: Data for calculation of mean and standard deviation of temperature:
Observation
numbers
Temperature 0C
(t
t 10 2
0
(t
t 10 4
( C)
0
1
2
3
4
5
6
7
8
9
10
90.68
90.83
90.79
90.64
90.63
90.94
90.60
90.68
90.76
90.65
-4
11
7
-8
-9
22
-12
-4
4
-7
16
121
49
64
81
484
144
16
16
49
Total
907.2
1040
The positive square root of s2 (q) is termed as estimated standard error of the mean.
The standard uncertainty u (q) associated with the input estimate q is the standard
error.
() ()
uq =sq
(4.4)
4.2.4
()
s2 q =
s 2p
(4.5)
The standard uncertainty is deduced from the value given by Eq. (4.4)
Example:
Table (4.1) is shown the data from a temperature measurement. We now
estimate different parameters as follows:
Mean Temperature:
t
j
j =1
= 90.72
t=
n
(4.6)
t = 90.72 C
(4.7)
Standard Deviation:
s (t ) =
1 n
(t j t
n 1 j =1
1
1040 10 4 = 10.75 10 2
9
(4.8)
()
st =
s (t ) 10.75 10
=
n
10
2
= 3.40 10 2 C
(4.9)
Standard uncertainty:
()
u t = 3.40 10 2 C
(4.10)
Degrees of freedom ()
= n 1 = 10 1 = 9
(4.11)
4.3
4.3.1
4.3.2
manufacturers specifications ;
The proper use of the available information for a Type B evaluation of standard
uncertainty of measurement calls for insight based on experience and general
knowledge. It is a skill that can be learned with practice. A well-based Type B
evaluation of standard uncertainty can be as reliable as a Type A evaluation of standard
uncertainty, especially in a measurement situation where a Type A evaluation is based
only on a comparatively small number of statistically independent observations. The
following cases must be discerned:
(a)
When only a single value is known for the quantity Xi, e.g. a single measured
value, a resultant value of a previous measurement, a reference value from the
literature, or a correction value, this value will be used for xi. The standard
uncertainty u (xi) associated with xi is to be adopted where it is given. Otherwise
it has to be calculated from unequivocal uncertainty data. If data of this kind are
not available, the uncertainty has to be evaluated on the basis of experience
taken as it may have been stated (often in terms of an interval corresponding to
expanded uncertainty).
(b)
Examples:
In cases, where the uncertainty is quoted to be particular multiple of standard deviation (), the
multiple becomes the specific factor (see Appendix A).
Case I:
A calibration certificate states that the mass of a given body of 10 kg is 10.000650 kg. The
uncertainty at 2 (at confidence level of 95 .45 %) is given by 300 mg. In such a case, the
standard uncertainty is simply,
u(m) = 300 / 2 = 150 mg
(4.12)
and estimated variance is
u2(m) = 0.0225 g2
(4.13)
Case II:
Suppose in the above example, the quoted uncertainty defines an interval having a 90% level
of confidence. The standard uncertainty is then
u(m) = 300 / 1.64 = 182.9 mg
(4.14)
Where we have taken 1.64 as the factor corresponding to the above level of confidence,
assuming the normal distribution unless otherwise stated.
Case III:
A calibration certificate states that the resistance of a standard resistor, Rs of nominal
value 10 is 10.000742 129 at 23 0 C and that the quoted uncertainty of 129
defines an interval having a level of confidence of 99%. The standard uncertainty of the
resistor may be taken as
u(Rs ) = 129 / 2.58 = 50
(4.15)
Therefore, in this case, specific factor is 2.58. The corresponding relative standard uncertainty
u(Rs )/ Rs = 5 10-6
(4.16)
(4.17)
Case IV:
A calibration certificate states that the length of a standard slip gauge (SG) of nominal value 50
mm is 50.000002 mm. The uncertainty of this value is 72 nm, at confidence level of 99.7 %
(corresponding to 3 times of standard deviation). The standard uncertainty of the standard slip
gauge is then given by
u(SG) = 72 nm /3 = 24 nm
(4.18)
(c)
If only upper and lower limits a + and a - can be estimated for the value of the
quantity Xi (e.g. manufacturers specifications of a measuring instrument, a
temperature range, a rounding or truncation error resulting from automated data
reduction), a probability distribution with constant probability density between
these limits (rectangular probability distribution) has to be assumed for the
possible variability of the input quantity Xi . According to case (b) above this
leads to
1
(a + + a - )
2
xi =
(4.19)
u 2 (x i ) =
1
(a + a - )2
12
(4.20)
for the square of the standard uncertainty . If the difference between the limiting values
is denoted by 2a , Eq. (4.20) yields
u 2 (x i ) =
1
(a )2
3
(4.21)
Examples:
The specifications of a dial type pressure gauge are as follows :
Range : 0 to 10 bar,
Scale : 1 division = 0.05 bar,
Resolution : division = 0.025 bar
Accuracy : 0.25 % Full Scale Deflection
Assuming that with the above specifications, there is an equal probability of the true
value lying anywhere between the upper (a + ) and lower (a ) limits. Therefore, for
rectangular distribution,
a=
Here ,
and
(a +
a )
2
(4.22)
(4.23)
(4.24)
(4.25)
u=
a
3
0.025
3
= 0.0144 bar
(4.26)
5.
5.1
For uncorrelated input quantities the square of the standard uncertainty associated with
the output estimate y is given by,
n
u 2 ( y ) = ui2 ( y )
(5.1)
i =1
(5.2)
where ci is defined as sensitivity coefficients associated with the input estimate xi i.e.
the partial derivative of the model function f with respect to Xi , evaluated at the input
estimates xi .
c i = (f /xi ) = (f /Xi ) at Xi = xi
(5.3)
5.2
The sensitivity coefficient ci, describes the extent to which the output estimate y is
influenced by variations of the input estimate xi. It can be evaluated from the function f
by Eq. (5.3) or by using numerical methods, i.e. by calculating the change in the output
estimate y due to a change in the input estimate xi of + u(xi) and -u(xi) and taking as the
value of ci the resulting difference in y divided by 2u (xi) . Sometimes it may be more
appropriate to find the change in the output estimate y from an experiment by repeating
the measurement at e.g. xi u (xi).
5.3
p X
i =1
(5.4)
corresponding sum or
y = pi x i
i =1
(5.5)
u 2 (y ) = p i2u 2 (x i )
i =1
(5.6)
5.4
f ( X 1, X 2 ,......., X N ) = c X ipi
(5.7)
i =1
the output
estimates
N
y = c X ipi
(5.8)
i =1
The sensitivity coefficients equal piy/xi in this case and an expression analogous to Eq.
(5.6) is obtained from Eq. (5.1), if relative standard uncertainties w(y) = u(y)/y and w (xi)
= u (xi) / xi are used,
w 2 (y ) =
p w (x )
i =1
2
i
(5.9)
6.
6.1
(6.1)
6.3
The reliability of the standard uncertainty assigned to the output estimate is determined
by its effective degrees of freedom (see Appendix B). However, the reliability criterion is
always met if none of the uncertainty contributions is obtained from a Type A evaluation
based on less than ten repeated observations.
6.4
If one of these conditions (normality or sufficient reliability) is not fulfilled, the standard
coverage factor k = 2 can yield an expanded uncertainty corresponding to a coverage
probability of less than 95 %. In these cases, in order to ensure that a value of the
expanded uncertainty is quoted corresponding to the same coverage probability as in
the normal case, other procedures have to be followed. The use of approximately the
same coverage probability is essential whenever two results of measurement of the
same quantity have to be compared, e.g. when evaluating the results of an
interlaboratory comparison or assessing compliance with a specification.
6.5
Even if a normal distribution can be assumed, it may still occur that the standard
uncertainty associated with the output estimate is of insufficient reliability. If, in this
case, it is not expedient to increase the number n of repeated measurements or to use
a Type B evaluation instead of the Type A evaluation of poor reliability, the method
given in Appendix B should be used.
6.6
For the remaining cases, i.e. all cases where the assumption of a normal distribution
cannot be justified, information on the actual probability distribution of the output
estimate must be used to obtain a value of the coverage factor k that corresponds to a
coverage probability of approximately 95 %.
7.
7.1
7.2
However, in cases where the procedure of Appendix A has been followed, the
additional note should read as follows: The reported expanded uncertainty in
measurement is stated as the standard uncertainty in measurement multiplied by
the coverage factor k which for a t-distribution with eff effective degrees of
freedom corresponds to a coverage probability of approximately 95 %. (See
Appendix B).
7.3
The numerical value of the uncertainty in measurement should be given to at most two
significant figures. The numerical value of the measurement result should in the final
statement normally be rounded to the least significant figure in the value of the
expanded uncertainty assigned to the measurement result. For the process of
rounding, the usual rules for rounding of numbers have to be used. However, if the
rounding brings the numerical value of the uncertainty in measurement down by more
than 5 %, the rounded up value should be used.
8.
8.1
8.2
A formal example of such an arrangement is given as Table (8.1) applicable for the
case of uncorrelated input quantities. The standard uncertainty associated with the
measurement result u(y) given in the bottom right corner of the table is the root sum
square of all the uncertainty contributions in the outer right column. Similarly, eff has to
be evaluated as mentioned in Appendix B.
Estimates Limits
xi
xi
Probability
Distribution
- Type A or
B
Standard
Sensitivity
Uncertainty coefficient
u(xi)
ci
Uncertainty
contribution
ui(y)
Degree
of
freedom
i
X1
X2
X3
x1
x2
x3
x1
x2
x3
-Type A or B
-Type A or B
-Type A or B
u(x1)
u(x2)
u(x3)
c1
c2
c3
u1(y)
u2(y)
u3(y)
1
2
3
XN
Y
xN
y
xN
-Type A or B
u(xN)
cN
uN(y)
uc(y)
N
eff
9.
(9.1)
Step 2
Step 3
Step 4
Step 5
Step 6
For input quantities for which the probability distribution is known or can be
assumed, calculate the expectation and the standard uncertainty u (xi) according
to paragraph 4.3.2 (b). If only upper and lower limits are given or can be
estimated, calculate the standard uncertainty u (xi) in accordance with paragraph
4.3.2(c).
Step 7
Calculate for each input quantity Xi the contribution ui (y) to the uncertainty
associated with the output estimate resulting from the input estimate xi according
to Eqs. (5.2) and (5.3) and sum their squares as described in Eq. (5.1) to obtain
the square of the standard uncertainty u(y) of the measurand.
Step 8
Step 9
()
Appendix A
Probability distribution
A.1
Normal distribution
P (x ) =
1
2
exp (x ) / 2 2 ,
2
-<x <+
(A.1)
Where is the mean and is the standard deviation. Figure (A.1) represents such a
distribution.
A.1.1
Table A.1: Confidence Level and the corresponding Coverage factor (k)
Confidence 68.27 %
level
Coverage 1.000
factor (k)
90 %
95 %
95.45 %
99 %
99.73 %
1.645
1.960
2.000
2.576
3.000
If based on available information, it can be stated that there is 50 % chance that the value of
input quantity Xi lies in the interval between a and a + and also it is assumed that the
distribution of Xi is normal, then the best estimate of Xi is:
(A.2)
If based on available information, it can be stated that there is 68% chance that the value of
input quantity Xi lies in the interval of a and a + ; and also it is assumed that the distribution of
Xi is normal, then the best estimate of Xi is :
x i = a, with u(x i ) = a
A.2
(A.3)
Rectangular distribution
The probability density function p(x) of rectangular distribution is as follows:
P (x ) =
1
, a < x < a + , where a = (a + a ) / 2
2a
(A.4)
E (X i ) = x i = (a + + a ) / 2
(A.5)
Var (X i ) = a 2 / 3 , where a = (a + a ) / 2
(A.6)
- a /3
+ a /3
A.3
E (X i ) = (a + + a ) / 2 ,
(A.7)
Var (X i ) = a 2 (1 + )2 /6
(A.8)
A.3.1
Triangular Distribution
When = 0, the symmetric trapezoidal distribution is reduced to a triangular
distribution. Figure (A.3) shows such a distribution. When the greatest concentration of
the values is at the center of the distribution, then one must use the triangular
distribution.
The expectation of Xi is given as,
E (X i ) = (a + + a ) / 2 ,
(A.9)
Var (X i ) = a 2 / 6
(A.10)
A.4
U-Shaped Distribution
This U-shaped distribution is used in the case of mismatch uncertainty in radio and
microwave frequency power measurements (shown in figure (A.4). At high frequency
the power is delivered from a source to a load, and reflection occurs when the
impedances do not match. The mismatch uncertainty is given by 2s L where s and
L are the reflection coefficients of the source and the load respectively. The standard
uncertainty is computed as:
u2 (xi) = (2 S L)2 / 2
(A.11)
Appendix B
B.2
B.3
Step 1
Step 2
Estimate the effective degree of freedom eff of the standard uncertainty u(y)
associated with the output estimate y from the Welch-Satterthwaite formula.
Step 3
Obtain the coverage factor k from the table of values of student t distribution.
If the value of eff is not an integer, it is truncated to the next lower integer and
the corresponding coverage factor k is obtained from the table.
eff
u4 (y)
= N 4
ui ( y )
i =1
(B.1)
where u i (y) (i = 1,2,3,..N) defined in Eqs. (5.1) and (5.2) , are the contributions to
the standard uncertainty associated with the output estimate y resulting from the
standard uncertainty associated with the input estimate x i which are assumed to be
mutually statistically independent , and the i is the effective degrees of freedom of the
standard uncertainty contributions u i (y).
Note :
The calculation of the degrees of the freedom for Type A and Type B of the evaluation
may be as follows:
Type A Evaluation
For the results of direct measurement (Type A evaluation), the degree of freedom is related to
the number of observations (n) as,
i = n - 1
(B.2)
Type B Evaluation
For this evaluation, when lower and upper limits are known
i
(B.3)
It is suggested that i should always be given when Type A and Type B evaluations of
uncertainty components are documented.
Where high precision measurements are undertaken, the accredited calibration laboratories
shall be required to follow ISO Guide to the expression of uncertainty in Measurement (1995).
Concerned laboratories should refer to Annexure G (with special emphasis on table G-2 )
and Annexure H for related examples.
However, assuming that i is not necessarily unrealistic, since it is a common practice to
chose a- and a+ in such a way that the probability of the quantity lying outside the interval a- to
a+ is extremely small.
Further interpretation on the above is given on page 29 & 30
Where :
u(xi)/ u(xi) is the relative uncertainty in the uncertainty
This is a number less than 1, but may for convenience be thought of as a percentage or a
fraction. The smaller the number, the better defined is the magnitude of the uncertainty.
For example, if relative uncertainty is 10%, i.e.
u(xi)/ u(xi) = 0.1
Then it can be shown that the number of degrees of freedom is 50. For a relative uncertainty of
25 % then = 8 and for relative uncertainty of 50 %, = is only 2.
Rather than become seduced by the elegance of mathematics, it is better to try to determine
the limits more definitely, particularly if the uncertainty is a major one. It is of the interest to note
that equation (1) tells us that when we have made 51 measurements and taken the mean, the
relative uncertainty in the uncertainty of the mean is 10%. This shows that even when many
measurements are taken, the reliability of the uncertainty is not necessarily any better than
when a type B assessment is make. Indeed, it is usually better to rely on prior knowledge
rather than using an uncertainty based on two or three measurements. It also shows why we
restrict uncertainty to two digits. The value is usually not reliable enough to quote to better than
1 % resolution.
National Accreditation Board for Testing and Calibration Laboratories
Doc. No: NABL 141
Issue No: 02
Once the uncertainty components have been combined, it remains to find the number of
degrees of freedom in the combined uncertainty. The degrees of freedom for each component
must also be combined to find the effective number of degrees of freedom to be associated
with the combined uncertainty. This is calculated using the Welch-Satterthwaite equation,
which is:
n
eff = [uc4 (y)/ { ui4 (y)/ i}]
1
. 2
Where:
eff
Table B.1: Student t-distribution for degrees of freedom . The t-distribution for defines an
interval -t p () to + t p () that encompasses the fraction p of the distribution. For p = 68.27 %,
95.45 %, and 99.73 %, k is 1, 2, and 3, respectively.
Degrees
Freedom
Fraction p in percent
()
68.27
90
95
95.45
99
99.73
1.84
6.31
12.71
13.97
63.66
235.80
1.32
2.92
4.30
4.53
9.92
19.21
1.20
2.35
3.18
3.31
5.84
9.22
1.14
2.13
2.78
2.87
4.60
6.62
1.11
2.02
2.57
2.65
4.03
5.51
1.09
1.94
2.45
2.52
3.71
4.90
1.08
1.89
2.36
2.43
3.50
4.53
1.07
1.86
2.31
2.37
3.36
4.28
1.06
1.83
2.26
2.32
3.25
4.09
10
1.05
1.81
2.23
2.28
3.17
3.96
11
1.05
1.80
2.20
2.25
3.11
3.85
12
1.04
1.78
2.18
2.23
3.05
3.76
13
1.04
1.77
2.16
2.21
3.01
3.69
14
1.04
1.76
2.14
2.20
2.98
3.64
15
1.03
1.75
2.13
2.18
2.95
3.59
16
1.03
1.75
2.12
2.17
2.92
3.54
17
1.03
1.74
2.11
2.16
2.90
3.51
18
1.03
1.73
2.10
2.15
2.88
3.48
19
1.03
1.73
2.09
2.14
2.86
3.45
20
1.03
1.72
2.09
2.13
2.85
3.42
25
1.02
1.71
2.06
2.11
2.79
3.33
30
1.02
1.70
2.04
2.09
2.75
3.27
1.000
1.645
1.960
2.000
2.576
3.000
Appendix C
Solved examples showing the applications of the method outlined here to
eight specific problems in different fields
XSTD +
(C.1)
Where YGUT is the micrometer reading [Gauge under test (GUT), XSTD is the gauge block size
and X is the error or the difference between the micrometer reading and gauge block size.
Uncertainty equation
STD
(C.2)
Measured results
Type A evaluation
Five readings are taken and the deviation from the nominal value is as follows.
Mean Deviation:
n
x = ( x j ) /n = 0.0006mm
(C.3)
j =1
Standard deviation :
s( x ) =
1 n
xj x
n 1 j =1
1
(
12 10 7 ) = 3 10 7 mm = 5.47 10 4 mm
4
(C.4)
Observation
numbers
1.
2.
3.
4.
5.
0.001
0.000
0.001
0.000
0.001
x
(mm)
0.0006
1.6
3.6
1.6
3.6
1.6
()
sx =
s 2 ( x ) 5.47 10 4
=
= 2.446 10 4 mm = 0.2446m
n
5
(C.5)
Standard uncertainty:
()
u x = 0.2446m
(C.6)
5 1 = 4
(C.7)
Type B evaluation
The uncertainty quoted in the gauge block calibration certificate is considered to be Type B
uncertainty of normal distribution.
Standard uncertainty (u1) due to the temperature measurement 1 0 C.
Standard thermal expansion coefficient of the gauge block is 11.510-6/ 0C
u 1 = 25 111.5 10-6 mm = 287.5 10-6 mm = 0.2875 m
(C.8 )
Standard uncertainty (u2) due to difference in temperature of micrometer and slip gauge
Assuming the temperature of the slip gauges and micrometer are the same but still it can have
a difference 1 0 C. Hence, again uncertainty component
u2 = 0.2875 m
(C.9)
Standard uncertainty (u3) due to difference in thermal expansion coefficient of the slip
gauge and micrometer
It is assumed that the difference in thermal expansion coefficient of standard slip gauge and
the micrometer screw is amounting to 20 %, hence the uncertainty component
[T = Tc Tref = 30 C],
u3 = 25 3 11.5 10-6 (20 /100) mm = 0.1725 m
(C.10)
(C.11)
(C.12)
u6 =
0.08
m = 0.046m
3
(C.13)
The sensitivity coefficients (ci) are 1 and degree of freedom is i = [Type B components ] in
all six cases .
Degrees of freedom (eff )
eff =
uc (y )
4
n (u c (y ))
4
i =1
Veff =
(0.531)4
(0.2446)4 + (0.165 )4 + (0.165 )4 + (0.099)4
(0.288 )4
(0.288 )4
(0.046)4
(C.14)
0.07950
0.00358
4
89
Combined uncertainty
Combined uncertainty [uc (YGUT)] is ,
0.531 m
(C.15)
(C.16)
Estimate
s
xi
(m)
Limits
x i
(m)
Probability
Distribution
Type A or B
- factor
Standard
Sensitivity Uncertainty Degree
Uncertainty coefficient contribution of
u(x i)
u i (y)
ci
freedom
(m)
(m)
i
u1
0.287
0.143
0.165
1.0
0.165
u2
0.287
0.143
0.165
1.0
0.165
u3
0.172
0.086
0.099
1.0
0.099
u4
0.5
0.25
0.288
1.0
0.288
u5
0.5
0.25
0.288
1.0
0.288
u6
0.08
0.04
0.046
1.0
0.046
Repeatability
u(x)
Rectangular
- Type B
- 3
Rectangular
- Type B
- 3
Rectangular
- Type B
- 3
Rectangular
- Type B
- 3
Rectangular
- Type B
- 3
Rectangular
- Type B
- 3
Normal
- Type A
- 5
0.547
1.0
0.2446
0.531
1.062
uc (YGUT )
Expanded
uncertainty
k=2
eff
= .
(C.17)
(C.18)
Reference standard (ms): The calibration certificate for the reference standard gives a
value of 10,000.005 g with an associated expanded uncertainty of 45 mg (coverage
factor k = 2)
2.
Drift of the value of the standard (mp): The drift of the value of the reference standard
is estimated from the previous calibrations to be zero within 15 mg.
3.
4.
Eccentricity and magnetic effect (mc): The variation of mass due to eccentric load and
magnetic effect is found to be 10 mg.
5.
Air buoyancy (A): The limit of air buoyancy correction is found to be 10 mg.
A reading
(g)
B reading
(g)
B reading
(g)
A Reading
(g)
1.
2.
3.
4.
5.
0.015
0.010
0.020
0.020
0.010
0.020
0.030
0.045
0.040
0.030
0.025
0.020
0.040
0.030
0.020
0.010
0.010
0.015
0.010
0.010
Observed
difference
(g)
0.010
0.015
0.025
0.020
0.015
Type A evaluation
Five observations of the difference in mass between the unknown mass (B) and the standard
(A) are obtained using the substitution method and the ABBA weighing sequences:
Arithmetic mean m = 0.017 g,
pooled estimate of standard deviation sp ( m) = 0.025 g (obtained from prior evaluation)
Standard uncertainty uA = u(m) = s(m ) =
25
mg = 11.18 mg,
5
Degrees of freedom = 5 1 = 4
Type B evaluation
1.
From the calibration certificate of the standard (A), the expanded uncertainty (U) is
certified as 45 mg with a coverage factor k = 2. the standard uncertainty
u (ms ) =
2.
u (m D ) =
3.
45
mg = 22.5mg
2
Assuming a rectangular
15
mg = 8.66mg
3
Due to eccentricity and magnetic effects, the variation of mass is found to be 10 mg.
Assuming a rectangular distribution, the standard uncertainty,
u (m c ) =
10
mg = 5.77mg
3
4.
u ( ) =
uB =
Assuming rectangular
10
mg = 5.77mg
3
(C.19)
Degrees of freedom (i): In all these four cases, the degree of freedom is i =
Degrees of freedom eff
v eff =
(27.8)4
(11.8)4 + (22.5)4 + (8.66)4 + (5.77 )4
4
(
5.77 )
+
598029.44
15623.1
4
~ 153.1
uc = u A2 + u B2 =
(11.18)2 + (25.45)2
= 27.8mg
(C.20)
Expanded uncertainty
U = ku c (m) = 2 27.8 mg = 55.6 mg
Reported result
The measured mass of the nominal 10 kg weight is 10.000025 kg 56 mg. The reported
expanded uncertainty of measurement is stated as the standard uncertainty in measurement
multiplied by the coverage factor k = 2, which for a normal distribution corresponds to a
coverage probability 99.5 %.
Table C.4: Uncertainty Budget:
Source of
Uncertainty
Estimates
xi
Xi
Limits
x i
ms
10.000005 45 mg
kg
m d
15 mg
7.5 mg
m c
10 mg
5 mg
10 mg
5 mg
Repeatabilit
y
Probability
Distribution
-Type A or B
- Factor
Standard
Sensitivity Uncertainty
Uncertainty coefficient contribution
u(x i)
(mg)
ci
Normal
- Type B
- 2
Rectangular
- Type B
- 3
Rectangular
- Type B
- 3
Rectangular
- Type B
- 3
Normal
- Type A
- 5
22.5
1.0
22.5
8.66
1.0
8.66
5.77
1.0
5.77
5.77
1.0
5.77
11.18
11.18
27.8
55.6
uc (m )
Expanded
uncertainty
k=2
u i (y)
(mg)
Degree
of
freedom
vi
Reference level: - The pressure generated at the reference level both at the secondary
standard (SPC) and test gauge (DWT) has been maintained constant.
Mathematical Model
The mathematical relationship can be modeled as:
PDWT = PSPC + P,
Where PDWT is the pressure as measured by the industrial DWT , PSPC
is the average
standard pressure as obtained from the standard gauge PSPC . The average sign indicate the
arithmetic mean of several repetitive measurement under identical condition.
P is the
difference between the two readings of the PDWT and the PSPC . Table (C.4) represents the data
where represents the standard uncertainty in each set of readings at a given pressure.
In the simplest case and also in this limited pressure range up to 100 MPa, it is normally
observed that the P is a linear function of PSPC , which is,
P = Po +
where P0 and S 1
S1
PSPC
(C.22)
SPC
O
S1
where the bracketed quantities are the standard uncertainties due to the repeatability of the
readings and the standard gauge, respectively. It may be mentioned here that the partial
derivatives of PDWT with PSPC and PDWT with Po is 1 , therefore, the sensitivity coefficients
also for both these cases equal to one. However partial derivative of PDWT with S1 is PSPC .
Therefore, we have to take into account the value of
standard uncertainty.
PSPC
Experimental Results
Table (C.4) shows the data as obtained from the experiment. As mentioned, we have taken
five readings (n = 5) of the same pressure point while increasing / decreasing the pressure
cycle. The PSPC is the average of these five readings. P is the error or the difference
between PDWT and PSPC . PW represents Piston and Weight hanger.
Type A evaluation of standard uncertainty
(1)
Repeatability
We have taken 5 repeatable readings at each and every pressure or n= 5. The maximum
standard uncertainty () from the table (C.5) is,
s(qk)
= 0.0026827 MPa
(C.24)
u1(s ) = 0.0026827
5 = 0.001199MPa
(C.25)
Degree of freedom i = 5 - 1 = 4
Table C.5: Increasing and Decreasing Pressure
Weight
used
Nominal
Pressure
PDWT
PW
PW , 7-8
PW,1
PW 1-2
PW , 1-3
PW,1-4
PW, 1-4 ,
6
PW , 1-5
PW , 1-6
PW, 1-9
PW , 1-9
PW , 1-6
PW, 1-5
PW, 1-4 ,
6
PW , 1-4
PW, 1-3
PW , 1-2
PW , 1
PW , 7-8
PW
Reading of SPC
(MPa)
1.0
5.0
10.0
20.0
30.0
40.0
45.0
0.9989
4.9965
9.9931
19.9861
29.9791
39.9694
44.9665
50.0
55.0
60.0
60.0
55.0
50.0
45.0
40.0
30.0
20.0
10.0
5.0
1.0
PSPC
(MPa)
0.9992
4.9967
9.993
19.9865
29.9787
39.9706
44.9654
(MPa)
0.9988
4.9964
9.9933
19.9867
29.9782
39.971
44.9659
0.9978
4.9971
9.9931
19.9871
29.9799
39.9701
44.9671
0.9991
4.9961
9.9935
19.9861
29.9787
39.9691
44.9669
0.9987
4.9965
9.9932
19.9865
29.9789
39.9700
44.9663
(kPa)
0.619
0.422
0.162
0.477
0.728
0.896
0.675
(kPa)
1.24
3.44
6.80
13.50
21.08
29.96
33.64
49.96
54.9548
59.9494
59.9497
54.9568
49.9641
44.9675
49.9611
54.9576
59.9532
59.9552
54.9576
49.9611
44.9654
49.9606
54.9556
59.9498
59.9549
54.9556
49.9606
44.9659
49.9604
54.9561
59.9512
59.9512
54.9561
49.9604
44.9671
49.9595
54.9545
59.9532
59.9532
54.9545
49.9595
44.9669
49.9603
54.9557
59.9513
59.9528
54.9561
49.9611
44.9665
0.626
1.364
1.819
2.682
0.780
0.911
1.000
39.68
44.28
48.64
47.16
43.88
38.86
33.44
39.9704
29.9799
19.9871
9.9939
4.9975
0.9997
39.9706
29.9787
19.9875
9.993
4.9967
0.9992
39.971
29.9782
19.9867
9.9933
4.9964
0.9988
39.9701
29.9801
19.9871
9.9931
4.9971
0.9978
39.9691
29.9787
19.9861
9.9935
4.9961
0.9991
39.9702
29.9791
19.9869
9.9933
4.9967
0.9989
0.603
0.943
0.425
0.392
0.467
0.116
29.76
20.88
13.10
6.64
3.24
1.08
(2)
Data Analysis
As mentioned in the mathematical model, there is a difference in pressure (P) at each
pressure of PSPC [(Eq. (C.22)] . Thus, this P can be fitted with PSPC
in a linear fitting
program.
In the present case, we have 20 data points as are shown in Table (C.5), the fitted
equation reduces to,
(C.26)
The different fitting parameters with standard uncertainty are shown in Table (C.6). From
the above Eq.(C.26), P is equal to 0.04682 MPa at 60 MPa. It is therefore P is
maximum at 60 MPa but reduces as we decrease the pressure.
Table C.6: Regression Output:
Source Parameters
Fitted Value
Po
U(Po)
(1 )
S1
u(S1 )
(1 )
Degrees of Freedom
-0.001304 (MPa)
0.0012762 (MPa)
0.000802
0.000014
19
The standard uncertainty in P is evaluated from Eq. (C.23) with the sensitivity
coefficient, as mentioned earlier, equals to PSPC .
( )
u2 ( ) = u (o ) + PSPC u (S1 )
2
(C.27)
At the maximum pressure [60 MPa], the standard uncertainty reduces to,
u2 ( ) =
(C.28)
u A = u1 (s ) + u ( )
2
(0.001199)2 + (0.001527)2
= 0.00194 MPa
(C.29)
a =
a+ a
= 0.01MPa
2
(C.30)
uB =
a
0.01
=
= 0.0057MPa
3
3
(C.31)
Degree of freedom i =
Table C.7: Summary of standard uncertainty components
Source of
uncertainty
(X i)
Estimates
( xi )
Limits
x i
(MPa)
Probability
Distribution
Type A or B
- factor
(MPa)
Standard
uncertainty
u(x i )
Sensitivity
coefficient
(MPa)
uB
0.01
0.01
Rectangular
- Type B
3
Normal
-Type A
uA
uc (PDWT )
Expanded
uncertainty
Uncertainty
contribution
u i (y)
Deg. of
freedom
(i )
(MPa)
0.00570
1.0
0.0057
0.00194
1.0
0.00194
19
0.00602
0.012
k = 2.00
uc (PDWT ) = u A2 + u B2
(C.32)
(C.33)
eff =
(uc )4
u14 u24 u B4
+
+
v1 v1
(0.00602)4
(0.001199)4 + (0.001527)4 + (0.0057)4
4
19
(C.34)
Expanded uncertainty
Using the students t-distribution table, k = 2.0 for a confidence level of approximately 95.45%.
Therefore, the expanded uncertainty is given by
k
U=
uc (PDWT)
0.012 MPa
Reporting of results
For the range 0- 60 MPa, the uncertainty U is 0.012 MPa which is approximately 0.02 % of
the full scale pressure. This is determined from a combined standard uncertainty uc = 0.00602
MPa and a coverage factor k = 2.00 based on students distribution for = degrees of
freedom and estimated to have a level of confidence of 95.45 %.
7.
Switch on the measuring equipment and let the auxiliary lamp warm up for 15 minutes.
Mount the standard lamp into the sphere center.
After burning in period, the indirect illumination by standard lamp Es is measured.
Turn the supply voltage down and switch off the standard lamp.
The switched on auxiliary lamp is moved into the sphere. It should remain switched on
always to avoid warm up period. The indirect illuminance EAS is measured.
The standard lamp is moved out of the sphere and the test lamp to be measured is
mounted into the sphere center with the auxiliary lamp still burning, the indirect
illuminance EAT is measured.
Turn on the test lamp to be measured. After burning in period the indirect illuminance
ET is measured.
From the above it is clear that the luminous flux of the test lamp is a function of the luminous
flux of the standard lamp S, the indirect illuminance from the auxiliary lamp EAS and EAT with
standard lamp inside the integrating sphere and test lamp inside the sphere respectively and
the indirect illuminance ES and ET produced by the standard lamp and the test lamp
respectively. The functional dependence of T, can be written as
T, = f (S, E S, E T, E AS, E AT)
(C. 35)
The luminous flux T, of the lamp can be calculated from the luminous flux S of the standard
lamp according to the following relation
T , = s
E s E AS
ET E AT
(C.36)
The factor E AS E AT considers the effect of different sizes and types of test lamps and the
standard lamps.
Since the quantities on the RHS of the Eq. (C.36) are in product form, the equation for the
combined standard uncertainty can be expressed as an estimated relative combined variances
u 2 ( T ) u 2 ( S ) u 2 (E S ) u 2 (ET ) u 2 (E AS ) u 2 (E AT )
=
+
+
+
+
2
2
T2
2S
E S2
ET2
E AS
E AT
(C.37)
If S for standard source is given to be 1045 m and the standard uncertainty is 9.12 m,
from the measurement of expectation value of ES / E T and E AS / E AT, the value T can be
calculated from Eq. (C.36).
Table C.8: Observations
S. No.
1.
2.
3.
4.
5.
6.
7.
8.
9.
10.
E AS
E AS
10.296
10.279
10.254
10.290
10.271
10.272
10.286
10.285
10.254
10.277
10.276
(E AS - E AS)2
400 10-6
9 10-6
484 10-6
196 10-6
25 10-6
16 10-6
100 10-6
81 10-6
484 10-6
441 10-6
For identical standard and test lamp of identical size, shape, same electrical parameters and
same colour temperature, the expectation value of ES / ET and EAS / EAT will be close to 1.
Hence the value of T will be almost equal to S and the uncertainty u (T) can be calculated
from Eq. (C.37). However, we will calculate value of T and standard uncertainty in u (T) by
the following example in which the uncertainties in the measurements of EAS and EAT and ES
and ET are calculated by statistical method and is an example of Type A evaluation of standard
uncertainty.
In the case study , for identical standard lamp and test lamp of identical shape and size , if the
expectation or the average values of the EAS , EAT , ES and ET are 10.276 lux , 10.20 lux, 81.14
lux and 83.76 lux , respectively , the value of T from Eq. (C.36) comes out to be
T = 1045
10.276 83.76
= 1086.6 m
10.20 81.14
(C.38)
Type A Evaluation
We will illustrate by an example the calculation of the standard uncertainty u(E AS ) and relative
standard uncertainty u(E AS) / E AS for one of the parameters, e.g. E AS by Type A evaluation.
(E
10
i =1
Variance
s (E AS
2
AS
E AS
(
)= E
10
AS
i =1
= 2326 10 6
E AS
9
(C.39)
2326 10 6
9
(C.40)
Standard Deviation
s (E AS ) = 16.08 10 3
(C.41)
( )
u E AS
16.08 10 3
=
= 5.08 10 3
10
(C.42)
( )
(C.43)
(C.44)
u E AS
= 5 10 4
E AS
Similarly the relative uncertainties for E AT, E S and E T, can be evaluated by Type A evaluation.
The degree of freedom in each case is also 9 as the total number of observations made in
each case are 10.
The values of the relative uncertainties are
( )
u E AT
E AT
( )
u ES
= 1.2 10 2 and degree of freedom in this case is = i = n - 1 = 10 1 = 9
ES
u (ET )
= 1.6 10 2 and degree of freedom in this case is = i = n - 1 = 10 1 = 9
ET
Type B Evaluation
The uncertainty in the value of S is 9.12 m and is taken from the certificate of the
calibration of the standard lamp. Assuming rectangular distribution, the standard uncertainty in
the value of S is :
u ( s ) =
9.12
= 5.3m
3
(C.45)
u ( s ) 5.3
=
= 5 10 3
1045
s
(C.46)
uc ( T )
for the value of T is calculated using
T
Eq. (C.37).
uc ( T )
= 25 10 6 + 144 10 6 + 256 10 6 + 0.36 10 6 + 0.25 10 6
T2
2
Therefore,
(C.47)
uc ( T )
= 2.06 10-2
T
u c ( T )
T
=
4
4
4
u ( S )
u (E S )
u (ET )
E
E
S
S
+
+ T +
9
9
4
eff
u (E AS )
E
AS +
9
4
u (E AT )
E
AT
9
(C.48)
eff =
[2.06 10 ]
] + [1.6 10 ] + [6 10 ] + [5 10 ]
2 4
[5 10 ] + [1.2 10
3 4
2 4
2 4
4 4
(C.49)
5 4
eff = 18.8 = 19
(C.50)
Estimates
xi
Xi
ES
u(ES )
ET
u(ET )
EAS
u(EAS )
EAT
u(EAT )
u(S )
Combined
uncertainti
es
uc(T )
Expanded
uncertainti
es
Limits
xi
Probability
Distribution
- Type A or B
- Factor
Standard
Uncertainty
u(xi )
Sensitivity Uncertainty
coefficient contribution
ci
ui(y)
Normal
- Type A
- 10
Normal
- Type A
- 10
Normal
- Type A
- 10
Normal
- Type A
- 10
Rectangular
- Type B
1.2 10 2
1.0
1.2 10 2
1.6 10 2
1.0
1.6 10 2
5 10 3
1.0
5 10 3
1.6 10 3
1.0
1.6 10 2
5.3 lm
1.0
5.3 m
22.4 m
19
46.8 m
19
k = 2.09
Deg. of
freedom
i
(C.51)
Where
T= Temperature measured.
D = Displayed temperature of the digital thermometer.
Correction = Correction due to the digital thermometer and Type K thermocouple
Uncertainty evaluation
The combined standard uncertainty (uc ) includes uncertainties of the repeatability of the
displayed readings, the digital thermometer and the thermocouple. This can be represented in
the equation below:
(C.52)
T in 0 C
500.1
500.0
501.1
499.9
499.9
500.0
500.1
500.2
499.9
500.0
Type A evaluation
(A) Standard uncertainty in the readings (u1)
Table (C.10) shows the data obtained from the experiment. Mean value,
T=
1 10
Ti = 500.02
10 i =1
(C.53)
where Ti are the 10 measurements taken as listed in Table (C.10). The temperature of the
chamber after taking into consideration the correction of the thermocouple is 500.50C.
The variance is calculated as follows:
s 2 (Ti ) =
1 n
Ti T
n 1 i =1
1
(0.096) = 0.0106 0 C 2
9
(C.54)
Standard deviation
s (Ti) = 0.103 0 C
(C.55)
()
u1 = s T =
s (Ti ) 0.103 0C
=
= 0.030 C
n
10
(C.56)
Type B evaluation
Standard uncertainty (u2)
From specifications, the uncertainty in the digital thermometer is 0.60 C. Assuming
rectangular distribution, the standard uncertainty in the digital thermometer (u2) is,
u2 =
0.6
= 0.35 0C
3
(C.57)
degree of freedom (i ) =
Standard uncertainty (u3)
From calibration report, with a confidence level of 99 % (k = 2.58), the uncertainty in the
thermocouple is 2.0 0 C.
u3 =
2.0
= 0.78 0 C
2.58
(C.58)
degrees of freedom (i ) =
Combined standard uncertainty
The value of the combined standard uncertainty is calculated using Eq. (C.52)
uc =
= 0.85 0 C
(C.59)
eff =
(0.85)4
(0.03)4
(C.60)
9
Expanded uncertainty
U = k uc = 2 0.85 = 1.7 0 C
(C.61)
Table C.11
Source of
Uncertainty
Estimates Limits
xI
xI
Xi
Digital
Thermometer
Thermocoupl
e
C
0.6
C
0.3
2.0
Repeatability
Combined
uncertainty
Expanded
uncertainties
Probability
Distribution
- Type A or B
- Factor
Rectangular
- Type B
- 3
Normal
- Type B
- 2.58
Normal
- Type A
10
Standard
Uncertainty
u(xi )
0
C
0.35
Sensitivity Uncertainty
coefficient contribution
ci
ui(y)
0
C
1.0
0.35
Deg. of
freedom
i
0.78
1.0
0.78
0.03
1.0
0.03
0.85
1.7
uc
U
k=2
Reporting of results
The temperature of the chamber was measured to be 500.5 0 C. The measurement uncertainty
is 1.7 0 C. The reported measurement uncertainty is estimated at a level of confidence of
approximately 95 % with a coverage factor k of 2.
(C.62)
(C.63)
The product VDC is extremely small and is neglected. The assumption is that the drift in
the values of is small and is also neglected and error of DMM in 1 volt range due to 1 count
is 1V and is also neglected. The precaution is that the interconnecting leads are coaxial
shielded and are kept very small. The reference plane of measurement (mid point of Tee) is
brought close to input plane of DMM (DUC). These precautions minimize the loading as well
as transmission errors. At frequencies up to 10 kHz, with above precautions taken, the error
contribution by above factors are very low ( 2to 3 10-6) and can be neglected.
The inputs are :
1.
The DC calibrator is regularly calibrated at intervals of six months. For the range of 1V
the uncertainty in the calibrator from its calibration certificate is 5.8 10-6 at 95 %
confidence level. Three months stability data from the manufacturers specifications is
5.0 10-6.
2.
The AC /DC transfer correction factor for the thermal converter is + 0.008 %. The
AC/DC transfer uncertainty is 0.01 % at 95 % confidence level.
Readings (V)
0.499986
0.499982
0.499991
0.499994
0.499993
Uncertainty evaluation
VAC = VDC + VDC + VDC
(C.64)
f
u = u 2 ( xi ),
i =1 xi
n
2
c
(C.65)
u2 (V)
u3 (V)
u4 (V)
c1 =
V AC
V AC
V AC
= 1, c 2 =
= 1, and c1 =
=1
V DC
V DC
(C.66)
Type A evaluation
Mean DC Voltage = 0.499989 V,
Standard deviation = 0.0000005 V,
Standard deviation of mean or standard uncertainty
()
sq =
0.000005
= 2.23 10 6V
5
(C.67)
Degrees of freedom
i = 5 1 = 4
(C.68)
Type B evaluation
1.
u1 (V ) =
a1
5.8
=
= 2.96 0.5 V = 1.48 V
1.96 1.96
(C.69)
Degrees of freedom =
2.
u 2 (V ) =
a2
3
5.0
3
(C.70)
Degrees of freedom =
3.
From AC/DC transfer at 95 % confidence level a3 = 100 10-6, distribution is normal and
coverage factor = 1.96
Standard uncertainty
u 3 (V ) =
a3
1 100
=
= 51.02 0.5 V = 25.5 V
1.96
1.96
(C.71)
Degrees of freedom =
(C.72)
(C.73)
eff =
[uc ]4
(u1 )4 + (u2 )4 + (u3 )4 + (u4 )4
eff =
(C.74)
[28.5]4
(1.48)4 + (1.44)4 + (2.23)4 + (25.5)4
(C.75)
Expanded uncertainty
For 95.45 % level of confidence, the coverage factor k = 2, thus
U = kuc (V) = 2 x 28.5 = 57V
(C.76)
Estimates Limits
xi
xi
u1
V
0.5
V
2.9
u2
0.0
2.5
u3
0.0
50.0
Repeatabilit
y
uc(Vac)
Expanded
uncertainty
Probability
Standard
Distribution Uncertainty
- Type A or B u(xi)
- Factor
V
Normal
1.48
-Type B
-1.96
Rectangular
1.44
-Type B
-3
Normal
25.5
-Type B
-1.96
Normal
-Type A
1.0
1.44
1.0
25.5
28.5
28.5
57.0
k=2.0
Reporting of result
The measured average AC Voltage corresponding to 0.500000 V indicated by the DMM,
VAC = 0.499989 (1 + 0.000008 ) V 57 V = 0.499993V 57 V
(C.77)
=
=
=
=
=
(C.78)
Method of measurement
Five separate measurements were taken which involved disconnection and reconnection of
both the unknown sensor and the standard sensor on a power transfer system. All
measurements are made in terms of voltage ratios that are proportional to calibration factor.
None of uncertainty contributions is considered to be correlated.
Mismatch uncertainty:
As the source is not perfectly matched and the phase relation of the reflection coefficients of
the source, the unknown and the standard sensors are not known, there is an uncertainty due
to mismatch for each sensor at the calibration frequency as well as reference frequency. The
corresponding limits of deviation is calculated from the well known formula:
Mismatch uncertainty =
(C.79)
Mismatch uncertainty =
(C.80)
G, s and x are the reflection coefficients for source, the standard and the unknown,
respectively.
(C.81)
C.82)
(C.83)
The long-term stability from the results of five annual calibrations was found to have limits not
greater than 0.4 % per year.
The values of reflection coefficients are themselves uncertain. This uncertainty has been
accounted for by adding it in quadrature with the actual measured values.
The instrumentation linearity uncertainty has been estimated to lie within 0.1 % from
measurements against a reference attenuation standard up to ratios of 2:1 at confidence level
of 95 %.
Reference source:
The ratio of power outputs of the reference source has been estimated to be 1.000 with
deviations 0.004.
Standard sensor:
The standard sensor was calibrated 6 months ago. The value of calibration factor from its
calibration certificate is
0.965 0.012
at confidence level of 95 %.
Uncertainty evaluation
Type A evaluation
The measured values of calibration factor for the unknown power sensor are shown in Table
(C.14).
The mean value is
K x = 0.9496 0.950,
(C.84)
Calibration Factor
0.958
0.946
0.951
0.950
0.943
s(K x ) = 0.0057
(C.85)
Standard uncertainty in u (s )K x is
u( K x ) =
0.0057
= 0.0025
5
(C.86)
Degree of freedom = 5 - 1 = 4.
Type B evaluation
1.
2.
(C.87)
Uncertainty in drift in calibration factor since its last calibration is 0.002. This is a
rectangular distribution and the standard uncertainty,
u ( Ds ) =
0.002
= 0.0012 , and
3
(C.88)
Degrees of freedom =
3.
u( REF ) =
0.004
= 0.0023, and
3
(C.89)
Degrees of freedom =
4.
Uncertainty due to the instrument linearity is 0.001. This is a normal distribution and
the standard Uncertainty u (DC)
u( DC ) =
0.001
= 0.0023, and
1.96
(C.90)
Degrees of freedom =
5.
This is U-shaped and the corresponding associated standard uncertainty figures are:
u (M s ) =
0.0008
= 0.00056 at 50MHz
2
(C.91)
u (M x ) =
0.0008
= 0.00056 at 50MHz
2
(C.92)
u (M s ) =
0.0108
= 0.0076 at 18GHz
2
(C.93)
u (M x ) =
0.012
2
= 0.0085 at 18GHz
(C.94)
(C.95)
u c (K x ) 0.0134
(C.96)
Estimates Limits
Probability
Distribution
- Type A or B
Standard
Sensitivity
Uncertainty Coefficient
u(xi)
ci
Uncertainty Degree of
Contribution Freedom
ui (y)
vi
Normal
-Type B
Rectangular
-Type B
Normal
-Type B
Rectangular
-Type B
0.0061
1.0
0.0061
0.0012
1.0
0.0012
0.0005
1.0
0.0005
0.0023
1.0
0.0023
xi
xi
Ks
0.965
0.012
Ds
0.002
0.002
DC
1.0
0.001
REF
1.0
0.004
0.0
0.0
0.0008
0.0008
-U shaped
-do-
0.00056
0.00056
1.0
1.0
0.00056
0.00056
0.0
0.0
0.0108
0.012
-do-do-Type B
Normal
-Type A
5
0.0076
0.0085
1.0
1.0
0.0076
0.0085
0.0025
1.0
0.0025
0.0134
0.0263
3301
Mismatch at
50 MHz
S
X
Mismatch at
18 GHz
S
X
Repeatabilit
y
uc(Kx)
Expanded
uncertainty
k = 1.96
(C.97)
Reported Result:
The calibration factor of unknown power sensor at 18 GHz is 0.450 0.027
The reported expanded uncertainty of measurement is stated as the combined standard
uncertainty multiplied by the coverage factor 1.96, which for a normal distribution corresponds
to a coverage probability or confidence level of 95 %.
(C.98)
errors due to loading and wire leads or connections are considered negligible, and
all input quantities are uncorrelated.
The calibrator is calibrated regularly at the interval of six months. For the range of 10 V,
the specifications are: Resolution = 10 V with the uncertainty at 99 % level of
confidence as (4.5 10-6 of output + 100 V),
Vx = Vs + Vx
(B)
(C.99)
4 digital multimeter specification are: for range 100 V full display is 99.99 with
resolution = 10 mV and the uncertainty at 99 % level of confidence as (10-5 of reading
+ 0.2 10-5 of full scale).
Observation:
Applied voltage from calibrator
10.00000 V
Even after repetition of the readings the multimeter reading indicating the same value or 1
due to digitizing process. This is due to better accuracy of the reference standard (calibrator).
f
u =
i =1 xi
N
2
c
2
u ( xi )
(C.100)
c1 =
V x
V x
= 1 and c2 =
=1
Vs
V x
(C.101)
Type A evaluation:
Even after repetition of the readings the multimeter reading indicating the same value or 1 due
to digitizing process. This is due to better accuracy of the reference standard (calibrator). In this
case, type A uncertainty can be assumed as negligible and the repeatability uncertainty can be
treated as type B uncertainty using the resolution error of the multimeter.
Type B evaluation:
(I) The uncertainty in applied voltage from the calibrator is
a1 =
(4.5 10 + 100) V
145 V
(C. 102)
At 99 % confidence level assuming normal distribution, coverage factor k = 2.58, the standard
uncertainty in applied voltage is
u1 (V) = a1/2.58 = (145 / 2.58) V = 56.20 V
(C.103)
Degree of freedom is = v1 =
(C.104)
(II)
a2 =
10
mV = 5mV
2
(C.105)
For rectangular distribution, the standard uncertainty due to the resolution uncertainty of
the multimeter is:
u 2V =
a2
3
5
3
mV = 2886.75 V
(C.106)
Degree of freedom is = v2 =
(C.107)
as v1 = and
(C.108)
v2 =
Estimates
xi
Vs
10.00 V
Vx
0.01 V
Repeatabilit
y
Limits
xi
5 mV
Probability
Distribution
- type A or B
Standard
Uncertainty
u(xi)
Sensitivity
coefficient
ci
Uncertainty
contribution
ui (y)
Degree of
Freedom
vi
Normal
Type B
- 2.5
56.2 V
1.0
56.2 V
Rectangular
- Type B
-3
2886.25 V
1.0
2886.75 V
Normal
- Type A
uc(Vx)
Expanded
Uncertainty
k=2
2.89 mV
5.78 mV
CORE GROUP