- •Цели и задачи дисциплины «Имитационное моделирование»
- •Тематика курса
- •Тема 1. Сущность, развитие и использование имитационного моделирования
- •Тема 2. Основные этапы построения имитационных моделей
- •Тема 3. Имитационная модель управления запасами
- •Тема 4. Метод Монте-Карло в имитационном моделировании
- •Тема 5. Генерирование случайных чисел и случайных величин во время машинной имитации
- •Тема 6. Планирование имитационных экспериментов
- •Тема 7. Общая и статистическая проверка результатов имитационных экспериментов
- •Тема 8. Планирование имитационных экспериментов в процессе исследования и оптимизации систем
- •Тема 1. Сущность, развитие и использование имитационного моделирования
- •1.1. Моделирование. Классификация видов моделирования.
- •Классификация видов моделирования:
- •1.2. Понятие имитационного моделирования.
- •1.3. Этапы развития имитационного моделирования.
- •1.4. Сферы применения имитационного моделирования.
- •1.5. Вопросы для самоконтроля по теме № 1.
- •Тема 2. Основные этапы построения имитационных моделей
- •2.1. Практическая реализация имитационного моделирования.
- •2.2. Основные этапы построения логической схемы имитационных моделей.
- •2.3. Вопросы для самоконтроля по теме №2.
- •Тема 3. Имитационная модель управления запасами
- •3.1. Задача оптимального управления запасами.
- •3.2. Детерминированная модель управления запасом.
- •3.3. Концептуальная модель оптимизации запасов.
- •3.4. Вопросы для самоконтроля по теме № 3.
- •Тема 4. Метод Монте-Карло в имитационном моделировании
- •4.1. Статистическое моделирование.
- •4.2. Метод Монте-Карло.
- •4.1. Вопросы для самоконтроля по теме № 4
- •Тема 5. Генерирование случайных чисел и случайных величин во время машинной имитации
- •5.1. Генератор случайных чисел.
- •5.2. Генерирование равномерной случайной последовательности чисел.
- •1. Табличный генератор.
- •2. Физический генератор.
- •3. Программный генератор
- •5.3. Имитация случайных событий.
- •5.4. Имитация дискретных случайных величин.
- •5.5. Имитация непрерывных случайных величин.
- •5.6. Вопросы для самоконтроля по теме № 5
- •Тема 6. Планирование имитационных экспериментов
- •6.1. Задачи планирования экспериментов.
- •6.2. Основные понятия и определение.
- •6.3. Управляемые и неуправляемые факторы.
- •6.3. Факторные планы.
- •Матрица планирования для
- •6.4. Вопросы для самоконтроля по теме № 6
- •Тема 7. Общая и статистическая проверка результатов имитационных экспериментов
- •7.1. Проверка однородности дисперсий.
- •7.2. Проверка адекватности модели.
- •7.3. Вопросы для самоконтроля по теме № 7.
- •Тема 8. Планирование имитационных экспериментов в процессе исследования и оптимизации систем
- •8.1. Планирование экспериментов при исследовании систем.
- •8.2. Планирование экспериментов для оптимизации систем.
- •8.3. Реализация имитационной модели средствами пакета имитационного моделирования дискретных систем gpss.
- •8.3.1. Структура системы моделирования gpss World
- •8.3.2. Объекты языка имитационного моделирования
- •Объекты gpss
- •8.3.3. Построение моделей с одноканальными устройствами, функционирующими в режиме занятия и освобождения устройства
- •8.3.4. Организация поступления транзактов в модель и удаления транзактов из нее
- •8.3.5. Разработка и эксплуатация моделей в gpss World. Создание объекта «Модель»
- •8.3.6. Стандартный отчет модели – структура файла report
- •8.3.7. Построение моделей.
- •8.4. Вопросы для самоконтроля по теме № 8
- •Используемая литература
- •Содержание.
6.2. Основные понятия и определение.
При выполнении экспериментальных исследований природа исследуемого процесса большей частью известна мало (в противном случае можно составить и исследовать математическую модель объекта, избегнув опытов, которые связаны со значительными затратами). Поэтому схема эксперимента напоминает схему кибернетического «черного ящика», свойства которого устанавливаются анализом зависимости измеренных эндогенных величин при определенных комбинациях уровней контролируемых входных величин.
Пусть y - эндогенная величина, а Х1, Х2, ..., Хn - контролируемые во время исследование (опытов) факторы. Тогда процесс, который изучается, может быть описан с помощью математической модели
у = f (Х1, Х2, ..., Хn)
Функцию f (Х1, Х2, ..., Хn) называют функцией (реакцией, поверхностью) отклика, а величины Х1, Х2, ..., Хn – факторами.
Рассмотрим определения основных категорий планирования экспериментов.
Отклик – эндогенная (обусловленная внутренними причинами) случайная величина, которая по предположению, зависит от факторов.
Факторы Х1, Х2, ..., Хn – переменные величины, которые по предположению, влияют на результаты экспериментов.
Функция отклика – математическая модель у = f (Х1, Х2, ..., Хn), которая представляет собой зависимость математического ожидания от факторов.
Поверхность отклика для функции одной переменной вырождается в линию. Для двух факторов поверхность отклика представляет собой поверхность в трехмерном пространстве. Эту поверхность для наглядности можно подать графически на плоском рисунке с помощью линий одинакового уровня (так, как это делается при создании топографических карт). Полученное изображение целевой функции является удобным для иллюстрирования особенностей поверхности отклика во время поиска экстремума.
6.3. Управляемые и неуправляемые факторы.
Факторы могут быть либо количественными, либо качественными. Количественные факторы, как правило, предполагают численные значения, тогда как качественные факторы обычно являют собой структурные допущения, которые не измеряются количественно. Следует заметить, что у некоторых факторов различие может быть не ясно.
При проведении имитационных экспериментов различают управляемые и неуправляемые факторы в зависимости от того, может ли руководство соответствующих реальных систем управлять ими.
В некоторых примерах из таблицы расположенной ниже возможность управления факторами зависит от ситуации; скажем, изменить резервный запас на фабрике можно при одних обстоятельствах и нельзя при других.
Таблица.
Примеры факторов и откликов
Возможные факторы |
Количест-венные? |
Качест-венные? |
Управ-ляемые? |
Неуправл-яемые? |
Возможные |
Контрольно-кассовый пункт в супермаркете |
|||||
Среднее время между прибытиями покупателей |
Да |
|
|
Да |
Задержка в очереди |
Среднее время обслуживания |
Да |
|
Да |
|
Время пребывания в системе |
Число проходов |
Да |
|
Да? |
Да? |
Длина очереди |
Наличие скоростных проходов |
|
Да |
Да |
|
Коэффициент занятости контролеров |
Стратегия добавления или снятия контролеров |
Да? |
Да? |
Да |
|
Частота добавления или снятия |
Производственная система |
|||||
Число станков |
Да |
|
Да? |
Да? |
Производительность, в деталях
|
Дисциплина обслуживания очереди |
|
Да |
Да? |
Да? |
Время пребывания детали в системе |
Резервный запас |
Да |
|
Да? |
Да? |
Коэффициент использования станков |
Скорости конвейеров |
Да |
|
Да? |
Да? |
Доходность |
Объединение станков по участкам |
|
Да |
Да? |
Да? |
|
(продолжение таблицы)
Возможные факторы |
Количест-венные? |
Качест-венные? |
Управ-ляемые? |
Неуправл-яемые? |
Возможные |
Сеть |
|||||
Интенсивность поступления сообщений |
Да |
|
|
Да |
Задержка сообщения |
Длительность сообщений |
Да |
|
|
Да |
Пропускная способность |
Число узлов |
Да |
|
Да? |
Да? |
Коэффициент использования узлов и линий |
Число линий |
Да |
|
Да? |
Да? |
Длина очередей |
Используемый протокол |
|
Да |
Да |
|
|
Режим эксплуатации |
|
Да |
Да |
|
|
Система управления запасами |
|||||
Среднее время между возникновением спроса |
Да |
|
|
Да |
Затраты на хранение |
Количество запрашиваемого товара |
Да |
|
|
Да |
Издержки, связанные с нехваткой товара |
Время получения заказа |
Да |
|
Да? |
Да? |
Коэффициент использования станков |
Количество товара, на которое подается заказ |
Да |
|
Да |
|
|
Поставки, выполнение которых было отложено, в сравнении с утраченными объемами продаж |
Да |
Да |
Да? |
Да? |
Поставки, выполнение которых было отложено, или утраченные поставки |
Рассматриваются преимущественно управляемые факторы имитационных экспериментов, так как они имеют наибольшее влияние на решения, которые могут быть приняты относительно реализации реальных систем. Неуправляемые факторы тоже представляют интерес при проведении имитационных экспериментов, поскольку возникает необходимость оценить, каким образом резкое увеличение интенсивности прибытия клиентов повлияет на загруженность системы. В имитационном моделировании, доступно управление всеми факторами, независимо от их практической управляемости.
Выдвигается ряд требований не только к отдельным факторам, а и ко всей совокупности факторов в целом. Прежде всего, факторы должны быть независимыми, то есть установление некоторого уровня одного фактора не должно зависеть от взятых значений других факторов. При проведении физических экспериментов важным является свойство совместимости факторов (любая комбинация их уровней может быть осуществима и безопасна).
При планировании важно выбрать наиболее удачную функцию отклика у = f (Х1, Х2, ..., Хn) - математическую модель процесса, который изучается. Поскольку истинное описание этой функции установить невозможно, то ее подают приближенно с помощью аппроксимирующего полинома – отрезка ряда Тейлора, в который раскладывается неизвестная зависимость в окрестности точки с нулевыми координатами:
где
Аппроксимирующий полином берут первого, второго и, реже, третьего порядка, причем степень может изменяться в зависимости от этапа эксперимента или специфики решаемой задачи. Коэффициенты полинома β из-за отсутствия истинного описания функции не могут быть найдены теоретически. Их определяют экспериментально, проводя опыты при некоторых фиксированных значениях факторов, причем, учитывая случайность исследуемого процесса или в результате ошибок при измерениях исходной величины, попытки дублируются.
После обработки результатов эксперимента исходный аппроксимирующий полином фактически заменяется уравнением регрессии:
где коэффициенты b0, bi, bij, bii - статистические оценки неизвестных теоретических коэффициентов β0, β i, βij, βii.
После выбора вида функции отклика необходимо установить границы областей определения факторов, их основные уровни и интервалы изменения. При выполнении экспериментов для каждого фактора выбирают два уровня, на которых можно изменять его значения. Выбор основных уровней (начальной точки факторного пространства) и интервалов изменения основывается на предшествующих знаниях относительно исследуемого процесса. Делают этот выбор, исходя из необходимости уменьшения числа попыток для решения поставленной задачи.
Итак, для любого i-го фактора определяются:
– основной уровень
i-го
фактора;
– верхний уровень
і-го
фактора;
–
нижний уровень
і-го
фактора;
–
интервал изменения.
Во время проведения экспериментов используются кодированные значения уровней факторов. При этом основной уровень берут равным нулю, верхний +1, нижний -1. Кодирование выполняют по формуле:
,
(i
= 1, 2, …, n)
Количество всех точек факторного пространства при двухуровневой системе изменения факторов, в которых необходимо экспериментально находить значения функции отклика, равняется 2n, где n – число факторов.
Если при проведении опытов возможны ошибки измерения эндогенной величины y или она является случайной величиной и оценивается с помощью математического ожидания (что является характерным, например, для задач управления запасами), то за расчетное значение функции отклика в некоторой точке факторного пространства берут среднее значение этого показателя, полученное обработкой результатов установленного заранее числа дублирований эксперимента. Число повторений опыта определяется статистическими методами в зависимости от взятой надежности и точности ожидаемых результатов.
Для определения коэффициентов b уравнения регрессии экспериментально находят значения эндогенной величины y в N точках факторного пространства (при двухуровневой системе смены факторов N ≤ 2), причем попытки в каждой точке дублируются k раз.
В общем случае число повторений опытов (параллельных попыток) в точках факторного пространства может быть разным, тем не менее, в имитационных моделях экономических систем это число целесообразно считать единым для всех точек факторного пространства, в которых будет производиться эксперимент.
Задача определения коэффициентов регрессии типична для регрессионного анализа. Поэтому для корректного использования разработанного в этой теории инструментария необходимо, чтобы выполнялись рассмотренные далее условия:
1. Результаты наблюдений y1, y2, ..., yn величины y в N точках факторного пространства представляют собой реализацию нормально распределенной случайной величины.
2. Дисперсии реализаций σ2{yi} (i = 1, 2, …, N) равны между собой, то есть дисперсия величины y не зависит от ее абсолютного значения.
3. Факторы Х1, Х2, ..., Хn - независимые величины, измеряемые с настолько маленькой погрешностью, которой можно пренебречь в сравнении с погрешностью определения величины у.
В случае имитационного моделирования третье условие выполняется всегда, а первое и второе должны проверяться с помощью специальных тестов.
