Федеральное Государственное Образовательное Учреждение Высшего Профессионального Образования Ставропольский Государственный Аграрный Университет Лекция. - презентация

1 Федеральное Государственное Образовательное Учреждение Высшего Профессионального Образования Ставропольский Государственный Аграрный Университет Лекция 13 на тему «Проблемно-ориентированные имитационные модели» по дисциплине «Имитационное моделирование экономических процессов»

2 Вопросы: 1. Модель «посещение пунктов местности коммивояжером» 2. Модель «стоянка маршрутного такси» 3. Модель «эффективность компьютеров в автоматизированной бухгалтерии» 4. Модель «минимизация производственных затрат» 5. Модель «Динамическое распределение ресурсов в сети под управлением UNIX» Контрольные вопросы.

3 1. Модель «Посещение пунктов местности коммивояжером» Предположим, что поставлена следующая задача: необходимо посетить все города, показанные на рис. 1, по кратчайшему маршруту, причем каждый город можно посетить только один раз. Коммивояжер находится в Санкт-Петербурге и имеет в своем распоряжении вертолет. Рассмотрим три варианта решения такой задачи.

4 Рис. 1. Развертка выбранного региона на плоскости (учебный пример)

5 Вариант 1. Коммивояжер, пользуясь только картой, строят маршрут по алгоритму «ближайшего непосещенного города». Чтобы не запутаться в большом числе частных маршрутов, он применяет сканирование на определенной широте. В данном случае он выполняет такие действия: из Санкт-Петербурга проводит маршрут в западном направлении (ближайший город - Рига); поскольку на данной широте не нашлось ближайшего города западнее, то далее маршрут строится южнее в восточном направлении (ближайший город ~ Москва); так как на данной широте не нашлось ближайшего города восточнее, далее маршрут строится южнее в западном направлении (ближайшая оптимальная последовательность городов - Плавск, Вильнюс, Балтийск);

6 поскольку на данной широте не нашлось ближайшего города западнее, то далее маршрут строится южнее в восточном направлении (ближайшая оптимальная последовательность – Минск, Яровщина, Клинцы); ввиду того что на данной широте не нашлось ближайшего города восточнее, далее маршрут строится южнее в западном направлении (ближайшая оптимальная последовательность - Новозыбков, Хойники, Пинск); остался последний не посещенный город южнее - Киев; туда проводится последний отрезок маршрута. Построенный маршрут (рис. 2) выглядит практически без изъянов: он действительно короткий, на нем нет петель (т.е. выполнено необходимое условие построения оптимального маршрута). Далее измеряем протяженность.

7 Рис. 2. Пример использования алгоритма «ближайшего непосещенного города»

8 Решение задачи: длина маршрута равна 4526,2 км. Вариант 2. Введем более жесткие условия: пункт старта - Санкт-Петербург; пункт финиша - Москва. В данном случае алгоритм «ближайшего непосещенного города» применить не удается: маршрут получится длиннее по сравнению с вариантом 1. Методы линейного и динамического программирования, метод «ветвей и границ» дают не лучшие решения при таком на первый взгляд несущественном усложнении задачи. Поэтому используем алгоритм «двух вертолетов». В результате создается маршрутное расписание Сам маршрут имеет два характерных участка, которые не позволяют его считать оптимальным, причем имеется петля (не выполнено необходимое условие оптимальности маршрута). Далее измеряем протяженность.

9 Решение задачи: длина маршрута равна 3338,2 км. На первый взгляд, это решение парадоксально: маршрут получился короче на 1188,0 км по сравнению с первым вариантом. Но объяснение такого эффекта простое: при большом числе населенных пунктов алгоритм «двух вертолетов» всегда дает решение, лучшее по сравнению с алгоритмом «ближайшего непосещенного города» по следующей причине: оптимизация осуществляется двунаправленным просмотром. Алгоритм «ближайшего непосещенного города» начинает запутывать маршрут и при большом числе пунктов «уводить» его от оптимального, так как в конец маршрута заглянуть невозможно. Если иметь средства автоматического отображения маршрута, то неоптимальные участки сразу видны. Поэтому возможна корректировка маршрута вручную.

10 Вариант 3. Используются условия варианта 2. Применяется алгоритм «двух вертолетов». Разрешена ручная корректировка неоптимальных участков. В результате получаем оптимальное маршрутное расписание. Далее измеряем протяженность. Решение задачи: длина оптимального маршрута равна 2807,0 км. Это на 531,2 км короче по сравнению с вариантом 2 (но без варианта 2 этот маршрут нельзя было бы построить). Далее необходимо разработать имитационную модель, которая должна обладать следующими свойствами: маршрутное расписание представляется в виде массива space, количество элементов которого - число населенных пунктов, а номер элемента массива - порядок посещения этого пункта; модель должна автоматически измерять длину маршрута (км) и время в пути (час).

11 Такая модель позволяет сравнивать эффективность различных алгоритмов. Структура модели показана на рис. 3. Узел «Генератор вызовов» поставляет транзакты в модель. Очередь 2 «Пункты вызовов» привязывает каждый транзакт к номеру пункта: если номер транзакта i, то он получает координаты из элемента i массива space, i=1,2,...,М. Функцию накопления эта очередь не выполняет.

12 Рис. 3. Схема моделирования маршрутного расписания коммивояжера

13 Маршрутное расписание в виде упорядоченных транзактов образуется в узле 3 «Расписание». Транзакты поступают на вход узла «Вертолет», который имитирует полет с очередному населенному пункту. По завершении полета к каждому пункту соответствующий транзакт удаляется из узла «Вертолет», а далее и из модели (терминатор 5 «Пункт посещен»). Узел «Вертолет» автоматически измеряет пройденный путь. Время полета фиксируется с помощью клапана (узел 6 «Время полета»). Этот клапан закрывается в начале имитации полета (операция hold из узла 2 queue при прохождении первого транзакта), а открывается в конце полета (операция rels из узла 5 term при удалении транзакта с номером М).

14 2. Модель «Стоянка маршрутного такси» Задача взаимного исключения возникает в следующих случаях: при создании сложных программных мультипроцессных средств; для моделирования логики экономического (логистического) процесса. Рассмотрим практическое решение такой задачи. Имеется остановка маршрутного такси. С определенными интервалами времени на остановку приходят пассажиры и подъезжают машины такси. Такси уезжает, когда в машину сядут 10 пассажиров. Если пришедший пассажир не обнаруживает стоящей машины такси, то он встает в очередь на посадку. Когда такси подъезжает на пустую остановку (нет ни одного человека), машина ждет пассажиров.

15 Принцип работы узла delet: вошедший в него транзакт семейства «такси» становится транзактом семейства «пассажиры», входящие следом. Как только в узел войдет заданное количество транзактов-пассажиров, они поглотятся, а транзакт «тележка» переходит в следующий узел; транзакты, собранные в нее, перестают существовать как заявки на обслуживание в системе. Параметры модели представлены в табл. 1. Таблица 1 Параметры имитационной модели п/п Параметр Среднее значение (минуты) 1Время моделирования: 10 ч 600,0 2Интервал прихода пассажиров 1,0 3Интервал приезда такси 10,0 Параметры имитационной модели

16 Опишем работу модели. Генератор 7 (рис. 4) создает транзакты, имитирующие пассажиров. Интервал прихода пассажиров в соответствии с теоремой о суперпозиции потоков событий имеет экспоненциальное распределение. Интервал генерации транзактов, имитирующих такси, имеет нормальное распределение (обоснование: маршрут такси состоит из множества отрезков, поэтому имеют место условия центральной предельной теоремы, даже если на маршруте несколько машин). Использование в описании генератора 8 нормального закона распределения интервала генерации означает, что время между приездами такси на остановку чаще оказывается ближе к своему среднему значению и реже - дальше от него (чем больше отклонение интервала от среднего, тем реже это бывает).

17 Рис. 4. Схема модели взаимного исключения и координации, процессов типа «производитель- потребитель»

18 Узлы 1 и 5 имитируют соответственно очереди пассажиров и такси. Ключ 6 в начале работы модели находится в открытом состоянии (по умолчанию), а ключ 2 закрывается при приходе в очередь 1 первого транзакта- пассажира. Это делается для того, чтобы в узел delet первым вошел транзакт-такси, а не пассажир (иначе пассажир станет «тележкой»). Как только в delet войдет первый транзакт-такси, ключ 6 закрывается, а ключ 2 одновременно с ним открывается. Теперь выходящие из генератора 6 транзакты-такси стоят в очереди 5, а транзакты-пассажиры заходят в узел delet (идет заполнение такси). Когда в узле delet накопятся 10 транзактов-пассажиров, транзакт-такси перейдет в терминатор (заполненное такси уедет). После этого ключ 6 открывается, чтобы в delet мог зайти следующий транзакт-такси, а ключ 2 в это же время закрывается, чтобы транзакты-пассажиры стояли в очереди 1 до прихода в delet транзакта-такси.

19 Таким образом, ключи 2 и 6 в модели всегда находятся в противоположных состояниях: если один открыт, то другой закрыт. Тем самым чередуется доступ к узлу delet разных типов транзактов - такси и пассажиров. Рассмотренный принцип «тележки» позволил эффективно решить задачу взаимного исключения и координации процессов «производитель-потребитель». Нетрудно заметить, что в модели нет ни одного узла обслуживания, но в очередях могут возникать существенные задержки.

20 3. Модель «Эффективность компьютеров в автоматизированной бухгалтерии» Рассмотрим в несколько упрощенном изложении процесс обработки бухгалтерской информации в фирме с помощью экономической информационной системы, работающей в бухгалтерии. Такие системы часто называются автоматизированным рабочим местом (АРМ). Далее попытаемся создать структурную схему АРМ бухгалтерии и имитационную модель, соответствующую этой схеме. В бухгалтерии выделим три подразделения (рис. 5): группа бухгалтеров в рабочей комнате бухгалтерии; архив, где документы, обработанные на компьютере, подшиваются в архивные папки; компьютерное подразделение, работающее в компьютерном зале бухгалтерии.

21 Рис. 5. Иллюстрация потока бухгалтерских документов

22 В рабочей комнате N бухгалтеров (N 1) готовят независимо друг от друга различные документы, которые направляются в компьютерный зал для дальнейшей подготовки или обработки. В этом зале имеется один столик, на котором каждый документ находится до того, как оператор возьмет его на обработку. На таком столике может возникнуть очередь документов, ожидающих обработки. Обработка выполняется оператором, находящимся за компьютером, в течение какого-то интервала времени. Если в компьютерном зале несколько компьютеров, объединенных в сеть, то за каждым компьютером должен быть свой оператор. Общее число компьютеров М(М 1). Обработанные (или распечатанные) на компьютере документы поступают в архив. Состав сотрудников архива нас не интересует.

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

24 сколько компьютеров необходимо для того, чтобы справиться с потоком документов в бухгалтерии ? сколько времени будет находиться документ, поступивший из рабочей комнаты, на столике в компьютерном зале, прежде чем его возьмет оператор на обработку? Рассмотрим три основные причины, по которым компьютеры не смогут справиться с обработкой документов. Введем обозначения: ta - средний интервал времени между двумя последовательными поступлениями документов из двери рабочей комнаты на столик в компьютерном зале; ts - средний интервал времени обработки документа оператором на компьютере.

25 Первая причина, по которой компьютеры не справятся с обработкой потока документов, это высокая интенсивность потока: ta ts / M. Будем предполагать, что владелец фирмы достаточно хорошо владеет арифметическими расчетами, чтобы выбрать столько компьютеров, сколько необходимо для обработки потока. Вторая причина заключается в том, что даже при выполнении соотношения ta > ts / М интервал поступления и интервал обработки - это случайные величины, которые имеют значительные отклонения от средних значений.

26 Третья причина - это ремонт компьютера в связи с неисправностями, возникающими в случайные моменты времени. Необходимый ремонт увеличивает разброс времени обработки документов и является дополнительной причиной снижения пропускной способности компьютера. Если допустить, что в зале один компьютер и количество бухгалтеров велико ( т.е. N >> 1), то можно считать поток документов простейшим, интервал поступления - случайной величиной, распределенной по экспоненциальному закону, а интервал обслуживания - случайной величиной, распределенной по любому закону со среднеквадратичным отклонением, обозначенным через.

27 При таких допущениях можно применить точную формулу Поллачека-Хинчина и оценить среднее время задержки документа на столике tq. Следует отметить, что при произвольном числе компьютеров N > 1 также можно найти аналитическое решение для опреде­ления времени tq. В нашем случае применение методов исследования операций становится невозможными по двум причинам: во- первых, число бухгалтеров не может быть большим и, во- вторых, присутствует фактор неисправности. Если попытаться создать математическую модель с помощью, например, аппарата вложенных цепей Маркова (метод Кендалла) или аппарата полумарковских процессов, то придется ввести большое число допущений, которые сделают погрешность метода при определении tq крайне большой (буквально «плюс-минус в несколько раз»).

28 В аспекте системного анализа можно представить модель «АРМ бухгалтерии» в виде «черного ящика» и выделить параметры четырех потоков (два входных и два выходных): входного потока документов; входного потока неисправностей; обработанного потока документов; потока устраненных неисправностей. Структурный анализ данной системы довольно просто закончить уже на втором слое. Рассмотрим возможную схему модели (рис. 6). Каждый бухгалтерский документ, исходящий из рабочей комнаты, будем считать неприоритетным транзактом. Тогда все бухгалтеры - это один генератор транзактов бесконечной емкости.

29 Рис Структурная схема модели «АРМ бухгалтерии»

30 Представим неисправности в виде особых транзактов, которые могут занять компьютер, прервав обработку документа - неприоритетного транзакта. Для этого необходимо ввести в рассмотрение второй генератор транзактов бесконечной емкости. Столик-очередь представим в виде узла типа очередь. Будем считать, что документы забирается со столика в хронологическом порядке по принципу «первым пришел - первым обслужен», т.е. оператор берет нижний документ, если вновь поступающие документы помещаются сверху на пачку необработанных документов. Документы обрабатываются на одном (или нескольких) компьютере. Процесс обработки представим в виде узла - обслуживающего прибора (узел 4 модели serv).

31 Обслуживающий прибор имеет М независимых каналов обслуживания (M 1). Логично допустить, что каждый канал - это компьютер со своим оператором, обрабатывающим документы. Канал занят, если соответствующий оператор взял документ и обрабатывает его. Архив представим в виде особого узла - терминатора, в который транзакты входят и в нем остаются (узел 5 модели term).

32 Рассмотрим путь транзакта-неисправности. Такой нематериальный транзакт выходит из второго генератора. Этот транзакт должен обладать абсолютным приоритетом по сравнению с транзактом-документом. Если на столике находятся несколько документов, то неисправность их «отодвигает», проникает в компьютер, и если на АРМ проходила обработка документа, то транзакт-неисправность прерывает эту обработку. Прерванный транзакт-документ после этого попадает в специальный накопитель (стек), где находится до тех пор, пока не будет отремонтирован компьютер. По завершении ремонта транзакт-неисправность будет направлен в особый узел модели - терминатор (узел 6 модели term), в котором он уничтожается, а прерванный транзакт-документ возвращается из стека, и обработка его будет продолжена.

33 Длительность ремонта компьютера всегда можно представить в виде суммы длительностей многочисленных элементарных операций - случайных величин, распределенных по произвольному закону. Поэтому в соответствии с центральной предельной теоремой будем считать, что эта длительность распределена по, нормальному закону (значение соответствующего параметра - norm). Если допустить, что процесс обработки оператором одного документа также можно представить в виде последовательности многочисленных элементарных операций и длительность каждой элементарной операции распределена по произвольному закону, то можно считать, что интервал обработки одного документа - это величина, распределенная по нормальному закону (значение соответствующего параметра - также norm).

34 Рассмотрим модель, которая относится к замкнутым имитационным моделям. Предположим, что мы открыли собственную швейную фабрику «Мануфактура Ко». На этой фабрике есть швейный цех, где имеются 50 рабочих мест со швейными машинами. Машины выходят из строя и требуют ремонта. Поэтому на фабрике есть ремонтный цех, в котором работают специалисты - наладчики машин. Среднее время безотказной работы одной машины 157,0 ч. Естественно допустить, что это время - случайная величина, распределенная по экспоненциальному закону в соответствии с предельной теоремой о суперпозиции потоков. Допущение основано на том, что машина - сложное устройство, состоящее из сотен деталей и узлов, которые могут выйти из строя. 4. Модель «Минимизация производственных затрат»

35 Среднее время ремонта машины равно 7,0 ч; среднеквадратичное отклонение времени ремонта - 3,0 ч. В соответствии с центральной предельной теоремой считаем, что время ремонта, состоящего из многих последовательных элементарных операций, распределено по нормальному закону. Сформулируем следующую задачу: сколько нужно арендовать резервных машин j (дополнительно к 50 собственным) и сколько необходимо нанять наладчиков i, чтобы сделать минимальными затраты на производство, связанные с ремонтом и наладкой машин? Все исходные данные для решения поставленной задачи сведены в табл. 3. В качестве денежной единицы выберем американский доллар (только для определенности).

36 Таблица 3 Параметры швейной фабрики п/п Параметр швейного производства Значение параметра, 1Название швейной фабрики«Мануфактура Ко» 2Заработная плата одного наладчика (в час)3,75 долл. 3Оплата за аренду одной швейной машины (в день)30,0 долл. 4Убыток из-за простоя одного рабочего места по причине неисправности швейной машины (в час) 20,0 долл. 5Среднее время безотказной работы одной швейной машины 157,0 ч 6Среднее время ремонта (наладки) одной машины 7,0 ч 7Среднеквадратичное отклонение времени ремонта (наладки) одной машины 3,0 ч 8Продолжительность одного рабочего дня 8,0 ч 9Число рабочих часов в одну неделю 40,0 ч 10Число рабочих недель в одном году 52 недели

37 Сначала попытаемся решить эту задачу без применения имитационной модели. Допустим, что необходимо j дополнительных машин, которые будем арендовать на каком-либо складе. Нетрудно представить, что каждая машина находится в двух состояниях: находится в швейном цехе и работает; находится в ремонтном цехе, неисправна и ремонтируется наладчиками. Средняя длительность такого цикла равна = 164. Поэтому введем в рассмотрение полезную загрузку одной машины р = 157/164 = 0,957. Составим простое уравнение, показывающее, сколько нужно иметь дополнительных машин, чтобы все 50 мест не простаивали: 50 = (50+j).

38 Решение этого уравнения элементарно: = 50(1-р) = 2,15 машин. Однако решение должно быть целочисленным. Поэтому округляем в большую сторону и получаем целое число у = 3. Время ремонта составляет 7 ч. Учитывая, что постоянно будут неисправны 3 машины, а рабочий день составляет 8 ч, наймем трех наладчиков, которые будут хорошо загружены в течение дня (два наладчика могут не справиться с потоком машин, поступающих в ремонт). Получили искомую величину: i = 3. Если поступим в соответствии с полученным решением и будем проводить натурное моделирование, собирая в течение длительного периода времени (не менее года) хронометрическую и финансовую информацию, то увидим, что решение, полученное выше, неверное. Основная ошибка заключалась в гипотезе о том, что машины находятся в двух состояниях. На самом деле таких состояний четыре (рис. 8).

39 Рис. 8. Иллюстрация к задаче оптимизации производственных затрат

40 Выбранные выше состояния (исправная работа и ремонт) осуществляются в двух цехах: швейном и ремонтном. В начальный момент, при организации производства, все машины были исправны. Введем в рассмотрение следующие переменные: Nowon - количество рабочих мест в швейном цехе и соответственно количество собственных швейных машин; Arend - число арендуемых дополнительных машин для замены вышедших из строя ; Men - количество наладчиков, производящих ремонт (наладку) швейных машин. Отметим все значимые для задержек состояния в табл. 4.

41 Таблица 4 Состояния машин Номер фразы Практический смысл фазы Тип узла в модели Число каналов 1Помещение, где находятся исправные швей­ные машины (их начальное число Arend=j ) queue- 2Швейный цех, в котором имеется Nowon рабочих мест для Nowon=50 швейных машин servNowon 3Помещение, где находятся неисправные швейные машины, ожидающие ремонта queue- 4Ремонтный цех. в котором работают специалисты- наладчики, число которых Меп=i senMen

42 В начальный момент решено арендовать Arend машин. Эти машины должны некоторое время находиться в каком-то помещении и «ждать», пока не выйдет из строя машина в цехе. Далее в процессе выхода из строя и восстановления машин в этом помещении будут находиться машины, число которых - случайная величина. Пребывание машин в этом помещении является состоянием 1. Пребывание машин в швейном цехе в исправности - состояние 2. Если какая-либо машина выходит из строя, то она поступит в ремонт только в том случае, если хотя бы один наладчик свободен (не занят ремонтом другой машины). Вероятность такого события, естественно, не равна единице. Поэтому неисправные машины в ожидании ремонта направляются в другое помещение. Пребывание в этом помещении является состоянием 3.

43 И, наконец, если подходит очередь ремонта неисправной машины, она поступает в ремонтный цех. Ремонт - состояние 4 каждой машины. Далее начнем строить имитационную модель. Наиболее рациональное решение заключается в том, чтобы выбрать машины в качестве транзактов. В табл. 4 приведены типы узлов, имитирующих соответствующие состояния. Количество каналов узла serv в состоянии 2 равно Nowon, а число каналов узла serv в состоянии 4 - это значение Men. Единственная проблема, с которой мы немедленно столкнемся, заключается в том, что в данной замкнутой схеме нет генераторов транзактов, В связи с этим непонятно, откуда тракзакты попадут в модель.

44 Узлы с номерами 1 (queue), 2 (serv), 3 (queue) и 4 (serv) имитируют пребывание машин в состояниях с соответствующими номерами. Эта замкнутая схема не требует дальнейших комментариев. Далее рассмотрим схему зарядки, которая состоит из трех узлов: узел 5 (генератор ag), узел 6 (транзактно- управляемый генератор creat) и узел 7 (терминатор term со вспомогательным оператором cneg). Модельное время, в течение которого будем осуществлять эксперимент, задается переменной Protime. Если мы хотим, чтобы это время составляло 3 года, а единицей измерения установим 1 ч, то это время равно 6240,0 ч. Узел ag генерирует первый транзакт через какой-либо заданный или случайный интервал времени tbeg. Этот интервал можно назвать временем подготовки производства. Сгенерированный первый транзакт (назовем его прародителем) войдет в узел creat. В результате из этого узла выйдут Nowon + Arena* транзактов, которые поступят в узел 1.

45 Рис. 9. Структурная схема модели для решений задачи минимизации производственных затрат

46 Первые Nowon транзактов без задержки пройдут этот узел и поступят в узел 2; причем каждый транзакт-машина займет свой канал, т.е. рабочее место. Остальные Arend транзактов останутся в узле 1. Что касается транзакта-прародителя, то он выйдет из узла creat и войдет в узел term. В этом узле перед уничтожением транзакт-прародитель выполнит операцию cheg для перенастройки генератора ag. Эта операция установит время генерации второго транзакта, равное Protime. Учитывая, что с начала запуска модели уже прошло время tbeg, момент генерации второго транзакта будет равен Protime + tbeg (т.e. после выполнения модели). Это означает, что второй транзакт не выйдет из генератора во время моделирования.

47 Нетрудно представить себе, что если число наладчиков мало или равно нулю, то машины в основном будут неисправны, а затраты на производство - велики из-за простоя рабочих мест. Если же число наладчиков очень велико, то будут большими затраты на заработную плату (рис. 10, а). Но при отсутствии наладчиков машины постепенно выйдут из строя, а затраты из-за простоя рабочих мест будут велики. Что касается резервных машин, то, если их число мало или равно нулю, затраты на производство будут велики из-за простоя рабочих мест. Если же число арендуемых машин очень велико, то будут большими затраты на их аренду (рис. 10,6). Поэтому можно допустить гипотезу о наличии какой-то поверхности оптимизации затрат (рис. 10,в). Однако, учитывая, что число наладчиков и количество машин - целью величины, эта поверхность включает в себя только точки с целочисленными координатами Men и Arend.

48 Необходимо найти минимальное значение на поверхности затрат Eij и соответственно значения координат iopt (переменная Men) и jopt (переменная Arend). He будем применять сложные методы поиска экстремальных значений на поверхности оптимизации затрат. Задача решается методом перебора, так как число возможных вариантов, которые необходимо сравнить, невелико. Однако в более сложных случаях требуется применять регрессионный анализ и строить поверхность отклика 2- го порядка.

49 Рис. 10. Интерпретация возможности минимизации затрат: а - зависимость затрат от числа арендуемых машин Arend при фиксированном числе наладчиков; б - зависимость затрат от числа наладчиков Men при фиксированном числе арендуемых машин; в - предполагаемая поверхность оптимизации

50 Первое, что нужно сделать при проведении имитационного эксперимента, - это привести все времена и темпы затрат к одним единицам измерения. С помощью табл. 3 получим табл. 5 со статьями ежедневных расходов на производство. В ней три строки. Обозначим значение, указанное в каждой строке, как Ак, к = 1,2,3. Таблица 5 Статьи ежедневных расходов Ак (к=1,2,3) Номе р строк и Вид ежедневных расходов Стоимость расходов, долл. 1Дневная заработная плата одного наладчика 30 2Дневная аренда одной машины 30 3Убыток из-за простоя одного рабочего места 160

51 5. Модель «Динамическое распределение ресурсов в сети под управлением UNIX» Локальная вычислительная сеть (ЛВС) имеет один процессор, один винчестер и несколько пользовательских терминалов. Пользователи посылают операционной системе OS UNIX запросы на выполнение своих прикладных задач, которые содержат некоторое количество операций ввода-вывода (I/O). Посланный запрос поступает на обработку процессором (CPU, число процессоров может быть 1, 2,4 и др.). По каждому запросу создается процесс, который начинает работать в ЛВС. Созданному процессу OS UNIX выделяет участок памяти (RAM) и затем несколько буферов I/O (ресурсы ЛВС). Процессы, требующие меньшего объема ресурса, имеют больший приоритет при его выделении.

52 Если ресурсов ЛВС недостаточно для выполнения очередного процесса, он ожидает освобождения ресурсов другими процессами. После этого процесс получает доступ к винчестеру (HDD) для выполнения одной операции I/O. Процесс, завершивший свою работу, освобождает занимаемые им ресурсы, возвращая их операционной системе. Параметры модели приведены ниже, схема модели показана на рис. 11.

53 Параметр модели Значение параметра Количество пользователей 50 Время обдумывания запроса пользователем (среднее), с Время обработки процесса в CPU (среднее), с 0,10 Время доступа HDD-вннчестера (операция I/O): минимальное, с 0,05 максимальное, с 0,15 наиболее вероятное, с 0,10 Объем (мощность) ресурса «RAM», Мбайт 64 Объем (мощность) ресурса «Буферы I/O» 10 Вероятность завершения процесса I/O0,20 Время моделирования, ч 8

54 Рис. 11. Модель работы ЛВС под управлением OS UNIX

55 После выделения ресурса и выполнения каждой операции ввода-вывода CPU определяет очередной этап обработки процесса. Если при очередной проверке состояния процесса необходимое количество операций I/O окажется выполненным, пользователю направляется сигнал о выполнении его запроса. Получив сигнал, пользователь обдумывает и посылает следующий запрос. Описание работы модели. В модели транзакт имеет несколько состояний: обдумываемый пользователем запрос; досланный пользователем запрос на организацию процесса ввода-вывода; созданный OS UNIX-V процесс I/O; процесс, получивший RAM; процесс, получивший буферы I/O; процесс, выполнивший хотя бы одну операцию I/O; сигнал пользователю о выполнении запроса.

56 Текущее состояние транзакта заносится в его параметр и является критерием перехода из сервера, имитирующего работу процессора. Изменение состояния транзакта осуществляется в узлах, имитирующих терминалы пользователей, работу CPU, обращение к HDD, выделение RAM и буферов I/O. Имена используемых параметров транзактов переопределены операторами #defme. Модель является замкнутой. С помощью схемы зарядки (узлы 9, 10 и 11) генерируется семейство транзактов количеством USERS (число терминалов пользователей). Сгенерированные транзакты поступают в сервер 8, имитирующий пользовательские терминалы. Число каналов сервера равно константе USERS. Время обслуживания - это время обдумывания пользователем запроса к OS UNIX. В соответствии с центральной предельной теоремой оно распределяется по нормальному закону.

57 Перед входом в узел состояние транзакта (параметр PROCESS_STATUS) получает значение RAM. Это значит, что транзакт, выйдя из сервера, будет играть роль запроса на создание процесса I/O и ему нужно будет выделить память. Из сервера транзакт переходит в очередь 1, имитирующую очередь к процессору, и далее в сервер 2, имитирующий процессор. Перед входом в узел на основании значения параметра PROCESS_STATUS принимается решение о дальнейшей обработке транзакта. Если состояние транзакта имеет значение HDD, это значит, что данный транзакт имитирует уже созданный процесс I/O. Если при этом установлен флажок доступа к винчестеру (параметр транзакта HDD_ACGESS равен 1), значит, процесс выполнил хотя бы одну операцию ввода-вывода и с вероятностью 0,2 является завершенным (значение генератора rundum больше 0,8). В этом случае транзакт отправляется в сервер 8 -OS UNIX посылает пользователю сигнал о выполнении его запроса.

58 Транзакт отправляется в очередь 3, имитирующую очередь запросов к винчестеру на выполнение операций I/O (процесс должен выполнить очередную операцию ввода- вывода), когда выполняется хотя бы одно из двух условий: значение HDD_ACCESS равно нулю (процесс не выполнил еще ни одной операции ввода-вывода); значение генератора rundum не превосходит 0,8 (процесс не завершен); если состояние транзакта имеет значение BUFF, это значит, что процесс получил память и теперь должен получить буферы I/O. В этом случае транзакт отправляется на «склад ресурсов» - узел attach 7, имитирующий выделение процессу буферов I/O.

59 Число каналов сервера 2, имитирующего процессор, равно единице, поэтому процессор одновременно обрабатывает один процесс (запрос). Время обслуживания - это время обработки ветви программы процесса между двумя прерываниями; оно имеет равномерный закон распределения. Если транзакт входит в узел CPU, а далее будет отправлен в сервер 8 (это ответ на терминал пользователя), то сбрасывается флажок доступа к винчестеру HDD_ACCESS (он получает значение нуль), а транзакт возвращает полученные им ресурсы (выполняются операции detach). В очереди 3 устанавливается флажок доступа к винчестеру HDD_ACCESS, и транзакт переходит в сервер 4, имитирующий обращение к винчестеру - операцию I/O. Число каналов сервера равно единице, поэтому одновременно выполняется одна операция I/O.

60 Время обслуживания в сервере есть время доступа к винчестеру; оно задано с помощью треугольного beta- распределения с минимальным временем доступа 0,05 с, максимальным - 0,15 с и наиболее вероятным - 0,1 с. Из этого сервера-винчестера транзакт переходит в очередь 1. На «складе ресурсов» 5 транзакту выделяется ресурс типа RAM. Перед входом в узел определяется объем запрашиваемой памяти (параметр транзакта RAM_AMOUNT) как случайная величина, составляющая какую-то долю от мощности - максимального объема памяти. Состояние транзакта PROCESS_STATUS получает значение BUFF, в результате процесс помечается как ожидающий выделения буферов I/O. В описании узла параметр RAM_AMOUNT указывает количество запрашиваемого ресурса. Флажок учета приоритета запроса установлен. Далее транзакт переходит к диспетчеру ресурсов 6, а оттуда в очередь 1.

61 На «складе ресурсов» 7 транзакту выделяется требуемое количество буферов I/O (параметр транзакта BUF_AMOUNT). Состояние транзакта получает значение HDD (процесс готов к выполнению операций I/O, доступ к винчестеру разрешен). В остальном узел 7 аналогичен узлу 5. Анализ результатов моделирования. За период моделирования к OS UNIX поступило запросов; столько же было образовано вычислительных процессов (счетчик входов узла 8). Каждый пользователь послал в среднем 500 запросов, и 7 пользователей на момент окончания моделирования обдумывают очередной запрос. В среднем пользователи потратили на обдумывание запросов 18 % времени работы с ЛВС. Было запросов от процессов на получение RAM (счетчик входов узла 5). Однако 41 запрос остался неудовлетворенным: из-за дефицита 41 транзакт находится в очереди attach на момент завершения моделирования.

62 Всего процессор находился в состоянии работы 71,3% времени, операции ввода-вывода заняли 44,6% времени. Среднее время ожидания процессора составило 0,05 с, а среднее время пребывания в очереди ввода-вывода - 0,01 с. Средняя загрузка ресурса RAM составляет 73,8%. На момент окончания моделирования остаток ресурса RAM составляет 9 Мбайт, а стоящие в очереди процессы требуют в общей сложности 2495 Мбайт памяти. Это объясняет большую задержку - 42,33 с при каждом обращении; данный ресурс - самое узкое место ЛВС. Средняя загруженность ресурса буферов I/O составляет 55,5%. Задержка в связи с возможным дефицитом сравнительно невелика - 0,38 с.

63 Контрольные вопросы: Если маршрут, построенный с помощью алгоритма «ближайшего непосещенного города» выглядит практически без изъянов, он действительно короткий, на нем нет петель (т.е. выполнено необходимое условие построения оптимального маршрута), то чем он плох? Какой узел имитационной модели может замерять пройденный путь? Как вертолет коммивояжера в модели ориентируется на местности? Можно ли в модели учитывать времена задержек коммивояжера в населенных пунктах? В постановке задачи взаимного исключения нет ни одного узла обслуживания, требующего на обслуживание какое-то время. Тогда почему возникают значительные задержки транзактов?

64 Каков принцип работы узла типа delet в модели «Стоянка маршрутного такси» при решении задачи взаимного исключения? Почему интервал прихода пассажиров в модели «Стоянка маршрутного такси» имеет экспоненциальное распределение? Почему интервал генерации транзактов, имитирующих такси в этой же модели, имеет нормальное распределение? Какая основная причина мешает в модели «Эффективность компьютеров в автоматизированной бухгалтерии» применить формулу Поллачека-Хинчина? Какое решение обойдется дешевле для бухгалтерии: купить более надежный компьютер или поставить еще один аналогичный с той же надежностью?

65 Когда интервал времени между двумя последовательными поступлениями документов из двери рабочей комнаты на столик в компьютерном зале распределен не по экспоненциальному, а по нормальному закону? По каким причинам компьютеры не смогут справиться с обработкой документов в модели «Эффективность компьютеров в автоматизированной бухгалтерии»? Почему аналитическое решение в модели «Минимизация производственных затрат» неправильное? Какие две основные тенденции показывают на возможность существования оптимального решения в этой модели? Почему поиск оптимального решения методом перебора начинается со значений i = 3, j = 3, с последующим увеличением i и j в модели «Минимизация производственных затрат»?

66 Для чего нужна схема зарядки в этой модели? Почему модель «Динамическое распределение ресурсов в сети под управлением Unix» нельзя исследовать как замкнутую стохастическую сеть? Что показывает «состояние транзакта» в параметре t iu0? Какие ресурсы влияют на быстродействие работы в вычислительной сети не меньше, чем быстродействие и число процессоров на сервере? Для чего нужна схема зарядки в модели «Динамическое распределение ресурсов в сети под управлением Unix»?