Метод построения нечеткой полумарковской модели функционирования сложной системы Текст научной статьи по специальности «Математика»

ное потоковое программирование и его реализация на компьютерных системах // Изв. РАН. ТиСУ. 2005. № 1. С. 75-96.

9. Kutepov V.P., Malanin V.M., Pankov N.A An approach to the development of programming software for distributed computing and information processing systems. ICSOFT-O8, International conference on software and data technologies, Porto, Portugal, 2008, pр. 83-90.

10. Бажанов С.Е., Кутепов В.П., Шестаков Д.А Язык функционального параллельного программирования и его реализация на кластерных системах // Изв. РАН. Программирование. 2005. № 5.

11. Бажанов С.Е., Воронцов М.Н., Кутепов В.П. Структурный анализ и планирование процессов параллельного выполнения функциональных программ // Изв. РАН. ТиСУ. 2005. № 6. С. 111-126.

12. Milner R. A calculus for communicating systems. Lecture notes in computing science. Springer - Verlag, New York, 1980. Vol. 92.

13. Хоар Ч. Взаимодействующие последовательные процессы. М.: Из-во «Мир», 1989. С. 240.

14. Журнал для разработчиков MSDN Magazine. М.: 2008, 11(83). URL: www.microsoft.com/rus/msdn/magazine (дата обращения: 16.04.2010).

15. Апериодические автоматы; под ред. В.И. Варшавского. М.: Наука, 1976. С. 424.

16. Jones S.T. Peyton. The implementation of functional programming languages. Prentice Hall, 1987.

17. Milner R. The standard ML core language. Polymorphism // The ML/LCF/ Hope Newsletter, October, 1985. Vol. 2. № 2.

УДК 004.8

МЕТОД ПОСТРОЕНИЯ НЕЧЕТКОЙ ПОЛУМАРКОВСКОЙ МОДЕЛИ ФУНКЦИОНИРОВАНИЯ СЛОЖНОЙ СИСТЕМЫ

(Работа поддержана РФФИ, грант 10-01-97506-р_центр_а)

Ю.Г. Бояринов, к.т.н.; В.В. Борисов, д.т.н.; В.И. Мищенко, д.т.н.; М.И. Дли, д.т.н. (Смоленский филиал Московского энергетического института (технического университета),

BYG@yandex. ги)

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

Ключевые слова: нечеткая полумарковская модель, нечеткий вывод, нечеткая функция, нечеткое отображение.

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

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

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

• неопределенность переменных, обусловленную экспертным характером значительной части информации или эвристическим описанием процессов;

• разнокачественность данных, а также их оценку с помощью различных шкал;

• нечеткость выделения и описания переменных или отдельных состояний, а также входных и выходных воздействий.

Одним из основных способов решения подобных проблем при анализе функционирования систем является использование нечетких полумарковских моделей [1-3].

Следует отметить определенные ограничения существующих нечетких полумарковских моделей:

• не всегда четко обосновано соответствие вложенных цепей полумарковской модели классу цепей Маркова;

• имеющиеся модели в общем случае используются только для прогнозирования характеристик пребывания моделируемой системы в заданном состоянии, что затрудняет решение задач управления;

• предлагаемые подходы к построению нечетких полумарковских моделей не учитывают наличие избыточности ресурса как необходимое условие функционирования сложной системы.

Использование полумарковских моделей

для анализа функционирования систем

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

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

На рисунке S1 - состояние применения системы по назначению (желаемое состояние); S2 - состояние проверки работоспособности системы (пригодности к использованию по назначению); Sз - состояние проверки непригодности (отказа) к использованию по назначению, а также восстановление до состояния пригодности (готовности) к использованию по назначению после обнаружения дефекта (отказа); Ротк - вероятность возникновения отказа системы, Ротк=Fотк(ТП); Fотк(ТП)= =Вер(^,тк<ТП) - функция распределения случайной величины ^,тк момента отказа, соответствующая значению ^,тк<ТП, то есть вероятности отказа в течение периода проверки системы ТП.

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

Полумарковский процесс задается с помощью матрицы F(t) условных функций распределения продолжительности пребывания в состояниях, матрицы W переходных вероятностей вложенной марковской цепи и начального состояния процесса, из которого он стартует:

F(t) = S2

Я,

0

^ (Т) Fз1(Т)

&

(Т) 0 0

я.

Яз ^з (Т) 0 0

я,

(0 1 - ^ (Гп ) ^ (Тп )>

W = я,

я,

1 0 0

1 0 0

3 ^ У

Ненулевые компоненты матрицы F(t) запишутся в следующем виде:

(t) =

Чз (t) =

F21(!) =

Fзl(t) =

0, t < ТП,

1, t > Тп ; 0, t < 0,

F0IK (t)

0 < t < тп

(ТП ) 1, t > Тп ;

0, ! < ТК,

! > Тк ; t < ТК + Тв.

t > Тк + Тв.

Здесь ТК, ТВ - математические ожидания про-должительностей проведения проверки и восстановления после достоверного обнаружения отказа системы.

Безусловные функции распределения Fi(t) определятся в соответствии с выражением Fi(t)=

3

= ^ Б , где му - вероятность перехода вло-j=l

женной марковской цепи за один шаг.

Для рассматриваемого примера получим:

0, t < 0,

F0IK (1),0 < t < Тп,

М > Тп;

[0,1 < тК,

11, Т > Тк;

[0,1 < тк + тв,

[1,1 > Тк + Тв. Среднее время пребывания системы в соответствующих состояниях находится следующим образом:

Ц (О = F2 (1) =

Fз (1) =

(1)

ш, = /[1 - Fi (t)] dt.

(2)

Тогда т1 = /[1 - FоIк (t)]dt; т2=Тк; тз=Тк+Тв. о

Финальное распределение вероятностей состояний вложенной марковской цепи определяется в результате решения системы уравнений в матричном виде P=PW, где Р - вектор-строка (Р1, Р2, Рз).

Эта система уравнений является линейно зависимой, так как для нее выполняется условие нормировки. Таким образом, для рассматриваемого случая получим систему уравнений [Р, = Р2 + Р3 , Р = Р, (1 - Ротк ),

Р3 = Р,Ротк ,

Р, + Р2 + Р3 = 1. Решение этого уравнения имеет следующий

вид:

Р1 = - ; Р2 = 1 Р<,ТК 1 2 2

Рз =

2

Ротк ^отк (ТП )

1 - Ртк (Тп ) ; 2 ;

(3)

22

В стационарном режиме распределение вероятностей состояний для полумарковского процесса в целом определяется из выражения Р.т.

2 Р]т]

]=1

Очевидно, что показателем эффективности системы, функционирование которой описывается с помощью приведенного полумарковского процесса, будет вероятность застать систему в состоянии S1 в произвольный момент. И эта вероятность является функцией периода и продолжительности контроля, то есть ТП и ТК, а также функцией распределения моментов отказа системы Fотк.

Выражение этой вероятности имеет вид

/ [1 - (t)]dt

П =

1 Тп

/ [1 - ^ (t)]dt + Тк + ^ (Тп )ТВ о

Таким образом, числитель данного выражения соответствует средней продолжительности пребывания системы в исправном состоянии S1, а знаменатель - сумме средних времен пребывания в состоянии исправности, контроля и восстановления. В этой полумарковской модели в отличие от марковской учтены периодичность и продолжительность контроля [4].

Способы введения нечеткости в полумарковские модели

Анализ функционирования системы во многом субъективный процесс, в ходе которого требуется учитывать как количественные, так и качественные переменные, плохо поддающиеся формализации. Это затрудняет использование традиционных полумарковских моделей. Кроме того, как отмечалось, сложность использования вероятностного подхода анализа функционирования систем обусловлена уникальностью и сложностью типизации оцениваемых ситуаций и отсутствием достаточных статистических данных.

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

Предлагаются два способа введения нечеткости в полумарковские модели.

1. Дополнение выражения для оценки распределения вероятностей состояний функциями принадлежностей времени пребывания системы в соответствующих состояниях.

В этом случае выражение для оценки распределения вероятностей состояний примет следую-Ц(т1 )Рт 1

щий вид: п =-

Р]т]

]=1

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

Данный способ введения нечеткости в полумарковскую модель заключается в следующем:

• вероятности состояний и времен пребывания системы в соответствующих состояниях заменяются на нечеткие числа (нечеткие множества);

• для задания нечетких переменных, характеризующих вероятности состояний и времен пребывания системы в соответствующих состояниях, используются нечеткие отображения, реализуемые в соответствии с одним из подходов (нечеткие продукции, нечеткие отношения, нечеткие функции);

• операции суммирования, произведения и деления над вероятностями состояний и времен пребывания системы в соответствующих состояниях заменяются на расширенные (например, по принципу обобщения Л. Заде или на основе интервального метода) операции суммирования, произведения и деления нечетких чисел соответственно.

Примечания:

1) использование указанных выше способов введения нечеткости в полумарковские модели зависит от характера используемой информации и особенностей решаемых задач анализа функционирования систем;

2) в предлагаемом методе построения и использования нечеткой полумарковской модели функционирования системы реализован второй способ введения нечеткости в полумарковскую модель.

Метод построения нечеткой полумарковской модели функционирования системы

Рассмотрим основные этапы предлагаемого метода.

Этап 1. Задание нечетких переменных р на базовых множествах значений вероятностей Pi состояний системы.

При этом значения вероятностей Pi состояний системы являются базовыми множествами, на ко-

Т

о

торых определены функции принадлежности нечетких множеств, число которых для каждого состояния системы может быть различным. Для наглядности при оценке будем использовать по три нечетких множества (термы {Низкий, Средний, Высокий}) для задания каждой нечеткой переменной ]: для р - Нр1, Срь Врь для Р2 - Нр2, Срг,

Врг; для Р3 - НР3, СР3, ВР3. С одной стороны, это не нарушает общности при рассмотрении предлагаемой нечеткой полумарковской модели (то есть при необходимости перечень термов для каждой

нечеткой переменной ] может быть изменен), с другой стороны, позволяет типизировать компоненты модели, что существенно упрощает задачу при ее построении и использовании.

Отметим, что, хотя возможно совпадение названий термов для различных нечетких переменных, на самом деле они могут характеризоваться разными параметрами.

Для формирования терм-множеств нечетких

переменных ] (в виде нечетких множеств) и построения их логико-лингвистических шкал целесообразно использовать типовые Ь^-функции, например, колоколообразного типа [5]. Так, для

показанной на рисунке 2 нечеткой переменной Р2 значение НР2 задается следующим образом:

(P?) = exp

«in

Г х2

1 (P? - а2

2

b,

где a2, b2 - параметры функции принадлежности колоколообразного типа.

Этап 2. Задание нечетких переменных mt на базовых множествах значений вероятностей mi времен пребывания системы в соответствующих состояниях.

Все значения вероятностей времен mi пребывания системы в соответствующих состояниях являются базовыми множествами, на которых определены функции принадлежности нечетких множеств нечетких переменных m j.

Как и для нечетких переменных Pj, задаются нечеткие множества для переменных inj: для mt

- Hmi, Cmi, Bmi; для m2 - Н2, Omz, В^г; для 1П3 -Нщз, Сщз, B„3 (аналогично этапу 1).

Этап 3. Задание нечетких отображений для нечетких переменных mj, характеризующих вероятности времен пребывания системы в соответствующих состояниях.

Исходя из выражений (1) и (2) зададим следующие нечеткие отображения для нечетких переменных inj: riij =cpi (ТП ),m2 = ср^ (ТК ),m3 = = ф3 (ТК, ТВ), где pp. - нечеткое отображение.

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

Рассмотрим пример реализации нечеткого отображения Рр3 для нечеткой переменной m3 на основе нечеткой продукционной модели MISO-ти-па (Multiple Inputs, Single Outputs) в предположении, что входные переменные ТК, ТВ также являются нечеткими и задаются соответствующими

нечеткими множествами: для ТК - Нк, СК, ВК; для ТВ - Нв, Св, Вв:

П1: ЕСЛИ -ГК есть Нк И ТВ есть Нв, ТО im3 есть Нтз, П2: ЕСЛИ ТК есть Нк И ТВ есть Св, ТО im3 есть Нтз,

П;: ЕСЛИ -ГК есть Ск И ТВ есть Св, ТО im3 есть Стз.

П9: ЕСЛИ -ГК есть Вк И ТВ есть Вв, ТО im3 есть Втз.

Для нечеткой переменной im3 в рамках сформированной нечеткой продукционной модели реализуется алгоритм нечеткого вывода Мамда-ни [5].

Этап 4. Задание нечетких отображений для нечетких переменных Pj, характеризующих вероятности состояний системы.

В соответствии с выражением (3) зададим следующие нечеткие отображения для нечетких

переменных Pj: Pj = ср (const), P2 = ф2 (ТП), P3 =cp3 (ТП), где cpi - нечеткое отображение.

Если входную переменную ТП представлять как нечеткую переменную ТП, то нечеткое отображение cp можно реализовать на основе нечеткой продукционной модели SISO-типа (Single Inputs, Single Outputs).

Рассмотрим пример реализации нечеткого отображения ф для нечеткой переменной P2 на основе нечеткой продукционной модели SISO-ти-па с учетом того, что нечеткая переменная ТП задается нечеткими множествами Нп, Сп, Вп:

П1: ЕСЛИ -ГП есть НП, ТО P2 есть BP2.

П2: ЕСЛИ ТП есть СП, ТО Р2 есть СР2.

П3: ЕСЛИ ТП есть ВП, ТО Р2 есть Н«.

Далее для нечеткой переменной Р2 в рамках сформированной нечеткой продукционной модели реализуется алгоритм нечеткого вывода Мамдани.

Если переменная ТП является четкой, то для реализации нечеткого отображения ф целесообразно использовать нечеткую функцию от этой четкой переменной ТП, формирующей образ ее четкой области определения в соответствующем

нечетком множестве : ТП ^ Р..

Следует отметить, что, поскольку понятия нечеткой функции четкой переменной и нечеткого отношения соответствуют друг другу в математическом смысле, нечеткая функция ф может быть

интерпретирована как нечеткое отношение Й и определена следующим образом:

У(1п , р!) е Тп х Р1, (Р,) = ^ (tП , р1).

Этап 5. Определение значений нечетких переменных, характеризующих распределение вероятностей состояний для нечеткого полумарковского процесса.

Для выполнения данного этапа предлагаемого метода введем необходимые понятия, определения и выражения для расчетов.

Нечеткое число X определяется как нечеткое множество, заданное на множестве действительных чисел R с функцией принадлежности ц XX (х) е [0,1], где xeR, и для которого выполняются следующие условия:

а) 3 xeR, такое, что ц ^ (х) = 1;

б) множества а-уровней нечеткого числа X представляют собой закрытые ограниченные интервалы Ха=[а, Ь]={х/ (х) >а};

в) 3 NeR, такое, что V xeR, если |х|>^ то Цх (х) =1.

Операции над нечеткими числами могут выполняться с использованием различных методов.

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

Пусть заданы два нечетких числа X и У, а также множества их а-уровней Xа=[а, Ь]={х/ (х) >а} и Уа=[с, d]={x/ (х) >а для каждого

значения ае[0, 1].

Интервальные операции над множествами

а-уровней нечетких чисел Ха=[а, Ь] и Уа=[с, d] вычисляются в соответствии со следующими выражениями [6].

Сложение:

XX а + X* = [а,Ь] + [c,d] = [а + с, Ь + d].

Умножение:

XXа хУа = [а, Ь]х[с, d] = [1шп{ас, Ьс, ad, bd},

тах{ас, Ьс, ad, bd}].

Деление:

XXа/V* = XXа X 1/= [а, Ь] X [1 / ё, 1 / С] =

= [тт{а/ё, Ь/ё, а/с, Ь / с},

тах{а/ё, Ь /ё, а/с, Ь / с}], если 0 ё [с, ё],

где =[1/а,1/с], о г уа.

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

По другому методу выполнения вычислений, основанному на принципе расширения Л. Заде, произвольная расширенная бинарная арифметическая операция * над нечеткими числами может быть задана следующим образом:

Ъ = XX * У (г) = V (^х (X) лц? (У)),

ъ=х х ? '

где XX, У,22 - нечеткие числа, * - обозначение произвольной расширенной бинарной арифметической операции над нечеткими числами [5].

Тогда, например, расширенная операция произведения ® над нечеткими числами Р(А = а4) и Р(В = Ь | А = а1) имеет следующий вид:

Р(А = а11В = Ь) =

=(Р(А=а) ® Р (В=Ь | А=а))«

^^(л^в^)(С) =

= с=(А=^) ^(А) =,) (В))

Примечание: в дальнейшем при осуществлении операций над нечеткими числами будем использовать интервальный арифметический метод над множествами а-уровней этих чисел.

Для задания правил вывода на основе нечеткой байесовой сети рассмотрим понятия нечеткой вероятностной меры и нечеткого вероятностного пространства.

Функция р : £^F(R) называется нечеткой

вероятностной мерой [7] на (□, в) в том и только в том случае, если:

а) 0Х ^ р (А) ^ 1х для VAe8;

б) Pf (□) =1Х и Pf (0) =0Х;

в) если А и В являются несовместными значениями из в (то есть А/пВ=0), то р (А) © рр (В) э

Э Р(А иБ) ;

г) если A и B являются событиями из 8, то 0x, если Pf (A) *Pf (B) < 0x, Pf (A) * Pf (B),

если 0x < Pf (A) * Pf (B) < 1, 1, если Pf (A)*Pf (B) > 1x.

Pf (A)*I>f (B) =

Здесь Q=Rm - некоторое универсальное множество, на котором задана переменная А; 8 - множество несовместных значений переменной А; F(R) - множество всех нечетких чисел, определенных на множестве действительных чисел R; 0Х и 1х - нечеткие числа 0 и 1.

Нечетким вероятностным пространством

называется тройка (□, 8, р).

Таким образом, в основе предлагаемого метода построения и использования нечеткой полумарковской модели лежит следующее выражение:

аi = (р ёмi) 7 ё мj) , (4)

где п 1 - нечеткая переменная, характеризующая распределение вероятностей состояний для нечеткого полумарковского процесса.

В заключение отметим, что в данной работе делается анализ функционирования систем на основе нечетких полумарковских моделей.

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

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

Предложен метод построения и использования нечеткой полумарковской модели функционирования системы.

Литература

1. Bhattacharyya M. Fuzzy Markovian decision process // Fuzzy Sets and Systems. Vol. 99, 1998, pp. 273-282.

2. Praba B., Sujatha R., Srikrishna S. Fuzzy reliability measures of fuzzy probabilistic semi-Markov model // Int. Journal of Recent Trend in Engineering. Vol. 2, No. 2, 2009, pp. 25-29.

3. Praba B., Sujatha R., Srikrishna S. A study on homogeneous fuzzy semi-Markov model // Applied Mathematical Sciences. Vol. 3, No 50, 2009, pp. 2453-2467.

4. Бояринов Ю.Г., Дли М.И., Круглов В.В. Оценка диапазона возможных значений вероятности пребывания в заданном состоянии марковской модели производственно-экономической системы // Программные продукты и системы. 2009. № 4. С. 88-91.

5. Круглов В.В., Борисов В.В. Искусственные нейронные сети. Теория и практика. М.: Горячая линия - ТЕЛЕКОМ, 2002.

6. Moore R. E. Interval analysis. New Jersey: Prentice Hall,

1966.

7. Kwakernaak H. Fuzzy random variables: definitions and theorems // Information Sciences. № 15(1), 1978, pp. 1-29.

УДК 004.415.2

ПРОГРАММНЫЙ ИНСТРУМЕНТАРИЙ ПРОЕКТИРОВАНИЯ БАЗ ЗНАНИЙ ДЛЯ ИНТЕЛЛЕКТУАЛЬНОГО УПРАВЛЕНИЯ

С.Г. Караткевич, кэ.н; Л.В. Литвинцева, кф.-м.н. (Научно-технологический парк «Дубна», г. Дубна, [email protected], s.ulyanov@ ntpdubna.ru);

О.Ю. Тятюшкина; П.Н. Григорьев

(Международный университет природы, общества и человека «Дубна», г. Дубна, [email protected], [email protected])

В статье дается анализ проблем реализации интеллектуальных систем управления. Описаны основные этапы работы с оптимизатором баз знаний. Эффективность процессов управления с оптимизатором баз знаний демонстрируется на конкретном примере модели динамического объекта управления («шарик-доска») в условиях неполной информации о параметрах структуры объекта управления.

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

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

Наряду с этим важными проблемами внедрения ИСУ являются также их слабая адаптивность

к изменениям параметров ОУ (вызванным, например, старением структуры ОУ или резким изменением внешней среды) и низкая робастность полученных законов управления в случае непредвиденных ситуаций управления. Актуальность решения данных проблем многократно диктовалась практическими задачами теории и систем управления и не раз отмечалась многими исследователями [1-3].

Для решения проблем были применены технологии интеллектуальных вычислений, вклю-