Анализ и решение задач выбора с параметрической нечеткостью Текст научной статьи по специальности «Математика»

УДК 510.22(075.8)

DOI: 10.14529/ mmp150402

АНАЛИЗ И РЕШЕНИЕ ЗАДАЧ ВЫБОРА С ПАРАМЕТРИЧЕСКОЙ НЕЧЕТКОСТЬЮ

М.Г. Матвеев

Для задач выбора, представленных моделями с параметрами в виде нечетких LR-чисел предложена методика решения, основанная на применении а-уровневого представления нечетких чисел, их дальнейшей модификации с помощью выпуклого линейного преобразования границ а-интервалов, сохраняющего основные характеристики нечеткости, предложенной алгебры модифицированных нечетких чисел и выпуклой линейной комбинации решений на границах промежутка изменения а. Достоинствами предложенной методики являются: ограниченность роста неопределенности при обработке нечеткой информации; сохранение естественной интерпретации промежуточных и конечных результатов вычислений; возможность организации вычислений

а

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

Ключевые слова: модели с параметрической нечеткостью; нечеткие числа LR-а

решения.

Введение и постановка задачи

Важными направлениями развития приложений теории нечетких множеств в управлении является нечеткая логика и нечеткая арифметика. Модели нечеткой логики можно рассматривать как альтернативу классическим аналитическим моделям теории управления, в частности теории принятия решений. При этом основные компоненты нечеткой логики: системы нечетких продукционных правил, нечеткие реляционные уравнения обладают рядом важных достоинств - относительной простотой синтеза системы управления в условиях нечеткой неопределенности, естественной ро-баСТНОСТЬЮ, возможностью учета человеческого опыта и интуиции. С другой стороны, модели и методы классической теории управления обладают несоизмеримо большими возможностями отображения разнообразных свойств систем управления (управляемость, устойчивость, оптимальность и т.п.) и обеспечивают направленность синтеза управления на их достижение. Нечеткая неопределенность может быть, при необходимости, включена в классические модели управления или принятия решений в виде параметров и переменных, оценки которых представляют собой приблизительные значения, получаемые обычно экспертным путем. Наиболее адекватным представлением таких оценок являются нечеткие числа, т.е. нечеткие подмножества числовой оси R с функцией принадлежности ß : R ^ [0, 1]. Известны несколько разновидностей нечетких чисел, а также соответствующие алгебры и алгоритмы, составляющие основу нечеткой арифметики и позволяющие производить необходимые вычисления

в условиях нечеткости параметров и получать решение задач в нечетком виде, например, [1-4]. Сочетание классических моделей и методов учета нечеткой неопределенности в данных представляется перспективным направлением развития теории управления. Разнообразные модели принятия решений часто используют представление неопределенности параметров в форме нечетких чисел, о чем свидетельствует большое количество публикаций по этой тематике, например, [4-8]. Широкое применение классических моделей управления и принятия решений с нечеткими параметрами сдерживается рядом существенных особенностей известных методов обработки нечетких чисел:

- операции над нечеткими числами могут повлечь неоправданное расширение неопределенности результата, делая его зачастую практически неприемлемым;

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

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

В работах [9, 10] были п р е дл оже н ы модифицированные треугольные нечеткие числа LR-типа и основные операции над ними, позволяющие преодолеть некоторые недостатки. Но правила выполнения операций над нечеткими числами в предложенной алгебре не совпадают с правилами выполнения аналогичных операций над действительными числами, что делает невозможным применение стандартных инструментальных программных средств решения задач управления, которые оперируют действительными числами. Для решения задач управления с нечеткими параметрами требуется либо разработка специального программного обеспечения, либо формирование такой алгебраической системы нечетких чисел, которая допускает использование стандартных пакетов (Maple, Matlab, Excel и т.п.). Очевидно, что второй путь менее затратный.

Кроме того, использование в работах [9, 10] а-уровневого подхода к решению широкого круга задач не гарантирует сохранение исходной формы представления нечетких чисел в полученном решении, что затрудняет его интерпретацию, а иногда делает ее невозможной. Это связано с проблемой устойчивости решений на различных а-уровнях, что обуславливает необходимость анализа устойчивости решения.

Цель настоящего исследования - подобрать алгебру изоморфную алгебре, предло-женнои в [9, 10], обеспечивающую возможность применения стандартных программных средств решения задач, а также разработать практические критерии для анализа устойчивости решения для некоторых типов задач.

1. Основные определения

Нечеткое треугольное чиело A LR-типа задается в виде тройки действительных чисел (xL,m,xR}, где т - мода, a xL, xR - соответственно левая и правая границы

носителя X. Функция принадлежности числа A определяется следующим образом

ßÄ{x) = <

x _ xxL

ßL =-r; x E lxL; m\,

л m — xL R

= ——; x E [m; xR], (1)

A m — xR

0 в остальных случаях.

Нечеткие числа можно представлять в виде объединения четких подмножеств Аа множества X, каждое го которых называется а-интервалом - Аа = {х Е X|^д(х) > а} а Е [0; 1]:

А= У аАа. (2)

ае [0;1]

Правая часть равенства (2) прбдст&влябт а

а

число п значений а, то равенство (2) позволяет свести решение исходной нечеткой задачи к совокупности п задач на соответствующих а-уровнях с параметрами, преда

а

практических задач, в качестве такой точки принимается середина интервала [7]. Границы а функциями обратными (1):

xL(a) = am + (1 — a)xL, xR(a) = am + (1 — a)x

, (1 (3)

В работе [9] была предложена алгебра P = (K; S\) на множестве K = {xr(a)} ; a E [0; 1] нечетких LL/RR-чжел. Элементами множества K

x(a) = a + ba (4)

ЯВЛЯЮТСЯ модифицированные нечеткие числа с треугольной функцией принадлежности, представленной только одной стороной треугольника, отображающей степень

b

кости равный длине носителя LL/RR - числа взятой с положительным (LL- число) или отрицательным (ДД-число) знаком; a E [0; 1] - уровневые значения функции принадлежности.

LL/RR-числа получаются из исходной экспертной оценки числовой величины LR-

(xL, m, xR)

a

x (a) = L (Xa) = \xL (a) + (1 — A) xR (a) , A E [0, 1]. (5)

Здесь A - управляющий параметр, значение которого определяется требуемыми свойствами решения задачи.

Операции Si над нечеткими числами (4) выполняются по следующим правилам

[9]:

xi(a) ф x2(a) = ai + a2 + (bi + b2)a, (6)

xi(a) ® x2(a) = aia2 + (aib2 + a2bi + bib2)a, (7)

х\(а) ф (—х2(а)) = а1 — а2 + (61 — Ь2)а, —х = —а — Ьа, (8)

Х1(а) ® х—1 (а) = а1/а2 + [(а2^ — аЬ)/^ + а2Ь2)]а,

Х-1(а) = а-1----г-а, а = 0, а = —Ь. (9)

^ а(а + Ь)

Так как в (6) - (9) операции линеины по переменной а, множество ^Сявляется замкнутым по всем операциям, то есть если хгЕ К, то хг * х^ Е К.

В работе [9] показано, что множество К с операциями (6) - (9) представляет собой алгебру модифицированных нечетких ЬЬ/ЯЯ -чисел с единственным нулевым элементом 0 = (0 + 0а) и единственным единичным числом 1 = 1 + 0а.

Алгебра Р = (К; Б1) позволяет применять операции независимо от расположения нечетких чисел на числовой оси, не допускает необоснованного расширения нечеткости, сохраняет треугольную форму нечеткого результата применения операций и естественные математические отношения в задаче, например, тождественность уравнения с нечеткими переменными понимается как равенство функций принадлежности левой и правой части уравнения.

Если по каким-либо причинам модификация нечетких чисел в ЬЬ/ЯЯ-числа недопустима, то в работе [10] предлагается алгебру Р = (К; Б1) применять отдельно к левой и правой части нечеткого ЬЯ-числа с треугольной функцией принадлежности. При этом все свойства результатов операций, достигнутые в [9], сохраняются.

Реализация операций (6) - (9) отличается от реализации аналогичных операций над действительными числами, что не позволяет применять при решении задач стандартные программные пакеты.

2. Построение подходящего изоморфизма и методики решения задач с нечеткими параметрами

Изоморфизмом алгебры Р = (К; Б1) на алгебру С = (М; Б2) будем называть отображение Г : К ^ М, удовлетворяющее условиям

ВД(к)) = ^(Г(к)), к Е К. (10)

^1(Г"1(т)) = Г"1^2(т)), т Е М. (11)

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

Пусть элементы множества М представлены в виде х(а) = х(0) + (х(1) — х(0))а, где х(0), х(1) - действительные числа. Сигнатура Б2 определяется выражением

х1(а) * х2(а) = ах(1) + (1 — а)х(0) = х(0) + (х(1) — х(0))а, (12)

где х(0) = х1(0) * х2(0); х(1) = х1(1) * х2(1); * Е {+, —, /} и все операции выполняются по правилам работы с действительными числами.

Предложение 1. Алгебры Р и С изоморфны.

МК

Очевидно, что любой элемент множества К -х(а) = а + Ьа получается из элементов

множества M простой подстановкой a = X(0), b = (X(1) — X(0)), в этом случае K = M, а отображение Г является отображением на себя. Множество K является замкнутым по операциям (6-9), поэтому для доказательства гомоморфизма (10) и обратного гомоморфизма (11) достаточно показать справедливость равенства S1(k) = S2(m). Поскольку результат операций в обеих алгебрах представляет собой линейную функцию от параметра а, достаточно показать, что эти результаты совпадают в двух точках, например, при а = 0 и а =1. Совпадение результатов легко проверяется выполнением операций по правилам (6-9) и (12) для двух нечетких чисел ai + bia и a2 + b2a.

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

Шаг 1. Преобразовать нечеткие параметры задачи, заданные как треугольные нечеткие числа в модифицированные нечеткие числа (4) с помощью преобразования

а

. . m — xL щее значение А для каждого нечеткого параметра, например А = —-- обеспечи-

XX^i _ XßL

вает минимальную потерю информации при переходе от LR-чисел к модифицированным LL/RR-числам [9].

Шаг 2. Решить задачу с модифицированными нечеткими параметрами как задачу с четкими параметрами при а=1. Обозначить решение - X(1).

а

решение - X(0).

Шаг 4. Построить модифицированное нечеткое решение как линейную комбинацию решений X(1) и X(0):

X(a) = X(0) + (X(1) — X(0))a.

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

Пример 1. Пусть эффективность некоторого производственного проекта описы-Ведется непрерывной функциональной зависимостью от изменения (увеличения или сокращения) объема выпуска товарной продукции - х, и эта зависимость имеет экстремум. Один из параметров этой зависимости определяется экспертным путем в виде нечеткого числа. Допустим, что необходимое условие экстремума приводит к решению квадратного уравнения X2 + Вx + 2 = 0, где В = (2; 3; 5).

Шаг 1. bL = 2 + a; bR = 5 — 2a. Пусть А выбирается из условия минимальных потерь экспертной информации, т.е. А=1/3. Тогда b(a) = 4 — а. Уравнение с модифицированным нечетким числом, которое предстоит решать, имеет вид

X2^) + (4 — а^(а) + 2 = 0.

Шаг 2. При а=1 решаем уравнение X2^) + 3X^) + 2 = 0. Получим X1 (1) = — 1, X2(1) = —2.

Шаг 3. При а=0 решаем уравнение X2^) + 4X^) + 2 = 0. Получим с точностью до двух знаков X1(1) = —0,59, X2(0) = —3, 41.

Шаг 4. Модифицированное нечеткое решение:

хг = -0, 59 - 0,41а; Х2 = -3, 41 + 1,41а. (13)

Нетрудно проверить, что решение (13) при подстановке в уравнение Х2(а) + (4 — а)х(а) + 2 = 0 превращает его в четкое тождество (в смысле равенства функций принадлежности).

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

Пример 2. Задача сетевого планирования проекта при нечетко заданных продол-ЖИТбЛЬНОСТЯХ работ. Пусть дан проект, заданный табл. 1. Соответствующий сетевой граф показан на рис. 1.

Таблица 1

Операция Предшествующие операции Нечеткая продолжительность

опе зации, г

г ь т г д

А 1 2 5

В 2 4 5

С А 3 7 9

Б В 4 6 9

Е В 9 10 12

Е С 4 5 6

С Б, Е 1 5 6

н Е, С 2 4 7

Требуется определить критические операции проекта и общее нечеткое время выполнения проекта. Примем для определенности единицу времени дни.

Необходимое решение можно получить как решение задачи линейного программирования с нечеткими параметрами:

tn - to ^ min;

tjs - tis > ts, Ув = 1,..., 8; I , (14)

T = E Tsi; i

si )

где t0, tn - моменты времени начала и окончания проекта; tis, tjs - моменты начала и окончания в-й операции; Ts - нечеткая продолжительность в-ой итерации; fsi -нечеткая продолжительность операций критического пути, т.е. операций для которых в результате решения выполняется равенство ¿jsi ¿is i — Ts i (в смысле равенства функций принадлежности левой и правой части уравнения). Такие ограничения в задачах линейного программирования принято называть дефицитными ограничениями; Т - продолжительность проекта.

Система неравенств в задаче (14) применительно к табл. 1 примет вид:

ti - ío > Ta (а); ¿2 - ¿o > tB (а); ¿5 - ¿i > te (а); ¿э - ¿2 > td (а); ¿4 - ¿2 > Te (а); ¿6 - ¿5 > tF (а); ¿6 - ¿4 > Tg (а); t7 - t6 > Тн (а).

Следует обратить внимание на то, что предложенная методика решения позволяет устранить проблему сравнения нечетких чисел в представленной системе неравенств, заменяя ее соответствующей системой сравнения действительных чисел. При решении следует предусмотреть фиктивную операцию, которой соответствует равенство Ьэ - t4 = 0.

Выполним необходимые шаги предложенной методики.

Шаг 1. Примем для всех нечетких чисел, определяющих продолжительность операций, параметр А=0,5. Соответствующие модифицированные числа примут вид:

Ta (а) = 3 - а; тв (а) = 3, 5 + 0,5а; fe (а) = 6 + а; т D (а) = 6, 5 - 0, 5а; тЕ(а) = 10, 5 - 0, 5а; тF(а) = 5; Tg (а) = 3, 5 + 1, 5а; тн (а) = 4, 5 - 0,5а.

Задача (13) принимает модифицированный вид:

тп(а) - ¿о (а) ^ min;

tjs^) - ^(а) > Ts (а), Ув = 1,..., 8; I (15)

т(а) = Е тsl(а) . I

si

Шаг 2. Решаем детерминированную задачу линейного программирования (15) а

следующий вид:

to(1) = 0; Ti(1) = 7; ¿2(1) = 4; Тэ(1) = 10; ¿4(1) = 14; T>(1) = 14; Тз(1) = 19; ¿7(1) = 23.

Критические операции: В; Е; G; Н. Продолжительность проекта: т(1) = 23 дня. а

¿o(0) = 0; ¿i(0) = 7, 5; ¿2(0) = 3, 5; ¿э(0) = 10; ¿4(0) = 14; ¿5(0) = 12, 5; tб(0) = 17, 5; t7(0) = 22.

Критические операции: В; Е; О; Н. Продолжительность проекта: т(0) = 22 дня.

Шаг 4. Модифицированное нечеткое решение задачи (15) будет ВКЛЮЧсХТЬ нечеткие интервалы моментов времени начала и окончания операций проекта:

°о(а) = 0; °Да) = 7, 5 — 0,5а; °Да) = 3, 5 + 0, 5а; ¿Да) = 10; ¿Да) = 14; ¿б(а) = 12, 5 + 1, 5а; ¿б(а) = 17, 5 + 1, 5а; ¿Да) = 22 + а.

Критический путь проекта составляют операции В; Е; О; Н. Результаты решения задачи (15) имеют естественную интерпретацию:

- значение критерия минимально по переменным ¿г (а) при всех значениях а; при

а

или влево рассматривается как мера риска;

- сумма нечетких продолжительностей операций критического пути равна (в смысле равенства функций принадлежности) нечеткому моменту времени окончания проекта, то есть нечеткая продолжительность проекта определяется выражением т(а) = 22 + а, т.е. возможно проект будет выполняться в интервале (22; 23) дня;

- общий резерв времени, например, операции А определяется выражением ¿Да) — Тд (а) = 4, 5 + 0, 5а, которое характеризует риск сокращения резерва времени в интервале (4,5; 5) дней. У критических операций резерва времени нет, например, общий резерв времени операции В равен ¿Да) — тв(а) = 3, 5 + 0, 5а — 3, 5 — 0, 5а = 0.

Выполнение этих естественных соотношении подтвер^кд&ет достоверность полученного решения.

Пример 3. В предыдущих примерах рассматривались нечеткие аналоги задач выбора, решением которых являются точки в п-мерном пространстве. Для того, чтобы проиллюстрировать нахождение решения в виде функции, будем рассматривать дифференциальное уравнение с нечеткими параметрами. Пусть имеется задача Коши с нечетким параметром

ш = А* л = (—22;0;1) ; (16)

у(г = 0) = 1.

Найдем частное решение этого уравнения, следуя предложенной методике. Шаг 1. Примем параметр А=0,5. Модифицированное представление нечеткого параметра уравнения (16) примет вид: 0(а) = —0, 5 — 0, 5а. Модифицированное представление дифференциального уравнения:

= 0, 5 — 0,5а)С(а); (17)

у(а, г = 0) = 1.

а

¿у(1)

¿г =—у(1);

0(1, г = 0) = 1.

Решение: у(1) = ехр(—г).

а

= — 0,50(0);

0(0, г = 0) = 1.

Решение: 0(0) = ехр(—0, Ы).

Шаг 4. Модифицированное нечеткое частное решение уравнения (17) п р еде т£ьв— ляет собой линейную комбинацию решений шагов 2 и 3:

y (а) = exp(—0, 5t) + (exp(—t) — exp(—0, 5Ь))а.

(18)

Подстановка решения (18) в уравнение (17) дает четкое тождество. Решение (18) представляет собой семейство экспонент (рис. 2). непрерывно зависящее от параметра а Е [0; 1]. Решение (18) будет устойчивым по Ляпунову для всех значений а.

Рис. 2. Устойчивое решение задачи Копти с нечеткими параметрами

В примере 3 решение имеет вид выпуклой линейной комбинации решений при а Е {0, 1} также как и в примерах 1 и 2, но является не нечетким числом, а нечеткой функцией двух переменных. В общем случае решение как модифицированная нечеткая функция будет иметь вид: y(t, а) = a(t) + b(t)a, а при фиксированном значении t

b(t)

t

t

3. Исследование проблемы устойчивости решения

K

нечетких чисел недостаточно для построения интерпретируемого решения задач с нечеткими параметрами. Если процедуру построения решения некоторой произвольной задачи обозначить как оператор X^) = Ф^; а(а)), где <1(а) Е K, то возникает вопрос - будет ли решение X^) Е K? То есть, будет ли нечеткое решение, полученное по предложенной методике, представимо как модифицированное нечеткое число? Подобные вопросы впервые рассматривались в работах [11. 12]. применительно к задаче линейного программирования с приближенно заданными параметрами. Небольшие изменения решения как результат небольших изменений параметров в этих работах

определялось как устойчивость решения. По аналогии свойство сохранения решением вида (4) будем называть устойчивостью решения.

Обсуждение проблемы устойчивости удобно начать с примеров.

Пример 1*. Вернемся к примеру 1 и введем новую нечеткую оценку параметра В = (1; 2; 5). Тогда для квадратного уравнения х2 + Вх + 2 = 0, а точнее его модифицированного аналога, полученного при А=0,5 и имеющего вид Х2(а) + (3 — а)Х(а) + 2 = 0 при а =1, решения не существует на множестве действительных чисел. При а=0 решение существует: Х1(0) = —0, 666 и Х2(0) = —3. Нетрудно проверить, что в точке а ~ 0,15 функция нечеткого решения имеет разрыв. По одну сторону от этой точки имеется непрерывное решение, по другую - решение отсутствует. Очевидно, что в этом случае предложенный метод получения решения, основанный на линейной

а

Пример 2*. Похожая ситуация может возникнуть и в примере 2. Поменяем параметры у операций Е и^в табл. 1 так, как это показано в табл. 2:

Таблица 2

Операция Предшествую щие операции Нечеткая продолжительность операций, Г

т ь т т д

А _ 1 2 5

В _ 2 4 5

С А 3 7 9

в В 4 6 9

Е В 6 7 9

Е С 4 5 12

С Б, Е 1 5 6

н Е, С 2 4 7

А

теперь нечеткое решение в виде ¿¿(а), г = 0; 1;...; 7, а также зависимость продола

этом легко убедиться, рассмотрев, например, оптимальное значение критерия задачи (14) в точках а Е [0; 0, 25; 0, 5; 0, 75; 1], представленное в табл. 3.

Таблица 3

Значение минимизируемого Критический путь

критерия

а=0 21,5 АСЕН

а=0,25 20,625 АСЕН

а=0,5 19,75 АСЕН

а=0,75 19,75 ВЕСН

а=1 20,0 ВЕСН

При а Е [0; 0, 5] функция изменяется линейно. На промежутке значений а от 0,5 до 0,75 функция решения имеет разрыв, так как критический путь изменяется на кон-

цах этого промежутка. Этот разрыв уже не позволяет интерпретировать результаты решения как в примере 2.

Примеры 1* и 2* показывают, как решение на промежутке изменения значений а теряют треугольную форму, что рассматривается как неустойчивость решений соот-ветствуюгцих задач. Можно утверждать, что решение большинства задач с нечеткими параметрами с помощью предложенной методики имеет интерпретируемый смысл только при условии устойчивости этого решения. В этой связи для практического использования предложенной методики актуально наличие простых, легко проверяемых условий устойчивости. С этой целью для решения в виде нечеткого числа или вектора нечетких чисел можно дать следующее определение устойчивости.

Определение 1. Задача Е(а(а),х) = 0, содержащая параметры, значения которых заданы нечеткими промежутками а° = [аД0); аД 1)], г = 1,...,п; устойчива по решению, если:

- она имеет решение для всех значений параметров а Е 01 х ... х 0°;

- для, любого заданного е > 0 найдется ё(е > 0) такое, что при выполнении условия \а1 — а2\ < ё, для всех а Е [0; 1], выполняется условие - для всех компонент х(а) решения задачи \х(а1) — х(а2)\ < £ (рис. 3).

Рис. 3. Иллюстрация к определению устойчивости

Иными словами, необходимым условием устойчивости является непрерывность

а

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

Определение 2. Задача, Е(а(а),х) = 0 имеет устойчивое решение если:

- она имеет решение для всех значений параметров на всех а-уровнях при а Е [0; 1];

- каждая компонента а(а)решения непрерывная функция на промежутке а Е [0; 1].

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

Предложение 2. Для того чтобы, нечеткое модифицированное решение х(а) = X(0) + (X(1) —x(0))— квадратного уравнения а(а)ж2+6(а)ж+с(а) = 0 было устойчивым, достаточно существования решений X(0), х(1).

Доказательство этого утверждения следует из простых рассуждений. Если существуют решения х(0), х(1), то необходимо следует существование неотрицательного дискриминанта D(a) = Ъ2(а) — 4а(а)Л(а) при а = 0 и а=1. Но в силу линейности операций (6) дискриминант - монотонная функция параметра а. Значит ИЗ ВЫ ПОЛ-нения неравенств D(0) > 0 и D(1) > 0 следует D(a) > 0 при всех а Е [0; 1], что и

а

, . —Ъ(а) ± (D(a))0'5

Решение х(а) = -——- непрерывная функция на промежутке а Е

2а (а)

[0; 1] и D(—) > 0 в силу линейности операций (6) - (9). Полученные результаты полностью соответствуют определению 2, что и доказывает утверждение 2.

Предложение 3. Нечеткое модифицированное решение хДа) = Xi (0) + (хД1) — Хг(0))а; i = 1,..., n задачи линейного программирования канонического типа, с нечеткими параметрами

n

Е aiXi ^ min;

i=l n

Е bijXi = Cj(а); j = 1, ...m; m < n;

i=1

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

J2hjXi = Cj(0) и J^bijXi = Cj(1).

ii

Действительно, если решения х(0) = (х1(0);...; Xn(0), х(1) = (х1(1);...; Xn(1)) существуют и принадлежат одноименным вершинам, то им соответствуют две системы уравнении относительно базисных переменных Xi, i = 1, ...,m

m n _

Y^bijXi = — E bijXi + ALj(0), j = 1,...,m;

i=1 i=m+1

mn

EbijXi = — E bijXi + Cj(1), j = 1, ...,m.

i=1 i=m+1

п - -fn\ rA1(0) An(0),

Представим решения этих уравнении в виде x(0) = (———;...; ———) и x(1) =

(———...; ———), где символом — обозначены определители правила Крамера. При —=0 аЕ [0; 1] будут существовать решения x(—) =

(————);...; ——-—). Непрерывность компонент вектор а решения X(а)oбycлoвлeнa линейностью операций (6) - (9). Поскольку каноническая форма задачи линейного программирования сводится к стандартной и общей форме, утверждение справедливо и для этих представлений задачи.

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

Пример 3*. В примере 3 изменим нечеткий параметр следующим образом: А

(-1; 0; 1, 6).

Уравнение с модифицированным параметром (А=0,5)

dy(a)

dt

(0, 3 - 0, 3a)y(a).

Дифференциальные уравнения при а=1 и а=0 примут вид:

dy (1)

0

dy(0)

_ ¿Ь _ ¿Ь

0,3у(а). Соответствующие решения будут иметь вид: у(1) = 1 и у(0) = ехр(0, 3Ь).

Очевидно, что первое решение устойчивое, второе не устойчивое по Ляпунову по

переменной Ь. Сформируем частное решение в виде выпуклой линейной комбинации

аа

нечеткой функции: у(Ь, а) = ехр(0, 3Ь) + (1 — ехр(0, 3Ь))а.

Рис. 4. Неустойчивое решение задачи Копти с нечеткими параметрами

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

¿У

Определение 3. Задача Кош,и — = а(а)у с начальными условиям,и у(Ь = 0) = у0,

аЬ

содержащая параметр, значения которого заданы нечетким промежутком а° = [аД0); аД1)], г =1, ...,п; устойчива по решению, если:

- она имеет решение для всех значений параметра а(а), а Е [0; 1];

- для, любого заданного е > 0 найдется ё(е) > 0 такое, что при выполнении условия \а1 — а2\ < ё, для, всех а Е [0; 1] и всех значений Ь, выполняется условие \у(Ь, а1) — у(Ь, а2)\ < е, где у(Ь, а) - реш,ение задачи.

Предложение 4. Задача Кошм, сформулированная в определении 3, будет иметь устойчивое решение, если будут существовать устойчивые по Ляпунову решения, аа

Это предложение достаточно очевидно, так как линейная форма нечеткого числа позволяет сохранить условия существования решения внутри промежутка [0; 1], если

существуют решения при а=1 и а=0, а устойчивость этих решений по Ляпунову обеспечит выполнение второго условия определения 3.

Заключение

В результате проведенных исследований получена обоснованная методика решения ряда задач выбора и управления, представленных моделями с параметрической нечеткостью. Методика состоит из следующих этапов:

1. В рассматриваемой модели каждый нечеткий параметр, заданный как нечеткое число LR-типа представляется в виде совокупности а-уровневых интервалов в соответствии с (2).

2. К границам интервалов применяется выпуклое линейное преобразование (5) с фиксированными параметрами А, в результате чего нечеткие параметры приобретают вид модифицированных нечетких чисел LL/RR-типа (4), отражающих основные, с точки зрения выбора, характеристики нечеткости.

3. В исходной модели нечеткие параметры переписываются в виде модифицированных чисел, и задача решается для границ промежутка изменения значений а G {0; 1}, как обычная задача с четкими параметрами. В результате получается два четких решения.

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

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

Предложенная методика обеспечивает процесс получения решения:

- с ограниченным расширением неопределенности, обусловленным применяемой для вычислений алгеброй;

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

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

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

Работа выполнена при финансовой поддержке гранта РФФИ № 13-08-00532. Литература

1. Fuzzy Theory Systems: Techniques and Applicatons / Ed. by Cornelius T. Leondes. - London:

Academic Press, 1999. - 1777 p.

2. Intelligent Systems for Information Processing: From Representation to Applications / Eds. by Bouchon-Meunier В., Foulloy L., Yager R.R. - Amsterdam: Elsevier, 2003. - 488 p.

3. Hanss, M. Applied Fazzy Arithmetic: An Introduction with Engineering Applications / M. Hanss. - Netherlands: Springer, 2005. - 256 p.

4. Нечеткие множества в моделях управления и искусственного интеллекта / А.Н. Авер-кин, И.З. Батыршин, Л.Ф. Блишун и др.; под ред. Д.А. Поспелова - М.: Наука. Глав, ред. физ-мат. лит., 1986. - 312 с.

5. Борисов, A.M. Принятие решений на основе нечетких моделей: примеры использования / A.M. Борисов, О.А. Крумберг, И.П. Федоров. - Рига: Знание, 1990. - 184 с.

6. Piegat, A. Fuzzy Modelling and Control / A. Piegat. - New York: Springer-Heidelberg, 2001. -371 p.

7. Chen, S.P. A Simple Approach to Fuzzy Critical Path Analysis in Project Networks / S.P. Chen, Y.J. Hsueh // Applied Mathematical Modeling. - 2008. - V. 32. - P. 1289-1297.

8. Bekheet, S. An Enhanced Fuzzy Multi Criteria Decision Making Model with a proposed Polygon Fuzzy Number/ S.S. Bekheet, A. Mohammed, H.A. Hefny // International Journal of Advanced Computer Science and Applications. - 2014. - V. 5, № 5. - P. 118-121.

9. Воронцов, Я.А. Алгебраические операции с нечеткими LR-числами с использованием преобразования L / Я.А. Воронцов, М.Г. Матвеев // Программная инженерия. - 2014. -№ 8. - С. 23-29.

10. Воронцов, Я.А. Арифметические операции над двухкомпонентыми нечеткими числами / Я.Ф. Воронцов, М.Г. Матвеев // Вестник ВГУ. Серия: Системный анализ и информационные технологии. - 2014. - № 2. - С. 75-82.

11. Ашманов, С.А. Линейное программирование / С.А. Ашманов. - М.: Наука, Гл. ред. физ.-мат. лит., 1981. - 340 с.

12. Агаян, Г.М. О задаче линейного программирования с приближенными данными / Г.М. Агаян, А.А. Рютин, А.Н. Тихонов (мл.) // Журнал вычислительной математики и математической физики. - 1984. - Т. 24, № 9. - С. 1303-1311.

Михаил Григорьевич Матвеев, доктор технических наук, профессор, кафедра «Информационные технологии управления:», Воронежский государственный университет (г. Воронеж, Российская Федерация), [email protected].

Поступила в редакцию 5 июня 2015 г.

MSC 06D72, 08A72 DOI: 10.14529/mmpl50402

Analyzing and Solving Problems of Decision Making with Parametric Fuzzy

M. G. Matveev, Voronezh State University, Voronezh, Russian Federation, mgmat veev@yandex. ru

The method of solution of decision making problems presented as models with parameters in the form of LR fuzzy numbers is proposed. This methodic is based on using

of a-level representation of fuzzy numbers, their subsequent modification by a convex linear transformation of the boundaries of a-intervals, preserving the basic characteristics of fuzziness, proposed algebra of modified fuzzy numbers and a convex linear combination of the boundaries of a-change interval. Bounded growth of uncertainty in fuzzy information processing, preservation of natural interpretation of intermediate and final results of calculations and the possibility of algorithm realization in software environments working

a

representation causes the problem of fuzzy solutions stability. We give the definition of stability for solutions in the form of el fllZZy point in n-dimensional space and in the form of a fuzzy function. For several kinds of problems we give a stability criteria, which is easily verified in practical calculations. We have examples of solving the problems with parametric fuzziness using the proposed method, confirming the validity of the results.

a

representation; algebra of fuzzy numbers; stability of fuzzy solution.

References

1. Fuzzy Theory Systems: Techniques and Applicatons. Ed. by Cornelius T. Leondes. London, Academic Press, 1999. 1777 p.

2. Intelligent Systems for Information Processing: From Representation to Applications. Amsterdam, Elsevier, 2003. 488 p.

3. Hanss M. Applied Fazzy Arithmetic: An Introduction with Engineering Applications. Netherlands, Springer, 2005. 256 p.

4. Averkin A.N., Batyrshin I.Z., Blishun A.F. etc. Nechetkie mnozhestva v modelyakh upravleniya i iskusstvennogo intellekta [Fuzzy Sets in Management Models and Artificial Intelligence]. Moscow, Nauka, 1986. 312 p.

5. Borisov A.M., Krumberg O.A., Fedorov I.P Prinyatie resheniy na osnove nechetkikh modeley: primery ispol'zovaniya [Decision-Making Based on Fuzzy Models: Examples of UseJ. Riga, Nauka, 1990. 184 p.

6. Piegat A. Fuzzy Modelling and Control. New York, Springer-Heidelberg, 2001. 371 p. DOI: 10.1007/978-3-7908-1824-6

7. Chen S.P., Hsueh Y.J. A Simple Approach to Fuzzy Critical Path Analysis in Project Networks. Applied Mathematical Modelling, 2008, vol. 32, pp. 1289-1297. DOI: 10.1016/j.apm.2007.04.009

8. Bekheet S.S., Mohammed A., Hefny H.A. An Enhanced Fuzzy Multi Criteria Decision Making Model with a Proposed Polygon Fuzzy Number. International Journal of Advanced Computer Science and Applications, 2014, vol. 5, no. 5, pp. 118-121. DOI: 10.14569/IJ ACS A.2014.050517

9. Vorontsov Y.A. Matveev M.G. Algebraic Operations with LR Fuzzy Numbers Using L Conversion. Software Engineering, 2014, no. 8, pp. 23-29. (in Russian)

10. Vorontsov Y.A. Matveev M.G. Arithmetic Operations on Two-Component Fuzzy Numbers. Proceedings of Voronezh State University. Series: System Analysis and Information Technologies, 2014, no. 2, pp. 75-82. (in Russian)

11. Ashmanov S.A. Lineynoe programmirovanie [Linear Programming]. Moscow, Nauka, 1981. 340 p.

12. Agayan G.M., Ryutin A.A., Tikhonov A.N. The Problem of Linear Programming with Approximate Data. USSR Computational Mathematics and Mathematical Physics, 1984, vol. 24, no. 5, pp. 14-19. DOI: 10.1016/0041-5553(84)90149-6

Received June 5, 2015