:= t0 + h,

t0 – момент поступления заявки (содержится в строке 0).

Полученное значение  запоминается в строке k0. Переход к оператору 6.

Оператор 6 присваивает признаку , соответствующему номеру занятого канала, значение 0, символизирующее занятость канала. Переход к оператору 7.

Оператор 7 исключает из массива Е0 номеров свободных каналов номер k0 занятого канала. Переход к оператору 8.

Оператор 8 добавляет номер k0 занятого канала к массиву Е1. Переход к оператору 9.

Оператор 9 формирует случайную величину продолжительности интервала между заявками в соответствии с плотностью распределения y(t). Переход к оператору 10.

Оператор 10. Сформированная датчиком случайных чисел с плотностью распределения y(t) случайная величина x добавляется к значению t0 и, таким образом, определяется момент поступления следующей заявки: t0:= t0+x. Возврат к блоку 2 внешнего контура, контролирующему общее число заявок, прошедших через систему.

Оператор 11 выполняет действия в случае, когда в момент поступления заявок все каналы системы заняты. При этом проверяется, заполнен ли буфер. Если не заполнен (число т содержащихся в буфере заявок меньше емкости буфера М), то переход к оператору 12, в противном случае – к оператору 13.

Оператор 12 увеличивает число заявок в буфере на единицу.

Оператор 13 реализует формирование случайной продолжительности ожидания заявки в соответствии с заданной плотностью распределения N(t).

Оператор 5. Сформированная оператором 12 случайная величина H используется для расчета момента времени освобождения места в очереди. Этот момент времени вычисляется по формуле

                                                           tn+m:= t0 +H,                                             (2.7)

t0 – момент поступления заявки (содержится в строке 0).

Полученное значение tn+m запоминается в строке n+m. Переход к оператору 9.

Оператор 15 увеличивает число заявок, получивших отказ (все каналы и буфер заняты), на единицу. Переход к оператору 9.

Рассмотрим теперь операции, реализуемые в модуле 1. Блок-схема модуля 1 приведена на рис. 2.4.

Рис. 2.3 - Блок-схема модуля 1

Модуль 1 начинает работать в случае, когда самое ранее из событий, отображаемых календарем, соответствует освобождению канала с номером r0.

Оператор 1 проверяет, есть ли хотя бы одна заявка, ждущая обслуживания в буфере. Если буфер не пуст (m¹0), то переход к оператору 2, в противном случае – к оператору 5.

Оператор 2 обеспечивает формирование случайной продолжительности h занятости канала r0 при обслуживании заявки, хранившейся в буфере. Переход к оператору 3.

Оператор 3 определяет момент окончания обслуживания каналом r0 заявки, взятой из буфера. Момент освобождения канала рассчитывается по формуле

                                                      :=  + h.                                       (2.8)

Переход к оператору 4.

Оператор 4 уменьшает число заявок, хранящихся в буфере и ожидающих освобождения какого-либо канала, на единицу. Возврат к оператору 2 внешнего контура.

Оператор 5 сдвигает массив заявок, ожидающих в очереди, на 1 позицию вверх.

Оператор 6 присваивает признаку -го значение 1. В результате этой операции строка r0, соответствующая освободившемуся, но не занятому каналу (буфер пуст), при очередном просмотре календаря не будет выделена (просматриваются только те строки, для которых cj=0). Если описанную операцию присваивания :=1 не выполнить, то при просмотре календаря та же строка r0 будет выбрана вновь (этой строке соответствует минимальное время наступления события) и процедура реализации модели зациклится. Переход к оператору 6.

Оператор 6 добавляет номер r0 к массиву свободных каналов. Переход к оператору 7.

Оператор 7 исключает номер r0 из массива занятых каналов.

Рассмотрим теперь операции, реализуемые в модуле 2. Блок-схема модуля 2 приведена на рис. 2.4.

Рис. 2.4 - Блок-схема модуля 2.

Оператор 1 очищает ячейку с номером n+r0 .

Оператор 2 сдвигает массив заявок, ожидающих в очереди, на 1 позицию вверх, начиная с номера n+r0+1

Оператор 3 уменьшает количество ожидающих заявок на 1.

Завершающим этапом работы имитационной модели является статистическая обработка результатов моделирования. После завершения работы модели в памяти остаются значения общего числа заявок N0, прошедших через систему, и числа заявок, получивших отказ – s.

2.6 Оптимизация параметров системы обслуживания

Данные, полученные в результате работы ИМ, могут быть использованы для подсчета критерия эффективности L функционирования СМО:

                                                                L = Пр – Затр,                                    (2.9)

где

Пр – средняя прибыль в единицу времени, получаемая в ходе работы СМО,

Затр – средние затраты в единицу времени, связанные с функционированием СМО.

При этом

                                              Пр = C0 (Tобс) (N0 – s),                                         (2.10)

                                            Затр = C1 s + Cэ (Tобс) n.                                      (2.11)

Тогда

                                      L = C0 (Tобс)(N0 – s) – C1 s - Cэ (Tобс) n.                      (2.12)

Полученное соотношение позволяет использовать имитационную модель для оптимизации СМО.

Проведем оптимизацию СМО с помощью метода Нелдера-Мида.

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

где

   - длина ребра симплекса, выбираемая, например равной 1.

В каждой из этих точек проведем серию имитационных экспериментов и, усреднив результаты в каждой, получим оценки средних значений функции отклика . Теперь, используя стандартную процедуру Нелдера-Мида, отыскивают «худшую» точку (если решается задача максимизации, то это точка, в которой значение функции отклика минимально).

Затем реализуется один из возможных вариантов деформирования многогранника (отражение, растяжение, сокращение или редукция), после чего в новой (или новых) точке выполняется имитационное моделирование и процедура продолжается.

Рассчитаем оптимальные параметры библиотечной системы обслуживания – число каналов обслуживания  и среднее время обслуживания .

Вершины начального симплекса:

Параметры имитационной модели:

Оптимизируемой функцией является (2.12)

Критерий останова:

Результат:

Значение критерия

3 Гражданская оборона

Защита населения от оружия массового поражения и при чрезвычайных ситуациях (ЧС) достигается максимальным осуществлением всех защитных мероприятий гражданской обороны, наилучшим использованием всех способов и средств защиты. Основными способами защиты населения при ЧС являются: укрытие населения в защитных сооружениях; рассредоточение в загородной зоне рабочих и служащих предприятий, учреждений и организаций, продолжающих свою деятельность в городах, а также эвакуация из этих городов всего остального населения; использование населением средств индивидуальной защиты (СИЗ). В данной дипломной работе рассмотрен вопрос об использовании СИЗ  в электронной промышленности.

Страницы: 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17