Форумы ANCHEM.RU: руководство по неопределенностям

Российский химико-аналитический портал

химический анализ и аналитическая химия в фокусе внимания ::: портал химиков-аналитиков ::: выбор профессионалов

карта портала ::: расширенный поиск

ANCHEM.RU » Форумы » 1. Аналитический форум ...

Форум химиков-аналитиков, аналитическая химия и химический анализ.

руководство по неопределенностям >>>

  Ответов в этой теме: 142
  Страница: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
  «« назад || далее »»

[ Ответ на тему ]

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

Garry пишет:
Я уже давно рассчитываю неопределенности методом Монте-Карло. Это очень просто, хотя и звучит страшненько. Достаточно иметь обычный Эксель, малость желания и понимания суть происходящего.
Есть соотвтствующий документ применения этого метода, который рекомендуется как 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 и по методу Монте-Карло с помощью Экселя.
Мой мейл где-то валяется здесь на форуме.

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

ЭКСАН, НПО, ООО
Научно-производственное объединение «ЭКСАН» было создано в 2002 г. группой сотрудников лаборатории газового анализа НИИМТ в г. Ижевске. Производит оборудования для экспрессного анализа содержания азота, водорода, кислорода, углерода и серы в металлах и сплавах.

Алеся
Пользователь
Ранг: 21

19.02.2009 // 8:56:49

FeaR пишет:
Нашел. Всем спасибо за помощь. Буду разбираться.

Ой, а не могли бы и со мной поделиться? Очень ручки чешутся посчитать что-нибудь. aleka{coбaчkа}inbox.ru Буду признательна.

Алеся
Пользователь
Ранг: 21

19.02.2009 // 9:01:29

Garry пишет:

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

Уважаемый Garry! Я желающая, можно мне выслать?! aleka{coбaчkа}inbox.ru Большое спасибо.

FeaR
Пользователь
Ранг: 197

19.02.2009 // 9:22:46

Алеся пишет:

FeaR пишет:
Нашел. Всем спасибо за помощь. Буду разбираться.
Ой, а не могли бы и со мной поделиться? Очень ручки чешутся посчитать что-нибудь. aleka{coбaчkа}inbox.ru Буду признательна.

Выслал. Проверяйте почтовый ящик.

Алеся
Пользователь
Ранг: 21

19.02.2009 // 9:52:05

FeaR пишет:

Алеся пишет:

FeaR пишет:
Нашел. Всем спасибо за помощь. Буду разбираться.
Ой, а не могли бы и со мной поделиться? Очень ручки чешутся посчитать что-нибудь. aleka{coбaчkа}inbox.ru Буду признательна.
Выслал. Проверяйте почтовый ящик.

Еще раз спасибо, и я буду разбираться.

Garry
VIP Member
Ранг: 1076

19.02.2009 // 10:11:23

Редактировано 1 раз(а)

DIV007 пишет:

Спасибо за развитие темы, интересно посмотреть пример
Буду очень признателен, div(собачка)tutl.by
Заранее благодарю.

Коллега, какая то проблема с Вашим ящиком. У меня The Bat матерится "Письмо не отправлено. Сервер сообщает: non-local recipient verification failed", а так вроде бы всем желающим выслал. Если что-то непонятно пишите... по возможности буду отвечать.
Есть некто в славном городе Харьков - проф. Захаров Игорь Петрович. Вот на его лекциях и семинарах основан принцип программирования Экселя на предмет составления бюджета неопределенности. Можно, конечно используя VBA вообще свести все манипуляции, кроме занесения данных по входным величинам к нажатию на одну кнопку. Экспериментально я уже это сделал, правда не в файлах для примера, которые разослал, а для своих расчетов неопределенности МВИ определения состава газовых смесей для одних заказчиков этого самого МВИ. Любознательные это могут сделать самостоятельно.


Ответов в этой теме: 142 Страница: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 «« назад \|\| далее »»	Ответ на тему

ПОСЛЕДНИЕ НОВОСТИ ANCHEM.RU: [ Все новости ]

Размещение рекламы / Контакты