ОЦЕНИВАНИЕ НЕОПРЕДЕЛЕННОСТИ ИЗМЕРЕНИЙ МЕТОДОМ МОНТЕ-КАРЛО - Студенческий научный форум

VIII Международная студенческая научная конференция Студенческий научный форум - 2016

ОЦЕНИВАНИЕ НЕОПРЕДЕЛЕННОСТИ ИЗМЕРЕНИЙ МЕТОДОМ МОНТЕ-КАРЛО

Жмылева (Николаева) А.Н. 1, Сибринин Б.П. 1
1Пензенский Государственный Университет
 Комментарии
Текст работы размещён без изображений и формул.
Полная версия работы доступна во вкладке "Файлы работы" в формате PDF
Содержание

Введение

3

1 Понятие «неопределенность измерений» и методы ее оценивания

5

  1.  
    1. Понятие «неопределенность измерений» и его связь с термином «погрешность измерений»

5

  1.  
    1. Требования к точности измерений и описание методик измерений

19

  1.  
    1. Основные этапы расчета погрешности и неопределенности измерений

22

  1.  
    1. Методы оценивания неопределенности измерений

24

   
  1. Анализ оценивания неопределенности измерений методом Монте-Карло и сравнение с методом оценивания неопределенности GUM

32

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

32

2.1.1 Возникновение метода Монте-Карло и его общая характеристика

32

2.1.2 Применение метода Монте-Карлодля оценивания неопределенности измерений

33

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

34

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

35

2.1.5 Интервал охвата для выходной величины

36

2.1.6 Адаптивная реализация метода Монте-Карло

37

2.1.7 Преимущества и недостатки метода Монте-Карло

39

  1.  
    1. Сравнение метода оценивания неопределенности GUM и метода Монте-Карло

42

   
  1. Сравнение результатов оценивания неопределенности измерений методами GUM и Монте-Карло

45

   

Заключение

54

Список использованных источников

55

   
   
   

Введение

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

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

В большинстве случаях оценивания неопределенности измерений метод Монте-Карло используется, когда:

- вклад разных составляющих неопределенности существенно неодинаков;

- распределение выходной величины не является нормальным или масштабированным смещенным t - распределением;

- модель достаточно сложная;

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

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

Цель данной работы состоит в изучении методов оценивания неопределенности измерений и применении метода Монте-Карло для оценивания неопределенности измерений.

Задачи работы:

- изучение различий между понятиями «погрешность измерений» и «неопределенность измерений»,

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

- сравнение результатов оценки неопределенности измерений методом Монте-Карло и расчетным методом, изложенным в ИСО/МЭК «Руководство по выражению неопределенности измерений (GUM)».

Автор аттестационной работы выражает глубокую благодарность руководителю – старшему преподавателю кафедры МСК ПГУ Сибринину Борису Петровичу за проявленное внимание и участие в решении проблем, которые возникали в процессе выполнения аттестационной работы.

1 ПОНЯТИЕ «НЕОПРЕДЕЛЕННОСТЬ ИЗМЕРЕНИЙ» И МЕТОДЫ ЕЕ ОЦЕНИВАНИЯ

1.1 Понятие «неопределенность измерений» и его связь с термином «погрешность измерений»

Считается, что термин «неопределенность» пришел на замену термину «погрешность». Однако, это не так. Термин «погрешность» остается существовать, и он входит в международный метрологический словарь VIM [12].

В 1993 г. был выпущен документ ИСО/МЭК «Руководство по выражению неопределенности измерений (GUM)» [15], в котором было введено понятие «неопределенность измерений».

В Руководстве большое внимание уделяется различию терминов «погрешность» и «неопределенность». Утверждается, что они не синонимы и представляют собой совершенно различные понятия. По определению авторов Руководства, погрешность измерения – отклонение результата измерения от истинного значения измерения. Поскольку истинное значение не может быть определено, то и погрешность, как идеализированное понятие не может быть известна точно. Поэтому, авторы постарались во - первых, по возможности отказаться, от использования при изложении понятий «погрешность» и «истинное значение измеряемой величины» в пользу понятий «неопределенность» и «оцененное значение измеряемой величины», во - вторых, переход от деления погрешностей по природе их проявления на «случайные» и «систематические» к другому делению – по способу оценивания неопределенностей измерений (по типу А – методами математической статистики и по типу В – другими методами).

Основным понятием, используемым в Руководстве, является понятие «неопределенность измерения». В разделе II дается два определения этого понятия:

- «в своем самом широком смысле «неопределенность измерения» означает сомнение относительно достоверности результата измерения»;

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

Обратимся к формулам, описывающим термины «неопределенность измерения» и «погрешность измерения».

При определении погрешности измерения, обозначаемой Δ, исходной является хорошо известная всем формула:

Δ= хизм - хист (1),

где хизми хист ‑ результат измерения и истинное значение измеряемой величины, соответственно.

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

Δхизм– хдейст (2)

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

Вопрос о том, насколько близким – один из наиболее трудных вопросов втеории погрешностей измерений.

Следует заметить, что в формуле (2) вместо знака строгого равенства появляется знак приближенного равенства «». С этого момента мы вынуждены рассматривать погрешность измерений (Δ), как сугубо неопределенную величину, для которой можно получать лишь более или менее приемлемые (по точности) оценки.

Основным количественным выражением неопределенности измерений является стандартная неопределенность.

В Руководстве не оперируют формулами вида (1) и (2), а измеряемая величина Y определяется как функция

Y = ƒ(X1, X2, …, XN), (3)

где X1, X2, …, XN – входные величины (непосредственно измеряемые или другие величины, влияющие на результат измерения);

N– число этих величин;

ƒ – вид функциональной зависимости.

Оценка измеряемой величины y вычисляется как функция оценок входных величин x1, x2, …, xN после внесения поправок на все известные источники, имеющие систематический характер:

y = f (x1, x2, …, xN) (4)

Затем вычисляются стандартные неопределенности входных величин u(xi) (i=1, 2,…, N).

Различают два типа вычисления стандартной неопределенности:

- вычисление по типу А – как среднеквадратическое (стандартное) отклонение закона распределения результата измерения с классической частотной интерпретацией;

- вычисление по типуВ – с использованием других способов.

Исходными данными для оценки стандартной неопределенности по типу А являются результаты многократных измерений: xi, ( i=1, 2,…, n)

В большинстве случаев наилучшей оценкой ожидаемого значения μx величины x, изменяющийся случайным образом, для которой получены n независимых наблюдений xi при одинаковых условиях измерения, является среднее арифметическое или среднее значение из nнаблюдений:

Экспериментальная дисперсия наблюдений, которая оценивает дисперсию σ2 распределения вероятностей x, определяется по формуле

Эта оценка дисперсии выборки и ее положительный квадратный корень sназываемый экспериментальным стандартным отклонением, характеризует изменчивость наблюдаемых значений xi или, точнее, их дисперсию относительно среднего значения . Наилучшая оценка σ2() дисперсии среднего значения определяется как:

(6).

Экспериментальная дисперсия среднего и экспериментальное стандартное отклонение среднего значения s(x), равное положительному квадратному корню из , количественно определяют, насколько хорошо оценивает математическое ожидание μx величины x. Таким образом, для входной величины Xi, определенной из nнезависимых наблюдений, стандартная неопределенность u(xi) ее оценки xi= есть u(xi) c , вычисленным согласно уравнению (5), т.е.

Исходными данными для оценки неопределенности по типу В являются:

- данные предшествовавших измерений величин, входящих в уравнение измерения;

- сведения о виде распределения вероятностей;

- данные, основанные на опыте исследователя или общих знаниях о поведении и свойствах соответствующих приборов и материалов;

- неопределенности констант и справочных данных;

- данные поверки, калибровки, сведения изготовителя о приборе и др.

Неопределенности этих данных обычно представляются в виде границ отклонения значения величины от ее оценки. Наиболее распространенный способ формализации неполного знания об измеряемой величине базируется на постулате равномерного закона распределения в указанных (нижней и верхней) границах (a-, a+). При этом стандартная неопределенность, оцененная по типуВ, равна:

Для симметричных границ (±a):

u = uB.

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

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

Суммарная стандартная неопределенность uс(y) – стандартная неопределенность результата измерений, полученного через значения других величин, равная положительному квадратному корню суммы членов, причем члены являются дисперсиями или ковариациями этих других величин, взвешенными в соответствии с тем, как результат измерений изменяется при изменении этих величин [2]. Суммарная дисперсия рассчитывается по формуле

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

В случае некоррелированных оценок (x1, x2, …, xN) суммарная стандартная неопределенность вычисляется по формуле

В случае коррелированных оценок x1 , x2 , …, xNuc(y) вычисляется по формуле

где ;

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

Следует отметить, что коэффициент корреляции является числом в диапазоне

-1 ≤ ≤ +1.

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

U =k×uc(y), (12)

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

Коэффициент охвата определяют по формуле:

k = t(νэфф),

где t(νэфф)– коэффициент Стьюдента при числе эффективных степеней свободы νэфф.

Значение νэфф находят по формуле Велча-Статтертвейта:

где νi= n-1 число степеней свободы для оценивания неопределенности по типу А,

n– число результатов измерений;

νi = ∞ для оценивания по типу В.

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

Y = y ± U = y±k×uc(y),

соответствующий выбранному уровню доверия, такому как 0,95 или 0,99. Равным образом, для заданного значения k хотелось бы иметь возможность четко указать уровень доверия, связанный с этим интервалом. Однако это нелегко осуществить на практике, поскольку требуется полное знание закона распределения вероятностей, характеризуемого результатом измерения y и его суммарной стандартной неопределенностью uc(y). Во многих практических случаях при вычислении неопределенностей измерений делают предположение о нормальности закона распределения возможных значений измеряемой величины и полагают: k =2 при уровне доверия p≈ 0,95 или k =3 при уровне доверия p≈ 0,99[1].

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

Составляющие неопределенности (погрешности) измерений рассмотрены в [6] и включают:

а) Методические составляющие:

1) составляющие, обусловленные неадекватностью выбранной модели объекта измерений его свойствам;

2) составляющие, обусловленные отклонением от номинальных значений параметров функции, связывающей измеряемую величину с величиной на входе средства измерений;

3) составляющие, обусловленные квантованием по уровню (при использовании средств измерений с аналого-цифровым преобразованием);

4) составляющие, обусловленные вычислительными алгоритмами;

б) Инструментальные составляющие:

1) основная погрешность средства измерений;

2) дополнительные погрешности средства измерений;

3) составляющая, обусловленная вариацией (гистерезисом) средства измерений;

4) составляющая, обусловленная взаимодействием средства измерений с объектом измерений;

5) динамическая составляющая, обусловленная инерционностью средства измерений;

6) составляющие, связанные с отбором и приготовлением проб веществ.

в) Составляющие, обусловленные действиями оператора (субъективные составляющие):

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

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

Результаты сопоставления количественных оценок погрешности измерений и неопределенности измерений представлены в таблице 1.

Сравнение способов оценки расширенной неопределенности и доверительных границ суммарной погрешности показывает, что:

- для большинства случаев получаемые оценки отличаются незначительно;

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

Таблица 1 – Сопоставление количественных оценок погрешности измерений и неопределенности измерений

Погрешность измерения

Неопределенность измерения

СКО, характеризующее случайную погрешность

где

Стандартная неопределенность, оцененная по типу А

где

СКО, характеризующее неисключенную систематическую погрешность

где k=1,1 при P=0,95 и

k=1,4при P=0,99 n>4

Стандартная неопределенность, оцененная по типу В

В случае неизвестного закона распределения вероятностей наибо-лее распространенный способ формализации неполного знания об измеряемой величине базируется на постулате вероятностного закона распределения (обычно равномерно-го) в указанных (нижней и верхней) границах (а+, а-). При этом стандартная неопределенность, оцененная по типу В, равна uB=

Для симметричных границ неисключенных систематических погрешностей (±а) uB=

СКО, характеризующее суммарную погрешность

Суммарная неопределенность

Окончание таблицы 1

Погрешность измерения

Неопределенность измерения

Доверительные границы погрешности

Упрощенный вариант

где t – квантиль распределения Стьюдента

Расширенная неопределенность

UР=kuc,

где k – коэффициент охвата

(k ≈2 на уровне доверия P=0,95;

k ≈3при уровне доверия P=0,99)

Интерпретация результата

Интервал (-ΔР; + ΔР) с вероятностью содержит погрешность измерений, что равносильно тому, что интервал (y-ΔР;y + ΔР)содержит истинное значение измеряемой величины

Интерпретация результата

Интервал (y-UР; y + UР) содержит большую долю (Р) распределения значений, которые могли бы быть обоснованно приписаны измеряемой величине

Примечание

1. Полная формула для вычисления доверительных границ погрешности

2. Полный вариант выбора коэффициента охвата k=tP(vэфф),

где эффективное количество степеней свободы

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

Указанные границы (отложенные от результата измерения) накрывают истинное значение измеряемой величины с заданной доверительной вероятностью (частотная интерпретация вероятности). В то же время расширенная неопределенность (аналогичный интервал) трактуется в Руководстве как интервал, содержащий заданную долю распределения значений, которые могли бы быть обоснованно приписаны измеряемой величине (субъективная интерпретация вероятности) [4].

Различие между традиционным подходом, использующим понятие «погрешность измерения», и подходом, изложенным в Руководстве, сводится к различию систем координат, относительно которых рассматривают значение измеряемой величины и результат измерений. При рассмотрении погрешности измерений начало системы координат совмещают со значением измеряемой величины, наблюдая рассеяние результата измерений; при рассмотрении неопределенности измерений – с результатом измерений, что и создает эффект рассеяния единственного значения измеряемой величины. Поскольку конкретному результату измерений соответствует совокупность возможных значений измеряемой величины, каждое из которых удовлетворяет условию: сумма возможного значения измеряемой величины и соответствующей ему реализации погрешности измерений должна быть равна результату измерений, закон распределения вероятностей возможных значений измеряемой величины определяется законом распределения вероятностей погрешности измерений.

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

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

На рисунке 1 изображена плотность вероятностей φ(aP), характеризующая распределение результата измерений aP, которую можно одновременно рассматривать и как плотность вероятностей φ(ΔИ) абсолютной погрешности измерений ΔИ, если перенести начало координат в точку, соответствующую истинному значению измеряемой величины Аист.

Рисунок 1 ‑ Плотность вероятностей φ(aP), характеризующая

распределение результата измерений aP

Представим наблюдателя, который всякий раз будет располагаться в точке числовой оси, соответствующей очередной реализации Аpi результата измерений, и наделим его способностью наблюдать на числовой оси истинное значение измеряемой величины Аист. В точке Ap1 он будет видеть Аист слева, в точке Ap2 ‑ справа, в точке Ap3 ‑ близко от себя, в точке Ap4 ‑ далеко. Он будет наблюдать Аист как случайную величину, закон распределения которой определяется законом распределения погрешности измерений.

Для более четкого понимания поместим наблюдателя в своего рода «инерциальную» систему координат, жестко привязанную к результату измерений, в которой наблюдатель остается «неподвижным». Чтобы подчеркнуть, что размер рассматриваемой измеряемой величины неизменен (фиксирован), наделим наблюдателя способностью видеть истинное значение измеряемой величины вместе с числовой осью истинных значений, которая в системе координат наблюдателя будет занимать, как это мы выяснили из рисунка 1, случайное положение. Тогда рисунок 1 трансформируется в рисунок 2, где Аистi – «реализации» истинного значения измеряемой величины в системе координат наблюдателя.

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

Для симметричных законов распределения свойство «зеркальности» не имеет значения. Иллюстрация «зеркальности» для несимметричного закона распределения приведена на рисунке 3 [8].

Рисунок 2 ‑ Плотность вероятностей φ(aист), характеризующая

распределение результата измерений aист

Рисунок 3 ‑ Несимметричный закон распределения

1.2 Требования к точности измерений и описание методик измерений

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

При описании требований к выражению погрешности и неопределенности измерений, выполненных с использованиемтеории шкал, применяют положения рекомендаций по межгосударственной стандартизации РМГ 83-2007 «ГСИ. Шкалы измерений. Термины и определения» с учетом особенностей конкретных шкал измерений.

Методы и средства измерений выбирают в соответствии с документами, относящимися к выбору методов и средств измерений данного вида, а при отсутствии таких документов – в соответствии с общими рекомендациями по метрологии МИ 1967-89 «ГСИ. Выбор методов и средств измерений при разработке методик выполнения измерений. Общие положения».

Если методика измерений предназначена для использования всфере государственного регулирования обеспечения единства измерений, то применяемые средства измерений и стандартные образцы должны быть утвержденных типов, внесенные в Государственный реестр СИ (стандартных образцов), поверенные, с не истекшим межповерочным интервалом; испытательное оборудование должно быть аттестовано.

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

Если полученное значение погрешности измерений выходит за заданные пределы, то погрешность измерений может быть уменьшена в соответствии с рекомендациями по межгосударственной стандартизации РМГ62-2003 «ГСИ. Обеспечение эффективности измерений при управлении технологическими процессами. Оценивание погрешности измерений при ограниченной исходной информации».

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

Планирование экспериментов по оценке характеристик погрешности методик измерений состава и свойств веществ и материалов и выбор способов экспериментальной оценки этих характеристик проводят в соответствии с ГОСТ Р ИСО 5725-1 ‑ ГОСТ Р ИСО 5725-6, неопределенности ‑ в соответствии с руководством ЕВРАХИМ/СИТАК «Количественное описание неопределенности в аналитических измерениях».

В документе, регламентирующем методику измерений, указывают [5]:

- наименование методики измерений;

- назначение методики измерений;

- область применения;

- условия выполнения измерений;

- метод (методы) измерений;

- допускаемую и (или) приписанную неопределенность измерений или норму погрешности и (или) приписанные характеристики погрешности измерений;

- применяемые средства измерений, стандартные образцы, их метрологические характеристики и сведения об утверждении ихтипов;

- операции при подготовке к выполнению измерений, в том числе по отбору проб;

- операции при выполнении измерений;

- операции обработки результатов измерений;

- требования к оформлению результатов измерений;

- процедуры и периодичность контроля точности получаемых результатов измерений;

- требования к квалификации операторов;

- требования к обеспечению безопасности выполняемых работ;

- требования к обеспечению экологической безопасности;

- другие требования и операции (при необходимости).

1.3 Основные этапы расчета погрешности и неопределенности измерений

Порядок расчета неопределенности измерений сводится к следующим этапам:

а) составление уравнения измерения, характеризующего измерительный процесс;

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

в) исследование возможности оценивания каждого источника по типу А или по типу В;

г) определение предполагаемого закона распределения для неопределенностей типа В;

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

е) расчет суммарной и расширенной неопределенности измерения.

Этапы расчета погрешности измерений:

а) составление уравнения измерения, характеризующего измерительный процесс;

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

в) формирование исходных данных для расчета погрешности;

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

д) объединение (суммирование) составляющих погрешности измерений, получение результирующей (суммарной) характеристики качества измерений [6].

1.4 Методы оценивания неопределенности измерений

Существует несколько методов оценивания неопределенности измерений, такие как метод, изложенный в ИСО/МЭК «Руководство по выражению неопределенности измерений (GUM)»; экспериментальный метод; метод Монте-Карло (подробнее данный метод рассмотрим в п. 2) и другие.

Процедура оценивания неопределенности измерений по GUM состоит из следующих этапов:

а) Выражают связь между измеряемой величиной Y и входными величинами Xi, от которых она зависит, в виде функциональной зависимости . Функция f должна содержать все величины, включая поправки и поправочные коэффициенты, которые могут существенно повлиять на неопределенность результата измерения;

б) Получают оценку xiвходной величины Xiлибо на основе статистического анализа ряда наблюдений, либо другими способами (к ним относятся величины, связанные с аттестованными эталонами, стандартными образцами веществ или материалов, и величины, значения которых указаны в справочниках);

в) Оценивают стандартную неопределенность u(xi) каждой входной оценки xi. При оценивании стандартной неопределенности по типу А используют формулу (7), по типу В – формулу (8).

г) Если среди входных величин есть коррелированные между собой, то оценивают их ковариации согласно формуле (11);

д) Рассчитывают результат измерения, т.е. находят оценку измеряемой величины у по функциональной зависимости f, используя в качестве аргументов Xi оценки xi, полученный на этапе б;

е) Определяют суммарную стандартную неопределенность uс(y) результата измерения у по стандартным неопределенностям и ковариациям входных оценок по формулам (9) и (10);

ж) Если требуется знать расширенную неопределенность U для определения интервала от у-U до у+U, в пределах которого, предположительно, находится большая часть распределения значений, которые можно с достаточным основанием приписать измеряемой величине Y, то суммарную стандартную неопределенность uс(y) умножают на коэффициент охвата k, обычно принимающий значения в диапазоне от 2 до 3, чтобы получить значение U по формуле (12);

и) Представляют результат измерения у вместе с его суммарной стандартной неопределенностью uс(y)или расширенной неопределенностью U.

Если мерой неопределенности результата измерения является суммарная стандартная неопределенность uс(y), то при представлении результата измерения необходимо:

- дать подробное определение измеряемой величины Y,

- привести оценку у измеряемой величины Y и суммарной стандартной неопределенности uс(y) с указанием единиц измерений,

- при необходимости указать относительную суммарную стандартную неопределенность uс(y)/|y|, у≠0.

Если мерой неопределенности результата измерения является расширенная неопределенность U, то при представлении результата измерения необходимо:

- дать подробное определение измеряемой величины Y,

- указать результат измерения в виде Y=yU с указанием единиц измерений для у и U,

- при необходимости указать относительную расширенную неопределенность U/|y|, у≠0,

- указать использованное для получения расширенной неопределенности значение k,

- указать приблизительный уровень доверия для интервала yU и пояснить, как он был определен.

Схематично последовательность оценивания неопределенности измерений по GUM представлена на рисунке 4 [2].

Рисунок 4 - Схема последовательности вычисления неопределенности

измерений

Экспериментальный метод оценивания неопределенности разделяется на два подхода:

а) использование данных внутрилабораторных исследований по разработке и оценки пригодности метода;

б) использование данных межлабораторных исследований.

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

- наиболее достоверную из имеющихся оценок общей прецизионности;

- наиболее достоверную имеющуюся оценку правильности и неопределенности (смещение нужно для введения поправки, неопределенность для оценки неопределенности);

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

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

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

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

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

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

Существует два случая прецизионности: повторяемость и воспроизводимость.

Прецизионность в условиях повторяемости:

- одна и та же методика,

- идентичные пробы,

- одинаковые условия,

- параллельные измерения (одновременно).

Прецизионнность в условиях воспроизводимости:

- одна и та же методика,

- разные условия (разные аналитики, партии реактивов одного типа, наборы мерной посуды, экземпляры средств измерений одного типа, лаборатории),

- разное время.

Количественными оценками прецизионности являются:

а) стандартное отклонение повторяемости – СКО результатов испытаний в одной лаборатории, полученных по методике в условиях повторяемости:

где - выборочное СКО результатов единичного анализа, полученных вi-ой лаборатории, которое определяется по формуле:

На основе расчета значений в m-ом образце проверяют гипотезу о равенстве генеральных дисперсий, используя критерий Кохрена:

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

для P=0,95.

в) стандартное отклонение воспроизводимости – СКО результатов анализа, полученных в условиях воспроизводимости:

где Xm – среднее арифметическое результатов анализа, полученных в L лабораториях, которое определяется по формуле:

г) предел воспроизводимости - допускаемое для принятой вероятности 95% абсолютное расхождение между наибольшими наименьшим из n результатов единичного анализа, полученных в условиях воспроизводимости:

где коэффициент k=1,2;…; 2,0, учитывающий условия проведения эксперимента, для P=0,95 .

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

а) Рассчитать оценку математического ожидания систематической погрешности методики анализа:

где - среднее значение результатов анализа; ‑ аттестованное значение m-го образца.

б) Проверить значимость вычисленных значений по критерию Стьюдента, вычисляя

где ‑ погрешность аттестованного значения m-го образца.

Значение определяется по формуле

Если условие выполняется, то Θпринимается равным нулю. Если условие не выполняется, то вводят дополнительное условие:

|Θ| ≤ ξσRm,

где ξ = 0,5; …; 1, в зависимости от условий, средств измерений и характеристик объекта анализа.

Стандартная неопределенность с использованием показателей правильности и прецизионности вычисляется согласно ISO 21748:2010 «Руководство по использованию оценок повторяемости, воспроизводимости и правильности при оценке неопределенности измерений».

Основная модель [9]:

где m – общее среднее;

В – лабораторная составляющая систематического сдвига в условиях повторяемости с нулевым средним значением и стандартным отклонением σL;

е – случайная погрешность лаборатории, распределенная нормально с нулевым средним значением и стандартным отклонением σW.

Исходя из этого выражения, неопределенность y будет равна:

где - оценка межлабораторного стандартного отклонения, равная значению ;

u(e) – оценка стандартного отклонения повторяемости, равная значению Sr, которое определяется по формуле

Так как – оценка дисперсии воспроизводимости, то

Если в модель включаются данные правильности, то она принимает следующий вид:

где ‑ опорное значение;

δ – систематический сдвиг метода.

Отсюда:

где - неопределенность, соответствующая сертифицированному значению , используемому для оценки правильности при совместном исследовании метода;

- оценка стандартного отклонения смещения, причем

где L – количество лабораторий;

n – количество повторений в каждой лаборатории.

Если y=f(x1, x2, … , xm)и существуют отклонения от номинальных значений xi, не учитываемые в процессе совместного исследования, то объединенная модель будет иметь следующий вид:

где - коэффициент чувствительности, определяемый по формуле

Неопределенность результата измерений в этом случае будет равна:

2 АНАЛИЗ ОЦЕНИВАНИЯ НЕОПРЕДЕЛЕННОСТИ ИЗМЕРЕНИЙ МЕТОДОМ МОНТЕ-КАРЛО И СРАВНЕНИЕС МЕТОДОМ ОЦЕНИВАНИЯ НЕОПРЕДЕЛЕННОСТИ GUM

2.1 Метод Монте-Карло

2.1.1 Возникновение метода Монте-Карло и его общая характеристика

В 1930-х годах Энрико Ферми (Италия), а затем и Джон фон Нейман и Станислав Улам в 1940-х в Лос-Аламосе (США) предположили, что связь между стохастическими (случайными) процессами и дифференциальными уравнениями можно использовать «в обратную сторону». Они предложили воспользоваться стохастическим подходом для аппроксимации многомерных интегралов в уравнениях переноса, возникших в связи с задачей движения нейтрона в изотропной среде.

Идея была развита Станиславом Уламом, который раскладывал пасьянс и заинтересовался вопросом: «Какова вероятность того, что пасьянс сложится?» Вместо того, что бы использовать обычные комбинации, он предположил, что можно поставить эксперимент множество раз и оценить вероятность, подсчитав число удачных попыток.

Годом рождения метода Монте-Карло является 1949 год, когда была опубликована статья Николаса Метрополиса и Станислава Улама «Метод Монте-Карло». Название метода происходит от названия города в княжестве Монако, широко известного своими многочисленными казино, ведь именно рулетка является одним из самых известных генераторов случайных чисел. Станислав Улам пишет в своей автобиографии «Приключения математика», что название было предложено Николасом Метрополисом в честь его дяди, который был азартным игроком.

При расчете неопределенности измерений метод Монте-Карло (ММК) используется как метод трансформирования распределений на основе моделирования случайных выборок из этих распределений [3].

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

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

В большинстве случаях оценивания неопределенности измерений ММК используется, когда:

- вклад разных составляющих неопределенности существенно неодинаков;

- распределение выходной величины не является нормальным или масштабированным смещенным t - распределением;

- трудно найти частные производные от функции измерения, как того требует закон трансформирования неопределенностей;

- модель достаточно сложная;

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

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

2.1.2 Применение метода Монте-Карло

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

Реализация метода Монте-Карло представлена на рисунке 5.

Рисунок 5 - Реализация ММК

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

Необходимо выбирать M ‑ количество оцениваний модели, которое необходимо произвести. Лучше всего выбирать достаточно большие значение М (например, превышающие в 104 раз) по сравнению с 1/(1-p). Тогда можно ожидать, что функция распределения G обеспечит приемлемое дискретное представление GY(η) (функции распределения для Y) вблизи границ 100 p %-ного интервала охвата для Y.

Для применения метода формируют М векторов xr,r=1,…, M в соответствии с плотностями распределения вероятностей gXii) для Nвходных величин Xi или, если это необходимо, из совместной плотности распределения gX (о).

Модель оценивается для каждого из М извлечений из функции плотности вероятности (ФПВ) для N входных величин. Конкретнее, необходимо обозначить M извлечений через x1, …, xM, где r-е значение xr состоит из случайных значений x1,r, …, xN,r, а xi,r ‑ случайное значение из ФПВ для Xi. Тогда значения модели можно представить в виде:

yr = f(xr), r = 1, …, M.

Дискретное представление G функции распределения GY(η) для выходной величины Y может быть получено следующим образом:

а) необходимо рассортировать значения модели yr, r = 1, …, M, полученные по ММК, в неубывающем порядке. Обозначить рассортированные значения модели как y(r), r = 1, …, M;

б) если необходимо, создать возмущения для любых дублирующих значений модели y(r) так, чтобы конечный полный набор y(r), r = 1, …, M, формировал строго возрастающую последовательность;

в) полученная последовательность y(r), r = 1, …, M определяет G.

Если выходная величина Y будет рассматриваться как входная величина при оценивании неопределенности другого измерения, то выборку из ее распределения легко получить случайным (равновероятным) выбором значений из y(r), r = 1, …, M.

Последовательность y(r) (или yr) может быть скомпонована в гистограмму (при подходящей ширине ячеек) и представлять собой частотное распределение, которое, при условии, что оно нормировано, чтобы иметь единичную площадь, обеспечивает приближение к ФПВ gY(η) для Y. Вычисления характеристик распределения обычно проводятся не в терминах этой гистограммы, разрешение которой зависит от выбора ширины ячеек, а в терминах G. Тем не менее, гистограмма может быть полезна для понимания природы ФПВ, например степени ее асимметрии.

В ряде случаев полезна аппроксимация GY(η) непрерывной функцией.

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

В качестве оценки выходной величины Y используется выборочное среднее:

а в качестве оценки ее стандартной неопределенности u(y) – выборочное стандартное отклонение

При некоторых особых обстоятельствах, когда одной из входных величин приписано t-распределение с числом степеней свободы менее трех, математическое ожидание и стандартное отклонение Y, описываемой ФПВ gY(η), могут не существовать. Формулы (14) и (15) могут не обеспечить получение содержательных результатов. Однако интервал охвата для Y может быть сформирован, т. к. G имеет смысл и может быть определена.

2.1.5 Интервал охвата для выходной величины

Пусть q = pM, если pM ‑ целое число. В противном случае в качестве q можно выбрать целую часть (pM+1/2). Тогда [ylow, yhigh] ‑ это100p % интервал охвата для Y, где для любого r из ряда r = 1, …, (M – q),ylow= y(r) и yhigh= y(r+q). Вероятностно симметричный 100p % интервал охвата можно получить, взявr = (M − q)/2, если (M − q)/2 — целое число, r=int[(M− q + 1)/2] ‑ в противном случае.

Для нахождения наименьшего 100p % интервала охвата задается определение r* ,чтобы, что для r = 1, …, M – q выполнялось неравенство:

y(r*+ q) – y(r*) ≤ y(r+q) – y(r).

Из-за стохастичности в ММК некоторые из протяженностей этихM-q интервалов могут быть короче, чем они будут в среднем, а некоторые ‑ длиннее. Таким образом, при выборе такой наименьшей протяженности (аппроксимация) к наикратчайшему 100p % интервалу охвата имеет тенденцию быть, по крайней мере, меньше, чем при вычислении из GY(η), приводя к тому, что типичная вероятность охвата меньше, чем 100p %. Однако для больших М этим отличием можно пренебречь.

2.1.6 Адаптивная реализация метода Монте-Карло

Суть адаптивной процедуры состоит в последовательном увеличении числа испытаний до тех пор, пока полученные числовые оценки статистических характеристик не станут установившимися. Численный результат считается установившимся, если соответствующее ему удвоенное стандартное отклонение станет меньше заданной точности вычисления стандартной неопределенности u(y) [3].

Пусть ndig ‑ число десятичных знаков, рассматриваемых как значащие в числовом представлении величины z. Предел погрешности вычисления δ, в таком случае, определяется следующим образом:

а) выражают значениеz в виде c*10l, где с ‑ это целое от ndig десятичных знаков, а l ‑ целое число;

б) δ определяют по формуле:

δ (15)

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

а) необходимо задать ndig;

б) задают M = max(J, 104), где J ‑ наименьшее целое, которое больше или равно 100/(1 - p);

в) задают h = 1 (счетчик итераций ММК);

г) проводят M испытаний методом Монте-Карло;

д) используют M полученных на выходе модели значений y1, …, yM, для вычислений h-й оценки y(h)величины Y, ее стандартной неопределенности u[y(h)], левой и правой границ 100p % интервала охвата;

е) еслиh = 1, то увеличивают h на единицу и возвращаются к шагу г);

ж) вычисляют выборочное стандартное отклонение sy, связанное со средним оценок y(1), …, y(h) величины Y, задающееся как:

и) аналогичным образом вычисляют выборочное стандартное отклонение для средних значений оценок u(y),,;

к) используют далее все hM доступные значения модели для вычисления u(y);

л) определяют предел погрешности δ, связанной с u(y);

м) если хотя бы одно из значений 2sy, 2su(y),или превышает δ, необходимо увеличить h на единицу и вернуться к шагу г);

н) если возврата к этапу г) не произошло, и значения всех вычисляемых значений можно считать установившимися, то на основе hM полученных значений выходной величины вычисляют y, u(y) и 100p% интервал охвата.

Обычно на этапе а) выбирается ndig=1 или ndig=2.

На этапе ж) y может рассматриваться как реализация случайной величины со стандартным отклонением sy.

В ситуациях, где интервал охвата не требуется, проверка стабилизации вычислений на этапе м) может быть основана вместо этого только на 2sy и 2su(y).

Альтернативный неадаптивный подход для симметричного интервала охвата при вероятности 95 %, основанный на использовании статистик биноминального распределения, состоит в следующем. Выбирают М = 105 или М= 106. Формируется интервал [y(r), y(s)], где для М = 105, r= 2420 и s=97581, а для М = 106, r= 24747 и s= 975254. Этот интервал ‑ 95 % статистический интервал охвата при уровне доверия 0,99, т. е. вероятность охвата будет не менее 95 % в, как минимум, 99 % применений ММК. Средняя вероятность охвата такого интервала будет (s− r)/(M + 1), что превышает 95 % на величину, которая становится меньше при возрастании M, а именно, 95,16 % для М = 105 и 95,05 % для М = 106.

Эти результаты могут быть распространены на другие значения вероятности охвата (и другие значения M).

В результате применения адаптивной процедуры должны быть определены:

а) оценка y величины Y,

б) соответствующая стандартная неопределенность u(y);

в) границы ylow и yhigh интервала охвата для Y, соответствующего заданной вероятности охвата.

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

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

2.1.7 Преимущества и недостатки метода Монте-Карло

Проанализируем недостатки и преимущества численного метода вычисления неопределенности на основе метода Монте-Карло.

Преимущества:

а) возможна оценка любых статистических характеристик результата измерения Y, а не только стандартного отклонения;

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

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

г) неопределенность входных величин может быть сколь угодно велика;

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

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

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

и) нет необходимости вычислять число эффективных степеней свободы по формуле Велча-Статерсвейта;

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

Недостатки:

а) необходимо иметь эффективные генераторы псевдослучайных чисел с длинным периодом;

б) сложные модели могут требовать много вычислительного времени для М реализаций;

в) коэффициенты чувствительности математической модели не могут быть получены;

г) математическая модель должна быть численно стабильной в связи с оцениванием не только в окрестности приписанного значения, но и на всех интервалах всех распределений входных величин [7].

Анализ недостатков:

а) любые современные математические программные пакеты содержат в себе ГПСЧ с длинным периодом (>> 106);

б) это не критично при расчете на ПЭВМ с помощью математических программных пакетов и, особенно, с использованием специализированного ПО, например, «Неопределенность 1.5» (демо-версия программы доступна на сайте http://www.novikov.biz.ua);

в) для работы алгоритма ММК не требуется знания коэффициентов чувствительности модели;

г) данный недостаток можно рассматривать, только если сама математическая модель измерения оценивалась численно.

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

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

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

GUM рассматривает ФПВ как фундамент оценивания неопределенности: в контексте закона распространения неопределенности он ясно ссылается на входные и выходную величины как описываемые или характеризуемые распределениями вероятностей.

Метод оценивания неопределенности GUM не определяет явно ФПВ для выходной величины. Однако на распределение вероятностей, используемое этим методом для определения характеристик выходной величины, иногда ссылаются как на «предусмотренное» или «следующее из» метода оценивания неопределенности GUM.

Так как метод оценивания неопределенности GUM явно использует только наилучшие оценки xi и связанные с ними стандартные неопределенности (а также ковариации и степени свободы, где это необходимо), он ограничен в информации, которую может предоставить о ФПВ для Y. По-существу, он ограничен предоставлением оценки y для Y и соответствующей стандартной неопределенности u(y), связанной с y, а также, возможно, связанных с ними (эффективных) степеней свободы.

Значения y и u(y) будут обоснованны для модели, которая линейна по X. Любая другая информация об Y, например интервалы охвата, получается с помощью дополнительных предположений о том, к примеру, что распределение для Y нормальное или масштабированное и сдвинутое t-распределение.

Некоторые преимущества ММК таковы:

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

б) улучшенная, в общем случае, оценка Y для нелинейных моделей;

в) улучшенная стандартная неопределенность, связанная с оценкой Y для нелинейных моделей, особенно когда Xi приписываются не гауссовские (например, асимметричные) ФПВ, без необходимости получать производные более высоких порядков;

г) предоставление интервала охвата, соответствующего оговоренной вероятности охвата, когда ФПВ для Y не может быть адекватно аппроксимирована нормальным распределением или масштабированным и сдвинутым t-распределением, т. е. когда центральная предельная теорема неприменима. Такая неадекватная аппроксимация может возникнуть, когда ФПВ, приписанная доминирующей Xi, не является распределением Гаусса или масштабированным и сдвинутым t-распределением, модель не линейна или ошибка аппроксимации, возникающая при использовании формулы Велча-Саттертвейта для эффективных степеней свободы, не является пренебрежимо малой;

д) при определении интервала охвата не требуется нахождение значения коэффициента охвата.

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

ММК имеет преимущества, когда:

а) линеаризация модели измеряемой величины не обеспечивает адекватного представления;

б) функция плотности вероятностей для выходной величины заметно отклоняется от распределения Гаусса или масштабированного и сдвинутого t-распределения, например, из-за явной асимметрии.

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

В случае б) результатом могут стать нереалистичные интервалы охвата (обобщение «расширенной неопределенности» в оценивании неопределенности GUM).

3 КОЛИЧЕСТВЕННОЕ СРАВНЕНИЕ РЕЗУЛЬТАТОВ ОЦЕНИВАНИЯ НЕОПРЕДЕЛЕННОСТИ ИЗМЕРЕНИЙ МЕТОДАМИ GUM И МОНТЕ-КАРЛО

Для количественного сравнения результатов оценивания неопределенности измерений методами GUM и Монте-Карло оценим неопределенность измерения объемного расхода природного газа в трубопроводе в рабочих условиях.

Объемный расход определяется по формуле из ГОСТ 8.586.1-2005(п. 5.1):

где d – диаметр отверстия сужающего устройства, м;

KШ – поправочный коэффициент, учитывающий шероховатость внутренней поверхности трубопровода;

KП – поправочный коэффициент, учитывающий притупление входной кромки диафрагмы;

С – коэффициент истечения;

Е – коэффициент скорости входа;

‑ коэффициент расширения среды;

‑ плотность среды, кг/м3;

Δp – перепад давления на диафрагме, Па.

Для упрощения расчетов предположим, что природный газ состоит исключительно из метана (т.е. примеси отсутствуют) и примем =0,668 кг/м3; KШ=1; KП=1; С =0,6050; Е=1,1426; =0,9964 (для следующих условий: t=20 oС; p=101,3 кПа).

Пусть , тогда .

Так как A=0,6619 (м3/кг)1/2, .

Рассчитаем неопределенность измерения методом GUM.

Базовый алгоритм оценивания неопределенности:

а) записываем модельное уравнение

б) входные величины равны:

x1 = (0,0800±0,0002) м,

x2 = (16,00±0,32)*103 Па;

в) вычисляем оценку результата измерения:

y=0,7578 м3/с;

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

- вычисление стандартной неопределенности по типу В проводим для симметричных границ по формуле

- вычисление вкладов неопределенности входных величин u(xi) и неопределенности измеряемой величины:

u(yi) =ciu(xi),

где сi – коэффициент чувствительности.

Коэффициенты чувствительности сi находят как частные производные выходной величины y по каждой из входной величин xi:

Коэффициенты чувствительности будут равны:

Вклад неопределенности измерения диаметра в неопределенность объемного расхода будет равен:

Вклад неопределенности измерения перепада давления в неопределенность объемного расхода будет равен:

3) Определяем суммарную стандартную неопределенность. Так как величины некоррелированы, то суммарная стандартная неопределенность рассчитывается по формуле:

4) Рассчитываем расширенную неопределенность по формуле:

где k – коэффициент охвата (в данном случае k=2).

е) Записываем результат измерения

Рассчитаем неопределенность измерения объемного расхода природного газа методом Монте-Карло.

а) Генерируем 2 массива случайных чисел (для входных величин d и Δp) объемом M=, подчиняющихся равномерным законам распределения с помощью программы STATGRAPHICS Plus v 5.0.

Фрагменты массивов полученных данных представлены на рисунках 6.1 и 6.2 (d– Col_1, Δp ‑Col_2).

Рисунок 6.1 – Массивы входных данных

Рисунок 6.2 – Массивы входных величин

б) Получаем массив оценки выходной величины qv (на рисунке 7 – Col_3). Сортировка массива по возрастанию и построение гистограммы выполнены в той же программе.

Рисунок 7 ‑ Массив оценки выходной величины qv

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

Рисунок 8 – Гистограмма распределения выходной величины

в) Вычисляем оценки параметров полученного распределения:

- математическое ожидание :

  • суммарная стандартная неопределенность:

  • расширенная неопределенность:

  • коэффициент охвата:

k=0,0182/0,0049=3,7.

г) Записываем результат измерения:

Сравним результаты измерения, полученные методами GUM и Монте-Карло:

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

Сформулируем основные преимущества оценивания неопределенности измерений методом Монте-Карло:

  • возможна оценка любых статистических характеристик результата измерения Y, а не только стандартного отклонения;

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

  • применимость для любых типов математических моделей результатов измерений, как линейных, так и нелинейных (не нужно определять необходимое количество членов ряда Тейлора для аппроксимации функции f);

  • неопределенность входных величин может быть сколь угодно велика;

  • нет необходимости знать или делать предположения о виде закона распределения выходной величины Y математической модели результата измерения;

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

  • нет необходимости оценивать коэффициенты чувствительности (частные производные первого порядка);

  • нет необходимости вычислять число эффективных степеней свободы по формуле Велча-Саттертвейта;

  • вычислительная сложность определяется числом М реализаций входных величин и временем вычисления функции f с последующей сортировкой массива данных.

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

Заключение

Результатом данной аттестационной работы стало выяснение различий между понятиями «погрешность измерений» и «неопределенность измерений», сравнение результатов оценки неопределенности измерений методом Монте-Карло и расчетным методом, изложенным в ИСО/МЭК «Руководство по выражению неопределенности измерений (GUM)».

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

Список использованных источников

  1. Захаров И. П. Теория неопределенности в измерениях. / И. П. Захаров, В. Д. Кукуш. ‑ Харьков: Консум, 2002. ‑ 256 с.

  2. РМГ 43-2001 ГСИ. Применение «Руководства по выражению неопределенности измерений».

  3. ГОСТ Р54500.3.1 – 2011/Руководство ИСО/МЭК 98-3:2008 «Неопределенность измерения. Часть 3. Руководство по выражению неопределенности измерения. Дополнение 1. Трансформирование распределений с использованием метода Монте-Карло».

  4. Артемьев Б.Г. Метрология и метрологическое обеспечение. М.: Стандартинформ, 2010.‑568 с.

  5. ГОСТ Р 8.563- 2009. ГСИ. Методики (методы) измерений

  6. ПМГ 96-2009.ГСИ. Результаты и характеристики качества измерений. Формы представления.

  7. Новиков В.В. Численные методы в вычислении неопределенности. / В. В. Новиков // Системи обробки інформації. – Харків, 2008. Вип.4 (71). С.126-128.

  8. Кузнецов В. П. Сопоставительный анализ погрешности и неопределенности измерений./ Кузнецов В. П.//Измерительная техника. ‑ 2003. ‑ №8. С. 18.

  9. ISO 21748:2010 «Руководство по использованию оценок повторяемости, воспроизводимости и достоверности при оценивании погрешностей измерений».

  10. Новиков В.В. Вычисление расширенной неопределенности / В. В. Новиков // Системи обробки інформації. – Харків, 2007, Вип.6 (64). С.73-77.

  11. Захаров И. П. Неопределенность измерений для чайников и … начальников: учеб. пособ. /И.П. Захаров, ‑ Харьков, 2013, ‑ 36 с.

  12. Международный словарь по метрологии. Основные и общие понятия и соответствующие термины: пер. с англ. и фр. / Всерос. науч.-исслед. ин-т метрологии им. Д. И. Менделеева, Белорус. гос. ин-т метрологии. Изд. 2-е, испр. ‑ СПб.: НПО «Профессионал», 2010. ‑ 82 с.

  13. РМГ 29-99 ГСИ. Метрология. Основные термины и определения.

  14. Novikov V.V. Numerical methods for uncertainty of measurements’ results calculation. // Конференція молодих учених із сучасних проблем механіки і математики імені академіка Я.С. Підстригача. Тези доповідей. – Львів, 2009. С. 147-148.

  15. ГОСТ Р 54500.3-2011 Неопределенность. Ч. 3. Руководство по выражению неопределенности измерения.

61

Просмотров работы: 12304