Ger t-анализ формирования технологических циклов управления космическими аппаратами

Таким образом, в данной статье был рассмотрен метод DEA, предложена и построена граница отставания, а также представлено разбиение множества объектов на два класса, что облегчает ЛПР классификацию объектов, так как предполагает существование всего двух классов: плохого и хорошего. Это удобно, когда выборка объектов достаточно велика. Если ЛПР требуется точное количественное значение эффективности, то он получает это решение из прямого решения задачи по методу DEA.

/

Рис. 3. Графики границы эффективности

Библиографический список

1. Моргунов, Е. П. Многомерная классификация на основе аналитического метода оценки эффективности сложных систем / Е. П. Моргунов ; НИИ систем упр., волновых процессов и технологий. Красноярск, 2003.

2. Coelli, T. An introduction to efficiency and productivity analysis / T. Coelli, D. S. Prasada Rao, G. E. Battese. Boston : Kluwer Academic Publishers, 1998.

Рис. 4. Классы объектов, полученные разбиением совокупности объектов с помощью разделяющей прямой

T. A. Rukavitsyna MODEL DEVELOPMENT OF THE DEA METHOD

DEA (Data Envelopment Analysis) method and its modification are considered. An inefficient frontier is suggested to be introduced. The intersection of the efficient frontier and inefficient frontier are investigated. A hyper plane going through the intersection of both frontiers is constructed to divide the sample into two classes.

Keywords: DEA, efficiency, efficientfrontier.

©Рукавицына Т. А., 2009

1

УДК 681.34

А. А. Ступина, Е. С. Разгулина

GERT-АНАЛИЗ ФОРМИРОВАНИЯ ТЕХНОЛОГИЧЕСКИХ ЦИКЛОВ УПРАВЛЕНИЯ КОСМИЧЕСКИМИ АППАРАТАМИ

Предложено использовать аппарат сетевого анализа, включая GERT-анализ, для решения задач оптимизации формирования технологических циклов управления (ТЦУ).

Ключевые слова: GERT-сети, GERT- анализ, ациклическая сетевая модель технологических циклов управления, оптимизационная модель, вероятностные характеристики ТЦУ, время реализации ТЦУ.

Решение задачи формирования технологических циклов времени ответа центра управления полетами (ЦУП) на за-

управления (ТЦУ) связано в первую очередь с необходи- прос от объекта. Ограничения на время реакции связыва-

мостью обеспечить взаимодействие с объектом (космичес- ются в этом случае с выполнением периодических действий

ким аппаратом) в соответствии с временными циклограм- в рамках технологического цикла управления.

мами в режиме реального времени. Более простая форма Естественной математической интерпретацией рас-

режима реального времени предполагает лимитирование пределенных, асинхронных, мультипроцессорных и муль-

типрограммных систем являются сетевые модели, которые позволяют отражать распределенность структуры, сетевой характер взаимосвязей между процессами и ресурсами, аппаратными и программными компонентами средств управления. В связи с этим для решения задач оптимизации формирования ТЦУ целесообразно привлекать аппарат сетевого анализа, включая GERT-анализ (GERT - аббревиатура от graphical evaluation and review technique). В последнее время GERT-сети получают все большее распространение для моделирования и оптимизации технических систем, к числу которых можно отнести ТЦУ космических аппаратов.

Общность методов построения моделей ТЦУ для анализа и проектирования, выражающаяся в использовании сетевых моделей, позволяет создать единые инструментальные средства оптимизации формирования ТЦУ.

Сетевое представление ТЦУ открывает широкие возможности для использования в моделях анализа реализуемости циклограмм математического аппарата, учитывающего вероятностную интерпретацию. Одна из таких возможностей связана с введением неопределенности в продолжительность реализации задач ТЦУ (аналогично использованию PERT- методики в сетевом анализе). Однако такая модель представляет собой типичную детерминированную сеть, для полного выполнения которой необходимо выполнение всех дуг, т. е. безоговорочное выполнение соответствующих операций ТЦУ Из этого условия следует, что в такую модель не могут быть включены операции с обратной связью, поскольку они представлены петлями, существование которых в свою очередь означает, что конечный узел операции должен быть выполнен раньше ее начального узла.

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

Согласно [1], узлы стохастической сети могут быть интерпретированы как состояния процесса, а дуги - как переходы из одного состояния в другое. Такие переходы можно рассматривать как реализацию обобщенных задач ТЦУ, характеризуемых плотностью распределения, или функцией массы, и вероятностью выполнения. В результате получается стохастическое графовое представление ТЦУ, где узлы являются входом и выходом для очередной задачи, а дуги характеризуют время выполнения реальной технологической операции.

Каждый внутренний узел стохастической сети выпол -няет две функции, одна из которых касается входа в узел, а другая - выхода. Обычно эти функции называют входной и выходной.

В [1] определен следующий тип входной функции: узел выполняется, если выполнена дуга, входящая в него, при условии, что в заданный момент времени может выполниться только одна дуга.

Для выхода там же определены два типа выходной функции: детерминированный выход и вероятностный

выход. Для детерминированной выходной функции характерна ситуация, когда все дуги, выходящие из узла, выполняются, если этот узел выполнен. У вероятностной выходной функции ровно одна дуга, выходящая из узла, выполняется, если узел выполнен. Выбор такой дуги может быть описан с помощью функции распределения вероятностей.

При анализе реализуемости ТЦУ стохастическая сеть определяется как сеть, которая работает только при выполнении некоторого подмножества дуг, при этом время выполнения каждой дуги (задачи ТЦУ) выбирается в соответствии с вероятностным распределением. В такой стохастической сети для выполнения узла нет необходимости в выполнении всех дуг, входящих в него. Поэтому в такой модели допускается существование циклов и петель.

Для более наглядного представления детерминированного случая можно обратится к простому ациклическому детерминированному ТЦУ, который имеет GERT-подобную узловую логику. Такая модель соответствует сети для проектирования (решения) ТЦУ Это указывает на то, что осуществляется выбор, т. е. принимается решение о том, какие задачи ТЦУ должны быть выполнены для минимизации некоторой целевой функции. Для учета вероятностных характеристик реализации ТЦУ вводится понятие случайных акций и рассматривается возможность многоразовой последовательной реализации ТЦУ до момента успешного завершения.

Пусть N - ациклическая сетевая модель ТЦУ с источниками и стоками (действия, соответствующие задачам ТЦУ, представляются дугами), где множество узлов обозначается V, а множество дуг - Е. Предположим, что N имеет только один исток, который обозначается через г и соответствует началу рассматриваемого ТЦУ Также предположим, что один из стоков N представляет собой успешное завершение ТЦУ Обозначение этого стока - 5. Оставшиеся стоки, если они есть, могут представлять собой различные виды неудачного завершения или прерывания процесса управления.

Определение 1. Ациклическую сетевую модель ТЦУ ЩУ, Е) только с одним истоком и со стоками назовем сетью для проектирования (решения) ТЦУ, если каждый узел i из N определен через входную характеристику СИ-. е {0, 1, ...,|Р(0|} и выходную характеристику СИ+.е{0, 1, ...,|£(0|}, где множество узлов обозначается V, а множество дуг - Е; |Р(0|, |£(.)| - мощность множества предшественников и последователей узлов I соответственно. Эти характеристики, формирующие узловую логику, имеют следующие значения:

а) узел активируется сразу же, как только входные действия СИ-. завершаются;

б) как только узел . активирован, то не более СИ+. выходных действий начинает выполняться. Если узел . не активируется, то ни одно выходное действие не выполняется. Иногда уместно заменить сочетание «не более» на «точно».

Два условия из определения 1 подразумевают, что каждое действие выполняется сразу как только это становится возможно.

Для источника г полагаем СИ-г = 0, т. е. он всегда активирован. Кроме того, СИ+. = 0 для . е £ где £ - множество стоков N.

Заметим, что если CH+. = 1, то узел i имеет OR-вход, а если CH-r = |P(i)|, то тогда i имеет AND-вход. И если в определении 1б «не более» заменяется на «точно», то CH+. = 1 соответствует вероятностному выходу, а CH+. = |S(i)| - детерминированному выходу. Если же заданная сеть N для формирования ТЦУ имеет множество источников R (|R |>1) и множество R'H R, R' Ф 0 активизируется в начале выполнения ТЦУ, то можно формально перевести N в соответствующую одноистоковую сеть.

Дуговые переменные обозначаются через w., при этом полагается, что w.. = 1, если (i, j) выполняется ((i, j) є E), и w.. = 0 в противном случае. Узловые переменные п.. = 1 (i є V), если i активируется, иначе п.. = 0, где ыг = 1, т. е. источник всегда активируется. Тогда условия узловой логики, введенные в определении 1, могут быть представлены в следующем виде:

X w. > CH-.u;

кєР

X wu < CH. + Mtu;

кєР

где M > I P( i) - CH-t |,

X w. < CH+.u;

.є8

Так как решающая модель ациклична, то каждая задача соответствующего ТЦУ либо выполняется только один раз, либо не выполняется вообще. Таким образом, каждая реализация ТЦУ может быть соотнесена с множеством выполняемых задач ТЦУ или с функцией w: E  {0, 1}; ((i, j) E), значения которой задаются как w((i, j)) = : w. = 1, если (і, j) выполняется, и 0, если иначе.

С другой стороны, если некоторая w-я реализация ТЦУ задана, то как узловые, так и дуговые переменные для этого случая специфицированы и можно говорить о допустимой реализации ТЦУ, если w удовлетворяет условиям узловой логики. Тогда e = {w: E  {0, 1}| w. удовлетворяет условиям (1)-(3); (i, j) є E } и e - множество всех допустимых реализаций ТЦУ.

Если в решающей сети для формирования ТЦУ обозначить вес дуги (і, j) є E как длительность d.. є R+ соответствующей задачи ТЦУ, то dW - длительность w-й реализации ТЦУ, т. е. время, требуемое для исполнения всех задач (i, j) при w. = 1. Необходимо минимизировать С^^при условиях: w активирует s; (w є є). Через є*: = {w є є | w активирует s} обозначим множество успешных реализаций ТЦУ Для є* Ф 0 d* = min dW соответствует минимальному значению целевой функции задачи.

Приняв, что каждая реализация ТЦУ начинается с задействования истока r в момент 0, считаем для w є e, что tw есть время активации узла j є V для w-й реализации, причем tJ = 0 и tw = ¥, если j не активируется в течение w-й реализации ТЦУ Для j є V\{r} имеем tw = min{t > 0 | существуют CH-, отличные от i є P(j) и такие, что w.. = 1 и tw + d.. < t}. 1 'J

j V ’

Кроме того, справедливо dW >tsw V, п'єс. Тогда dW > tsw, если сток s активируется, пока некоторые действия (i,

s * w

j) при w.. = 1 все еще выполняются. Так как ґ. = mint -

J w j wü

самое раннее из возможных времен задействования узлаj в течение какой-либо возможной реализации ТЦУ, то te = min max {(tw + d ) w },

s ms‘ (i,j)є£ . Г iJ

Е„={(/, .) е Е | =1},

где Ем - множество задач, выполняемых в течение м-й реализации ТЦУ

Учитывая вышесказанное, можно утверждать, что минимальная длительность успешной реализации ТЦУ равна самому раннему моменту задействования стока 5: й* = Р для е* ф 0. Таким образом, можно найти величину йР, вычисляя минимально возможные моменты задействования а е V).

Рассматривая моменты I. как компоненты вектора временной развертки (ВВР) ТЦУ, которые удовлетворяют условию (/. - /. - й..)м. > 0, (/,у) е Е; (. > 0, . е И{г}; / = 0, можно утверждать, что для некоторой м-й реализации ТЦУ (м е е *) моменты I м отвечают этим ограничениям, а самые ранние моменты Iе удовлетворяют им для соответствующей минимальной м е е*. Соответствующая оптимизационная задача имеет вид: минимизировать тах {(/.м + й.) м. } при выполнении условий (1)-(3) и ограничений ВВР (и = 0).

По сравнению с минимизацией затрат на выполнение ТЦУ здесь имеются не только добавочные ограничения, но и более сложная целевая функция. Кроме того, дополнительно к узловым и дуговым переменным существуют моменты /.; (/ е V), которые также являются переменными данной оптимизационной задачи.

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

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

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

Процедура, используемая для анализа вероятностных характеристик ТЦУ, базируется на стохастических сетях типа GERT и совмещает теорию потоков в графах, функции генерации момента и PERT-анализ для получения результата. Стохастическое представление ТЦУ в виде GERT-сети позволяет получить достаточное количество полезной информации о временных характеристиках реализации ТЦУ Используя неравенство Чебышева, можно показать пределы, в которых будет изменяться фактическое время реализации ТЦУ, а также получить более сильные утверждения. Кроме того, если реальный показатель не соответствует этим оценкам, то строится ряд критериев для проверки гипотез, позволяющих определить перспективные характеристики времени реализации ТЦУ

Рассмотренная GERT- модель является своеобразной альтернативой традиционным методам определения директивных времен реализации ТЦУ, поскольку система моделей позволяет учесть случайные отклонения и неопределенность, возникающие непосредственно во время выполнения каждой отдельной задачи ТЦУ Следовательно, в полученный результат уже включены все случайные колебания и нет необходимости вносить в него дополнительные поправки, не считая тех, которые соответствуют нештатным ситуациям или авариям.

Анализ моделей формирования ТЦУ показывает, что они имеют, как правило, алгоритмически заданный критерий, а часто - и ограничения. Очевидно, что при такой постановке задачи имеет смысл говорить только об эвристических процедурах случайного поиска.

Формально указанные оптимизационные модели могут быть сведены к задаче псевдобулевой оптимизации вида

Определим инверсию р (1 < . < п) как подстановку на множестве ВП. И пусть X, У е В2П. Тогда направлением

F(X)  extr,

XES

К(Х, У) с ВП назовем множество таких точек 2 є В' которые можно представить в виде

2 = п. -п. -...-п. X,

, Pjt, такое, что

£ = |X е В^ : X = т, т < п }

на решение которой ориентированы алгоритмы метода изменяющихся вероятностей (МИВЕР) [2].

Алгоритмы МИВЕР не укладываются в традиционную схему алгоритмов случайного поиска для непрерывной оптимизации, которая в общем случае включает два этапа:

- первый этап - в соответствии с некоторым правилом и вероятностным распределением на множестве всех возможных направлений поиска случайным образом выбирается одно из направлений;

- второй этап - по заданному алгоритму выбирается длина шага перемещения поиска в выбранном направлении.

Для реализации этой схемы в пространстве булевых переменных необходимо определить понятия направления в В2П.

где п. е Щ(X,У), ЩХ, У) = pk,Р2

•)ш

У = п. -п. -...%. X (к < п).

¡1 12 1к

В соответствии с этим определением для любой точки X е В2П существует (2п - 1) направлений, каждое из которых однозначно определяется любой точкой У е В2П либо множеством Щ(А", У). Длину шага по направлению К(X, У) определяем, используя метрику Хэмминга.

Для введенных таким образом понятий направления и длины шага общая схема построения алгоритмов направленного поиска выглядит следующим образом.

1. Случайным образом выбирается начальная точка поиска X0. В соответствии с заданным алгоритмом выбирается точка У и формируется множество инверсий Щ(X0, У), определяющих выбираемое направление.

2. По заданному правилу определяется длина шага I < саМЩ^,^). Формируется множество инверсий Щ(X0,2) с Щ(^0, У), саМЩ^, 2)) = I. Осуществляется переход в точку

2 =п. -п. -... -п. ,п. е Щ(X0,2).

¡1 12 Ц ’ 1] ' > '

Критерием остановки работы алгоритма, как и при случайном поиске на множестве Я”, может быть ограничение по числу итераций, определение точки локального минимума и т. п. Предложенная схема четко подразделяется на два этапа: определение направления поиска и выполнение шага в выбранном направлении.

Приведенные выше процедуры направленного случайного поиска были реализованы программно и прошли апробацию при решении задач формирования ТЦУ Численные результаты показали, что в ряде случаев алгоритмы направленного поиска позволяют получить лучшие результаты по сравнению с алгоритмами МИВЕР.

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

Сепарабельные псевдобулевые функции определяются следующим образом:

А^) = Ахг •••,х”) = X &(X^

1< ] <п

где Д: В 2 ”Я \

Оценкой сверху числа необходимых вычислений значений целевой функции, когда поиск начинается в точке (1, •, 1) (самой дальней от минимума), будет п2 - п + 2 вычислений, а оценкой в среднем (по выбору начальной точки поиска) будет величина (п2 + 4)/2 - 1/2п. Данные оценки позволяют сравнить эффективность локального поиска на классе сепарабельных псевдобулевых функций с эффективностью других подходов, например эволюционных алгоритмов.

Рассмотренный выше алгоритм локального поиска реализует информационную сложность класса слабо немонотонных псевдобулевых функций. При оптимизации произвольных функций, в том числе и унимодальных, он не дает никаких гарантий за исключением полного перебора.

Библиографический список

1. Филлипс, Д. Методы анализа сетей / Д. Филлипс, А. Гарсия-Диас. М. : Мир, 1984.

2. Kovalev, I. Optimization Problems when Realizing the Spacecrafts Control / I. Kovalev // Advances in Modeling and Control. 1998. Vol. 52, №> 1, 2. P. 63-70.

A. A. Stupina, E. S. Razgulina

GERT-ANALYSIS OF THE TECHNICAL CONTROL CYCLES FORMATION FOR SPACE CRAFTS

In this paper authors recommend to use network analysis method (including GERT analysis) to optimize formation of technical control cycles (ТСС).

Keywords: GERT-networks, GERT-analysis, acyclic network model for TCC, optimization model, TCC relative parameters, time for TCC accomplishment.

© Ступина А. А., Разгулина Е. С., 2009