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


ANCHEM.RU » Форумы » 1. Аналитический форум ...
  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
Журнал Масс-спектрометрия Журнал Масс-спектрометрия
Оригинальные научные статьи, обзоры, учебные и некоторые справочные материалы по всем разделам теории и практики масс-спектрометрии.
Алеся
Пользователь
Ранг: 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
  «« назад || далее »»

Ответ на тему


ААС, ИСП-АЭС, ИСП-МС - прямые поставки в 2022 году

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


ЖУРНАЛ ЛАБОРАТОРИИ ЛИТЕРАТУРА ОБОРУДОВАНИЕ РАБОТА КАЛЕНДАРЬ ФОРУМ

Copyright © 2002-2022
«Аналитика-Мир профессионалов»

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