8. Моделирование стохастических процессов методом статистических испытаний8.1. Основные положенияРанее мы познакомились с методами построения математических моделей, дающих возможность установить аналитическую (формульную) связь между заданными условиями операции и результатом (исходом) операции, характеризующимся одним или несколькими параметрами – показателями эффективности. Если в ход операции вмешиваются случайные факторы, то она представляет собой случайный процесс, а показатель эффективности – вероятность какого-то события или же математическое ожидание какой-то случайной величины. Иногда удается построить аналитическую модель случайного процесса (например, систему дифференциальных уравнений для вероятностей состояния или алгебраических уравнений для предельных вероятностей состояния) и связать заданные условия операции с ее исходом аналитическими зависимостями. Однако, это удается далеко не всегда - главным образом, в тех случаях, когда случайный процесс, протекающий в рассматриваемой системе, марковский или близок к марковскому. Для произвольных потоков событий, переводящих систему из состояния в состояние, аналитические решения получены только для отдельных частных случаев, а в общем случае удовлетворительных методов математического описания соответствующих процессов не существует. В тех случаях, когда построение аналитической модели явления по той или иной причине трудно осуществимо, применяется другой метод моделирования, известный под названием метода статистических испытаний или, иначе, метода Монте-Карло. Существо метода. Вместо того, чтобы описывать случайное явление с помощью аналитических зависимостей, производится «розыгрыш» - моделирование случайного явления с помощью некоторой процедуры, дающей случайный результат. Так же как в жизни, когда конкретное осуществление процесса складывается каждый раз по иному, так и в результате «розыгрыша» мы получаем один экземпляр – одну реализацию случайного явления. Произведя такой «розыгрыш» очень большое число раз, мы получим статистический материал – множество реализаций случайного явления – который можно обработать обычными методами математической статистики. Нередко такой прием оказывается проще, чем попытки построить аналитическую модель явления и исследовать зависимость между его параметрами на этой модели. Для сложных операций, в которых участвует большое число элементов (машин, систем, людей, коллективов) и в которых случайные факторы сложным образом взаимодействуют между собой, метод статистических испытаний, как правило, оказывается проще аналитического. В сущности, методом «розыгрыша» может быть решена любая вероятностная задача: однако оправданным он становится только в случае, когда процедура «розыгрыша» проще, а не сложнее применения аналитических методов. Пример 1. Решается задача: по некоторой цели производится 4 независимых выстрела, каждый из которых попадает в нее с р=0,5. Для поражения цели требуется не менее двух попаданий. Определить вероятность поражения цели.^ Аналитический способ. Вероятность поражения цели вычисляется через вероятность противоположного события. Вероятность непоражения цели = сумме вероятностей ни одного попадания и ровно одного попадания; вероятность ни одного попадания 0,54; вероятность одного попадания равна:, следовательно.Розыгрыш. Будем моделировать процедуру стрельбы с помощью другой, тоже случайной процедуры. Будем бросать 4 монеты: условимся - герб попадание, решка – промах. Цель поражена - не менее двух гербов. Розыгрыш в нашем случае – бросание четырех монет, результат этого опыта - «поражение или не поражение цели». Повторим такой опыт (бросание 4 монет) очень много раз подряд. Тогда согласно теореме Бернулли, частота поражения цели почти наверняка будет мало отличаться от вероятности этого события W; значит, если мы бросим четыре монеты большое число раз N, мы почти наверняка получим число близкое к W, т.е. к 0,688. В данном случае определение вероятности W розыгрышем было несравненно трудней, чем аналитическим расчетом. Пример 2. Анализ поведения дискретного объекта (дискретные входные и выходные переменные) - «задача о пьяном прохожем или задача о случайном блуждании». Прохожий решил прогуляться, стоя на углу улиц. Пусть вероятность того, что, достигнув очередного перекрестка, он пойдет на север, юг, восток и запад, одинакова. Какова вероятность того, что пройдя 10 кварталов, прохожий окажется не далее 2 кварталов от места, где он начал прогулку. Обозначим его местонахождение на каждом перекрестке двумерным вектором (x1, x2) («выход»), где x1 – направление с востока на запад и x2– направление с севера на юг. Каждое перемещение на один квартал к востоку (x1 + 1), а каждое перемещение на один квартал к западу (x1 – 1) (x1 - дискретная переменная). К северу x2 + 1, к югу x2 – 1. Начальное положение (0,0). Если в конце прогулки абсолютные значения х1 и х2 будут больше 2, то будем считать, что он ушел дальше двух кварталов в конце прогулки протяженностью в 10 кварталов. Т.к. вероятность движения нашего прохожего в любом из 4 направлений по условию одинакова и равна 0,25, то можно оценить его передвижение с помощью таблицы случайных чисел. Условимся, что если случайное число (СЧ) лежит в пределах от 0 до 24, пьяный пойдет на восток и мы увеличим х1 на 1; если от 25 до 49, то он пойдет на запад и х1-1; если от 50 до 74, он пойдет на север и x2 + 1; если от 75 до 99, то на юг и x2 – 1.Блок-схема поведения прохожего.Нужно провести достаточно большое число «машинных опытов», чтобы получить достоверный результат. Другими методами такую задачу решить практически невозможно. В литературе этот метод получил название метода имитационного моделирования, а также машинного, статистического, вероятностного, Монте-Карло или метода машинной имитации. В принципе имитационное моделирование можно осуществлять на широком спектре устройств, начиная с аналоговых ЭВМ и кончая листом бумаги с карандашом, однако, как правило, ориентируются на использование ЭВМ. Метод имитационного моделирования может рассматриваться как своеобразный экспериментальный метод. Отличие от обычного эксперимента заключается в том, что в качестве объекта экспериментирования выступает имитационная модель, реализованная в виде программы на ЭВМ. При таком экспериментировании с моделью (в отличие от «решения» модели при аналитическом, например, моделировании) могут быть применены статистические методы. Заметим, что методом статистических испытаний можно находить не только вероятности событий, но и средние значения (математические ожидания) случайных величин. При этом применяется закон больших чисел (теорема Чебышева). Согласно этой теореме, при большом числе опытов среднее арифметическое наблюденных значений случайной величины почти наверняка мало отличается от её математического ожидания. Аналогичным образом могут быть найдены не только математические ожидания, но и дисперсии интересующих нас случайных величин. Метод Монте-Карло есть метод математического моделирования случайных явлений, в которых сама случайность непосредственно включается в процесс моделирования и представляет собой его существенный элемент. Каждый раз, когда в ход операции вмешивается тот или другой случайный фактор, его влияние имитируется с помощью специально организованного «розыгрыша» или жребия. Таким образом, строится одна реализация случайного явления, представляющая собой как бы результат одного «опыта». При большом числе реализаций интересующие нас характеристики случайного явления (вероятности, математические ожидания) находится так же, как они находятся из опыта. Моделирование случайных явлений методом Монте-Карло имеет общие черты с процессом набора опыта отдельными людьми и человеческими коллективами. И тут, и там каждая отдельная реализация случайна; устойчивые закономерности обнаруживаются лишь при многократном наблюдении явления, при обширном опыте. Большое число реализаций, требующееся при применении метода Монте-Карло, делает его вообще громоздким и трудоемким.^ 8.2. Способы организации единичного жребия Основным элементом, из совокупности которых складывается монте-карловская модель, является одна случайная реализация моделируемого явления, например: один «обстрел» цели», один «день работы» транспорта, одна «эпидемия» и т.п. Реализация представляет собой как бы один случай осуществления моделируемого случайного явления (процесса) со всеми присущими ему случайностями. Она разыгрывается с помощью специально разработанной процедуры или алгоритма, в котором важную роль играет собственно «розыгрыш» или бросание жребия». Каждый раз, когда в ход моделируемого процесса вмешивается случайность, её влияние учитывается не расчетом, а бросанием жребия. Предположим, что в ходе моделируемого процесса наступил момент, когда его дальнейшее развитие (а значит и результат) зависит от того, появилось ли на данном этапе событие А или не появилось (например: произошло ли попадание в цель, обнаружен ли некоторый объект, исправна ли некоторая аппаратура и т.д). Тогда нужно «бросанием жребия» решить вопрос: появилось событие А или не появилось? Для этого нужно привести в действие некоторый случайный механизм розыгрыша (бросить игральную кость, несколько монет или выбрать число из таблицы случайных чисел) и условиться о том, какой результат жребия означает появление, а какой – непоявление события А). Ниже мы увидим, что розыгрыш всегда можно организовать так, чтобы событие А имело любую наперед заданную вероятность. Кроме событий, появляющихся случайным образом, на ход и исход операции могут так же влиять разные случайные величины (время, координаты и т.д.). С помощью жребия можно разыграть значения любой случайной величины или совокупность значений нескольких случайных величин. Условимся называть единичным жребием любой элементарный опыт, в котором решается один из вопросов:Произошло или не произошло событие А?Какое из возможных событий А1,А2,…Аkпроизошло?Какое значение приняла случайная величина Х?Какую совокупность значений приняла система случайных величин Х1,Х2,…Хk? Рассмотрим способы организации всех разновидностей единичного жребия. При любой организации жребия должен быть пущен в ход какой-то механизм случайного выбора. Механизмы могут быть самыми разнообразными, однако любой из них может быть заменен стандартным механизмом, позволяющим решить одну задачу: получить случайную величину, распределенную с постоянной плотностью от 0 до 1. Условимся для краткости называть такую случайную величину «случайное число от 0 до 1» и обозначать R.^ Появилось или нет событие А?Пусть вероятность события А равна p: Р(А)=р. Выберем с помощью стандартного механизма случайное число R и будем считать, что если оно меньше р, событие А произошло, если больше р – не произошло. Действительно: если R – случайное число от 0 до 1, то , где f(r)=1 при 01 или.^ Какое из нескольких возможных событий появилось? Пусть имеется полная группа несовместных событий: А1,А2,…Аkс вероятностями р1,р2,…рk. Т.к. события несовместны и образуют полную группу, то р1+р2+…+рk=1. Разделим весь интервал от 0 до 1 на k участков длиной р1,р2,…рk. Если случайное число R, выданное стандартным механизмом, попало, например, на участок р3, это означает, что появилось событие А3. ^ 3. Какое значение приняла случайная величина?Пусть нам требуется «разыграть» значение случайной величины Х, имеющей известный закон распределения. Случай, когда величина Х дискретна (т.е. имеет отдельные значения х1,х2,…хkс вероятностями р1,р2,…рk) рассматривать не будем, т.к. он сводится к предыдущему пункту 2. Рассмотрим случай, когда случайная величина Х непрерывна и имеет заданную непрерывную функцию распределения F(x) (рис.)Докажем следующее утверждение: если взять на оси ординат случайное число ^ R (от 0 до 1) и найти то значение Х, при котором F(x)=R, то полученная случайная величина Х будет иметь функцию распределения F(x). Действительно, возьмем случайную величину Х и найдем её функцию распределения, т.е. вероятность Р(Х , что и требовалось доказать.Т.о. розыгрыш значения случайной величины ^ Х с заданной функцией распределения F(x) сводится к следующей процедуре: получить случайное число R от 0 до 1 и в качестве значения Х взять X=F-1(R), где F-1 – функция, обратная по отношению к F. Пример 1. Случайная величина Х распределена по показательному закону с плотностью f(x)=le-lx (x>0). Построить процедуру единичного жребия для получения значения Х. По заданной плотности f(x) находим функцию распределения: (x>0).График F(x) дан на рисунке.Графически значение случайной величины ^ Х можно разыграть так: взять случайное число от 0 до 1 на оси ординат и найти соответствующее ему значение абсциссы Х. Это же можно сделать расчетом, если написать: R=1-e-lX и решать это уравнение относительно Х (т.е. найти обратную по отношению к F функцию). Имеем e-lX=1-R, -lX=ln(1-R)ÞX= ln(1-R). Эту формулу можно упростить: вспомним, что если R – случайное число от 0 до 1, то (1-R)- также случайное число от 0 до 1, поэтому можно взять X =lnR. Т.о. процедура розыгрыша сводится к следующему: взять случайное число от 0 до 1, прологарифмировать его при натуральном основании, изменить знак и разделить на l. Пример 2. Розыгрыш значения случайной величины, распределенной по нормальному закону (короче – «нормальной») с математическим ожиданием mxи средним квадратичным отклонением. Плотность распределения случайной величины Х имеет вид: . Удобнее применить не общее правило, а поступить иначе: перейти от Х к другой (т.н. «нормированной») случайной величины Z=, разыграть значение этой величины, а затем уже по ней найти ^ Х. Это удобно потому, что mz=0, и придется только один раз и навсегда найти обратную функцию. Легко показать, что значение нормальной случайной величины Х с характеристиками mx,разыгрывается по формуле:, где Ф-1 – функция, обратная функции Лапласа. Есть и другой способ, основанный на центральной предельной теореме теории вероятностей. Согласно этой теореме, при сложении достаточно большого числа независимых случайных величин, сравнимых по своим дисперсиям, получается случайная величина, распределенная приближенно по нормальному закону, причем этот закон тем ближе к нормальному, чем больше случайных величин складывается (для большинства прикладных задач достаточно складывать 6 случайных величин от 0 до 1). В результате получается следующая процедура.^ 4. Какую совокупность значений примет система случайных величин? Пусть имеется система случайных величин: Х1, Х2, …Хn с совместной плотностью распределения f(х1,х2,…,хn. Если случайные величины независимы, то f(х1,х2,…,хn)=f1(x1)f2(x2) …fn(xn) и розыгрыш совокупности значений системы х1, х2,…,хn сводится к тому, чтобы разыграть каждую из них в отдельности, т.е. организовать n единичных жребиев типа, описанного в п.3. Если случайные величины зависимы, то f(х1,х2,…,хn)=f1(x1)f(x2/x1)f(x3/x1x2)…, где каждая последующая плотность распределения берется условная, при условии, что предыдущие случайные величины приняли определенные значения. При розыгрыше последовательности значений случайных величин получается сначала значение х1 случайной величины Х1; это значение берется в качестве аргумента в условной плотности f(x2/x1); разыгрывается значение х2 случайной величины Х2, оба значения х1, х2 берутся в качестве аргументов в условной плотности f(x3/x1x2) и т.д.^ 8.3. Современное содержание терминов «имитация», «имитационная модель».До сих пор термины «имитация», «имитационное моделирование», «имитационная модель», «имитационный эксперимент» охарактеризовывались с разных точек зрения. Соберем воедино все эти характеризации и постараемся выработать целостное представление о содержании этих терминов. Итак, термины «имитация» и «имитационный эксперимент» появились сначала в теории вероятностей и математической статистике как способ вычисления статистических характеристик интересующих нас случайных величин посредством воспроизведения реализаций соответствующего случайного процесса с помощью его математической модели. Воспроизведение реализаций случайного процесса и есть то, что естественно называть имитационным экспериментом, поскольку реальные эксперименты с измерением интересующих нас случайных величин как бы заменяются их имитацией с помощью математической модели данного процесса. Вскоре после начала использования методов прикладной математики в управлении экономикой, планировании, исследовании операций, проектировании термины «имитация», «имитационный эксперимент» приобрели в этих областях смысл, не совпадающий с их первоначальной трактовкой. Этими терминами стали обозначать способ выбора рационального управления сложным процессом (рационального плана, рациональной конструкции проектируемого изделия), состоящий в следующем. Некоторым образом разрабатываются варианты управлений (планов, конструкций). Затем эти варианты сравниваются. Для этого при каждом таком варианте процесс (функционирование проектируемого изделия) воспроизводится с помощью его математической модели. Сравнение может происходить по некоторым формальным критериям, а может носить неформальный характер, причем чем сложнее используемая модель, чем больше она содержит реальных факторов, влияющих на принятие решений, тем более естественна неформальная оценка сравниваемых результатов. Математические модели, ориентированные на такое их использование, получили название имитационных, процесс их составления стал называться имитационным моделированием, а каждая акция воспроизведения процесса (функционирования проектируемого изделия) - имитационным экспериментом. Если изучаемый процесс является случайным (если процесс достаточно сложен, то почти неизбежно он случаен), для сравнения, о котором шла речь выше, необходимо выполнять то, что в теории вероятностей и математической статистике называется имитацией, т. е. вычислять статистические характеристики этого случайного процесса путем набора необходимого для этого количества реализаций (предполагается, что аналитическими средствами вычислить статистические характеристики нельзя) и именно эти статистические характеристики и сравнивать. Даже если дело обстоит именно таким образом, специалисты в области управления, планирования, проектирования, исследования операций, произнося слово «имитация», будут иметь в виду не способ вычисления характеристик случайных процессов путем набора статистики (для них это некоторая необходимая техническая деталь), а то, что альтернативные варианты управлений (планов, конструкций проектируемого объекта) являются внешними по отношению к модели процесса, задаются «извне» ее, а не являются ее продуктом. Они будут называть воспроизведения процесса имитацией и тогда, когда процесс детерминирован. Для них термин «имитация» несет смысловую нагрузку противопоставления термину «оптимизация», в то время как для специалистов в области теории вероятностей и математической статистики термин «имитация» несет оттенок противопоставления аналитическим методам расчета статистических характеристик случайного процесса. Сопоставление оптимизационных задач с реальным содержанием задач планирования, управления, проектирования приводило к попыткам улучшить модели, лежащие в основе оптимизационных задач, что влекло за собой их усложнение, появление вместо одного критерия оптимальности нескольких или же вообще отказ от оптимизации в рамках усложнившихся моделей и использование их в режиме вариантных расчетов с задаваемыми извне модели вариантами планов (управлений, конструкций проектируемого изделия). Поскольку в этих усложненных моделях присутствовали, как правило, случайные факторы, то получение обоснованных результатов требовало вычисления статистических характеристик, т. е. имитации в том смысле, в котором ее понимают в теории вероятностей и математической статистике. Очень скоро специалисты в области использования математических методов в планировании, управлении, исследовании операций, проектировании присвоили новое содержание терминам «имитация», «имитационная модель», «имитационный эксперимент» — то, о котором говорилось выше. Конечно, содержание этих гуманитарных терминов нельзя очертить достаточно четко: каждый специалист имеет право понимать под ними то, что ему заблагорассудится. Тем не менее, суммируя все то, о чем говорилось, сопоставляя различные трактовки этих терминов, даваемые специалистами в разных сферах, можно сказать то, что уже говорилось ранее. Именно, имитационная модель — это модель, обладающая качествами из следующего набора: «сложность» модели, наличие в ней случайных факторов, описание процесса, развивающегося во времени, невозможность получения результатов без ЭВМ, предназначенность модели для использования ее в режиме вариантных расчетов, т. е. для сравнения путем выполнения имитационных экспериментов, заданных заранее, «извне модели» вариантов планов, управлений, конструкций. Ни одно из перечисленных качеств не является обязательным для того, чтобы именовать модель имитационной. Например, специалист по оптимизации назовет даже очень простую детерминированную модель имитационной, если в ее рамках никакая оптимизационная задача не решается, а она используется в режиме вариантных расчетов для сравнения заранее сформированных альтернативных вариантов управлений. Тем не менее, перечисленные качества в совокупности дают целостное представление о том, какой смысл имеет понятие «имитационная модель» в современной русскоязычной научной литературе. В англоязычной литературе терминам «имитация», «имитационная модель», «имитационный эксперимент» приблизительно соответствует термин «simulation». Наиболее известна трактовка термина «имитация», которую дал Р.Шеннон, определив ее как «процесс конструирования модели реальной системы и постановки экспериментов на этой модели с целью понять поведение системы либо оценить (в рамках ограничений, накладываемых некоторым критерием или совокупностью критериев) различные стратегии, обеспечивающие функционирование данной системы». Аналогичное определение этому термину дает Т. Нейлор: «численный метод проведения на цифровых вычислительных машинах экспериментов с математическими моделями, описывающими поведение сложных систем в течение продолжительных периодов времени». Обе приведенные трактовки термина «simulation» укладываются в то содержание терминов «имитация», «имитационная модель», «имитационный эксперимент», которое было очерчено выше. Первым результатом «столкновения» оптимизационных задач с реальным содержанием задач планирования, управления, проектирования в реальных производственных структурах стало, как указывалось выше, усложнение моделей. В рамках усложненных моделей происходил отказ от «жесткой» оптимизации и появление в работе с моделью режима, который получил название имитационного. Именно, вместе с оптимальными управлениями (планами) определялись и близкие к ним, предоставлялась возможность выбрать критерий оптимизации из некоторого множества, изменить параметры, ограничения, назначить часть позиций плана (часть управлений), а оптимизацию провести по остальным и т. д. Описанные возможности потребовали интерактивного режима работы с ЭВМ и соответствующего программного обеспечения. Таким образом возникли интерактивные проблемно-ориентированные оптимизационные системы. Часто усложнение модели приводило к отказу от оптимизации в ее рамках потому, что на оптимизацию уже не хватало вычислительных ресурсов, а еще чаще потому, что включение в модель большого количества факторов, которые необходимо учитывать при реальном управлении, приводило к невозможности постановки в ее рамках оптимизационных задач на выбор управлений (планов). Такие модели стали использовать в режиме вариантных расчетов, как об этом говорилось выше. Это использование также потребовало интерактивного режима работы с ЭВМ, применения средств визуализации информации и вообще разнообразных средств манипулирования информацией. Так возникли проблемно-ориентированные имитационные системы. Очень скоро выяснилось, что «чистый» режим имитации при работе с имитационной моделью во многих случаях почти столь же непрактичен, как и «чистая» оптимизация, о чем много говорилось выше. Например, среди многочисленных позиций плана лишь несколько являются предметом особого интереса для эксперта, занимающегося составлением данного плана. Назначение всех остальных - рутинная работа. Если это назначение выполняется «вручную», то от использования ЭВМ для эксперта мало толку по сравнению с обычной технологией составления плана. Существенный выигрыш получается тогда, когда рутинная работа по назначению позиций плана, мало интересующих эксперта, автоматизируется на основе оптимизации. Позиции же плана действительно существенные «доверять» оптимизации нельзя. Очень полезно получить оптимальные значения этих позиций, однако лишь для анализа ситуации, как «информацию для размышления». Таким образом, в проблемно-ориентированных имитационных системах возникает оптимизационный режим работы. Он обеспечивается совокупностью упрощенных моделей изучаемого процесса вместе с алгоритмами, вычисляющими в рамках этих моделей оптимальные управления.^ 8.4. Приемы построения и эксплуатации имитационных моделейПри создании имитационных моделей в настоящее время используется два подхода: дискретный и непрерывный. Выбор подхода в значительной мере определяется свойствами объекта-оригинала и характером воздействия на него внешней среды. Метод статистического моделирования (метод Монте-Карло) – можно рассматривать как частный случай дискретных вероятностных имитационных моделей. При использовании дискретного подхода к созданию имитационных моделей обычно применяются абстрактные системы (математические схемы) трех основных типов: автоматные системы, системы массового обслуживания и агрегативные системы. В случае непрерывного подхода моделируемый объект независимо от его природы формализуется в виде непрерывной абстрактной системы, между элементами которой циркулируют потоки той или иной природы. Структура такой системы представляется графически в виде диаграммы (схемы) потоков. Основными элементами непрерывной системы рассматриваемого типа являются абстрактные «бункеры» (емкости, резервуары), а также элементы задержки. Под имитационной моделью объекта – оригинала, в общем случае, мы можем понимать определенную систему, состоящую из отдельных подсистем (элементов, компонентов) и связей между ними, причем функционирование (изменение состояний) и внутреннее изменение всех элементов модели под действием связей может быть алгоритмизировано тем или иным образом, так же как и взаимодействие системы с внешней средой. Тогда имитация функционирования системы сводится к пошаговому воспроизведению на ЭВМ процесса функционирования всех её элементов с учетом их взаимодействия и воздействий внешней среды. В имитационном моделировании могут быть алгоритмизированы и воспроизведены процессы функционирования и взаимодействия самых различных элементов абстрактной системы - дискретных и непрерывных, вероятностных и детерминированных, выполняющих функцию обслуживания, задержки и других. В качестве имитационной модели объекта при такой постановке выступает программа на ЭВМ (вместе с обслуживающими, сервисными программами), написанная либо на универсальных языках высокого уровня, либо с применением специализированных языков имитационного моделирования (GPSS/320 – для дискретных систем типа систем массового обслуживания, Q-GERT язык сетевого моделирования для дискретных систем, SIMPL/1, GASP IV и др.). До настоящего времени разработка и применение имитационных моделей все ещё в большей степени искусство, чем наука. Имитационное моделирование наиболее эффективно использовать на высоких уровнях иерархии, при рассмотрении взаимодействия во времени большого числа сложных объектов, что характерно для высшей степени технологии и производственных процессов. ^ Дискретные имитационные модели. Цель имитационного моделирования состоит в воспроизведении поведения исследуемой системы на основе результатов анализа наиболее существенных взаимосвязей между её элементами. Результаты исследования имитационной модели, как правило, представляют собой оценки функциональных характеристик той системы, поведение которой имитируется. Так, например, при имитационном моделировании любой СМО практический интерес могут представить такие показатели, как средняя продолжительность обслуживания заявки, средняя длина очереди, доля времени простоя и т.д. Первый шаг к созданию имитационной модели состоит в описании реально существующей системы с использованием характеристик основных событий. Событие определяется как точка во времени, в которой происходят изменения характеристик системы. Обычно изменения имеют место в тех случаях, когда кончается один процесс (или несколько процессов) и начинаются другие. Для получения требуемых результатов моделирования достаточно наблюдать систему в те моменты, когда происходят события. Для иллюстрации рассмотрим пример СМО с одним каналом (СМО с ожиданием). Оценка характера функционирования: среднее время пребывания заявки в очереди, средняя длина очереди и доля времени простоя системы. Эти характеристики могут менять свои значения либо в момент поступления дополнительного требования на обслуживание, либо при завершении обслуживания (возможны различные ситуации). Можно получить необходимую информацию, наблюдая различные условия, которые возникают при наступлении того или иного события. Для эксплуатации любой имитационной модели необходимо выбрать единицу времени. В зависимости от природы моделируемой системы такой единицей может быть минута, месяц и т.п. (для аэропорта крупного города – минута, небольшого города – час). Допустим, что надо моделировать работу системы в течение Т единиц времени. Работа начинается с данными, относящимися к нулевому моменту времени и отмечаются соответствующие события на шкале времени в хронологическом порядке. Т.о., модель функционирует, перепрыгивая от одного события к другому, непосредственно за ним следующему. Каждое событие сопровождается корректировкой протокола, отражающей возможные изменения в показателях функционирования.Резкие переходы (скачки), совершаемые моделью при переходе от одного события к другому, указывают на то, что процесс протекает в дискретном времени, откуда появилось название «дискретное моделирование». В случае дискретного моделирования между реальным временем и временем работы модели нет ничего общего (время функционирования модели обычно значительно меньше реального). Пример. Пусть мы хотим моделировать работу одноканальной СМО с ожиданием. Поступление требований в этой системе подчинено пуассоновскому распределению со средним 3 клиента/час, а время обслуживания равно 0,2 часа с вероятностью 0,5 или 0,6 часа с вероятностью 0.5. Клиенты обслуживаются согласно дисциплине «первым пришел - первым обслуживается», длина очереди, а также источник поступления клиентов не ограничены. Предположим, что в начальный момент моделирования клиентов нет. Для пуассоновского входного потока промежутки времени между требованиями имеют экспоненциальное распределение и могут быть получены из формулы:p=-(1/l) lnR = -(1/3) lnR.Время обслуживания определяется как:.В одноканальной системе возможны события только двух типов: поступление клиентов и их уход (окончание обслуживания). Действия, вызываемые этими событиями можно охарактеризовать следующим образом.^ Событие, связанное с поступлением клиента. 1. Генерация момента времени, в который поступает следующее требование на обслужив