УДК 519.876.2
С.П. Вовк, Л. А. Гинис МОДЕЛИРОВАНИЕ ПЕРЕХОДОВ МЕЖДУ ЭТАЛОННЫМИ СИТУАЦИЯМИ В СЛОЖНЫХ СИСТЕМАХ В УСЛОВИЯХ НЕОПРЕДЕЛЕННОСТИ*
Рассмотрены вопросы моделирования процедуры выбора решения по управлению динамическими подсистемами сложных организационных систем. Предложен способ нахождения переходов между эталонными ситуациями в многоуровневых сложных системах, основанный на сравнении нечетких интервалов. Система описана нечеткой ситуационной сетью. Рассмотрены два способа описания признаков в эталонных ситуациях, определяющие переходы между эталонными ситуациями: описание признаков с помощью нечетких чисел и с помощью нечетких интервалов. Определены проблемы и приведены способы поиска решений в нечеткой иерархической системе управления.
Ложная система; управляющее решение; неопределенность; нечеткие числа; нечеткие интервалы; эталонная ситуация; нечеткий вывод.
S.P. Vovk, L.A. Ginis
MODELLING OF TRANSITIONS AMONG ETALON SITUATIONS
IN DIFFICULT SYSTEMS IN THE CONDITIONS OF UNCERTAINTY
The paper considers the problem of decision-making modeling for controlling of dynamic subsystems in complex organizational systems. The way of finding of transitions among etalon situations in the multilevel difficult systems, based on comparison of fuzzy intervals is offered. The system is described by the fuzzy situational network. There are considered two ways of factors description in etalon situations that are define transitions between etalon situations. Those ways are description by means of fuzzy numbers and by means of fuzzy intervals. Decisions for known problems of search of decisions in a fuzzy hierarchical control system are defined and resulted.
Difficult system; the operating decision; uncertainty; fuzzy numbers; fuzzy intervals; a reference situation; fuzzy inference.
Проблема разработки методов моделирования и исследования функционирования сложных систем в последнее десятилетие стала весьма актуальна. К классу таких систем, как правило, относят крупные энергетические или производственные комплексы с автоматизированным управлением, социально-экономическую систему, вычислительные комплексы поддержания работы глобальных сетей и т. п. Отличительной и неотъемлемой характеристикой сложной системы является наличие у нее эмерджентных свойств. Измерение таких свойств актуальный вопрос, но количественно их измерить весьма затруднительно, поэтому говорят об изучении косвенных проявлений системы.
К классу сложных систем относятся различные организационные системы. Особый интерес представляют организационные системы, состоящие из нескольких подсистем, которые могут действовать согласованно или конфликтно в целях достижения наилучшей (в некотором смысле) общей для системы динамики нарастания некоторого необходимого качества и выстроены иерархично. При моделировании выбора решения по управлению такими динамическими подсистемами приходится принимать во внимание интересы всех сторон, влияющих своими действиями на реальную динамику накопления качества, поскольку осуществляя их они оказывают активное или пассивное противодействие выполнению мероприятий, намеченных другой стороной. В таких подсистемах учет интересов сторон представляет слабо-
* Работа поддержана РФФИ, проект № 11-01-00011а.
структуризованную проблему, т.е. каким-либо образом можно описательно представить критерии, но установление зависимости между ними - принципиально невозможно, только на основе имеющейся у исследователя объективной информации.
Для исследования таких сложных систем предлагается создавать иерархию взаимосвязанных моделей - систему моделей [1]. Система моделей включает наиболее агрегированные (обобщенные) и детализированные модели, определяет условия их взаимодействия: выходные переменные агрегированных моделей рассматриваются как экзогенные, так и эндогенные параметры детализированных моделей. На каждом уровне системы определяется набор вершин, среди которых имеются и эталонные. Определение переходов между эталонными ситуациями зависит от способа описания признака в эталонной ситуации.
Рассмотрим два способа описания признаков в условиях неопределенности. Признаки в эталонных ситуациях описываются нечеткими числами. Набор эталонных ситуаций 5,^ не содержит нечетко равных ситуаций при заданном
пороге равенства 1;пс. в отличие от набора типовых ситуаций шага управления 5 1. Поскольку плохо определенных ситуаций нет, то в один класс эквивалентности попадают нечетко равные между собой ситуации (которые можно считать одной ситуацией) в пределах достоверности, ограничиваемой порогом нечеткого равенства 1;пс, м(5;, 5.) > ^пс. Для идентификации типовых ситуаций эталонными выполняются следующие операции [2]:
1) определяется степень равенства /и(Б1,5;)ситуаций 5 в с помощью
формулы М~,~,) = &(г<~,~.)&КМ(~-,5;)) = &(М (у),// (у))
уеУ уеУ
где М(М5;(у,),М5.(у,)) = &, В(Мм;(у„)(ТкР),Мм.(у„)(ТкР)),
)) =■
ММ(ур)(Ю (Уp)(TР'), если /Мр)^^1 ^«с^
и М1 (УР)(Тк ^(1—ос^'пД
1 если ,1пс) или Мм у.Ас);
^(Ур^к Чгс’Чгс; (ур)^к .
2) определяется ситуация , имеющая максимальную степень равенства с ^ т.е. М(50,~) = тахм(~0,~).
5 е550
С учетом дерева управляющих воздействий для иерархической системы определяются переходы из некоторой вершины иерархии, к которой типовая ситуация оказалась наиболее близкой, во все возможные эталонные ситуации следующего уровня иерархии.
Определение переходов между эталонными состояниями под воздействием управляющих воздействий разной силы, допустимых с учетом желаемого уровня достижения цели “11-12”, позволяет представить совокупность сценариев с интервальными оценками исходов для представителей разных классов в виде нечеткой ситуационной сети управления и подробно описаны в [3].
Признаки в эталонных ситуациях описываются нечеткими интервалами. В этом случае сравнение нечетких множеств производится не путем определения степени равенства, а с помощью следующей процедуры [4]. Опишем процедуру сравнения нечетких множеств и Ск на единичном интервале. Необходимо пройти следующие шаги.
1. Представить С0; и СО на единичном интервале, для чего нечеткие числа,
описываемые этими множествами, умножить на масштабирующий коэффициент 1/ (1).
2. Построить а-уровневые множества С0аи Ска на основе СО и СО .
3. Определить средние величины М (%) элементов множества Ска с помощью следующих формул:
♦ при построении термов с использованием статистических данных в случае
О < *11 < * 12 < *21 * *22 * ... * *„1 * *„2 < 1,
* + *
„ I(*^)(* г2 - * ,1)
М С ) = ^{* г1 < V < * г2 }, М С ) = г= „ 2-------------;
I (* г2 - * г1)
г =1
♦ при построении термов с помощью метода парных сравнений
м (СЬ =1 £С‘г.
Пы
4. Определить средние величины М С) элементов множества С0а аналогично п. 3.
Ршах
5. Определить функцию упорядоченности рС) = | МСг )йр, где Ршах -
О
максимальная степень принадлежности элементов нечеткого множества
6. Определить функцию упорядоченности Б( С'') аналогично п. 5.
7. Если Б( СО )<Р( Ск ), то множество СО больше, чем СО .
В сложной системе в условиях неопределенности возникает необходимость решения следующих проблем:
1) представителем какого класса по способностям к приобретению нужного качества является субъект;
2) наблюдается ли отклонения от эталонной промежуточной цели в рассматриваемый момент времени;
3) какое из максимизирующих управлений применимо для устранения отклонения от эталонной промежуточной цели;
4) моделирование силы управляющего воздействия, позволяющего осуществить переход из ситуации в ситуацию;
5) определение возможности события, описывающего признак уч в момент ^. Для поиска решений перечисленных проблем в нечеткой иерархической системе управления, где отдельные вершины представлены нечеткими интервалами с границами на разных шкалах, неприменим подход, основывающийся на определении степени нечеткого равенства нечетких чисел. Поэтому предлагается применять второй подход, основанный на сравнении нечетких интервалов. Приведем пример решения для каждой проблемы.
Проблема 1. Выяснить, по исходу взаимодействия является ли 12 представителем класса к с “высокими” способностями по приобретению нужного качества, описываемого эталонной ситуацией Ск . В текущий момент субъект находится в
ситуации С О .
Для представления двух сравниваемых нечетких интервалов, представленных на шкалах различной размерности, воспользуемся следующими функциями принадлежности термов лингвистической переменной «рейтинг»:
Со={0/<0,6>; 0,02/<0,7>; 0,16/<0,8>; 0,39/<0,9>; 1/<1>; 1/<1,33>; 1/<1,5>;
0,27/<1,66>} и Ск ={0 /<0,5>; 0,36/<0,6>; 0,39/<0,7>; 0,5/<0,8>; 1/<0,9>; 1/<1>; 1/<1,2>; 0,27/<1,33>}.
Для сравнения нечетких множеств С 0 и Ск предварительно представим их
на единичном интервале: (~0={0/<0,36>; 0,02/<0,42>; 0,16/<0,48>; 0,39/<0,54>; 1/<0,6>; 1/<0,8>; 1/<0,9>; 0,27/<1>} и Ск = {0 /<0,375>; 0,36/<0,45>; 0,39/<0,525>; 0,5/<0,6>; 1/<0,675>; 1/<0,75>; 1/<0,9>; 0,27/<1>}.
Построим на основе С 0 уровневые множества М0а:
С0а= {0,36; 0,42; 0,48; 0,54; 0,6; 0,8; 0,9; 1} для 0<а< 0,3;
С0а= {0,54; 0,6; 0,8; 0,9} для 0,3<а< 0,5;
С0а= {0,6; 0,8; 0,9} для 0,5<а< 1.
М0а= 5,1/ 8= 0,6375 для 0<а< 0,3;
М0а=2,84/ 4= 0,71 для 0,3<а< 0,5;
М0а=2,3/ 3= 0,7666 для 0,5<а< 1.
Определим функцию упорядоченности этого множества 0,3 0,5 1
Б(С0) = 10,6375^а+ 10,7Ша+ 10,7666йа=
0 0,3 0,5
= 0,3*0,6375+0,2*0,71+0,5*0,7666= 0,7165.
Построим на основе Ск уровневые множества Мка:
Ска= {0,375; 0,45; 0,525; 0,6; 0,675; 0,75; 0,9; 1} для 0<а< 0,3;
Ска= {0,45; 0,525; 0,6; 0,675; 0,75; 0,9} для 0,3<а< 0,5;
Ска= {0,675; 0,75; 0,9} для 0,5<а< 1.
Мка=5,275/ 8= 0,6594 для 0<а< 0,3;
Мка=3,9/ 6= 0,65 для 0,3<а< 0,5;
Мка= 2,325/ 3= 0,775 для 0,5<а < 1.
Функция упорядоченности
0,3 0,5 1
Б(Ск) = 10,6594^а+ 10,65^а+ 10,775йа =
0 0,3 0,5
= 0,3* 0,6594+0,2* 0,65+0,5* 0,775 = 0,7153.
Поскольку Б(С0)> Р(Ск), то 12 является представителем класса к с “высокими ” способностями по приобретению нужного качества.
Проблема 2. Выяснить, имеется ли отклонение от промежуточной цели по признаку уч (рейтинг) для 12 как представителя класса к с “высокими ” способностями по приобретению нужного качества. С 0 представим на единичном интервале: С0={0/<0,4>; 0,31/<0,5>; 1/<0,6>; 1/<0,75>; 0,82/<0,8>; 0,29/<1>}.
С0а= {0,4; 0,5; 0,6; 0,75; 0,8; 1} для 0<а< 0,3;
С0а= {0,5; 0,6; 0,75; 0,8} для 0,3<а< 0,5;
С0а= {0,6; 0,75; 0,8} для 0,5<а< 1.
М0а= 4,05/ 6= 0,675 для 0<а < 0,3;
М0а= 2,65/ 4= 0,6625 для 0,3<а < 0,5;
М0а= 2,15/ 3= 0,7166 для 0,5<а< 1.
0,3 0,5 1
Б(С0) = 10,675^а + 10,6625^а + 10,7166^а=
0 0,3 0,5
= 0,3*0,675+0,2*0,6625+0,5*0,7166=0,7153.
Поскольку эталонная ситуация Ск осталась неизменной и взята из проблемы 1,
для сравнения нечетких множеств С 0 и Ск ={0/<0,375>; 0,36/<0,45>; 0,39/<0,525>;
0,5/<0,6>; 1/<0,675>; 1/<0,75>; 1/<0,9>; 0,27/<1>} воспользуемся ранее найденной функцией упорядоченности Р(Ск)= 0,7153.
А учитывая, что Р(С0)=Б(Ск), то отклонения от эталонной промежуточной цели не наблюдается и для дальнейшего взаимодействия к 12 как представителю класса к с “высокими ” способностями по приобретению нужного качества следует применять управления со “слабой” силой воздействия.
Проблема 3. Результат взаимодействия 12 описывается нечетким множеством С0={0/<0,4>; 0,31/<0,5>; 0,67/<0,6>; 1/<0,7>; 1/<0,9>; 0,82/<1>; 0,29/<1,2>}. Это
одна из ситуаций множества Парето-оптимальных решений ю3=[0,81; 0,82].
Выяснить, имеется ли отклонение от промежуточной цели по признаку уч (рейтинг) для 12 как представителем класса к с “высокими ” способностями по приобретению нужного качества.
С0 представим в виде нечеткого множества на единой шкале:
С0={0/<0,3>; 0,31/<0,4>; 0,67/<0,5>; 1/<0,6>; 1/<0,75>; 0,82/<0,8>; 0,29/<1>}.
С0а= {0,3; 0,4; 0,5; 0,6; 0,75; 0,8; 1} для 0<а< 0,3;
С0а= {0,4; 0,5; 0,6; 0,75; 0,8} для 0,3<а< 0,5;
С0а= {0,5; 0,6; 0,75; 0,8} для 0,5<а< 1.
М0а= 4,35/ 7= 0,6214 для 0<а< 0,3;
М0а= 3,05/ 5= 0,61 для 0,3<а< 0,5;
М0а= 2,65/ 4= 0,6625 для 0,5<а< 1.
0,3 0,5 1
Б(С0) = 10,6214^а + 10,6Ша + 10,6625^а=0,3* 0,6214+0,2* 0,61 +0,5* 0,6625 0 0,3 0,5
= 0,1864+ 0,122+ 0,3312= 0,6396.
Поскольку эталонная ситуация Ск осталась неизменной и взята из проблемы 1, для сравнения нечетких множеств С0 и Ск ={0 /<0,375>; 0,36/<0,45>;
0,39/<0,525>; 0,5/<0,6>; 1/<0,675>; 1/<0,75>; 1/<0,9>; 0,27/<1>} воспользуемся ранее найденной функцией упорядоченности Б(Ск)= 0,7153.
Поскольку Р(С0)<Б(Ск), то имеет место отклонение от эталонной промежуточной цели и для дальнейшего взаимодействия с 12 требуется применить максимизирующее управление с тем, чтобы на следующем шаге попытаться вернуть субъект к эталонной ситуации для представителя класса к с “высокими ” способностями по приобретению нужного качества.
Проблема 4. Признак уч эталонной ситуации, описывающей состояние представителя класса с “высокими” способностями в момент ^.1, представим нечетким исходом.
% = сД‘={0/< Т^ 1 0,31/<Т1 >; 1/<Т/°'м) >; 0,82/<Т/10 >; 0,29/<
<Г/(гм) >}= {0/<0,7>; 0,31/<0,8>; 1/<0,93>; 0,82/<1>; 0,29/<1,2>}.
Этот признак для эталонной ситуации в ^ описывается Ск‘>= й2 ={0/<>;0,31/<Т^1"1 >;1/<Т3Ф1'} >;1/<Т4Ф,'} >;0,27/<Т5Ф1'} >}=
= {0/<0,7>; 0,31/<0,8>; 1/<0,93>; 1/<1,2>; 0,27/<1,32>}.
Выполнить моделирование силы управляющего воздействия, позволяющего
осуществить переход из Ск’1 в С к*1. Сила управляющего воздействия определяется с помощью формулы Ф*(1н)= Ск* 1 х Ск’. Опишем нечеткие множества с помощью термов. Для заданных ситуаций она будет описывается матрицей Фг(1>1):
Т 5 (Г 3 ) Г 1 т 5 (Г 3 ) Г 2 Т 5 (Г 3 ) Г 3 Т 5 (Г 3 ) Г 4 гр 5 ( Г з ) Г 5
Т5(Г 3 -1 ) Г1 0 0 0 0 0
Т5(Гз-1) Г 2 0 0,31 0,31 0,31 0,27
Т5(Гз-1) Г3 0 0,31 1 1 0,27
Т5(Гз-1) Г4 0 0,31 0,82 0,82 0,27
Т5 (Г 3-1) Г 5 0 0,29 0,29 0,29 0,27
Проблема 5. Задано нечеткое множество, описывающее признак эталонной ситуации в ^-1 следующим образом: С^’ ~‘ = £~1={0/<Т^1 -1)>; 0,31/<Т^1 -1) >; 1/<
<Т/10 >; 0,82/<Т4?(11) >; 0,29/<Т/(11) >} = {0/<0,7>; 0,31/<0,8>; 1/<0,93>;
0,82/<1>; 0,29/<1,2>} и сила управляющего воздействия Фг(1н) из проблемы 4.
Ск*1' ={0/<Т^1) >; 0,5/<Т2Ф 1) >; 1/<Т3Ф1) >; 1/<Т4Ф 1) >; 1/<Т/(*1 )>}= {0/<0,7>;
0,5/<0,8>; 1/<0,93>; 1/<1,2>; 1/<1,33>}
Определить возможность события, описывающего признак уч в момент ^. [5].
Возможность события С 1 к при наличии сведений С/11 ° Ф‘({ 11) обознача-
ется РОББ (Ск 1 / Ск 11 ° Ф 11 ) и может быть определена по формуле
РО88 (Ск ‘ /с;- о ф —) = sup тіп(ц .(ТТ”),Ис „ )),
т; з етм 5 з Ск оФ Ск
где Ц^(Т^) = sup тіп^,,- лт?1-0)^^,- Ж3,^)). Имеем:
Ц = (Гі5^)) = пшх (0; 0; 0; 0; 0) = 0;
Ск оФ 3 1,1 т5(,1)
Ц ,. () = Г5'*) = тах (0; 0,31; 0,31; 0,31; 0,29) = 0,31;
Ск оФ <Г,-1) 2 Г5(.1) 4 '
Цз ФГ(Г,1) = Т35(Г.) = трс (0; 0,31; 1; 0,82; 0,29) = 1;
с к о ф 31 ° г 5 (Г,-1)
ЦСг, фг(г,-1) = Т4д(Г) = пик (0; 0,31; 1; 0,82; 0,29) = 1;
с к о ф 3 1 ^ т 5 (Гз -1)
U (= T5q(tj) = max (0; 0,27; 0,27; 0,27; 0,27) = 0,27.
С^ о Ф (tj-1) 5 Tq(tj-i)
Тогда композиция представляет собой: Ckhl о = {0/< Tiq(tj) >; 0,31/<
T2q(tj) >; 1/<T3q((j) >; 1/<T4q(tj) >; 0,27/<T5q((j) >} = {0/<0,7>; 0,31/<0,8>; 1/<0,93>; <1/<1,2>; 0,27/<1,33>}. А возможность события определяется как
POSS (Cktj / Ckh 1 о 0l{(i-l)) = max (0; 0,31; 1; 1; 0,27) = 1.
TNq(tj)
В статье рассмотрен способ нахождения максимизирующего управления, который позволяет работать с вершинами нечеткой ситуационной сети и текущим состоянием управляемой системы, описанным в виде нечетких интервалов. Оптимальность найденного решения обеспечивается за счет того, что при выполнении прогноза учитывается полная предыстория процесса управления.
Предложенный подход, использующий нечеткое моделирование при анализе сложных систем и сложноформализуемых процессов, в зависимости от специфики организационных систем предоставляет эффективные методы и средства для проведения прогнозных исследований в иерархических многоуровневых системах с целью принятия оптимального решения.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Боженюк А.В., Гинис Л.А. Об использовании нечетких баз и антибаз при анализе нечетких когнитивных карт. - НАНУ, Институт проблем искусственного интеллекта, «Искусственный интеллект». - Украина, Донецк: ТПТТТТ «Наука i осв1та». - 2004. - № 4. - C. 276-285.
2. Мелихов А.Н., Берштейн Л.С., Коровин С.Я. Ситуационные советующие системы с нечеткой логикой. - М.: Наука, 1990. - 272 с.
3. Вовк С.П. Ситуационное управление и нечеткие игры в моделировании организационных систем. - Таганрог: Изд-во ТРТУ, 2002, - 96 с.
4. Борисов А.Н., Крумберг О.А., Федоров И.П. Принятие решений на основе нечеткого описания системы исходов // Принятие решений на основе нечетких моделей. - Рига: Зинатне, 1990. - С. 126-132
5. Vovk S.P., Ginis Ь.А. Modelling and forecasting of transitions between levels of hierarchies in Difficult formalized systems // European Researcher. - 2012. - Vol. (20), № 5-1. - P. 541-545.
Статью рекомендовал к опубликованию д.т.н., профессор В.П. Карелин.
Гинис Лариса Александровна - Федеральное государственное автономное образовательное учреждение высшего профессионального образования «Южный федеральный университет»; e-mail: gla@sfedu.ru; 347928, г. Таганрог, пер. Некрасовский, 44; тел.: 88634371743; кафедра информационно-аналитических систем безопасности; к.п.н.; доцент.
Вовк Светлана Павловна - e-mail: vovk61@list.ru; кафедра информационно-аналитических систем безопасности; к.т.н.; доцент.
Ginis Larisa Aleksandrovna - Federal State-Owned Autonomy Educational Establishment of Higher Vocational Education “Southern Federal University”; e-mail: gla@sfedu.ru; 44, Nekrasovskiy, Taganrog, 347928, Russia; phone: +78634371743; the department of information-analytical systems of safety; cand. of eng. sc.; associate professor.
Vovk Svetlana Pavlovna - e-mail: vovk61@list.ru; the department of information-analytical systems of safety; cand. of eng. sc.; associate professor.