Иногда на практике результаты имитационного моделирования не удовлетворяют ни одному из этих требований. Тем не менее их выполнение гарантирует наличие корректных способов сбора наблюдений над имитационной моделью.
Рассмотрим сначала вопрос о стационарности распределений (первое требование). Результаты наблюдений над моделью зависят от продолжительности периода имитации. Начальный период неустойчивого поведения модели (системы) обычно называется переходным. Когда результаты имитационного эксперимента стабилизируются, говорят, что система работает в установившемся режиме. Продолжительность переходного периода определяется в значительной степени начальными характеристиками модели, и невозможно предсказать, когда наступит установившийся режим. В общем случае, чем длиннее продолжительность прогона модели, тем выше шанс достичь установившегося состояния.
Рассмотрим теперь второе требование, состоящее в том, что наблюдения над имитационной моделью должны иметь нормальное распределение. Это требование можно выполнить, если привлечь центральную предельную теорему (см. раздел 12.4.4), утверждающую, что распределение среднего выборки является асимптотически нормальным независимо от распределения генеральной совокупности, из которой взята выборка. Центральная предельная теорема, таким образом, есть главное средство удовлетворения требования о нормальности распределения.
Третье требование касается независимости наблюдений. Природа имитационного эксперимента не гарантирует независимости между последовательными наблюдениями над моделью. Однако использование выборочных средних для представления отдельных наблюдений позволяет смягчить проблему, связанную с отсутствием независимости. Для этого, в частности, следует увеличивать интервал времени имитации для получения выборочного среднего.
Понятия переходного и установившегося состояний имеют силу в ситуациях, именуемых незаканчивающейся имитацией, т.е. имитацией, применяемой к системам, которые функционируют бесконечно долго. При заканчивающейся имитации (например, работа банка, если он обычно работает восемь часов в день) переходное поведение является частью нормального функционирования системы и, следовательно, не может игнорироваться. Единственным выходом в такой ситуации является увеличение, насколько это возможно, числа наблюдений.
Обсудив "подводные камни" имитационного эксперимента и средства, с помощью которых их можно обойти, рассмотрим теперь три наиболее общих метода сбора информации в процессе имитационного моделирования: метод подынтервалов, метод повторения и метод циклов.
18.6.1. Метод подынтервалов
На рис. 18.10 проиллюстрирована идея метода подынтервалов. Предположим, что имитация длится на протяжении Т единиц времени (т.е. длина прогона модели равна Т) и требуется получить п наблюдений. В соответствии с методом подынтервалов необходимо сначала "обрезать" информацию, относящуюся к переходному периоду, а затем разделить остаток результатов имитации на п равных подынтервалов (групп). Среднее значение искомой величины (например, длины очереди или времени ожидания в очереди) внутри каждого подынтервала используется затем в качестве единственного наблюдения. Отбрасывание начального переходного периода означает, что статистические данные, собранные на протяжении этого периода, не используются.
Преимущество данного метода состоит в том, что влияние переходных (нестационарных) условий уменьшается, в частности, на те данные, которые собраны в конце времени имитации. Недостаток заключается в том, что последовательные группы с общей границей являются коррелированными, что приводит к невыполнению предположения о независимости. Влияние корреляции может быть уменьшено путем увеличения интервала времени для каждой группы.
Переходный
период Группа 1 Группа 2
Группа п
Время имитации
Рис. 18.10. Иллюстрация к методу интервалов
Пример 18.6.1
На рис. 18.11 показано изменение длины очереди в системе обслуживания с одним сервисом (модель простой очереди) как функции времени. Период имитации составляет Т= 35 часов, а длина переходного периода оценивается в 5 часов. Необходимо получить 5 наблюдений, т.е. п = 5. Соответствующая длина интервала времени для каждой группы равна (35 - 5)/5 = 6 часов.
Длина
очереди Q
Переходный | | | |
период Группа 1 Группа 2 | , Группа 3 , | Группа 4 | Группа 5 , |
~*---* . ,J L 2=оГ | т*-»1 Lj3=n | **-** А4=6Г> | **-Н |
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 | .....1 | . 1 Г. . 1 | Л5=15Ц 1 1 1 1 1 1 |
J I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I
5 10 15 20 25 30 35
Время имитации Рис. 18.11. Изменение длины очереди
Пусть <2, представляет среднюю длину очереди в группе /. Так как длина очереди является переменной, зависящей от времени, то
Q,=-, / = 1,2, ...,5,
где А-, - площадь под кривой длины очереди, t - длина интервала времени для группы. В рассматриваемом примере t = 6 часов.
Анализ данных, приведенных на рис. 18.11, приводит к следующей таблице.
Наблюдение/ | | | | | |
| | | | | |
| 2,33 | 1,67 | 1,83 | 1,00 | |
Выборочное среднее | = 1,87 | | Выборочная дисперсия = 0,35 | |
Выборочные среднее и дисперсию можно использовать, если это необходимо, для вычисления доверительного интервала.
Вычисление выборочной дисперсии в этом примере основано на использовании следующей хорошо известной формулы:
5>.-*)2
s2=-.
Эта формула является лишь приближением точного значения дисперсии, так как не учитывает эффекта автокорреляции между последовательными группами. Точную формулу можно найти в работе [2].
18.6.2. Метод повторения
В данном методе каждое наблюдение представляется независимым прогоном (имитацией) модели, в котором переходный период не учитывается, как показано на рис. 18.12. Вычисление средних величин выборки для каждой группы проводится точно так, как и в методе подынтервалов. Единственное отличие в том, что в данном случае стандартная формула для дисперсии применима, так как группы не коррелированы между собой.
Преимуществом этого метода является то, что каждый имитационный прогон модели определяется своей последовательностью случайных чисел из интервала [О, 1], что действительно обеспечивает статистическую независимость получаемых наблюдений. Недостаток состоит в том, что все наблюдения могут оказаться под сильным влиянием начальных переходных условий. Этот недостаток можно смягчить, увеличив длину прогона модели.
Рис. 18.12. Иллюстрация к методу повторения