ПАРАМЕТРИЧЕСКИЙ α -УРОВНЕВЫЙ МЕТОД λ -ПРОДОЛЖЕНИЯ ДЛЯ ЗАДАЧИ НЕЧЕТКОГО ЛИНЕЙНОГО ПРОГРАММИРОВАНИЯ Текст научной статьи по специальности «Математика»

УДК 519.688

DOI: 10.18101/2304-5728-2018-1-34-51

ПАРАМЕТРИЧЕСКИЙ а -УРОВНЕВЫЙ МЕТОД Л -ПРОДОЛЖЕНИЯ ДЛЯ ЗАДАЧИ НЕЧЕТКОГО ЛИНЕЙНОГО ПРОГРАММИРОВАНИЯ

© Шаталова Алевтина Юрьевна

аспирант,

Кубанский государственный университет

Россия, 350040, г. Краснодар, ул. Ставропольская, 149

E-mail: al-shatalova@yandex.ru

© Лебедев Константин Андреевич

доктор физико-математических наук, доцент

заведующий кафедрой,

Кубанский государственный университет

Россия, 350040, г. Краснодар, ул. Ставропольская, 149

E-mail: klebedev.ya@yandex.ru

В статье рассмотрен известный метод четырех задач нечеткого линейного программирования, который обобщен до параметрического а -уровневого метода Л -продолжения. Предложенный метод позволяет сократить количество рассмотренных случаев с четырех до двух, позволяет получать дополнительные решения на подмножествах квадрата а - Л. Преимущества предложенного подхода состоят в том, что они включают в себя четыре описанные ранее задачи как частные случаи. Метод позволяет получать более общие и гибкие решения, дает возможность влиять на свойства линейной задачи оптимизации: меру принадлежности, чувствительность, устойчивость. В статье также рассмотрен численный пример, показывающий работу нового метода и реализованный в программной среде МаШСАБ.

Ключевые слова: нечеткое линейное программирование; оптимизационные задачи; параметрическое линейное программирование с нечеткими данными; принцип расширения; метод продолжения; нечеткие множества; нечеткая логика.

1. Введение

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

При разработке инструментария для нечеткого подхода используются различные методы. Задачи, в которых нечетко описано множество допустимых выборов (альтернатив) и (или) функции полезности, принимают вид задач нечеткого математического программирования [1]. Существует два различных подхода к их решению:

1) подход Беллмана — Заде — задача формулируется как задача выполнения нечеткой цели при нечетких ограничениях, причем решение задается пересечением нечетких множеств цели и ограничений [2].

2) выбираются эффективные альтернативы в понимании Парето [1]. Необходимо отметить, что нахождение полного решения некоторых

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

В работах [3; 4] представлены методы линейной оптимизации с нечеткими данными, позволяющие находить оптимальные решения задач линейного программирования с коэффициентами в виде нечетких чисел путем сведения их к параметрическим задачам линейного программирования.

В [3] задачи рассматривались с коэффициентами в виде нечетких треугольных чисел, имеющих левую (Ь), правую (Я) стороны и середину (С), поэтому в линейной задаче выражения для целевой функции и ограничений могут быть записаны в форме левой (Ь) или правой (Я) стороны (рис. 1) с использованием среднего значения (С) [5].

м

Рис. 1. Треугольное число а = (аЬ,аС,аК)

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

а = (aL,aC,ciR), b = (bL,bC,bR)

аС < ЬС — сравнение, аналогичное сравнению действительных чисел;

aL<bR — «слабое» отношение;

aR<bL — «сильное» отношение. Данные обозначения вводятся в [3].

Таким образом, возникает два расширения бинарного отношения нера-

min

венства «меньше или равно» (<): «слабое» — А < В . (А' <ВИ) и «силь-

шах

ное» — А < В, (А <В ), где An В — линейные выражения относительно неизвестных с нечеткими треугольными коэффициентами (в частности, треугольное число). Так как задача линейного программирования содержит отношение неравенства (или равенства) в ограничениях, а целевая функция может быть записана в форме левой (L) или правой (R) стороны коэффициентов, то метод задает четыре связанные между собой параметрические задачи. В статье [5] рассмотрен пример применения данного подхода к экономической задаче оптимального инвестирования.

В данной статье, в отличие от [3; 5; 6], предлагается усовершенствовать метод: ввести дополнительный параметр Ле [0,1], и вместо четырех связанных между собой задач с а е [0,1] рассматривать две, но с двумя параметрами а,/1е[0Д]. Параметр Л назовем параметром «гибкости», а сам метод — параметрическим а -уровневым методом Л -продолжения задачи нечеткого линейного программирования.

Построена модель и создана программа с помощью пакета MathCAD, позволяющая решить задачу линейной оптимизации при нечетких коэффициентах в ограничениях и целевой функции. Параметрический а -уровневый метод Л -продолжения задачи нечеткого линейного программирования позволяет получить решение задачи при различных значениях а,/1 е [0,1].

2. Метод параметризации задачи

В работе [3] приведен параметрический а -уровневый метод задачи нечеткого линейного программирования с использованием а -эффективного решения, где значение а является уровнем эффективности. Метод позволяет для различных значений ае[0,1] и

min

четырех разных задач с использованием «слабого» — А < В и

шах

«сильного» — А < В отношений построить семейство задач, связанное с классическим линейным программированием.

Параметр ае[0,1] имеет значения в отрезке от 0 до 1 (ае[0Д]). Используя метод погружения задачи в семейство задач с параметром продолжения Ае[0Д], можно показать, что методы, изложенные в [3], позволяют решать задачи только при Л = 0 или Л = 1, что в системе

координат (рис. 2а) отображается двумя отрезками Ы^-КЬ на верхней и нижней сторонах квадрата, а другая область квадрата недоступна. Предлагаемый здесь метод позволяет рассматривать точки и в середине области квадрата, например, со значением а = 0,7 и А = 0,5 (рис. 26), в промежутке между слабым и сильным отношениями.

X

1Я

Ей

И.

О

>

о

1 £Г

Рис. 2а. Вид области определения параметров а - 0; 1,Л = 0; 1

о

>

о

1 а

Рис. 26. Вид области определения параметров а е (0; 1), А = (0; 1)

Задача линейного программирования с нечеткими данными имеет вид

[3]:

/(х) = (с,х)->тш(тах), х е X = {х е Еп : Ах < Ь], (1)

где с — вектор с элементами в виде нечетких треугольных чисел с (/ = 1, п ), А — матрица, состоящая из нечетких треугольных чисел ау = ^ауС• а1/И) (I = ■ / = \.п) с функциями принадлежности вида

г-а}- а*-г

¡иа (0 = тах{0,тт(—-—1—-;—--—)} на вещественном носителе /?.

"г ~ау А _

Ь — вектор с элементами в виде нечетких треугольных чисел bj (/ = 1, т ) со своими функциями принадлежности на носителе /?.

Обозначим через ЬЫ лево-право-стороннее: {Ах)[л < Ь'1: а через ЫЬ право-лево-стороннее отношение: {Лх)ил < Ь'А. где А' . А" . Ьь, 6й вычисляются по формулам, вывод которых очевиден из рис. 3 [3]:

0 а 1 ^

Рис. 3. Построение функций принадлежности треугольного числа

а = {аЬ,аС ,аК)

Аь =аЬ + (аС - аЬ)а ,

=а^-(а^-аС)а , (2)

Ъь =Ы + (ЬС-Ы)а ,

а = {аЬ,аС,аК), Ь = (ЬЬ,ЬС,ЬЬ) — треугольные числа (нижние индексы опущены). Параметр ае[0,1] можно трактовать как «показатель эффективности», т. к. при а = 1 коэффициенты принимают значения

действительных чисел, а при а = 0 — все значения из сечения А0 треугольного числа.

Коэффициенты целевой функции также представляются в виде нечетких чисел вида: с = (cL,cC,cR).

Из равенств (2) видно, что при а е [0,1] коэффициенты ограничений принимают значения на некотором множестве в зависимости от выбранного а и при (cTx)La —> min или (с'х)" —> min задача линейного программирования в этом случае делится на четыре возможных варианта (табл. 1) [3]. Понятно, что задачу можно рассматривать и на максимум (max) целевой функции, что не приведет к принципиальным отличиям.

Таблица 1

Вид целевой функции Вид ограничений

(cTx)La min(max) (Ах)ьа<Ъ* , (1 задача)

(Ax)Ra<bLa, (2 задача)

(c''xt min(max) (Ах)1'; <//\ (3 задача)

(Ах)1<%, (4 задача)

Данный подход с четырьмя задачами позволяет получить оптимальные решения для четырех случаев (четырехсторонний вариант) на двух сторонах квадрата (рис. За). Множество точек во внутренней области квадрата не может быть использовано. Развитие метода заключается в погружении задачи в семейство задач с параметром Л, который дает возможность решить задачу и для любой точки квадрата, любого сочетания значений а и 1 из квадрата [0,1] х [ОД] (рис. 26). Например, для точки (а = {).! . /1 = 0,5) или серии задач при (а = 0,7 , Л е [ОД]) или (Л = 0,5, а е [ОД]) или вообще на любом подмножестве из области квадрата.

Предложим метод погружения исходной проблемы с 4 задачами в семейство задач с параметром продолжения. Запишем ограничения задачи в виде:

Л(Ах% + (1 -1) • (Ах)яа < ЛЪка + (1 -1) • Ъьа ,

где Л е [ОД].

При Л = 0 получается «сильное» ограничение (Ах)^ < Ь'А. а при Л = \ — «слабое» (Ах)'^ < //'. Также и целевая функция получает два случая с Я = 0 и 1 = 1. Параметр Л позволяет непрерывно, через промежуточные задачи, переходить от задачи ЬЫ (Л = 1) к ЫЬ (Л = 0). Этот метод (таблица 2) решения задачи линейного программирования назовем параметрическим а -уровневым методом ^-продолжения задачи нечеткого линейного программирования.

Таблица 2

Вид целевой функции

Вид ограничений

X(cTxfa + (1 - X) • (cTx)Ra min(max)

X(Ax)La+(\-X).(Ax)Ra< XbRa+{\-X)-bLa

X(cTxfa + (1 - X) ■ (cTx)La min(max)

a

С помощью выбора X можно оценить уровень влияния факторов отношения LR — RL. Ясно, что, задавая промежуточные значения X в промежутке от нуля до единицы, получим возможность влиять на характеристики задачи: разрешимость, устойчивость, чувствительность к изменению параметров.

3. Пример

Сравним метод четырех задач [3] и предложенный метод на простом примере

/(х) = CjXj + с2х2 +1 —> max ,

Пусть Cj, с2, а}. а2, h] . b2 — треугольные нечеткие числа:

с1=(2, 3, 4), с2=( 1, 2, 3), Oj=(l, 2, 3), а2=(2, 3, 4), 6, =(11, 12, 13), 62=(23,

24, 25). Имеем четыре задачи согласно таблице 1. Все четыре задачи решены в разработанной программе в среде МаШСАБ для четырех угловых случаев ЬЫ, ЫЬ (рис. 1., табл. 1), для дискретных значений а = 0; 0,1; 0,4; 0,6; 0,8; 1,0 . Выпишем все четыре задачи.

(3)

1) fL(x) = (2 + a)Xj + (1 + а)х2 +1 —»max

(4)

2) fL(x) = (2 + a)Xj + (l + a)x2 +1 —»max

(5)

(3-a)Xj <(11 + a), (4 -a)x2 < (23 + a).

3) fR (х) = (4 - а)хх + (3 - а)х2 +1 —> шах

(6)

4) fR (х) = (4 - а)хх + (3 - а)х2 +1 —> шах

(7)

(3 -а)хх < (11 +а), (4 -а)х2 < (23 + а).

Например, при а = 0,9 задача (4) имеет оптимальное значение /ь = 35,258 (х1 = 6,368, х2 = 8,31). При а = 0,9 задача (5) имеет оптимальное значение /^=32,082 (х1 =5,667, х2 =7,71). На рис. 4 показаны оптимальные значения /ь (х) для задач (4) и (5) при изменении а и X -1.

Аналогично для задач (6)-(7) построим соответствующие кривые (рис. 5), показывающие изменения оптимальных значений задач (6)-(7) от изменения уровня а , (Я = 1).

0,2 0,4 0,6 0,8 1

_(4) ¿J

Рис. 4. Изменение оптимальных значений (верхняя кривая — для задачи (4), нижняя — для (5))

Зэ

а

а

(6) Р)

Рис. 5. Изменение оптимальных значений

(верхняя кривая — для задачи (6), нижняя — для (7))

Из рис. 4-5 видно, что при а = 1 (все коэффициенты становятся «четкими») значение целевой функции у всех задач одинаково f(x) = 35 . Также видно, что в задачах (4) и (7), когда для «слабого» отношения целевая функция строится по левосторонней формуле (задача (4)) и для «сильного» отношения — по правосторонней формуле (задача (7)), изменение целевой функции в меньшей степени зависит от а, чем в задачах (5) и (6).

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

Так как в работе используются только треугольные числа, то результатом сложения левых частей ограничений задачи (3) являются также треугольные числа [5]:

Тогда левые части ограничений задачи (3) имеют вид:

g = an(a)x1+... + ai„x„,ieM, L,=(gf(x),gf(x),gf(x)),

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

¡uL (х) = п

-}}, te R, i еМ .

с соответствующими функциями принадлежности:

t — х Зх —t (х) = max{0,min{-L,—!-}}, feR,

/ \ fг\ . ./ ^х^ л I"*

[iL (х) = max{0,min{--,—=-}}, feR.

При а = 0,9, (Я = 1 ) задача (4) имеет оптимальные значения Xj = 6,368, х2 = 8,31, функция принадлежности имеет вид:

гп • /"6,368 19,104 — i ц(x) = max{0,min{ ' , ' }},teR, (8)

6,368 6,368

tn . J-16,62 33,24-1 /iL;(x) = max{0,min{ , }}, te R.

o,31 O, 31

Задача (5) при a = 0.9 < Л = I ) имеет оптимальные значения Xj =5,667, х2 = 7,71, и функция принадлежности для левой части будет иметь вид:

fn • J-5,667 17,001 — i /i (x) = max{0,min{ ' , ' }}, te R, (9)

5,667 5,667

in • Д~ 15,42 30,84-i

u, (x) = тахШ,т1щ-,-ieR.

7,71 7,71

Строится функция принадлежности нечеткого отношения /; . На рис.

6 для примера показаны функции принадлежности левой и правой частей первого ограничения задачи (3) для случаев (8) (рис. 6а) и (9) (рис. 66) соответственно. На рис. 6а показан вид функций принадлежности задачи (3) при а = 0,9, Л = 1 для случая (4) с функциями принадлежности (10). На рис. 66 — для случая (5) с функциями принадлежности (11).

Рис. 6а. Вид функций принадлежности задачи (3) для случая (4)

/

/| '1 \

Л

Mi, Mb,

Рис. 66. Вид функций принадлежности задачи (3) для случая (5)

mili шах

Функции принадлежности /; нечетких отношений < и < задаются формулами:

ц^ (Ll,bl) = (sup{min{^ (t),iuh (z)} \ t<z}, t,ze R), (10)

/"min (4 A) = (supímin^ (t),¡ub¡ (z)} \t<z},t,zeR).

цт (A A) = (inf{max{l - nh {t) ,1 - ^ (z)} \t < z}, t,z eR), (11)

(L2,b2) = (inf{max{l - ^ (t),l - ^ (z)} \t < z}, t,z eR).

где j.ih , i = 1,2 — функции принадлежности правых частей ограничений. На рис. 7 для примера показаны функции принадлежности /;

нечеткого отношения неравенства первого ограничения задачи (3) для случаев (10) и (11) соответственно. На рис. 7а показан вид функций принадлежности /; нечеткого отношения задачи (3) для первого ограничения при а = 0,9, Л = 1 для случая (4) с функцией принадлежности (10), а на рис. 76 — для случая (5) с функциями принадлежности (11).

Л «

Рис. 7а. ¡Ll~ нечеткого отношения задачи (3) для случая (4)

Рис. 76. ¡Ll~ нечеткого отношения задачи (3) для случая (5) На рис. 7а sup/i (Ц,Ь{) «0,83 , а на рис. 76 inf ju {Ц,Ь{) «0,85 . Это

mili

показывает, что < можно трактовать как более «слабое» требование

шах

ограничения по сравнению с < . Допустимое множество решений X для задачи (3) имеет функцию принадлежности вида [3]:

цх(х) = min((Ц,^),/и£min (L2,b2),tи£Ш,Х (Ц,6,),¡u^ (L2,b2)) . Так как целевая функция / — это треугольное число вида / = (2Xj + х2 +1; 3Xj + 2х2 +1; 4x¡ + Зх2 +1) с функцией принадлежности вида:

. J-2x,-xn-l 4x+3x0+l-f „ ¡лг =max{0,min{-!-=-,—!-=-}}, te R,

I эс 2 ^С^ I »Х- 2

то удовлетворяющее (оптимальное) решение X * имеет функцию принадлежности следующего вида:

¡их,(х) = min((Ц,6,),¡u^ (L2,b2),(АЛ).^ (4 АXHf) • Для задач (4) и (5) /;; имеет вид соответственно:

f. . í - 22,046 51,402-f

=max{0,min{-,-}}, íeR,

f 14,678 14,678

f. . í-20,004 46,798 -1 fif = max{0,min{ , }}, íeR.

13,377 13,377

Значения /;; меняются в зависимости от найденных оптимальных xt и х2 для каждого а . При требуемом уровне эффективности а = 0,9 (Я = 1 ) inf /i (х) ~ 0,15, который можно снижать с помощью уменьшения значений а .

Используя таблицу 2, выпишем две задачи для параметрического а -уровневого метода Л -продолжения задачи нечеткого линейного программирования. Расширенные условия задачи будут иметь двойной вид:

Л((2 + а)х1 +(1 + а)х2 +1) + (1-Я)((4-а)х1 +(3-а)х2 +1) ^-тах (12) ЯД + аЦ + (1-Я)(3-а)х1 <Я(13-а) + (1-Я)(11 + а), Л(2 + а)х2 + (1 - Л)(4-а)х2 < Л(25 - а) + (1 - Я)(23 +а), где а,Яе [0,1].

Я^-аЦ +(3-а)х2 +1) + (1-Я)((2 + а)х1 + (1 + а)х2 +1) ^-тах(13) ЯД + аЦ + (1-Я)(3-а)х1 <Я(13-а) + (1-Я)(11 + а), Я(2 + а)х2 + (1 - Я)(4 -а)х2 < Л(25 - а) + (1 - Я)(23 +а), где а,Яе [ОД].

При Л = 1 задачи (12) и (13) сводятся к задачам (4) и (5) (рис. 7а). При Л = 0 задачи (12) и (13) сводятся к задачам к (6) и (7). На рис. 8 показано решение задачи (12)—(13) при Л = \ (рис. 8а) и Я = 0,3 (рис. 86). При а = 0 и Л = 1 разность значений целевой функции задач (12)и(13) равна 26 (рис. 7а), а при а = 0и Я = 0.3 эта разность составляет всего 6 единиц (рис. 76). Это показывает, что разработанный метод дает широкие возможности влиять на свойства задачи, на значения целевой функции, задавая величину Я.

Значение Я влияет и на функции принадлежности для множества допустимых решений X задачи нечеткого линейного программирования, а также для удовлетворяющего (оптимального) решения X* задачи (3). На рис. 8а показано сравнение оптимальных значений задач (12)—(13) при а = 0; ОД; 0,4; 0,6; 0,8; 1 при Я = 1, рис. 86 — Я = 0,3 (верхняя кривая — для задачи (12), нижняя — для (13)).

а 0,3 ас 1

а

(4)

Рис. 8а. Сравнение значений задач (12)—(13) при Я = 1

Рис. 86. Сравнение значений задач (12)—(13) при Л = 0,3

Так как при а = 0,9, Л = 1 задача (12) имеет оптимальные значения х1 = 6,368, х2 = 8,31, то функция принадлежности имеет вид (8), а а = 0,9, Л = 0,3 (х1 = 5,863, х2 = 7,882 ), функция принадлежности имеет вид:

г« • ,'-5,863 17,589-Г„ ^(х) = тах{0,1шп{ ' , ' }}, ГеЯ, (14)

5,ооЗ 5,ооЗ

,п • ,'-15,768 31,528 — ^ №ь2(х) = тах{0,тт{ ' 7~882 ^

4. Обсуждение результатов

Сравнив функции принадлежности левой части первого ограничения задачи (12) при а = 0,9, Л = \ (рис. 9а) (совпадает с задачей (4)) и при а = 0,9, Я = 0,3 (рис. 96), можно увидеть, что во втором случае серединное значение левой части ближе к серединному значению правой части. На рис. 9а показан вид функций принадлежности левой и правой частей первого ограничения задачи (12) при а = 0,9, Л = 1; рис. 96 — при а = 0,9, Л = 0,3 ; (внешняя кривая — функция принадлежности левой части первого ограничения задачи (12), внутренняя — функция принадлежности правой части первого ограничения задачи (12)).

Рис. 9а. Функция принадлежности задачи (12) при а = 0,9, Л = 1

5 10 1.5 20

Рис. 96. Функция принадлежности задачи (12) при а = 0,9, Л = 0,3

Аналогичным образом можно сравнить функции принадлежности удовлетворяющего (оптимального) решения X* задачи (3) для случаев (12) и (13). На рис. 10а показан вид функций принадлежности /; нечеткого отношения задачи (12) для первого ограничения при а = 0,9, Л = 1; рис. 106 — при а = 0,9, Л = 0,3 .

Рис. 10а. задачи (12) при а = 0,9, Л = 1

Рис. 106. ¡uz задачи (12) при а = 0,9, Л = 0,3

На рис. 10а sup/i (Ц,Ь{) «0,85, а на рис. 86 sup/i (Ц,Ь{) «0,98.

Приведенные результаты во втором случае, с точки зрения экспертного понимания природы этой задачи, свидетельствуют о большей степени «четкости» получаемых данных. Таким образом, при а = 0,9, Л = 0,3 степень нечеткости снижена и задача близка к четкой постановке.

5. Заключение

В данной статье рассмотрен известный [3] метод четырех задач нечеткого линейного программирования, который обобщен до параметрического а -уровневого метода Л -продолжения. Предложенный метод позволяет сократить количество рассмотренных случаев с четырех до двух, позволяет получать дополнительные решения на подмножествах квадрата а - Л . Преимущество такого подхода состоит в том, что он включает в себя четыре задачи как частные случаи. Новый подход позволяет получать более общие и гибкие решения, дает возможность влиять на свойства линейной задачи оптимизации: меру принадлежности, чувствительность, устойчивость. В статье также рассмотрен численный пример, показывающий работу нового метода и реализованный в программной среде MathCAD. Вышеизложенные исследования показывают, что нечеткий подход обобщает известную методику [1], но при этом усложняет расчеты. Полученные результаты нуждаются в интерпретации, которые возможны за счет анализа функций принадлежности и сферы применимости метода, экспертного понимания природы таких задач, например в экономической области, финансовых инвестиционных задачах.

Литература

1. Зайченко Ю. П. Исследование операций. Нечеткая оптимизация. Киев: Выща школа, 1991. 191 с.

2. Беллман Р., Заде Л. Принятие решений в расплывчатых условиях // Вопросы анализа и процедуры принятия решений. М.: Мир, 1976. С. 172— 215.

3. Задачи линейной оптимизации с неточными данными / М. Фидлер [и др.]. М.; Ижевск: Изд-во ин-та компьютерных исследований, 2008. 288 с.

4. Стародубцев И. Ю. Решение задачи линейного программирования с нечеткими параметрами // Технические науки — от теории к практике: материалы VI Междунар. науч.-практ. конф. Новосибирск: Изд-во СибАК, 2012.

5. Иегат А. Нечеткое моделирование и управление: пер. с англ. 2-е изд. М : БИНОМ, 2013.798 с.

6. Шаталова А. Ю., Лебедев К. А. Нечеткое линейное программирование в задаче оптимального финансирования инвестиционных проектов, максимизирующей получаемый предприятием доход // Международный журнал прикладных и фундаментальных исследований. 2015. Ч. 1, № 9. С.35-36.

7. Ортега Дж., Рейнболдт В. Итерационные методы решения нелинейных систем уравнений со многими переменными. М.: Мир, 1975. 559 с.

PARAMETRIC a -LEVEL METHOD OF Я -CONTINUATION FOR FUZZY LINEAR PROGRAMMING PROBLEM

Alevtina Yu. Shatalova

Research Assistant,

Kuban State University

149 Stavropolskaya St., Krasnodar 350040, Russia

E-mail: al-shatalova@yandex.ru

Konstantin A. Lebedev

Dr. Sci. (Phys. and Math.), A/Prof.,

Kuban State University

149 Stavropolskaya St., Krasnodar 350040, Russia

E-mail: klebedev.ya@yandex.ru.

The article considers the well-known method of four fuzzy linear programming problems, which is generalized to the parametric a -level method of A -continuation. The proposed method makes it possible to reduce the number of cases considered from four to two, and obtain additional solutions on subsets of the square a - Л. The advantage of the proposed approach is that one task includes the previous four tasks as special cases. This method allows us to obtain more general and flexible solutions and to influence the properties of a linear optimization problem: degree of membership, sensitivity, stability. We also give the numerical example of the new method implementation in the programming support environment MathCAD.

Keywords: fuzzy linear programming; optimization problems; parametric linear programming with fuzzy data; expansion principle; the method of continuation; fuzzy sets; fuzzy logic.

References

1. Zaichenko Yu. P. Issledovanie operatsii. Nechetkaya optimizatsiya [Operations Research. Fuzzy Optimization], Kiev: Vishcha Shkola, 1991. 191 p.

2. Bellman R., Zadeh L. A. Decision Making in a Fuzzy Environment. Management Science. 1970. No. 17B (4). Pp. 141-164.

3. Fiedler M., Nedoma J., Ramik J., Röhn J., Zimmermann K. Linear Optimization Problems with Inexact Data. Springer US, 2006.

4. Starodubtsev I. Yu. Reshenie zadachi lineinogo programmirovaniya s nechetkimi parametrami [A Solution of the Problem of Linear Programming with Fuzzy Parameters]. Tekhnicheskie nauki - ot teorii k praktike — Technical Sciences: Bedside Approach. Coll. art. VI Int. sci.-pract. conf. Novosibirsk: SibAK Publ., 2012.

5. Piegat A. Fuzzy Modeling and Control. Germany: Physica-Verlag GmbH Heidelberg, 2010.

6. Shatalova A. Yu., Lebedev K. A. Nechetkoe lineinoe programmirovanie v zadache optimal'nogo finansirovaniya investitsionnykh proektov, mak-simiziruyushchei poluchaemyi predpriyatiem dokhod [Fuzzy Linear Programming in the Problem of Optimal Financing of Investment Projects, Maximizing the Enterprise's Income]. Mezhdunarodnyi zhurnal prikladnykh i fundamen-tal'nykh issledovanii —International Journal of Applied and Fundamental Research. 2015. No. 9 (Part 1). Pp. 35-36.

7. Ortega J. M., Rheinboldt W. Iterative Solution of Nonlinear Equations in Several Variables. New York: Academic Press, 572 p.