Слайд 2
Метод динамики средних
Поставим задачу следующего вида:
N элементов и
M состояний.
Пусть элементы системы блуждают по своим состояниям: k=1…..M.
Найти среднюю численность k-го состояния (мат. ожидание) и дисперсию k-го состояния.
Мk(t)=?
Dk(t)=?
Введём в рассмотрение случайную величину xk(t) – численность элементов, находящихся в k-ом состоянии.
Тогда x1(t)+x2(t)+…+xM(t)=N
Слайд 3
Метод динамики средних
Xk(i)(t) = 1, если i-ый элемент
находится в состоянии k;
Xk(i)(t) = 0, в противном случае.
Т.к.
элементы однородны, состояния независимы, то случайная величина xk(t) имеет биномиальное распределение (распределение Бернулли).
Таким образом, получаем связь между вероятностью k- го состояния элемента в произвольный момент времени с математическим ожиданием численности этих состояний по всем элементам.
Т.е. достаточно определить вероятности Pk(t) для одного элемента.
Слайд 4
Метод динамики средних
Алгоритм:
Описать состояния одного элемента системы.
Нарисовать граф
состояний для одного элемента.
Получить для каждого состояния его среднюю
численность.
Составить дифференциальные уравнения динамики средних. При этом каждый член в уравнении равен произведению интенсивности данного перехода на среднюю численность состояния до перехода. Исходящий поток имеет знак минус, входящий – плюс.
Решить систему дифференциальных уравнений относительно средних численностей.
Вычислить значения дисперсий и средних квадратических отклонений.
Слайд 5
Метод динамики средних
S0
S2
S3
S1
Слайд 6
Системы массового обслуживания
Слайд 7
СМО. Основные понятия
Под СМО понимают системы, на вход
которых подается случайный поток однотипных заявок (событий), которые обрабатывают
однотипными каналами или устройствами.
Условно систему массового обслуживания можно представить следующим образом:
Очередь на обслуживание заявки
Обслуженные заявки
Слайд 8
СМО. Основные понятия
Каналом обслуживания называется техническое или какое-либо
другое устройство, обеспечивающее выполнение какой-либо заявки.
В теории массового обслуживания
работа канала характеризуется тем временем, которое затрачивается на обслуживание одной заявки.
В общем случае оно считается случайным и часто время обслуживания в теории массового обслуживания принимают подчиняющимся пуассоновскому закону распределения.
Тогда на выходе канала обслуживания будет появляться обслуженные заявки с интервалом времени, подчиняющимся пуассоновскому закону распределения.
Также, если на входе системы существует очередь и каждая заявка обслуживается независимо от других и не покидает систему во время обслуживания, а только после его завершения, то интервалы времени между обслуженными заявками будут также независимы и подчинены пуассоновскому закону распределения.
Слайд 9
Структура СМО
Структура системы задается:
потоком заявок,
количеством приборов (каналов) в
системе,
длительностями обслуживания,
числом мест ожидания,
в сложных системах – алгоритмом прохождения
заявки по различным каналам системы.
Слайд 10
Структура СМО
Однородный поток заявок полностью описывается случайными моментами
их поступления в СМО, то есть ему соответствует случайный
процесс.
В случае параллельных каналов и обслуживания заявок на любом свободном канале, СМО является многоканальной. При единственном канале – одноканальной.
Если обслуживание идет последовательно на разных приборах, то выделяют:
многофазные СМО, если одновременно на разных приборах могут обслуживаться несколько заявок;
многоэтапные СМО, если одновременно возможно обслуживание только на одном приборе.
Слайд 11
Дисциплина обслуживания
Дисциплина обслуживания – это порядок распределения заявок
между свободными каналами, кроме того это закон образования очереди,
закон поведения заявок в очереди и порядок обслуживания заявок.
Различают следующие варианты порядка обслуживания:
прямой;
инверсионный;
случайный;
с разделением времени;
с разделением процессора;
пакетный
и др.
Слайд 12
Порядок обслуживания
Прямой порядок обслуживания – обслуживание в порядке
поступления – FIFO (First In First Out).
Инверсионный или стековый
порядок обслуживания – обслуживание заявки, поступившей в систему последней – LIFO (Last In First Out).
Случайный порядок обслуживания – любая заявка, из имеющихся в очереди, может быть выбрана на обслуживание с равной вероятностью.
Режим разделяемого времени – вся длительность обслуживания разбивается на этапы, после завершания некоторого этапа тем или иным образом выбирается одна из заявок для обслуживания её очередного этапа.
Слайд 13
Порядок обслуживания
Режим равномерного разделения процессора – процессор (т.е.
канал) может одновременно обслуживать несколько заявок. Но при этом
интенсивность (т.е. скорость) обслуживания каждой заявки уменьшается во столько же раз, сколько одновременно заявок обрабатывается процессором.
Пакетное обслуживание. Задается некоторое правило формирования пакетов – например, в пакет, который поступает на обслуживание, включают все заявки, которые находились в системе на момент освобождения канала от предыдущего пакета.
Слайд 14
Порядок обслуживания
Кроме того, заявкам, тем или иным способом,
может назначаться приоритет. Это обычно требуется делать при наличии
разнородных требований в системе.
Приоритетные правила приписывают каждому типу заявок соответствующее значение приоритета. Способы задания могут быть следующие:
относительный;
абсолютный;
чередующийся;
в порядке возрастания длительности обслуживания;
в порядке убывания длительности обслуживания
и др.
Слайд 15
Порядок обслуживания
Относительный приоритет – при завершении обслуживания очередной
заявки, из очереди следующей выбирается заявка с наивысшим приоритетом.
Абсолютный
приоритет – при поступлении в СМО заявки с более высоким приоритетом, чем заявка, которая обслуживается в данный момент, происходит прерывание обслуживания. Поступившая заявка начинает обслуживаться прибором(каналом), а прерванная заявка:
может быть возвращенной в очередь с последующим дообслуживанием;
может быть возвращенной в очередь и затем быть обслуженой заново;
может быть отброшена.
Слайд 16
Порядок обслуживания
Чередующийся приоритет. За всеми заявками того же
типа, что и обслуживающаяся в данный момент, закрепляется наивысший
приоритет. После того, как все такие заявки будут обслужены, из очереди некоторым образом выбирается следующее требование, например по правилу относительных приоритетов.
При обслуживании в первую очередь требований с наименьшим временем обслуживания, можно получить оптимальное обслуживание, а именно, - минимальное среднее время пребывания в системе требований при отсутствии задержек и потерь из-за переориентации и изменения порядка обслуживания. Однако, чаще всего неизвестно заранее, какова будет длительность обслуживания.
Слайд 17
Характериситки СМО
Работа СМО оценивается с помощью параметров эффективности:
среднее
число заявок, обслуженных в единицу времени;
среднее число занятых каналов;
средняя
длина очереди;
среднее время ожидания в очереди;
среднее число заявок, получивших отказ в обслуживании и т.д.
Также для анализа могут быть интересны следующие вероятности:
вероятность застать систему в свободном состоянии;
вероятность потери заявки из-за занятости системы;
вероятность ожидания заявкой начала обслуживания в очереди при поступлении в систему;
вероятность нахождения в систему определенного количества заявок.
Слайд 18
Задачи СМО
Можно выделить несколько основных классов задач:
выяснение условий
существования стационарных распределений искомых характеристик;
нахождение значений определенных вероятностей, распределений
вероятностей и их числовых характеристик.
оценка характеристик СМО при помощи реальных наблюдений или имитационного моделирования;
сравнение и оптимизация СМО по различным критериям за счет подбора значений числовых параметров системы, порядка обслуживания, структуры и т.д.
исследование асимптотического поведения характеристик СМО при стремлении каких-либо парметров к некоторым граничным значениям.
Слайд 19
Замкнутые и разомкнутые СМО
Разомкнутые системы: на вход поступает
извне некоторый поток заявок, причем источник заявок в
систему не входит и анализу не подвергается. Разомкнутые системы массового обслуживания бывают с отказами и без отказов: если система свободная, она обслуживается, если нет, то получает отказ на обслуживание.
В замкнутой системе число источников заявок ограничено и интенсивность поступления заявок зависит от состояния источников, обусловленных работой самой системы. При поступлении заявки на обслуживание заявка либо обслуживается, либо, если система занята, уходит из очереди. Очередь может быть ограниченная и не ограниченная. В замкнутой СМО заявка рано или поздно обслуживается, то есть никогда не получает отказа.
Слайд 20
Разомкнутые СМО
На вход поступает поток заявок на обслуживание
с интенсивностью λ.
На выходе:
λ0 - интенсивность потока обслуженных заявок;
λн
- интенсивность потока необслуженных заявок.
λ0
λн
Слайд 21
Одноканальная разомкнутая СМО
S0
S1
Слайд 22
Многоканальная СМО с отказами
все каналы однотипны;
время обслуживания заявки
случайно и образует простейший поток интенсивностью μ;
поступающий поток заявок
интенсивностью λ будем считать простейшим;
если все каналы заняты, то очередная заявка не будет удовлетворена.
Слайд 23
Многоканальная СМО с отказами
Уравнения Колмогорова для финальных вероятностей:
Слайд 24
Многоканальная СМО с отказами
Приведенная интенсивность заявок - среднее
число поступивших заявок за среднее время обслуживания одной заявки.
α
= λ/μ
Тогда, решая уравнения, получим формулы Эрланга:
Слайд 25
Многоканальная СМО с отказами
Если все каналы заняты –
заявка получит отказ. Поэтому вероятность отказа:
Вероятность обслуживания заявки:
Абсолютная пропускная
способность – среднее число заявок, обслуженных в единицу времени:
Среднее число загруженных каналов:
Коэффициент загрузки одного канала:
Слайд 26
Многоканальная СМО с бесконечной очередью
Допущения:
все каналы однотипны;
время обслуживания
заявки случайно и образует простейший поток интенсивностью μ;
поступающий поток
заявок интенсивностью λ будем считать простейшим;
если все каналы заняты, то очередная заявка ставится в конец очереди;
система может иметь бесконечное число состояний:
S0 – все каналы свободны;
S1 – занят один канал;
Sn – заняты все каналы;
Sn+m – заняты все каналы и m заявок в очереди;
и т.д.
Слайд 27
Многоканальная СМО с бесконечной очередью
Преобразуя уравнения Колмогорова для
конечных вероятностей аналогично СМО с отказами получим:
Слайд 28
Многоканальная СМО с бесконечной очередью
Если α>=n, то длина
очереди стремится к бесконечности.
Среднее число загруженных каналов:
Коэффициент загрузки одного
канала:
Среднее число заявок в очереди:
Среднее время ожидания в очереди:
Слайд 29
Замкнутая одноканальная СМО
S0
S1
S2
Sk
Sn
Слайд 31
Имитационное моделирование
Ситуации, в которых рекомендуется применять имитационное моделирование:
не
существует законченной постановки задачи исследования;
аналитические методы являются математически трудоемкими
и трудно реализуемыми;
одной из задач является наблюдение за поведением отдельных составляющих системы в течение некоторого периода времени;
процесс или явление невозможно наблюдать в реальной обстановке;
Слайд 32
Имитационное моделирование
необходимо управлять протеканием процессов в системе замедляя
или ускоряя их в процессе моделирования (имитации);
осуществляется подготовка специалистов
и освоение новой техники;
при проверке гипотез и изучении новых ситуаций в сложных системах, о которых мало что известно;
основная задача моделирования – изучение последовательности событий, возникающих в проектируемой системе, поиск узких мест и оптимизация проекта.
Слайд 33
Имитационное моделирование
Суть имитационного подхода – процесс функционирования сложной
системы представляют в виде алгоритма, реализуемого на компьютере.
Как и
при аналитическом подходе, при постановке задачи выделяются гипотезы о функционировании как всей системы, так и отдельных элементов.
Имитационный подход позволяет максимально использовать всю известную информацию о системе.
Слайд 34
Имитационное моделирование
Моделируется не только система (структура), но и
время.
Обычно используют один из трех следующих представлений времени:
реальное время
моделируемой системы;
модельное время, по которому обеспечивается синхронизация событий в системе;
машинное время, соответствующее затратам ресурсов времени компьютера в процессе моделирования.
Слайд 35
Имитационное моделирование
Построение имитационной модели (имитатора) можно разбить на
этапы, похожие на этапы построения математической модели вообще, но
с некоторыми отличиями.
Этапы построения имитационной модели:
содержательное описание объекта моделирования;
концептуальное описание объекта моделирования;
формальное описание объекта моделирования;
Слайд 36
Имитационное моделирование
описание имитатора;
программирование и отладка имитационной модели;
испытание и
исследование модели;
применение модели;
анализ результатов.
Слайд 37
Имитационное моделирование
Содержательное описание объекта моделирования
Формулируются основные цели моделирования,
определяется объект имитации, границы и ограничения модели, выбираются критерии
эффективности для сравнения различных вариантов системы
Слайд 38
Имитационное моделирование
Концептуальное описание объекта моделирования
Выдвигаются основные гипотезы и
допущения относительно модели. Формулируется упрощенное алгоритмическое описание объекта моделирования.
Проводится декомпозиция системы.
Таким образом имеется:
описание объекта моделирования;
параметры и переменные моделирования;
критерии эффективности;
список методов обработки и способов представления результатов моделирования.
Слайд 39
Имитационное моделирование
Формальное описание объекта моделирования
Формальное представление поведения компонентов
системы и их взаимодействия с помощью аппроксимации функциональными ззависимости,
алгоритмического описания или смешанного подхода.
Производится предварительная проверка построенного формального описания на соответствие объекту.
Выбираются средства программирования и аппаратное обеспечение.
Слайд 40
Имитационное моделирование
Описание имитатора
Перевод формального описания в описание имитатора.
Решаются вопросы о моделировании времени и об основных моментах
функционирования имитатора – начальные условия, условия остановки, сбор и сохранение логов и промежуточных данных, планирование.
Слайд 41
Имитационное моделирование
Программирование и отладка имитационной модели
Техническое задание на
ПО. Проектирование ПО, кодирование, тестирование и отладка.
Создание технической документации
– руководства, средства помощи и справки и т.д.
Слайд 42
Имитационное моделирование
Испытание и исследование модели
Проверка правильности алгоритма в
ходе работы программы-имитатора, проверка адекватности. При необходимости – калибровка
модели.
Оценка точности и устойчивости результатов, чувствительности критериев эффективности.
Слайд 43
Имитационное моделирование
Применение модели
Составление плана экспериментов – максимум информации
при минимальных временных и вычислительных затратах.
Проведение экспериментов и сбор
результатов.
Анализ результатов
Анализ результатов с целью поддержки принятия решений, выработки рекомендаций по проектированию или модификации системы, или другого использования полученной информации.
Применение различных средств представление результатов – графики, диаграммы, средства компьютерной графики.
Слайд 44
Имитационное моделирование
Дискретно- событийное моделирование (discrete-event simulation (DES))
Агентное моделирование
(Agent-based modeling)
Системная динамика (System dynamics)
Слайд 45
Имитационное моделирование
При построении имитационной модели необходимо использовать языки
программирования. Часто, имитационные модели сложных технологических систем разрабатывают с
использованием универсальных языков программирования (C, C++, Java, C# и т.д.).
Альтернативой является использование методов основанных на специализированных языках и системах имитационного моделирования GPSS, SIMAN, SLAM, SIMSCRIPT, SIMULA, GASP.
Слайд 46
Имитационное моделирование
Специализированные языки моделирования описывают поведение системы от
события к событию, означающим начало или окончание технологической операции.
Основные преимущества использования специализированных языков компьютерной имитации:
упрощение и ускорение создания имитационных моделей;
автоматическое формирование определенных типов данных, необходимых в процессе имитационного моделирования;
без специального на то указания пользователя собирается множество статистических данных, описывающих поведение модели;
возможность конструирования сложных имитационных моделей пользователям, не являющимися профессиональными программистами;
Слайд 47
Имитационное моделирование
программы имитационных моделей на специализированных языках моделирования
близки к описаниям моделируемых систем на естественном языке (например,
при моделировании СМО последовательность обработки заявки будет иметь вид: занять прибор (команда SEIZE), задержать на время обслуживания (команда ADVANCE), освободить прибор (команда RELEASE));
направленность на моделирование систем без получения аналитических закономерностей процессов, так как техпроцесс отображается не системой уравнений, а взаимодействием отдельных динамических элементов во времени и пространстве.
Слайд 48
СМО
Многие специализированные системы и языки имитационного моделирования используют
в качестве основы Системы Массового Обслуживания.
Слайд 49
СМО
Пример имитатора одноканальной СМО с отказами.
Имеется источник заявок,
который характеризуется:
общим числом заявок за всё время наблюдения Nmax;
числом
заявок N(t), сгенирированных к моменту времени t;
законом распределения интервалов времени Δt между появлением заявок;
временем появления текущей заявки tn.
Слайд 50
СМО
Канал обработки характеризуется следующими параметрами:
состоянием S(t) в момент
времени t – свободен или занят;
числом Nw обслуженных заявок;
суммарным
временем Tw нахождения в занятом состоянии;
производительностью Δw, т.е. законом распределения времен выполнения заявок определенного типа;
временем Wn окончания обслуживания n-й заявки.
Слайд 51
СМО
Параметры эффективности.
Вероятность отказа можно оценить как частость отказов.
Pотк(t)
= Wотк(t) = (N(t)-Nw(t))/N(t)
Коэффициент загрузки канала:
Kз(t) = Tw(t)/t
Задается число
Nmax заявок, которые должны быть сгенерированы и статистические параметры этих заявок, имитирующие случайный процесс поступления в реальной системе.
Чем больше число заявок, тем больше точность оценок.
Также задаются статистические характеристики процесса обработки заявок, имитирующие случайный процесс обработки в реальной системе.
С увеличением времени наблюдения за системой, т.е. числа заявок, параметры эффективности будут стремится к предельным значениям.
Слайд 54
Системная динамика
Этот подход был разработан и предложен Джеем
Форрестером в конце 1950х как “исследование информационных обратных связей
в промышленной деятельности с целью показать как организационная структура, усиления (в политиках) и задержки (в принятии решений и действиях) взаимодействуют, влияя на успешность предприятия”.
Приложения системной динамики также включают социальные, урбанистические, экологические системы.
Слайд 55
Системная динамика
Основными компонентами модели системной динамики являются:
накопители (и
их уровни), stocks, (например, материальных объектов, знаний, людей, денег);
потоки,
flows, характеризующие изменение накопителей;
связи, через которые потоки и накопители влияют друг на друга;
вспомогательные переменные, для задания влияний параметров.
Слайд 56
Системная динамика
СД абстрагируется от отдельных объектов и событий
и предполагает “агрегатный” взгляд на процессы, концентрируясь на политиках,
этими процессами управляющих.
Моделирование с помощью методологии СД – это представление структуры и поведения системы как множества взаимодействующих положительных и отрицательных обратных связей и задержек.
Причинные связи показывают как процесс А влияет на процесс Б. При этом процесс Б тоже, возможно, влияет на процесс А, но через длинную цепочку причинно-следственных связей.
Только изучение динамики всей системы со всеми ее связями может привести к корректному пониманию процессов развития системы.
Слайд 57
Системная динамика
Классическая модель распространения нового продукта Франка Басса:
Потенциальные
клиенты (Potential Adopters) становятся клиентами (Adopters) со скоростью распространения
(Adoption Rate), которая зависит от рекламы и “устной рекламы”, т.е. общения клиентов с не-клиентами.
Влияние рекламы моделируется следущим образом:
Некий постоянный процент потенциальных клиентов (Potential Adopters * Advertising Effectiveness).
Все контактируют со всеми. Количество контактов человека в единицу времени обозначается как Contact Rate.
Если клиент общался с потенциальным клиентом, последний становится клиентом с вероятностью Adoption Fraction.
Таким образом, в единицу времени все клиенты обратят
Adopters * Contact Rate * Adoption Fraction * [Potential Adopters / (Potential Adopters + Adopters)] потенциальных клиентов в клиентов.
Слайд 59
Модель распространения нового продукта Франка Басса
d(Потенциальные клиенты)/dt =
-Продажи
d(Клиенты)/dt = Продажи
Продажи = Продажи из-за рекламы + Продажи
из-за устной рекламы
Продажи из-за рекламы = Эффективность рекламы * Потенциальные клиенты
Продажи из-за устной рекламы = Частота контактов * Эффективность устной рекламы * Потенциальные клиенты * Клиенты / (Потенциальные клиенты + Клиенты)
Слайд 62
Системная динамика
Математически, системно-динамическая модель – это система дифференциальных
уравнений.
Имитационное моделирование с использованием системной динамики можно представить как
последовательность шагов:
определить границ задачи (и системы);
выявить основные накопители и потоки, которые меняют уровни этих накопителей;
выявить источники информации, которые влияют на потоки;
выявить основные обратные связи;
нарисовать диаграмму, связывающую накопители, потоки и источники информации;
записать уравнения, которые определяют потоки;
оценить параметры и начальные условия (экспертной оценкой, статистической оценкой, используя внешние источники информации);
провести моделирование.
Слайд 64
Агентное моделирование
Исследование поведения децентрализованных агентов, которое определяет поведение
системы в целом. Моделирование снизу вверх.
В основе модели лежит
набор основных элементов, из взаимодействия которых рождается обобщенное поведение системы.
«Возникающее» поведение (emergent behavior) представляет собой результат взаимодействия элементов системы. Соответственно, в рамках данного подхода к моделированию возникает необходимость корректно отобразить механизм поведения и взаимодействия элементов системы – т.н. «агентов».
Агент представляет собой «некую сущность, которая обладает активностью, автономным поведением, может принимать решения в соответствии с некоторым набором правил, может взаимодействовать с окружением и другими агентами, а также может изменяться»
Слайд 65
Агентное моделирование
По одному из распространенных представлений агент должен
обладать следующими характеристиками :
Агент является «идентифицируемым», т.е. представляет собой
конечного индивидуума с набором определенных характеристик и правил, определяющих его поведение и правила принятия решений. Агент автономен и может независимо действовать и принимать решения по взаимодействию с другими агентами.
Агент находится в определенной среде, позволяющей ему взаимодействовать с другими агентами. Агент может коммуницировать с другими (контактировать при определенных условиях и отвечать на контакт).
Агент имеет определенную цель, влияющую на его поведение.
Агент гибок и обладает способностью самообучения с течением времени на основе собственного опыта. В ряде случаев агент может даже изменять правила поведения на основе полученного опыта.
Слайд 66
Агентное моделирование.
Допущения
Поведение (его ключевые аспекты) агентов может
быть описано.
Механизмы взаимодействия могут быть описаны.
Система может быть построена
снизу вверх.
Примеры
Люди, группы, организации.
Другие живые организмы с социальной структурой.
Механизмы, роботы и их системы.
Слайд 67
Агентное моделирование. Пример правил
Простое моделирование стаи.
Сплоченность – двигаться
к средней позиции ближайших соседей.
Разделение – двигаться, чтобы избежать
давки.
Ориентация - двигаться в среднем направлении ближайших соседей
http://www.red3d.com/cwr/boids/
Слайд 69
Агентное моделирование
Нет централизованного управления:
работой всей системы;
моделированием системы;
изменением состояний
модели/системы.
Агенты соединяются с соседними объектами по средствам топологии взаимодействия.
Существует
много различных видов топологий.
Топологии могут быть статическими и динамическими.
Слайд 70
Агентное моделирование. Топологии взаимодействия
Евклидово пространство (2д, 3д)
Дискретная сетка
Сетевая
структура
ГИС – геологическая информационная система
…
Слайд 71
Агентное моделирование
http://www.anylogic.ru/screentshots
Слайд 72
Агентное моделирование
http://www.anylogic.ru/screentshots
Слайд 73
Агентное моделирование
http://www.anylogic.ru/screentshots
Слайд 74
Когда использовать агентное моделирование?
Существует естественное представление компонентов системы
как агентов:
есть решения и поведения, которые могут быть определены
дискретно (с ограничениями);
адаптация и изменение поведения агентов являются важными;
необходимо обучение и динамическое стратегическое поведение агентов;
важны динамические связи между агентами, связи могут формироваться и исчезать;
образование структур из агентов, адаптация и обучение важны на уровне структур;
есть пространственная составляющая для поведения и взаимодействия агентов.
Может потребоваться масштабирование вверх до произвольного уровня.
Структурные изменения процесса – результат моделирования, а не входные данные.
Слайд 76
Метод Монте-Карло
Интеграл заменяется математическим ожиданием функции g(ξ) случайной
величины ξ, равномерно распределенной в произвольной области Ω n-мерного
пространства, причем для приближенной оценки математического ожидания используется усреднение по достаточно большой выборке значений функции g(ξk).
Слайд 77
Метод Монте-Карло
Принципиальную основу метода Монте-Карло составляет усиленный закон
чисел в форме Колмогорова:
Чтобы среднее арифметическое независимых реализаций случайной
величины ξ сходилось с вероятностью единица к ее математическому ожиданию, необходимо и достаточно, чтобы это математическое ожидание существовало.
Слайд 82
Статистический ансамбль
Ансамбль – много копий (может даже бесконечно
много) некоторой системы, каждая из которых представляет собой состояние,
в котором система может находится.
Т.е. по сути, статистический ансамбль – распределение вероятностей состояния системы.
Слайд 84
Метод Монте-Карло
Усреднение по ансамблю можно свести к усреднению
вдоль марковской цепи.
Задача состоит в организации такого случайного блуждания
в конфигурационном пространстве, при котором выполняются 2 условия:
для полученной марковской цепи должно существовать стационарное распределение вероятности;
вероятность появления различных состояний должна совпадать с функцией плотности вероятности для данного ансамбля.
Слайд 87
Метрополис-Гастингс
Slide by Anton S. Kaplanyan
Karlsruhe Institute of Technology
Слайд 88
Метрополис-Гастингс
Slide by Anton S. Kaplanyan
Karlsruhe Institute of Technology
Слайд 89
Метрополис-Гастингс
Slide by Anton S. Kaplanyan
Karlsruhe Institute of Technology
Слайд 90
Метрополис-Гастингс
Slide by Anton S. Kaplanyan
Karlsruhe Institute of Technology
Слайд 91
Метрополис-Гастингс
Slide by Anton S. Kaplanyan
Karlsruhe Institute of Technology
Слайд 92
Метрополис-Гастингс
Slide by Anton S. Kaplanyan
Karlsruhe Institute of Technology
Слайд 93
Метрополис-Гастингс
Slide by Anton S. Kaplanyan
Karlsruhe Institute of Technology
Слайд 94
Slide by Anton S. Kaplanyan
Karlsruhe Institute of Technology
Метрополис-Гастингс
Слайд 95
Метрополис-Гастингс
Slide by Anton S. Kaplanyan
Karlsruhe Institute of Technology
Слайд 96
Метод Монте-Карло. Примеры
Вычисление многомерных интегралов.
Моделирование всех возможных конформаций
цепных молекул (полимеров и др.)
Получение полной фазовой диаграммы системы,
моделирование фазовых переходов (классический пример – модель Изинга)
Слайд 97
Метод Монте-Карло в естественных науках
Квантовый метод Монте-Карло
Диффузный метод
Монте-Карло (DMC)
Вариационный метод Монте-Карло (VMC)
Метод Монте-Карло в статистической физике
Вычисление
многомерных интегралов для вычисления макроскопических свойств системы. Например для некоторого свойства A:
- вектор всех степеней свободы системы
Слайд 98
Метод Монте-Карло в статистической физике
Вместо того, чтобы моделировать
динамическое поведение системы, производится генерация всех возможных состояний в
соответствии со статистикой Больцмана.
Для определения нового состояния системы из текущего используется Марковская цепь.
Каждый проход (итерация) метода – переход в новое состояние.
