Статистическое моделирование систем на эвм

В практике моделирования систем информатики наиболее часто приходится иметь дело с объектами, которые в процессе своего функционирования содержат элементы стохастичности или подвергаются стохастическим воздействиям внешней среды. Поэтому основным методом получения результатов с помощью имитационных моделей таких стохастических систем является метод статистического моделирования на ЭВМ, использующий в качестве теоретической базы предельные теоремы теории вероятностен. Возможность получения пользователем модели результатов статистического моделирования сложных систем в условиях ограниченности машинных ресурсов существенно зависит от эффективности процедур генерации псевдослучайных последовательностей на ЭВМ, положенных в основу имитации воздействий на элементы моделируемой системы.

ОБЩАЯ ХАРАКТЕРИСТИКА МЕТОДА СТАТИСТИЧЕСКОГО МОДЕЛИРОВАНИЯ

На этапе исследования и проектирования систем при построении и реализации машинных моделей (аналитических и имитационных) широко используется метод статистических испытаний (Монте-Карло), который базируется на использовании случайных чисел, т. е. возможных значений некоторой случайной величины с заданным распределением вероятностей. Статистическое моделирование представляет собой метод получения с помощью ЭВМ статистических данных о процессах, происходящих в моделируемой системе. Для получения представляющих интерес оценок характеристик моделируемой системы Я с учетом воздействий внешней среды ? статистические данные обрабатываются и классифицируются с использованием методов математической статистики [10, 13, 18).

Сущность метода статистического моделирования. Таким образом, сущность метода статистического моделирования сводится к построению для процесса функционирования исследуемой системы 5 некоторого моделирующего алгоритма, имитирующего поведение и взаимодействие элементов системы с учетом случайных входных воздействий и воздействий внешней среды ?, и реализации этого алгоритма с использованием программно-технических средств ЭВМ.

Различают две области применения метода статистического моделирования: 1) для изучения стохастических систем; 2) для решения детерминированных задач. Основной идеей, которая используется для решения детерминированных задач методом статистического моделирования, является замена детерминированной задачи эквивалентной схемой некоторой стохастической системы, выходные характеристики последней совпадают с результатом решения детерминированной задачи. Естественно, что при такой замене вместо точного решения задачи получается приближенное решение и погрешность уменьшается с увеличением числа испытаний (реализаций моделирующего алгоритма) N.

В результате статистического моделирования системы 5 получается серия частных значений искомых величин или функций, статистическая обработка которых позволяет получить сведения о поведении реального объекта или процесса в произвольные моменты времени. Если количество реализаций N достаточно велико, то полученные результаты моделирования системы приобретают статистическую устойчивость и с достаточной точностью могут быть приняты в качестве оценок искомых характеристик процесса функционирования системы 5.

Теоретической основой метода статистического моделирования систем на ЭВМ являются предельные теоремы теории вероятностей [2, 13]. Множества случайных явлений (событий, величин) подчиняются определенным закономерностям, позволяющим не только прогнозировать их поведение, но и количественно оценить некоторые средние их характеристики, проявляющие определенную устойчивость. Характерные закономерности наблюдаются также в распределениях случайных величин, которые образуются при сложении множества воздействий. Выражением этих закономерностей и устойчивости средних показателей являются так называемые предельные теоремы теории вероятностей, часть из которых приводится ниже в пригодной для практического использования при статистическом моделировании формулировке. Принципиальное значение предельных теорем состоит в том, что они гарантируют высокое качество статистических оценок при весьма большом числе испытаний (реализаций) N Практически, приемлемые при статистическом моделировании количественные оценки характеристик систем часто могут быть получены уже при сравнительно небольших (при использования ЭВМ) N.

Неравенство Чебышева. Для неотрицательной функции? (?) случайной величины? и любого АГ> 0 выполняется неравенство.

В частности, если ^ (О = —^)² и К=к²о² (где х среднее арифметическое;

о — среднее квадратическое отклонение), то.

Теорема Бернулли. Если проводится N независимых испытаний, в каждом из которых некоторое событие А осуществляется с вероятностью р, то относительная частота появления события т/ЛГ при N-*00 сходятся по вероятности к р, т. с. при любом ?>0.

где т — число положительных исходов испытания.

Теорема Пуассона. Если проводится N независимых испытаний и вероятность осуществления события А в /-м испытании равна р_ь то относительная частота появления события т/И при ЛГ-*оо сходится по вероятности к среднему из вероятностей />" т. е. при любом I > О.

Теорема Чебышева. Если в N независимых испытаниях наблюдаются значения …, х" случайной величины то при N-*0о среднее арифметическое значений случайной величины сходится по вероятности к ее математическому ожиданию а, т. е. при любом е>0.

Обобщенная теорема Чебышева. Если …" — независимые случайные величины с математическими ожиданиями а₁,…, а" и дисперсиями а, о%, ограниченными сверху одним и тем же чистом, то при М-*сс среднее арифметическое значений случайной величины сходится по вероятности к среднему арифметическому их математических ожиданий:

Теорема Маркова. Выражение (4.6) справедливо и для зависимых случайных величин {_1э…, если только.

Совокупность теорем, устанавливающих устойчивость средних показателей, принято называть законом больших чисел.

Центральная предельная теорема. Если — независимые одинаково рас пределенные случайные величины, имеющие математическое ожидание а и дисперсию о², то при ЛГ-*оо закон распределения суммы? х, неограниченно приближается к нормальному:

Здесь интеграл вероятностей.

Теорема Лапласа. Если в каждом из N независимых испытаний событие А появляется с вероятностью р. то.

где т — число появлений события А в N испытаниях. Теорема Лапласа является частным случаем центральной предельной теоремы.

Примеры статистического моделирования. Статистическое моделирование систем на ЭВМ требует формирования значений случайных величин, что реализуется с помощью датчиков (генераторов) случайных чисел. Не останавливаясь пока на способах их реализации для целей моделирования на ЭВМ, поясним сущность метода статистического моделирования следующими примерами.

Пример 4.1. Необходимо методом статистического моделирования найти оценки выходных характеристик некоторой стохастической системы ?_я, функционирование которой описывается следующими соотношениями: х = 1 —е"*— входное воздействие, е~^ф — воздействие внешней среды, где Я и <�р — случайные величины, для которых известны их функции распределения. Целью моделирования является оценка математического ожидания М[у] величины у. Зависимость последней от входного воздействия х и воздействия внешней среды.

V имеет вид у = х²+г².

Рис. 4.1. Структурная схема системы.

В качестве оценки математического ожидания А/[у], как следует из приведенных теорем теории вероятностей, может выступать среднее арифметическое, вычисленное по формуле.

где уI — случайное значение величины у; N — число реализаций, необходимое для статистической устойчивости результатов.

Структурная схема системы 5* показана на рис. 4.1.

Здесь элементы выполняют следующие функции: вычисление В.:

возведение в квадрат К:

суммирование С:

извлечение квадратного корня И:

Схема алгоритма, реализующего метод статистического моделирования для оценки М [у] системы приведена на рис. 4.2. Здесь ЬЛ и /7 функции распределения случайных величин Л и q> N — заданное число реализаций; Isi — номер текущей реализации; LATsX (_tFI/ я <�р,;

ВИД ГЕН ВРМ[…] — процедуры ввода исходных данных, генерации псевдослучайных последовательностей и выдачи результатов моделирования соответственно.

Таким образом, данная модель позволяет получить методом статистического моделирования на ЭВМ статистическую оценку математического ожидания выходной характеристики М у) рассмотренной стохастической системы Точность и достоверность результатов взаимодействия в основном будут определяться числом реализаций N.

Пример 4.2. Необходимо методом статистического моделирования найти оценку площади фигуры (рис. 4.3), ограниченной осями координат, ординатой а=1 и кривой у —/(а); при этом для определенности предполагается, что 0 </(«).

Рис. 4.2. Схема моделирующего алгоритма системы.

Таким образом, данная задача является чисто детерминированной и ее аналитическое решение сводится к вычислению определенного интеграла, т. е. искомая площадь фигуры.

Для решения этой детерминированной задачи методом статистического моделирования необходимо предварительно построить адекватную по выходным характеристикам стохастическую систему оценки характеристик которой будут совпадать с искомыми в данной детерминированной задаче. Вариант структурной схемы такой системы показан на рис. 4.4, где элементы выполняют следующие функции:

Рис. 4.3. Геометрическая интерпретация оценки площади фигуры.

Рис. 4.4. Структурная схема системы 5_С

вычисление В₂: S=*h'IN.

Система S_D функционирует следующим образом: получается пара независимых случайных чисел интервала (0, 1), определяется координата точки (x_it х,+|), показанной на рис. 4.3, вычисляется ордината у,=/(х;) и проводится сравнение величин уi и причем если точка (х,•, x_f+i) попала в площадь фигуры (в том числе и на кривую f (x)), то исход испытания считается положительным 1 ив итоге можно получить статистическую оценку площади фигуры 5ф по заданному числу реализаций ЛЛ Логическая схема моделирующего алгоритма вероятностной системы S_D представлена на рис. 4.5. Здесь Кау*/(в) — заданная функция (табличная кривая); N —- заданное число реализаций; /а/ — номер текущей реализации; ХЫх& и Д7/ах,_{+ 1}; Н1шЬ•, 5=5; 5Я = А'"? А, — суммирующая ячейка.

Таким образом, построение некоторой стохастической системы Яр позволяет методом статистического моделирования получить оценки для детерминированной задачи.

Пример 4.3. Необходимо методом статистического моделирования решить следующую задачу. Проводится 5″ 10 независимых выстрелов по мишени, причем вероятность попадания при одном выстреле задана и равна р. Требуется оценить вероятность того, что число попаданий в мишень будет четным, т. е. 0, 2, 4, 6, 8, 10.

Данная задача является вероятностной, причем существует ее аналитическое решение:

Рис. 4.5. Схема моделирующего алгоритма системы 5₀

Рис. 4.6. Структурная схема системы г

В качестве объекта статистического моделирования можно рассмотреть следующую вероятностную систему, структура которой представлена на рис. 4.6, где элементы выполняют такие функции: анализ А.:

Выходным воздействием в данной системе является событие четного числа попаданий в мишень в серии из десяти выстрелов. В качестве оценки выходной характеристики необходимо при числе испытаний (серий выстрелов), равном найти вероятность четного числа попаданий:

Логическая схема алгоритма статистического моделирования для оценки искомой характеристики такой системы Р (у) приведена на рис. 4.7. Здесь Рвр — заданная вероятность попадания в мишень при одном выстреле; N — заданное число реализаций;

о данном моделирующем алгоритме посте ввода исходных данных в реализации операторов цикла происходит обращение к генератору случайных чисел, т. е. получаются значения X/ случайной величины, равномерно распределенной в интервале (0, 1). Вероятность попадания случайной величины в интервал (0, р), где 1, равна длине этого отрезка, т. е. Р{х/<�р} Поэтому при каждом моделировании выстрела полученное случайное число X/ сравнивается с заданной вероятностью р и при Х (<�р регистрируется «попадание в мишень», а в противном случае — «промах». Далее моделируются серии из десяти испытаний каждая, подсчитывается четное число «попаданий» в каждой серии и находится статистическая опенка искомой характеристики Р (у).

Таким образом, подход при использовании статистического моделирования независимо от природы объекта исследования (будет ли он детерминированным или стохастическим) является общим, причем при статистическом моделировании детерминированных систем (система Зн в примере 4.2) необходимо предварительно построить стохастическую систему, выходные характеристики которой позволяют оценить искомые.

Отметим, что во всех рассмотренных примерах не требуется запоминания всего множества генерируемых случайных чисел, используемых при статистическом моделировании системы 5. Запоминается только накопленная сумма исходов и общее число ре;

Рис. 4.7. Схема моделирующего алгоритма системы ализаций.

Это немаловажное обстоятельство вообще является характерным при реализации имитационных моделей методом статистического моделирования на ЭВМ.