- •Введение Понятие модели. Имитационная модель. Основные характеристики сложной системы.
- •1. Классификация моделей
- •2. Структура моделей
- •3. Схема взаимодействия компонентов системы между собой
- •4. Последовательные этапы процесса имитации
- •5. Представление исходных данных для имитации
- •6. Моделирующий алгоритм
- •7. Принципы построения моделирующих алгоритмов для сложных систем
- •8. Организация статистического моделирования систем на эвм
- •8.1. Общая характеристика метода статистического моделирования (метод Монте-Карло)
- •8.2. Алгоритм метода статистических испытаний
- •8.3. Псевдослучайные числа и процедура их генерации
- •8.4. Моделирование испытаний в схеме случайных событий
- •8.5. Формирование возможных значений св
- •8.6. Формирование реализаций случайных векторов
- •8.7 Определение необходимого числа реализаций
- •8.8. Особенности фиксации и статистической
- •8.9. Случайный процесс
- •8.10. Особенности использования критериев согласия в методах регрессионного и корреляционного анализа при обработке результатов моделирования и их интерпретации
- •8.10.1. Критерий Пирсона ( критерий 2 )
- •8.10.2. Критерий Колмогорова
- •8.10.3. Критерий Смирнова
- •8.10.4. Критерий Стьюдента
- •9. Динамическое моделирование
- •9.1 Основные теоретические положения
- •9.1.1. Основные этапы построения динамической модели
- •9.1.2. Структура динамической модели
- •9.1.3. Математическое описание динамической модели
- •9.1.4. Запаздывания
- •9.1.5. Процесс принятия решения
- •9.2. Пример анализа системы методом динамического моделирования
- •10. Регрессионный и корреляционный анализ
- •10.1. Моделирование систем массового обслуживания (смо)
- •10.2. Описание q -схем с использованием марковских случайных процессов (сп)
- •10.3. Уравнение Эрланга и формула Эрланга
- •10.4. Правила составления ду
- •10.5. Моделирование смо с помощью метода статистических испытаний
- •10.6. Формирование входного потока ( 3 -ий блок )
- •10.7. Подалгоритм выбора канала
- •10.8. Подалгоритм выбора заявки из очереди на обслуживание
- •10.9. Подалгоритм моделирования сбоев
- •10.10. Агрегаты, основные понятия
- •10.11. Процесс функционирования агрегата
- •10.12. Представление смо в виде агрегата
- •11. Регрессионный и корреляционный анализ
- •11.1. Регрессионный анализ
- •11.2. Корреляционный анализ
9.2. Пример анализа системы методом динамического моделирования
Рассматривается работа цеха, имеющего две конвейерные линии, изготовляющие детали типов А и В соответственно, и отдел сборки, выпускающий готовое изделие С. Изделие С состоит из определенного количества деталей А и В. Конвейерные линии получают заготовки со склада с определенным темпом. Ресурсы заготовок на складе ограничены, но возможны дополнительные случайные или неслучайные поступления. Сборочный отдел осуществляет контроль количества поступающих деталей А и В и регулирует производительность конвейерных линий с целью увеличения выпуска готовых изделий и улучшения ритмичности предприятия в целом. Выпуск готовых изделий будет максимальным в том случае, если в сборочном отделе не произойдет накапливания деталей типа А или В.
Рисунок 9.4 - структура цеха, выпускающего изделие С.
Методом динамического моделирования можно определить следующие характеристики рассматриваемой системы:
– число изделий, выпущенных за определенный промежуток времени Т, если ресурс заготовок на складе задан и известен характер их пополнения;
– оптимальный интервал принятия решения Δt, обеспечивающий при заданных ресурсах выпуск максимального количества готовых изделий;
– функцию принятия решения, изменяющую темпы потоков деталей в линиях и оптимизирующую производство по какому-либо критерию (например, для обеспечения постоянного времени выпуска 1-го изделия);
– графики выпуска готовых изделий и динамику изменения ресурсов с учетом пополнения и т. д.
Схема динамической модели исследуемой системы приведена на рисунке.
Рисунок 9.5 - схема динамической модели.
Пусть в начальный момент времени t = 0 в сборочном отделе имеется деталей типа А и деталей типа В. Введем некоторый коэффициент μ, характеризующий величину уровней деталей типа А и В в сборочном отделе, следующим образом:
(9.3)
где (t), (t) – уровень деталей типов А и В в момент времени t в сборочном отделе; П(A), П(В) – комплектность изделия, т. е. количество деталей типа А и В, необходимых для изготовления готового изделия.
Очевидно, что коэффициенты μ изменяются от 0 до 1. Если оба коэффициента больше или равны 1, то сборочный отдел выдает готовое изделие. Будем считать, что определение величины μ ведется только в конце интервала Δt, т. е. в момент принятия решения. Коэффициенты μ и μсравниваются, и если Δμ = (μ – μ) > Δμ, то в силу вступает обратная связь, регулирующая задержку на той или иной линии.
Допустим, что μ > μ и Δμ > Δμ, тогда темпы в линииА должны быть уменьшены за счет увеличения задержек, а в линии В – наоборот.
Темп потока в i-й линии (i = 1, 2) на выходе ј-го звена (ј = 1, 2) Х, в некоторый момент времени t + kΔt (k = 0, 1, 2,...) определяется отношением уровня y(t + kΔt) к величине задержки изделий в этом звене D, т. е.
(9.4)
Введение переменной задержки деталей в j-м звене позволяет не только представить динамику процесса выпуска изделий, но и вести управление темпами Х. Это осуществляется путем обратного воздействия сборки на величину задержки и соответственно на величину iј-гo темпа.
Минимальная и максимальная величины темпов выпуска деталей устанавливается за счет введения определенных границ D и D, в пределах которых может изменяться задержка. Величина задержки в ј-м звене в момент времени t + kΔt может быть найдена из уравнения
(9.5)
где α – коэффициент, характеризующий величину управляющего воздействия (функцию решения). Величина α может быть постоянной или определяться из некоторого уравнения, связывающего постоянные и переменные параметры системы, что позволяет уменьшить колебания ритма всей системы. На практике величину α при заданных конкретных условиях и критериях оценки можно найти путем динамического моделирования, варьируя исходные данные и параметры исследуемой системы.
В выражения (9.3)–(9.5) входят переменные уровней деталей в момент времени t +kΔt, которые должны быть вычислены в первую очередь по формуле
(9.6)
где k = 1, 2, 3, … .
Интервал времени Δt может быть выбран таким образом (при опре- деленной величине Δμдоп), что обратная связь включается в момент окончания интервала Δt не каждый раз.
Обратная связь, открывающая «клапан» Dпоступ, включается только при выполнении условия y(t + kΔt) ≤ 0, 1 от исходного y ∕ t=0 , причем поступление можно считать мгновенным или с определенным темпом