CC BY
55РЕШЕНИЕ НЕЧЕТКОЙ МНОГОКРИТЕРИАЛЬНОЙ ЗАДАЧИ ОПТИМИЗАЦИИ В УСЛОВИЯХ РИСКА
Т. Ф. Бекмуратов, Д. Т. Мухамедиева
Центр разработки программных продуктов и аппаратно-программных комплексов при Ташкентском университете информационных технологий, 100084, Ташкент, Узбекистан
УДК 519.71(575.1)
Рассматриваются подходы к решению нечеткой многокритериальной задачи оптимизации в условиях риска.
Ключевые слова: нечеткое множество, многокритериальная задача, оптимизация, риск, адаптивный метод, принятие решения.
An approach to solving the optimization problem and the solution of the fuzzy multiobjective optimization problem in terms of risk.
Key words: fuzzy set, multicriteria task, optimization, risk, adaptive method, decision making.
Введение. Существует большое количество работ, в которых исследуется проблема принятия решений в условиях риска [1-5]. Широко используются минимаксный подход, оптимизация ожидаемой полезности, минимизация среднего ущерба или вероятности неблагоприятного события, модели стохастического программирования и др. Изучению вопросов формализации рисков в экономическом плане (оценка риска инвестиции, принятия проектов) посвящены работы А. О. Недосекина [6]. В этих работах предлагается своего рода конструктор, который может быть использован любым экспертом по своему усмотрению. Идея применения нечетких множеств при проведении финансового анализа предприятий предложена в [6] как способ избежать неопределенности, когда она имеет не только статистический, но и лингвистический характер.
Однако данных постановок недостаточно для принятия решений в нечетких условиях, когда невозможно ориентироваться на средние показатели эффективности решений, поскольку они оправданы в случае многократно повторяющихся ситуаций, в то время как рисковые ситуации уникальны: они могут возникнуть в любой момент или не возникнуть никогда. Последние ситуации характеризуются возможностью крайне мало вероятных, но исключительно больших потерь, способных создать проблему выживания рассматриваемой системы. Ясно, что такие традиционные показатели риска, как дисперсия, в данном случае неадекватны.
Для оценки рисков в нечетких условиях предлагается дополнить систему ограничений стандартной задачи принятия решений набором ограничений по возможным потерям, а именно построить для избранных сценариев модель их последствий (ущербов) как функций управляющих параметров и накладывать экспертные ограничения на приемлемый уровень относительного ущерба для каждого сценария.
1. Подходы к решению оптимизационной задачи. Очевидно, что найти идеальный вариант критерия "максимальная доходность — минимальный риск" удается лишь в очень
Подходы к решению оптимизационной задачи
Подход
Модель
Подход "максимум выигрыша": из всех вариантов выбирается тот, который позволяет получить максимальное значение выигрыша Г при приемлемом для лиц, принимающих решения (ЛПР), значении Кпрлт,
Подход "оптимальная вероятность": из возможных решений выбирается такое решение, при котором обеспечивается максимум математического ожидания выигрыша М(Г) при приемлемом для ЛПР риске К
Сочетание подходов "оптимальная вероятность" и "оптимальная колеблемость": выбирается вариант, обеспечивающий минимальное значение коэффициента вариации СУ (Г) при приемлемом для ЛПР риске К
Подход "минимум риска": из всех возможных вариантов выбирается тот, который позволяет получить ожидаемый выигрыш, т. е. предельно допустимое значение Г при минимальном риске
Подход "максимальная доходность — минимальный риск": выбирается вариант, обеспечивающий максимум Г при минимуме К в случае удовлетворения заданным ограничениям
Г — тах,
К Кnp.don,
с к
з
М (Г) — тах, Е ХК3 С к
з
СУ (Г) — тах,
Е ХК3 С К
з
Г Гnp.don,
К — тт,
Е ХКз С к
з
Г — тах, К — тт,
ЕХк с к
редких случаях. Поэтому предлагаются следующие подходы к решению этой оптимизационной задачи (см. таблицу). В таблице К = (у\у Е Кт, у < д} — заданное выпуклое подмножество пространства Кт
УЛкК = < У\У Е Кт, у = ^\аъз ,хг Е X, а^ Е К у с Кт, г = !,...,п; з = 1,...,т
2. Решение нечеткой многокритериальной задачи оптимизации. В многокритериальных задачах сложно оценить решение задачи с учетом всех критериев. Наиболее распространенными являются метод аддитивной свертки и оценка лицом, принимающим решение (ЛПР) [7-9].
Предлагается применять нечеткие методы при оценке альтернативных решений. Решение многокритериальной задачи оптимизации включает следующие этапы [8-11]:
— формирование целевой функции в нечеткой постановке;
— определение значений критериев оценки в нечетком виде;
— разработка функций принадлежности для критериев;
— определение базы правил и (или) базы предпочтений для критериев;
— вычисление значений целевой функции;
— дефазификация (приведение к четкому виду) целевой функции.
Обозначим через Г(X, Л) операцию свертки частных критериев оптимальности, где Л Е Ба С К — вектор весовых множителей (Л = {Л^}, г = 1,^); = (Аг\Аг > 0, ЕАг = 1, г Е [1 : в]} — множество допустимых значений этого вектора.
Задача параметрического программирования с в независимыми параметрами Л = {Лг}, г = 1,5, или 5-задача параметрического программирования, в матричном виде записывается следующим образом:
Г(X, Л) = (ао + Л'Ь^ + еЛ ^ ех1г, ^ХК С К> Л е Я
Здесь X = {х-}, ] = 1,п — решение 5-задачи параметрического программирования; а0,Ь,е — коэффициенты, являющиеся нечеткими подмножествами числовой оси базовых множеств А0,В0,Е0 соответственно, представляемыми обычно в виде нечетких множеств с заданными функциями принадлежностей Цщ(а0) (а0 С А0), ^ь(Ь) (Ь С Б) и це(е) (е С Е).
Для решения задачи параметрического программирования при нечетких исходных данных (коэффициентов а0 ,Ь,е ) предлагается использовать три подхода:
1. Получение с использованием различных операций дефазификации нечетких множеств а0,Ь,е (интегрирования, суммирования, осреднения и др.) нечетких оценок значений коэффициентов а0,Ь,е [8]. Введя эти значения в 5-задачу параметрического программирования вместо нечетких коэффициентов и записав ограничения в виде соответствующих неравенств, исходную задачу сведем к виду
Г(X, Л) = (а0 + Л' Ь)Х + еЛ —► ех1г(шт), ^ХК < д, Л е Я3. (1)
Заметим, что в силу нечеткости описания коэффициентов а0 и Ь оценка любого решения х(Л) е X (и соответственно значения функции Г(X, Л) при х = х(Л)) представляет собой нечеткое подмножество числовой оси базового множества Х.
2. Сведение решения исходной задачи к решению задач линейного программирования для каждого дискретного а-уровня [9]. В результате нечеткие ограничения записываются в интервальном виде
Р =) га(ац)х1 + аа(а12)х2 + ... + Са(аы)хп С Са(Ы), г =1,т,а =1,р,
1'3
> ] = 1,п.
Здесь X = {х-}, ] = 1,п — решение многокритериальной задачи параметрического программирования на каждом дискретном а-уровне; аа(аг,-), оа(Ьг) — интервальные значения коэффициентов аг- и Ьг на каждом дискретном а-уровне.
3. Решение задачи многокритериальной оптимизации адаптивным методом [7]. Каждая итерация этих методов включает фазу анализа, выполняемую ЛПР, и фазу расчетов, выполняемую системой многокритериальной оптимизации.
Прямой адаптивный метод решения многокритериальной задачи, который исследуется в данной работе, основан на предположении существования функции предпочтений ЛПР Г(X, Л) = (а0 + Л'Ь^ + еЛ, которая определена на множестве Дх допустимых значений вектора варьируемых параметров X и выполняет отображение этого множества на множество действительных чисел К. При этом задача многокритериальной оптимизации сводится к задаче выбора вектора X* е Дх (X* = {х*}, ] = 1, п), такого что
min F(X, Л) = F(X*, Л), X е Dx. (2)
При каждом фиксированном векторе метод скалярной свертки сводит решение задачи (1) к решению однокритериальной задачи глобальной условной оптимизации (2).
Заметим, что в случае аддитивной свертки F(X, Л) вектор X принадлежит множеству эффективных по Парето векторов [9].
Указанное обстоятельство позволяет полагать, что в данном случае функция предпочтений ЛПР определена не на множестве Dx, а на множестве Da'- F : Л ^ R. В результате многокритериальная задача сводится к задаче выбора вектора Л* е Da, такого что
minF(X, Л) = F(X, Л*), Л е Da. (3)
л
Поскольку обычно s « n, переход от задачи (1) к задаче (3) важен с точки зрения уменьшения вычислительных затрат.
Вектор Л* е Dл находится с помощью нечетких правил вывода:
k / в \
U I П V = " с весом wjp\ ^ F(X, Л) = F(X, Л*). p=1 Дг=1 ) /
Здесь Фi,jp — лингвистический терм, которым оценивается переменная \i в строке с номером jp; Фjp — весовой коэффициент правила с порядковым номером jp; F(X, Л) = F(X, Л*) —
выход нечеткого правила. Величину Ф будем считать лингвистической переменной со значе-
*
ниями от "Очень-очень плохо" до "Отлично". Ядро нечеткой переменной Ф обозначим Ф [8],
*
так чтобы значению переменной Ф "Очень-очень плохо" соответствовало ф = 1, а значению
*
"Отлично" — Ф = l.
В результате многокритериальная задача сводится к задаче отыскания вектора Л* е Dл,
*
обеспечивающего максимальное значение дискретной функции Ф(Л):
*
тахф(А) = Ф(А*) = Ф*, Л е Бл. Л
Каждая входная переменная имеет собственные функции принадлежности нечетким термам Фур .
Функция принадлежности элемента \ терму Фур имеет вид
азр(\ ) = _1_
1 1 + [(А - ь?)/сГ]2'
где Цр, с1р — параметры функции принадлежности.
Общая схема рассматриваемого метода является итерационной и включает следующие перечисленные ниже основные этапы.
Этап 1. Случайно последовательно генерируются з векторов А1, А2,... , А3, и для каждого из этих векторов выполняются следующие действия: 1) решается многокритериальная задача
шт Г(X, Л) = Г(X*, Л), X € Бх;
X
(4)
2) ЛПР предъявляет найденное решение X*, а также соответствующие значения всех частных критериев оптимальности /1 (X*), /2(X*),... , /s(X*);
3) ЛПР оценивает эти данные и вводит в задачу соответствующее значение своей функ-
ции предпочтений ф(Лг).
Этап 2. На основе всех имеющихся значений А1, А2,... , Аs вектора Л и соответствующих оценок функции предпочтений выполняются следующие действия:
1) строится функция Г1 (X, Л), аппроксимирующая функцию Г(X, Л) в окрестности то-
3) с найденным вектором Л* решается задача вида (4) — находятся вектор параметров и соответствующие значения частных критериев оптимальности, а затем предъявляются ЛПР, которое оценивает указанные данные и вводит в систему соответствующее значение своей функции предпочтений Г (Л*).
Этап 3. На основе всех имеющихся в системе значений А1, А2, ... , Аs вектора Л и соответствующих оценок функции предпочтений Г(X, Л1) выполняется аппроксимация функции Г(X, Л) в окрестности точек Л1; Л2,... , Л&, Л1, строится функция Г2(X, Л) по схеме первого этапа, до тех пор пока ЛПР не примет решение о прекращении вычислений.
Входами системы нечеткого вывода являются значения весов частных критериев оптимальности — нечеткие термы Аг,, г = 1, 2,..., к € [1 : в]. Выходной переменной системы нечеткого вывода является лингвистическая переменная, ядро которой принимает значения
Совокупность значений указанных нечетких входных переменных, выходных лингвистических переменных, а также правил нечетких продукций образуют нечеткую базу знаний.
3. Настройка нечетких баз знаний с применением генетических алгоритмов оптимизации. Алгоритм настройки параметров функций принадлежности И = (Ь1,Ь2,... ,ЪЯ), С = (с1, с2, ... , оя) и W = (ш1, ш2, ... , ) включает следующие этапы.
1. Формирование исходной популяции. Для реализации генетического алгоритма следует задать способ кодирования нечетких моделей. Сведем неизвестные параметры W, В, С в один вектор (см. рисунок):
5 = В, С) = (^1, шг, ... , Ьц} сп,..., Ьц1} сщ, Ьп1, Сп1,..., Ьпгъ спь1).
Здесь N — общее число строк в нечеткой базе знаний; ¡г — количество термов-оценок входной переменной Аг; ¡1 + ¡2 + ... + ¡п = д, г = 1, п; д — общее число термов.
2. Скрещивание. Поскольку операция скрещивания является основной операцией генетического алгоритма, его производительность в первую очередь зависит от производительности используемой операции скрещивания. В результате скрещивания двух хромосом-родителей 51 и 52 получаются хромосомы-отпрыски СН1 и Ск2 путем обмена генов относительно (п + 1)-й точки скрещивания.
3. Мутация. Каждый элемент вектора 5 может подвергнуться операции мутации с вероятностью Рт. Обозначим мутацию элемента з через Ми(в):
*
А
1, 2, . . . , ¡.
Xi Xn Веса правил
r л
bii С11 bi c11, bni Cni К Cn1n w1 Wn
Терм 11
Терм И1 Терм щ Терм н1„
Код нечеткой модели — хромосома
Mu (Wj) = RANDOM ([;w, w]) ,j = 1, N,
Mu (bip) = RANDOM ([xi, x]) ,
Mu (ap) = RANDOM ([c,,^]) ,
где w,W — нижняя и верхняя границы интервала возможных значений весов правил [w, W] С [0,1]; С, (c-i)^ — интервал возможных значений коэффициента концентрации-растяжения функций принадлежности термов-оценок входной переменной x,, (c-i)^ С (0, ,i = 1,n; RANDOM (£]) обозначает операцию нахождения равномерно распределенного на интервале ^ случайного числа.
Разработано программное обеспечение в среде DELPHI-7 для решения практических задач многокритериальной оптимизации и получены результаты решения оптимизационной задачи [10, 11].
Заключение. Таким образом, показана целесообразность объединения метода нечеткого вывода и генетических алгоритмов в задачах с неопределенной или лингвистической информацией, а также в задачах, для которых характерны интуитивные решения. Использование предложенного метода позволяет существенно улучшить качество решения многокритериальных задач оптимизации с нечетко заданными параметрами и критериями. В дальнейшем планируется изучение различных гибридных методов применительно к оптимизационным задачам, а также методов автоматического формирования базы нечетких правил, что позволит значительно усовершенствовать процесс автоматизации.
Список литературы
1. Норкин В. И. Об измерении и профилировании катастрофических рисков // Кибернетика и систем. анализ. 2006. № 6. С. 80-94.
2. Ermoliev Yu. M., Ermolieva T. Yu. Catastrophic risk management: flood and seismic risks case studies / Applications of stochastic programming / Ed. by S. W. Wallace, W. T. Ziemba. Philadelphia: MPS-SIAM, 2005. P. 425-444.
3. Михалевич В. С., Кнопов П. С., Голодников А. Н. Математические модели и методы оценки риска на экологически опасных производствах // Кибернетика и систем. анализ. 1994. № 2. С.121-138.
4. Кнопов П. С., Марьянович Т. П. О некоторых актуальных проблемах оценки риска сложных систем в условиях недостаточной информации // Кибернетика и систем. анализ. 2003. № 4. С. 125-137.
5. Сергиенко И. В., Яненко В. М., Атоев К. Л. Общая концепция управления риском экологических, техногенных и социогенных катастроф // Кибернетика и систем. анализ. 1997. № 2. С. 65-87.
6. Недосекин А. О. Нечетко-множественный анализ риска фондовых инвестиций. СПб.: Сезам, 2002.
7. Lee C. C. Fuzzy logic in control systems: fuzzy logic controller. 1 // IEEE Trans. Systems, Man Cybernet. 1990. V. 20, N 2. P. 419-435.
8. Алтунин А. Е. Модели и алгоритмы принятия решений в нечетких условиях / А. Е. Алту-нин, М. В. Семухин. Тюмень: Тюмен. гос. ун-т, 2000.
9. Орловский С. А. Проблемы принятия решений при нечеткой исходной информации. М.: Наука, 1981.
10. Бекмуратов Т. Ф. Систематизация интеллектуальных систем поддержки принятия решений // Пробл. информатики и энергетики. 2003. № 4. С. 25-35.
11. Мухамедиева Д. Т. Задачи нечеткого параметрического программирования в случае зависимости от многих параметров коэффициентов целевой функции // Вопр. вычисл. и прикл. математики (Ташкент). 2004. Вып. 114. С. 81-87.
Бекмуратов Тулкун Файзиевич — акад. АН РУз, д-р техн. наук, гл. науч. сотр.
Центра разработки программных продуктов и аппаратно-программных комплексов при Ташкентском университете информационных технологий;
тел. (+99871) 262-71-53; e-mail: bek.tulkun@yandex.ru;
Мухамедиева Дилноз Тулкуновна — д-р техн. наук, вед. науч. сотр.
Центра разработки программных продуктов и аппаратно-программных комплексов при Ташкентском университете информационных технологий;
тел. (+99871) 262-71-55, e-mail: dilnoz134@rambler.ru
Дата поступления — 24.09.12
