Российский химико-аналитический портал | химический анализ и аналитическая химия в фокусе внимания ::: портал химиков-аналитиков ::: выбор профессионалов |
|
ANCHEM.RU » Форумы » 1. Аналитический форум ... |
руководство по неопределенностям >>>
|
Garry VIP Member Ранг: 1076 |
18.02.2009 // 11:12:07
Редактировано 1 раз(а) Я уже давно рассчитываю неопределенности методом Монте-Карло. Это очень просто, хотя и звучит страшненько. Достаточно иметь обычный Эксель, малость желания и понимания суть происходящего. Есть соотвтствующий документ применения этого метода, который рекомендуется как Suplements к традиционному GUM. Обзывается он JCGM 101:2008. Evaluation of measurement – Supplement to the «Guide to the expression of uncertainty in measurement – propagation of distributions using a Monte-Carlo method». Я наверное, являюсь пропагандистом этого метода, поскольку знаю на своей шкуре насколько он облегчает все телодвижения по этим расчетам. Суть метода сводится к тому, что вам нужно "нарисовать" перед собой уравнение расчета конечной велины (выходной величины), т.е. той которую вы вписываете в протоколы измерений. Так, вот все компоненты этой формулы стоящие после "=" являются входными величинами, которые имеют свой тип распределения неопределенности. Например, масса, отмериваемый объем - это неопределенность типа "Б" и распределена она равномерно. Коэффициенты регрессии калибровочного уравнения, площадь хроматографического пика - неопределенности типа "А" (нормальное распределение). Их неопределенности в виде СКО арифметического среднего нужно оценить из реальных экспериментальных данных. Т.е. для оценивания площади нужно 10 раз подряд продозировать одну и ту же концентрацию стандарта и определить СКО площади пика. Для оценивания коэффициентов калибровки существуют формулы для получения оценок. Далее в Экселе составляете так называемый бюджет неопределенности, в котором вписываете в столбик входные величины, затем в столбик их неопределенности для "А" - СКО среднего арифметического, а для "Б" данные погрешности из паспотра на пипетку, колбу и т.д. деленное на корень из 3. Открываете следующий лист в той же книге Экселя где используя генератор случайных чисел (модуль анализа должен быть установлен в Экселе ! если его нет, то должны установить из офисного пакета) генерируете массивы данных для каждой из входных величин в соответствии с типом распределения и неопределенностями из первого листа. Эксель это позволяет делать элементарно. Количество генерируемых данных должно быть 10000. Далее после того как сгенерированы все массивы данных в следующем столбце прописываете формулу - уравнение расчета выходной величины и размножаете ее на 10000 ячеек. В этом столбце Эксель используя данные из столбцов входных величин формирует массив данных выходной величины в котором учтены все неопределенности входных величин. Далее копируете значения ! массива выходной величины в следующий столбец и ранжируете данные по возрастанию. Находите интерквантильный промежуток путем вычитания значения из ячейки 9750 и ячейки 250. Эта величина будет вашей расширенной неопределенностью. Стандартную суммарную неопределенность находят простым вычислением СКО массива выходной величины, а коэффициент охвата находят делением расширенной неопределенности на стандартную суммарную. Вот и все. Всем желающим могу выслать пример расчета и по GUM и по методу Монте-Карло с помощью Экселя. Мой мейл где-то валяется здесь на форуме. |
||
ANCHEM.RU Администрация Ранг: 246 |
|||
skaerta Пользователь Ранг: 79 |
18.02.2009 // 11:54:47
Если не затруднит kaerta(собачка)mail.ru Заранее благодарна |
||
FeaR Пользователь Ранг: 197 |
18.02.2009 // 12:27:45
d_phil_b{coбaчkа}list.ru, пожалуйста |
||
DIV007 Пользователь Ранг: 58 |
18.02.2009 // 16:31:58
Спасибо за развитие темы, интересно посмотреть пример Буду очень признателен, div(собачка)tutl.by Заранее благодарю. |
||
Alex75 Пользователь Ранг: 459 |
18.02.2009 // 22:13:37
To FeaR: Огромная благодарность. To Garry: Игорь Николаевич, если можно, то и со мной поделитесь методом расчета - dimmu2000{coбaчkа}mail.ru |
||
Каталог ANCHEM.RU Администрация Ранг: 246 |
|
||
Алеся Пользователь Ранг: 21 |
19.02.2009 // 8:56:49
Ой, а не могли бы и со мной поделиться? Очень ручки чешутся посчитать что-нибудь. aleka{coбaчkа}inbox.ru Буду признательна. |
||
Алеся Пользователь Ранг: 21 |
19.02.2009 // 9:01:29
Уважаемый Garry! Я желающая, можно мне выслать?! aleka{coбaчkа}inbox.ru Большое спасибо. |
||
FeaR Пользователь Ранг: 197 |
19.02.2009 // 9:22:46
Выслал. Проверяйте почтовый ящик. |
||
Алеся Пользователь Ранг: 21 |
19.02.2009 // 9:52:05
Еще раз спасибо, и я буду разбираться. |
||
Garry VIP Member Ранг: 1076 |
19.02.2009 // 10:11:23
Редактировано 1 раз(а) Коллега, какая то проблема с Вашим ящиком. У меня The Bat матерится "Письмо не отправлено. Сервер сообщает: non-local recipient verification failed", а так вроде бы всем желающим выслал. Если что-то непонятно пишите... по возможности буду отвечать. Есть некто в славном городе Харьков - проф. Захаров Игорь Петрович. Вот на его лекциях и семинарах основан принцип программирования Экселя на предмет составления бюджета неопределенности. Можно, конечно используя VBA вообще свести все манипуляции, кроме занесения данных по входным величинам к нажатию на одну кнопку. Экспериментально я уже это сделал, правда не в файлах для примера, которые разослал, а для своих расчетов неопределенности МВИ определения состава газовых смесей для одних заказчиков этого самого МВИ. Любознательные это могут сделать самостоятельно. |
|
||
Ответов в этой теме: 142
|
|
ЖУРНАЛ | ЛАБОРАТОРИИ | ЛИТЕРАТУРА | ОБОРУДОВАНИЕ | РАБОТА | КАЛЕНДАРЬ | ФОРУМ |
Copyright © 2002-2022 «Аналитика-Мир профессионалов» |
Размещение рекламы / Контакты |