ФЕДЕРАЛЬНОЕ АГЕНТСТВО ПО ТЕХНИЧЕСКОМУ РЕГУЛИРОВАНИЮ И МЕТРОЛОГИИ
ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ УНИТАРНОЕ ПРЕДПРИЯТИЕ
«УРАЛЬСКИЙ НАУЧНО-ИССЛЕДОВАТЕЛЬСКИЙ ИНСТИТУТ МЕТРОЛОГИИ»
(ФГУП «УНИИМ»)

РЕКОМЕНДАЦИЯ

ГОСУДАРСТВЕННАЯ СИСТЕМА ОБЕСПЕЧЕНИЯ ЕДИНСТВА
ИЗМЕРЕНИЙ

УСТАНОВЛЕНИЕ ПРОСЛЕЖИВАЕМОСТИ
АТТЕСТОВАННЫХ ЗНАЧЕНИЙ СТАНДАРТНЫХ
ОБРАЗЦОВ

МИ 3174-2009

Екатеринбург
2009

ПРЕДИСЛОВИЕ:

1. РАЗРАБОТАНА: Федеральным государственным унитарным предприятием «Уральский научно-исследовательский институт метрологии» Ростехрегулирования.

ИСПОЛНИТЕЛИ: Д.П. Налобин канд. хим. наук, Е.В. Осинцева канд. хим. наук, А.В. Скутина.

2. УТВЕРЖДЕНА: ФГУП «УНИИМ».

3. ЗАРЕГИСТРИРОВАНА: ФГУП «ВНИИМС».

4. ВВЕДЕНЫ ВПЕРВЫЕ.

Содержание

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

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

3. Термины, определения и обозначения. 3

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

5. Способы установления прослеживаемости аттестованного значения стандартного образца. 5

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

7. Прослеживаемость при установлении аттестованного значения стандартного образца по стандартному образцу. 8

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

9. Прослеживаемость при установлении аттестованного значения стандартного образца посредством межлабораторного эксперимента. 15

10. Прослеживаемость при установлении аттестованного значения стандартного образца по процедуре приготовления. 17

Приложение А. Статистические таблицы (справочное) 17

Приложение Б. Оценивание неопределенности, обусловленной способом установления аттестованного значения государственного стандартного образца состава калия двухромовокислого 1-го разряда с применением первичного метода (справочное) 18

Библиография. 24

 

 

Государственная система обеспечения единства измерений

Установление прослеживаемости аттестованных значений
стандартных образцов

МИ 3174-2009

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

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

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

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

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

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

ГОСТ 8.315-97. Государственная система обеспечения единства измерений. Стандартные образцы состава и свойств веществ и материалов. Основные положения.

ГОСТ Р 8.563-96. Государственная система обеспечения единства измерений. Методики выполнения измерений.

ГОСТ Р ИСО 5725-1-2002. Точность (правильность и прецизионность) методов и результатов измерений. Часть 1. Основные положения и определения.

ГОСТ Р ИСО 5725-2-2002. Точность (правильность и прецизионность) методов и результатов измерений. Часть 2. Основной метод определения повторяемости и воспроизводимости стандартного метода измерения.

ГОСТ Р ИСО 5725-3-2002. Точность (правильность и прецизионность) методов и результатов измерений. Часть 3. Промежуточные показатели прецизионности стандартного метода измерения.

ГОСТ Р ИСО 5725-6-2002. Точность (правильность и прецизионность) методов и результатов измерений. Часть 6. Использование значений точности на практике.

ГОСТ Р 50779.10-2000 (ИСО 3534.1-93). Статистические методы. Вероятность и основы статистики. Термины и определения.

РМГ 53-2002. Государственная система обеспечения единства измерений. Стандартные образцы. Оценивание метрологических характеристик с использованием эталонов и образцовых средств измерений.

РМГ 29-99. Государственная система обеспечения единства измерений. Метрология. Основные термины и определения.

Р 50.2.056-2007. Государственная система обеспечения единства измерений. Образцы материалов (веществ) стандартные. Термины и определения.

Р 50.2.058-2007. Государственная система обеспечения единства измерений. Оценивание неопределенностей аттестованных значений стандартных образцов.

3. Термины, определения и обозначения

3.1. В настоящих рекомендациях использованы термины из ГОСТ 8.315, ГОСТ Р 8.563, ГОСТ Р 50779.10 (ИСО 3534-1), ГОСТ Р ИСО 5725-1, ГОСТ Р ИСО 5725-3, РМГ 29, Р 50.2.058, Р 50.2.056 и термины со следующими определениями:

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

3.1.2. стандартная неопределенность: Неопределенность результата измерений, выраженная в виде стандартного отклонения [1].

3.1.3. способ установления аттестованного значения СО: Метрологически обоснованная процедура установления аттестованного значения СО.

3.1.4. неопределенность, обусловленная способом установления аттестованного значения СО (неопределенность от способа): Составляющая неопределенности аттестованного значения СО, обусловленная неопределенностями результатов измерений при установлении аттестованного значения СО.

3.1.5. стандартная неопределенность, обусловленная способом установления аттестованного значения СО (стандартная неопределенность от способа): Неопределенность, обусловленная способом установления аттестованного значения СО, выраженная в виде стандартного отклонения.

3.1.6. прослеживаемость: Свойство результата измерения или значения эталона, заключающееся в возможности его соотнесения с принятыми реперами - обычно национальными или международными эталонами - посредством непрерывной цепи калибровок или сличений с установленными неопределенностями [3].

3.1.7. прослеживаемость аттестованного значения СО: Свойство аттестованного значения СО, заключающееся в возможности его соотнесения с принятыми реперами, использованными при установлении аттестованного значения СО, с установленными неопределенностями.

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

Примечание - Номенклатура принятых в настоящих рекомендациях реперов приведена в 5.2.

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

3.1.10. стандартная неопределенность, обусловленная прецизионностью: Неопределенность, обусловленная прецизионностью выраженная в виде стандартного отклонения.

3.1.11. неопределенность, обусловленная использованными реперами (неопределенность от использованных реперов): Составляющая неопределенности от способа установления аттестованного значения СО, обусловленная неопределенностями реперов, использованных при установлении аттестованного значения СО.

3.1.12. стандартная неопределенность, обусловленная использованными реперами (стандартная неопределенность от использованных реперов): Неопределенность, обусловленная использованными реперами, выраженная в виде стандартного отклонения.

3.1.13. суммарная стандартная неопределенность: Стандартная неопределенность результата измерений, полученного через значения других величин, равная положительному квадратному корню суммы членов, причем члены являются дисперсиями или ковариациями этих других величин, взвешенные в соответствии с тем, как результат измерений изменяется при изменении этих величин [1].

3.1.14. расширенная неопределенность: Величина, определяющая интервал вокруг результата измерений, в пределах которого, как можно ожидать, находится большая часть распределения значений, которые с достаточным основанием могли быть приписаны измеряемой величине [1].

3.1.15. метод измерений: Логическая последовательность операций, описанная в общем виде, которая применяется при выполнении измерений; [2].

3.1.16. эмпирические методики выполнения измерений (эмпирические МВИ): МВИ, в которых единицы измерения результатов зависят от применяемого метода [2].

3.1.17. рациональные методики выполнения измерений (рациональные МВИ): МВИ, в которых единицы измерения результатов не зависят от применяемого метода [2].

3.1.18. первичный метод измерений: Метод измерений, обладающий наивысшими метрологическими свойствами, операции которого полностью описываются и истолковываются в единицах СИ, и результаты которого принимаются без сравнения с эталоном той же величины [3].

Примечание

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

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

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

СО - стандартный образец;

МВИ - методика выполнения измерений;

u - стандартная неопределенность;

U - расширенная неопределенность;

uchar - стандартная неопределенность, обусловленная способом установления аттестованного значения СО;

uprec - стандартная неопределенность, обусловленная прецизионностью;

utr - стандартная неопределенность от использованных реперов;

uс - суммарная стандартная неопределенность;

uc,char - суммарная стандартная неопределенность, обусловленная способом установления аттестованного значения СО;

А - аттестованное значение СО;

 - аттестуемая характеристика СО;

σr - стандартное отклонение повторяемости;

σI(z) - стандартное отклонение промежуточной прецизионности при различии по фактору z;

R - диапазон значений результатов измерений;

I - значение выходного сигнала средства измерений;

sc(y) - стандартное отклонение не исключенной систематической погрешности МВИ.

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

4.1. В соответствии с Р 50.2.056, [4 - 6] одним из важных характеристик стандартного образца (СО) является прослеживаемость аттестованного значения СО к единице СИ или принятым реперам.

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

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

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

4.3. Существуют СО, аттестованные значения которых не могут иметь строгую прослеживаемость к единицам СИ. В этих случаях устанавливают прослеживаемость аттестованного значения СО до соответствующих реперов, например [5, 6]:

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

- до согласованного значения результатов измерений, полученных в строгом соответствии с эмпирической МВИ, или до согласованных эталонов, четко описанных и принятых всеми заинтересованными сторонами.

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

5. Способы установления прослеживаемости аттестованного значения стандартного образца

5.1. Имеются следующие способы установления прослеживаемости аттестованного значения СО:

1) использование для установления аттестованного значения СО методики выполнения измерений (МВИ), основанной на первичном методе (далее - прослеживаемость при установлении аттестованного значения СО первичным методом);

2) использование для установления аттестованного значения СО другого СО, имеющего установленную прослеживаемость (далее - прослеживаемость при установлении аттестованного значения СО по стандартному образцу);

3) использование для установления аттестованного значения СО МВИ, предусматривающей градуировку по СО, прослеживаемость аттестованного значения которого установлена (далее - прослеживаемость при установлении аттестованного значения СО с использованием градуировки);

4) установление аттестованного значения СО, выполняемое с привлечением нескольких лабораторий (далее - прослеживаемость при межлабораторном способе установления аттестованного значения СО), предусматривающее применение:

- одной эмпирической МВИ,

- одной рациональной МВИ, прослеживаемость которой установлена при ее аттестации,

- нескольких МВИ, основанных на различных физических и химических принципах;

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

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

5.2. В качестве реперов в данных рекомендациях приняты:

1) значения эталонов;

2) аттестованные значения СО, прослеживаемость которых установлена;

3) результаты измерений по МВИ, прослеживаемость которых установлена при аттестации МВИ;

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

5) стандартные справочные данные.

5.3. Неопределенность от способа установления аттестованного значения СО имеет два источника неопределенности:

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

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

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

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

5.4. Суммарную стандартную неопределенность, обусловленную способом установления аттестованного значения СО uс,char, оценивают по формуле

(5.1)

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

Одним из способов установления прослеживаемости аттестованного значения СО является использование при установлении аттестованного значения СО методики выполнения измерений (МВИ), основанной на первичном методе.

6.1. Значение аттестуемой характеристики аттестуемого СО , определяемое по МВИ, основанной на первичном методе, устанавливают по уравнению

(6.1)

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

Аттестованное значение À аттестуемой характеристики при определенных условиях проведения измерений I (w1, ... ,wm) получают по результатам измерения входных величин

(6.2)

где х1, х2, ... xm - измеренные или приписанные значения входных величин,

w1wk - параметры, характеризующие условия выполнения измерений.

6.2. Условия проведения измерений:

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

- для МВИ, используемой при установлении аттестованного значения СО, должно быть известно значение стандартного отклонения промежуточной прецизионности σI(z), оцененного в соответствии с ГОСТ Р ИСО 5725-3;

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

Значение i-го результата измерения аттестуемой характеристики представляют в виде

(6.3)

где εi - случайная составляющая погрешности i-го результата измерения значения аттестуемой характеристики, полученной в условиях промежуточной прецизионности.

Предполагают, что математическое ожидание M(εi) = 0 и дисперсия

6.3. Для установления аттестованного значения СО выполняют в соответствии с МВИ измерения аттестуемой характеристики в пробах материала СО, получая результаты измерений A1, … ,An.

Вычисляют:

- среднее значение

(6.4)

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

(6.5)

- диапазон значений результатов измерений R (разность между наибольшим Аmax и наименьшим Amin результатом измерения)

R = Amax - Amin,

(6.6)

- значение критического диапазона для уровня вероятности 95 %

(6.7)

где f(n) - значение коэффициента критического диапазона, зависящего от числа результатов измерений n, возможные значения f(n) приведены в таблице А.1 приложения A; σI(z) - стандартное отклонение промежуточной прецизионности МВИ.

Сравнивают полученный диапазон R с критическим диапазоном, проверяют выполнение неравенства

(6.8)

Если неравенство (6.8) выполняется, то все результаты считают приемлемыми и аттестованное значение СО оценивают как среднее арифметическое значение, вычисляемое по формуле (6.4)

(6.9)

а значение стандартной неопределенности, обусловленной прецизионностью, оценивают по формуле

(6.10)

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

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

6.5. Стандартную неопределенность от используемого репера utr оценивают с учетом значений стандартной неопределенности входных величин u(хi) и коэффициентов чувствительности ci по формуле

(6.11)

где r(xi, xj) - коэффициент корреляции входных величин xi и xj.

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

(6.12)

Для статистически независимых входных величин формула (6.11) преобразуется к виду

(6.13)

Суммарную стандартную неопределенность способа установления аттестованного значения СО uc,char оценивают по формуле (5.1).

6.6. Пример установления аттестованного значения с применением первичного метода СО состава калия двухромовокислого 1-го разряда с применением метода кулонометрического титрования приведен в приложении Б.

7. Прослеживаемость при установлении аттестованного значения стандартного образца по стандартному образцу

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

7.2. Для аттестованного значения используемого при испытаниях СО A0 должна быть определена стандартная неопределенность u(A0) и число степеней свободы .

7.3. Измерения, проводимые при установлении аттестованного значения СО, должны выполняться по МВИ, стандартное отклонение повторяемости σr которой должно быть оценено в соответствии ГОСТ Р ИСО 5725-2.

7.4. Стандартное отклонение повторяемости используемой МВИ и стандартная неопределенность аттестованного значения используемого СО u(A0) должны удовлетворять следующим условиям

(7.1)

где Uдоп - допускаемое значение расширенной неопределенности аттестованного значения аттестуемого СО.

7.5. Измерение значения аттестованной характеристики проводят в пробах материала каждого СО. Получают в условиях повторяемости 2n результатов y0i и y1i (результаты с индексом 0 для проб материала используемого при испытаниях СО, а результаты с индексом 1 для проб материала аттестуемого СО). Очередность получения результатов y0i и y1i определяют случайным образом.

Минимальное количество пар результатов nmin определяют, исходя из соотношений между стандартной неопределенностью u(A0), стандартным отклонением повторяемости σr МВИ и Uдоп из неравенства

(7.2)

7.6. Результат измерения аттестуемой характеристики аттестуемого СО представляют в виде

yi = Ai + ∆c + εi,

(7.3)

где Ai - значение аттестуемой характеристики аттестуемого СО,

с - систематическое смещение результатов измерений относительно аттестованного значения А01,

εi - случайная составляющая погрешности i-го результата измерения значения аттестуемой характеристики, полученной в условиях повторяемости.

Относительно величин, входящих в модель (7.3) принимают следующие предположения:

- математическое ожидание равно нулю

M(εi) = 0,

(7.4)

- дисперсия εi (D(εi)) равна квадрату стандартного отклонения повторяемости МВИ

(7.5)

- математическое ожидание результата измерения yi (M(yi)) при выполнении предположения (7.4) равно

(7.6)

где М(А) - математическое ожидание результатов измерений аттестуемой характеристики СО.

При измерении аттестуемой характеристики в пробах используемого СО математическое ожидание значений аттестованной характеристики СО в уравнении (7.6) приравнивают аттестованному значению СО

(7.7)

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

(7.8)

Из уравнений (7.6), (7.7) и (7.8) получают уравнение для оценивания систематического смещения результатов измерений

(7.9)

Стандартную неопределенность оценивания систематического смещения результатов измерений () аттестованной характеристики используемого СО относительно аттестованного значения используемого СО (А0) - стандартную неопределенность от используемого репера utr - оценивают по формуле

(7.10)

Оценивают относительное смещение результатов измерений при значении аттестуемой характеристики равной А0

(7.11)

Систематическое смещение среднего значения результатов измерений  аттестуемой характеристики аттестуемого СО оценивают по формуле

(7.12)

_________

1 Систематическое смещение результатов измерений включает систематическую погрешность МВИ и систематическую погрешность лаборатории

7.7. После выполнения всех измерений вычисляют диапазон результатов измерений

(7.13)

 

(7.14)

где y0,max, y0,min и y1,max, y1,min - наибольшие и наименьшие результаты измерения из всех результатов измерений y0i и y1i, соответственно.

7.8. Вычисляют значение критического диапазона CR0,95(n) по формуле

(7.15)

где f(n) - значение коэффициента критического диапазона, зависящего от числа результатов измерений, устанавливают по таблице А.1 приложения А.

Сравнивают диапазоны R0 и R1 с критическим диапазоном CR0,95(n), т.е. проверяют выполнение неравенств

(7.16)

7.9. Если неравенства (7.16) выполняются, все полученные результаты измерений принимают и используют для оценивания аттестованного значения.

В противном случае выясняют причину большого рассеяния результатов измерений, либо исключают часть результатов измерений, либо повторяют процедуру измерения в соответствии с пунктами 7.5 - 7.8.

7.10. В случае выполнения неравенств (7.16) по полученным результатам вычисляют:

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

(7.17)

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

(7.18)

 

(7.19)

- коэффициент корреляции результатов измерений

(7.20)

7.11. Проверяют статистическую гипотезу о равенстве среднего значения  и аттестованного значения А0. Для проверки этой гипотезы вычисляют отношение

(7.21)

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

(7.22)

Сравнивают отношение  с квантилем распределения Стьюдента с числом степеней свободы

(7.23)

7.12. Если неравенство (7.23) выполняется, то гипотезу о равенстве среднего значения  и аттестованного значения  принимают (т.е. систематическое смещение результатов ∆c = 0) и аттестованное значение аттестуемого СО À оценивают по формуле

(7.24)

Суммарную стандартную неопределенность от способа установления аттестованного значения СО uc,char оценивают по формуле

(7.25)

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

7.13. При невыполнении неравенства (7.23) гипотезу о равенстве среднего значения  и аттестованного значения A0 используемого СО отвергают, систематическое смещение результатов оценивают по формуле (7.12) и аттестованное значение СО À аттестуемого СО оценивают по формуле

(7.26)

Стандартную неопределенность от используемого репера  и стандартную неопределенность, обусловленную прецизионностью uprec, оценивают по формулам:

(7.27)

 

(7.28)

Суммарную стандартную неопределенность от способа установления аттестованного значения СО uc,char в этом случае оценивают по формуле

(7.29)

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

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

8.1. Преобразование градуировочной характеристики

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

(8.1)

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

.

(8.2)

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

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

Примечание - В работе [8] показано, что количество экземпляров СО в градуировочном комплекте не влияет на величину относительной стандартной неопределенности от использованных реперов, если между аттестованными значениями СО комплекта существует статистическая зависимость. Это будет, например, в случае, когда материалы СО комплекта получены путем разбавления исходного вещества или аттестованные значения СО комплекта определены по одной МВИ. В этом случае следует, очевидно, ограничиться минимальным числом СО для определения градуировочной характеристики, т.е. использовать для характеризации СО один или два СО из комплекта, наиболее близких по значениям выходных сигналов к выходным сигналам от материала аттестуемого СО.

8.2. Градуировка по одному стандартному образцу

8.2.1. Аттестация СО по МВИ, предусматривающей градуировку по одному СО, имеющему установленную прослеживаемость, возможна если:

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

- функциональная зависимость между значениями аттестуемой характеристики СО Ai и величиной выходного сигнала средства измерения Ii, имеет вид

.

(8.3)

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

8.2.2. Для оценивания аттестованного значения А аттестуемого СО с помощью используемого для исследований СО попеременно n раз измеряют выходные сигналы от проб материала используемого СО I0i и аттестуемого СО Iai. Обработку результатов измерений проводят следующим образом. Вычисляют:

- для каждой i-й пары выходных сигналов отношение ϑi

(8.4)

- среднее значение отношений выходных сигналов

(8.5)

- стандартное отклонение среднего отношения выходных сигналов

(8.6)

8.2.3. Аттестованное значение СО оценивают по формуле

(8.7)

Стандартные неопределенности, обусловленные используемыми реперами utr и прецизионностью uprec, оценивают по формулам

(8.8)

 

(8.9)

Суммарную стандартную неопределенность, обусловленную способом установления аттестованного значения СО uc,char, вычисляют по формуле

(8.10)

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

8.3. Градуировка по двум стандартным образцам

8.3.1. Аттестация СО по МВИ, предусматривающей градуировку по двум СО, имеющим установленную прослеживаемость, возможна если:

- имеется два СО, близких по составу к материалу аттестуемого СО, аттестованные значения первого СО A1, второго СО - А2, значение A1 должно быть меньше А2 (А1 < А2), значение аттестуемой характеристики аттестуемого СО должно находится в интервале [А1, А2];

- для любого значения аттестуемой характеристики из интервала [А1, А2] зависимость между значениями аттестуемой характеристики СО Ai и величиной выходного сигнала используемого средства измерения Ii имеет вид

Ai = b + a·Ii.

(8.11)

8.3.2. Для оценивания аттестованного значения  получают от используемых при испытаниях двух СО попеременно n раз выходные сигналы I1i и I2i и выходные сигналы от проб аттестуемого СО Iai. Обработку результатов проводят следующим образом. Вычисляют:

- для каждого i-го измерения выходных сигналов значение отношения

(8.12)

- значение аттестуемой характеристики Аi

(8.13)

8.3.3. После выполнения всех измерений вычисляют:

- средние значения отношений выходных сигналов

(8.14)

- стандартные неопределенности отношений выходных сигналов

(8.15)

 

(8.16)

8.3.4. Аттестованное значение СО оценивают по формуле

(8.17)

Суммарную стандартную неопределенность, обусловленную способом установления аттестованного значения СО uc,char, оценивают с учетом составляющих:

(8.18)

 

(8.19)

где  - коэффициент корреляции отношений  и  и r(A1, A2) - коэффициент корреляции аттестованных значений СО A1 и A2.

8.3.5. Из определения параметров ϑ1i и ϑ2i по уравнению (8.12) следует равенство

ϑ1i + ϑ2i = 1,

(8.20)

из которого для средних значений получают уравнение

(8.21)

Следовательно, коэффициент корреляции r(ϑ1, ϑ2) = -1 и их стандартные отклонения равны. С учетом этого уравнение (8.19) может быть представлено

(8.22)

где стандартную неопределенность  оценивают по уравнению

(8.23)

8.3.6. В зависимости от способов установления аттестованных значений двух используемых СО A1 и A2 возможны два варианта оценки суммарной стандартной неопределенности uc,char.

Если аттестованные значения A1 и A2 установлены по результатам независимых измерений, то коэффициент корреляции r(А1, А2) равен нулю и уравнение (8.18) имеет вид

(8.24)

где u(А1) и u(А2) - стандартная неопределенность аттестованных значений А1 и А2 первого и второго используемых СО соответственно.

Если между аттестованными значениями A1 и А2 существует функциональная зависимость, то коэффициент корреляции r(A1, A2) = 1 и уравнение (8.18) имеет вид

(8.25).

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

9.1. При установлении аттестованного значения СО по результатам межлабораторного эксперимента в зависимости от применяемых МВИ и вида измеряемой величины возможны следующие ситуации:

1) измеряемая величина определяется методом, измерения во всех лабораториях выполняют по одной эмпирической МВИ;

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

3) измеряемая величина не зависит от метода, измерения выполняют по нескольким МВИ, основанным на различных физических и химических принципах.

9.2. В случае перечисления 1) из 9.1 согласованность результатов, полученных в разных лабораториях, оценивание аттестованного значения СО и стандартной неопределенности способа установления аттестованного значения СО производят в соответствии с разделом 7.2.3 Рекомендации Р 50.2.058.

9.3. В случае перечисления 2) из 9.1 для применяемой МВИ должны быть известны характеристики прецизионности, установленные в соответствии с ГОСТ Р ИСО 5725-2, и стандартное отклонение не исключенной систематической погрешности МВИ sc(y).

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

Аттестованное значение аттестуемой характеристики СО оценивают средним арифметическим значением

(9.1)

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

Стандартную неопределенность, обусловленную прецизионностью, оценивают по формуле

(9.2)

Стандартную неопределенность от используемых реперов оценивают по значению стандартной неопределенности от реперов МВИ

(9.3)

9.4. В случае перечисления 3) из 9.1 в межлабораторном эксперименте принимают участие р лабораторий, которые применяют q МВИ (q ≥ 2), основанных на различных физических или химических принципах. Если лаборатория применяет в межлабораторном эксперименте несколько МВИ, то результаты измерений этой лаборатории по каждой МВИ представляют как результаты отдельной лаборатории.

Предполагают, что каждая j-я МВИ имеет оцененные в соответствии с ГОСТ Р ИСО 5725-2 характеристики прецизионности (стандартное отклонение повторяемости и стандартное отклонение воспроизводимости σRj и стандартные отклонения не исключенной систематической погрешности МВИ sj,c(y)).

После представления всех результатов для каждой j-й МВИ вычисляют:

- среднее значение в i-й лаборатории  по j-й МВИ

(9.4)

где  - k-й результат, полученный в i-й лаборатории по j-й МВИ;

- общее среднее значение по всем результатам, полученным по j-й МВИ

(9.5)

где pj - число лабораторий, представивших результаты по j-й МВИ,

- стандартное отклонение средних значений по j-й МВИ

(9.6)

Оценку согласованности результатов, полученных в разных лабораториях по каждой МВИ, согласованность средних результатов, полученных по всем МВИ проводят в соответствии с подразделом 7.2.4 Р 50.2.058.

Аттестованное значение СО рассчитывают, как среднее арифметическое значение

(9.7)

где q - число МВИ, используемых при характеризации СО.

Стандартные неопределенности аттестованного значения СО, обусловленные прецизионностью и использованными реперами, оценивают соответственно по формулам

(9.8)

 

(9.9)

9.5. Суммарную стандартную неопределенность, обусловленную способом установления аттестованного значения СО uc,char для описанных случаев оценивают по формуле (5.1).

10. Прослеживаемость при установлении аттестованного значения стандартного образца по процедуре приготовления

10.1. Аттестуемую характеристику СО  выражают в виде зависимости от двух групп входных величин

(10.1)

где X1, ... ,Xn - значения аттестуемой характеристики в исходных веществах, используемых для приготовления материала СО;  - количества исходных веществ.

10.2. Аттестованное значение СО À оценивают по оценкам входных величин, используя зависимость (10.1) по уравнению

(10.2)

где x1, … ,xn - оценки значения аттестуемой характеристики в исходных веществах, имеющие прослеживаемость к выбранными по 5.2 реперами; q1, ... ,qn - измеренные при приготовлении материала СО количества исходных веществ (массы, объёмы), определенные с применением поверенных средств измерений.

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

10.3. Стандартную неопределенность utr оценивают по стандартным неопределенностям входных величин и коэффициентов чувствительности, описывающих изменение аттестуемой характеристики при изменении i-ой входной величины, по уравнению

(10.3)

10.4. Стандартную неопределенность uрrес оценивают по стандартным отклонениям случайных погрешностей средств измерений, используемых при измерении количеств исходных веществ по уравнению

(10.4)

где s(qi) - стандартное отклонение случайной погрешности средства измерений количества i-го исходного вещества при приготовлении материала СО.

Приложение А

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

Статистические таблицы

Таблица А.1 - Значение коэффициентов для вычисления критического диапазона

n

f(n)

n

f(n)

2

2,8

25

5,2

3

3,3

26

5,2

4

3,6

27

5,2

5

3,9

28

5,3

6

4,0

29

5,3

7

4,2

30

5,3

8

4,3

31

5,3

9

4,4

32

5,3

10

4,5

33

5,4

11

4,6

34

5,4

12

4,6

35

5,4

13

4,7

36

5,4

14

4,7

37

5,4

15

4,8

38

5,5

16

4,8

39

5,5

17

4,9

40

5,5

18

4,9

45

5,6

19

5,0

50

5,6

20

5,0

60

5,8

21

5,0

70

5,9

22

5,1

80

5,9

23

5,1

90

6,0

24

5,1

100

6,1

Таблица А.2 - Квантили распределения Стьюдента t0,95(v)

v

t0,95(v)

v

t0,95(v)

v

t0,95(v)

1

12,706

11

2,201

22

2,074

2

4,303

12

2,179

24

2,064

3

3,182

13

2,160

26

2,056

4

2,776

14

2,145

28

2,048

5

2,571

15

2,131

30

2,042

6

2,447

16

2,120

32

2,037

7

2,365

17

2,110

34

2,032

8

2,306

18

2,101

36

2,028

9

2,262

19

2,093

38

2,024

10

2,228

20

2,086

40

2,021

Приложение Б

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

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

Б.1. Способ установления аттестованного значения ГСО состава калия двухромовокислого 1-го разряда

Установление массовой доли основного вещества в материале Государственного стандартного образца состава калия двухромовокислого 1-го разряда (далее ГСО-КД) выполняют методом косвенных измерений с использованием установки высшей точности кулонометрического титрования УВТ 14-А-81 по реакции окисления-восстановления. Аттестуемая характеристика ГСО-КД массовая доля калия двухромовокислого, %.

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

Fe3+ + e → Fe2+

(катод)

Cr6+ + 3Fe2+ → Cr3+ + 3Fe3+

(раствор)

Массу прореагировавшего бихромата калия определяют по количеству электричества, прошедшего через раствор от начала титрования до достижения точки конца титрования (ТКТ), согласно закону Фарадея

(Б.1)

где mэксп - масса определяемого вещества, г;

М - молекулярная масса определяемого вещества, для бихромата калия равная 294,1846 г/моль;

Q - количество электричества, Кл;

z - зарядовое число иона, участвующего в реакции, равное для бихромат-иона 6;

F - постоянная Фарадея, принятая равной 96485,3383 Кл/моль.

Количество электричества Q, Кл, рассчитывают по формуле

Q = I·t,

(Б.2)

где I - значение генераторного тока, проходящего через электролит, А,

t - время прохождения тока, с.

Поскольку значения М, z, F для данной реакции величины постоянные, то формулу (1) можно представить в виде

(Б.3)

Силу генераторного тока I, А, рассчитывают по падению напряжения U на термостатированной катушке с сопротивлением R по формуле

(Б.4)

В значение R вводят температурную поправку согласно паспорту на катушку сопротивления. Время измеряют частотомером электронно-счетным Ф 5041, который включается и выключается синхронно с включением и выключением тока.

Для титрования используют электролитическую ячейку типа «двойной защитной трубки». Точку конца титрования (ТКТ) определяют амперометрическим методом с использованием в качестве индикаторных электродов двух поляризуемых пластинчатых платиновых электродов.

Ниже представлен алгоритм оценивания неопределенности установления аттестованного значения ГСО состава калия двухромовокислого (далее ГСО-КД).

Б.2. Уравнение измерений

Значение массовой доли основного вещества в ГСО калия двухромовокислого 1-го разряда по методу кулонометрического титрования рассчитывают по уравнению

(Б.5)

где : А - массовая доля основного вещества в ГСО, %.

Для входных величин в формуле (Б.1) приняли следующие обозначения:

x1 = U1 - значение напряжения при предварительном титровании, В;

x2 = U2 - значение напряжения при основном титровании, В;

x3 = U3 - значение напряжения при заключительном титровании, В;

x4 = t'1 - общее время предварительного титрования, с;

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

x6 = t2 - время основного титрования, с;

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

x8 = M - молекулярная масса калия двухромовокислого, г/моль (М = 294,1846 г/моль);

x9 = F - постоянная Фарадея, Кл/моль (F = 96485,3383 Кл/моль);

x10 = m - масса пробы с учетом поправки на подъемную силу воздуха, г;

x11 = R - сопротивление катушки, Ом.

z - зарядовое число бихромат-иона (z = 6);

Множитель 100 переводит значение А из долей в проценты.

Б.3. Оценивание стандартной неопределенности, обусловленной прецизионностью

Для оценивания аттестованного значения СО выполнили в соответствии с МВИ 11 измерений массовой доли основного вещества в пробах материала СО. Результаты измерений приведены в таблице Б.1.

Таблица Б.1 - Результаты измерения массовой доли основного вещества

Номер результата i

Результат Ai

Номер результата i

Результат Ai

1

99,966

1

99,997

2

99,964

8

99,964

3

99,987

9

99,997

4

99,968

10

99,967

5

99,987

11

99,968

6

99,971

-

-

По полученным результатам вычисляют:

- по формуле (7.4) среднее значение

(Б.6)

- по формуле (7.5) стандартное отклонение

(Б.7)

Находим наибольшее значение Аmах = 99,997 и наименьшее значение Amin = 99,964. Вычисляют диапазон значений результатов

(Б.8)

Стандартное отклонение повторяемости, оцененное при характеризации МВИ, σr = 0,015. По таблице А.1 находят для n = 11 значение функции f(11) = 4,6 и определяют по формуле (7.6) критический диапазон

CR0,95(11) = f(11)·σr = 4,6·0,015 = 0,069.

(Б.9)

Диапазон значений результатов W меньше критического диапазона CR0,95(11), следовательно все результаты считают приемлемыми и в качестве аттестованного значения СО принимают среднее значение

(Б.10)

Стандартную неопределенность от прецизионности uprec вычисляют по формуле (7.8)

(Б.11)

Б.4. Оценивание стандартной неопределенности от используемых реперов

Б.3.1. Неопределенность от используемых реперов определяют композицией составляющих неопределенности, обусловленных взвешиванием пробы (m), измерением напряжения генерирования (U1, U2, U3) на катушке сопротивления (R), измерением времени титрования (t1, t1, t2, t3), определением молекулярной массы (М) и константы Фарадея (F).

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

Б.3.2. Стандартную неопределенность от используемых реперов рассчитывают по формуле

(Б.12)

где xj - оценка j-ой входной величины;

u(xj) - стандартная неопределенность j-ой входной величины;

cj - коэффициенты чувствительности, описывающие изменение выходной величины А (массовой доли основного вещества) при изменении j-ой входной величины xj.

Б.3.3. Коэффициенты чувствительности cj для каждой входной величины xj определяют как частные производные функции формулы (Б.1). Для расчета численных значений коэффициентов чувствительности используют соответствующие экспериментальные данные, усредненные по серии измерений МДОВ в ГСО

(Б.13)

 

(Б.14)

 

(Б.15)

 

(Б.16)

 

(Б.17)

 

(Б.18)

 

(Б.19)

 

(Б.20)

 

(Б.21)

 

(Б.22)

 

(Б.23)

Б.3.4. Стандартную неопределенность u(xi) для входных величин x4, x6, x10, x11,… для которых в документации на средства измерения, применяемые для измерения этих величин, заданы границы значений погрешности , оценивают исходя из предположения о равномерном распределении вероятности погрешностей по формуле

(Б.24)

Стандартная неопределенность измерения напряжения генерирования u(xi) для входных величин х1, х2, х3, является комбинацией двух составляющих, обусловленных:

а) погрешностью компаратора напряжений Р3003;

б) погрешностью нормального элемента Х480, используемого для калибровки Р3003.

Предполагая равномерное распределение вероятности, неопределенность измерения напряжения генерирования оценивают по формуле

(Б.25)

Примечание - характеристика погрешности компаратора напряжений Р3003, ∆i Р3003, зависит от величины измеряемого напряжения.

Стандартная неопределенность u(xi) для входных величин x5 и x7, значения которых находят в результате обработки экспериментальной кривой титрования, является комбинацией двух составляющих неопределенности:

1) составляющей u1(хi), обусловленной погрешностью частотомера Ф 5041 и рассчитываемой по формуле (Б.24);

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

Таким образом, стандартную неопределенность u(t1) оценивают по формуле

(Б.26)

Оценивание стандартной неопределенности u2(xi) для амперометрического метода индикации точки конца титрования (ТКТ) производят следующим образом.

Экспериментальные данные для построения кривой титрования Iинд = f(t) разбивают на две группы - точки до и после возрастания индикаторного тока. Критерием отнесения текущей точки к группе возрастания индикаторного тока является превышение величины индикаторного тока по сравнению со значением в предыдущей точке более, чем на 0,050 мкА. Используя положения метода наименьших квадратов, группы точек до и после возрастания индикаторного тока описывают с помощью линейных уравнений в виде

Iинд = ax + b

(Б.27)

Приравнивая уравнения прямых до и после возрастания индикаторного тока, находят значение времени в ТКТ по формуле

(Б.28)

где а1, b1 - коэффициенты в уравнении прямой до роста индикаторного тока;

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

x - время в точке конца титрования, полученное расчетным методом.

Примечание - в формулах (Б.22) и (Б.23) индексы при х опущены, эти формулы одинаково справедливы для x5 и для х7.

Рассчитывают средние квадратические отклонения коэффициентов а и b, S(a) и S(b) соответственно, по формулам

(Б.29)

 

(Б.30)

где Sy - среднее квадратическое отклонение значений индикаторного тока, полученных экспериментально, , от рассчитанных по уравнению (Б.28), ;

l - число точек, использованных для построения прямой.

Sy вычисляют по формуле

(Б.31)

Выполняют вычисления по формулам (Б.31) - (Б.31) для каждого результата измерений МДОВ в серии, получают n значений Si(a1), Si(a2), Si(b1) и Si(b2). Стандартную неопределенность нахождения коэффициента a1, u(a1), принимают за среднее квадратическое значение, рассчитываемое по формуле

(Б.32)

Аналогично оценивают стандартные неопределенности нахождения коэффициентов а2, b1 и b2 - u(a2), u(b1) и u(b2) соответственно. Полученные значения затем подставляют в формулу (Б.32).

Коэффициенты чувствительности c(a1), с(а2), с(b1) и с(b2) рассчитывают как частные производные формулы (Б.32)

(Б.33)

 

(Б.34)

 

(Б.35)

 

(Б.36)

Для нахождения численных значений коэффициентов чувствительности используют средние по серии измерений МДОВ значения а1, a2, b1 и b2.

Стандартную неопределенность u2(xi) рассчитывают по формуле

(Б.37)

Стандартная неопределенность u(x8), обусловленная неопределенностью определения молекулярной массы анализируемого вещества, рассчитывают по значениям неопределенности оценки атомных масс элементов, входящих в состав молекулы, согласно справочным данным [9], по формуле

(Б.38)

где U(ArK) - неопределенность определения атомной массы калия, равная 0,0001 [9];

U(АrCr) - неопределенность определения атомной массы хрома, равная 0,0006 [9];

U(АrO) - неопределенность определения атомной массы кислорода, равная 0,0003 [9].

Стандартную неопределенность u(F), обусловленную неопределенностью определения константы Фарадея, согласно справочным данным [10], оценивают по формуле

(Б.39)

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

Таблица Б.2 - Стандартные неопределенности входных величин и коэффициенты чувствительности

Входная величина

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

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

Коэффициент чувствительности ci

Вклад в utr (ciui)

U1

0,07072

2,94Е-05

4,50Е+00

1,32Е-04

t′1

577,5

1,73Е-04

1,18Е-03

2,04Е-07

t1

307,2

6,31Е-02

-1,18Е-03

-7,45Е-05

U2

1,50612

3,61Е-05

6,61Е+01

2,38E-03

t2

3971,2

1,15Е-03

2,51Е-02

2,89Е-05

U3

0,05937

9,74Е-07

4,63Е+00

4,51Е-06

t3

278,5

9,40Е-02

9,88Е-04

9,29Е-05

m

0,305450

8,66Е-06

3,27Е+02

2,83Е-03

R

9,99977

5,77Е-04

-1,00Е+01

-5,77Е-03

M

294,1846

1,40Е-03

3,40Е-01

4,76Е-04

F

96485,3383

4,79Е-03

-1,04Е-03

-4,96Е-06

Стандартная неопределенность от прецизионности uprec

0,004

Стандартная неопределенность от используемых реперов utr

0,007

Суммарная стандартная неопределенность установления аттестованного значения СО

0,008

Расширенная неопределенность, U (k = 2)

0,016

Библиография

1. International - Vocabulary of Basic and General Terms in Metrology. ISO Geneva, 1993.

2. Руководство по выражению неопределенности измерения. Перевод с английского под редакцией В.А. Слаева. - ВНИИМ. - С-Пб, 1999.

3. EURACHEM/CITAC Guide Traceability in Chemical Measurement. A guide to achieving comparable results in chemical measurement (Руководство ЕВPAXИМ/СИTАК Прослеживаемость в химических измерениях. Руководство по достижению сопоставимых результатов химического анализа. Пер. с англ./под редакцией Л.А. Конопелько. - ВНИИМ. - С-Пб, 2005).

4. ISO Guide 30 Terms and definitions used in connection with reference materials.

5. ISO Guide 34 General requirements for the competence of reference material producers.

6. ISO Guide 35 Reference Materials - General and statistical principles for certification.

7. EURACHEM/CITAC. Guide Quantifying Uncertainty in Analytical Measurement; Количественное описание неопределенности в аналитических измерениях. Руководство ЕВРАXИМ/СИТАК/Пер. с англ., под ред. Л.А. Конопелько. - СПб: ВНИИМ, 2002.

8. Налобин Д.П., Осинцева Е.В. Способы установления прослеживаемости аттестованных значений стандартных образцов (продолжение).//Стандартные образцы. 2008 №3.С 10-15.

9. «Atomic weights of the elements 1995. (Technical Report)», 1995, IUPAC.

10. Peter J. Mohr, and Barry N. Taylor, CODATA Recommended Values of the Fundamental Physical Constants: 2002, Rev. Mod. Phys. vol. 77 (1) 1-107. 2005.