Количественная оценка неопределенности
Количественная оценка неопределенности ( UQ ) — это наука о количественной характеристике и оценке неопределенностей как в вычислительных, так и в реальных приложениях. Он пытается определить, насколько вероятны определенные результаты, если некоторые аспекты системы точно не известны. Примером может служить прогнозирование ускорения человеческого тела при лобовом столкновении с другим автомобилем: даже если бы скорость была точно известна, небольшие различия в изготовлении отдельных автомобилей, в том, насколько сильно был затянут каждый болт и т. д., приведет к различным результатам, которые можно предсказать только в статистическом смысле.
Многие проблемы в области естественных наук и техники также изобилуют источниками неопределенности. Компьютерные эксперименты по компьютерному моделированию являются наиболее распространенным подходом к изучению проблем количественной оценки неопределенности. [1] [2] [3] [4]
Источники
[ редактировать ]Неопределенность может присутствовать в математических моделях и экспериментальных измерениях в различных контекстах. Одним из способов классификации источников неопределенности является рассмотрение: [5]
- Параметр
- Это происходит из-за параметров модели, которые являются входными данными для компьютерной модели (математической модели), но чьи точные значения неизвестны экспериментаторам и не могут контролироваться в физических экспериментах, или чьи значения не могут быть точно определены статистическими методами . Некоторыми примерами этого являются локальное ускорение свободного падения в эксперименте с падающим объектом, различные свойства материалов в анализе методом конечных элементов для техники и неопределенность множителя в контексте оптимизации макроэкономической политики .
- Параметрический
- Это происходит из-за изменчивости входных переменных модели. Например, размеры заготовки в процессе производства могут не соответствовать заданным и указанным в инструкциях, что может привести к изменениям в ее характеристиках.
- Структурная неопределенность
- Это явление, также известное как неадекватность модели, предвзятость модели или несоответствие модели, возникает из-за отсутствия знаний физики, лежащей в основе проблемы. Это зависит от того, насколько точно математическая модель описывает истинную систему для реальной ситуации, учитывая тот факт, что модели почти всегда являются лишь приближениями к реальности. Одним из примеров является моделирование процесса падения объекта с использованием модели свободного падения; сама модель неточна, поскольку всегда существует трение воздуха. В этом случае, даже если в модели нет неизвестного параметра, все равно ожидается несоответствие между моделью и истинной физикой.
- алгоритмический
- Также известна как числовая неопределенность или дискретная неопределенность. Этот тип возникает из-за числовых ошибок и численных аппроксимаций при реализации компьютерной модели. Большинство моделей слишком сложны, чтобы их можно было точно решить. Например, метод конечных элементов или метод конечных разностей можно использовать для аппроксимации решения уравнения в частных производных (которое вносит числовые ошибки). Другими примерами являются численное интегрирование и усечение бесконечной суммы, которые являются необходимыми приближениями при численной реализации.
- Экспериментальный
- Эта ошибка, также известная как ошибка наблюдения, возникает из-за изменчивости экспериментальных измерений. Экспериментальная неопределенность неизбежна, и ее можно заметить, многократно повторяя измерения с использованием одних и тех же настроек для всех входных данных/переменных.
- Интерполяция
- Это происходит из-за отсутствия доступных данных, собранных в результате компьютерного моделирования и/или экспериментальных измерений. Для других входных настроек, для которых нет данных моделирования или экспериментальных измерений, необходимо интерполировать или экстраполировать, чтобы спрогнозировать соответствующие ответы.
Алеаторика и эпистемика
[ редактировать ]Неопределенность иногда подразделяют на две категории: [6] [7] широко распространен в медицинских целях. [8]
- Алеаторик
- Алеаторическая неопределенность также известна как стохастическая неопределенность и отражает неизвестные, которые различаются каждый раз, когда мы проводим один и тот же эксперимент. Например, один выстрел стрелы из механического лука, который точно повторяет каждый запуск (с одинаковым ускорением, высотой, направлением и конечной скоростью), не попадет в одну и ту же точку цели из-за случайных и сложных вибраций стрелы. знание которых не может быть определено в достаточной степени, чтобы исключить возникающий в результате разброс точек удара. Аргумент здесь, очевидно, заключается в определении слова «нельзя». Тот факт, что мы не можем проводить достаточные измерения с помощью наших доступных в настоящее время измерительных устройств, не обязательно исключает существование такой информации, которая переместила бы эту неопределенность в категорию ниже. Алеаторика происходит от латинского alea или кости, что означает азартную игру.
- Эпистемическая неопределенность
- Эпистемическая неопределенность также известна как систематическая неопределенность и связана с вещами, которые в принципе можно знать, но не делают на практике. Это может быть связано с неточными измерениями, с тем, что модель не учитывает определенные эффекты или с тем, что определенные данные были намеренно скрыты. Примером источника этой неопределенности может быть сопротивление в эксперименте, предназначенном для измерения ускорения силы тяжести вблизи земной поверхности. Обычно используемое гравитационное ускорение 9,8 м/с² игнорирует влияние сопротивления воздуха, но сопротивление воздуха для объекта можно измерить и включить в эксперимент, чтобы уменьшить возникающую неопределенность при расчете гравитационного ускорения.
- Комбинированное возникновение и взаимодействие алеаторической и эпистемической неопределенности
- Алеаторическая и эпистемическая неопределенность также могут возникать одновременно в одном термине. Например, когда экспериментальные параметры демонстрируют алеаторическую неопределенность, и эти экспериментальные параметры вводятся в компьютерное моделирование. Если затем для количественной оценки неопределенности суррогатная модель , например, гауссовский процесс или полиномиальное расширение хаоса , изучается из компьютерных экспериментов, этот суррогат демонстрирует эпистемическую неопределенность, которая зависит от алеаторической неопределенности экспериментальных параметров или взаимодействует с ней. [4] Такую неопределенность уже нельзя классифицировать исключительно как алеаторическую или эпистемическую, она представляет собой более общую неопределенность, связанную с выводами.
В реальных приложениях присутствуют оба вида неопределенностей. Количественная оценка неопределенности предполагает явное выражение обоих типов неопределенности по отдельности. Количественная оценка алеаторических неопределенностей может быть относительно простой, при этом традиционная (частотная) вероятность является самой базовой формой. такие методы, как метод Монте-Карло Часто используются . Распределение вероятностей может быть представлено его моментами (в гауссовском случае достаточно среднего и ковариации , хотя, вообще говоря, даже знание всех моментов сколь угодно высокого порядка все еще не определяет функцию распределения однозначно), или, в последнее время, как такие методы, как Кархунен-Лёве и полиномиальные разложения хаоса. Чтобы оценить эпистемическую неопределенность, предпринимаются усилия, чтобы понять (недостаток) знаний о системе, процессе или механизме. Эпистемическая неопределенность обычно понимается через призму байесовской вероятности , где вероятности интерпретируются как показатель того, насколько уверен рациональный человек в отношении конкретного утверждения.
Математическая перспектива
[ редактировать ]В математике неопределенность часто характеризуется с точки зрения распределения вероятностей . С этой точки зрения эпистемическая неопределенность означает отсутствие уверенности в том, каково соответствующее распределение вероятностей, а алеаторическая неопределенность означает отсутствие уверенности в том, какой случайная выборка, будет полученная из распределения вероятностей.
Типы проблем
[ редактировать ]Существует два основных типа проблем при количественной оценке неопределенности: одна — это прямое распространение неопределенности (когда различные источники неопределенности распространяются через модель для прогнозирования общей неопределенности в реакции системы), а другая — обратная оценка неопределенности модели. и неопределенность параметров (когда параметры модели калибруются одновременно с использованием тестовых данных). Исследования первой проблемы получили широкое распространение, и для ее решения было разработано большинство методов анализа неопределенности. С другой стороны, последняя проблема привлекает все большее внимание в сообществе инженеров-проектировщиков, поскольку количественная оценка неопределенности модели и последующие прогнозы истинного ответа(ов) системы представляют большой интерес при проектировании устойчивых систем.
Вперед
[ редактировать ]Распространение неопределенности — это количественная оценка неопределенностей в выходных данных системы, распространяемых из неопределенных входных данных. Основное внимание уделяется влиянию на результаты параметрической изменчивости, указанной в источниках неопределенности. Целями анализа распространения неопределенности могут быть:
- Чтобы оценить моменты низкого порядка выходных данных, т.е. среднее значение и дисперсию .
- Оценить надежность выходных данных. Это особенно полезно в проектировании надежности , где выходные данные системы обычно тесно связаны с ее производительностью.
- Оценить полное распределение вероятностей выходных данных. Это полезно в сценарии оптимизации полезности , когда для расчета полезности используется полное распределение.
Обратный
[ редактировать ]Учитывая некоторые экспериментальные измерения системы и некоторые результаты компьютерного моделирования на основе ее математической модели, количественная оценка обратной неопределенности оценивает несоответствие между экспериментом и математической моделью (что называется коррекцией смещения ) и оценивает значения неизвестных параметров в модели, если они существуют. любые (что называется калибровкой параметра или просто калибровкой ). Обычно это гораздо более сложная проблема, чем прямое распространение неопределенности; однако это имеет большое значение, поскольку обычно реализуется в процессе обновления модели. Существует несколько сценариев количественной оценки обратной неопределенности:
Только коррекция смещения
[ редактировать ]Коррекция смещения количественно определяет неадекватность модели , то есть несоответствие между экспериментом и математической моделью. Общая формула обновления модели для коррекции смещения:
где обозначает экспериментальные измерения как функцию нескольких входных переменных , обозначает реакцию компьютерной модели (математической модели), обозначает аддитивную функцию невязки (также известную как функция смещения), а обозначает экспериментальную неопределенность. Цель состоит в том, чтобы оценить функцию невязки , и в качестве побочного продукта полученная обновленная модель . Доверительный интервал прогноза предоставляется вместе с обновленной моделью в качестве количественной оценки неопределенности.
Только калибровка параметров
[ редактировать ]Калибровка параметров оценивает значения одного или нескольких неизвестных параметров в математической модели. Общая формулировка обновления модели для калибровки:
где обозначает реакцию компьютерной модели, которая зависит от нескольких неизвестных параметров модели. , и обозначает истинные значения неизвестных параметров в ходе экспериментов. Цель состоит в том, чтобы либо оценить , или придумать распределение вероятностей это включает в себя лучшее знание истинных значений параметров.
Коррекция смещения и калибровка параметров
[ редактировать ]Он рассматривает неточную модель с одним или несколькими неизвестными параметрами, а формулировка обновления модели объединяет их вместе:
Это наиболее полная формулировка обновления модели, включающая все возможные источники неопределенности, и ее решение требует наибольших усилий.
Выборочные методологии
[ редактировать ]Было проведено много исследований для решения проблем количественной оценки неопределенности, хотя большинство из них связано с распространением неопределенности. За последние одно-два десятилетия также был разработан ряд подходов к решению обратных задач количественного определения неопределенности, которые оказались полезными для большинства задач малого и среднего масштаба.
Прямое распространение
[ редактировать ]Существующие подходы к распространению неопределенности включают вероятностные подходы и невероятностные подходы. В основном существует шесть категорий вероятностных подходов к распространению неопределенности: [9]
- Методы, основанные на моделировании: моделирование Монте-Карло , выборка по важности , адаптивная выборка и т. д.
- Общие методы, основанные на суррогатах: при неинтрузивном подходе суррогатная модель изучается, чтобы заменить эксперимент или моделирование дешевым и быстрым приближением. Методы, основанные на суррогатах, также могут использоваться полностью байесовским способом. [10] [4] [11] [12] Этот подход оказался особенно эффективным, когда стоимость выборки, например, дорогостоящего моделирования, непомерно высока.
- Методы, основанные на локальном разложении: ряд Тейлора , метод возмущений и т. д. Эти методы имеют преимущества при работе с относительно небольшой входной изменчивостью и выходными данными, которые не выражают высокую нелинейность. Эти линейные или линеаризованные методы подробно описаны в статье Распространение неопределенности .
- Методы, основанные на функциональном расширении: расширение Неймана, ортогональное или расширение Карунена-Лоэва (KLE), с полиномиальным расширением хаоса (PCE) и вейвлет-разложением в качестве особых случаев.
- Методы на основе наиболее вероятной точки (MPP): метод надежности первого порядка (ФОРМ) и метод надежности второго порядка (СОРМ).
- Методы, основанные на численном интегрировании: полное факторное численное интегрирование (FFNI) и уменьшение размерности (DR).
Для невероятностных подходов интервальный анализ [13] Нечеткая теория , теория возможностей и теория доказательств являются одними из наиболее широко используемых.
Вероятностный подход считается наиболее строгим подходом к анализу неопределенности в инженерном проектировании из-за его соответствия теории анализа решений. Его краеугольным камнем является вычисление функций плотности вероятности для выборочной статистики. [14] Это можно строго выполнить для случайных величин, которые можно получить путем преобразования гауссовых переменных, что приводит к точным доверительным интервалам.
Обратная неопределенность
[ редактировать ]Частотный
[ редактировать ]В регрессионном анализе и наименьших квадратов задачах легко доступна стандартная ошибка , оценок параметров которую можно расширить до доверительного интервала .
Байесовский
[ редактировать ]существует несколько методологий количественной оценки обратной неопределенности В рамках байесовской структуры . Наиболее сложным направлением является решение задач как с коррекцией смещения, так и с калибровкой параметров. Проблемы таких проблем включают не только влияние неадекватности модели и неопределенности параметров, но также отсутствие данных как компьютерного моделирования, так и экспериментов. Распространенной ситуацией является то, что входные настройки не совпадают в экспериментах и симуляциях. Другая распространенная ситуация заключается в том, что параметры, полученные в результате экспериментов, вводятся в моделирование. Для дорогостоящего моделирования часто необходима суррогатная модель , например, гауссовский процесс или полиномиальное расширение хаоса , определяющая обратную задачу для поиска суррогатной модели, которая лучше всего аппроксимирует моделирование. [4]
Модульный подход
[ редактировать ]Подходом к количественной оценке обратной неопределенности является модульный байесовский подход. [5] [15] Модульный байесовский подход получил свое название от процедуры, состоящей из четырех модулей. Помимо имеющихся на данный момент данных, априорное распределение необходимо задать неизвестных параметров.
- Модуль 1: Моделирование гауссовского процесса для компьютерной модели
Чтобы решить проблему отсутствия результатов моделирования, компьютерная модель заменяется моделью гауссовского процесса (GP).
где
- размерность входных переменных, а — размерность неизвестных параметров. Пока заранее определен, , известные как гиперпараметры модели GP, необходимо оценивать с помощью оценки максимального правдоподобия (MLE) . Этот модуль можно рассматривать как обобщенный метод кригинга .
- Модуль 2: Моделирование гауссовского процесса для функции несоответствия
Аналогично первому модулю функция невязки заменяется на модель ГП.
где
Вместе с априорным распределением неизвестных параметров и данными компьютерных моделей и экспериментов можно получить оценки максимального правдоподобия для . В то же время, из Модуля 1 также обновляется.
- Модуль 3: Апостериорное распределение неизвестных параметров
Теорема Байеса применяется для расчета апостериорного распределения неизвестных параметров:
где включает все фиксированные гиперпараметры в предыдущих модулях.
- Модуль 4: Прогнозирование экспериментального ответа и функции несоответствия
Полный подход
[ редактировать ]Полностью байесовский подход требует, чтобы не только априорные значения неизвестных параметров но также априорные значения для других гиперпараметров должно быть назначено. Это следует за следующими шагами: [16]
- Выведите апостериорное распределение ;
- Интегрировать выйти и получить . Этот единственный шаг завершает калибровку;
- Прогнозирование экспериментального ответа и функции невязки.
Однако данный подход имеет существенные недостатки:
- В большинстве случаев является весьма трудноразрешимой функцией . Следовательно, интеграция становится очень проблематичной. Более того, если априорные значения для других гиперпараметров не выбраны тщательно, сложность численного интегрирования возрастает еще больше.
- На этапе прогнозирования прогноз (который должен, по крайней мере, включать ожидаемое значение реакций системы) также требует численного интегрирования. Цепь Маркова Монте-Карло (MCMC) часто используется для интегрирования; однако это вычислительно дорого.
Полностью байесовский подход требует огромного количества вычислений и пока может оказаться непригодным для решения самых сложных ситуаций моделирования. [16]
Известные проблемы
[ редактировать ]Теории и методологии распространения неопределенности гораздо лучше разработаны по сравнению с количественной обратным определением неопределенности. Для последнего остается нерешенным ряд трудностей:
- Проблема размерности: вычислительные затраты резко возрастают с увеличением размерности задачи, т. е. количества входных переменных и/или количества неизвестных параметров.
- Проблема идентификации: [17] Множественные комбинации неизвестных параметров и функции несоответствия могут дать один и тот же экспериментальный прогноз. Следовательно, разные значения параметров не могут быть различены/идентифицированы. Эту проблему можно обойти в байесовском подходе, где такие комбинации усредняются. [4]
- Неполный ответ модели: относится к модели, не имеющей решения для некоторых комбинаций входных переменных. [18] [19]
- Количественная оценка неопределенности входных величин: важные события отсутствуют в имеющихся данных или критические количества не идентифицированы аналитиками, например, из-за ограничений существующих моделей. [20]
- Мало внимания уделяется влиянию решений, сделанных аналитиками. [21]
См. также
[ редактировать ]- Компьютерный эксперимент
- Необходимы дальнейшие исследования
- Количественная оценка маржи и неопределенностей
- Вероятностные числа
- Байесовская регрессия
- Байесовская вероятность
Ссылки
[ редактировать ]- ^ Сакс, Джером; Уэлч, Уильям Дж.; Митчелл, Тоби Дж.; Винн, Генри П. (1989). «Планирование и анализ компьютерных экспериментов» . Статистическая наука . 4 (4): 409–423. дои : 10.1214/ss/1177012413 . JSTOR 2245858 .
- ^ Иман, Рональд Л.; Хелтон, Джон К. (1988). «Исследование методов анализа неопределенности и чувствительности компьютерных моделей». Анализ рисков . 8 (1). Уайли: 71–90. Бибкод : 1988РискА...8...71И . дои : 10.1111/j.1539-6924.1988.tb01155.x . ISSN 0272-4332 .
- ^ Уокер, МЫ; Харремоэс, П.; Ротманс Дж.; ван дер Слейс, JP; ван Асселт, MBA; Янссен, П.; Крайер фон Краусс, депутат парламента (2003 г.). «Определение неопределенности: концептуальная основа управления неопределенностью в поддержке принятия решений на основе моделей». Комплексная оценка . 4 (1). Издательство Swets & Zeitlinger: 5–17. Бибкод : 2003IntAs...4....5W . дои : 10.1076/iaij.4.1.5.16466 . hdl : 1874/386032 . ISSN 1389-5176 .
- ^ Jump up to: а б с д и Ранфтл, Саша; фон дер Линден, Вольфганг (13 ноября 2021 г.). «Байесовский суррогатный анализ и распространение неопределенности» . Форум физических наук . 3 (1): 6. arXiv : 2101.04038 . дои : 10.3390/psf2021003006 . ISSN 2673-9984 .
- ^ Jump up to: а б Кеннеди, Марк С.; О'Хаган, Энтони (2001). «Байесовская калибровка компьютерных моделей» . Журнал Королевского статистического общества, серия B (статистическая методология) . 63 (3): 425–464. дои : 10.1111/1467-9868.00294 .
- ^ Дер Кюрегян, Армен; Дитлевсен, Уве (2009). «Алеаторный или эпистемический? Имеет ли это значение?». Структурная безопасность . 31 (2): 105–112. doi : 10.1016/j.strusafe.2008.06.020 .
- ^ Маттис, Герман Г. (2007). «Количественная оценка неопределенности: современное вычислительное представление вероятности и приложения». Экстремальные техногенные и природные опасности в динамике сооружений . Серия «Безопасность НАТО через науку». стр. 105–135. дои : 10.1007/978-1-4020-5656-7_4 . ISBN 978-1-4020-5654-3 .
- ^ Абхая Индраян, Медицинская биостатистика , второе издание, Chapman & Hall/CRC Press, 2008, страницы 8, 673.
- ^ Ли, Ш.; Чен, В. (9 мая 2008 г.). «Сравнительное исследование методов распространения неопределенности для задач типа черного ящика». Структурная и междисциплинарная оптимизация . 37 (3). ООО «Спрингер Сайенс энд Бизнес Медиа»: 239–253. дои : 10.1007/s00158-008-0234-7 . ISSN 1615-147X . S2CID 119988015 .
- ^ Карденас, IC (2019). «Об использовании байесовских сетей в качестве подхода к метамоделированию для анализа неопределенностей при анализе устойчивости склонов». Геориск: оценка и управление рисками для инженерных систем и опасных геологических процессов . 13 (1): 53–65. Бибкод : 2019GAMRE..13...53C . дои : 10.1080/17499518.2018.1498524 . S2CID 216590427 .
- ^ Ранфтл, Саша; Мелито, Джан Марко; Бадели, Вахид; Рейнбахер-Кестингер, Алиса; Эллерманн, Катрин; фон дер Линден, Вольфганг (31 декабря 2019 г.). «Количественная оценка байесовской неопределенности с использованием данных разной точности и гауссовских процессов для импедансной кардиографии расслоения аорты» . Энтропия . 22 (1): 58. Бибкод : 2019Entrp..22...58R . дои : 10.3390/e22010058 . ISSN 1099-4300 . ПМЦ 7516489 . PMID 33285833 .
- ^ Ранфтл, Саша; Мюллер, Томас Стефан; Виндбергер, Урсула; фон дер Линден, Вольфганг; Бренн, Гюнтер (03 мая 2021 г.), Данные и коды для «Байесовского подхода к реологическим неопределенностям крови в аортальной гемодинамике» , номер документа : 10.5281/zenodo.5237189 , получено 12 января 2022 г.
- ^ Жолен, Л.; Киффер, М.; Дидрит, О.; Уолтер, Э. (2001). Прикладной интервальный анализ . Спрингер. ISBN 1-85233-219-0 .
- ^ Арнаут, Л.Р. Неопределенность измерений в реверберационных камерах - I. Примеры статистики. Технический отчет TQE 2, 2-й. ред., с. 3.1, Национальная физическая лаборатория, 2008.
- ^ Марк К. Кеннеди, Энтони О'Хаган, Дополнительные сведения о байесовской калибровке компьютерных моделей , Шеффилд, Университет Шеффилда: 1–13, 2000 г.
- ^ Jump up to: а б Баярри, MJ ; Бергер, Дж. О.; Лю, Ф. (01 марта 2009 г.). «Модуляризация в байесовском анализе с упором на анализ компьютерных моделей» . Байесовский анализ . 4 (1). Институт математической статистики: 119–150. дои : 10.1214/09-ba404 . ISSN 1936-0975 .
- ^ Арендт, Пол Д.; Апли, Дэниел В.; Чен, Вэй ; Лэмб, Дэвид; Горсич, Дэвид (28 сентября 2012 г.). «Улучшение идентификации при калибровке модели с использованием множественных ответов». Журнал механического дизайна . 134 (10). ASME International: 100909. doi : 10.1115/1.4007573 . ISSN 1050-0472 .
- ^ Карденас, IC (2019). «Об использовании байесовских сетей в качестве подхода к метамоделированию для анализа неопределенностей при анализе устойчивости склонов». Геориск: оценка и управление рисками для инженерных систем и опасных геологических процессов . 13 (1): 53–65. Бибкод : 2019GAMRE..13...53C . дои : 10.1080/17499518.2018.1498524 . S2CID 216590427 .
- ^ ван ден Эйнден, AP; Швекендик, Т; Хикс, Массачусетс (2021). «Метамоделирование для анализа геотехнической надежности с использованием зашумленных и неполных моделей» . Геориск: оценка и управление рисками для инженерных систем и опасных геологических процессов . 16 (3): 518–535. дои : 10.1080/17499518.2021.1952611 . S2CID 238819106 .
- ^ Карденас, И.; Авен, Т.; Флагж, Р. (2022). «Решение проблем количественной оценки неопределенности. Случай оценки геологических опасностей» . Геосци. Модель Дев. Обсуждать . 16 (6): 1601–1615. doi : 10.5194/gmd-16-1601-2023 . HDL : 11250/3105739 .
- ^ Карденас, И.; Авен, Т.; Флагж, Р. (2022). «Решение проблем количественной оценки неопределенности. Случай оценки геологических опасностей» . Геосци. Модель Дев. Обсуждать . 16 (6): 1601–1615. doi : 10.5194/gmd-16-1601-2023 . HDL : 11250/3105739 .