Руководстве по выражению неопределенности измерения gum

Утвержден и введен в действие

Приказом Федерального

агентства по техническому

регулированию и метрологии

от 22 ноября 2013 г. N 1665-ст

НАЦИОНАЛЬНЫЙ СТАНДАРТ РОССИЙСКОЙ ФЕДЕРАЦИИ

НЕОПРЕДЕЛЕННОСТЬ ИЗМЕРЕНИЯ

ЧАСТЬ 3

РУКОВОДСТВО

ПО ВЫРАЖЕНИЮ НЕОПРЕДЕЛЕННОСТИ ИЗМЕРЕНИЯ

ДОПОЛНЕНИЕ 2

ОБОБЩЕНИЕ НА СЛУЧАЙ ПРОИЗВОЛЬНОГО ЧИСЛА ВЫХОДНЫХ ВЕЛИЧИН

Uncertainty of measurement. Part 3. Guide to the expression

of uncertainty in measurement. Supplement 2.

Extension to any number of output quantities

ISO/IEC Guide 98-3:2008/Supplement 2:2011

Uncertainty of measurement — Part 3: Guide to the expression

of uncertainty in measurement (GUM:1995) — Supplement 2:

Extension to any number of output quantities

(IDT)

ГОСТ Р 54500.3.2-2013/

Руководство ИСО/МЭК 98-3:2008/

Дополнение 2:2011

Группа Т80

ОКС 17.020

Дата введения

1 сентября 2014 года

Предисловие

1. Подготовлен Федеральным государственным унитарным предприятием «Всероссийский научно-исследовательский институт метрологии им. Д.И. Менделеева» (ФГУП ВНИИМ) и Автономной некоммерческой организацией «Научно-исследовательский центр контроля и диагностики технических систем» (АНО «НИЦ КД») на основе собственного аутентичного перевода на русский язык международного документа, указанного в пункте 4.

2. Внесен Техническим комитетом по стандартизации ТК 125 «Применение статистических методов».

3. Утвержден и введен в действие Приказом Федерального агентства по техническому регулированию и метрологии от 22 ноября 2013 г. N 1665-ст.

4. Настоящий стандарт идентичен международному документу Руководство ИСО/МЭК 98-3:2008/Дополнение 2:2011 «Неопределенность измерения. Часть 3. Руководство по выражению неопределенности измерения. Дополнение 2. Обобщение на случай произвольного числа выходных величин» [ISO/IEC Guide 98-3:2008/Supplement 1:2011 «Uncertainty of measurement — Part 3: Guide to the expression of uncertainty in measurement (GUM:1995) — Supplement 2: Extension to any number of output quantities»].

При применении настоящего стандарта рекомендуется использовать вместо ссылочных международных документов соответствующие им национальные стандарты, сведения о которых приведены в дополнительном Приложении ДА.

5. Введен впервые.

Правила применения настоящего стандарта установлены в ГОСТ Р 1.0-2012 (раздел 8). Информация об изменениях к настоящему стандарту публикуется в ежегодном (по состоянию на 1 января текущего года) информационном указателе «Национальные стандарты», а официальный текст изменений и поправок — в ежемесячном информационном указателе «Национальные стандарты». В случае пересмотра (замены) или отмены настоящего стандарта соответствующее уведомление будет опубликовано в ближайшем выпуске ежемесячного информационного указателя «Национальные стандарты». Соответствующая информация, уведомление и тексты размещаются также в информационной системе общего пользования — на официальном сайте Федерального агентства по техническому регулированию и метрологии в сети Интернет (gost.ru).

Введение

В «Руководстве по выражению неопределенности измерений» (GUM) [JCGM 100:2008] рассматриваются в основном одномерные модели измерений, включающие в себя единственную скалярную выходную величину. Однако на практике часто встречаются измерительные задачи с двумя и более выходными величинами. Примеры таких задач имеются в GUM для случаев электрических измерений с тремя выходными величинами [JCGM 100:2008 (раздел H.2 приложения H)] и температурных измерений с двумя выходными величинами [JCGM 100:2008 (раздел H.3 приложения H)]. В настоящем стандарте рассматриваются многомерные модели измерения, включающие в себя произвольное число выходных величин. В большинстве случаев выходные величины коррелированны, поскольку зависят от общих входных величин. В настоящем стандарте рассматривается обобщение способа оценивания неопределенности по GUM [JCGM 100:2008 (раздел 5)], позволяющее получить оценки выходных величин, а также стандартные неопределенности и ковариации, соответствующие этим оценкам. Входные и выходные величины модели измерения могут быть действительными или комплексными.

Дополнение 1 к GUM [JCGM 101:2008] рассматривает трансформирование распределений [JCGM 100:2008, 5] при заданной модели измерения как основу для выражения неопределенности измерения и реализацию данной процедуры посредством метода Монте-Карло [JCGM 100:2008 (раздел 7)]. Как и в GUM, в нем рассмотрены только модели с единственной скалярной выходной величиной [JCGM 101:2008 (раздел 1)]. Настоящий стандарт рассматривает обобщение метода Монте-Карло с целью получения дискретного представления совместного распределения вероятностей для выходных величин многомерной модели. Такое дискретное представление служит основой для получения оценок выходных величин, их стандартных неопределенностей и ковариации. Использование метода Монте-Карло является альтернативой способу оценивания неопределенности по GUM, особенно в ситуациях, когда последний не способен обеспечить достоверные результаты измерений вследствие того, что (а) линеаризация модели приводит к существенному искажению результатов измерения или (б) распределение вероятностей для выходной величины (или величин) не может быть описано многомерным нормальным распределением.

Настоящий стандарт устанавливает также метод определения области охвата для выходных величин многомерной модели, являющейся аналогом интервала охвата в случае одномерной модели, для заданной вероятности охвата. Рассматриваются области охвата в форме эллипсоидов или прямоугольных параллелепипедов. Применение численных процедур расчета неопределенности измерения с использованием метода Монте-Карло дает возможность приближенного построения областей охвата наименьшего объема.

1. Область применения

Настоящий стандарт является дополнением к «Руководству по выражению неопределенности измерений» (GUM) (JCGM 100) и распространяется на модели измерения с произвольным числом входных и выходных величин. Входящие в модель измерения величины могут быть действительными и/или комплексными. Рассмотрено два подхода к использованию таких моделей. Первый представляет собой обобщение способа оценивания неопределенности по GUM. Второй — использование метода Монте-Карло для трансформирования распределений. Использование метода Монте-Карло дает возможность получить достоверные результаты в ситуациях, когда условия применимости первого подхода не выполняются.

Способ оценивания неопределенности по GUM применим, когда информацию о входных величинах можно представить в виде их оценок (например, полученных измерением), связанных с этими оценками стандартных неопределенностей и, при необходимости, ковариаций. Использование соответствующих формул и процедур позволяет на основе указанной информации получить оценки, а также соответствующие им стандартные неопределенности и ковариации для выходных величин. Эти формулы и процедуры применимы к моделям измерения, для которых выходные величины (a) выражены непосредственно как функции от выходных величин (функции измерения) или (b) могут быть получены решением уравнений, связывающих входные и выходные величины.

В целях упрощения формулы, применяемые в настоящем стандарте, даны в матричной форме записи. Дополнительным преимуществом такой формы записи является ее приспособленность к реализации на многих языках программирования и в системах, которые поддерживают матричную алгебру.

Способ оценивания неопределенности измерения с применением метода Монте-Карло основывается на (i) присвоении входным величинам модели измерения соответствующих распределений вероятностей [JCGM 101 (раздел 6)], (ii) определении дискретного представления совместного распределения вероятности для выходных величин и (iii) получения из этого дискретного представления оценок выходных величин, их стандартных неопределенностей и ковариаций. Данный подход является обобщением метода Монте-Карло, установленного в JCGM 101 применительно к моделям с единственной скалярной выходной величиной.

Применение вышеуказанных подходов позволяет получить при заданной вероятности охвата область охвата для выходных величин многомерной модели — аналог интервала охвата для одномерной модели с единственной скалярной выходной величиной. Рассматриваемые в настоящем стандарте области охвата имеют формы гиперэллипсоидов (далее — эллипсоидов) и прямоугольных гиперпараллелепипедов (далее — параллелепипедов) в многомерном пространстве выходных величин. В случае применения метода Монте-Карло приведена также процедура приближенного построения области охвата наименьшего объема.

Применение стандарта иллюстрировано подробными примерами.

Настоящий стандарт служит дополнением к GUM и должен быть использован вместе с ним и с Дополнением 1 к GUM (соответственно, JCGM 100 и JCGM 101). Настоящий стандарт предназначен для тех же пользователей, что и два вышеуказанных документа (см. также JCGM 104).

2. Нормативные ссылки

В настоящем стандарте использованы нормативные ссылки на следующие документы:

JCGM 100:2008 Оценивание данных измерений. Руководство по выражению неопределенности измерения (GUM) [JCGM 100:2008, Evaluation of measurement data — Guide to the expression of uncertainty in measurement (GUM)]

JCGM 101:2008 Оценивание данных измерений. Дополнение 1 к «Руководству по выражению неопределенности измерения». Трансформирование распределений с использованием метода Монте-Карло (JCGM 101:2008, Evaluation of measurement data — Supplement 1 to the «Guide to the expression of uncertainty in measurement» — Propagation of distributions using a Monte Carlo method)

JCGM 104:2009 Оценивание данных измерений. Введение к «Руководству по выражению неопределенности измерения» и сопутствующим документам (JCGM 104:2009, Evaluation of measurement data — An introduction to the «Guide to the expression of uncertainty in measurement» and related documents)

JCGM 200:2008 Международный словарь по метрологии. Основные и общие понятия и связанные с ними термины (VIM) [JCGM 200:2008, International Vocabulary of Metrology — Basic and general concepts and associated terms (VIM)].

3. Термины и определения

В настоящем стандарте применены термины по JCGM 100 и JCGM 200, некоторые из которых (при необходимости, модифицированные) приведены в настоящем разделе, а также следующие термины с соответствующими определениями (обозначения, использованные в настоящем стандарте, приведены в Приложении D).

3.1. Действительная величина (real quantity): величина, числовое значение которой является действительным числом.

3.2. Комплексная величина (complex quantity): величина, числовое значение которой является комплексным числом.

Примечание. Комплексная величина может быть представлена двумя действительными величинами в форме алгебраической

или тригонометрической

где символ «^T» обозначает транспонирование;

i — мнимая единица, i² = -1;

Z_R и Z_I — соответственно действительная и мнимая части ;

Z_r и — соответственно модуль и аргумент .

3.3. Векторная величина (vector quantity): совокупность величин, упорядоченных в виде элементов матрицы с одним столбцом.

3.4. Действительная векторная величина (real vector quantity): векторная величина, элементами которой являются действительные величины.

Пример — Действительная векторная величина , состоящая из N элементов (действительных чисел) X₁, …, X_N может быть представлена в виде матрицы размерности N x 1 (матрицы-столбца):

3.5. Комплексная векторная величина (complex vector quantity): векторная величина, элементами которой являются комплексные величины.

Пример — Комплексная векторная величина , состоящая из N элементов (комплексных чисел) Z₁, …, Z_N может быть представлена в виде матрицы размерности N x 1 (матрицы-столбца):

3.6. Векторная измеряемая величина (vector measurand): векторная величина, подлежащая измерению.

Примечание. Данное определение модифицировано по отношению к JCGM 200 (словарная статья 2.3).

3.7. Модель (измерения) (measurement model): математическое соотношение между всеми величинами, используемыми для получения результата измерения.

Примечание 1. Данное определение модифицировано по отношению к JCGM 200 (словарная статья 2.48).

Примечание 2. В общем виде модель измерения имеет вид уравнения h(Y, X₁, …, X_N) = 0, где Y — выходная величина модели измерения, являющаяся одновременно измеряемой величиной, значение которой должно быть получено на основе информации о входных величинах X₁, …, X_N.

Примечание 3. Если модель измерения содержит две и более выходные величины, то она включает в себя более одного уравнения.

3.8. Многомерная модель (измерения) (multivariate measurement model): модель измерения с произвольным числом выходных величин.

Примечание 1. В общем случае многомерная модель измерения имеет вид уравнений

h₁(Y₁, …, Y_m, X₁, …, X_N) =

= 0, …, h_m(Y₁, …, Y_m, X₁, …, X_N) = 0,

где Y₁, …, Y_m — m выходных величин, в совокупности составляющих измеряемую величину, значения которых должны быть получены на основе информации о входных величинах многомерной модели X₁, …, X_N.

Примечание 2. Общий вид многомерной модели измерения может быть представлен также в векторной форме

где и — матрицы размерности m x 1.

Примечание 3. В случае одной выходной величины, т.е. m = 1 (см. примечание 1), модель измерения называют одномерной.

3.9. Многомерная функция (измерения) (multivariate measurement function): функция, определяющая зависимость выходных величин от входных величин в многомерной модели измерения.

Примечание 1. Данное определение модифицировано по отношению к JCGM 200 (словарная статья 2.49).

Примечание 2. Если уравнения, входящие в модель измерения , могут быть разрешены в явном виде , где — входные величины, а — выходные величины модели измерения, то — многомерная функция измерения. В более общем случае под можно понимать алгоритм, посредством которого устанавливается однозначное соответствие значений выходных величин значениям входных величин .

Примечание 3. В случае одной выходной величины, т.е. m = 1 (см. примечание 2), функцию измерения называют одномерной.

3.10. Модель (измерения) с действительными величинами (real measurement model): модель измерения (в общем случае многомерная), в состав которой входят только действительные величины.

3.11. Модель (измерения) с комплексными величинами (complex measurement model): модель измерения (в общем случае многомерная), в состав которой входят комплексные величины.

3.12. Модель многоступенчатого измерения (multistage measurement model): модель измерения (в общем случае многомерная), состоящая из последовательности подмоделей, связанных между собой таким образом, что выходные величины подмодели одной ступени являются входными величинами подмодели следующей ступени.

Примечание. Зачастую потребность в определении области охвата для выходных величин (на основе их совместного распределения) имеет место только на заключительном этапе измерения.

Пример — Измерение, включающее в себя процедуру калибровки, может рассматриваться как двухступенчатое. Для первой подмодели значениями входных величин являются передаваемые от эталонов и соответствующие им показания средства измерений, а выходными величинами — параметры калибровочной функции (градуировочной характеристики). Эта подмодель определяет способ определения выходных величин по входным величинам, например решением системы уравнений, получаемых при применении метода наименьших квадратов. Входными величинами второй подмодели являются параметры калибровочной функции и новое показание средства измерений, а выходной величиной — измеряемая величина, для получения значения которой было применено средство измерений.

3.13. Функция (совместного) распределения (вероятностей) (joint distribution function): функция, дающая для каждого значения значение вероятности того, что каждый элемент X_i случайной векторной переменной будет меньше или равен .

Примечание. Функцию распределения случайной переменной обозначают , где

3.14. Плотность (совместного) распределения (вероятностей) (joint probability density function): неотрицательная функция , удовлетворяющая условию

3.15. Маргинальная плотность распределения (вероятностей) (marginal probability density function): плотность распределения элемента X_i случайной векторной переменной с плотностью совместного распределения , которая имеет вид

Примечание. Если все элементы X_i, i = 1, …, N, составляющие случайную переменную , независимы, то .

3.16. Математическое ожидание (expectation): характеристика случайной величины X_i, являющейся элементом случайной векторной переменной с плотностью совместного распределения , которая имеет вид

Примечание 1. Данное определение модифицировано по отношению к JCGM 101 (словарная статья 3.6).

Примечание 2. Математическим ожиданием случайной векторной величины является — матрица размерности N x 1.

3.17. Дисперсия (variance): характеристика случайной величины X_i, являющейся элементом случайной векторной переменной с плотностью совместного распределения , которая имеет вид

Примечание. Данное определение модифицировано по отношению к JCGM 101 (словарная статья 3.7).

3.18. Ковариация (covariance): характеристика двух случайных величин X_i и X_j, являющихся элементами случайной векторной переменной с плотностью совместного распределения , которая имеет вид

где — плотность совместного распределения случайных величин X_i и X_j.

Примечание 1. Данное определение модифицировано по отношению к JCGM 101 (словарная статья 3.10).

Примечание 2. Ковариационной матрицей случайной векторной величины является симметричная положительно полуопределенная матрица размерности N x N, элементами которой являются ковариации Cov(X_i, X_j), i = 1, …, N, j = 1, …, N. Некоторые операции с использованием налагают более строгое ограничение в виде положительной определенности этой матрицы.

3.19. Корреляция (correlation): характеристика двух случайных величин X_i и X_j, являющихся элементами случайной векторной переменной с плотностью совместного распределения , которая имеет вид

Примечание. Величина Corr(X_i, X_j) имеет размерность единица.

3.20. Ковариационная матрица (оценок) (measurement covariance matrix): связанная с оценкой действительной векторной величины размерности N x 1 симметричная положительно полуопределенная матрица размерности N x N, на главной диагонали которой расположены квадраты стандартных неопределенностей, соответствующих оценкам элементов векторной величины, а остальные члены матрицы представляют собой ковариации между парами соответствующих оценок элементов векторной величины.

Примечание 1. Термин и определение модифицированы по отношению к JCGM 101 (словарная статья 3.11).

Примечание 2. Ковариационная матрица размерности N x N, соответствующая вектору оценок векторной величины , имеет вид

где u(x_i, x_i) = u²(x_i) — дисперсия (квадрат стандартной неопределенности) оценки x_i;

u(x_i, x_j) — ковариация между x_i и x_j. Если элементы X_i и X_j вектора некоррелированны, то u(x_i, x_j) = 0.

Примечание 3. В JCGM 101 ковариационная матрица называется матрицей неопределенности.

Примечание 4. При работе с ковариационными матрицами могут возникать некоторые вычислительные трудности. Например, ковариационная матрица , соответствующая оценке , может не быть положительно определенной (это зависит от того, каким образом была рассчитана матрица ). Как следствие, для такой матрицы не будет существовать разложение Холецкого, часто применяемое в численных методах вычислений (см. [7] и Приложение B). Более того, дисперсия для линейной комбинации элементов , которая предположительно должна иметь небольшое положительное значение, может оказаться отрицательной. Для таких ситуаций разработаны методы «коррекции» , после применения которых полученная матрица будет положительно определена, и, соответственно, для нее будет существовать разложение Холецкого, а дисперсия линейной комбинации элементов будет всегда положительна. Один из таких методов приведен в [27], а его принцип состоит в следующем. Выполняют спектральное разложение матрицы , представляя ее в виде

где — матрица, столбцы которой являются ортонормированными собственными векторами матрицы , а — диагональная матрица, на главной диагонали которой расположены соответствующие собственные значения . Строят новую диагональную матрицу , заменяя в матрице элементы, меньшие, чем d_min, на d_min, где d_min равно произведению наибольшего элемента на единичную ошибку округления компьютера, применяемого при вычислениях. Тогда «корректированная» ковариационная матрица, применяемая для последующих вычислений, будет иметь вид

Примечание 5. Некоторые операции с использованием требуют, чтобы данная матрица была положительно определенной.

3.21. Корреляционная матрица (оценок) (correlation matrix): связанная с оценкой действительной векторной величины размерности N x 1 симметричная положительно полуопределенная матрица размерности N x N, членами которой являются корреляции между парами соответствующих оценок элементов векторной величины.

Примечание 1. Корреляционная матрица размерности N x N, соответствующая вектору оценок векторной величины , имеет вид:

где r(x_i, x_i) = 1, а r(x_i, x_j) — корреляция между x_i и x_j. Если элементы X_i и X_j вектора некоррелированны, то r(x_i, x_j) = 0.

Примечание 2. r(x_i, x_j) называют также коэффициентом корреляции.

Примечание 3. Корреляционная матрица и ковариационная матрица (см. 3.20) связаны между собой соотношением

где — диагональная матрица размерности N x N с диагональными элементами u(x₁), …, u(x_N). Элементы матрицы могут быть представлены в виде

u(x_i, x_j) = r(x_i, x_j)u(x_j)u(x_j).

Примечание 4. Корреляционная матрица будет положительно определенной/сингулярной в том и только в том случае, если соответствующая ей ковариационная матрица будет положительно определенной/сингулярной. Некоторые операции с использованием требуют, чтобы данная матрица была положительно определенной.

Примечание 5. При представлении численных значений недиагональных элементов корреляционной матрицы часто достаточно округлять их с точностью до трех знаков после запятой. Однако если корреляционная матрица близка к сингулярной, то, чтобы избежать вычислительных сложностей при использовании корреляционной матрицы среди прочих исходных данных в оценивании неопределенности измерения, число сохраняемых десятичных знаков необходимо увеличить. Это число зависит от характера последовательных вычислений, но в качестве ориентировочного значения рекомендуется брать его равным числу десятичных знаков, необходимых для представления наименьшего собственного значения корреляционной матрицы с двумя значимыми десятичными знаками. Так для корреляционной матрицы размерности 2 x 2 собственные значения и равны соответственно 1 + |r| и 1 — |r|, где r — недиагональный элемент корреляционной матрицы, и, значит, таким наименьшим собственным значением будет 1 — |r|. Если заранее известно, что корреляционная матрица является сингулярной, то округление к меньшему по модулю снижает риск того, что после операции округления корреляционная матрица не окажется положительно полуопределенной.

3.22. Матрица (коэффициентов) чувствительности (sensitivity matrix): матрица частных производных первого порядка функций, описывающих модель измерения с действительными величинами, по входным или входным величинам в точке оценок этих величин.

Примечание. В случае модели с N входными и m выходными величинами матрицы чувствительности в отношении входных величин и выходных величин имеют размерности соответственно m x N и m x m.

3.23. Интервал охвата (coverage interval): интервал, построенный на основе имеющейся информации и содержащий значение скалярной случайной величины с заданной вероятностью.

Примечание 1. Данное определение модифицировано по отношению к JCGM 101 (словарная статья 3.12).

Примечание 2. Вероятностно симметричный интервал охвата для скалярной величины представляет собой интервал охвата, для которого вероятность того, что значение случайной величины меньше наименьшего значения (нижней границы) интервала охвата, равна вероятности того, что значение случайной величины больше наибольшего значения (верхней границы) интервала [см. JCGM 101 (словарная статья 3.15)].

Примечание 3. Наименьший интервал охвата представляет собой интервал охвата, имеющий наименьшую длину среди всех возможных интервалов охвата для данной случайной величины с одинаковой вероятностью охвата [см. JCGM 101 (словарная статья 3.16)].

3.24. Область охвата (coverage region): область, определенная на основе имеющейся информации и содержащая значение векторной случайной величины с заданной вероятностью.

3.25. Вероятность охвата (coverage probability): вероятность того, что значение случайной величины находится в границах интервала охвата или области охвата.

Примечание 1. Данное определение модифицировано по отношению к JCGM 101 (словарная статья 3.13).

Примечание 2. Вероятность охвата иногда называют уровнем доверия [JCGM 100 (6.2.2)].

3.26. Наименьшая область охвата (shortest coverage region): область охвата, имеющая наименьший объем среди всех возможных областей охвата для данной случайной величины с одинаковой вероятностью охвата.

Примечание. В случае скалярной случайной величины наименьшая область охвата совпадает с наименьшим интервалом охвата. Для случайной величины, описываемой вектором в двумерном пространстве, наименьшая область охвата представляет собой поверхность с наименьшей площадью из всех, имеющих ту же вероятность охвата.

3.27. Многомерное нормальное распределение (вероятностей) <1> (multivariate Gaussian distribution): распределение вероятностей векторной случайной величины размерности N + 1, плотность совместного распределения которого имеет вид:

Примечание. — математическое ожидание , — ковариационная матрица , которая должна быть положительно определена.

———————————

<1> Многомерное нормальное распределение называют также многомерным распределением Гаусса.

3.28. Многомерное t-распределение (multivariate t-distribution): распределение вероятностей векторной случайной величины размерности N x 1, плотность совместного распределения которого с параметрами , и имеет вид:

где — гамма-функция, z > 0.

Примечание 1. Многомерным t-распределением описывается векторная случайная величина размерности N x 1, удовлетворяющая соотношению , где — векторная случайная величина размерности N x 1, имеющая нормальное распределение с нулевым математическим ожиданием и положительно определенной ковариационной матрицей размерности N x N, а W — скалярная случайная величина, имеющая (распределение хи-квадрат) с степенями свободы.

Примечание 2. Плотность t-распределения нельзя представить в виде произведения N плотностей распределения элементов вектора даже в том случае, когда — диагональная матрица. В общем случае между элементами вектора существует статистическая зависимость. Например, при N = 2, v = 5 и — единичной матрице размерности 2 x 2 вероятность того, что X₁ > 1 составляет 18%, в то время как условная вероятность того, что при X₂ > 2 значение X₁ будет превышать единицу, составляет 26%.

4. Соглашения и условные обозначения

В настоящем стандарте использованы следующие соглашения и условные обозначения.

4.1. В GUM [JCGM 100 (пункт 4.1.1, примечание 1)] для экономии условных обозначений один и тот же символ (прописная буква) используется для:

(i) физической величины, которая предполагает наличие единственного истинного значения;

(ii) случайной переменной, ассоциированной с этой физической величиной.

Примечание. Случайная переменная выполняет разные роли при оценивании неопределенности по типу A и B. При оценивании неопределенности по типу A случайная переменная представляет собой «… возможный результат наблюдения величины». При оценивании неопределенности по типу B вероятность распределения случайной переменной характеризует имеющиеся знания о возможных значениях этой величины.

Эта двойственность обозначений в большинстве случаев не вызывает неудобств.

В настоящем стандарте (так же, как и в JCGM 101) в случае входных величин, неопределенность которых оценивают по типу A, один и тот же символ (прописная буква) использован для трех понятий, а именно:

a) физическая величина;

b) случайная переменная, для которой получают результаты наблюдений;

c) случайная переменная, распределение вероятности которой ассоциируют с имеющимися знаниями о возможных значениях физической величины.

Два последних понятия, относящиеся к случайной величине, в GUM (JCGM 100) не разделяются, что может явиться источником недоразумений. Так рассматриваемая в настоящем стандарте и в JCGM 101 процедура оценивания неопределенности с использованием метода Монте-Карло может быть неправильно истолкована как реализация процедуры, изложенной в JCGM 100 (пункт 4.1.4, примечание 1). В действительности же, хотя указанные процедуры схожи в том, что в обеих получают выборку значений выходной величины для данной модели измерения из соответствующего распределения, сами распределения в общем случае будут разными. В JCGM 100 (пункт 4.1.4, примечание 1) это частотное распределение, т.е. случайная величина интерпретируется в смысле перечисления b), тогда как в методе Монте-Карло это распределение случайной величины, интерпретируемой в смысле перечисления c). Для большинства измерительных задач подход, предложенный в JCGM 100 (пункт 4.1.4, примечание 1), не рекомендуется (см. [2]).

4.2. Для входных величин модели измерения в настоящем стандарте принято обозначение X₁, …, X_N или в виде матрицы размерности N x 1 (символ «^T» обозначает транспонирование).

4.3. Для выходных величин модели измерения в настоящем стандарте принято обозначение Y₁, …, Y_m или в виде матрицы размерности m x 1.

4.4. Если Y_i могут быть выражены через в явном виде, то модель измерения имеет вид

(1)

где — многомерная функция измерения.

Другая форма записи для той же модели (см. 3.9) имеет вид

где являются составляющими .

4.5. Если Y_i не выражены в явном виде через , то модель измерения имеет вид

(2)

или в другой форме записи (см. 3.8)

4.6. Оценку обозначают в виде — матрицы размерности N x 1. Ковариационную матрицу, соответствующую , обозначают в виде — матрицы размерности N x N (см. 3.20).

4.7. Оценку обозначают в виде — матрицы размерности m x 1. Ковариационную матрицу, соответствующую , обозначают в виде — матрицы размерности m x m.

Примечание. в случае многомерной модели с m выходными величинами является аналогом дисперсии u²(y) для y в случае одномерной модели измерения, рассматриваемой в JCGM 100 и JCGM 101. В JCGM 100 u(y) обозначается как u_c(y), где подстрочный индекс «c» применительно к стандартной неопределенности обозначает «суммарная». Как и в JCGM 101, в настоящем стандарте использование подстрочного индекса «c» в данном контексте рассматривается как излишнее [см. JCGM 101 (пункт 4.10)].

4.8. Если оценки выходных величин предполагается использовать по отдельности, то каждая из этих величин может рассматриваться как выходная в соответствующей одномерной модели измерения. Если же, например, для последующих расчетов эти оценки должны быть использованы совместно, то должны быть приняты во внимание корреляции между ними.

4.9. Стандартную неопределенность, соответствующую x, обозначают u(x). Если контекст исключает возможность ошибочного истолкования, то может применяться сокращенная форма записи u_x. Данная форма записи не рекомендуется, если при x имеется индекс или иной знак, например x_i или .

4.10. Под можно понимать как «оценки входных величин», так и «оценку входной величины (векторной)». В настоящем стандарте преимущественно используется последнее определение (то же самое справедливо для выходных величин).

4.11. Как указано в 4.2 — 4.10, величина в общем случае обозначается с помощью прописной буквы, а ее оценка или некоторое фиксированное значение величины (такое, как математическое ожидание) соответствующей строчной буквой. Данное правило удобно для общего анализа, но зачастую не подходит для обозначения величин в конкретных приложениях из-за устоявшейся практики использования для конкретных физических величин специальных обозначений, например T для температуры и t для времени. Поэтому в некоторых примерах настоящего стандарта используются иные обозначения: физическая величина обозначается ее общепринятым символом, а ее математическое ожидание или оценка этим же символом с циркумфлексом («крышкой»). Например, амплитуда переменного тока (пример 1 из 6.2.2) обозначается I, а оценка I — [см. JCGM 101 (пункт 4.8)].

4.12. Настоящий стандарт отступает от обозначений, часто используемых для обозначения плотностей распределения вероятностей и функций распределения. В JCGM 100 одно и то же обозначение f использовано как для функции измерения, так и для плотности распределения вероятностей, что создает неоднозначность понимания. Поскольку в настоящем стандарте моделям уделено особое внимание, для плотности распределения вероятностей и функции распределения вместо обозначений f и F использованы соответственно g и G. Применяемые подстрочные индексы соответствуют случайной величине, о которой идет речь. Обозначение f оставлено для описания функции измерения (в скалярной или векторной форме).

4.13. Плотность распределения может быть поставлена в соответствие как скалярной (X), так и векторной величине. В случае скалярной величины плотность распределения для X обозначается как , где — переменная, принимающая возможные значения величины X. Здесь X рассматривается как случайная переменная с математическим ожиданием E(X) и дисперсией V(X).

4.14. В случае векторных величин плотность распределения для обозначается как , где — переменная, принимающая возможные значения величины . Здесь рассматривается как случайная переменная с ожиданием и ковариационной матрицей .

4.15. Аналогично, в случае скалярных величин (Y) плотность распределения обозначается как , а в случае векторных величин .

4.16. Для обозначения десятичной дроби используется запятая <1>.

———————————

<1> В оригинале на английском языке в данном подразделе указывается на использование в качестве десятичного знака точки вместо запятой.

5. Основные принципы

5.1. Общие положения

5.1.1. В GUM [JCGM 100 (пункт 4.1)] измерение моделируется функцией, связывающей действительные входные величины X₁, …, X_N и действительную выходную величину Y в виде формулы (1), т.е. , где — действительная векторная входная величина. Это одномерная функция измерения для действительных величин (см. 3.9, примечание 3).

5.1.2. На практике не все измерения могут быть смоделированы с помощью функции измерения с одной скалярной выходной величиной. В реальных измерительных задачах могут иметь место:

a) несколько выходных величин Y₁, …, Y_m (которые совместно обозначаются действительной векторной выходной величиной , для которых формула (1) принимает вид ;

b) более общий вид модели измерения в виде формулы (2), т.е. .

5.1.3. Кроме того, некоторые или все элементы и, соответственно, элементы могут представлять собой комплексные величины. Если каждую такую комплексную величину представить в виде двух составляющих (действительная и мнимая часть или модуль и аргумент комплексного числа), то, в принципе, без нарушения общности модель измерения может рассматриваться как модель с действительными величинами. Однако в большинстве случаев вид алгоритмов, работающих с комплексными величинами, проще, чем если бы модель включала только действительные величины [14]. Применение моделей измерения с комплексными величинами позволяет записать закон трансформирования неопределенностей в компактном матричном виде (см. 6.4 и Приложение A).

5.1.4. В настоящем стандарте модели, указанные в 5.1.2 и 5.1.3, рассматриваются в более общем виде.

5.2. Основные этапы оценивания неопределенности

5.2.1. Основные этапы оценивания неопределенности включают в себя формулировку измерительной задачи, трансформирование распределений и получение окончательного результата:

a) формулировка измерительной задачи включает в себя:

1) задание выходной величины (измеряемой векторной величины);

2) выявление входных величин, составляющих векторную входную величину , от которых зависит ;

3) составление модели измерения, определяющей взаимосвязь с в виде функции измерения [см. формулу (1)] или в более общем виде [см. формулу (2)];

4) приписывание распределений вероятностей (нормального, прямоугольного и т.д.) входным величинам X_i (элементам вектора ) или совместного распределения вероятностей входным величинам, не являющимся независимыми, на основе имеющейся о них информации,

b) трансформирование распределений предусматривает определение плотности совместного распределения выходной величины на основе плотностей распределения входных величин X_i и используемой модели измерения,

c) получение окончательного результата предполагает использование плотности распределения для определения:

1) оценки математического ожидания в виде ;

2) ковариационной матрицы , соответствующей ;

3) области охвата, содержащей с заданной вероятностью p (вероятность охвата).

5.2.2. Формулировку измерительной задачи осуществляет метролог. Рекомендации по выбору плотности распределения [стадия 4) этапа a) в 5.2.1] для некоторых общих случаев приведены в JCGM 101 и в 5.3. Этапы трансформирования распределений и получения окончательных результатов [б) и в) в 5.2.1], для которых приведены подробные указания, не требуют дополнительной метрологической информации и могут быть выполнены с любой требуемой вычислительной точностью для поставленной задачи.

Примечание. Как только этап постановки задачи a) в соответствии с 5.2.1 выполнен, тем самым плотность распределения вероятностей для выходной величины формально полностью определена. Однако вычисление математического ожидания, стандартного отклонения и области охвата может потребовать применения численных методов, обладающих некоторой степенью приближения.

5.3. Функции плотности вероятности для входных величин

5.3.1. Общие положения

Руководство по выбору плотностей распределения для входных величин X_i на этапе формулировки измерительной задачи приведено в JCGM 101 (раздел 6) для некоторых общих случаев. Однако единственным многомерным распределением, рассмотренным в JCGM 101, является многомерное нормальное распределение JCGM 101 (пункт 6.4.8). Это распределение приписывают входной величине , если доступная информация об включает в себя только оценку и соответствующую ковариационную матрицу . В 5.3.2 рассматривается еще одно многомерное распределение — t-распределение. Его применяют, если единственной доступной информацией о величине является выборка наблюдений (предполагаемых независимыми) векторной величины из многомерного нормального распределения с неизвестными математическим ожиданием и ковариационной матрицей (см. также 6.5.4).

5.3.2. Многомерное t-распределение

5.3.2.1. Предположим, что для векторной величины размерностью N x 1, имеющей многомерное нормальное распределение с неизвестными математическим ожиданием и ковариационной матрицей размерностью N x N, доступны n независимых наблюдений, n > N. Пусть — искомое значение . Тогда, выбирая в качестве априорных распределений для и соответствующие неинформативные распределения и используя теорему Байеса, получим, что совместным распределением для (или распределением, приписываемым ) будет многомерное t-распределение с степенями свободы [11], где

Примечание. При наличии соответствующих оснований в качестве априорных распределений могут быть взяты другие распределения, что может привести к другому значению числа степеней свободы для или даже к другому типу распределения для .

5.3.2.2. Плотность распределения, полученного для , имеет вид

где — гамма-функция аргумента z.

5.3.2.3. Математическим ожиданием и ковариацией будут соответственно

где определено только для (что соответствует n > N + 1).

5.3.2.4. Чтобы сформировать случайное выборочное значение из , возьмем N случайных выборочных значений Z_i, i = 1, …, N, из стандартного распределения Гаусса N(0, 1) и одно выборочное значение w из с степенями свободы. Тогда

где — нижняя треугольная матрица размерности N x N в разложении Холецкого [13].

Примечание. Матрица может быть определена, например, как в [13].

5.3.3. Построение многомерных функций плотности распределения

Когда входные величины X₁, …, X_N коррелированны, то обычно доступной о них информацией является вид плотности распределения для каждой из этих величин (например, для одной — нормальное, для другой — прямоугольное и т.п.), оценки x₁, …, x_N, используемые в качестве математических ожиданий, стандартные неопределенности u(x₁), …, u(x_N), используемые в качестве стандартных отклонений, и ковариации, соответствующие парам x_i. Построить по маргинальным распределениям X₁, …, X_N совместную плотность распределения для можно, зная их копулу. Однако вышеуказанной исходной информации может соответствовать множество копул, поэтому вид построенной совместной плотности распределения будет не единственным.

5.4. Трансформирование распределений

5.4.1. В левой части рисунка 1 показан пример модели измерения с N = 3 взаимно независимыми входными величинами и m = 2 выходными величинами . Функция измерения — . Величинам X_i, i = 1, 2, 3, приписаны плотности распределения , а характеризуется совместной плотностью распределения . В правой части рисунка 1 показан пример, в котором X₁ и X₂ взаимно зависимы и характеризуются совместной плотностью распределения 

Рисунок 1. Трансформирование распределений для модели

с N = 3 входными величинами и m = 2 выходными величинами,

когда входные величины X₁, X₂ и X₃ взаимно независимы

(слева) и когда X₁ и X₂ взаимно зависимы (справа)

5.4.2. Выходная величина может сама служить основой для получения следующей величины, например, . Тогда будет рассматриваться как входная величина в модели измерения, описываемой, например, функцией измерения и имеющей вид

Так, может представлять собой набор эталонов массы, а — суммы некоторых из них.

5.4.3. Объединение функций измерения и для двух подмоделей позволяет получить зависимость непосредственно от входных величин . Однако в ряде измерительных задач желательно сохранить разбиение на подмодели, если они относятся к функционально разным этапам. Совокупность двух подмоделей представляет собой пример модели многоступенчатого измерения (см. 3.12).

5.4.4. Случай, когда на финальном этапе многоступенчатого измерения с применением многомерных подмоделей имеется единственная выходная скалярная величина, может быть рассмотрен с применением JCGM 101.

5.5. Получение итоговой информации

5.5.1. Оценка выходной величины рассматривается как математическое ожидание . Ковариационная матрица , соответствующая , — как ковариационная матрица .

5.5.2. Для вероятности охвата p область охвата R_Y для получают решением уравнения

Примечание 1. Некоторым величинам могут быть поставлены в соответствие случайные переменные с распределениями, у которых математического ожидания и ковариационной матрицы не существует (см., например, 5.3.2). Однако область охвата для существует всегда.

Примечание 2. В общем случае существует более одной области охвата для заданной вероятности охвата p.

5.5.3. Прямого многомерного аналога вероятностно симметричного 100p%-ного интервала охвата, рассмотренного в JCGM 101, не существует. Однако существует аналог наименьшего 100p%-ного интервала охвата — это 100p%-ная наименьшая область охвата.

5.6. Способы трансформирования распределений

5.6.1. Трансформирование распределений осуществляют несколькими способами:

a) аналитическими методами, обеспечивающими получение математического представления плотности распределения для ;

b) применением закона трансформирования неопределенностей, основанного на замене функции измерения ее аппроксимацией рядом Тейлора с членами разложения первого порядка [обобщение подхода, изложенного в JCGM 100 (пункт 5.1.2)];

c) численными методами [см. JCGM 100 (пункт G.1.5)], в том числе с использованием метода Монте-Карло (ММК).

Примечание 1. Аналитические методы превосходят все прочие с той точки зрения, что они не используют приближений. Однако они применимы только в простых случаях. Такие методы в настоящем стандарте не рассматриваются, за исключением примеров, где они используются для сравнения.

Примечание 2. Метод Монте-Карло в настоящем стандарте используется для получения распределения векторной выходной величины, а не в качестве метода имитационного моделирования. При оценивании неопределенности на этапе трансформирования распределений решаемая задача является детерминированной, поэтому в имитационном моделировании физического случайного процесса нет необходимости.

5.6.2. В законе трансформирования неопределенностей оценка для и соответствующая ковариационная матрица подвергаются преобразованию посредством линеаризованной модели измерения. В настоящем стандарте данная процедура рассматривается для моделей разных типов.

5.6.3. В левой части рисунка 2 показан обобщенный закон трансформирования неопределенностей для модели измерения с N = 3 взаимно независимыми входными величинами и m = 2 выходными величинами . Оценкой является с соответствующими стандартными неопределенностями u(x₁), u(x₂) и u(x₃). Оценкой является с соответствующей ковариационной матрицей . В правой части рисунка 2 тот же закон показан для случая, когда X₁ и X₂ взаимно зависимы и имеют ковариацию u(x₁, x₂) оценок x₁ и x₂.

Рисунок 2. Обобщенный закон трансформирования

неопределенностей для N = 3 взаимно независимых величин

X₁, X₂ и X₃ и m = 2 взаимно зависимых выходных

величин (слева) и тот же закон,

но для взаимно зависимых и X₁, и X₂ (справа)

5.6.4. В методе Монте-Карло совместное распределение вероятностей для , представленное в цифровом виде, трансформируется с помощью модели измерения для того, чтобы получить дискретное представление совместного распределения вероятности для , на основе которого затем получают окончательные результаты измерения.

6. Способ оценивания неопределенности по GUM

6.1. Общие положения

6.1.1. В настоящем стандарте способ оценивания неопределенности через трансформирование неопределенностей, рассмотренный в JCGM 100 (пункты 6.2 и 6.3) для моделей измерения вида , обобщен на более широкий класс моделей с многими выходными переменными. Хотя непосредственно в JCGM 100 такие модели не рассматриваются, для их изучения могут быть применены те же самые основные принципы трансформирования оценок входных величин и соответствующих им неопределенностей в оценки выходных величин и соответствующих им неопределенностей. Для математического представления указанных процедур вместо сумм величин с подстрочными индексами, как это сделано в JCGM 100, удобнее использовать компактную матрично-векторную форму записи, хорошо приспособленную для современных пакетов программ и языков программирования.

6.1.2. Для применения закона трансформирования неопределенностей используется та же информация о входных величинах, что и для одномерной модели измерения, рассмотренной в JCGM 100:

a) оценка входной величины ;

b) ковариационная матрица , соответствующая , содержащая ковариации u(x_i, x_j), i = 1, N, j = 1, …, N, соответствующие x_i и x_j.

6.1.3. Описание трансформирования неопределенностей, приведенное в 6.2 и 6.3, распространяется на модели с действительными величинами, включая случаи комплексных величин, представленных парами действительных составляющих. Трансформирование неопределенностей в случае моделей с комплексными величинами рассматривается в 6.4 (см. также 5.1.3).

6.1.4. Способ получения области охвата для векторной выходной величины описан в 6.5.

6.2. Трансформирование неопределенностей для многомерных моделей измерения с явным видом функциональной зависимости

6.2.1. Общие положения

6.2.1.1. Многомерная модель измерения с явным видом функциональной зависимости между выходной величиной и входной величиной имеет вид

где обозначает многомерную функцию измерения.

Примечание. Аргументами отдельных функций могут быть разные подмножества . При этом каждый элемент должен являться аргументом как минимум одной функции .

6.2.1.2. При заданной оценке для оценка для имеет вид

6.2.1.3. Ковариационная матрица размерности m x m, соответствующая , имеет вид

где Cov(y_j, y_j) = u²(y_j), и определяется по формуле

(3)

где — матрица чувствительности размерности m x N, определяемая по формуле

где все производные берутся в точке ([19], страница 29).

6.2.2. Примеры

Пример 1 — Активное и реактивное сопротивления элемента цепи [JCGM (раздел H.2)]

Активное R и реактивное X сопротивления элемента цепи определяют путем измерения амплитуды V изменяющегося по гармоническому закону напряжения на его клеммах, амплитуды I проходящего через элемент переменного тока и фазового сдвига между напряжением и силой тока. Двумерная модель измерения для R и X, выраженных через V, I и , имеет вид

(4)

В обозначениях, принятых в настоящем стандарте, N = 3, m = 2, и .

Оценку активного и реактивного сопротивлений получают по формуле (4) в точке — оценке входной величины .

Ковариационную матрицу размерности 2 x 2, соответствующую y, рассчитывают по формуле (3), где — матрица чувствительности размерности 2 x 3, получаемая вычислением

в точке , а — ковариационная матрица размерности 3 x 3, соответствующая .

Примечание. В JCGM 100 реактивное сопротивление обозначено X. Это обозначение использовано и в настоящем стандарте. Реактивное сопротивление X является элементом векторной выходной величины , и его не следует путать с векторной входной величиной .

Пример 2 — Коэффициент отражения, измеряемый микроволновым рефлектометром (подход 1)

Комплексный коэффициент отражения , измеряемый микроволновым рефлектометром, например, таким, который используют для определения повреждения кабельных линий, описывается моделью с комплексными величинами в виде

(5)

где — комплексный неисправленный коэффициент отражения, а , и — комплексные коэффициенты, полученные при градуировке (калибровке) рефлектометра [10, 16, 26].

В обозначениях настоящего стандарта, описывая комплексные величины через их действительные и мнимые части, получаем N = 8, m = 2, и .

Оценку комплексного коэффициента отражения в виде его действительной и мнимой частей вычисляют по формуле (5), подставляя в нее оценку входной величины .

Ковариационную матрицу размерности 2 x 2, соответствующую , рассчитывают по формуле (3), где — матрица чувствительности размерности 2 x 8, получаемая при вычислении производных

в точке , а — ковариационная матрица размерности 8 x 8, соответствующая .

Пример 3 — Калибровка эталонов массы

Этот пример описывает модель многоступенчатого измерения (см. 3.12, 5.4.2 и 5.4.3).

Набор из q эталонов массы со значениями калибруют сличением с эталоном килограмма с использованием компаратора массы, калибровочной гири для определения калибровочной функции компаратора и ряда вспомогательных приборов, таких как термометр, барометр и гигрометр для определения поправок на выталкивающую силу воздуха. Эталон килограмма и гиря имеют массы m_R и m_S соответственно. Калибровку проводят в соответствии с подходящей методикой измерений посредством достаточного числа k сличений между наборами эталонов с получением видимых, т.е. заметных при измерениях в воздухе, разностей . Вычисляют соответствующие поправки на выталкивающую силу воздуха . Разности масс в вакууме получают из подмодели , в которой .

Оценку масс обычно получают решением по методу наименьших квадратов переопределенной системы уравнений , где — матрица размерности k x q с элементами, равными единице, минус единице или нулю в соответствии с тем, какие эталоны массы включены в сличение, с учетом неопределенностей, соответствующих оценке величины . В этом случае формула для определения имеет вид

(6)

где ковариационную матрицу размерности q x q, соответствующую , получают по формуле , а — ковариационная матрица размерности k x k, соответствующая . Более подробное описание подмодели, а также процедура получения через — ковариационной матрицы, соответствующей оценке величины , — приведено в [3].

Многомерная модель измерения для этого примера имеет вид

где — функция измерения. В принятых обозначениях настоящего стандарта N = k, m = q и .

Примечание. С вычислительной точки зрения для получения оценки предпочтительнее использовать не формулу (6), а алгоритм, основанный на ортогональном разложении матриц (см. [13]).

6.3. Трансформирование неопределенностей для многомерных моделей измерения с неявным видом функциональной зависимости

6.3.1. Общие положения

6.3.1.1. Многомерная модель измерения с неявным видом функциональной зависимости между выходной величиной и входной величиной имеет вид

6.3.1.2. При заданной оценке величины оценку величины получают решением системы уравнений

(7)

Примечание. Систему уравнений (7) относительно обычно решают численными методами, например методом Ньютона [12] или одной из его модификаций, задавая начальное значение корня и последовательно приближаясь к решению.

6.3.1.3. Ковариационную матрицу размерности m x m, соответствующую , получают решением системы уравнений

(8)

где — матрица чувствительности размерности m x m, содержащая частные производные , l = 1, …, m, j = 1, …, m, а — матрица чувствительности размерности m x N, содержащая частные производные , l = 1, …, m, i = 1, …, N. Производные вычисляются в точках и .

Примечание 1. Ковариационная матрица в формуле (8) не определена, если матрица является вырожденной (сингулярной).

Примечание 2. Формулу (8) получают аналогично формуле (3) с использованием правила дифференцирования неявной функции.

6.3.1.4. Из формулы (8) следует, что решение относительно ковариационной матрицы может быть записано в виде

(9)

где

(10)

матрица размерности m x N, сформированная из коэффициентов чувствительности.

6.3.1.5. Процедура расчета матрицы приведена в Приложении B. Применение для этих целей непосредственно формул (9) и (10) не рекомендуется вследствие неустойчивости соответствующих им алгоритмов вычислений.

6.3.2. Примеры

Пример 1 — Давления, задаваемые грузопоршневым манометром

Давление p, задаваемое грузопоршневым манометром, определяют в соответствии с уравнением преобразования, имеющим вид

(11)

где m_w — полная приложенная масса (груза и поршня), и — плотности воздуха и приложенного груза соответственно, g_l — локальное значение ускорения свободного падения, A₀ — эффективная площадь манометра при нулевом давлении, — коэффициент деформации поршневой пары манометра, — коэффициент теплового расширения, — отклонение от нормальных условий по температуре (20 °C) [17].

Пусть p₁, …, p_q обозначают давления уравновешивания для приложенных масс, соответственно, m_{w, 1}, …, m_{w, q} и температурных отклонений .

В обозначениях, принятых в настоящем стандарте, N = 6 + 2q, m = q, , .

Модель измерения, определяющая вид зависимости между и , имеет вид

j = 1, …, q. (12)

Оценку величины p_j получают решением уравнения (12) при заданных оценках A₀, , , , m_{w, j}, , и g_l. Однако полученные оценки , j = 1, …, q, имеют соответствующие ковариации, т.к. все они зависят от одних и тех же случайных величин A₀, , , , и g_l.

Ковариационную матрицу размерности q x q, соответствующую , вычисляют по формуле (8), где — матрица чувствительности размерности q x q, содержащая частные производные , l = 1, …, q, j = 1, …, q; — матрица чувствительности размерности q x (6 + 2q), содержащая частные производные , l = 1, …, q, i = 1, …, 6 + 2q, (все производные берут в точках и ), а — ковариационная матрица размерности (6 + 2q) x (6 + 2q), соответствующая .

Примечание 1. В данном примере выражение зависимости Y_j (или, что то же самое, p_j) через может быть задано в явном виде как решение квадратного уравнения. Однако числовой алгоритм нахождения корня квадратного уравнения не всегда будет устойчив. Более того, иногда в уравнение преобразования включают дополнительные члены, представляющие собой степени p более высоких порядков. В таких случаях получение явной функции измерения не всегда возможно.

Примечание 2. Рассматриваемая в данном примере модель измерения может быть представлена разными способами. Например, вместо зависимости, описываемой формулой (12), может быть использована модель в виде сравнения с нулем разности между левой и правой частями уравнения (11). От выбора модели зависит эффективность и устойчивость численного решения.

Примечание 3. Могут быть рассмотрены более полные модели давления для грузопоршневого манометра, которые включают, например, поправки, учитывающие эффекты поверхностного натяжения.

Примечание 4. Функции измерения имеют в качестве своих аргументов не все входные величины. Так, в выражение для j-й функции входят только A₀, , , , m_{w, j}, , и g_l.

Пример 2. Коэффициент отражения, измеренный микроволновым рефлектометром (подход 2)

Другой подход к задаче, описанной в примере 2 из 6.2.2, заключается в выражении зависимости между входной величиной и выходной величиной через двумерную модель измерения, имеющую вид

(13)

где и соответственно действительная и мнимая части выражения

Преимущество этого подхода состоит в том, что вычисление производных и, следовательно, коэффициентов чувствительности производится более прямым способом.

Оценку комплексного коэффициента отражения находят в результате подстановки в формулы (13) и численного решения полученных уравнений.

Ковариационную матрицу размерности 2 x 2, соответствующую , вычисляют по формуле (8), где — матрица чувствительности размерности 2 x 2, содержащая частные производные , l = 1, 2, j = 1, 2; — матрица чувствительности размерности 2 x 8, содержащая частные производные , l = 1, 2, i = 1, …, 8 (все производные вычисляют в точках и ), а — ковариационная матрица размерности 8 x 8, соответствующая .

Пример 3 — Калибровка рефлектометра

Калибровку рефлектометра (см. пример 2 из 6.2.2) обычно проводят, измеряя неисправленный коэффициент отражения при применении эталонов с заданными значениями коэффициента отражения . Часто в этих целях используют три эталона, что позволяет получить систему из трех совместных уравнений:

j = 1, 2, 3. (14)

Разделение выражений в левой части уравнения (14) на действительную и мнимую части приведет к получению шести совместных линейных уравнений, решение которых позволяет найти действительную и мнимую части коэффициентов , и калибровочной функции при заданных значениях действительной и мнимой частей неисправленных коэффициентов отражения и коэффициентов отражения для эталонов.

В обозначениях, принятых в настоящем стандарте, N = 12, m = 6, и .

Входные и выходные величины связаны между собой посредством многомерной модели измерения, в которой и , j = 1, 2, 3, — соответственно действительная и мнимая части левой части уравнения (14).

Оценку калибровочных коэффициентов получают, подставляя оценки для и в уравнения (14) и решая эти уравнения численно.

Ковариационную матрицу размерности 6 x 6, соответствующую , вычисляют по формуле (8), где — матрица чувствительности размерности 6 x 6, содержащая частные производные , l = 1, …, 6, j = 1, …, 6; — матрица чувствительности размерности 6 x 12, содержащая частные производные , l = 1, …, 6, i = 1, …, 12 (все производные вычисляют в точках и ), а — ковариационная матрица размерности 12 x 12, соответствующая .

Примечание 1. При наличии программы обработки данных, работающей с комплексными величинами, разделение уравнений модели измерения на действительную и мнимую части необязательно. Эти уравнения могут быть решены непосредственно для , и .

Примечание 2. Каждое j-е уравнение модели измерения включает только четыре входные величины: W_{j, R}, W_{j, I}, и .

6.4. Трансформирование неопределенности для моделей с комплексными величинами

В Приложении A приведен компактный алгоритм вычисления частных производных многомерных комплексных функций измерения первого порядка, которые необходимо знать при распространении закона трансформирования неопределенностей на модели с комплексными величинами. Данный алгоритм может быть применен для многомерных моделей измерения с комплексными величинами общего вида.

Пример — Коэффициент отражения, измеряемый микроволновым рефлектометром (подход 3)

Рассмотрим вновь пример 2 из 6.2.2.

Комплексная выходная величина и комплексная входная величина связаны между собой моделью измерения, описываемой формулой (5). С учетом результатов, приведенных в Приложении A, матрицу чувствительности размерности 2 x 8 получают в виде

где

в точке оценки для величины . Например, для

(15)

использование результатов Приложения A дает

где Q_R и Q_I — соответственно действительная и мнимая части выражения в правой части формулы (15).

Ковариационную матрицу размерности 2 x 2, соответствующую , где

вычисляют по формуле (A.1) Приложения A, где — ковариационная матрица размерности 8 x 8, соответствующая .

6.5. Область охвата векторной выходной величины

6.5.1. Общие положения

6.5.1.1. В некоторых областях метрологии, например, связанных с измерениями электрических величин, для дальнейшего использования результатов измерения выходную величину удобно оставить в векторной форме вместе с поставленным ей в соответствие совместным распределением вероятностей. Такое представление результата измерения позволяет в максимальной степени сохранить всю полученную информацию о выходной величине.

6.5.1.2. Если же результат измерения представляют в виде полученной оценки выходной величины , соответствующей ей ковариационной матрицы и вероятности охвата p, то это требует определения области охвата R_y в m-мерном пространстве, которая содержит с вероятностью p.

6.5.1.3. Если доступная информация о выходной величине содержит только значения оценок и , то в соответствии с принципом максимума энтропии случайной величине, ассоциированной с измеряемой величиной и характеризующей имеющиеся на данный момент представления о ее возможных значениях, приписывают многомерное нормальное распределение [см. JCGM 101 (пункт 6.4.8)].

Примечание. Такой выбор функции распределения согласуется с используемым в способе оценивания неопределенности по GUM представлением о нормальном распределении скалярной выходной величины Y в случаях, когда число степеней свободы, связанное с оценкой Y, бесконечно.

6.5.1.4. В общем случае, как только получено совместное распределение для выходной величины , появляется возможность определить вероятность охвата для некоторой заданной области охвата или, наоборот, построить область охвата для заданной вероятности охвата. Эта задача не вызывает затруднений, если совместное распределение является многомерным нормальным распределением (см. 6.5.2, 6.5.3 и 6.5.4). Для других распределений возможно получение приближенных решений с точностью, приемлемой для практических задач, благодаря использованию численных методов, таких как метод Монте-Карло (см. раздел 7).

6.5.1.5. В 6.5.2 рассматривается метод определения области охвата для двумерной величины, что потом позволит распространить его на общий многомерный случай (см. 6.5.3). Также рассматривается определение области охвата для случая, когда оценка выходной величины получена усреднением результатов наблюдений этой величины, представляемых как случайная независимая выборка из многомерного нормального распределения (см. 6.5.4).

6.5.2. Двумерный случай

6.5.2.1. На примере двумерной модели измерения можно продемонстрировать все отличия в определении многомерной области охвата от получения одномерного интервала охвата. Рассмотрим точку в прямоугольной системе координат, где Y₁ является абсциссой точки, а Y₂ ее ординатой. Пусть измерение обеих координат проведено с использованием одного средства измерений. Полученная информация об будет включать в себя оценки y₁ и y₂ ее координат, стандартные неопределенности u(y₁) и u(y₂), соответствующие этим оценкам, и ковариацию u(y₁, y₂), наличие которой в данном случае обусловлено использованием одного и того же средства измерений для каждой координаты.

6.5.2.2. Согласно способу оценивания неопределенности по GUM при наличии информации о выходной величине , указанной в 6.5.2.1, ей приписывают плотность распределения в виде двумерного нормального распределения (см. 6.5.1.3), где

6.5.2.3. Из возможных форм областей охвата рассматриваются две:

a) эллипс с центром в точке [19, страница 38], описываемый формулой

(16)

где k_p — постоянная, которую находят из условия, что интеграл от по площади, ограниченной эллипсом, равен p. При данной форме области охвата в учет принимается взаимная зависимость между Y₁ и Y₂. Если характеризируется нормальным распределением, то величина

(17)

имеет распределение хи-квадрат с двумя степенями свободы. Отсюда следует, что представляет собой квантиль уровня p указанного распределения, т.е.

где имеет распределение хи-квадрат с двумя степенями свободы. Для вероятности охвата p = 0,95 k_p = 2,45 (см. 6.5.3);

b) прямоугольник с центром в точке со сторонами, параллельными осям и равными длинам по отдельности определенных интервалов охвата для Y₁ и Y₂. Интервал охвата для Y₁ находят из маргинального распределения Y₁, определяемого формулой

(18)

и он не зависит от имеющейся информации о величине Y₂. То же самое справедливо для интервала охвата для Y₂. Интервалы охвата y_j +/- k_qy(y_j), j = 1, 2, определяют для вероятности охвата

q = 1 — (1 — p)/2 = (1 + p)/2. (19)

Данный прямоугольник представляет собой область охвата для , соответствующий вероятности охвата, не меньшей p [5]. Если характеризуется нормальным распределением, то маргинальное распределение (18) для Y₁, равно как и распределение для Y₂, также будет нормальным. Отсюда следует, что k_q представляет собой квантиль уровня q стандартного нормального распределения, т.е.

q = Pr(|Z| <= k_q),

где Z имеет нормальное распределение N(0,1). Для вероятности охвата p = 0,95 q = 0,975 и k_q = 2,24 (см. 6.5.3).

Примечание 1. При выполнении условий применимости способа оценивания неопределенностей по GUM область охвата в виде эллипса, указанная в перечислении a), является наименьшей 100p%-ной областью охвата.

Примечание 2. Если Y₁ и Y₂ взаимно независимы, то q в формуле (19) может быть заменено на q = p^1/2.

Пример 1 — Рассмотрим двумерную величину , характеризуемую нормальным распределением , где

На рисунке 3 слева показаны 95%-ные эллиптическая и прямоугольная области охвата для , полученные методами a) и b). Также на рисунке показано 1000 точек, представляющих случайную выборку из указанного распределения вероятности. Интервал охвата в форме эллипса будет иметь наименьшую площадь при заданной вероятности охвата. Для данной случайной выборки 950 из 1000 точек содержатся внутри эллиптической области охвата площадью 26,6 квадратных единиц и 953 — внутри прямоугольной области охвата площадью 28,4 квадратных единиц.

Рисунок 3. Эллиптические и прямоугольные области охвата

для двумерной величины , описываемой нормальным

распределением, элементы которой Y₁ и Y₂ взаимно

независимы (слева) и коррелированны (справа)

Пример 2 — Рассмотрим двумерную величину , характеризуемую нормальным распределением , где

На рисунке 3 справа показаны 95%-ные эллиптическая и прямоугольная области охвата для , полученные методами a) и b). В отличие от примера 1, элементы Y₁ и Y₂ величины являются коррелированными. Для случайной выборки из указанного распределения 957 из 1000 точек содержатся внутри эллиптической области охвата площадью 11,8 квадратных единиц и 972 — внутри прямоугольной области охвата площадью 40,1 квадратных единиц, что показывает, что вероятность охвата для прямоугольной области превышает 0,95. Область в форме прямоугольника, построенная согласно методу b) без учета коррелированности элементов выходной величины и распределения точек выборки на плоскости, может рассматриваться как представление области охвата для , не соответствующее данной измерительной задаче. Прямоугольник со сторонами, параллельными осям эллипса, будет иметь меньшую площадь и мог бы рассматриваться как более подходящая область охвата, но неудобство его применения состоит в том, что для него теряется ясный физический смысл соответствия формы построенной области охвата элементам выходной величины.

6.5.2.4. Другие примеры областей охвата для выходных величин двумерной модели измерения приведены в разделе 9.

6.5.3. Многомерный случай

Если размерность вектора выходной величины более двух, то способы построения областей охвата становятся менее наглядными, но по существу ничем не отличаются от описанных в перечислениях a) и b) в 6.5.2.3 для двумерного случая. Необходимо построить такую область в m-мерном пространстве, которая содержала бы с вероятностью p. Как и в 6.5.2.3 [перечисления a) и b)], рассматриваются формы двух видов:

a) эллипсоид в m-мерном пространстве, описываемый формулой (16). Если характеризируется многомерным нормальным распределением, то величина, задаваемая формулой (17), имеет распределение хи-квадрат с m степенями свободы. Таким образом, представляет собой квантиль уровня p указанного распределения, т.е.

где имеет распределение хи-квадрат с m степенями свободы. Значения коэффициента охвата k_p для вероятности охвата p = 0,95 и разных значений m приведены в таблице 1;

b) параллелепипед в m-мерном пространстве с центром в точке с ребрами, параллельными осям элементов Y_j, j = 1, …, m, выходного вектора в прямоугольной системе координат и равными длинам по отдельности определенных интервалов охвата для Y_j. Интервалы охвата определяют для вероятности охвата q = 1 — (1 — p)/m. Данный параллелепипед представляет собой область охвата для , соответствующий вероятности охвата, не меньшей p [5]. Для каждого элемента Y_j, j = 1, …, m, интервал охвата вычисляют через его маргинальное распределение. Если характеризуется нормальным распределением, то маргинальное распределение для каждого Y_j также будет нормальным. Отсюда следует, что, как и в 6.5.2.3 [перечисление b)], k_q представляет собой квантиль уровня q стандартного нормального распределения. Значения коэффициента охвата k_q для вероятности охвата p = 0,95 и разных значений m приведены в таблице 2.

Таблица 1

Коэффициенты охвата для областей охвата в форме m-мерного

эллипсоида, соответствующие вероятности охвата p = 0,95

Таблица 2

То же, что в таблице 1, но для областей охвата

в форме m-мерного параллелепипеда

Примечание 1. Для одномерного случая (m = 1) выражение (16) упрощается до

что дает граничные точки интервала охвата для Y

Для вероятности охвата p = 0,95 k_p = 1,96 (см. таблицу 1).

Примечание 2. Если Y_j, j = 1, …, m, взаимно независимы, то q может быть заменено на q = p^1/m.

6.5.4. Область охвата для оценки в виде выборочного среднего из многомерного нормального распределения

Рассмотрим n векторов размерностью m x 1 каждый, n > m, соответствующих повторяющимся наблюдениям многомерной величины . Предположим, что могут быть интерпретированы как реализации независимых случайных векторов , каждый из которых характеризуется многомерным нормальным распределением с математическим ожиданием и ковариационной матрицей . Определим среднее и ковариационную матрицу

размерности m x 1 и m x m соответственно. Тогда случайная переменная

будет иметь распределение Фишера F_{m, n — m} (называемое также F-распределением) с m и n — m степеней свободы [19, пункт 3.5.2.1].

Примечание. Аналогом этого результата для одномерной величины будет следующее утверждение: для независимых случайных переменных Y₁, …, Y_n, каждая из которых характеризуется одномерным нормальным распределением с математическим ожиданием и дисперсией , величина (n — 1)^1/2 имеет t-распределение с n — 1 степенями свободы, где

Пример — Рассмотрим n = 12 повторных парных наблюдений объемных долей микроклина (A₁) и биотита (A₂) в одном тонком разрезе гранита G-2 [4, 25]. На рисунке 4 значения парных наблюдений изображены в виде точек. На нем также построена 95%-ная эллиптическая область охвата для математического ожидания величины A размерности 2 x 1. Оценка математического ожидания и соответствующая ей ковариационная матрица имеют вид

а 95-я перцентиль распределения F_{2, 10} равна 4,10. 95%-ная область охвата для A — это эллипс, описываемый формулой

Небольшое число наблюдений в данном примере не позволяет сделать содержательные выводы о том, насколько справедливы исходные допущения, чтобы считать построенную область охвата достоверной.

Рисунок 4. Двенадцать пар выборочных значений и 95%-ная

эллиптическая область охвата вокруг их выборочного среднего

7. Метод Монте-Карло

7.1. Общие положения

7.1.1. В настоящем разделе рассматривается применение метода Монте-Карло для трансформирования распределений (соответствующая процедура описана в 7.1.7 и представлена в виде диаграммы на рисунке 5).

Рисунок 5. Этапы трансформирования распределений

и получения результатов оценивания неопределенности

методом Монте-Карло для случая явной зависимости

выходных величин от входных величин

7.1.2. Метод Монте-Карло позволяет реализовать общий подход к получению дискретного приближенного представления функции распределения для [18, страница 75]. Суть подхода состоит в получении повторных выборок из плотности распределения для X_i (или совместной плотности распределения для ) и вычислении для каждого выборочного значения векторной выходной величины.

7.1.3. Поскольку содержит максимально полную информацию о , любые характеристики , такие как математическое ожидание, дисперсия и ковариация, а также области охвата могут быть рассчитаны из полученного приближения . В общем случае достоверность получаемых характеристик возрастает с увеличением числа выборок.

7.1.4. Полученные в соответствии с 7.1.2 значения выходной величины рассматриваются как независимая выборка из совместного распределения вероятности для . Математические ожидания, дисперсии (и высшие моменты), а также ковариации могут быть рассчитаны непосредственно по этим выборочным значениям. Определение областей охвата требует предварительного анализа полученных значений (см. 7.7).

7.1.5. Пусть , r = 1, …, M, обозначает выборочные значения выходной величины (см. 7.1.4). Выборка позволяет получить приближенные значения математического ожидания и дисперсии величины . Как правило, в качестве моментов величины [включая и ] принимают соответствующие выборочные моменты. Обозначим число векторов в выборке , для которых каждый их элемент не превосходит соответствующий элемент некоторого вектора размерности m x 1. Тогда вероятность может быть приближенно представлена отношением . Таким образом, выборка векторов позволяет получить дискретное представление функции распределения .

7.1.6. Приближение является первым результатом применения метода Монте-Карло и представляет собой матрицу размерности m x M

7.1.7. Процедура применения метода Монте-Карло для трансформирования распределений в случае явной зависимости через и заранее заданного числа испытаний M (в противном случае см. 7.8) показана в виде диаграммы на рисунке 5 и включает в себя следующие этапы:

a) выбирают число испытаний M (см. 7.2);

b) формируют в каждом из M испытаний N-мерный вектор входных величин X_i, элементами которого являются случайные выборочные значения из распределений для X_i или совместного распределения для (см. 7.3);

c) рассчитывают для каждого выборочного значения вектора входной величины вектор выходной величины , получая таким образом выборку векторов выходной величины объемом M (см. 7.4);

d) формируют представление функции распределения в виде ряда M значений векторной выходной величины (см. 7.5);

e) на основе вычисляют оценку величины и ковариационную матрицу , соответствующую (см. 7.6);

f) на основе строят соответствующую область охвата для для заданной вероятности охвата p (см. 7.7).

Примечание. Выборочное среднее для M векторных выходных величин имеет математическое ожидание и дисперсию . Таким образом, расхождение между и его оценкой в среднем будет пропорционально M^-1/2.

7.1.8. Эффективность метода Монте-Карло при определении , и области охвата для зависит от адекватного выбора числа испытаний M [этап a) в 7.1.7]. Рекомендации по выбору достаточного числа испытаний M и по другим вопросам реализации метода Монте-Карло приведены в [8] (см. также 7.2 и 7.8).

7.2. Число испытаний метода Монте-Карло

7.2.1. Для применения метода Монте-Карло необходимо выбрать число испытаний M, т.е. объем выборки векторной выходной величины. Это число может быть выбрано заблаговременно (до проведения испытаний), но тогда будет исключена возможность управления точностью результатов, полученных с помощью данного метода. Причиной этому служит то, что число испытаний, необходимое для получения результата вычисления с заданной точностью, зависит от формы плотности распределения выходной величины и от заданного значения вероятности охвата. Кроме того, метод вычисления является стохастическим по своей природе, поскольку зависит от случайной выборки.

7.2.2. Поскольку нельзя заранее гарантировать, что выбранное значение M обеспечит достаточную точность приближения, можно использовать процедуру адаптивного выбора, уточняя значение M в процессе испытаний. Адаптивная процедура, установленная в 7.8, позволяет оптимальным образом получить значение M, соответствующее заданной точности вычислений.

Примечание. Для сложной модели, например, требующей получения решения методом конечных элементов, применение большого числа испытаний может оказаться невозможным. В этом случае рекомендуется приближенно представить плотность распределения выходной величины нормальным распределением (как в GUM). Это позволяет использовать относительно небольшое число испытаний M, например 50 или 100, а полученные по результатам испытаний выборочное среднее и выборочные ковариации принять, соответственно, в качестве оценок и . Для описания и построения области охвата используют плотность нормального распределения . Хотя уменьшение числа испытаний неизбежно ухудшает свойства метода в части аппроксимации распределения выходной величины, оно все же позволяет учесть нелинейность модели измерения.

7.3. Получение выборок из распределений вероятности

7.3.1. Для применения метода Монте-Карло формируют M векторов , r = 1, …, M, в соответствии с плотностями распределения вероятностей для N входных величин X_i или, если это необходимо, из совместной плотности распределения ФПВ .

7.3.2. Рекомендации по формированию выборки для наиболее распространенных распределений (равномерного, нормального, многомерного нормального и t-распределения) приведены в JCGM 101 (пункт 6.4 и приложение C). Однако возможно получение выборок, соответствующих и другим распределениям [см. JCGM 101 (раздел C.2)]. Некоторые распределения могут быть аппроксимированы распределениями, полученными в результате применения метода Монте-Карло при предыдущих вычислениях неопределенности [см. JCGM 101 (пункт 6.5 и приложение D)].

7.3.3. Процедура формирования выборки для многомерного t-распределения описана в 5.3.2.4.

Примечание. Для достоверности результатов применения метода Монте-Карло необходимо, чтобы генераторы псевдослучайных чисел, используемые для формирования выборок из заданных распределений, обладали соответствующими свойствами. В JCGM 101 (пункт C.3.2) приведены некоторые тесты сформированных выборок на случайность.

7.4. Вычисление векторной выходной величины

7.4.1. Векторную выходную величину вычисляют для каждого из M выборочных значений N-мерного вектора входной величины, полученных из соответствующих распределений. Если обозначить M выборочных значений вектора входной величины размерности N x 1 через , где r-й вектор содержит x_{1, r}, …, x_{r, N}, и x_{i, r} — выборочное значение из распределения для X_i, то соответствующие значения векторов выходной величины для модели измерения с явным видом функциональной зависимости вычисляют по формуле

r = 1, …, M.

Примечание. При использовании закона трансформирования неопределенностей, когда аналитические выражения производных функции измерения по входным величинам известны точно, значения этих производных и значения функции измерения берут в точке оценок входных величин. Если аналитические выражения для производных неизвестны и для их оценок используют приближение в виде конечных разностей, то получают оценки только для функции измерения. Согласно рекомендации GUM [JCGM 100 (примечание 2 к пункту 5.1.3)] значения функции измерения берут в точках оценок входных величин, а также в точках, отстоящих по обе стороны от этих оценок на расстоянии одной стандартной неопределенности (варьируя по очереди для каждой входной величины). В методе же Монте-Карло значения функций измерения получают при варьировании входных величин в окрестности их оценок, т.е. в точках, которые могут отстоять от этих оценок на несколько стандартных отклонений. Поскольку в методе Монте-Карло значения функции измерения получают в разных точках, может возникнуть вопрос о свойствах вычислительной процедуры, в частности, о ее устойчивости и (в случае применения адаптивной процедуры) сходимости. При возникновении сомнений пользователю следует убедиться в том, что метод дает достоверные оценки выходной величины для достаточно больших окрестностей оценок входных величин. Однако следует ожидать, что вопросы устойчивости и сходимости численного метода могут стать критическими только в исключительных случаях.

7.4.2. Если X_i являются зависимыми величинами, то в 7.4.1 распределения для X_i следует заменить на совместное распределение для .

7.4.3. В случае модели с неявным видом функциональной зависимости в процедуре, описанной в 7.4.1, значения векторной выходной величины получают в результате решения уравнений

r = 1, …, M.

7.5. Дискретное представление функции распределения для выходной величины

Дискретное представление функции распределения для векторной выходной величины формируется из M значений векторной выходной величины, полученных согласно 7.4. В общем случае это представление является матрицей размерности m x M, r-й столбец которой является r-м значением векторной выходной величины. Для одномерной модели — вектор-строка.

Примечание 1. j-я строка матрицы является приближенным дискретным представлением функции маргинального распределения для Y_j.

Примечание 2. Для m = 1, 2, 3 возможна графическая интерпретация приближения . Пусть m = 2. Для r = 1, …, M построим в плоскости (Y₁, Y₂) точку с координатами, соответствующими двум элементам в столбце r матрицы . При достаточно большом M плотность нанесенных точек в любой локальной области плоскости будет приблизительно пропорциональна плотности распределения в этой области.

Примечание 3. Построенное приближение позволяет извлечь разнообразную информацию, в частности, относительно моментов высших порядков. Однако в способе оценивания неопределенности по GUM для трансформирования неопределенностей и последующего оценивания неопределенности измерения [JCGM 100 (пункт 0.4)] для получения оценки и соответствующей ей ковариационной матрицы необходимо знать моменты только первых двух порядков.

Примечание 4. Если величина будет использоваться в качестве входной величины на следующем этапе многоступенчатого измерения, то на этом этапе выборку входных величин можно получить из уже имеющихся значений , r = 1, …, M (или, что то же самое, из столбцов матрицы ), последовательно осуществляя случайный выбор (с равной вероятностью) из этих значений [JCGM 101:2008, 6.5].

7.6. Оценка выходной величины и соответствующей ей ковариационной матрицы

Среднее арифметическое и ковариационную матрицу

принимают, соответственно, в качестве оценки величины и ковариационной матрицы , соответствующей .

Примечание. Независимо от того, является ли модель измерения линейной или нелинейной, при M, стремящемся к бесконечности, стремится по вероятности к , если последнее существует.

7.7. Область охвата для векторной выходной величины

7.7.1. Общие положения

Теоретически существует сколь угодно много 100p%-ных областей охвата для . В настоящем стандарте рассматриваются области охвата трех видов, которые, как ожидается, могут найти наибольшее практическое применение. Однако в отдельных обстоятельствах может оказаться предпочтительным использование области охвата специальной формы. Независимо от вида области охвата построение последней основано на полученном приближении , т.е. на наборе M точек , случайным образом выбранных из распределения вероятностей для при реализации метода Монте-Карло (см. 7.5). Рассматриваемые здесь области охвата имеют вид:

a) эллипсоида. Данная область охвата будет близка к наименьшей области охвата в случае, если распределение вероятностей для близко к нормальному;

b) параллелепипеда. Область охвата данного вида допускает простую интерпретацию, но зачастую ее объем существенно превосходит объем наименьшей области охвата;

c) наименьшей области охвата, которая в общем случае не имеет какой-либо заданной геометрической формы и определяется в некотором приближении, зависящем от M.

7.7.2. Область охвата в форме эллипсоида

Уравнение 100p%-ной области охвата для в форме эллипсоида имеет вид:

(20)

где определяет положение центра эллипсоида в пространстве выходных величин, — форму эллипсоида (степень вытянутости и ориентацию в пространстве), а k_p — его размер (объем). Процедура построения области охвата в форме эллипсоида по формуле (20) при известных оценках и (см. 7.6) требует нахождения k_p и состоит в следующем:

a) трансформируют точки в (r = 1, …, M) по формуле

r = 1, …, M, (21)

где — нижняя треугольная матрица размерности m x m разложения Холецкого ;

b) сортируют трансформированные точки в порядке возрастания значения d_r, где

r = 1, …, M;

c) используют упорядоченные для определения коэффициента охвата k_p, значение которого определяют из условия: d_r < k_p для всех с номером r, не превышающим pM;

d) строят 100p%-ную область охвата для , границы которой определены уравнением (20).

Примечание 1. Основы описанной процедуры изложены в [1], где массив векторных данных сортируют по метрике

где — статистика параметра положения, а — статистика параметра масштаба (дисперсии). Формулы, приведенные в процедуре, получены заменами и . Использование в процедуре трансформирования точек необходимо для того, чтобы сделать ковариационную матрицу, соответствующую трансформированным точкам, единичной, т.е. элементы векторов выборки становятся некоррелированными. Это позволяет сопоставлять (сортировать) трансформированные точки по их расстоянию от оценки . Из способа получения точек следует, что они характеризуют разброс значений случайной величины вокруг центральной точки . Область охвата в виде эллипсоида с центром в по своему смыслу должна содержать 100p% общего числа точек .

Примечание 2. Приемлемость полученной области охвата для представления неопределенности измерения зависит от конкретной измерительной задачи. Полученное представление может стать неприемлемым, если распределение точек в пространстве выходной величины плохо согласуется с описанием этой величины многомерным нормальным распределением.

Примечание 3. Матрица , полученная на этапе a), может быть определена, например, способом, описанным в [13] (см. пример ниже).

Примечание 4. Процедура построения области охвата учитывает существование зависимости между элементами вектора .

Примечание 5. Для сложной модели, например, требующей получения решения методом конечных элементов, применение большого числа испытаний M может оказаться невозможным. В этом случае рекомендуется приближенно представить плотность распределения выходной величины нормальным распределением. Это позволяет использовать относительно небольшое число испытаний M, например 50 или 100, а полученные по результатам испытаний выборочное среднее и выборочные ковариации принять, соответственно, за оценки и . Для описания и построения области охвата используют плотность нормального распределения с математическим ожиданием и ковариационной матрицей . Хотя уменьшение числа испытаний неизбежно ухудшает свойства метода в части аппроксимации распределения выходной величины, оно все же позволяет учесть нелинейность модели измерения. В таблице 1 приведены значения k_p для вероятности охвата p = 0,95 в зависимости от размерности m вектора выходных величин при допущении нормальности распределения .

Пример — Рассмотрим модель измерения

Y₁ = X₁ + X₃, Y₂ = X₂ + X₃, (22)

в которой входным величинам X₁ и X₂ приписано нормальное распределение N(0; 0,1) каждой, X₃ приписано прямоугольное распределение R(-(3 x 1,9)^1/2, (3 x 1,9)^1/2), и все входные величины являются независимыми. Математическими ожиданиями входных величин X_i являются x_i = 0, i = 1, 2, 3, а их дисперсиями — u²(x_i) = 0,1, i = 1, 2 и u²(x₃) = 1,9. Тогда, как в примере 2 из 6.5.2.3, имеет оценку математического ожидания и соответствующую ей ковариационную матрицу

Трансформированные точки на этапе a) формируют с использованием L^-1, имеющей следующий вид (с точностью до третьего знака после запятой):

На рисунке 6 слева показаны 1000 выборочных точек для распределения вероятности , определяемого моделью измерения [формула (22)] и приведенными выше распределениями вероятности для входных величин X_i, i = 1, 2, 3. 95%-ная эллиптическая область охвата для , полученная в предположении, что характеризуется двумерным распределением Гаусса [как и на этапе a) в 6.5.2.3], показана пунктирной линией. Эта область имеет площадь 11,8 квадратных единиц, k_p = 2,45 и содержит 968 из 1000 точек. Эллиптическая область охвата, определенная на основе 1000 выборочных точек в соответствии с процедурой 7.7.2, показана сплошной линией. Эта область имеет площадь 10,6 квадратных единиц, k_p = 2,33, что немногим меньше, чем область, полученная в предположении нормальности распределения выходной величины. Как и должно следовать из способа построения области охвата, эта область содержит ровно 950 точек.

Рисунок 6. Эллиптические области охвата, построенные

в соответствии с процедурами 6.5.2.3 a) и 7.7.2

для примера из 7.7.2 (слева) и прямоугольные области

охвата, построенные в соответствии с процедурами

6.5.2.3 b) и 7.7.3 для примера из 7.7.3 (справа)

Модель измерения [формула (22)] рассматривается более подробно в разделе 9, где приведены также другие примеры построения областей охвата для двумерных выходных величин.

7.7.3. Область охвата в форме параллелепипеда

100p%-ная область охвата для в форме параллелепипеда имеет вид

y_j +/- k_qu(y_j), j = 1, …, m, (23)

где определяет положение центра параллелепипеда в пространстве выходных величин, а и k_q определят его размер (объем). Процедура построения области охвата в форме параллелепипеда по формуле (23) при известных оценках и требует нахождения k_q и состоит в следующем:

a) трансформируют точки в (r = 1, …, M) по формуле (21), но где теперь является диагональной матрицей размерности m x m с диагональными элементами u(y₁), …, u(y_m);

b) сортируют трансформированные точки в порядке возрастания значения d_r, где d_r теперь определяют по формуле

r = 1, …, M;

c) используют упорядоченные для определения коэффициента охвата k_q, значение которого определяют из условия: d_r < k_q для всех с номером r, не превышающим pM;

d) строят 100p%-ную область охвата для , границы которой определены формулой (23).

Примечание. Процедура построения области охвата учитывает существование зависимости между элементами вектора .

Пример — На рисунке 6 (справа) показаны те же 1000 точек, что и на рисунке 6 (слева). 95%-ная прямоугольная область охвата для , полученная в предположении, что Y_j, j = 1, 2, характеризуется нормальным распределением N(y_j; u²(y_j)), а корреляция оценок y₁ и y₂ [как и на этапе b) в 6.5.2.3] не учитывается, показана пунктирной линией. Эта область имеет площадь 40,1 квадратных единиц, k_q = 2,24 и содержит все 1000 точек. Прямоугольная область охвата, определенная на основе 1000 выбранных точек в соответствии с процедурой 7.7.3, показана сплошной линией. Эта область имеет площадь 25,5 квадратных единиц, k_q = 1,78, что немногим меньше, чем область в предположении нормальности распределения выходной величины. Как и должно следовать из способа построения области охвата, эта область содержит ровно 950 точек.

7.7.4. Наименьшая область охвата

Процедура построения наименьшей 100p%-ной области охвата включает в себя следующие этапы:

a) в m-мерном пространстве выходных величин строят исходный параллелепипед, предположительно охватывающий большую долю значений, которые могут быть приписаны векторной выходной величине ;

b) равномерной m-мерной сеткой делят исходный параллелепипед на ячейки в виде подобных равных плотно уложенных малых параллелепипедов;

c) для каждого малого параллелепипеда подсчитывают число попавших в него выборочных точек ;

d) долю точек , попавших в малый параллелепипед, рассматривают как приближенное значение вероятности попадания в данный параллелепипед значения случайной векторной величины ;

e) все малые параллелепипеды располагают в порядке уменьшения приписанной им вероятности;

f) суммируют вероятности, приписанные малым параллелепипедам, начиная с первого и последовательно прибавляя вероятность каждого последующего вплоть до того, пока полученная сумма не превысит или не станет равной p;

g) совокупность малых параллелепипедов, внесших вклад в получение указанной суммы, рассматривают в качестве приближения наименьшей области охвата.

Примечание 1. Данная процедура, заимствованная из [20], заключается в разбиении пространства выходных величин на ряд малых ячеек (параллелепипедов), аппроксимации вероятности попадания значения случайной выходной величины в каждую из таких ячеек долей попавших в нее точек выборки и объединении минимального числа ячеек, в совокупности содержащих не менее 100p% из M выборочных значений векторной выходной величины, в конфигурацию, рассматриваемую как наименьшая область охвата.

Примечание 2. Параллелепипед, построенный на этапе a), должен включать в себя все точки .

Примечание 3. Число малых параллелепипедов, используемых для разбиения на этапе b) пространства выходных величин, ограниченного большим параллелепипедом, построенным на этапе a), влияет на качество полученного приближения.

Примечание 4. В большинстве случаев качество приближения улучшается с ростом M. Чтобы получить достаточно хорошее приближение, особенно в случае, когда размерность пространства выходных величин m более двух или трех, может потребоваться очень большое число испытаний M.

Примечание 5. Построенная в соответствии с вышеописанной процедурой наименьшая область охвата может оказаться несвязной, особенно, если M недостаточно велико.

Примечание 6. На этапе d) вероятность, приписанная каждому малому параллелепипеду, определяется относительной частотой попадания в него выборочного значения выходной величины. Использование более сложной аппроксимации (см. [23]) может позволить улучшить связность области охвата и сгладить ее границы.

Примечание 7. Приписанные на этапе d) вероятности могут оказаться одинаковыми для двух или нескольких малых параллелепипедов. В этом случае упорядочение малых параллелепипедов, выполняемое на этапе e), будет не единственным. Разным возможным вариантам упорядочения могут соответствовать разные приближения области наименьшего охвата.

Примечание 8. В двумерном случае (m = 2) этапы от a) до d) входят в начальные шаги типичного алгоритма визуализации распределения для (см. 9.1.6).

Пример — Рассмотрим ту же самую задачу, что и в примере 7.7.2. На рисунке 7 показано приближение для наименьшей 95%-ной области охвата, полученное с использованием вышеуказанной процедуры при разбиении исходного большого прямоугольника на 10 x 10 малых прямоугольников (слева) и 100 x 100 малых прямоугольников (справа). Область охвата на рисунке 7 слева построена по 1000 точкам, случайно выбранным из распределения для . Она имеет площадь 11,3 квадратных единиц и содержит 955 точек. Область охвата на рисунке 7 справа построена по 1000000 точек, полученных случайной выборкой из распределения для . Она имеет площадь 9,4 квадратных единиц и содержит 950074 точки. Для сравнения на каждом рисунке сплошной линией показана 95%-ная эллиптическая область охвата для , построенная в соответствии с процедурой из 7.7.2.

Рисунок 7. Эллиптические области охвата, построенные

по процедуре из 7.7.2, и приближения наименьших областей

охвата, построенные по процедуре из 7.7.4 для разбиения

10 x 10 и 1000 точек выборки (слева) и для разбиения

100 x 100 и 1000000 точек выборки (справа)

7.8. Адаптивная процедура метода Монте-Карло

7.8.1. Общие положения

7.8.1.1. Эффективность метода Монте-Карло при определении оценки выходной величины , соответствующей ковариационной матрицы и области охвата для зависит от числа испытаний M. Значение M может быть выбрано заранее, как в 7.2. Другой подход состоит в использовании адаптивной процедуры метода Монте-Карло, согласно которой число испытаний последовательно увеличивают до тех пор, пока полученные числовые оценки искомых величин не станут установившимися. Численный результат считается установившимся, если соответствующее ему удвоенное стандартное отклонение станет меньше заданной точности вычисления (см. 7.8.2).

7.8.1.2. Целью адаптивной процедуры, описанной в 7.8.3, является получение в среднем с заданной вычислительной точностью значений следующих величин:

a) оценки величины ;

b) вектора стандартных неопределенностей, связанных с оценками;

c) положительно определенной матрицы размерности m x m коэффициентов корреляции r_ij = r(y_i, y_j), связанных с парами оценок;

d) коэффициента охвата k_p, определяющего 100p%-ную область охвата для в форме m-мерного эллипсоида.

Примечание 1. То, что выполнение требования к точности вычислений может быть гарантировано не безусловно, а только в среднем, обусловлено природой случайности, используемой в методе Монте-Карло.

Примечание 2. Как правило, чем больше вероятность охвата p, тем большее число испытаний метода Монте-Карло необходимо для определения k_p с заданной точностью вычислений.

Примечание 3. Стандартные неопределенности и корреляционная матрица вместе определяют ковариационную матрицу , соответствующую оценкам (см. примечание 3 к 3.2.1 и 7.8.2.4).

Примечание 4. Если требуется построить область охвата иной формы, нежели эллипсоид в m-мерном пространстве, то проверку устойчивости метода выполняют для параметров, характеризующих область охвата данной формы. Например, для области охвата в форме m-мерного параллелепипеда оценка точности вычисления может быть выполнена для коэффициента охвата k_q. В этом случае процедура, описанная в 7.8.3, должна быть модифицирована соответствующим образом.

Примечание 5. Если в представлении результата измерения не требуется указывать область охвата, то процедуру увеличения числа испытаний останавливают после получения установившихся численных значений оценки , соответствующих стандартных неопределенностей и корреляционной матрицы .

Примечание 6. Матрица может не быть положительно определенной в случае существования зависимостей между выходными величинами (для примера см. 9.4, в частности, примечание 2 к 9.4.2.3).

7.8.2. Точность вычисления числовых значений

7.8.2.1. Если обозначить через n_dig число существенных значащих цифр в числовом представлении величины Z, то предел погрешности вычисления значения Z определяют следующим образом:

a) представляют значение Z в виде c x 10^l, где c — целое число, состоящее из n_dig значащих цифр, l — целое число;

b) определяют по формуле

7.8.2.2. При проверке получения установившейся оценки y_j величины Y_j и стандартной неопределенности u(y_j), соответствующей y_i, при реализации адаптивной процедуры метода Монте-Карло, описанной в 7.8.3, предел погрешности вычисления , j = 1, …, m, определяют по числу существенных значащих цифр в числовом представлении u(y_j).

7.8.2.3. При проверке получения установившейся оценки матрицы корреляционных коэффициентов r_ij, соответствующей оценке , предел погрешности вычисления определяют по числу существенных значащих цифр в числовом представлении — наибольшем собственном значении матрицы (см. примечание 3 к 3.21).

7.8.2.4. Матрица занимает ключевое место в процедурах последовательного оценивания неопределенности. Обычно эта процедура связана с получением значения скалярной величины Q, представляющей собой некоторую комбинацию выходных величин Y_i, т.е.

Используя формулу

(см. примечание 3 к 3.21), стандартную неопределенность u(q), соответствующую оценке

величины Q, можно представить в виде

где

7.8.2.5. Точность вычисления u(q) зависит от точности вычисления и , а последняя зависит от , и, следовательно, от .

7.8.2.6. При оценивании неопределенности в целях последующих расчетов, которые включают в себя более сложные преобразования, например, вычисления по методу наименьших квадратов, связанные с явным или неявным обращением матрицы , необходимо использовать другие критерии останова адаптивной процедуры вычислений. Так в случае применения метода наименьших квадратов критерий может быть основан на расчете (см. 7.8.2.3) по числу существенных значащих цифр в числовом представлении наименьшего собственного значения матрицы (см. примечание 5 к 3.21). При этом устойчивость численного алгоритма решения задачи с использованием метода наименьших квадратов зависит от числа обусловленности матрицы . Поэтому для уменьшения числа испытаний в методе Монте-Карло при постановке измерительной задачи ее следует определять через параметры, позволяющие сделать число обусловленности настолько малым, насколько это возможно.

7.8.2.7. Точность вычисления , используемая для проверки получения установившегося значения коэффициента охвата k_p, определяют по числу существенных значащих цифр в числовом представлении k_p.

7.8.2.8. При последующих расчетах неопределенности, основанных на использовании (см. 7.5) в качестве приближения функции распределения , необходимо убедиться, что такое дискретное представление пригодно для этих расчетов, особенно если их целью является получение области охвата для некоторой величины. Более подробное рассмотрение данного вопроса выходит за рамки настоящего стандарта.

7.8.3. Описание адаптивной процедуры

Практическая реализация адаптивной процедуры метода Монте-Карло с последовательным увеличением числа испытаний состоит в следующем:

a) задают в качестве n_dig небольшое положительное целое число (см. 7.8.2);

b) задают M = max(J, 10⁴), где J — наименьшее целое, большее или равное 100/(1 — p);

c) задают h = 1 (счетчик итераций метода Монте-Карло);

d) выполняют M испытаний методом Монте-Карло (см. 7.3 и 7.4);

e) используют M полученных на выходе модели значений векторной выходной величины для вычислений очередных, h-х значений , , и , являющихся соответственно оценкой величины , соответствующей стандартной неопределенностью, соответствующей корреляционной матрицей и коэффициентом охвата для 100p%-ной области охвата;

f) если h <= 10, то увеличивают h на единицу и возвращаются к этапу d);

g) для j = 1, …, m вычисляют выборочные стандартные отклонения средних значений полученных в результате итераций оценок величины Y_j, по формуле

h) аналогичным образом вычисляют выборочные стандартные отклонения средних значений элементов вектора , а также для и ;

i) используют все hM значений векторной выходной величины для вычисления , и k_p;

j) для j = 1, …, m определяют пределы погрешности вычисления для u(y_j), как описано в 7.8.2.1 и 7.8.2.2;

k) определяют предел погрешности для матрицы коэффициентов корреляции, как описано в 7.8.2.1 и 7.8.2.3;

l) вычисляют численную точность , соответствующую k_p, как описано в пунктах 7.8.2.1 и 7.8.2.7;

m) если для любого j = 1, …, m или [или ] превосходит , или превосходит , или превосходит , то увеличивают h на единицу и возвращаются к этапу d);

n) если возврата к этапу d) не произошло, то считают все вычисленные оценки установившимися и используют все hM значений векторной выходной величины для вычисления , и коэффициент охвата k_p для 100p%-ной области охвата, как описано в 7.6 и 7.7.

Примечание 1. Чтобы алгоритм адаптивной процедуры был сходящимся, должны существовать математическое ожидание и ковариационная матрица .

Примечание 2. Выбор M на этапе b) является произвольным, но практика показала его пригодность.

Примечание 3. Использование в процедуре минимального числа итераций, равного 10, позволяет предотвратить преждевременный останов алгоритма и, кроме того, делает более обоснованным предположение, сделанное в примечании 6. Аналогичное изменение, будучи внесенным в адаптивную процедуру метода Монте-Карло для моделей с единственной скалярной выходной величиной [JCGM 101 (пункт 7.9.4)], также позволит улучшить результаты применения этой процедуры для ряда задач.

Примечание 4. На этапе g) y_j рассматривается как реализация случайной переменной со стандартным отклонением .

Примечание 5. Стандартные отклонения, полученные на этапах g) и h), уменьшаются по закону h^-1/2 (см. примечание к 7.1.7).

Примечание 6. Коэффициент 2, используемый на этапе m), основан на представлении выборочных средних случайными, нормально распределенными величинами и соответствует вероятности охвата, приблизительно равной 95%.

Примечание 7. В [28] рассматриваются некоторые улучшения адаптивной процедуры, изложенной в JCGM 101 (подраздел 7.9).

8. Проверка результатов оценивания неопределенности

по GUM сравнением с методом Монте-Карло

8.1. Способ оценивания неопределенности по GUM (далее обозначаемый GUF) во многих случаях позволяет получить удовлетворительные результаты. Однако не всегда можно сразу определить, соблюдены ли все условия для его применения [см. JCGM 101, пункты 5.7 и 5.8]. Обычно гораздо проще оценить неопределенность с использованием метода Монте-Карло (при наличии соответствующего программного обеспечения), чем выяснить, выполнены ли все условия оценивания по GUM [8]. При наличии сомнений в обоснованности применения способа оценивания по GUM полученные с его помощью результаты нуждаются в проверке, а поскольку диапазон условий, при которых может быть применен метод Монте-Карло, значительно шире, чем для метода по GUM, то для такой проверки рекомендуется сопоставить результаты оценивания по GUM с результатами оценивания методом Монте-Карло. Если сравнение подтвердит обоснованность применения GUM, то способ оценивания неопределенности по GUM можно будет применять в будущем для схожих задач. В противном случае следует рассмотреть возможность замены на другой способ оценивания неопределенности, включая тот же метод Монте-Карло.

8.2. Для сравнения двух методов необходимо вначале:

a) применить способ оценивания неопределенности по GUM для получения (i) оценки величины , (ii) стандартной неопределенности , соответствующей , (iii) корреляционной матрицы , соответствующей , и (iv) коэффициент охвата , определяющего 100p%-ную область охвата для в виде m-мерного эллипсоида;

b) применить адаптивную процедуру метода Монте-Карло (см. 7.8.3), чтобы получить аналогичные оценки , , и .

8.3. Задача процедуры сравнения состоит в том, чтобы определить, согласуются ли между собой результаты, полученные способом оценивания неопределенности по GUM и методом Монте-Карло, в рамках заданной точности вычислений. Для этого:

a) задают в качестве n_dig небольшое положительное целое число (см. 7.8.2);

b) для j = 1, …, m рассчитывают пределы погрешности вычисления для u(y_j), как указано в 7.8.2.1 и 7.8.2.2;

c) рассчитывают предел погрешности вычисления для матрицы коэффициентов корреляции, как указано в 7.8.2.1 и 7.8.2.3;

d) рассчитывают предел погрешности вычисления для k_p, как указано в 7.8.2.1 и 7.8.2.7;

e) сравнивают оценки, соответствующие стандартные неопределенности, коэффициенты корреляции, а также коэффициенты охвата, полученные с использованием способа оценивания неопределенности по GUM и метода Монте-Карло, чтобы определить, обеспечивает ли первый из указанных способов требуемое число правильных цифр в числовой записи полученных результатов. Для этого определяют

j = 1, …, m,

j = 1, …, m,

т.е. абсолютные разности соответствующих численных результатов. Тогда, если для всех j = 1, …, m и не больше чем , не больше чем , а не больше чем , то результат сравнения считают положительным, а способ оценивания неопределенности по GUM — успешно прошедшим проверку.

Примечание 1. Выбор вероятности охвата и формы области охвата влияет на результат сравнения. Поэтому сравнение двух способов оценивания выполняют только для заданных вероятности охвата p и формы области охвата.

Примечание 2. В тех ситуациях, когда построение области охвата не требуется, проверку проводят только на основании полученных значений , и . Если форма области охвата должна быть отлична от эллипсоида, то проверку проводят для соответствующей величины, определяющей объем области охвата заданной формы. Например, если область охвата должна иметь вид m-мерного параллелепипеда, то сравнивают полученные двумя способами значения коэффициента охвата k_q с учетом требуемой точности вычисления .

Примечание 3. При проверке применимости способа оценивания неопределенности по GUM сравнением его результатов с результатами метода Монте-Карло последние должны быть получены при достаточном числе испытаний M. Если для сравнения применяют адаптивную процедуру метода Монте-Карло, то достаточно задать пределы погрешности вычислений в адаптивной процедуре в пять раз меньшими соответствующих пределов погрешности в процедуре проверки [см. JCGM 101 (пункт 8.2)] или, другой вариант, задать число значащих цифр при числовом представлении величин, вычисляемых с помощью адаптивной процедуры, на единицу большим, чем используется в процедуре проверки.

9. Примеры

9.1. Иллюстрации положений настоящего стандарта

9.1.1. В первом примере (см. 9.2) рассматривается линейная модель измерения, в которой входные величины могут быть общими для всех выходных величин или влиять только на некоторые из них. Для частных случаев данного примера существуют аналитические решения.

9.1.2. Во втором примере (см. 9.3) рассматривается нелинейная модель преобразования декартовых координат (действительной и мнимой части комплексной величины) в полярные координаты (модуль и аргумент комплексной величины). Для этого примера также в ряде случаев имеются аналитические решения [6].

9.1.3. В третьем примере (см. 9.3) рассматривается более сложная нелинейная модель. Он аналогичен примеру из GUM, связанному с одновременным измерением активного и реактивного сопротивлений [JCGM 100 (раздел H.2)]. Пример иллюстрирует обработку ряда одновременных независимых наблюдений векторной величины.

9.1.4. Четвертый пример (см. 9.5) посвящен измерению температуры с использованием термометра сопротивления. Этот пример демонстрирует обработку данных для одномерной и многомерной моделей измерения.

9.1.5. Многие из рисунков, используемых в примерах, для их лучшего восприятия даны в цветном исполнении. На контурных графиках каждому уровню контурной линии соответствует свой цвет. Если рисунок состоит из двух и более графиков, то для каждого из таких графиков один и тот же цвет использован для одних и тех же значений уровня, за исключением особо оговоренных случаев. Если для сравнения результатов используется два и более рисунка, как это имеет место при сопоставлении результатов, полученных способом оценивания по GUM и методом Монте-Карло, то соответствующие графики на этих рисунках изображены в одних и тех же границах осей, за исключением случаев, когда между этими результатами имеется существенное различие.

9.1.6. Поскольку первичными выходными данными для метода Монте-Карло являются M выборочных векторов для выходной величины , собранные в матрицу размерности m x M (см. 7.1.6), зачастую желательно представить эти данные в виде приближения соответствующей функции плотности распределения и изобразить эту функцию в виде контурного графика. Рисунки настоящего раздела показывают контурные графики для случая двумерной выходной величины, m = 2. С ростом числа испытаний M контурные линии выборочного распределения должны все больше приближаться к линиям распределения для , что требует соответствующего сглаживания [22, 24]. Некоторые контурные графики, приведенные в настоящем разделе, построены исходя непосредственно из приближения соответствующей плотности распределения. Для других применен соответствующий алгоритм сглаживания контуров. На одном из рисунков (рисунок 10) для демонстрации эффекта сглаживания показаны сглаженные и несглаженные контурные линии.

9.2. Аддитивная модель

9.2.1. Постановка задачи

В этом примере рассматривается аддитивная (линейная) двумерная модель измерения (см. пример в 7.7.2)

Y₁ = X₁ + X₃, Y₂ = X₂ + X₃ (24)

для трех разных примеров сочетаний плотностей распределения для входных величин X_i, рассматриваемых как независимые. Из трех входных величин X₁, X₂ и X₃ величина X₃ описывает фактор, влияющий на обе выходные величины Y₁ и Y₂, тогда как каждая из величин X₁ и X₂ описывает факторы, влияющие только на одну из выходных величин — Y₁ и Y₂ соответственно. В первом примере (см. 9.2.2) все являются плотностями нормального распределения с нулевым математическим ожиданием и единичным стандартным отклонением. Второй пример (см. 9.2.3) идентичен первому, за исключением того, что является плотностью равномерного распределения X₃ также с нулевым математическим ожиданием и единичным стандартным отклонением. Третий пример (см. 9.2.4) идентичен второму, за исключением того, что стандартное отклонение X₃ равно трем, что демонстрирует доминирующее влияние фактора, соответствующего данной входной величине, на результат измерения.

9.2.2. Вычисления и результаты (пример 1)

9.2.2.1. В данном примере каждая входная величина X_i описывается стандартным нормальным распределением, т.е. оценки X_i имеют вид x_i = 0, i = 1, 2, 3, с соответствующими стандартными неопределенностями u(x_i) = 1. Результаты, полученные с применением способа оценивания неопределенности по GUM (см. раздел 6) и методом Монте-Карло (см. раздел 7), показаны в таблице 3 и на рисунках 8 — 10. Некоторые данные в таблице с целью облегчения их сравнения представлены в виде чисел с четырьмя значащими цифрами.

Таблица 3

Результаты измерения способом оценивания неопределенности

по GUM (GUF) и методом Монте-Карло (ММК) для аддитивной

модели [формула (24)], с входными величинами X_i,

описываемыми стандартным нормальным распределением (9.2.2)

9.2.2.2. Способ оценивания неопределенности по GUM, обобщенный на случай нескольких выходных величин, дает оценку величины . Ковариационная матрица для данной оценки, имеющая вид

получена по формуле (3), , где в соответствии с условиями настоящего примера

Коэффициент корреляции, соответствующий оценкам y₁ и y₂ (см. 3.21), будет r(y₁, y₂) = 0,5. 95%-ные области охвата для в форме эллипса и прямоугольника определяются, соответственно, коэффициентами охвата k_p = 2,45 (таблица 1) и k_q = 2,24 (таблица 2). Эти результаты приведены в строке для GUF таблицы 3. В соответствии со способом оценивания неопределенности по GUM выходная величина описывается двумерным нормальным распределением .

9.2.2.3. Применение метода Монте-Карло с числом испытаний M, равным 10⁵, 10⁶ и 10⁷, дает результаты, приведенные в следующих строках таблицы для ММК. В строках для адаптивного ММК указаны значения M, равные 0,35 x 10⁶ и 0,45 x 10⁶ и представляющие собой число испытаний при двух применениях адаптивной процедуры метода Монте-Карло (см. 7.8.3) для n_dig, равного трем (см. примечание 3 в 8.3). Все числовые значения, полученные с применением адаптивной процедуры, в частности, для y₁, y₂, u(y₁), u(y₂), r(y₁, y₂), k_p и k_q, являются установившимися.

9.2.2.4. Полученная аналитически плотность распределения для представляет собой двумерное нормальное распределение с и , приведенными в 9.2.2.2.

9.2.2.5. На рисунке 8 слева показана плотность двумерного нормального распределения для , полученная способом оценивания неопределенности по GUM (совпадающая с аналитическим решением). Контурные линии равной вероятности имеют форму эллипсов и определяются уравнениями

для различных значений k. На рисунке 8 справа показаны контурные линии приближения плотности распределения для , полученные с применением метода Монте-Карло с M = 10⁷ испытаниями. Это приближение представляет собой отношение числа точек выборки из распределения , попавших в ячейку сетки разбиения плоскости , к общему числу испытаний M, приписанное области плоскости , занимаемой этой ячейкой (внутри данной ячейки плотность распределения считают постоянной). Контурные линии приближенного распределения построены для тех же значений k, что были использованы для построения эллиптических контурных линий на рисунке 8 слева. Алгоритм сглаживания контурных линий при этом не применялся (см. 9.1.6).

Примечание. На рисунке 8 и других контурных графиках приведена цветовая шкала, показывающая соответствие цвета контурной линии определенному уровню (вероятности).

Рисунок 8. Контурные линии совместных плотностей

распределения выходных величин в аддитивной модели

измерения [формула (24)], полученных способом оценивания

неопределенности по GUM (слева) и методом Монте-Карло

(справа) без сглаживания контуров при стандартном

нормальном распределении входных величин (9.2.2)

9.2.2.6. На рисунке 9 показаны маргинальная плотность распределения для Y₁, полученная способом оценивания неопределенности по GUM, а также приближение плотности распределения той же величины, полученное методом Монте-Карло с числом испытаний M = 10⁷. Приближение представлено в виде распределения частот (гистограммы). Эти две плотности распределения визуально практически не различимы. Аналогичный результат будет получен и для Y₂.

Рисунок 9. Маргинальная плотность распределения Y₁

для тех же условий, что и на рисунке 8 (9.2.2)

9.2.2.7. На рисунке 10 слева показаны контурные линии приближения плотности распределения для , полученные методом Монте-Карло с 0,45 x 10⁶ испытаниями. Эти линии значительно менее «гладкие», чем те, что изображены на рисунке 8 (справа) и получены для гораздо большего числа испытаний. На рисунке 10 справа приведен пример сглаженных контурных линий (см. 9.1.6).

Примечание. В большинстве случаев гладкость контурных линий для приближения плотности может быть обеспечена выбором большего числа испытаний в методе Монте-Карло и более частой сеткой разбиения плоскости . Аналогичный эффект имеет место в случае единственной скалярной выходной величины, когда для обеспечения гладкости огибающей кривой построенного приближения плотности распределения (гистограммы) необходимо большое значение M и большое число классов гистограммы.

Рисунок 10. Контурные линии совместных плотностей

распределения выходных величин в аддитивной модели

измерения [формула (24)], полученных с использованием

адаптивной процедуры Монте-Карло без сглаживания (слева)

и со сглаживанием (справа) для тех же условий,

что и на рисунке 8 (9.2.2)

9.2.2.8. Для сопоставления результатов, полученных способом оценивания неопределенности по GUM и с помощью адаптивной процедуры метода Монте-Карло, применена процедура проверки, описанная в разделе 8, при значении n_dig, равном двум. В этом случае точности вычисления будут определяться значениями

Адаптивная процедура метода Монте-Карло была применена дважды. В обоих случаях продемонстрирована хорошая согласованность результатов, чего и следовало ожидать, поскольку в данном примере полностью соблюдены все условия применимости способа оценивания неопределенности по GUM. Исключением мог стать только результат вычисления k_q, поскольку расчет значения этой величины способом оценивания неопределенности по GUM не учитывает взаимную зависимость между Y₁ и Y₂ (см. 6.5.2.3). Однако и для k_q результаты, полученные двумя методами, показали хорошую согласованность в пределах заданной точности вычисления.

9.2.3. Вычисления и результаты (пример 2)

9.2.3.1. Этот пример измерительной задачи идентичен описанному в 9.2.2, за тем исключением, что X₃ характеризуется равномерным распределением с нулевым математическим ожиданием и единичным стандартным отклонением. Оценки X_i и соответствующие стандартные неопределенности остались неизменными: x_i = 0, u(x_i) = 1, i = 1, 2, 3. Результаты, полученные с применением способа оценивания неопределенности по GUM (см. раздел 6) и методом Монте-Карло (см. раздел 7), показаны в таблице 4 и на рисунках 11 и 12.

Таблица 4

Результаты измерения для тех же условий, что в таблице 3,

за исключением того, что X₃ описывается равномерным

распределением (9.2.3)

Рисунок 11. То же, что на рисунке 8 для примера с входной

величиной X₃, описываемой равномерным

распределением (9.2.3)

Рисунок 12. То же, что на рисунке 9 для примера с входной

величиной X₃, описываемой равномерным

распределением (9.2.3)

9.2.3.2. Способ оценивания неопределенности по GUM дает абсолютно одинаковую двумерную плотность распределения для (см. рисунок 11, слева) независимо от того, распределена ли входная величина X₃ по нормальному (как в 9.2.2) или равномерному (как в настоящем примере) закону, поскольку в данной измерительной задаче вид этой плотности распределения полностью определен оценками входных величин и их стандартными отклонениями, а они в обоих примерах одинаковы. Поскольку модель измерения изначально линейна, ее дополнительной линеаризации для применимости закона трансформирования неопределенностей не требуется, и в этом смысле способ оценивания неопределенности по GUM никаких приближений не использует. Тем не менее, результаты, полученные способом оценивания неопределенности по GUM для оценки величины и соответствующей ковариационной матрицы , не будут совпадать с аналитическим решением, поскольку последнее зависит не только от оценок входных величин и их стандартных отклонений, но и от вида их распределений, в данном примере от распределения X₃.

9.2.3.3. На рисунках 11 (справа) и 12 показано, каким образом негауссовость распределения X₃ влияет на приближения плотности распределения и маргинальной плотности распределения для Y₁, полученные методом Монте-Карло.

9.2.3.4. Проверка применимости способа оценивания неопределенности по GUM, описанная в разделе 8 (см. также 9.2.2), была проведена сравнением результатов, полученных с его помощью и с помощью адаптивной процедуры Монте-Карло. Из таблицы 4 видно хорошее совпадение результатов для y₁, y₂, u(y₁), u(y₂) и r(y₁, y₂), но недостаточно хорошее для k_p и k_q.

9.2.4. Вычисления и результаты (пример 3)

9.2.4.1. Этот пример измерительной задачи идентичен описанному в 9.2.3, за исключением того, что стандартное отклонение для X₃ теперь равно не единице, а трем. Результаты, полученные с применением способа оценивания неопределенности по GUM (см. раздел 6) и методом Монте-Карло (см. раздел 7), показаны в таблице 5 и на рисунках 13 и 14.

Таблица 5

Результаты измерения для тех же условий, что в таблице 4,

за исключением того, что стандартное отклонение для X₃

равно трем (9.2.3)

Рисунок 13. То же, что на рисунке 11 для примера с входной

величиной X₃, имеющей стандартное отклонение,

равное трем (9.2.4)

Рисунок 14. То же, что на рисунке 12 для примера с входной

величиной X₃, имеющей стандартное отклонение,

равное трем (9.2.4)

9.2.4.2. Способ оценивания неопределенности по GUM дает оценку величины . Ковариационная матрица для данной оценки, имеющая вид

получена по формуле (3), где в соответствии с условиями настоящего примера

Коэффициент корреляции для оценок y₁ и y₂ будет r(y₁, y₂) = 0,9. 95%-ные области охвата для в форме эллипса и прямоугольника определяются, соответственно, коэффициентами охвата k_p = 2,45 (таблица 1) и k_q = 2,24 (таблица 2). Эти результаты приведены в строке для GUF таблицы 5. В соответствии со способом оценивания неопределенности по GUM выходная величина описывается двумерным нормальным распределением , показанным в виде контурного графика на рисунке 13 (слева). В этом примере наблюдается более сильная корреляция между Y₁ и Y₂ по сравнению с примерами, рассмотренными ранее, поскольку увеличен относительный вклад общего фактора, описываемого через X₃, в неопределенности, связанные с оценками Y₁ и Y₂.

9.2.4.3. Из рисунков 13 (справа) и 14 видно большее влияние (в сравнении с результатами в 9.2.3) распределения, описывающего X₃, на приближения плотности распределения для и маргинальной плотности распределения для Y₁, полученные методом Монте-Карло.

9.2.4.4. Проверка применимости способа оценивания неопределенности по GUM, описанная в разделе 8 (см. также 9.2.2 и 9.2.3), была проведена сравнением результатов, полученных с его помощью и с помощью адаптивной процедуры Монте-Карло. Из таблицы 5 видно хорошее совпадение результатов для y₁, y₂, u(y₁), u(y₂) и r(y₁, y₂), но плохое для k_p и k_q. Полученное способом оценивания неопределенности по GUM значение k_p приблизительно на 7%, а значение k_q приблизительно на 20% больше тех, что получены методом Монте-Карло.

9.3. Преобразование системы координат

9.3.1. Постановка задачи

9.3.1.1. Комплексная величина может быть представлена в алгебраической форме

X₁ + iX₂,

где X₁ = Z_R и X₁ = Z_I — соответственно действительная и мнимая части Z, или в тригонометрической форме

где R и — соответственно модуль и аргумент . Эти две формы записи соответствуют представлениям точки в комплексной плоскости соответственно с декартовой и полярной системой координат. Преобразованию первой системы во вторую можно поставить в соответствие модель измерения вида

tanY₂ = X₂/X₁ (25)

с входной величиной и выходной величиной .

Примечание. Формула (25) описывает двумерную модель измерения для выходных величин R и . Поскольку модуль R всегда неотрицателен, то он однозначно может быть определен как положительный квадратный корень из R². Если аргумент определять через так называемый арктангенс двух переменных (в языках программирования обозначаемый «atan2»), то при заданном отношении Z_|/Z_R значение аргумента также будет однозначно определено на интервале . Таким образом, формула (25) допускает использование ее в качестве двумерной модели измерения.

9.3.1.2. Исходными данными для расчета неопределенности являются оценки x₁ и x₂ величин X₁ и X₂, полученные из измерительной системы, и соответствующие им стандартные неопределенности u(x₁) и u(x₂) и ковариации u(x₁, x₂) = ru(x₁)u(x₂), где r = r(x₁, x₂) — коэффициент корреляции [см. JCGM 100 (пункт 5.2.2)]. При наличии таких данных [см. JCGM 101 (пункт 6.4.8.1)] приписывают двумерное нормальное распределение с математическим ожиданием и ковариационной матрицей соответственно

Предполагается, что размерность X₁ и X₂ равна единице.

9.3.1.3. Оценки величины и соответствующей ковариационной матрицы определены для разных сочетаний величин x₁, x₂, u(x₁), u(x₂) и r(x₁, x₂).

9.3.1.4. Рассмотрено шесть примеров, в каждом из которых x₂ взято равным нулю, а u(x₁) = u(x₂) = u_x = 0,010. В первых трех примерах рассматривается случай некоррелированных входных величин, r(x₁, x₂) = 0 для оценок соответственно x₁ = 0,001, x₁ = 0,010 и x₁ = 0,100 (см. 9.3.2). В остальных трех примерах использованы те же оценки для X₁, но при сильной корреляции входных величин: r(x₁, x₂) = 0,9 (см. 9.3.3). На рисунке 15 [где не соблюден принцип использования на каждом графике одного цвета для одинаковых уровней (см. 9.1.5)] показаны контурные линии совместной плотности распределениях для примера 1 [x₁ = 0,001 и r(x₁, x₂) = 0] и примера 4 [x₁ = 0,001 и r(x₁, x₂) = 0,9]. Для остальных примеров графики плотности распределения получают переносом контурных линий вдоль оси X₁ таким образом, чтобы их центрам соответствовала координата x₁ = 0,010 (примеры 2 и 5) или x₁ = 0,100 (примеры 3 и 6).

Рисунок 15. Контурные линии совместных плотностей

распределения входных величин в модели преобразования

системы координат для примеров 1 (слева)

и 4 (справа) (9.3.1.4)

9.3.2. Вычисления и результаты (случай нулевой ковариации)

9.3.2.1. Общие положения

9.3.2.1.1. Оценивание неопределенности выполняется применением трансформирования распределений (a) аналитически (для целей сравнения), (b) с использованием обобщенного способа оценивания неопределенности по GUM и (c) с использованием метода Монте-Карло.

9.3.2.1.2. Совместная плотность распределения для и маргинальные плотности распределения для Y₁ и Y₂ могут быть получены аналитически в случае независимых X₁ и X₂, описываемых нормальными распределениями и соответственно (см. раздел C.2).

9.3.2.1.3. Согласно способу оценивания неопределенности по GUM величине соответствует двумерное нормальное распределение , при этом оценку для получают решением уравнений

tany₂ = x₂/x₁,

а ковариационную матрицу получают применением закона трансформирования неопределенностей (см. разделы 6 и C.3).

9.3.2.1.4. Метод Монте-Карло применен с числом испытаний M = 10⁷ (см. раздел 7).

9.3.2.2. Пример 1 (x₁ = 0,001)

9.3.2.2.1. Результаты, полученные аналитически, способом оценивания неопределенности по GUM и методом Монте-Карло для входной оценки x₁ = 0,001 и коэффициента корреляции r(x₁, x₂) = 0, приведены в таблице 6, строки 1 — 3.

Примечание. Приведенные в первой строке таблицы 6 значения y₁ и u(y₁) получены в результате выполненных с высокой точностью численных расчетов определенных интегралов, представляющих собой заданные в виде формул через маргинальную плотность распределения Y₁ выражения для E(Y₁) и V(Y₁) [см. формулу (C.2) в Приложении C]. Аналогично значения y₂ и u(y₂) рассчитаны численно по формулам для E(Y₂) и V(Y₂), включающим в себя маргинальное распределение Y₂ [см. формулу (C.3) в Приложении C)]. Нетрудно показать, что cov(Y₁, Y₂) = 0, и, следовательно, r(y₁, y₂) = 0.

Таблица 6

Результаты измерений для модели преобразования системы

координат для случая нулевой ковариации между оценками

входных величин (9.3.2.2.1, 9.3.2.3.1 и 9.3.2.4.1)

9.3.2.2.2. На трех верхних графиках рисунка 16 показаны плотности распределения для , полученные аналитически, способом оценивания неопределенности по GUM и методом Монте-Карло. Последняя хорошо согласуется с аналитическим решением, в то время как результат оценивания по GUM отличается от них очень сильно, вплоть до того, что для него использована другая шкала соответствия цвета и уровня (см. 9.1.5). Более того, способ оценивания неопределенности по GUM дает ненулевую вероятность появления физически невозможных значений выходных величин в областях , и .

Рисунок 16. Совместные плотности вероятности для ,

полученные аналитически (сверху), оцениванием по GUM

(в центре слева) и методом Монте-Карло (в центре справа)

и маргинальные плотности распределения для Y₁ (внизу слева)

и Y₂ (внизу справа), полученные аналитически (пунктирная

линия), оцениванием по GUM (сплошная линия) и методом

Монте-Карло (гистограмма) для модели преобразования

координат с x₁ = 0,001 и r(x₁, x₂) = 0

(9.3.2.2.2 и 9.3.2.2.3)

9.3.2.2.3. На двух нижних графиках рисунка 16 изображены маргинальные плотности распределения для выходных величин Y₁ и Y₂, полученные из совместных распределений для . Маргинальные плотности распределения, полученные методом Монте-Карло (показаны каждая в форме гистограммы или распределения частот), согласуются с теми, что получены аналитически (показаны пунктирными линиями), в то время как маргинальные распределения, полученные способом оценивания неопределенности по GUM (показаны сплошными линиями), очень сильно отличаются от аналитических решений. Для удобства представления все графики маргинальной плотности распределения для Y₂ ограничены областью физически возможных значений , т.е. интервалом , хотя полученное способом оценивания по GUM нормальное распределение N(0, 10²) для Y₂ простирается далеко за границы этого интервала. На этих графиках также хорошо видно, что согласно способу оценивания неопределенности по GUM физически невозможным значениям и присвоена ненулевая вероятность (см. 9.3.2.2.2).

9.3.2.3. Пример 2 (x₁ = 0,010)

9.3.2.3.1. Результаты, полученные аналитически, способом оценивания неопределенности по GUM и методом Монте-Карло для входной оценки x₁ = 0,010 и коэффициента корреляции r(x₁, x₂) = 0, приведены в строках 4 — 6 таблицы 6. На рисунке 17 показаны контурные графики совместных плотностей распределения для , полученных с использованием трех вышеуказанных методов, а также маргинальные плотности вероятности для Y₁ и Y₂, полученные из соответствующих совместных распределений.

Рисунок 17. То же, что на рисунке 16, но для x₁ = 0,010

(9.3.2.3.1)

9.3.2.3.2. Видно, что результаты, полученные с использованием метода Монте-Карло, согласуются с аналитическим решением. В то время как результаты, полученные с использованием способа оценивания неопределенности по GUM, от аналитического решения отличаются, хотя и не так существенно, как в примере с x₁ = 0,001. Относительные разности между стандартными неопределенностями u(y₁) и u(y₂), определенные с использованием способа оценивания неопределенности по GUM и аналитическим методом, составляют приблизительно 25% и 10% (в сравнении с 40% и 470% для примера с x₁ = 0,001).

9.3.2.4. Пример 3 (x₁ = 0,100)

9.3.2.4.1. Результаты, полученные аналитически, оцениванием неопределенности по GUM и методом Монте-Карло для входной оценки x₁ = 0,100 и коэффициента корреляции r(x₁, x₂) = 0, приведены в строках 7 — 9 таблицы 6 и на рисунке 18.

Рисунок 18. То же, что на рисунке 16, но для x₁ = 0,100

(9.3.2.4.1)

9.3.2.4.2. Видно, что результаты, полученные с использованием способа оценивания неопределенности по GUM и методом Монте-Карло, согласуются с аналитическим решением. Показанные на рисунке 18 маргинальные распределения, полученные тремя методами, практически неразличимы. Оценки, соответствующие стандартные неопределенности и соответствующие коэффициенты корреляции совпадают с точностью до двух значащих цифр.

9.3.3. Вычисления и результаты (случай ненулевой ковариации)

9.3.3.1. Оценивание неопределенности выполняется применением трансформирования распределений с использованием (a) обобщенного способа оценивания неопределенности по GUM (см. разделы 6 и C.3) и (b) метода Монте-Карло с числом испытаний M = 10⁷ (см. раздел 7).

9.3.3.2. Результаты, полученные для входных оценок x₁ = 0,001, x₁ = 0,010 и x₁ = 0,100 (примеры 4, 5 и 6 соответственно) и коэффициента корреляции r(x₁, x₂) = 0,9, приведены в таблице 7. На рисунках 19 — 21 показаны совместные плотности распределения для и маргинальные плотности распределения для Y₁ и Y₂, полученные указанными двумя методами для всех трех примеров.

Таблица 7

Результаты измерений для модели преобразования системы

координат для случая ненулевой ковариации между оценками

входных величин (9.3.3.2)

Рисунок 19. Совместные плотности вероятности для ,

полученные оцениванием по GUM (вверху слева) и методом

Монте-Карло (вверху справа), и маргинальные плотности

распределения для Y₁ (внизу слева) и Y₂ (внизу справа),

полученные оцениванием по GUM (сплошная линия) и методом

Монте-Карло (гистограмма) для модели преобразования

координат с x₁ = 0,001 и r(x₁, x₂) = 0,9 (9.3.3.2)

Рисунок 20. То же, что на рисунке 19,

но для x₁ = 0,010 (9.3.3.2)

Рисунок 21. То же, что на рисунке 19,

но для x₁ = 0,100 (9.3.3.2)

9.3.3.3. Для случаев x₁ = 0,001 и x₁ = 0,010 результаты, полученные способом оценивания неопределенности по GUM и методом Монте-Карло, плохо согласуются друг с другом. В частности, метод Монте-Карло дает в этих двух примерах маргинальные плотности распределения для Y₂ с двумя пиками, что сильно отличается от одномодального нормального распределения, полученного оцениванием по GUM. Эти пики находятся вблизи точек и , что соответствует углам ориентации главной оси эллиптических контурных линий плотности распределения для [см. рисунок 15 (справа)].

9.3.3.4. Для случая x₁ = 0,100 результаты, полученные способом оценивания неопределенности по GUM и методом Монте-Карло, согласуются намного лучше.

9.3.4. Обсуждение результатов

9.3.4.1. Для обоих случаев (с нулевой и ненулевой ковариацией) по мере удаления оценки x₁ от нуля результаты, полученные способом оценивания неопределенности по GUM и методом Монте-Карло, начинают приближаться друг к другу.

9.3.4.2. Для оценок x₁ = 0,001, x₁ = 0,010 и, вообще говоря, для всех значений x₁, близких к нулю, эффект ненулевой ковариации значительно изменяет результаты, полученные с помощью метода Монте-Карло.

9.3.4.3. Численные данные, представленные в таблицах 6 и 7, соответствуют показанным на рисунках совместным и маргинальным плотностям распределения. В некоторых случаях такие данные могут быть не подходящими или недостаточными для описания распределения, характеризующего выходную величину. Так, для примера с входной оценкой x₁ = 0,001 и ненулевой ковариацией r(x₁, x₂) = 0,9 (рисунок 19) маргинальная плотность распределения для Y₂, полученная с помощью метода Монте-Карло, существенно бимодальна, между тем оценка y₂ величины Y₂ находится между модами распределения в той области значений Y₂, где плотность вероятности мала.

Примечание. Для выходной величины, имеющей многомерное нормальное распределение, вектор математического ожидания и соответствующая ковариационная матрица описывают это распределение исчерпывающим образом.

9.3.4.4. Для входных оценок x₁ = 0,001, x₁ = 0,010 и, вообще говоря, для всех значений x₁, близких к нулю, определение интервалов охвата в форме эллипсов и прямоугольников не является подходящим.

9.4. Одновременное измерение активного и реактивного сопротивлений

9.4.1. Постановка задачи

9.4.1.1. Активное R и реактивное X сопротивления элемента цепи определяют путем измерения амплитуды V изменяющегося по гармоническому закону напряжения на его клеммах, амплитуды I проходящего через элемент переменного тока и фазового сдвига между напряжением и силой тока. Входными величинами модели являются V, I и , а выходными — указанные сопротивления R и X, а также модуль полного импеданса элемента цепи Z, Z² = R² + X².

9.4.1.2. Применение закона Ома дает трехмерную модель измерения

(26)

связывающую входную величину с выходной величиной .

Примечание 1. В настоящем примере в целях упрощения не учитываются систематические эффекты, которые могут оказывать влияние на оценки V, I и .

Примечание 2. Аналогичный пример рассмотрен в JCGM 100 (раздел H.2), где реактивному сопротивлению присвоено обозначение X. Такое же обозначение используется в настоящем примере. Реактивное сопротивление X является элементом векторной выходной величины , и его не нужно путать с — векторной входной величиной.

9.4.1.3. В одинаковых условиях проведено n = 6 повторных независимых наблюдений , …, входной величины . Результаты этих наблюдений приведены в таблице 8.

Примечание. В примере из JCGM 100 (раздел H.2) число наблюдений было равно пяти, и результаты тех наблюдений приведены в первых пяти строках таблицы 8. Однако для определения ковариационной матрицы (см. 9.4.2.5) необходимо как минимум 6 наблюдений. Поэтому в таблицу 8 было добавлено шестое наблюдение, полученное как среднее арифметическое значение первых пяти наблюдений. Для настоящего примера несущественно, каким способом было получено шестое наблюдение, но выбор в качестве шестого наблюдения среднего значения выборки позволяет сохранить это среднее значение неизменным.

Таблица 8

Данные n = 6 одновременных независимых наблюдений входных

величин для модели измерения активного и реактивного

сопротивлений (9.4.1.3)

9.4.2. Вычисления и результаты

9.4.2.1. Результат измерения в настоящем примере должен быть представлен в виде оценки величины и соответствующей ковариационной матрицы . Этот результат должен быть получен на основе модели измерения, заданной формулами (26), и данных наблюдений , …, , приведенных в таблице 8. Решение данной задачи на основе способа оценивания неопределенности по GUM приведено в JCGM 100 (пункт H.2.3). В настоящем примере помимо данного способа рассматриваются модификация способа оценивания неопределенности по GUM на основе применения многомерного t-распределения (см. 5.3.2) для входных величин, а также применение метода Монте-Карло.

9.4.2.2. Для того, чтобы применить способ оценивания неопределенности по GUM, необходимо знать оценку входной величины , определяемую как выборочное среднее по наблюдениям, данные которых приведены в таблице 8 [JCGM 100 (пункт 4.2)]:

Ковариационная матрица , соответствующая , содержит дисперсии, соответствующие средним значениям [JCGM 100 (пункт 4.2)], и ковариации, соответствующие каждой паре средних значений [JCGM 100 (пункт 5.2.3)], и вычисляется по формуле:

где — матрица сумм квадратов и произведений. Оценки входных величин и соответствующие стандартные неопределенности приведены в таблице 9, а коэффициенты корреляции, соответствующие парам таких оценок, — в таблице 10.

Таблица 9

Оценки входных величин при одновременном измерении

активного и реактивного сопротивлений и соответствующие

стандартные неопределенности (9.4.2.2)

Таблица 10

Коэффициенты корреляции, соответствующие парам оценок

входных величин , при одновременном измерении

активного и реактивного сопротивлений (9.4.2.2)

Примечание. В JCGM 100 (раздел H.2) для данной задачи рассмотрены два способа оценивания неопределенности измерения, причем принцип, положенный в основу второго способа, изложен в JCGM 100 в примечании к пункту 4.1.4. В настоящем стандарте этот второй способ не рассматривается по причинам, указанным в 4.1.

9.4.2.3. Согласно способу оценивания неопределенности по GUM, оценку выходной величины рассчитывают на основе оценки по формулам (26)

Ковариационную матрицу , соответствующую , вычисляют по формуле (3) из 6.2.1.3, , где — матрица чувствительностей, имеющая вид

(27)

Результаты, полученные с применением способа оценивания неопределенности по GUM, приведены в первой строке (метод GUF) таблицы 11.

Примечание 1. В последнем столбце таблицы 11 приведены значения для 1 — r(X, Z), поскольку коэффициент корреляции r(X, Z) близок к единице (см. 3.21, примечание 5).

Примечание 2. При имеющейся в модели, описываемой формулами (26), зависимости между выходными величинами

Z² = R² + X², (28)

ковариационная матрица теоретически должна быть сингулярна. Вследствие ошибок округления корреляционная матрица, элементы которой приведены в первой строке таблицы 11, не является ни сингулярной, ни положительно определенной. Но из-за малости стандартных неопределенностей для оценок входных величин эта сингулярность не оказывает практически значимого влияния на полученные результаты [20, раздел 4].

Таблица 11

Результаты одновременного измерения активного и реактивного

сопротивлений (9.4.2.3, 9.4.2.4 и 9.4.2.5)

9.4.2.4. В предположении, что данные, приведенные в таблице 8, являются единственной доступной информацией о входных величинах и что каждое наблюдение можно рассматривать как выборку из одного и того же многомерного нормального распределения, входную величину описывают многомерным t-распределением t_v ), с степенями свободы, где — выборочное среднее, определенное в 9.4.2.2 (см. 5.3.2). Результаты, полученные с применением M = 10⁶ испытаний метода Монте-Карло, показаны во второй строке (метод ММК) таблицы 11.

9.4.2.5. Величина , описываемая многомерным t-распределением (см. 9.4.2.4) имеет ковариационную матрицу

Учет этого обстоятельства позволяет несколько улучшить процедуру оценивания неопределенности по GUM, рассмотренную в 9.4.2.3, заменив в формуле (3), , матрицу на . Результаты, полученные с применением такого альтернативного подхода, приведены в третьей строке (метод Альтернативный GUF) таблицы 11 (см. также [15]).

Примечание 1. — ковариационная матрица, полученная с учетом дополнительной информации, рассмотренной в 9.4.2.4.

Примечание 2. Ковариационная матрица для , полученная таким способом, может быть определена только при выполнении условия . Именно по этой причине для настоящего примера потребовалось не менее n = 6 одновременных наблюдений для N = 3 входных величин (см. примечание к 9.4.1.3).

9.4.2.6. Ковариационная матрица , полученная из распределения для , связана с ковариационной матрицей , используемой в способе оценивания неопределенности по GUM, соотношением

9.4.2.7. Расхождение результатов, полученных при применении метода Монте-Карло и альтернативного способа оценивания неопределенности по GUM, пренебрежимо мало. Это означает, что функции измерения в формуле (26) могут быть линеаризованы с хорошей степенью приближения в окрестности оценок входных величин.

9.4.2.8. В JCGM 100 рассматривается возможность описания выходной величины t-распределением с использованием формулы Уэлча-Саттертуэйта [см. JCGM 100, формула (G.2b)] для расчета числа эффективных степеней свободы. Однако применение этой формулы предполагает независимость входных величин, стандартным отклонениям которых соответствуют конечные числа степеней свободы. Данное условие в настоящем примере не выполняется.

9.5. Измерение температуры с использованием термометра сопротивления

9.5.1. Основное

В настоящем примере рассматривается измерение температуры промышленным платиновым термометром сопротивления путем сравнения сопротивления термометра с эталонным сопротивлением в схеме измерительного моста. Если измерению подлежит конкретная температура, то для этой цели используют одномерную модель измерения (см. 9.5.2), а если нескольких температур — то многомерную модель (см. 9.5.3). В примере рассматривается обработка данных в рамках указанных моделей измерения способом оценивания неопределенности по GUM.

9.5.2. Измерение одной температуры

9.5.2.1. Температуру измеряют сравнением сопротивления термометра сопротивления с эталонным сопротивлением R_S измерительного моста. На интервале температур от 0 до 30 °C сопротивление термометра может быть приближенно описано квадратичной функцией его температуры :

(29)

где R₀, A и B — параметры, определяемые при градуировке (калибровке) термометра. Оценки R₀, A и B и соответствующие стандартные неопределенности приведены в таблице 12, а коэффициенты корреляции для пар таких оценок — в таблице 13.

Таблица 12

Оценки входных величин 

и соответствующие стандартные неопределенности

при измерении температуры (9.5.2.1, 9.5.2.2 и 9.5.2.3)

Таблица 13

Коэффициенты корреляции, соответствующие парам оценок

входных величин при измерении температуры

(9.5.2.1, 9.5.2.2 и 9.5.2.3)

9.5.2.2. Оценка величины R_S и соответствующая стандартная неопределенность, получаемые при калибровке, приведены в таблице 12. R_S не зависит от параметров R₀, A и B.

9.5.2.3. Посредством измерительного моста определяют отношение сопротивлений

(30)

Полученное значение r и соответствующая стандартная неопределенность приведены в таблице 12. Отношение сопротивлений не зависит от параметров R₀, A и B термометра сопротивления и эталонного сопротивления R_S. Таким образом, в таблице 13 содержатся все ненулевые коэффициенты корреляции входных величин.

9.5.2.4. Модель измерения температуры , полученная из формул (29) и (30), имеет вид

(31)

В обозначениях, принятых в настоящем стандарте, N = 5, m = 1, , и

Примечание. Модель измерения, определяемая формулой (31), может быть преобразована к явному виду путем решения квадратного уравнения относительно . Однако такое преобразование включает в себя процедуру вычитания близких чисел, что может привести к потере точности вычисления, и, кроме того, оно усложняет вычисление коэффициентов чувствительности.

9.5.2.5. Оценку температуры получают на основе полученного отношения сопротивлений r, подставляя данные из таблицы 12 в уравнение (31) и находя решение этого уравнения. Искомая оценка равна .

9.5.2.6. Стандартную неопределенность , соответствующую оценке y, вычисляют с использованием формулы (8) из 6.3.1.3, . Вычисление матриц чувствительности и по формулам

с подстановкой оценок входных величин, приведенных в таблице 12, и соответствующей оценки выходной величины дает

c_y = 0,389 Ом·°C^-1,

Элементы ковариационной матрицы , соответствующей оценкам входных величин, вычисляют из стандартных неопределенностей в таблице 12 и коэффициентов корреляции в таблице 13, что дает .

9.5.3. Измерение нескольких температур

9.5.3.1. Термометр сопротивления, эталонное сопротивление и измерительный мост, описанные в 9.5.2, используют для измерения отношений сопротивлений r_{1, …}, r₁₀, соответствующих десяти разным температурам , …, .

9.5.3.2. Оценки входных величин R₀, A, B и R_S и соответствующие стандартные неопределенности приведены в таблице 12, а оценки r₁, …, r₁₀ — в таблице 14. Единственными ненулевыми коэффициентами корреляции, соответствующими парам оценок, по-прежнему остаются коэффициенты, приведенные в таблице 13. Предполагается, что отношения сопротивлений независимы между собой. Данное предположение можно считать справедливым, если преобладающий вклад в неопределенность измерения отношения сопротивлений вносят случайные эффекты.

Таблица 14

Оценки отношения сопротивлений и соответствующие

им стандартные неопределенности при измерении нескольких

температур (9.5.3.2)

9.5.3.3. Из формулы (31) следует, что связь отношения сопротивлений r_j с соответствующей температурой задается уравнением

j = 1, …, 10. (32)

В обозначениях, принятых в настоящем стандарте, N = 14, m = 10, , и

Примечание. Модель измерения, описываемая формулой (32), может быть приведена к явному виду (см. 9.5.2.4).

9.5.3.4. Оценки температуры (см. таблицу 15) получают подстановкой оценок, приведенных в столбцах с первого по четвертый таблицы 12, а также в таблице 14, в уравнения (32) и решением этих уравнений.

Таблица 15

Оценки выходных величин и соответствующие стандартные

неопределенности при измерении нескольких температур

(9.5.3.4 и 9.5.3.5)

9.5.3.5. Ковариационную матрицу , соответствующую , вычисляют с использованием формулы (8) из 6.3.1.3, , где и — матрицы чувствительности, получаемые из матриц и в точках оценок входных и выходных величин; — диагональная матрица размерности 10 x 10 с диагональными элементами ; — матрица размерности 10 x 14 вида

где

матрица размерности 10 x 4, а — диагональная матрица размерности 10 x 10 со всеми диагональными элементами, равными (-R_S). Ковариационную матрицу вычисляют из стандартных неопределенностей, приведенных в столбцах 1 — 4 таблицы 12 и в таблице 13. Стандартные неопределенности для оценок температур и коэффициенты корреляции для пар этих оценок, полученные из матрицы , приведены в таблицах 15 и 16, соответственно.

Таблица 16

Коэффициенты корреляции, соответствующие парам оценок

выходных величин , при измерении нескольких температур

(9.5.3.5)

9.5.3.6. Результаты, приведенные в таблице 15 и на рисунке 22, показывают, как стандартная неопределенность изменяется вместе с оценкой температуры . Неопределенность измерения температуры минимальна вблизи нуля и резко возрастает при температурах выше 25 °C. Этот эффект обусловлен тем, что градуировка термометра сопротивления была выполнена при опорных значениях температуры 0, 15, 20 и 25 °C, причем температура 0 °C создавалась при помощи ледяной ванны, в условиях которой стандартная неопределенность была в три раза меньшей, чем в условиях масляной ванны, используемой при градуировке на трех других опорных значениях.

Примечание. Отрезки прямых линий, соединяющих точки на рисунках 22 и 23, использованы в целях большей наглядности.

Рисунок 22. Стандартная неопределенность , соответствующая

оценке температуры (9.5.3.6)

Рисунок 23. Коэффициент корреляции, соответствующий паре

оценок и температур и (9.5.3.7)

9.5.3.7. На рисунке 23 график, построенный по данным последнего столбца таблицы 16, показывает, как изменяется коэффициент корреляции, соответствующий паре оценок и температур по Цельсию и в зависимости от , j = 1, …, 9. Коэффициент корреляции имеет максимум при и стремится к нулю при увеличении абсолютной разности . Пример демонстрирует, что величины, измеряемые одним и тем же средством измерения, могут иметь высокую степень корреляции.

Приложение A

(справочное)

ПРОИЗВОДНЫЕ МНОГОМЕРНЫХ ФУНКЦИЙ ИЗМЕРЕНИЯ

С КОМПЛЕКСНЫМИ ВЕЛИЧИНАМИ

A.1. В настоящем приложении рассматривается компактный алгоритм вычисления частных производных первого порядка функции измерения в многомерной модели измерения с комплексными величинами

где

и

— обозначает комплексную величину X_i,R + iX_i,I, X_i,R и X_i,I — действительные величины, а i² = -1. Аналогичные представления справедливы для комплексных величин и .

A.2. Пусть обозначает ковариационную матрицу размерности 2N x 2N, соответствующую оценке величины . имеет вид

где

ковариационная матрица размерности 2 x 2, соответствующая оценкам и комплексных величин и , соответственно.

A.3. Ковариационная матрица

размерности 2m x 2m, где

соответствующая оценке

величины , определяется обобщенным законом трансформирования неопределенностей

(A.1)

A.4. — матрица чувствительности размерности 2m x 2N, получаемая вычислением

в точке , где — матрица размерности 2 x 2 частных производных первого порядка действительных и мнимых частей , соответствующих действительной и мнимой частям :

A.5. Для произвольной комплексной скалярной величины рассмотрим матричное представление размерности 2 x 2 вида [14]:

Тогда можно представить как

Данное представление является основой для расчетов частных производных первого порядка комплексных величин по .

Приложение B

(справочное)

ВЫЧИСЛЕНИЕ КОЭФФИЦИЕНТОВ ЧУВСТВИТЕЛЬНОСТИ И КОВАРИАЦИОННЫХ

МАТРИЦ ДЛЯ МНОГОМЕРНЫХ МОДЕЛЕЙ

B.1. Если измерительная задача может быть выражена в терминах линейной алгебры [13], то вычислительно устойчивый алгоритм определения матрицы , являющейся решением уравнения (8), состоит в следующем:

a) для матрицы выполняют разложение Холецкого , в результате чего получают матрицу ;

b) матрицу представляют в виде произведения , где — ортогональная матрица, а — верхняя треугольная матрица;

c) матрицу представляют в виде произведения , где — нижняя треугольная, а — верхняя треугольная матрицы;

d) решают матричное уравнение относительно ;

e) решают матричное уравнение относительно ;

f) вычисляют матрицу ;

g) вычисляют матрицу ;

h) вычисляют матрицу ;

i) матрицу приводят к треугольному виду ;

j) вычисляют .

B.2. Указанная процедура может быть проверена методами элементарной матричной алгебры (см. [7]).

Приложение C

(справочное)

ПРЕОБРАЗОВАНИЕ СИСТЕМЫ КООРДИНАТ

C.1. Основные положения

В настоящем приложении рассматриваются некоторые аспекты задачи преобразования системы координат (см. 9.3). В разделе C.2 приведен аналитический вывод совместной плотности распределения для в случае, когда X₁ описывается нормальным распределением , а X₂ — нормальным распределением и при этом X₁ и X₂ взаимно независимы.

C.2. Аналитическое решение для особого случая

C.2.1. Предположим, что имеет плотность распределения , и — взаимно-однозначное преобразование значений величины в значения величины . Тогда [19, страница 35] имеет плотность распределения

(C.1)

где — детерминант матрицы Якоби ,

рассматриваемый как функция . Предполагается, что нигде не равен нулю или бесконечности.

Примечание 1. Формулу (C.1) иногда называют формулой замены переменных.

Примечание 2. В случае одномерной величины (N = 1) преобразование переменных , где f(.) — дифференцируемая и монотонная функция, дает следующую плотность распределения для Y [21, страницы 57 — 61]:

C.2.2. Для задачи преобразования системы координат, рассмотренной в 9.3, со значениями , со значениями и

Таким образом,

и

Из этого следует, что при 

C.2.3. Рассмотрим случай, когда X₁ описывается нормальным распределением , а X₂ — нормальным распределением и X₁ и X₂ взаимно независимы. Тогда

C.2.4. Маргинальное распределение для Y₁ = R имеет вид

(C.2)

где

а I₀ — модифицированная функция Бесселя первого рода нулевого порядка.

Примечание 1. Полученное распределение представляет собой распределение Райса с параметрами y₁ и u_x.

Примечание 2. Если y₁ = 0, то данное распределение является распределением Рэлея с параметром u_x.

Примечание 3. Если u_x = 1, то данное распределение представляет собой нецентральное распределение хи-квадрат с двумя степенями свободы и параметром нецентральности .

C.2.5. Маргинальное распределение для имеет вид

(C.3)

где

а

дополнительная функция ошибок.

C.2.6. Если, кроме того, выполнено условие x₁ = x₂ = 0, тогда

и, следовательно, Y₁ и Y₂ взаимно независимые величины с маргинальным распределением Рэлея с параметром u_x для Y₁

и маргинальным равномерным распределением на интервале от до для Y₂

C.3. Применение способа оценивания неопределенности по GUM

C.3.1. Для задачи преобразования системы координат, рассмотренной в 9.3, модель измерения может быть записана как двумерная модель

Y₂ = f₂(X₁, X₂) = arctan(X₂/X₁),

при этом подразумевается, что Y₁ >= 0 и . Входные величины X₁ и X₂ имеют оценки x₁ и x₂, соответствующие стандартные неопределенности u(x₁) и u(x₂) и ковариацию u(x₁, x₂).

C.3.2. Из 6.2.1.2 следует, что оценки величин Y₁ и Y₂ имеют вид

y₂ = arctan(x₂/x₁).

C.3.3. Матрицу чувствительности размерности 2 x 2 получают, вычисляя

в точке X₁ = x₁, X₂ = x₂. Таким образом, при условии получаем

C.3.4. Из 6.2.1.3 следует, что

является ковариационной матрицей, соответствующей оценкам , с u(y₂, y₁) = u(y₁, y₂) и

u²(y₁) = u²(x₁)cos²y₂ + u²(x₂)sin²y₂ + 2u(x₁, x₂)cosy₂siny₂,

u(y₁, y₂) = (u²(x₂) — u²(x₁))(siny₂cosy₂)/y₁ +

+ u(x₁, x₂)(cos²y₂ — sin²y₂)/y₁,

C.3.5. В рамках способа оценивания неопределенности по GUM приписывают двумерное нормальное распределение , по которому могут быть построены области охвата для при заданной вероятности охвата p (см. 6.5).

C.3.6. Рассмотрим случай, когда u(x₁) = u(x₂) = u_x и u(x₁, x₂) = 0 (см. C.2.3). Тогда

u(y₁, y₂) = 0, 

с , характеризуемым двумерным нормальным распределением, как в C.3.5. Т.е. в этом случае Y₁ и Y₂ являются независимыми величинами, и двумерное нормальное распределение для Y₁ и Y₂ распадается на два одномерных распределения N (y₁, u²(y₁)) и N(y₂, u²(y₂)), соответственно.

Примечание. Напротив, при аналитическом решении (см. C.2) в случае, когда не характеризуется двумерным распределением Гаусса, выполнения условий u(x₁) = u(x₂) и u(x₁, x₂) = 0 недостаточно для независимости Y₁ и Y₂. Чтобы выходные величины были независимы, необходимо также, чтобы x₁ = x₂ = 0 (см. C.2.6).

Приложение D

(справочное)

ОСНОВНЫЕ ОБОЗНАЧЕНИЯ

Приложение ДА

(справочное)

СВЕДЕНИЯ О СООТВЕТСТВИИ ССЫЛОЧНЫХ МЕЖДУНАРОДНЫХ ДОКУМЕНТОВ

НАЦИОНАЛЬНЫМ СТАНДАРТАМ РОССИЙСКОЙ ФЕДЕРАЦИИ

Таблица ДА.1

БИБЛИОГРАФИЯ

МЕЖГОСУДАРСТВЕННЫЙ СТАНДАРТ

НЕОПРЕДЕЛЕННОСТЬ ИЗМЕРЕНИЯ

Часть 1

Введение в руководства по выражению неопределенности измерения

Uncertainty of measurement. Part 1. Introduction to guides on the expression of uncertainty in measurement

ОКС 17.020

Дата введения 2018-09-01

Предисловие

Цели, основные принципы и основной порядок проведения работ по межгосударственной стандартизации установлены в ГОСТ 1.0-2015 «Межгосударственная система стандартизации. Основные положения» и ГОСТ 1.2-2015 «Межгосударственная система стандартизации. Стандарты межгосударственные, правила и рекомендации по межгосударственной стандартизации. Правила разработки, принятия, обновления и отмены»

Сведения о стандарте

1 ПОДГОТОВЛЕН Межгосударственным техническим комитетом по стандартизации МТК 125 «Статистические методы в управлении качеством продукции» на основе собственного перевода на русский язык англоязычной версии международного документа, указанного в пункте 5

2 ВНЕСЕН Федеральным агентством по техническому регулированию и метрологии (Росстандарт)

3 ПРИНЯТ Межгосударственным советом по стандартизации, метрологии и сертификации (протокол от 14 июля 2017 г. N 101-П)

За принятие проголосовали:

4 Приказом Федерального агентства по техническому регулированию и метрологии от 12 сентября 2017 г. N 1064-ст межгосударственный стандарт ГОСТ 34100.1-2017/ISO/IEC Guide 98-1:2009 введен в действие в качестве национального стандарта Российской Федерации с 1 сентября 2018 г.

5 Настоящий стандарт идентичен международному документу ISO/IEC Guide 98-1:2009* «Неопределенность измерения. Часть 1. Введение в выражение неопределенности измерения» («Uncertainty of measurement — Part 1: Introduction to the expression of uncertainty in measurement», IDT).

________________

* Доступ к международным и зарубежным документам, упомянутым в тексте, можно получить, обратившись в Службу поддержки пользователей. — .

Международный документ разработан Рабочей группой JCGM/WG 1 Объединенного комитета по руководствам в метрологии (как JCGM 104:2009) и одобрен национальными комитетами Международных организаций по стандартизации ISO и IEC.

Официальные экземпляры международного стандарта, на основе которого подготовлен настоящий межгосударственный стандарт, и международных стандартов, на которые даны ссылки, имеются в Федеральном агентстве по техническому регулированию и метрологии.

Наименование настоящего стандарта изменено относительно наименования указанного международного документа для приведения в соответствие с ГОСТ 1.5 (подраздел 3.6).

При применении настоящего стандарта рекомендуется использовать вместо ссылочных международных стандартов и документов соответствующие им межгосударственные стандарты, сведения о которых приведены в дополнительном приложении ДА

6 ВВЕДЕН ВПЕРВЫЕ

7 ПЕРЕИЗДАНИЕ. Июнь 2018 г.

Информация об изменениях к настоящему стандарту публикуется в ежегодном информационном указателе «Национальные стандарты», а текст изменений и поправок — в ежемесячном информационном указателе «Национальные стандарты». В случае пересмотра (замены) или отмены настоящего стандарта соответствующее уведомление будет опубликовано в ежемесячном информационном указателе «Национальные стандарты». Соответствующая информация, уведомление и тексты размещаются также в информационной системе общего пользования — на официальном сайте Федерального агентства по техническому регулированию и метрологии в сети Интернет (www.gost.ru)

Предисловие к международному документу ISO/IEC Guide 98.1:2009

В 1997 г. семью международными организациями, подготовившими в 1993 г. «Руководство по выражению неопределенности измерения» (GUM) и «Международный словарь по метрологии. Основные и общие понятия и связанные с ними термины» (VIM), был образован Объединенный комитет по руководствам в метрологии (JCGM), возглавляемый директором Международного бюро мер и весов (МБМВ), который принял на себя ответственность за указанные два документа от Технической консультативной группы по метрологии ИСО (ИСО/ТАГ 4).

Учредителями JCGM помимо МБМВ являются Международная электротехническая комиссия (МЭК), Международная федерация клинической химии и лабораторной медицины (МФКХ), Международное сотрудничество по аккредитации лабораторий (ИЛАК), Международная организация по стандартизации (ИСО), Международный союз теоретической и прикладной химии (ИЮПАК), Международный союз теоретической и прикладной физики (ИЮПАП) и Международная организация по законодательной метрологии (МОЗМ).

В рамках JCGM созданы две Рабочие группы (РГ). Задачей РГ 1 «Выражение неопределенности измерения» является содействие использованию Руководства (GUM), подготовка дополнений к Руководству и иных документов, способствующих его широкому применению. Задачей РГ 2 «Рабочей группы по Международному словарю основных и общих терминов в метрологии (VIM)» является пересмотр VIM и содействие его применению. Более подробную информацию о деятельности JCGM можно найти на сайте www.bipm.org.

Настоящий документ был подготовлен РГ 1 на основе детальных обзоров, подготовленных организациями — членами JCGM.

Настоящий документ является частью серии документов JCGM под общим названием «Оценивание данных измерений», включающей в себя:

— JCGM 100:2008 Оценивание данных измерений. Руководство по выражению неопределенности измерения (GUM) (см. раздел 2 настоящего стандарта);

— JCGM 101:2008 Оценивание данных измерений. Дополнение 1 к «Руководству по выражению неопределенности измерения». Трансформирование распределений с использованием метода Монте-Карло (см. раздел 2 настоящего стандарта);

— JCGM 102 Оценивание данных измерений. Дополнение 2 к «Руководству по выражению неопределенности измерения». Модели с произвольным числом выходных величин;

— JCGM 103 Оценивание данных измерений. Дополнение 3 к «Руководству по выражению неопределенности измерения». Моделирование;

— JCGM 104 Оценивание данных измерений. Введение к «Руководству по выражению неопределенности измерения» и сопутствующим документам (настоящий стандарт);

— JCGM 105 Оценивание данных измерений. Понятия и основные принципы;

— JCGM 106 Оценивание данных измерений. Роль неопределенности измерения в оценке соответствия;

— JCGM 107 Оценивание данных измерений. Применение метода наименьших квадратов.

Введение

Данные о неопределенности измерения должны всегда приниматься во внимание при оценке соответствия результата измерения его целям. Покупатель в овощной лавке не будет возражать, если при покупке килограмма фруктов весы покажут отклонение от истинного значения в пределах, допустим, двух граммов. В то же время размеры деталей гироскопов, используемых в системах навигации воздушных судов, контролируют до миллионных долей.

Неопределенность измерения — это общее понятие, связанное с любым измерением, которое используют при необходимости принятия обоснованных решений в разных областях практической деятельности и теоретических исследований. По мере наблюдаемого ужесточения допусков в технологических процессах роль неопределенности измерений при оценке соответствия этим допускам все более возрастает. Центральную роль неопределенность измерения играет также при оценке качества и в стандартах качества.

Измерения присутствуют практически во всех видах человеческой деятельности, включая промышленность, торговлю, науку, здравоохранение, обеспечение безопасности и охрану окружающей среды, помогая принимать обоснованные решения. Знание неопределенности измерения позволяет сопоставлять результат измерения с установленными требованиями при оценке соответствия, находить вероятность принятия неправильного решения и с ее учетом управлять возникающими рисками.

Настоящий документ служит введением в концепцию неопределенности измерения, в GUM и сопутствующие документы, указанные в предисловии. Для оценивания неопределенности используется вероятностный подход. Аббревиатуры, использованные в настоящем документе, приведены в приложении А.

В последующих изданиях JCGM 200 (VIM) предполагается дать четкое разграничение в применении термина «погрешность» к величине и к значению величины. То же самое относится к термину «показание». Поскольку в действующем издании JCGM 200:2008 такого разграничения нет, то данный вопрос рассматривается в настоящем документе.

1 Область применения

Настоящий документ подготовлен Объединенным комитетом по руководствам в метрологии (JCGM) с целью продвижения идей оценивания неопределенности измерения, изложенных в «Руководстве по выражению неопределенности измерения» (GUM), и в качестве вводного руководства по применению дополнений к GUM (далее при ссылках — JCGM 100:2008 (GUM)), включая JCGM 101:2008, а также другим документам, разрабатываемым JCGM (см. [3], [4], [5], [6], [7]).

Как и JCGM 100:2008, настоящий документ, в первую очередь, рассматривает выражение неопределенности измерения хорошо определенной величины, характеризуемой единственным истинным значением (JCGM 200:2008 (VIM), словарная статья 2.11, примечание 3) и называемой измеряемой величиной (JCGM 200:2008 (VIM), словарная статья 2.3). В JCGM 100:2008 (GUM) приведены обоснования, почему не рекомендуется использовать термин «истинное значение», однако в настоящем документе этот термин рассматривается для предотвращения возможных неясностей или путаницы с его применением.

Дополнения к GUM и другие сопутствующие документы разрабатываются JCGM с целью оказать помощь в понимании принципов, установленных в GUM, и расширить сферу его применения. Дополнения к GUM вместе с другими сопутствующими документами создают область применения концепции неопределенности измерения, существенно превышающую ту, что установлена GUM.

Настоящий документ знакомит с понятием неопределенности измерения, с GUM и дополнениями к GUM, а также документами, поддерживающими GUM. Он ограничивается преимущественно вопросами измерения величин, которые могут быть охарактеризованы непрерывными переменными, такими как длина, температура, время, количество вещества.

Настоящий документ распространяется на следующие сферы деятельности (но не ограничивается ими):

— наука;

— промышленность;

— деятельность калибровочных и испытательных лабораторий в промышленности, а также в сферах здравоохранения, обеспечения безопасности и охраны окружающей среды;

— деятельность органов по аккредитации, а также органов контроля, надзора и оценки соответствия.

Настоящий документ может быть использован при проектировании изделий, поскольку установление характеристик изделий с учетом последующих требований к контролю и связанными с ним измерениями позволит избежать завышенных технологических требований при их производстве. Применение настоящего документа в сфере высшего образования позволит включать в программы по различным дисциплинам разделы по неопределенности измерения. Результатом должна стать лучшая подготовленность специалистов к пониманию концепции неопределенности измерения и применению ее в разных измерительных задачах, что, в конечном итоге, послужит улучшению качества измерений в целом.

Настоящий документ, дополнения к GUM и другие сопутствующие документы следует использовать совместно с Международным словарем по метрологии (VIM, далее при ссылках — JCGM 200:2008 (VIM)), а также с международными стандартами ISO 3534-1:2006, ISO 3524-2:2006 и ISO 3534-3:1999, в которых определены термины, используемые в математической статистике и теории вероятностей (включая прикладную статистику и планирование экспериментов), и показано их место в структуре понятий в соответствии с нормативной терминологической практикой. Последнее важно с учетом того обстоятельства, что теоретической основой оценивания данных измерений и неопределенности измерений является математическая статистика и теория вероятностей.

2 Нормативные ссылки

В настоящем стандарте использованы нормативные ссылки на следующие стандарты и документы*:

________________

* Таблицу соответствия национальных стандартов международным см. по ссылке. — .

JCGM 100:2008, Evaluation of measurement data — Guide to the expression of uncertainty in measurement (GUM) (Оценивание данных измерений. Руководство по выражению неопределенности измерения (GUM))

JCGM 101:2008, Evaluation of measurement data — Supplement 1 to the «Guide to the expression of uncertainty in measurement» — Propagation of distributions using a Monte Carlo method (Оценивание данных измерений. Дополнение 1 к «Руководству по выражению неопределенности измерения». Трансформирование распределений с использованием метода Монте-Карло)

JCGM 200:2008, International Vocabulary of Metrology — Basic and general concepts and associated terms, 3rd Edition (Международный словарь по метрологии. Основные и общие понятия и связанные с ними термины, 3-е издание)

ISO 3534-1:2006, Statistics — Vocabulary and symbols — Part 1: General statistical terms and terms used in probability (Статистика. Словарь и условные обозначения. Часть 1. Общие термины по статистике и термины, используемые в теории вероятностей)

ISO 3534-2:2006, Statistics — Vocabulary and symbols — Part 2: Applied statistics (Статистика. Словарь и условные обозначения. Часть 2. Прикладная статистика)

ISO 3534-3:1999, Statistics — Vocabulary and symbols — Part 3: Design of experiments (Статистика. Словарь и условные обозначения. Часть 3. Планирование эксперимента)

_______________

Данный международный стандарт заменен на международный стандарт ISO 3534-3:2013 «Statistics — Vocabulary and symbols — Part 3: Design of experiments». Однако для однозначного соблюдения требования настоящего стандарта, выраженного в датированной ссылке, рекомендуется использовать только указанное в этой ссылке издание.

3 Понятие неопределенности измерения

3.1 Цель измерения состоит в получении информации об интересующей величине, называемой измеряемой величиной (JCGM 200:2008 (VIM), словарная статья 2.3). Измеряемой величиной может быть объем сосуда, разность потенциалов на клеммах батареи или массовая концентрация свинца в колбе с водой.

3.2 Абсолютно точных измерений не существует. При проведении измерения его результат зависит от измерительной системы (JCGM 200:2008 (VIM), словарная статья 3.2), методики измерения, квалификации оператора, внешних условий и других факторов [1]. Так, если измерять одну и ту же величину несколько раз одним способом и в одинаковых условиях, то, как правило, при достаточной разрешающей способности измерительной системы, позволяющей различать близкие показания (JCGM 200:2008 (VIM), словарная статья 4.1), эти показания (полученные значения измеряемой величины [JCGM 200:2008 (VIM), словарная статья 2.10]) всякий раз будут разными. Показания рассматривают как мгновенные значения соответствующей случайной переменной (показываемой величины).

3.3 Разброс показаний позволяет судить о качестве проведенного измерения. Их среднее должно обеспечить значение оценки (ISO 3534-1:2006, словарная статья 1.31) истинного значения величины (JCGM 200:2008 (VIM), словарная статья 2.11), которая в общем случае будет более достоверной, чем отдельное показание. Разброс показаний и их число дают некоторую информацию в отношении среднего значения как оценки истинного значения величины. Однако эта информация в большинстве случаев не будет достаточной.

3.4 Измерительная система может давать показания, которые рассеяны не вокруг истинного значения величины, а вокруг некоторого другого, смещенного значения. Разницу между смещенным значением и истинным значением величины иногда называют значением систематической погрешности (JCGM 200:2008 (VIM), словарная статья 2.17). Возьмем для примера домашние весы в ванной. Предположим, что в отсутствие нагрузки они показывают не ноль, а некоторое отличное от нуля значение. Тогда вне зависимости от числа повторных измерений массы встающего на весы человека влияние этого смещения будет неизменно присутствовать в среднем значении показаний. В большинстве случаев систематическая погрешность, рассматриваемая как величина, — это составляющая погрешности, которая остается постоянной или зависит определенным образом от какой-то другой величины.

3.5 Существует два вида погрешности измерения: систематическая и случайная (JCGM 200:2008 (VIM), словарная статья 2.19). Систематическая погрешность [значение оценки которой называют смещением при измерении (JCGM 200:2008 (VIM), словарная статья 2.18)] проявляется в том, что полученное значение измеряемой величины содержит сдвиг. Случайная погрешность проявляется в том, что при повторении измерения полученное значение измеряемой величины в большинстве случаев будет отличаться от предыдущего. Случайность заключается в том, что последующие значения измеряемой величины нельзя точно предсказать по предыдущим (если бы такая возможность существовала, то в результат измерений можно было бы внести соответствующую поправку). В общем случае каждый из видов погрешности может быть обусловлен действием нескольких факторов.

3.6 Для каждого проведенного измерения необходимо решить, как наилучшим образом представить информацию, которую удалось получить об измеряемой величине. Указание значений систематических и случайных погрешностей наряду с наилучшей оценкой измеряемой величины — это тот подход, который часто использовался до разработки GUM. GUM предложило другой подход к пониманию измерения, в частности, к тому, как выражать качество результата измерения. Вместо представления результата измерения в виде наилучшей оценки измеряемой величины вместе с информацией о систематической и случайной погрешностях (в форме «анализа погрешностей»), GUM рекомендует выражать результат измерения как наилучшую оценку измеряемой величины вместе с соответствующей неопределенностью измерения.

3.7 Одним из основных исходных положений подхода GUM является утверждение о возможности охарактеризовать качество измерения, исходя из единообразного обращения с систематической и случайной погрешностями, с предложением метода, как это сделать (см. 7.2). Этот метод возвращает к исходной информации, какой она была до применения «анализа погрешностей», и подводит под нее вероятностную основу с помощью концепции неопределенности измерения.

3.8 Другое базовое положение GUM состоит в утверждении, что нельзя установить, насколько хорошо известно единственное истинное значение величины, а можно только сформулировать степень нашей уверенности в том, что оно известно. Таким образом, неопределенность измерения можно представить через степень уверенности. Такая неопределенность будет отражать неполноту знания об измеряемой величине. Понятие «уверенности» очень важно, так как оно перемещает метрологию в сферу, где результат измерения должен рассматриваться и численно определяться в терминах вероятностей, которые выражают степень доверия.

3.9 Все сказанное выше касается прямого измерения величины, которое встречается довольно редко. Так весы в ванной комнате могут преобразовывать измеренное растяжение пружины в оценку измеряемой величины — массы человека на весах. Соотношение между растяжением данной пружины и массой определяют с помощью калибровки (JCGM 200:2008 (VIM), словарная статья 2.39) весов.

3.10 Соотношение, подобное тому, что описано в 3.9, устанавливает правило преобразования численного значения некоторой величины в соответствующее значение измеряемой величины. Это правило обычно называют моделью измерений (JCGM 200:2008 (VIM), словарная статья 2.48) или просто моделью. На практике встречаются измерения разных видов, и им соответствуют разные правила преобразования или модели. Даже одному конкретному виду измерений может соответствовать несколько моделей. Так для бытовых измерений может быть достаточна простая модель (например, в виде прямо пропорциональной зависимости массы на весах от растяжения пружины). Тогда как для научных целей или на производстве для получения более точных результатов могут использоваться более сложные модели взвешивания, учитывающие дополнительные факторы, например выталкивающую силу воздуха. Как правило, определение измеряемой величины зависит от ряда других величин, таких как температура, влажность, смещение объекта, которые также необходимо измерять.

3.11 Если условия измерений несколько отличаются от заданных, то в величины, входящие в модель, должны быть внесены поправки, соответствующие значениям систематической погрешности (JCGM 200:2008 (VIM), словарная статья 2.17). Если поправку можно оценить, то соответствующую величину следует скорректировать на полученное значение оценки [см. JCGM 100:2008 (GUM) (3.2.4)]. Это внесет дополнительную неопределенность в результат измерения, даже если значение оценки, как это часто случается на практике, будет равно нулю. Примерами источников систематических погрешностей, возникающих при измерениях высоты, могут быть отклонение средства измерений от вертикали или отличие от предписанного значения температуры окружающей среды. Ни угол отклонения средства измерений, ни температуру окружающей среды нельзя узнать точно, но можно получить некоторую информацию о возможных значениях этих величин, например, что угол отклонения от вертикали не может превышать 0,001° или что температура окружающей среды во время измерения отличается от предписанной не более чем на 2°С.

3.12 Величина, входящая в модель измерения, может зависеть от времени, например, если она отражает распад радионуклида с определенной скоростью. В этом случае соответствующая временная зависимость должна быть включена в модель, чтобы дать возможность соотнести измеряемую величину с временем проведения измерения.

3.13 Зачастую модель измерения предполагает использование помимо результатов наблюдений входящих в нее случайных величин также данных иной природы, в частности, физических констант, известных с некоторой точностью. Примерами таких констант могут служить физические характеристики материалов, например модуль упругости или удельная теплоемкость. Также в модель в качестве значений оценок величин могут быть включены данные, заимствованные из справочников, сертификатов о калибровке и других аналогичных источников.

3.14 Составляющие модели, необходимые для определения измеряемой величины, называют входными величинами модели измерений (JCGM 200:2008 (VIM), словарная статья 2.50). Саму модель, определяющую правило преобразования входных величин, часто называют функциональной зависимостью [см. JCGM 100:2008 (GUM) (4.1)]. Выходной величиной модели измерений (JCGM 200:2008 (VIM), словарная статья 2.51) является измеряемая величина.

3.15 Формально, связь выходной величины, обозначаемой Y, в отношении которой требуется получить информацию, с входными величинами, обозначаемыми , …, , информация о которых доступна, часто представляют моделью [см. JCGM 100:2008 (GUM), (4.1.1)] в виде функции измерения (JCGM 200:2008 (VIM), словарная статья 2.49)

. (1)

3.16 В общем виде модель измерения (см. JCGM 200:2008 (VIM), примечание 1 к словарной статье 2.48) может быть представлена формулой

. (2)

Предполагается, что для модели, задаваемой формулой (2), существует способ вычисления Y по данным , …, и что получаемое при этом значение Y единственно.

3.17 Истинные значения входных величин ,…, неизвестны. В подходе, принятом GUM, , …, ассоциируют со случайными переменными (ISO 3534-1:2006, словарная статья 2.10) с соответствующими распределениями вероятностей [см. JCGM 100:2008 (GUM) (3.3.5), а также ISO 3534-1:2006, словарную статью 2.11]. Эти распределения, принимаемые на основе имеющихся знаний об , …, , описывают вероятности нахождения истинных значений входных величин в разных интервалах. Иногда входные величины (все или некоторые) могут быть связаны между собой, и для их описания используют совместные распределения. В настоящем документе рассматриваются, преимущественно, независимые случайные переменные, однако полученные выводы могут быть легко обобщены и на случай взаимосвязанных величин.

3.18 Если из сертификатов, отчетов, документации изготовителей, анализа данных измерений и других источников известны значения оценок , …, соответствующих входных величин , …, , то ассоциированные с , …, распределения вероятностей должны иметь значения , …, в качестве своих математических ожиданий [см. JCGM 101:2008 (3.6), а также ISO 3534-1:2006, словарную статью 2.12]. Для каждого значения оценки i-й входной величины существует ассоциированная с ней стандартная неопределенность (JCGM 200:2008 (VIM), словарная статья 2.30), обозначаемая и определяемая как стандартное отклонение [см. JCGM 101:2008 (3.8), а также ISO 3534-1:2006, словарную статью 2.37] входной величины, ассоциированной с . Значение оценки понимают как наилучшее для данной входной величины в том смысле, что будет меньше, чем математическое ожидание квадрата отклонения от любого другого значения.

3.19 Принцип использования всей доступной информации для установления распределения вероятностей, характеризующих входящую в модель величину, справедлив как для каждой входной величины , так и для выходной величины Y. В последнем случае распределение вероятностей определяют на основе функциональной зависимости (1) или (2) и известных распределений вероятностей для . Данный способ получения распределения вероятностей для Y известен как трансформирование распределений [см. JCGM 101:2008 (5.2)].

3.20 Априорное знание об истинном значении выходной величины Y также может быть использовано соответствующим образом. Так в отношении измерений на домашних весах в ванной комнате априорными будут сведения, что масса человека на весах положительна, и что измеряют массу именно человека, а не, например, автомобиля. Учет такой дополнительной информации может помочь обоснованно выбрать распределение вероятностей для Y с меньшим стандартным отклонением, что, соответственно, даст меньшую стандартную неопределенность, ассоциированную со значением оценки Y ([2], [13], [24]).

4 Основные понятия и принципы

4.1 Основные понятия и принципы теории вероятностей, которые положены в основу концепции неопределенности измерения, изложенной в разделе 3, представлены в [4].

4.2 Неопределенность измерения определяют как (JCGM 200:2008 (VIM), словарная статья 2.26) «неотрицательный параметр, характеризующий рассеяние значений величины, приписываемых измеряемой величине на основании используемой информации«.

Это определение согласуется с положениями, изложенными в 3.8, а также в 3.17-3.20.

4.3 При вычислении неопределенности используются два представления распределения вероятностей [см. JCGM 101:2008 (3.1), а также ISO 3534-1:2006, словарную статью 2.11] случайной переменной X:

— через функцию распределения [см. JCGM 101:2008 (3.2), а также ISO 3534-1:2006, словарную статью 2.7], дающую для любого значения ее аргумента вероятность того, что Х меньше или равна этому значению;

— через функцию плотности вероятностей [см. JCGM 100:2008 (3.3), а также ISO 3534-1:2006, словарную статью 2.26], являющуюся производной от функции распределения.

4.4 Информацию о каждой входной величине в модели измерений, как правило, представляют в виде наилучшего значения оценки и ассоциированной с ней стандартной неопределенностью (см. 3.18). Если для любых i и j и связаны между собой (зависимы), то соответствующая информация должна быть отражена в виде меры тесноты этой связи, выражаемой через ковариацию (ISO 3534-1:2006, словарная статья 2.43) или корреляцию случайных переменных. Если и не связаны между собой (независимы), то соответствующая ковариация будет равна нулю.

4.5 Оценивание данных измерения в контексте модели измерений (1) или (2) — это использование имеющейся информации о входных величинах , …, для получения ассоциированных с ними распределений вероятностей и последующего вывода распределения вероятностей, ассоциированного с выходной величиной Y. Последнее распределение, таким образом, можно рассматривать как результат оценивания данных измерения.

4.6 Информация о входной величине в модели измерений может быть получена из повторных показаний (оценивание неопределенности по типу А) [см. JCGM 100:2008 (GUM) (4.2), а также JCGM 200:2008 (VIM), словарную статью 2.28] или из обоснованных суждений на основе имеющихся данных о возможных значениях этой величины (оценивание неопределенности по типу В) [см. JCGM 100:2008 (GUM) (4.3), а также JCGM 200:2008 (VIM), словарную статью 2.29].

4.7 При оценивании неопределенности по типу A (JCGM 200:2008 (VIM), словарная статья 2.28) часто делают предположение, что распределение, наилучшим образом соответствующее входной величине X в условиях имеющихся повторных независимых показаний, это распределение Гаусса (ISO 3534-1:2006, словарная статья 2.50). В таком случае Х характеризуется математическим ожиданием, наилучшей оценкой которого является среднее арифметическое показаний, и стандартным отклонением, равным стандартному отклонению среднего арифметического. Если неопределенность оценивают по малому числу показаний (являющихся мгновенными реализациями величины, распределенной по нормальному закону), то соответствующим распределением будет t-распределение (ISO 3534-1:2006, словарная статья 2.53). На рисунке 1 показаны плотности вероятности для распределения Гаусса (сплошная линия) и t-распределения с четырьмя степенями свободы (пунктирная линия). Сказанное выше не будет справедливо, если показания нельзя рассматривать как независимые.

4.8 При оценивании неопределенности по типу В (JCGM 200:2008 (VIM), словарная статья 2.29) часто единственной доступной информацией является то, что Х лежит в определенном интервале [а, b]. Информация такого вида может быть формализована в виде прямоугольного распределения вероятностей [см. JCGM 100:2008 (GUM) (4.3.7), а также ISO 3534-1:2006, словарную статью 2.60] с границами а и b (рисунок 2). Если бы о рассматриваемой величине была доступна информация иного рода, то распределение вероятностей должно было быть согласовано с этой имеющейся информацией [26].

4.9 После того, как составлена модель измерения, и входные величины , …, описаны через соответствующие распределения вероятностей, распределение вероятностей для измеряемой величины Y полностью определено (см. также 3.19). Математическое ожидание Y используется в качестве оценки измеряемой величины, а стандартное отклонение Y — в качестве стандартной неопределенности, ассоциированной с этой оценкой.

4.10 На рисунке 3 показано трансформирование двух разных прямоугольных распределений вероятностей для входных величин и в симметричное трапецеидальное распределение вероятностей для измеряемой величины Y в случае аддитивной функции измерения .

4.11 Часто необходимо знать интервал, содержащий Y с заданной вероятностью. Такой интервал, называемый интервалом охвата (JCGM 200:2008 (VIM), словарная статья 2.36), может быть получен из распределения вероятностей для Y. Заданную вероятность называют вероятностью охвата (JCGM 200:2008 (VIM), словарная статья 2.37).

4.12 Для установленной вероятности охвата существует множество интервалов охвата, среди которых различают:

a) вероятностно симметричный интервал охвата [см. JCGM 101:2008 (3.15)], для которого вероятности (в сумме равные единице за вычетом вероятности охвата) расположения значения величины справа или слева от интервала равны;

b) наименьший интервал охвата [см. JCGM 101:2008 (3.16)], протяженность которого является наименьшей из всех интервалов охвата, имеющих ту же вероятность охвата.

4.13 На рисунке 4 показано усеченное и масштабированное распределение Гаусса (в виде спадающей кривой) с граничными точками наименьшего (сплошные вертикальные линии) и вероятностно симметричного (пунктирные вертикальные линии) 95%-ных интервалов охвата для величины, с которой ассоциировано это распределение. Распределение асимметрично, поэтому указанные два интервала охвата различаются между собой (особенно заметно различие в граничных точках справа). Левая граничная точка наименьшего интервала охвата точно совпадает с нулем — наименьшим возможным значением для этой величины. Для данного примера вероятностно симметричный интервал охвата на 15% протяженней наименьшего интервала охвата.

Рисунок 1 — Распределение Гаусса (сплошная линия) и t-распределение с четырьмя степенями свободы (пунктирная линия) (для случайной переменной размерности D размерность плотности распределения будет )

Рисунок 2 — Прямоугольное распределение на интервале [-0,10; 0,10] (для случайной переменной размерности D размерность плотности распределения будет )

Рисунок 3 — Трансформирование распределений для аддитивной функции измерения при прямоугольных распределениях вероятностей для входных величин

Рисунок 4 — 95%-ные интервалы охвата: наименьший (сплошные вертикальные линии) и вероятностно симметричный (пунктирные вертикальные линии) для величины с усеченным масштабированным распределением Гаусса (для случайной переменной размерности D размерность плотности распределения будет )

4.14 Коэффициенты чувствительности , …, [см. JCGM 100:2008 (GUM) (5.1.3)] показывают, как на значение оценки y величины Y будут влиять небольшие изменения в значениях оценок , …, входных величин , …, . Для функции измерения (1) равен частной производной первого порядка от по в точке , и т.д. Если функция измерения линейна:

, (3)

то при независимых , …, изменение значения на приведет к изменению значения у на . То же самое будет справедливо в некотором приближении для большинства моделей, описываемых формулами (1) и (2) (см. 7.2.4). Сравнение значений для разных i позволяет оценить вклад каждой входной величины в стандартную неопределенность , ассоциированную с у.

4.15 Стандартную неопределенность , ассоциированную со значением оценки у выходной величины Y, получают суммированием не самих значений , а их квадратов, т.е.

. (4)

Формула (4) будет справедлива в некотором приближении для большинства моделей измерения, определяемых формулами (1) и (2).

4.16 Если входные величины взаимозависимы, то формулу (4) следует дополнить слагаемыми, содержащими ковариации [см. JCGM 100:2008 (GUM) (5.2.2)], которые могут увеличить или уменьшить значение .

_______________

Пункт 4.17, относящийся к использованию десятичного разделителя только в отношении англоязычной версии документа, исключен.

5 Этапы оценивания неопределенности

5.1 Основные этапы оценивания неопределенности включают в себя формулировку измерительной задачи и вычисления. Последнее включает в себя трансформирование распределений вероятностей и получение окончательного результата.

5.2 Этап формулировки измерительной задачи (см. раздел 6) включает в себя:

a) определение выходной величины Y (измеряемой величины);

b) выявление входных величин, от которых зависит Y;

c) составление модели измерения, определяющей соотношение Y с входными величинами;

d) приписывание распределений вероятностей (нормального, прямоугольного и т.д.) входным величинам (или совместного распределения вероятностей входным величинам, не являющимся независимыми) на основе имеющейся информации.

5.3 Этап вычислений (см. раздел 7) состоит из трансформирования по данной модели измерения распределений вероятностей для входных величин в распределение вероятностей для выходной величины Y и использования этого распределения для получения:

a) математического ожидания Y, принимаемого как значение оценки у величины Y;

b) стандартного отклонения величины Y, принимаемого как стандартная неопределенность , ассоциированная с у [см. JCGM 100:2008 (GUM) (E.3.2)];

c) интервала охвата, содержащего Y с заданной вероятностью охвата.

6 Составление модели измерений

6.1 Этап формулировки измерительной задачи при оценивании неопределенности включает в себя разработку модели измерений, учет соответствующих поправок и других воздействий, если это необходимо. В некоторых областях измерений выполнение данного этапа может представлять сложность. Он также включает в себя использование доступной информации для описания входных величин модели через распределения вероятностей. В [6] приведено руководство по разработке и применению модели измерений. Приписывание распределений вероятностей входным величинам модели измерений рассмотрено в JCGM 101:2008 и в [5].

6.2 В начале составляют модель, связывающую выходную величину с входными величинами. В некоторых задачах выходных величин может быть более одной (см. 6.5). Модель формируют на основе теоретических и/или эмпирических знаний с учетом специфики измерительной задачи (например, измерения электрических параметров, линейных размеров, температуры, массы). Затем модель дополняют другими входными величинами, посредством которых описывают эффекты случайного и систематического влияния на результат измерения. Руководство по учету дополнительных входных величин приведено в [6].

6.3 Класс моделей, рассматриваемых в [6], более широк, чем в GUM, и включает в себя классификацию по следующим признакам:

a) по виду входящих в модель величин: действительные или комплексные;

b) по виду модели: в виде функции измерений [формула (1)] или в общем виде [формула (2)];

c) по числу выходных величин: одна или более (см. 6.5).

Комплексные величины, указанные в перечислении а), используются, главным образом, в измерениях электрических величин, в акустике и оптике. Для функции измерений, указанной в перечислении b), выходная величина выражается непосредственно в виде формулы, в которую входят величины, в то время как модель измерения в общем виде представляет собой уравнение, которое необходимо решить относительной выходной величины (см. 6.5).

6.4 Разнообразные варианты применения [6] проиллюстрированы на примерах из разных областей метрологии. Кроме того, в этом документе приведено руководство по разным аспектам численного анализа в связи с рассматриваемыми примерами. Документ также включает в себя рассмотрение вопросов замены переменных таким образом, чтобы устранить или уменьшить корреляцию входящих в модель величин.

6.5 В GUM и JCGM 101:2008 рассматриваются, в основном, модели измерений в виде функций измерения, в которых есть только одна выходная величина Y. Однако существует множество измерительных задач, в которых необходимо рассматривать несколько выходных величин , …, , зависящих от одних и тех же входных величин. Приведенные в [6] примеры включают в себя случаи, когда а) выходная величина является комплексной и представлена в виде действительной и мнимой частей (или амплитуды и фазы); b) выходные величины представляют собой параметры калибровочной характеристики; с) выходные величины описывают геометрию поверхности объекта. Хотя подобные вопросы затрагиваются в GUM при рассмотрении примеров одновременного измерения активного и реактивного сопротивления [JCGM 100:2008 (GUM) (раздел Н.2)] и калибровки термометра [JCGM 100:2008 (GUM) (раздел Н.3)], специальному анализу в GUM они не подвергаются.

6.6 Этап формулировки измерительной задачи при оценивании неопределенности для случая с более чем одной измеряемой величиной мало отличается от аналогичного этапа для модели измерения с единственной измеряемой величиной. Он включает в себя разработку модели и приписывание распределений вероятностей входным величинам на основе доступной информации. Как и для модели измерений с одной выходной величиной, существует оценка каждой входной величины и стандартной неопределенности, ассоциированные с этой оценкой (и возможные ковариации, ассоциированные с парами оценок). Но так как в общем случае каждая выходная величина зависит от всех входных величин, то в дополнение к определению оценок этих выходных величин и ассоциированных с ними стандартных неопределенностей необходимо будет оценивать ковариации, ассоциированные со всеми парами выходных оценок.

6.7 Эквивалентом функции измерения (1) для произвольного числа m выходных величин являются формулы

, , (5)

для m функций , …, . Схематично формула (5) изображена рисунком 5.

Рисунок 5 — Функция измерения с тремя входными величинами , и и двумя выходными величинами и

6.8 В [6] рассматриваются также модели многоступенчатого измерения, в которых выходные величины предшествующих ступеней становятся входными величинами для последующих ступеней. Типичным примером модели многоступенчатого измерения может служить построение и применение калибровочной характеристики (JCGM 200:2008 (VIM), словарная статья 2.39) (см. рисунок 6):

a) параметры калибровочной характеристики оценивают, сравнивая размеры единицы измерения, переданные от эталонов, с соответствующими показаниями измерительной системы. Стандартные неопределенности, ассоциированные с полученными значениями измеряемой величины и значениями показаний, являются источниками стандартных неопределенностей для значений оценок параметров калибровочной характеристики и, в общем случае, ковариаций для оценок этих параметров;

b) полученное измерительной системой показание по калибровочной характеристике преобразуют в значение измеряемой величины. Для этого используется функция, обратная калибровочной характеристике. Стандартные неопределенности и ковариации, ассоциированные со значениями оценок параметров калибровочной характеристики, вместе со стандартной неопределенностью, ассоциированной со значением очередного показания, являются источниками для расчета стандартной неопределенности, ассоциированной с полученным значением измеряемой величины.

7 Трансформирование распределений и вычисление значений оценок

7.1 Общие положения

7.1.1 Этап вычислений включает в себя процедуру, известную как трансформирование распределений [см. JCGM 101:2008, (5.2)], которая может быть реализована следующими способами:

а) в виде используемого в GUM закона трансформирования неопределенностей с описанием случайной переменной, ассоциированной с выходной величиной Y, распределением Гаусса или t-pacпределением (см. 7.2);

Рисунок 6 — Модель двухступенчатого измерения, включающего построение калибровочной характеристики и ее применение к показаниям измерительной системы

b) в виде аналитического вывода формы распределения вероятностей для Y методами математического анализа (см. 7.3);

c) с помощью метода Монте-Карло, в котором приближенную функцию распределения для Y получают численным моделированием, генерируя случайные значения из распределений вероятностей для входных величин и преобразуя их в значения измеряемой величины посредством модели измерений (см. 7.4).

7.1.2 Для конкретной задачи оценивания неопределенности измерений может быть использован любой из способов, перечисленных в 7.1.1 (или какой-нибудь иной способ), причем способ а) является в большинстве случаев приближенным, способ b) — точным, а способ с) дает решение с числовой точностью, которую можно контролировать.

7.1.3 Применение способов а) и с) к функциям измерения для общеупотребительных моделей измерения с любым числом входных величин рассматривается в 7.5.

7.2 Способ расчета неопределенности по GUM

7.2.1 Способ оценивания неопределенности по GUM [см. JCGM 100:2008 (GUM) (3.4.8) и (5.1)] (схематично показанный на рисунке 7) для получения значения оценки у выходной величины Y и ассоциированной с ней стандартной неопределенности использует:

a) наилучшие значения оценок входных величин ;

b) стандартные неопределенности , ассоциированные с ;

c) коэффициенты чувствительности (см. 4.14).

7.2.2 Способ, указанный в 7.2.1, несколько видоизменяется [см. JCGM 100:2008 (GUM) (5.2)], если входные величины являются взаимозависимыми (на рисунке 7 такая модификация не показана). Если случайная переменная, ассоциированная с выходной величиной Y, имеет распределение Гаусса, то это позволяет построить интервал охвата для Y с заданной вероятностью охвата [см. JCGM 100:2008 (GUM) (раздел G.2)]. Если каждому соответствует конечное число степеней свободы [ISO 3534-1:2006, словарная статья 2.54], то по ним можно рассчитать число эффективных степеней свободы для , а выходную величину Y ассоциировать с t-распределением.

7.2.3 Для многих измерительных ситуаций способ расчета неопределенности по GUM [см. JCGM 100:2008 (GUM) (раздел 5)] позволяет получить достоверные результаты. Если функция измерения линейна относительно входных величин и эти величины распределены по нормальному закону, то способ оценивания неопределенности по GUM дает точные результаты [см. JCGM 101:2008 (5.7)]. Но даже если указанные условия не соблюдаются, данный способ может достаточно хорошо работать на практике [см. JCGM 101:2008 (5.8)].

7.2.4 Однако существуют измерительные ситуации, при которых способ оценивания неопределенности по GUM нельзя считать удовлетворительным. Так будет, в том числе, если:

a) функция измерения нелинейна;

b) распределения вероятностей для входных величин асимметричны;

c) , …, , дающие вклад в неопределенность (см. 4.14), не являются величинами приблизительно одного порядка [см. JCGM 100:2008 (GUM) (G.2.2)];

Рисунок 7 — Способ расчета неопределенности по GUM (левая часть рисунка, выделенная пунктирной линией, относится к получению значения оценки у и ассоциированной с ней стандартной неопределенности , остальная — к получению интервала охвата для Y)

d) распределение вероятностей для выходной величины либо асимметрично, либо существенно отличается от нормального распределения или t-распределения.

Иногда заранее трудно решить, позволяет ли данная измерительная задача применять способ оценивания неопределенности по GUM.

7.2.5 Использование способа оценивания неопределенности по GUM усложняется при нахождении частных производных (или их численных приближений) для сложной модели измерений, что является необходимым для применения закона трансформирования неопределенностей, особенно, если необходимо рассчитывать производные высших порядков [см. JCGM 100:2008 (GUM) (раздел 5)]. В таких случаях более подходящим и удобным для применения является метод Монте-Карло (см. 7.4).

7.3 Аналитический вывод

7.3.1 Аналитические методы, с помощью которых может быть получена алгебраическая формула для распределения вероятностей выходной величины, не содержат никаких приближений, но могут быть применены только в сравнительно простых случаях. В [8], [12] показаны возможности применения таких методов. В число измерительных задач, для которых возможен аналитический вывод, входят те, где выходная величина является линейной функцией входных величин [см. формулу (3)], которые все распределены либо по нормальному закону, либо по прямоугольному закону в одних и тех же границах. Пример для двух входных величин (N=2) с прямоугольными распределениями вероятности, которые дают трапецеидальное распределение выходной величины (см. [10]), показан на рисунке 3.

7.3.2 Часто аналитический вывод возможен для случаев, когда модель измерения включает в себя только одну входную величину (N=1) (см. [25, с.57-61]). Такие случаи возникают при преобразовании единиц измерения, например, из линейных в логарифмические (см. [10, с. 95-98]).

7.3.3 Преимуществом аналитического вывода является то, что он дает возможность понять суть измерения, показывая зависимость распределения вероятностей выходной величины от параметров распределений вероятностей входных величин.

7.4 Метод Монте-Карло

7.4.1 JCGM 101:2008 устанавливает подробное руководство по методу Монте-Карло как способу трансформирования распределений [см. JCGM 101:2008 (5.9)]. Для метода Монте-Карло существует меньше ограничений по применению, чем для способа оценивания неопределенности по GUM [см. JCGM 101:2008 (5.10)]. Схематично метод показан на рисунке 8. В JCGM 101:2008 приведены примеры сравнения метода Монте-Карло со способом оценивания неопределенности по GUM [см. JCGM:2008 101 (раздел 9)].

7.4.2 JCGM 101 устанавливает адаптивную процедуру для метода Монте-Карло, в которой число испытаний определяется автоматически с использованием меры сходимости всего процесса в целом [см. JCGM 101:2008 (7.9)].

Рисунок 8 — Оценивание неопределенности методом Монте-Карло (левая часть рисунка, выделенная пунктирной линией, относится к получению значения оценки у и ассоциированной с ней стандартной неопределенности , остальная — к получению интервала охвата для Y)

7.4.3 В JCGM 101:2008 показано, как метод Монте-Карло может быть применен, чтобы решить, приемлемо ли применение способа оценивания неопределенности по GUM для каждого конкретного случая [см. JCGM 101:2008 (раздел 8)].

7.5 Модели измерения с произвольным числом выходных величин

7.5.1 В случае измерений с использованием моделей с произвольным числом выходных величин способы оценивания ассоциированных с ними неопределенностей и ковариаций, установленные как в GUM, так и в JCGM 101:2008, требуют соответствующего обобщения. Суть такой модификации показана в GUM на ряде примеров [см. JCGM 100:2008 (GUM) (F.1.2.3)].

7.5.2 В [5] установлено, что закон трансформирования неопределенностей, составляющий основное содержание способа оценивания неопределенности по GUM для модели измерения с одной выходной величиной, может быть кратко представлен в матричной форме. Преимущество матричного представления состоит в том, что оно удобно для программной реализации метода, а также легко допускает обобщение на другие модели измерения.

7.5.3 Указанное обобщение использовано в [5] для случая функции измерения с произвольным числом выходных величин. Аналогичное обобщение рассматривается в указанном документе и для случая модели измерения, представленной в общем виде (см. 3.16).

7.5.4 Документ [5] также рассматривает применение для модели измерения с произвольным числом выходных величин метода Монте-Карло. В нем дается вывод распределений вероятностей дискретного вида для выходных величин. На основе этих распределений получены формулы для значений оценок выходных величин, стандартных неопределенностей, ассоциированных с этими оценками, и ковариаций, ассоциированных с парами этих оценок.

7.5.5 Требования к представлению результата измерения могут включать в себя, помимо указания значений оценок выходных величин вместе с ассоциированными стандартными неопределенностями и ковариациями, указание области, содержащей выходные величины с заданной вероятностью (охвата). Такие области естественно было бы рассматривать как обобщения для вероятностно симметричного интервала охвата и наименьшего интервала охвата. Но если для наименьшего интервала охвата его пространственный аналог существует (хотя его построение и сопряжено со значительными трудностями), то этого нельзя сказать в отношении области, аналогичной вероятностно симметричному интервалу охвата, которая не может быть определена единственным образом.

7.5.6 В ряде случаев целесообразно указывать приближенную область охвата, имеющую простую геометрическую форму. С этой точки зрения рассматриваются две формы области охвата. Первая вытекает из ассоциирования выходных величин с совместным распределением Гаусса, например, на основе использования центральной предельной теоремы [см. JCGM 100:2008 (GUM) (раздел G.2)]. Тогда наименьшая область охвата будет иметь вид многомерного эллипсоида. Другой формой является многомерный параллелепипед. В [5] приведены способы построения наименьших областей охвата указанных двух форм.

8 Применение неопределенности измерения при оценке соответствия

8.1 Оценка соответствия — важный аспект управления качеством производства, метрологического надзора, проверки соответствия требованиям безопасности и санитарным нормам. Так при контроле качества деталей на производстве принимают решения о соответствии деталей техническим условиям. Аналогичные вопросы возникают при проверке соответствия законодательно установленным нормативам (например, по выбросам, уровню радиации, содержанию химических веществ, наличию следов допинга). Руководство, в котором рассматриваются подобные вопросы, приведено в [7] (см. также [18]).

8.2 Измерение является неотъемлемой частью оценки соответствия, когда необходимо решить, соответствует ли выходная (измеряемая) величина установленному требованию. Для единственной величины такое требование обычно принимает вид границ, определяющих интервал допустимых значений величины. При отсутствии неопределенности полученное значение измеряемой величины, лежащее в пределах границ, считают соответствующим требованиям, в противном случае — не соответствующим. Наличие неопределенности измерения влияет на процедуру контроля и делает необходимым установление баланса рисков производителя и потребителя.

8.3 Возможные значения контролируемой величины Y представляют распределением вероятностей. Можно рассчитать вероятность, с которой Y соответствует установленным требованиям при данном распределении вероятностей и заданных границах допустимых значений.

8.4 Из-за неполного знания величины Y (что отражает ее распределение вероятностей) существует риск ошибочного решения при определении соответствия установленным требованиям. Ошибочные решения бывают двух типов: когда значение величины признано соответствующим требованиям, но на самом деле им не является, и когда значение величины признано несоответствующим, но на самом деле установленным требованиям удовлетворяет. Связанные с этим риски относят, соответственно, к риску потребителя и риску производителя (см. [7]).

8.5 Риски ошибочного решения в части соответствия или несоответствия установленным требованиям можно уравновесить, выбирая интервал приемки для полученных значений измеряемой величины таким образом, чтобы минимизировать потери, связанные с ошибочными решениями [19]. В [7] рассматривается задача вычисления вероятности соответствия и вероятностей ошибочных решений указанных двух типов для заданных распределений вероятностей и заданных границ интервала приемки. Выбор границ интервала приемки зависит от последствий принятия ошибочных решений.

8.6 Хотя сказанное в 8.3 и 8.5 справедливо для любых распределений вероятностей, в [7], в основном, рассматривается случай нормального распределения как наиболее характерного для практики.

9 Применение метода наименьших квадратов

9.1 Руководство по применению метода наименьших квадратов (известного также как подгонка методом наименьших квадратов) для задач по оцениванию данных в метрологии представлено в [3]. В таких задачах часто используется некоторое теоретическое соотношение между независимой и зависимой переменными. Это соотношение составляет основу для подгонки кривой под имеющиеся данные посредством подбора параметров теоретической зависимости. Входные величины в соответствующей модели измерений — это зависимые и независимые переменные, для которых получены данные измерений. Выходные величины — это искомые параметры зависимости. Способ, которым выходные величины получают из входных величин посредством метода наименьших квадратов, определяет модель измерения.

9.2 Применительно к калибровке (см. 6.8) значение измеряемой величины независимой переменной в большинстве случаев получают от эталона. Значение зависимой переменной будет показанием, полученным измерительной системой для соответствующего значения независимой переменной. Установленная в [3] процедура подгонки кривой, частным случаем которой является калибровочная характеристика, получаемая в процессе калибровки, является обобщением обычного метода наименьших квадратов.

9.3 Измерительная задача состоит в том, чтобы оценить параметры (а иногда и число этих параметров) поданным, представляющим собой набор пар из полученного значения измеряемой величины и соответствующего показания. Эти пары вместе с ассоциированными стандартными неопределенностями и, когда уместно, ковариациями, составляют исходные данные для процедуры подгонки.

9.4 Типичные измерительные задачи, к которым может быть применено руководство [3], включают в себя: а) подгонку линейных и нелинейных зависимостей, включая случай неточно известных значений независимой переменной; b) выбор модели из некоторого класса для оценки параметров физического процесса. Применение [3] не ограничено в самом строгом смысле задачами подгонки кривой. Это руководство может также быть использовано для обработки данных, например, в задачах свертки [21], согласования фундаментальных констант [22] и оценивания данных ключевых сличений [9].

9.5 Задачи, указанные в 9.4, перечисление а), предполагают, что после оценивания методом наименьших квадратов параметров калибровочной характеристики и ассоциированных с ними стандартных неопределенностей и ковариаций измерительная система будет далее использоваться для проведения измерения, в ходе которого оценки параметров калибровочной характеристики вместе со значением полученного показания используют для оценивания измеряемой величины. Стандартную неопределенность, ассоциированную со значением оценки измеряемой величины, вычисляют с использованием стандартных неопределенностей и ковариаций для параметров калибровочной характеристики и стандартной неопределенностью, ассоциированной с показанием измерительной системы.

9.6 В [3] особо подчеркивается, что постановку и решение измерительной задачи методом наименьших квадратов следует осуществлять с учетом структуры неопределенности, т.е. с учетом стандартных неопределенностей для зависимых и независимых переменных и ковариаций для пар этих переменных.

9.7 Задачи, указанные в 9.4, перечисления а) и b), редко предполагают подгонку к значениям только одной выходной величины. Чаще случается так, что выходных величин несколько, поэтому соответствующие математические выражения удобнее представлять в матричной форме. В [3] матричный формализм использован максимально широко, что облегчает программирование алгоритма вычислений и соответствует потребностям практики измерений (см. также 7.5).

Приложение А
(справочное)

Используемые сокращения

В таблице А.1 представлены сокращения, используемые в настоящем документе

Таблица А.1

Приложение ДА
(справочное)

Сведения о соответствии ссылочных международных стандартов и документов межгосударственным стандартам

Таблица ДА.1

_______________

В Российской Федерации данный международный стандарт был введен как Р 50.1.040-2002 «Статистические методы. Планирование экспериментов. Термины и определения».

Приложение ДБ
(справочное)

Дополнительные замечания к межгосударственным стандартам, вводящим международные руководства в области неопределенности измерения

ДБ.1 Общие замечания к серии межгосударственных стандартов ГОСТ 34.100

ДБ.1.1 Серия межгосударственных стандартов ГОСТ 34.100 вводит документы, разрабатываемые рабочей группой JCGM/WG 1 «Рабочая группа по выражению неопределенности измерения», входящей в состав объединенного комитета JCGM «Объединенный комитет по руководствам в метрологии» при Международном бюро мер и весов (см. «Предисловие к международному документу ISO/IEC Guide 98.1:2009» настоящего стандарта).

ДБ.1.2 Документы, разрабатываемые JCGM/WG 1, устанавливают общий единообразный подход к оценке точности измерений через концепцию неопределенности измерений и включают в себя как методы вычисления неопределенности измерения в разных измерительных задачах, так и учет неопределенности измерения при применении результатов измерения.

ДБ.1.3 Концепция неопределенности измерения разработана для выражения качества результата измерения взамен концепции погрешностей измерений с целью придания методической корректности используемым теоретико-вероятностным моделям.

В концепции погрешностей измерений результат измерения представляют в виде суммы истинного значения и погрешности, которая, в свою очередь, является суммой систематической и случайной составляющих. При этом для оценки точности измерения обычно используют один из двух способов: консервативный (оценка сверху) и теоретико-вероятностный. Выбор того или иного способа оценивания определяется конкретной измерительной задачей и дальнейшим использованием результата измерения. Каждый из этих подходов имеет ограничения в применении.

ДБ.1.4 При консервативном способе оценивания границы суммарной погрешности определяются арифметическим суммированием границ ее составляющих. Главный недостаток консервативного способа — слишком широкие границы суммарной погрешности, особенно в случае большого числа составляющих. Консервативный подход может найти применение в измерительных задачах, где необходимо обеспечить нахождение истинного значения измеряемой величины в установленных границах наверняка.

ДБ.1.5 При теоретико-вероятностном подходе для описания результата измерения используется случайная переменная, математическое ожидание которой совпадает с истинным значением измеряемой величины или смещено относительно него на величину систематической погрешности. Это дает возможность в условиях ограниченного числа повторных наблюдений измеряемой величины строить для нее точечные и интервальные оценки.

В теории погрешностей использована частотная интерпретация вероятности, наблюдения рассматриваются как выборка из заданной генеральной совокупности, оценки измеряемой величины и характеристик погрешности являются статистиками. В качестве интервальной оценки используется построенный на основе статистик доверительный интервал, соответствующий заданной доверительной вероятности.

Главным ограничением использования частотного подхода является невозможность его корректного распространения на задачу оценивания систематических погрешностей. Подход, основанный на «рандомизации» систематических погрешностей, применим лишь в отдельных случаях. В результате в рамках частотного подхода невозможно указать в общем виде правило построения доверительного интервала погрешности, особенно при наличии нескольких влияющих факторов, каждый из которых может описываться своей генеральной совокупностью и для которых могут быть получены свои выборки наблюдений. При отсутствии строгих математических методов метрологам часто приходилось обращаться к инженерным (эмпирическим) процедурам определения доверительных интервалов без оценки качества получаемых результатов.

_______________

Примером такой инженерной процедуры является способ оценивания доверительных границ погрешности в РМГ 43-2001 «Применение «Руководства по выражению неопределенности измерений».

ДБ.1.6 Введение в метрологическую практику концепции неопределенности измерения «Руководством по выражению неопределенности измерения (GUM)», опубликованным в 1993 г. (см. «Предисловие к международному документу ISO/IЕС Guide 98.1:2009» настоящего стандарта), явилось попыткой дать математически строгий единый подход к оценке составляющих неопределенности, обусловленных как случайными, так и систематическими факторами, при заданных условиях измерительной задачи. Однако GUM не смог в полной мере решить эту задачу, он появился как внутренне противоречивый документ, использующий одновременно частотную и байесовскую концепции вероятности. Единая процедура вывода, наиболее корректно и последовательно описанная в JCGM 101:2008, основана на отказе от частотной интерпретации вероятности при оценке точности измерения в пользу субъективного представления о вероятности. Если в частотном подходе понятие случайной переменной использовано для описания результата/погрешности измерения, то в субъективном подходе случайная переменная использована для описания возможных значений измеряемой величины. При этом получение распределения вероятностей, ассоциированного с измеряемой величиной, осуществляется на основе:

— составления для данной измерительной задачи модели измерений, связывающей измеряемую величину (выходную величину) со всеми значимыми влияющими величинами (входными величинами модели);

— приписывания входным величинам распределений вероятностей (в общем случае, совместных), исходя из имеющейся информации об этих величинах и их наблюдений (при наличии);

— преобразования совместного распределения входных величин в распределение выходной величины согласно правилам преобразования случайных переменных.

В отличие от теории погрешностей (на основе частотного подхода) концепция неопределенности (на основе субъективной вероятности) не имеет принципиальных ограничений в получении окончательного результата измерения в виде функции распределения, ассоциированной с измеряемой величиной, что позволяет вычислить интервал вероятности (охвата) для любой заданной вероятности. Однако во многих измерительных задачах аналитическое решение задачи преобразования плотностей вероятностей невозможно. В этом случае точное решение (в пределах точности вычислений) всегда может быть получено числовым методом Монте-Карло (см. JCGM 101:2008).

ДБ.1.7 При наличии выборки наблюдений одной или нескольких входных величин (например, показываемой величины — см. JCGM 104:2009, пункт 3.2) входное распределение для этой величины получают применением теоремы Байеса. Поэтому переход от концепции погрешностей к концепции неопределенности может рассматриваться как переход от частотного (объективного) подхода в интерпретации вероятностей к байесовскому (субъективному).

Примечание — Существует широкий круг измерительных задач, в которых получают только одно наблюдение для входной величины. Однако и в этом случае возможно формальное применение теоремы Байеса, поэтому концепцию неопределенности измерения можно связывать с байесовским подходом без потери общности.

ДБ.1.8 Важными характеристиками результатов измерений в обоих подходах являются интервальные оценки, которые, однако, имеют разное содержание. В частотном подходе это доверительный интервал, неявно предполагающий возможность проведения неограниченной серии измерений и гарантирующий накрытие истинного значения измеряемой величины в заданной доле р таких измерений. В байесовском подходе это интервал охвата, содержащий с вероятностью q значение измеряемой величины.

Примечание 1 — Часто, задавая р=q, пытаются провести количественное сопоставление получаемого доверительного интервала с интервалом охвата. Однако необходимо иметь в виду, что подобные попытки некорректны ввиду сопоставления разных величин.

Примечание 2 — Встречающееся в литературе утверждение, что оба подхода дают одинаковые интервальные оценки, несмотря на их разную интерпретацию, в общем случае неверно. Равенство оценок имеет место только в отдельных измерительных задачах, хотя к ним, например, относится часто встречающийся случай, когда можно обоснованно предположить наличие одной доминирующей влияющей величины, распределенной по нормальному закону. Для данной задачи, действительно, доверительный интервал (наименьший) совпадет с интервалом охвата (наименьшим), поскольку центральная статистика, используемая для построения доверительного интервала, подчиняется тому же t-распределению, которое после операций сдвига и масштабирования дает апостериорное распределение для измеряемой величины (при условии задания неинформативных априорных распределений для математического ожидания и дисперсии нормального распределения) в байесовском подходе.

ДБ.1.9 Разница между частотным и байесовским подходами наглядно проявляется в том, насколько в рамках данного подхода легко получить ту или иную характеристику результата измерения. Частотный подход основан на получении оптимальных точечных оценок (статистик), по которым потом можно построить (не всегда) доверительный интервал. Распределение случайной погрешности, характеризующей качество измерений, может быть получено только в отдельных частных случаях. В байесовском подходе ситуация противоположная. В первую очередь, получают распределение вероятностей случайной величины, ассоциированной с измеряемой величиной. На его основе всегда есть возможность построить интервал охвата. Точечную оценку получают из распределения вероятностей после принятия каких-либо дополнительных допущений.

Примечание — В зависимости от целей измерений точечной оценкой могут служить разные параметры полученного распределения для измеряемой величины, такие как математическое ожидание, медиана или мода.

ДБ.1.10 Достоинством байесовского подхода, а значит и концепции неопределенности измерений, является наличие формализованной процедуры учета априорной информации разного рода (в том числе, о возможных или наиболее вероятных значениях измеряемой величины) при получении результата измерений.

Сопоставление концепций погрешности и неопределенности измерения проиллюстрировано на рисунке ДБ.1.

Примечание — Вопросительные знаки на схеме частотного подхода показывают, что получение оценки данной характеристики затруднено или невозможно. Если особенности измерительной задачи позволяют получить распределение погрешности, то доверительный интервал может быть рассчитан. Обратное утверждение, вообще говоря, неверно.

Рисунок ДБ.1 — Обобщенная схема получения результата измерения в рамках частотного и байесовского подходов

ДБ.1.11 В рамках байесовского подхода решением измерительной задачи в общем случае является распределение, ассоциированное с измеряемой величиной, которое, в общем случае, индивидуально для каждой измерительной задачи и в наиболее полном виде описывает всю собранную при решении данной задачи информацию.

В целях сокращения объема передаваемых данных и удобства их хранения в документах, разрабатываемых JCGM/WG 1, основным способом представления результата измерения принят интервал охвата (или область охвата в случае многомерной измеряемой величины). При этом, однако, следует помнить, что за областью охвата всегда стоит распределение соответствующей случайной переменной и, главное, во многих практических приложениях результатов проведенного измерения необходимо знать не интервал охвата, а распределение, из которого оно получено. Поэтому, как правило, желательно сохранять результат измерения в виде распределения вероятностей случайной переменной, ассоциированной с измеряемой величиной.

ДБ.2 Дополнительные замечания к настоящему стандарту

ДБ.2.1 Настоящий стандарт является введением международного документа JCGM 104:2009, который, в свою очередь, основан на плане работ JCGM/WG 1. Следует и меть в виду, что хронологически документы, заявленные к разработке, разрабатываются и публикуются не в соответствии с их обозначениями. Так, на настоящий момент опубликовано пять таких документов: JCGM 100:2008, JCGM 101:2008, JCGM 102:2011, JCGM 104:2009 (настоящий стандарт) и JCGM 106:2012. Остальные запланированные документы еще не начали разрабатываться или их разработка находится на ранней стадии.

ДБ.2.2 В разделе 8 рассматривается документ JCGM 106:2012, который должен установить общий подход к использованию неопределенности измерения при оценке соответствия измеряемой величины установленному критерию. Особенностью концепции неопределенности измерения является использование теоретико-вероятностных методов на основе субъективной (байесовской) интерпретации вероятности в отличие от теории погрешностей измерений, где при описании результата измерения используется частотная интерпретация вероятности. Выбранный байесовский подход при вычислении неопределенности измерения позволил обосновать единые правила вычисления и суммирования составляющих неопределенности. Поэтому необходимо обратить внимание, что в документе JCGM 106:2012 наблюдается отход от методологии неопределенности измерения, решение о соответствии принимается не на основании распределения вероятностей, ассоциированного с измеряемой величиной, а на основе распределения показаний средства измерений, которое имеет частотную интерпретацию. Таким образом, результат измерения рассматривается не как индивидуализированное распределение, полученное для конкретного измерения (см. ДБ.1.11), а как выборочное значение из заданной генеральной совокупности, характеризующей измерительную систему (см. JCGM 106:2012, п.6.2.4). Такая постановка задачи оценки соответствия может рассматриваться только как частный случай, а указанная интерпретация вероятностей соответствует не байесовскому, а частотному подходу.

Библиография

Электронный текст документа
и сверен по:

, 2018