Учебно-методический комплекс по учебной дисциплине Контроль и испытания продукции для направления специальности 1-54 01 01-01 Метрология, стандартизация и сертификация (машиностроение и приборостроение)
.pdf461
● данные, полученные в результате опыта, или общие знания
оповедении и свойствах соответствующих материалов и приборов;
●спецификация изготовителя;
● данные, которые приводятся в свидетельствах о калибровке
исертификатах;
●справочные данные.
Корректное использование фонда доступной информации для
оценивания стандартной неопределенности по типуВ требует интуиции,
основанной на опыте и общих знаниях, и является мастерством, которое приходит с практикой. Имеющуюся информацию и знания или даже предположения о величинах Xi необходимо описать с помощью распределения вероятностей, чтобы затем определить оценки величин и их стандартные отклонения. При этом в метрологической практике чаще всего используются распределение Гаусса (нормальное), прямоугольное
(равномерное), треугольное, трапецеидальное.
Если бы лаборатория обладала неограниченными ресурсами, то все составляющие неопределенности оценивались бы статистическими методами на основе экспериментальных данных. Однако в некоторых измерительных ситуациях, когда не представляется возможным организовать достаточное количество независимых наблюдений, оценка стандартной неопределенности по типуВ может быть такой же надежной, как и оценка по типу А.
Две входные величины могут быть независимы или связаны между собой, то есть, взаимозависимы или коррелированы. С позиций концепции неопределенности измеряемая величина трактуется как скаляр, в то же время ряд связанных измеряемых величин, определенных одновременно в том же самом измерении, требует замены скалярной измеряемой величины и ее дисперсии на векторную измеряемую величину и ковариационную матрицу.
Значительная ковариация между двумя входными величинами может наблюдаться в случае, если «при их определении используют один и тот же измерительный прибор, физический эталон измерения или справочные
462
данные, имеющие значительную стандартную неопределенность» [2, с. 23].
Например, если поправка на температуру, необходимая для оценки хi одной входной величины Хi,получается с помощью некоторого термометра и такая же поправка на температуру, необходимая для оценки xj входной величины
Xj, тоже получается с помощью этого же термометра, то две входные величины Хi и Xj могут быть значительно коррелированы.
Ковариация двух случайных переменных xи z является мерой их взаимной зависимости и определяется по формуле:
cov(x,z) cov(z, x) (x x )(z z ) p(x, z)dxdz |
(5.10) |
|
Ковариация cov(x,z) может быть оценена с помощью дисперсии s(xi,zi),
полученной из n независимых пар xiи zi одновременных наблюдений xи z:
|
|
|
1 |
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
||||
s(xi , zi |
) |
(xi |
x |
)(zi |
|
z |
) |
|
|
||||||||||||
|
|
|
|
||||||||||||||||||
|
|
|
n 1 i 1 |
(5.11) |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||
Оцененная ковариация средних значений x и z определяется как: |
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
s(xi , zi ) |
|
|
|
|
|
|
|
|
||||
|
s( х ,z ) |
|
|
|
|
|
|
|
|||||||||||||
|
n |
(5.12) |
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||
Степень корреляции |
между xи z |
характеризуется оцененным |
|||||||||||||||||||
коэффициентом корреляции: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
r(x, z) |
|
u(x, z) |
|
|
|
|
|
|
|
|||||||||||
|
u(x )u(z) |
, |
(5.13) |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где r(x; z) =r(z; x) и –1 ≤r(x; z) ≤ +1.
Ковариация, связанная с оценками двух входных величин Хи Z может устанавливаться равной нулю или рассматриваться как пренебрежимо малая,
если:
а) величины Х и Zявляются независимыми друг от друга, например,
если они в различных, независимых один от другого экспериментах многократно, но не одновременно наблюдались или если они представляют
(описывают) результирующую величину различных, независимых друг от друга проведенных исследований;
б) одна из входных величин Х или Z может рассматриваться как константа;
463
в) исходя из наших знаний и предположений просто не имеется
никаких оснований для существования корреляции между входными величинами Хи Z.
Расчет значения измеряемой или выходной величины Yзаключается
в нахождении значения ее оценки у.Полученное значение оценки
принимается за результат измерения. Оценку выходной величины получают
из уравнения модели, заменяя входные величины Xi их оценками xi:y =f(x1,x2,...,xn).
При этом предполагается, что значения входных величин являются лучшими оценками входных величин, что были внесены поправки на влияния и эффекты, значимые для данной модели. Если это не так, то необходимые поправки должны вводиться в модель в качестве отдельных
входных величин.
Приведенные уравнения базируются на аппроксимации функции модели y =f(x1,x2,...,xN)рядом Тейлора первого порядка. Это справедливо для
линейных функций. При значительной |
нелинейности функции f(х)в ряд |
|||||
Тейлора нужно включать члены более высокого порядка. |
|
|||||
С учетом |
ковариации |
выражение |
для |
суммарной |
||
неопределенностизапишется в виде: |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
n |
n 1 |
n |
|
|
|
uc ( y) ci2 u 2 (xi ) 2 ci c j u(xi , x j ) |
|
|
|||
|
|
i 1 |
i 1 j i 1 |
, |
(5.14) |
|
|
|
|
|
|
||
|
|
|
|
|
||
где u(xi,xj) – оцененная ковариация |
xiи xj; сi, |
сi |
коэффициенты |
|||
чувствительности, |
показывающие, |
как |
выходная |
оценка |
y изменяется |
|
с изменением значений входных оценок |
x1,x2,...xn, в большинстве случаев |
|||||
коэффициенты чувствительности рассчитываются как частные производные:
c f
i |
xi |
|
|
|
(5.15) |
||
|
|
|
|
Таким образом, выражение (5.14) запишется в виде: |
|
||
|
n |
|
f |
|
2 |
|
n 1 n |
|
|
|
uc ( y) |
|
|
|
u |
2 |
(xi ) 2 |
|
|||
|
|
|
|
|
|
|||||
|
|
|
|
|
|
x j |
||||
|
i 1 |
|
xi |
|
|
i 1 j i 1xi |
|
|||
u(xi , x j )
, |
(5.16) |
Если выразить взаимозависимость величин через корреляции, выражение (5.14) примет вид:
|
n |
|
f |
|
2 |
|
n 1 n |
|
|
|
|
uc ( y) |
|
|
|
u |
2 |
(xi ) 2 |
|
u(xi )u(x j )r(xi , x j ) |
|||
|
|
|
|
|
|
||||||
|
|
|
|
|
|
x j |
|||||
|
i 1 |
|
xi |
|
|
i 1 j i 1xi |
|
|
|||
Для некоррелированных величин суммарная неопределенность рассчитывается из выражений:
464
коэффициент
(5.17)
стандартная
|
|
n |
|
|
|
|
|
|
|
|
uc ( y) |
ci2 u 2 (xi ) |
|
|
(5.18) |
||||
|
|
i 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
f |
|
2 |
|
|
|
|
|
uc ( y) |
|
|
2 |
|
|
|
|
|
|
|
u |
(xi |
) |
|
|
|||
|
|
xi |
|
|
|
||||
|
|
i 1 |
|
|
|
|
|
(5.19) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Коэффициенты |
чувствительности |
|
иногда |
определяются |
|||||
экспериментальным |
путем с помощью |
измерения изменения выходной |
|||||||
величины Y, вызванного изменением в входной величины Xi, поддерживая при этом остальные входные величины неизменными. В этом случае знание
функции у= f(х) сводится к эмпирическому разложению в ряд Тейлора первого порядка, основанного на измеренных коэффициентах чувствительности.
В случае коррелированных входных величин знак u(x,z) обязательно должен приниматься во внимание.
Если функция модели f является суммой или разностью
некоррелированных |
входных |
величин |
|
Xi |
с коэффициентами |
чувствительности, равными множителям рi: |
|
|
|
||
|
|
|
|
|
|
|
|
N |
|
|
|
|
y f ( X1 , X 2 ,..., X N ) pi |
X i |
, |
(5.20) |
|
|
|
i 1 |
|
||
|
|
|
|
|
|
|
|
|
|
|
|
то выражение для суммарной неопределенности выходной величины запишется в виде:
(5.21)
Если функция модели f является произведением или отношением некоррелированных входных величин Xi со степенями рi
|
|
|
465 |
|
|
|
|
|
N |
|
|
y f ( X1 , X 2 |
,..., X N ) X ipi |
, |
(5.22) |
|
i 1 |
||
|
|
|
|
то можно записать выражение для расчета относительной суммарной
неопределенности (в этом случае |
|
нет |
|
|
необходимости в вычислении |
||||
коэффициентов чувствительности): |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
u c ( y) |
|
N |
|
u(xi |
) |
2 |
|
|
|
|
|
pi |
|
|
||||
|
|
|
|
|
|||||
|
y |
i 1 |
|
xi |
|
|
(5.23) |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Бюджет неопределенности (табл. 5.1) служит для обобщения и анализа |
|||||||||
вкладов каждого источника неопределенности в суммарную стандартную неопределенность с целью определения точности измерительного процесса,
корректировки модели измерения или поиска способов уменьшения значений неопределенностей от некоторых вкладов.
Вклад в суммарную неопределенностьui(y)(i=1, 2,...,n) рассчитывается
по формуле: |
|
ui ( y) ci u(xi ) |
(5.24) |
Для прямого измерения сi = 1, поэтому вклады входных величин будут численно равны их стандартным неопределенностям:ui(y) =u(хi).
Следует обратить внимание, что в то время как стандартная неопределенность u(xi) всегда положительна, вклад в неопределенность ui(y)
в зависимости от знака коэффициента чувствительности может принимать положительное или отрицательное значение. В случае некоррелированных величин этот знак не играет роли, так как в этом случае при расчете суммарной стандартной неопределенности вклады в неопределенность u(xi)
возводятся в квадрат.
Расширенную неопределенность U получают путем умножения суммарной стандартной неопределенности uс(y) на коэффициент охватаk:
U p(y) k uc ( y) |
(5.25) |
|
|
В случае указания расширенной неопределенности результат измерения выражается в виде интервала Y =y ±U (y –U y y +U), который содержит большую часть распределения вероятностей, характеризуемого
|
|
466 |
результатом измерения и его суммарной |
стандартной неопределенностью, |
|
и является вероятностью охвата или уровнем доверия этого интервала. |
||
При |
выборе значения коэффициента охвата k необходимо полное |
|
знание о |
распределении вероятностей |
выходной величины. При этом |
вконцепции неопределенности используются положения Центральной Предельной теоремы, в соответствии с которой можно допустить, что распределение вероятностей, характеризуемое результатом измерения и его суммарной стандартной неопределенностью, может считаться нормальным
вслучаях, когда:
а) оценку y измеряемой величины Y получают из оценок хi числа n3
3 входных величин, которые описываются «хорошо ведущими себя» распределениями, такими как нормальное, прямоугольное, треугольное;
б) |
стандартные |
неопределенности u(xi) этих |
оценок дают |
|||
сопоставимые вклады |
в суммарную неопределенность |
uc(y), связанную |
||||
с оценкой выходной величины у; |
|
|
|
|
||
в) |
оценки неопределенностей |
входных |
величин |
являются |
||
достаточно надежными. |
|
|
|
|
|
|
Следовательно, |
в случаях, когда |
выходной |
величине |
может |
||
приписываться нормальное распределение вероятностей, для расчета расширенной неопределенности Up =kpuc(y), которая обеспечивает уровень доверия р, можно использовать для kp значения из нормального закона распределения. На практике, а также при проведении калибровокпринимают k = 2 для интервала, имеющего уровень доверия р = 95 % и k = 3 для интервала, имеющего уровень доверия р = 99 %.
При оценивании стандартной неопределенности по типуА из малого числа (n < 10) повторных наблюдений неопределенность такой оценки
(«неопределенность неопределенности») может достигать 50 % при n = 3.
Когда выполняется условия центральной предельной теоремы, но не выполняется условие надежности распределение вероятностей результата измерения описывается распределением Стьюдента (t-распределением)
467
с эффективными степенями свободы νeff.
свободы eff для стандартной uc(y)осуществляется с помощью уравнения
Оценка эффективных степеней неопределенности измерения
Велча-Саттерсвейта:
|
eff |
|
uс |
4 |
( y) |
|
|
|||
|
|
|
|
4 |
|
|
|
|
||
|
|
|
ui ( y) |
|
|
|||||
|
|
|
N |
|
|
|
|
|
|
|
|
|
|
i 1 |
|
i |
|
|
, |
(5.26) |
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
где i эффективная степень свободы вклада в неопределенность ui(y). |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
N |
|
|
|
|
|
|
|
|
νeff ≤ |
i |
. |
|
|
(5.27) |
||||
|
i 1 |
|
|
|
|
|||||
|
|
|
|
|
|
|||||
|
|
|
|
|
||||||
|
Критерий надежности |
в общем полностью выполняется, если вклад |
||||||||
в неопределенность, обусловленный стандартной неопределенностью
входной величины, оцененной по типу А, определялся из числа повторных наблюдений не менее 10. Для значения стандартной неопределенности измерения u(xi), которое определяется по методу оценивания неопределенности по типуА, степень свободы рассчитывается как i =n –1. Cтепени свободы для стандартной неопределенности измерения, значение
которой определено по методу оценивания неопределенности по типу В,
могут приниматься i .Коэффициент охвата k определяют из табл. 5.2.
Эта таблица базируется наt-распределении, которое установлено для вероятности охвата р = 95,45 % и р = 99,73 %. Если effявляется не целым числом, то его уменьшают до ближайшего целого числа.
Таблица 5.1 – Бюджет неопределенности
ВеличинаX |
Единицаизмерения |
хОценка |
Интервалот–адо+а |
Тип оцениваниянеопредел енности |
Распределение вероятностей |
Стандартнаянеопределенностьu(x |
Количествостепеней свободыν |
Коэффициентчувстви тельностис |
вВкладсуммарную неопределенность u |
Процентныйвклад, % |
|
|
|
|
|
|
) |
|
|
|
|
|
|
|
|
|
|
i |
|
|
|
|
i |
|
|
|
|
|
|
|
i |
|
|
|
|
i |
|
|
|
|
|
|
|
|
|
|
|
|
(y) |
|
|
|
|
i |
|
|
|
|
|
X1 |
x1 |
u(x1) |
с1 |
u1(y) |
Х2 |
х2 |
u(x2) |
с2 |
u2(y) |
.... |
|
|
|
|
ХN |
хn |
u(xn) |
сn |
un(y) |
Y |
y |
|
|
uc(y) |
468
Таблица 5.2 – Значения эффективных степеней свободы и коэффициентов
Стьюдента
eff |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
10 |
20 |
50 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k95 |
13,97 |
4,53 |
3,31 |
2,87 |
2,65 |
2,52 |
2,43 |
2,37 |
2,28 |
2,13 |
2,05 |
2,00 |
|
|
|
|
|
|
|
|
|
|
|
|
|
k99 |
235,8 |
19,21 |
9,22 |
6,22 |
5,51 |
4,90 |
4,53 |
4,28 |
3,96 |
3,42 |
3,16 |
3,00 |
|
|
|
|
|
|
|
|
|
|
|
|
|
Когда вклад неопределенности какой-либо входной величины
в суммарную стандартную неопределенность является доминирующим, то распределение вероятностей выходной величины предполагается аналогичным распределению вероятностей доминирующего вклада.
Для сравнения вкладов источников неопределенности в суммарную стандартную неопределенность используют понятие процентного вклада.
Процентный вклад неопределенности i-ой входной величины в суммарную
стандартную неопределенность вычисляется как:
u2 |
y |
(5.28) |
|
u2 |
y 100% |
||
i |
|
|
|
c
Сумма всех процентных вкладов должна быть равна 100 %.
При записи результата измерения рекомендуется применять следующие способы. В качестве примера рассмотрим эти способы записи для эталона массы msc номинальным значением 100 г. Если мерой неопределенности измерения является суммарная стандартная неопределенностьuс(y):
1)«ms = 100,02147 г (с суммарной стандартной неопределенностью)
uс(y) = 0,35 мг»
или «100,02147 г; 0,35 мг»;
или «100,02147 г; 3,5·10-6», где число после «..» без указания единиц величин является относительной стандартной неопределенностью uc(y)/│y│;
2) «ms = 100,02147(35) г», где цифры в скобках являются численным значением суммарной стандартной неопределенности uc(y),
соответствующие последним цифрам приведенного результата;
469
3) «ms = 100,02147(0,00035) г», где число в скобках является
численным значением суммарной стандартной неопределенности uc(y),
выраженной в единицах величин результата измерения;
или |
«100,02147 г |
(0,35 мг)», |
где |
число |
в скобках |
является |
численным |
|
значением суммарной стандартной неопределенности |
uс(y), |
выраженной |
||||||
в указанных единицах величин; |
|
|
|
|
|
|
||
или |
«100,02147 г |
(3,5·10-6)», |
где |
число |
в скобках |
без |
указания единиц |
|
величин являются относительной стандартной неопределенностью uс(y)/│y│;
Если мерой неопределенности измерения является расширенная неопределенность U, то лучше всего указать результат в виде y U.При
этомнаиболее полной будет следующая форма записи:
«ms = (100,02147 ± 0,00079) г», где число следующее за знаком ±, является численным значением расширенной неопределенности U =kuс(y), причем U
определено из |
суммарной |
стандартной неопределенности |
uс(y) = 0,35 мг |
и коэффициента |
охватаk = |
2, основанного на нормальном |
распределении, |
и определяет интервал, оцененный как имеющий уровень доверия р = 95 %»;
а более краткой «100,02147 г ± 0,00079 г (k =2, р = 95 %)»
или «100,02147 г ± 0,79 мг (k =2, р = 95 %)», где в скобках указано значение коэффициента охвата и уровня доверия.
Указание неопределенности измерения в свидетельстве о калибровки может быть представлено в следующей форме. В свидетельстве о калибровке
должен быть указан полный результат измерения, состоящий из оценки
у измеряемой величины и связанной с нею расширенной неопределенности измерения U в форме y U. При этом рекомендуется использовать фразу следующего содержания:
«Указанная расширенная неопределенность является произведением
стандартной неопределенности измерения и коэффициента охвата |
k = 2, |
и соответствует при нормальном распределении вероятности |
охвата |
приблизительнор = 95 %». |
|
470
В случаях, где имеет смысл нахождение степеней свободы, эта фраза
должна звучать следующим образом:
«Указанная расширенная неопределенность является произведением стандартной неопределенности измерения и коэффициента охвата k =XX,
и соответствует при t-распределении с eff =YY эффективными степенями свободы вероятности охвата приблизительно р = 95 %».
Значения оценки у и ее стандартной неопределенности uc(y) или расширенной неопределенности U не следует давать с избыточным числом
цифр. Обычно достаточно привести их с двумя значащими цифрами, хотя в некоторых случаях может быть необходимо сохранить дополнительные цифры для того, чтобы избежать погрешности округления в следующих расчетах.
При сообщении окончательных результатов иногда может быть уместным округлить неопределенности в сторону увеличения, а не до ближайшей цифры.
Пример:uc(y) = 10,47 мОм можно округлить до 11 мОм. Однако здравый смыслдолжен возобладать, и значение, такое как uc(y) = 28,05 кГц,
следует округлить до 28 кГц.
Выходные и входные оценки должныокругляться так, чтобы соответствовать своим неопределенностям.
Пример: еслиу = 10,05762 Ом с uc(y) = 27 мОм, то у следует округлить до 10,058. Коэффициенты корреляции должны даваться с точностью до третьей цифры, если их абсолютные значения близки к единице.
