Использование стохастического описания реализации производственных процессов позволяет решать задачу нахождения математического ожидания и стандартного отклонения директивного времени на их реализацию в условиях неопределенности. Для сложного производственного или технологического процесса директивное время рассматривается как случайная величина, описываемая подходящей функцией распределения с конечным математическим ожиданием и дисперсией [3]. Для получения дисперсионных оценок необходимы некоторые предположения, касающиеся стохастических характеристик каждого элемента производственной структуры в стандартных и, если необходимо, аварийных условиях.
Согласно [4], узлы стохастической сети могут быть интерпретированы как состояния процесса, а дуги - как переходы из одного состояния в другое. Такие переходы можно рассматривать как реализацию обобщенных операций производства, характеризуемых плотностью распределения, или функцией массы, и вероятностью выполнения. Таким образом, получается стохастическое графовое представление производственных процессов, где узлы являются входом и выходом для операций. Дуги характеризуют время выполнения реальной производственной операции.
Каждый внутренний узел стохастической сети выполняет две функции, одна из которых касается входа в узел, а другая - выхода. Обычно эти функции называют входной и выходной. В [4] определен следующий тип входной функции:
узел выполняется, если выполнена дуга, входящая в него, при условии, что в заданный момент времени может выполняться только одна дуга.
Для выхода там же определены два типа выходной функции:
детерминированный выход и вероятностный выход. Для детерминированной выходной функции характерна ситуация, когда все дуги, выходящие из узла, выполняются, если этот узел выполнен. Для вероятностной выходной функции - ровно одна дуга, выходящая из узла, выполняется, если узел выполнен. Выбор такой дуги может быть описан с помощью функции распределения вероятностей.
Электронный журнал ИССЛЕДОВАНО В РОССИИ 2572 Процедура, предлагаемая для анализа вероятностных характеристик процесса реализации производственных операций, базируется на стохастических сетях типа GERT (Graphical Evaluation and Review Technique) [5] и совмещает теорию потоков в графах, функции генерации момента и PERT - анализ [4] для получения результата.
Можно установить следующие соотношения между сетями PERT-типа, графами с потоками и стохастическими сетями.
1. Сети PERT-типа соответствуют стохастическим GERT-сетям, в которых все узлы являются детерминированными "И"-узлами (And - Deterministic).
2. Графы с потоками соответствуют стохастическим сетям с простым мультипликативным параметром (все аддитивные параметры, такие как время, установлены равными нулю). Вероятностная интерпретация для мультипликативного параметра исключена из рассмотрения.
Рассмотрим реализацию производственного процесса (программы), заданную сетью G=(N, A), которая содержит только GERT-узлы, образующие множество N.
Пусть время выполнения операции (i, j) (для стохастической интерпретации реализации) есть случайная величина Hij. По определению (i, j) может быть выполнена только в том случае, если выполнен узел i. Поэтому для изучения вопросов, связанных с выполнением этой операции, необходимо знать условную вероятность (в дискретном случае) случайной величины Hij при условии, что узел i выполнен. Это, в свою очередь, позволит нам провести исследования, связанные с выполнением всей сети. В частности, можно определить моменты распределения времени выполнения сети, с помощью которых будут вычислены математическое ожидание и дисперсия времени выполнения сети.
В [4] описаны некоторые наиболее важные функции распределения и указаны соответствующие производящие функции моментов и первые (математические ожидания) и вторые центральные моменты.
Пусть pij - вероятность того, что операция (i, j) будет выполнена при условии, что узел i выполнен. Для случайной величины Hij определим W-функцию как Wij (s) = pij M (s), где Mij - условная производящая функция моментов случайной ij величины Hij [4].
С помощью этого преобразования всегда можно определить сеть G', структура которой идентична структуре сети G, только вместо двух параметров дуг pij и Hij присутствует один параметр Wij.
Согласно системе GERT-моделей можно включить в описание в качестве параметра дуги время выполнения соответствующей операции производственного процесса. Однако в действительности можно рассматривать любой характерный параметр, обладающий аддитивностью по дугам любого пути. Если времена выполнения сети G представляются независимыми случайными величинами, то для G' справедливы правила эквивалентности (с вычислительной точки зрения):
а) для случая, когда G' состоит из двух последовательных дуг;
б) G' - из двух параллельных ветвей;
в) G' - из одной ветви и одной петли.
Электронный журнал ИССЛЕДОВАНО В РОССИИ 2573 Перечислим основные этапы решения поставленных задач моделирования посредством стохастических сетей типа GERT.
1. Перевести качественное описание рассматриваемых производственных процессов в GERT-сетевую модель.
2. Собрать необходимые структурные данные для описания связей в GERTсети (включая характеристики операций).
3. Применить топологическое уравнение для определения эквивалентной функции (или функций) GERT-сети.
4. Вычислить через эквивалентную функцию две следующие характеристики функционирования сети:
- вероятность выполнения конкретного узла;
- функции генерации момента для времени, связанного с узлом, если он выполняется.
5. Сделать выводы относительно реализации производственных процессов на основе информации, полученной в пункте 4.
В нашем случае важно рассмотреть подход к минимизации времени и максимизации прибыли при анализе реализуемости с учетом стохастической реализации производственного процесса. Первоначально рассмотрим простой ациклический детерминированный процесс, который имеет "GERT-подобную узловую логику" [6]. Такую модель будем называть сетью для планирования.
Термин "планирование/решение" показывает, что осуществляется процесс выбора, т.е. принимается решение о том, какие операции должны быть выполнены для оптимизации некоторой целевой функции. Это приводит к задаче комбинаторной оптимизации, частным случаем которой является, например, "decision CMP" - метод критического пути для случая, когда присутствуют только узлы двух типов AND и OR [7].
Для учета вероятностных характеристик реализации вводится понятие случайных акций и рассматривается возможность многоразовой последовательной реализации до момента успешного завершения.
Пусть N - ациклическая сетевая модель с источниками и стоками (действия, соответствующие операциям процесса, представляются дугами), где множество узлов обозначается V, а множество дуг - E. Предположим, что N имеет только один исток, который обозначается через r и соответствует началу рассматриваемого процесса. Предполагается также, что один из стоков N представляет собой успешное завершение и обозначается s. Оставшиеся стоки, если они есть, могут представлять собой различные виды неудачного завершения или прерывания. Будем использовать следующие Определения, введенные в [8].
Определение (1) Ациклическую сетевую модель - N(V,E) только с одним истоком и со стоками назовем сетью для планирования/решения, если каждый узел i из N определен через входную характеристику X 0,1,..., P(i) и выходную i + характеристику X 0,1,..., S(i), где множество узлов обозначается V, а множество i дуг - E; P(i), S(i) - мощность множества предшественников и последователей узлов I Электронный журнал ИССЛЕДОВАНО В РОССИИ 2574 соответственно. Эти характеристики, формирующие узловую логику, имеют следующие значения:
(1а) - узел активируется сразу же, как только входные действия Xi- завершаются;
(1б) - как только узел i активирован, то не более Xi+ выходных действий начинает выполняться, если узел i не активируется, то ни одно выходное действие не выполняется.
Иногда уместно заменить термин "не более" на "точно". Два условия из Определения (1) подразумевают, что каждое действие выполняется сразу, как только это становится возможно.
Для источника r полагаем X = 0, т.е. он всегда активирован. Кроме того, r Xi+ = 0 для i S, где S - множество стоков N.
Нужно отметить, что, во-первых, если Xi- =1, тогда узел i имеет OR-вход, и, если Xi- = P(i), то тогда i имеет AND-вход. И если "не более" заменяется на "точно" в (1б), то Xi+ =1 соответствует вероятностному выходу, а Xi+ = S(i) соответствует детерминированному выходу. Во-вторых, если данная сеть N для решения имеет множество источников R (|R|>1) и множество R' R, R' 0 активизируется в начале выполнения процесса, то можно формально перевести N в соответствующую одноистоковую сеть следующим образом. Введем новый единый источник r0 и для каждого i R введем вспомогательную дугу< r0,i >.
Кроме того, установим Xr- = Xr+ = 0. Определим 0 для Xi- =.
1 для i R' i R / R' Введем дуговые переменные ( E ):
i, j если выполняется;
1, wij = 0, иначе;
и узловые переменные (i V ):
i если активируется;
1, uij = 0, иначе;
где ur =1, т.е. источник всегда активируется.
Тогда два условия узловой логики, введенные в Определении (1) могут быть переписаны в следующем виде:
wki Xi-ui (i V /{r}); (1) kP (i) wki < Xi- + Miui, (2) kP (i) где Mi > P(i) - Xi- (i V /{r});
Электронный журнал ИССЛЕДОВАНО В РОССИИ 2575 + X ui (i V / S). (3) wij i jS(i ) Если в (1) ui =1, то это значит, что в результате активации узла i выполняется по крайней мере Xi+ входящих воздействий. Если ui =0 в (2), то это значит, что узел i не активирован, так как менее чем Xi+ входящих действий выполнено. Таким образом, оба эти неравенства вместе соответствуют (1а). Полагая ui =1 и ui =0 в (3), обеспечиваем выполнение второго условия узловой логики.
Определение (2) Так как сетевая модель ациклична, то каждая операция соответствующего процесса либо выполняется только один раз, либо не выполняется вообще. Таким образом, каждая реализация процесса (или реализация сетевой модели) может быть соотнесена с множеством выполняемых действий сети или с функцией w : E {0,1}; ( i, j E), значения которой задаются как i, j выполняется если 1, w( i, j )=: wij = 0, иначе.
С другой стороны, если некоторая реализация сети w задана, то, как узловые, так и дуговые переменные для этого случая также специфицированы. Будем говорить о допустимой реализации, если w удовлетворяет условиям узловой логики. Тогда e*= {w : E {0,1}wij -удовлетворяет (1)-(3); E } и e* - множество всех допустимых реализаций.
Итак, стохастическое представление моделей реализуемости в виде GERTсети позволяет получить достаточное количество полезной информации о временных характеристиках протекания производственных процессов. Используя неравенство Чебышева, можно показать пределы, в которых будет изменяться фактическое время реализации процессов, а также получить более сильные утверждения. Кроме того, если реальный показатель не соответствует этим оценкам, то, как показано в [8], можно построить ряд критериев для проверки гипотез, позволяющих определить перспективные характеристики времени реализации процессов.
В заключение необходимо отметить, что рассмотренная GERT-сетевая модель, является своеобразной альтернативой традиционным методам определения директивных времен реализации производственных процессов. При использовании традиционных методов предполагается, что время выполнения каждой отдельной операции постоянно. После суммирования этих времен в полученный результат вносится некоторая поправка с целью учесть случайные колебания или устранить неустойчивость действительных времен. Система же GERT-моделей позволяет включать случайные отклонения и неопределенность, возникающие непосредственно во время выполнения каждой операции производственного процесса. Следовательно, в полученный результат уже включены все случайные колебания и нет необходимости вносить в него дополнительные поправки, не считая тех, которые соответствуют аварийным производственным ситуациям.
Электронный журнал ИССЛЕДОВАНО В РОССИИ 2576 Литература 1. Дегтерев А.С., Нейман Г.А. Моделирование задачи оптимизации загрузки технологического оборудования // Экономика и финансы, 2002, № 20 (22), с. 46-49.
2. Дегтерев А.С., Ерыгин Ю.В., Лобков К.Ю. Оптимизация портфеля инновационных проектов на машиностроительном предприятии ОПК в условиях конверсии // Конверсия в машиностроении, 2004, № 3 (64), с. 83-85.
3. Корячко В.П., Шибанов А.П., Шибанов В.А. Численный метод нахождения закона распределения выходной величины GERT-сети // Информационные технологии, 2001, № 7,.с. 16-21.
4. Филлипс Д., Гарсиа-Диас А. Методы анализа сетей. М.: Мир, 1984.
5. Шибанов А.П. Нахождение закона распределения выходной величины GERT-сети большой размерности // Информационные технологии, 2002, № 1, с. 4245.
6. Antamoshkina O., Kovalev I. Modelling, Optimization and Computer-Realization of Control Cyclograms. SAA-Press. Krasnoyarsk, 1996.
7. Neumann K. Stochastic Project Networks. Springer-Verlag, 1990.
8. Kovalev I. System of Multi-Version Development of Spacecrafts Control Software. Pro Universitate Verlag, Sinzheim, 2001.
Книги по разным темам