Применение обобщенного логического критерия для аппроксимации области эффективности в многокритериальных задачах оптимизации Текст научной статьи по специальности «Математика»

Применение обобщенного логического критерия для аппроксимации области эффективности в многокритериальных задачах оптимизации

1 2 Д.Е. Шапошников , И.В. Костина

1 Нижегородский государственный университет им. Н.И. Лобачевского 2Нижегородский государственный технический университет им. Р.Е. Алексеева

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

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

Задачи многокритериального выбора широко распространены в различных отраслях деятельности, и, следовательно, в различных предметных областях [1-5]. При этом под задачей выбора понимается ситуация необходимости одной (наиболее подходящей) альтернативы из множества возможных (допустимых) вариантов. Большинство таких задач являются многокритериальными (задачи многокритериального выбора -МКВ), то есть имеющиеся альтернативы оцениваются сразу по нескольким критериям [6-8].

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

Существуют подходы к задаче многокритериального выбора, учитывающие качественную информацию о предпочтениях ЛПР [10]. Основным недостатком этих методов является требование к ЛПР о вводе

}J

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

Одним из наиболее распространенных способов решения задач МКВ является применение обобщенного критерия оптимальности, в состав которого входят весовые коэффициенты важности w е Dw, отражающие

представление лица, принимающего решение (ЛПР) об относительной важности частных критериев оптимальности. При этом важность критериев здесь понимается в смысле аксиоматической теории важности [11], что позволяет считать: если известна дополнительная информация вида " i -й критерий не менее важен, чем j-й критерий (Qi > Qj), то для весовых

коэффициентов wi и w j справедливо соотношение:

Q. ^ Qj « Wi > Wj . (1)

Весовые коэффициенты могут быть назначены ЛПР различными способами [10, 13, 14]. Но все известные методы требуют полных ответов (чаще всего, с числовыми оценками) на вопросы, которые выбранный метод задает.

В ряде работ [10, 12] был предложен и развит подход, при котором весовые коэффициенты важности частных критериев являются неконтролируемыми факторами, их значения могут быть различными для различных вариантов выбора, и значения их могут быть вычислены по принципу гарантированного результата:

х* = arg min max F(q(x), w). (2)

xeD weDW

В качестве обобщенного критерия оптимальности могут быть использованы различные функции [10], в частности, аддитивный и

}J

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

^шах((х^™) = шах\WjVj (х)}. (3)

1<} <п У -1

Таким образом, исходная задача (2) при использовании данного обобщенного логического критерия формулируется так:

j (x)). (4)

*

x = arg mm max max \w.q. (x

xeD weDw 1<j<n J J

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

Dw = jwwj > wo > 0, j = 1n; f4wl = rJ. (5)

В дальнейшем будем считать, что в общем случае (при отсутствии качественной информации о предпочтениях частных критериев) весовые коэффициенты wl, i = 1, n, нормированы относительно положительного параметра R > 0 и могут принимать численные значения не меньше некоторой неотрицательной величины w0 .

Можно предположить, что от ЛПР получена дополнительная качественная информация, устанавливающая для некоторых L пар частных

2

критериев (необязательно для всех Cn) предпочтение l-го критерия над j-м критерием на всем множестве D допустимых вариантов:

(6)

„ = { - О.}, l = 1,...,L < ^.

2

Информация (6) является качественной, так как из нее следует, что 1-й критерий важнее _/-го критерия, но нельзя сказать, во сколько раз. Тогда качественная информация (6) в соответствии с соотношением (1) позволяет

уточнить определение области допустимых значений (5) весовых

коэффициентов следующим образом:

^ =

wwj > > 0, } = 1,п; 2wi = Я; (wi > wJ■ j , I = 1,Ь

/=1 е'

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

В некоторых случаях ЛПР имеет возможность уточнить информацию о взаимной предпочтительности частных критериев и соответствующих весовых коэффициентах wi и Wj, связанных отношением ^ QJ^ с

помощью параметра 4 > 1:

П2

wwj > Wo > 0, J = 1,п; 2wi = Я; ^ > ,4.4 > 1,I = 1,Ь

г =1 '

. (8)

Будем считать, что дополнительная качественная информация (6) является непротиворечивой, если область ВМ! не пуста.

Для данного типа обобщенного критерия и приведенной области допустимых значений может быть получено аналитическое решение задачи вычисления весовых коэффициентов важности частных критериев оптимальности по принципу гарантированного результата [15]:

М>( х) = тах ^тах ( ^х) X (9)

где ^тах (w,Q(•х)Х = тах (x)}, (10)

1<г< N

а область допустимых значений Dw определяется выражением (8).

Предположим, что качественная информация от ЛПР представлена в виде графа О(I,О),где I - множество вершин, соответствующих частным критериям, О - множество ребер, соединяющих г -ю вершину с ■ -й тогда и только тогда, когда выполняется соотношение Qi У Qj.

Введем дополнительные обозначения:

I, i е {1,...,N} - множество вершин графа G(I,О), из которых имеется путь в вершину i, включая ее саму; П - мощность множества I;

p — произвольный путь в графе G(I,О) : p = 1П(p)}, где ^ — дуга, входящая в путь, п(p) — число дуг в пути p; - множество всех путей из вершины i в вершину к;

"' =

maxП", Pk

pep lep для всех i, k е I, (11)

1, P,k = 0,

где " - уточняющий коэффициент для соответствующей дуги l графа

G (I, П);

" = max "k . (12)

kel k

Учитывая введенные обозначения, можно сформулировать следующую теорему.

Теорема.

Оптимальным решением задачи (8)-(10) является вектор w* с компонентами, определяемыми выражениями

\r'"r)/ • Wo, i e Ir,

keIr (13)

'Wo, i e I \Ir,

N

где R' = R - w0> 0, (14)

i=1

величины "k и "" определяются из соотношений (11) и (12) соответственно, а индекс r определяется из соотношения

qr = max qk, (15)

1<k < N

*

W; =

где чк = тах

/е/.

а (*)■

(г)

ке/„

гг

+ # ■ Щ

(16)

Доказательство теоремы.

Доказательство проведем индукцией по слоям.

1. Пусть в графе предпочтений имеется только один первый слой (£ = 1). В этом бинарные отношения предпочтения на множестве частных критериев отсутствуют, и, следовательно, соотношения (13)-(16) сводятся к выражению:

Чг = тах {(К - (N - 1)щ )■ (х)} . (17)

/е/

Поскольку информация о предпочтениях отсутствует, /, = /, п/ = 1 для всех / = 1,2,..,N. Из (13) видно, что сформулированная теорема справедлива: \Я - (N -1) щ0, / = г,

*

Щ0, I ф г.

2. Пусть теорема справедлива для к слоев, то есть для любой вершины / из слоев 1,...,к имеем соотношение:

Ч/Що, Nk) ^ ЧгЩо, Nk) для всех / е -1к,

где Jk - множество вершин, расположенных на первых к слоях графа предпочтений, Nk = \Jk\ - количество вершин на первых к слоях графа, щ

вычисляются по формулам (13)-(16).

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

// = и /2 и... и / и / = ~/ и / ,

^ 1 2 т / ^

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

Оценим значение обобщенного критерия /тах(т, Q( х)):

^тах( Т Q( х)) = тах

тах ^к (х\

V

Для первого значения тах (х)

кеТу

сформулированную теорему, положив

Я = Я - т у, Т = Ту, N = Ык.

Тогда для выражения (21) можем записать:

^тахК ^х)) = тах1 та~х (x), | <

(18)

можем применить

кеТ

< тах( (Я, Т, #), т^у) = /Т)

т0 < < тах

)

1<< т..

Рассмотрим два возможных случая.

1. qr (Я, Т, N > В этом случае решением задачи (9) будет вектор

весовых коэффициентов Т , обеспечивающих максимум выражения

qr (Я, Т, N = тах

ге1„

Qг (Х)"

Г Г

Я - Ту - Ту

VV

кеТ.

4

у У / ,

Ш

ке1„

(19)

Выражение (19) линейно относительно ту и достигает максимального значения, когда Ту = т0. Тогда

/ (Ту ) = / К ) = qr (Я, Т0, Ык +1). (20)

2. qr (R, W, N) < WjQj. В этом случае весовой коэффициент для вершины j должен обеспечивать максимум выражения wjQj. Это выражение линейно относительно wj и достигает максимального значения при

W /}

Wj = W = max < 'U f. Тогда

1<,<mj { / Gk, J

f (Wj) = f (W) = q j (R, wo, Nk +1). (21)

Выражение (21) получается подстановкой значения wj = W в (16). Таким образом, учитывая выражения (20) и (21), получаем

Fmax К Q(x)) = max(qr (R, w0, Nk + 1), qj (R, ^ Nk +1)).

Если qr(R,w0,Nk +1)>qj(R,w0,Nk +1), то условие (14) выполняется

для первых k +1 слоев графа предпочтений, следовательно, соотношения (13)-(16) справедливы и для k +1 слоя.

Если qr (R, w0, Nk +1)< qj (R, w0, Nk +1), то в качестве индекса

r принимаем индекс j -той вершины. В этом случае выражение (13) имеет место, так как для всех k е I j выполняется

=.max G),

1<,<m,

Wj =(R'Gj)/ TGkJ + Gr W0, i е Ir

kelj

^ =^>0, I е 1 \ ,

а величины R определяются из соотношений (14) и (12) соответственно.

Теорема доказана.

Следствие 1.

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

А

w е Е

N

N

wi > Wo > 0, ] = 1, N; 2 wi = Я; wl > ^+1,I = 1, N -1, £ > 1

]=1

В этом случае решением задачи (9) является вектор весовых коэффициентов w* с компонентами, определяемыми выражениями:

Г(Г)/ 2 4кГ + 4'N ■ Wo, I = \г,

w = <

а; Nwo,

к=1

N

I = г +1, N,

где Я' = 1 - Wo24N > 0,

I=1

а индекс г определяется из соотношения

д, = тах дк,

1<к < N

(22)

(23)

(24)

где дк = тах

1<1<к

О, (х)'

^+ 4' * ■ Wo

2 4/

] =1

(25)

Следствие 2.

При отсутствии уточняющей информации от ЛПР, то есть в случае, когда 41 = 1, I = 1, Ь, область допустимых значений весовых коэффициентов приобретает вид:

А,, =

w е Е

N

N

Wj > w0 > 0,} = 1, N; 2 w I = Я; wl > w, I = 1, Ь

]=1

В этом случае решением задачи (2.1) является вектор весовых коэффициентов w* с компонентами, определяемыми выражениями:

w =«

Я - (N - пг) ■ w0

п

Wг

I е I,

I е I \ I,

(26)

где индекс г определяется из соотношения

*

дг = шах Чк,

1<к< N

где Чк =

'я - (N - пк) • V

пг

Ок (х)

(27)

(28)

Выражения (26)-(28) могут быть получены из выражений (13)-(16),

если положить все величины и , определяемые формулами (11) и (12)

соответственно, равными единице. Выражения для весовых коэффициентов (26)-(28) аналогичны выражениям, полученным в [10].

Следствие 3.

Для случая линейной упорядоченности частных критериев оптимальности при отсутствии дополнительной уточняющей информации = 1,1 = 1,Ь) область допустимых значений весовых коэффициентов имеет вид:

=

щ е Е

N

N

Щу > щ0 > 0, у = 1, N; ^ Щ = Я; Щ > Щ +1, г = 1, N -1

у =1

Тогда оптимальным решением задачи (9) является вектор весовых коэффициентов ч>* с компонентами

*

Щ =

Я - (N - г) • щ

0

г = 1, г,

г

(29)

"0, г = г +1, N,

где индекс г определяется из соотношения

Чг = шах Чк,

1<к < N

где Чк =

Я - (N - к) • Щ0 к

Ои (х).

(30)

(31)

Выражения (29)-(31) могут быть получены из выражений (22)-(25), если положить все величины <Е,\к, определяемые формулой (11), равными единице. Аналогичный результат приведен в [10].

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

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

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

Построим схему аппроксимации области эффективности задачи МКВ, основанную на предложенной методике вычисления весовых коэффициентов важности частных критериев. Поскольку в предложенной схеме учитывается дополнительная качественная информация, которая уточняет предпочтения ЛПР на множестве частных критериев. В качестве варьируемых параметров схемы используются уточняющие коэффициенты ^, I = 1, Ь.

Схема построения аппроксимации подобласти эффективности:

- задаются начальные значения для уточняющих коэффициентов:

41 = 1, l = й;

- выбирается шаг сетки А4, (например А4 = 0.1);

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

4/;

- для каждого набора вычисляются весовые коэффициенты важности щ , и для них находится значение обобщенного логического критерия;

- находится точка критериального пространства Од, являющаяся

оптимальной при вычисленном наборе весовых коэффициентов щ ;

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

Основные характеристики предложенной схемы:

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

- строится аппроксимация подобласти области Парето (определяется предпочтениями ЛПР);

- предпочтения ЛПР определяются по описанной выше методике;

- для аппроксимации применяется схема параметрического использования уточняющих коэффициентов 4/;

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

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

Приведем пример построения схемы аппроксимации в стандартной тестовой задаче многокритериальной оптимизации DTLZ 4.

В общем виде задача формулируется следующим образом: min f (x), min f (x),min f3 (x),

xeD xeD xeD

fl( x) = (1 + g )cos f2( x) = (1 + g)coS f3( x) = (1 + g)sin

nx

2

nxi

cos

7Ш

sin

V

2

nx2 2

(32)

та

l 2

а

где g = £(( - 0.5)2, (33)

I=3

а - задаваемый параметр степени,

0 < хг < 1, I = 1, п - область допустимых значений варьируемых параметров, п - число варьируемых пара метров.

Рассмотрим данную задачу для параметров степени а = 1 и количества варьируемых параметров п = 2.

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

Задача 1.

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

01 ^ 02, 01 ^ 03. (34)

Граф предпочтений для данного случая представлен на рис.1.

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

а

/ лд / / \

Рисунок 1. - Граф предпочтений для (34)

Таблица №1

Вычисление величин / и // для (34)

1 * , Р*

1 / =/1, / =/2 тах //2 }

2 - 1

3 - 1

Значения /1 и /2 являются варьируемыми параметрами схемы аппроксимации. Задается сетка значений с шагом А/ = 0.1. На каждой итерации в области критериев Бд выделяется и отображается эффективное

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

Аппроксимация, полученная для задачи (32)-(33) при наличии дополнительной качественной информации от ЛПР (34) представлена на рис.2.

Рис. 2. - Аппроксимация области эффективности задачи (32)-(33) при наличии отношений предпочтения (34)

Задача 2.

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

02 ^ й, 02 ^ въ- (35)

Рис. 3 - Граф предпочтений для (35)

Дополнительные коэффициенты <Е,\к и вычисляются аналогично первой задаче. Результаты вычислений приведены в таблице 2.

Таблица №2

Вычисление величин 6к и <Е,1 для (35)

1 £к , Рк 6

1 - 1

2 шах{б, 62}

3 - 1

Значения и 62 задаются в узлах сетки с шагом А^ = 0.1, аппроксимация строится аналогично первой задаче.

Аппроксимация, полученная для задачи (32)-(33) при наличии дополнительной качественной информации от ЛПР (35) представлена на рис.4.

С

0 2

О О

О

<3

Рис. 4 - Аппроксимация области эффективности задачи (32)-(33) при наличии отношений предпочтения (35)

Задача 3.

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

(36)

0 у ^ \

Рис. 5 - Граф предпочтений для (36)

Дополнительные коэффициенты Ц' и Ц вычисляются аналогично предыдущим задачам. Результаты вычислений приведены в таблице 3.

Таблица №3

Вычисление величин ^ и Ц для (36)

1 4 , Р 41

1 - 1

2 - 1

3 431 = 41, 432 = 42 ^Х^Ь 12 }

Значения 41 и 12 задаются в узлах сетки с шагом А4 = 0.1, аппроксимация строится аналогично предыдущим задачам.

Аппроксимация, полученная для задачи (32)-(33) при наличии дополнительной качественной информации от ЛПР (36) представлена на рис.6.

0.8 0 6

Рис. 6 - Аппроксимация области эффективности задачи (32)-(33) при наличии отношений предпочтения (36).

Литература

1. Батищев Д.И. Поисковые методы оптимального проектирования. - М.: Советское радио, 1975. 286 с.

2. Батищев Д.И. Методы оптимального проектирования. М.: Радио и связь, 1984. 248 с.

3. Батищев Д.И., Шапошников Д.Е. Многокритериальный выбор с учетом индивидуальных предпочтений. ИПФ РАН. Нижний Новгород, 1994. 92 с.

4. Batishchev D.I., Anuchin V.F., Shaposhnikov D.E. The Use of the Qualitative Information on the Importance of Particular Criteria for the Computation of Weighting Coefficients. // Multiobjective Problems of Mathematical Programming / Lecture Notes in Economic and Mathematical Systems, v.35. Springer-Verlag, 1991. P. 2-7

5. Касимова (Костина) И.В., Шапошников Д.Е. Использование обобщенного критерия максимального риска для вычисления весовых коэффициентов важности при решении многокритериальных задач // Материалы XVIII международной научно-технической конференции «Информационные системы и технологии ИСТ-2012». Нижний Новгород, изд-во Нижегородского государственного технического университета, 2012. С. 7879.

6. Бронштейн Е.М. Аппроксимация выпуклых множеств многогранниками // Современная математика. Фундаментальные направления. Т22, 2007. С. 5-37.

7. Гермейер Ю.Б. Введение в теорию исследования операций. М.: Наука, 1971. 384 с.

8. Демьянов В.Ф., Малоземов В.Н. Введение в минимакс. М.: Наука, 1972. 368 с.

9. Подиновский В.В., Ногин В. Д. Парето-оптимальные решения многокритериальных задач. М.: Наука, 1982. 256 с.

10. Многокритериальная оптимизация: Математические аспекты / Березовский Б.А., Барышников Ю.М., Борзенко В.И., Кемпнер Л.М. М.: Наука, 1989. 128 с.

11. Подиновский В.В. Об относительной важности критериев в многокритериальных задачах принятия решений // Многокритериальные задачи принятия решений. М.: Машиностроение, 1978. С. 48-92

12. Каменев Г.К., Лотов А.В., Рябиков А.И. Использование параллельных вычислений при аппроксимации многомерной границы Парето в задачах многокритериальной оптимизации // Доклады Пятой Международной конференции «Параллельные вычисления и задачи управления». М.: ИПУ РАН, 2010. С. 241-264.

13. Климова О.Н., Ногин В.Д. Учет взаимно зависимой информации об относительной важности критериев в процессе принятия решений // Журнал

вычислительной математики и математической физики, т. 46, №12. СПб., 2006. С. 2179-2191

14. Березкин В.Е, Каменев Г.К., Лотов А.В. Гибридные адаптивные методы аппроксимации невыпуклой многомерной границы Парето // Журнал вычислительной математики и математической физики, 2006, т. 46, №11. - С. 2009-2023

15. Лоскутов А.Б., Солнцев Е.Б., Петрицкий С.А., Терентьев П.В. Методика интегральной оценки уровня энергоэффективности непромышленных объектов // Инженерный вестник Дона, 2014, №3. URL: ivdon.ru/ru/magazine/archive/n3y2014/2477.

16. Земцов А.Н., Болгов Н.В., Божко С.Н. Многокритериальный выбор оптимальной системы управления базы данных с помощью метода анализа иерархий // Инженерный вестник Дона, 2014, №2. URL: ivdon.ru/ru/magazine/archive/n2y2014/2360

17. Лотов В.А., Поспелова И.И. Многокритериальные задачи принятия решений: учебное пособие. М: МАКС Пресс, 2008. 197 с.

18. Nelyubin A.., Podinovskiy V. Algorithmic decision rule using ordinal criteria importance coefficients with a first ordinal metric scale // Comput. Math. Math. Phys. 2012. № 1. P. 43-59.

19. Podinovskiy V. Using interval information on relative criteria tradeoffs in analyzing multicriterial decision making problems // Autom. Remote Control. 2010. Т. 71. № 8. P. 1648-1660.

20. Лазарев Е.А., Мисевич П.В., Шапошников Д.Е. Бикритериальная модель сети передачи данных // Системы Управления И Информационные Технологии. 2011. № 3.2 (45). С. 255-258.

References

1. Batishchev D.I. Poiskovye metody optimal'nogo proektirovaniya [Search methods of optimal design]. Moscow: Sovetskoe radio, 1975. 286 p.

2. Batishchev D.I. Metody optimal'nogo proektirovaniya [Methods of optimal design]. Moscow: Radio i svyaz', 1984. 248 p.

3. Batishchev D.I., Shaposhnikov D.E. Mnogokriterial'nyy vybor s uchetom individual'nykh predpochteniy [Multicriterion Choice with Taking Into Account Individual Preferences]. IPF RAN, Nizhniy Novgorod, 1994. 92 p.

4. Batishchev D.I., Anuchin V.F., Shaposhnikov D.E. Lecture Notes in Economic and Mathematical Systems, v.35. Springer-Verlag, 1991. pp. 2-7

5. Kasimova (Kostina) I.V., Shaposhnikov D.E. Materialy XVIII mezhdunarodnoy nauchno-tekhnicheskoy konferentsii «Informatsionnye sistemy i tekhnologii IST-2012». Nizhniy Novgorod, izd-vo Nizhegorodskogo gosudarstvennogo tekhnicheskogo universiteta, 2012. pp. 78-79.

6. Bronshteyn E.M. Sovremennaya Matematika. Fundamental'nye Napravleniya. V. 22, 2007. pp. 5-37.

7. Germeyer Yu.B. Vvedenie v teoriyu issledovaniya operatsiy [The Introduction of the Operations Research Theory]. Moscow: Nauka, 1971. 384 p.

8. Dem'yanov V.F., Malozemov V.N. Vvedenie v minimaks [The Introduction of the MiniMax]. Moscow: Nauka, 1972. 368 p.

9. Podinovskiy V.V., Nogin V.D. Pareto-optimal'nye resheniya mnogokriterial'nykh zadach [Pareto-optimal Solutions of the Multiobjective Problems]. Moscow: Nauka, 1982. 256 p.

10. Berezovskiy B.A., Baryshnikov Yu.M., Borzenko V.I., Kempner L.M. Mnogokriterial'naya optimizatsiya: Matematicheskie aspekty [Multicriteria Optimization: Mathematical Aspects] / Berezovskiy B.A., Baryshnikov Yu.M., Borzenko V.I., Kempner L.M. Moscow: Nauka, 1989. 128 p.

11. Podinovskiy V.V. Ob otnositel'noy vazhnosti kriteriev v mnogokriterial'nykh zadachakh prinyatiya resheniy [On the Relative Preferences of Criteria Importance

in Problems of Multicriteria Decision Making] / Mnogokriterial'nye zadachi prinyatiya resheniy. Moscow: Mashinostroenie, 1978. pp. 48-92

12. Kamenev G.K., Lotov A.V., Ryabikov A.I. Doklady Pyatoy Mezhdunarodnoy konferentsii «Parallel'nye vychisleniya i zadachi upravleniya». Moscow: IPU RAN, 2010. pp. 241-264

13. Klimova O.N., Nogin V.D. Computational Mathematics and Mathematical Physics, V. 46, №12. 2006. pp. 2179-2191

14. Berezkin V.E, Kamenev G.K., Lotov A.V. Computational Mathematics and Mathematical Physics, 2006, V. 46, №11. - pp. 2009-2023

15. Loskutov A.B., Solntsev E.B., Petritskiy S.A., Terent'ev P.V. Inzenernyj vestnik Dona (Rus), 2014, №3. URL: ivdon.ru/ru/magazine/archive/n3y2014/2477.

16. Zemtsov A.N., Bolgov N.V., Bozhko S.N. Inzenernyj vestnik Dona (Rus), 2014, №2. URL: ivdon.ru/ru/magazine/archive/n2y2014/2360

17. Lotov V.A., Pospelova I.I. Mnogokriterial'nye zadachi prinyatiya resheniy [Multicriteria tasks of Decision Making]. Uchebnoe posobie. Moscow, MAKS Press, 2008. 197 p.

18. Nelyubin A.., Podinovskiy V. Comput. Math. Math. Phys. 2012, № 1. P. 4359.

19. Podinovskiy V. Autom. Remote Control. 2010. V. 71. № 8. pp. 1648-1660.

20. Lazarev E.A., Misevich P.V., Shaposhnikov D.E. Sistemy Upravleniya I Informatsionnye Tekhnologii. 2011, № 3.2 (45). pp. 255-258.