ФИЗИКО-МАТЕМАТИЧЕСКИЕ НАУКИ
УДК 519.816
DOI 10.34822/1999-7604-2021-1-51-62
МЕТОД АППРОКСИМАЦИОННОЙ МАТРИЦЫ ФОРМИРОВАНИЯ ВЕСОВ ОБЪЕКТОВ В МНОГОКРИТЕРИАЛЬНЫХ ЗАДАЧАХ ВЫБОРА
В. П. Корнеенко
Институт проблем управления им. В. А. Трапезникова Российской академии наук, Москва, Россия E-mail: vkorn@ipu.ru
В статье представлен оптимизационный метод формирования количественных весов объектов (важности критериев, приоритетов альтернатив) по исходной экспертной матрице суждений в многокритериальных задачах выбора. Поскольку матрицу парных сравнений можно рассматривать как некоторое возмущение мультипликативной матрицы, то предлагаемый метод базируется на аппроксимации исходной матрицы парных сравнений мультипликативной матрицей по матричному критерию минимума расстояний между матрицами. Проведен сравнительный анализ эффективности нового метода с методом аналитической иерархии Т. Саати по критерию близости к исходной матрице суждений мультипликативных матриц, элементы которых восстановлены по найденным нормированным элементам весов объектов. На примере решения задачи формирования весов важности критериев оценена точность приближенного решения метода аналитической иерархии. Данный метод математически обоснован и в связи с простотой вычисления может быть рекомендован вместо метода анализа иерархий Т. Саати при решении прикладных многокритериальных задач.
Ключевые слова: многокритериальный выбор, нормированные веса объектов, мультипликативная матрица, матричный критерий.
APPROXIMATION MATRIX METHOD FOR WEIGHTS FORMATION OF OBJECTS IN MULTICRITERIA PROBLEMS
V. P. Korneenko
V. A. Trapeznikov Institute of Control Sciences, Russian Academy of Sciences, Moscow, Russia E-mail: vkorn@ipu.ru
The article presents an optimization method for the formation of quantitative weights of objects (importance of criteria, priorities of alternatives) according to the original expert judgment matrix in multicriteria selection problems. Since the pairwise comparison matrix can be considered as some perturbation of the multiplicative matrix, the proposed method is based on the approximation of the original pairwise comparison matrix by a multiplicative matrix according to the matrix criterion of minimum distances between matrices. A comparative analysis of the effectiveness of the new method with the well-known method of the analytical hierarchy of T. Saaty is carried out according to the criterion of closeness to the original judgment matrix of multiplicative matrices, the elements of which are reconstructed from the found normalized elements of the weights of objects. On the example of solving the problem for weights forming of the importance of criteria, the accuracy of the approximate solution of the analytical hierarchy method is estimated. This method is mathematically substantiated and, due to its computational simplicity, can be recommended instead of T. Saaty's hierarchy analysis method in solving applied multicriteria problems.
Keywords: multicriteria choice, normalized weights of objects, multiplicative matrix, matrix criterion.
Введение
При решении задач многокритериального выбора на множестве объектов (альтернатив, управленческих решений, вариантов), представленных в виде отношений предпочтений, возникает проблема их экспертного оценивания по количественной шкале отношений и аддитивной свертке критериев.
Для решения указанной проблемы предложены подходы и методы, базирующиеся на результирующих отношениях предпочтения для сужения множества недоминируемых альтернатив, а также на парных сравнениях объектов (решений, критериев) [1-3], которые не всегда позволяют выделить единственную альтернативу или управленческое решение без привлечения дополнительной информации.
Однако при решении прикладных задач, связанных с измерением объектов в экспертных шкалах, а также формированием локальных весов критериев, представленных в виде иерархического дерева, обычно используются экспертные методы оценки и ранжирования объектов [4]. Прямые методы экспертного оценивания весов критериев нашли применение в методике планирования посредством относительных показателей технической оценки -PATTERN (Planning Assistanse Throgh Technical Relevance Numbe), - в которой экспертам предлагается оценить по количественной шкале нормированные локальные веса критериев на каждом уровне иерархии для вычисления глобальных весов перемножением локальных весов по ветвям многоуровневого дерева критериев [5].
Другим экспертным подходом к назначению «весов» конечному набору сравниваемых объектов на основе матрицы парных сравнений является метод анализа иерархий (МАИ) (the Analytic Hierarchy Process) Т. Саати [6-7], который в настоящее время прочно вошел в теорию и практику задач многокритериального выбора.
В соответствии с МАИ экспертами формируется так называемая матрица парных
сравнений (суждений) объектов V = Vj , i, j = 1, n, в шкале отношений, а затем находится
правый собственный вектор w = (w^,..., этой матрицы, отвечающий максимальному
собственному значению. За искомый принимается вектор, элементы которого нормированы по сумме элементов правого собственного вектора.
Поскольку вычисление векторов весов объектов (критериев и альтернатив) выполняется численным (приближенным) методом, то, отдавая себе в этом отчет, Т. Саати ввел специальный числовой показатель consistency index - индекс совместимости матрицы суждений V и мультипликативной матрицы W, полученной на основе значений собственного вектора, в виде произведения Адамара [8]:
SJ. = \erVoWTe,er =(l,l,..., l), п
где V - исходная матрица суждений, а W =
w-
ч
wj
мультипликативная квадратная
матрица, элементы которой определяются по нормированным элементам правого собственного вектора м> матрицы суждений V.
При этом имеет место соотношение
VIЦ^к . . —
мл- =---= —, /,7 = 1, п.
Индекс 8.1. характеризует степень доверия к полученным с помощью МАИ результатам и трактуется как своеобразная мера отклонения исходной возмущенной матрицы суждений V от мультипликативной Ж.
В работе Т. Саати [7, с. 76] показано, что если выполнить матричное произведение
Amax
Адамара, то индекс совместимости примет вид S .1. =-, где A max > 1 - максимальное
n
собственное число матрицы V. При достаточно близком приближении к единичному значению индекса матрица парных сравнений V «близка» к мультипликативной матрице W. Если же индекс согласованности превышает некоторое «пороговое» значение, то сделать вывод о близости указанных матриц нельзя, поэтому и применять МАИ в таких случаях не рекомендуется.
В. Д. Ногин утверждает, что «по значению индекса совместности можно лишь опосредованно судить о величине итоговой "модельной" ошибки: точно она никогда и никем не может быть определена. Такова специфика данного эвристического подхода» [9, с. 1262].
В статье предлагается более эффективный, чем метод анализа иерархий Т. Саати, метод аппроксимационной матрицы (МАМ) формирования оптимальных весов объектов по матричному критерию расстояния к исходной матрице суждений.
Постановка задачи аппроксимации матрицы суждений
В многокритериальных прикладных задачах механизм агрегирования обычно представлен в виде аддитивной свертки оценок объектов (альтернатив, вариантов) а, l = 1, nA, по критериям с весами важности в виде [3]:
F (У, = ^=i w( fj ) х j,
где у() = fj(а) - оценка объекта аг е А = {а'11 = 1, пА} в результирующей шкале по критерию fj, ] = 1, п; ^^ = w( fj) - количественный (нормированный) вес критерия.
Рассмотрим постановку формирования весов объектов по критерию близости к исходной матрице парных сравнений мультипликативной матрицы. Пусть в качестве исходных данных имеем: V = у ] - исходная экспертная матрица суждений относительной важности
объектов (критериев) ^ е Г, /', j = 1, п; Жп — множество мультипликативных квадратных матриц п-го порядка над полем действительных чисел.
За меру близости ё(V, Ж) между исходной V = у ] и мультипликативной Ж = []
матрицей, где Ж еЖп, примем квадрат евклидовой матричной /2-нормы, равной разности этих матриц [8]:
d(v, ж)=V - we- WS- (i)
Тогда математическая постановка задачи выбора мультипликативной матрицы Ж, аппроксимирующей исходную матрицу V = [vj ], т. е. наиболее приближающейся к исходной
экспертной матрице суждений, сводится к минимизации показателя (1) в виде:
■ i/ ■ л(Уа - w и) ^ min, (2)
,=1 ^j =Л ч у; w^ (2)
при условии мультипликативности между элементами матрицы W = [wj ]:
wij = wikwkj, Vi, j k =1 n. (3)
В силу условия (3) аналитически не удается получить решение исходной задачи одним из классических методов оптимизации, например методом множителей Лагранжа.
Воспользуемся свойствами мультипликативной матрицы и решим данную проблему сведением исходной задачи к эквивалентной.
Особые линейные свойства мультипликативной матрицы
Для решения исходной задачи (2)-(3) необходимо найти элементы мультипликативной W = [Wy ], Vi, j = 1, n. матрицы, обеспечивающие минимум квадратичному критерию
d(V, W) (1) при условии (3). Оказывается, для мультипликативной матрицы справедлива взаимосвязь между элементами столбцов матрицы и элементами правого собственного вектора. Для этого рассмотрим два утверждения.
В работе Б. Г. Миркина [10, с. 183-184] доказано, что «матрица B = \bj], Vi, j = 1, n
является сверхтранзитивной», если существует положительный вектор х = (хъ ..., xN), такой,
г. X - _t 1
что = —, и при этом вектор х является точкой равновестия процесса: q =— Bq ,
xk Л
t = 1,2,..., который в пределе приводит к собственному вектору q = lim q . Покажем, что
t
мультипликативная матрица обладает еще рядом свойств, которые окажутся полезными при сведении исходной, аналитически не решаемой задачи (2)-(3) в рамках классических методов оптимизации к эквивалентной.
Теорема 1 (О взаимосвязи между нормированными элементами столбцов мультипликативной матрицы и элементами правого собственного вектора).
1. Между элементами Wj мультипликативной матрицы W = \Wj ], Vi, j = 1, n и любой
Т
парой компонент (wh Wj) правого собственного вектора w = (wl7wn) справедливо биек-
w■ w-
тивное отображение —\—> wt ■ , ставящее каждому отношению — в однозначное соответ-
wj wj
w-
ствие элемент w^ матрицы Wu обратно wif i—> —, при этом справедливо равенство:
Wj
W- —
wj =—, Vh j =1 n. (4)
j wj
2. Нормированные элементы Wy столбцов Wj=(w^j,...,wnj)T,j=\,n, совпадают между собой и равны нормированному правому собственному вектору, т. е.:
wj =
w
\nJ J
Wi
= w, V/ = 1, n, (5)
\wnj
М>и
где м>у =----нормированный элемент ¡-го столбца ;
№
ма =-—--нормированный элемент /-» строки правого собственного вектора м>
3. Мультипликативная матрица имеет один базисный столбец.
Доказательство. 1. Так как по условию матрица Ж = [и у ] мультипликативная, покажем, что если для любой пары чисел из {м?1,..., ип} выполняется равенство (4), то в этом слу-
чае выполняется условие мультипликативности между элементами матрицы Ж =
Ж
т. е.
Ж
существует взаимно однозначное отображение между элементами и у о- —L.
тт - Ж Жк щ Жк
Действительно, если щ. = —L и Жу = —к, то имеем = —L х к
wj ж Ж
Жк ' Ж Ж
т. е. для всех г, j, к = 1, п выполняется условие мультипликативности.
Т
2. Пусть для элементов вектора м? = (м^,..., м?п) справедливо равенство (4). Покажем,
у1кик)
Ж
=Ж
Ж жу щ
у
—■ 1
что м? = ("и>1з ...,м>п) - правый собственный вектор матрицы Ж = Убедимся, что справедливо Ж\\> = Ам>:
Ж
Ж
, щ у=1, п.
Жм> =
г Ж1 / Ж Ж1 / Ж ... Ж1 / жп Ж | ( пж ^ М>2 / "N>1 1^2 / 1^2 ... М?2 / М?п
КЖ / Ж ип / и2 ... Ж / Ж
п У
Ж,
К ип У
/7114
К пит У
Ж
= П
Ж,
Ж
К т У
= ПМ'
откуда Атак = п - максимальное собственное число.
Для произвольных столбцов и'/( и (1 < к, < п) матрицы Ж, элементы которой удовлетворяют условию мультипликативности щ (3), выполним нормировку элементов.
Так как = и ^уд, то и ^ = и д / , откуда имеем для любых номеров столбцов к, д\
Щ
Жк = ■
1к
Жд / Жд
Ж
гд
Ж
гд
=1 Жк Жк
1=
= Жд, V/ = 1, п,
Ж
=1г г д
т. е. нормированные компоненты столбцов матрицы совпадают между собой. С другой сто-
Ж
P0ны, Жгд =
г д
Ж / Жд
Ж
= ж, V/ = 1, п.
Ж
■ 1 и гд У ■ 1 Ж / ид У ■ 1Ж 1=1 гд =1 г д ¿—<1=1
Таким образом, нормированные компоненты столбцов матрицы равны также и нормированным элементам правого вектора мультипликативной матрицы.
3. Так как матрица Ж =
Ж
Ж
с помощью линейных преобразований приводится к виду
(щ/ Ж щ/ ... щ/ ип Л и2 / и2 / Ж2 ... Ж2 /
Ж / Щ Ж / Ж и / и
К'г п' 1 п' 2 п' п
■П ¡=1 ик
г1/ и1 1/ и2 ... 1/ ипЛ 1 / Ж1 1 / Ж2 ... 1 / Жп
п ' п У
1 / Ж1 1 / Ж2 ... 1 / Жп
■П 1=1 ик
г1/ и1 1/ и2 ... 1/ ипЛ 0 0... 0
пУ
0...
0
0
V
/
то ранг мультипликативной матрицы равен единице: = 1, т. е. мультипликативная матрица имеет один базисный столбец (базисную строку). Теорема доказана.
Теорема 2 (О приведении элементов столбцов матрицы к целочисленным значениям
и равенству правому собственному вектору). Пусть элементы д^, /, к = 1, п мультипликативной матрицы Ж = [д^], V/, к = 1, п в общем случае представлены целыми и рациональными числами в виде неправильных дробей.
Тогда, если элементы д1к любого столбца дк, к = \,п, матрицы умножить на наименьшее общее кратное щ знаменателей рациональных элементов столбца, то полученные в результате этой процедуры целочисленные элементы = щд^, Vk = 1, п, столбцов матрицы совпадут между собой по любому номеру строки и будут равны целочисленным элементам правого собственного вектора м> = (м>ъ ...,м>пУГ, соответствующего
собственному максимальному числу = п:
г1к
V гпк у
V жп
1 < к < п.
(6)
Доказательство. Так как элементы любой мультипликативной обратно-симметричной матрицы в соответствии с утверждением теоремы 1 можно привести к виду (4), а именно:
ж
д1к =—-, где ж - положительные целые числа - элементы правого собственного вектора,
жк
то наименьшее общее кратное натуральных чисел для знаменателей рациональных элементов такой мультипликативной матрицы будет равно пк = ж к, Vk = 1, п.
В результате чего для к-го столбца дк =
V жк
ж
\Т
ж
находим
к у
(м>А
п$к =пкх V жп / пк ) = м>кх . жп / пк ) — V жп у
= м>,\/к = 1, п,
т. е. убедились в справедливости (6). Очевидно, что м> - правый собственный вектор матрицы Ж, что и требовалось доказать. Теорема доказана.
Пример 1. Рассмотрим мультипликативную матрицу:
Ж =
Г 1 2 1/2 1 1/3 2/3
3 5 ^ 3/2 5/2 1 5 / 3
ч 1/5 2/5 3/5 1 у
(7)
Для столбцов матрицы наименьшие общие кратные составят:
п1 = ж = 2 х 3 х 5 = 30, п = ж = 3 х 5 = 15, п3 = ж = 2 х 5 = 10, п4 = 2 х 3 = 6.
Откуда придем к матрице:
Ж =
(30 30 30 301
15 15 15 15
10 10 10 10
6 6 6 6
Легко убедиться, что правый собственный вектор = (30,15,10, 6) матрицы Ж (7),
соответствующий собственному числу X = 4, и нормированные элементы столбцов матрицы совпадают между собой и равны элементам нормированного правого вектора матрицы:
=
Г 30 15 10 6 ^ 61' 61' 61' 61
Очевидно, что исходная матрица представима также и в виде Ж =
Ж
.Ж
Ж =
(30/30 30/15 30/10 30/61
15/30 15/15 15/10 15/6
10/30 10/15 10/10 10/6
6/30 6/15 6/10 6/6
Сведение исходной задачи к эквивалентной
Так как в соответствии с (5) элементы /-х строк нормированных вектор-столбцов мультипликативной матрицы равны /-ой нормированной компоненте правого собственного вектора, а именно:
м>л=--- = м>у=--- = = м>и
(8)
где = '' то
ж
Ж
имеем: =
ч
/ уп ' у к=1
Ж
ж
Убедимся, что условие мультипликативности (3) для нормированного элементов строк (8)
ж ж
матрицы выполняется: мг^мг^- = х = = м^-, V/, у, к = 1, п.
м>к
Таким образом, нахождение элементов аппроксимирующей матрицы Ж = [ и у ] задачи (2)-(3) равносильно нахождению элементов нормированного вектор-столбца:
В качестве целевого показателя задачи аппроксимации примем квадрат разности между нормированными элементами исходной и мультипликативной матрицы в виде:
т
Тогда математическая постановка исходной задачи сведется к нахождению нормированного правого вектор-столбца аппроксимирующей матрицы и обеспечивающего минимум критерия (9):
^п -г—\п ~ ч2
(10)
где Уу =
уч
- нормированные элементы матрицы V = [у, ] суждений , = 1, п.
г =1*4
Обоснованность перехода от исходной возмущенной матрицы V = [у, ] к нормированной V = \Уу ] базируется на том, что для исходной и нормированной матрицы сохраняются
ранг и величина отношения между элементами одного и того же вектор-столбца.
Теорема 3 (Об оптимальном решении). Оптимальным решением задачи (10) является
(
нормированный вектор м?* =
~ * УWnJ
, компоненты которого принимаются за коэффициенты
важности критериев м>* = и'(.//) и вычисляются по формуле:
М* = - X '=Л, =1, п
п
доставляя при этом минимум показателю (9). Элементы оптимальной аппроксимирующей матрицы Ж =
М;,
определяются по формуле: м>*- =
Я;
За оценку приближения аппроксимационной матрицы к исходной матрице суждений примем евклидову матричную норму:
а (V, ) =V - = £ 'IX и*, - мМ* )2.
(11)
Доказательство тривиально и базируется на необходимых и достаточных условиях существования оптимального решения функции многих переменных.
Алгоритм формирования оптимальных весов критериев
Шаг 1. Нормирование элементов исходной матрицы парных сравнений важности критериев по сумме элементов X П-\*ч столбцов исходной матрицы суждений:
УУ =
г=1
- нормированные элементы матрицы суждений V =
у
=У7
Шаг 2. Вычисление для каждой строки нормированной матрицы V = среднего значения, принимаемого за нормированные веса критериев, а именно:
1 п
V;,
суждений
=-
X П=1Уу, где М* = М (/) =!> п
Шаг 3. Восстановление элементов мультипликативной матрицы по оптимальным
( Л
нормированным весам критериев: Ж
МАМ
Ж
Ж'
V 1 У
Шаг 4. Оценка близости между элементами исходной матрицы парных сравнений и оптимальной мультипликативной матрицей по формуле ё(V, Ж*) (11).
Сравнение эффективности предлагаемого метода с методом анализа иерархии Т. Саати
Сравним погрешность вычисления весов важности объектов (критериев, альтернатив) методом Т. Саати с оптимизационным методом аппроксимационной матрицы формирования весов объектов в многокритериальных задачах. В качестве мультипликативных матриц рас-
смотрим квадратные матрицы W]
Г ~ Л ж
МАИ
Ж V 1
и Ж
(
ж
МАМ
г,/=1 ,п.
"Ж-
V 1 У
в которых эле-
I, У=1, га-
менты определяются по вектору приоритетов, найденных по методу анализа иерархий и методу аппроксимационной матрицы.
Для этого воспользуемся данными из примера покупки дома семьей со средними доходами, приведенными в работе Т. Саати [7, с. 41-44]. Проблема заключается в выборе дома из трех имеющихся альтернатив {А, В, С} по восьми факторам, выступающим в качестве критериев задачи многокритериального выбора.
В табл. 1 приведены экспертные оценки парного сравнения важности критериев, а также вектор приоритетов, рассчитанный по максимальному собственному числу в соответствии с ме-
8 811
тодом анализа иерархий Т. Саати ( Л.шяx = 8,811, индекс совместимости £ .I. = —-« 1,101 ).
8
Таблица 1
Исходная матрица V суждений и вектор приоритетов по МАИ
Фактор (критерии) /1 /2 /3 /4 /5 /6 /7 /8 Вектор приоритетов, п>МАИ
Размер, / 1 5 3 7 6 6 1/3 1/4 0,175
Транспорт, / 1/5 1 1/3 5 3 3 1/5 1/7 0,062
Окружение, / 1/3 3 1 6 3 4 1/2 1/5 0,103
Возраст, /4 1/7 1/5 1/6 1 1/3 1/4 1/7 1/8 0,019
Двор, /5 1/6 1/3 1/3 3 1 1/2 1/5 1/6 0,034
Удобства, / 1/6 1/3 1/4 4 2 1 1/5 1/6 0,041
Состояние, / 3 5 2 7 5 5 1 1/2 0,221
Финансы, / 4 7 5 8 6 6 2 1 0,348
Примечание: составлено по [7, с. 42].
Нормированные элементы матрицы суждений и найденный по этим данным вектор весов важности критериев по методу аппроксимационной матрицы (МАМ) представлены в табл. 2.
Нормированная матрица V
Таблица 2
1=1
-— суждений и вектор приоритетов по МАМ
и
Фактор (критерии) /1 /2 /з /4 /5 /б /7 /8 и* МАМ
Размер, Л 0,11 0,23 0,25 0,17 0,23 0,23 0,07 0,10 0,174
Транспорт, / 0,02 0,05 0,03 0,12 0,11 0,12 0,04 0,06 0,068
Окружение, / 0,04 0,14 0,08 0,15 0,11 0,16 0,11 0,08 0,108
Окончание табл. 2
Фактор (критерии) fx f2 f3 f4 f5 f6 f7 f8 w МАМ
Возраст, /4 0,02 0,01 0,01 0,02 0,01 0,01 0,03 0,05 0,021
Двор, ¡5 0,02 0,02 0,03 0,07 0,04 0,02 0,04 0,07 0,038
Удобства, /6 0,02 0,02 0,02 0,10 0,08 0,04 0,04 0,07 0,047
Состояние, / 0,33 0,23 0,17 0,17 0,19 0,19 0,22 0,20 0,212
Финансы, /8 0,44 0,32 0,41 0,20 0,23 0,23 0,44 0,39 0,333
Примечание: авторская разработка.
Для оценки точности восстановим мультипликативные матрицы парных отношений по векторам приоритетов:
^маи = (0Д75; 0,062; 0,103; 0,019; 0,034; 0,041; 0,221; 0,348), (12)
^мам = (0,174; 0,068; 0,108; 0,021; 0,038; 0,047; 0,212; 0,333), (13)
полученных по методу вычисления собственного вектора и методу аппроксимации матрицы попарных сравнений по критерию минимума расстояния, и сравним результаты.
В табл. 3-4 представлены мультипликативные матрицы ЖМЛИ, Жмлм, сформированные из нормированных весов векторов приоритетов м'мш (12) и (13).
Таблица 3
Мультипликативная матрица ^^мди' сформированная из нормированных весов вектора приоритетов и'мдд
Фактор fx f2 f3 f4 f5 f6 f7 f
Размер, / 1,00 2,82 1,70 9,21 5,15 4,27 0,79 0,50
Транспорт, /2 0,35 1,00 0,60 3,26 1,82 1,51 0,28 0,18
Окружение, /3 0,59 1,66 1,00 5,42 3,03 2,51 0,47 0,30
Возраст, /4 0,11 0,31 0,18 1,00 0,56 0,46 0,09 0,05
Двор, /5 0,19 0,55 0,33 1,79 1,00 0,83 0,15 0,10
Удобства, /6 0,23 0,66 0,40 2,16 1,21 1,00 0,19 0,12
Состояние, / 1,26 3,56 2,15 11,63 6,50 5,39 1,00 0,64
Финансы, /8 1,99 5,61 3,38 18,32 10,24 8,49 1,57 1,00
Примечание: авторская разработка.
Таблица 4
Мультипликативная матрица Щ^ЛИ' сформированная из нормированных весов вектора приоритетов и'мдд
Фактор fx f2 f3 f4 f5 f6 f7 f8
Размер, /1 1,00 2,54 1,62 8,37 4,63 3,70 0,82 0,52
Транспорт, /2 0,39 1,00 0,64 3,30 1,83 1,46 0,32 0,21
Окружение, / 0,62 1,57 1,00 5,17 2,86 2,29 0,51 0,32
Возраст, /4 0,12 0,30 0,19 1,00 0,55 0,44 0,10 0,06
Двор, /5 0,22 0,55 0,35 1,81 1,00 0,80 0,18 0,11
Удобства, /6 0,27 0,69 0,44 2,26 1,25 1,00 0,22 0,14
Состояние, / 1,22 3,10 1,98 10,21 5,65 4,52 1,00 0,64
Финансы, /8 1,91 4,86 3,10 16,02 8,86 7,08 1,57 1,00
Примечание: авторская разработка.
Находим значения нормы разности между исходными матрицами и мультипликативными матрицами отношений, сформированных из значений вектора приоритетов (важности объектов):
¿маи =|V-ИМАИ||Ч200,66 «14,2; ¿мам =\\У-ЖМАМ|| = ,[ш~6\ «11,8.
Определим точность решения ¿маи по методу Т. Саати по отношению к оптимальному ¿мам , полученному в рамках оптимизационной задачи по критерию (11). По формуле 8-приближения находим, что
£ =
_ |ИМАИ ИМАМ|
d
МАМ
14 2 -118
х 100 % = -^ х 100 % « 20,3 %.
11 8
Сравнение методов по точности получения весов объектов для исходной возмущенной матрицы 8-го порядка представлено на рисунке.
Расстоя н и е мул ьтипл и кативных матри ц до исходной матрицы суждений
16,0
14,0 12,0 10,0 8,0 6,0 4,0 2,0 0,0
14,2
11,8
Метод аппроксимационной матрицы
Метод анализа иерархий
Рисунок. Сравнение методов по точности получения весов объектов
Примечание: составлено автором.
Отсюда следует, что оценка погрешности равна 20,3 % и решение по МАИ, отличающееся от оптимального на эту величину, не может считаться удовлетворительным.
Таким образом, метод нахождения приоритетов важности критериев и объектов по собственному правому вектору матрицы парных сравнений Т. Саати следует отнести не к точным, как декларирует автор МАИ, а к приближенным.
Заключение
В работе исследована проблема формирования количественных весов объектов по матрице парных сравнений в шкале отношений. Поскольку в методе анализа иерархий Т. Саати весовой вектор определяется с помощью полиномов, то задача нахождения корней полиномов степени n > 5 неразрешима в радикалах (доказано Э. Галуа в 30-х гг. XIX в.).
Поэтому в методе Т. Саати для числа объектов не менее 5 процедура нахождения собственных значений матрицы степеней превосходства важности критериев или предпочтений альтернатив осуществляется с применением численных методов нахождения корней полинома, реализованных в пакете Expert Choice [11].
С другой стороны, из-за обратной симметричности элементов матрицы суждений
V = [Vj ], i,j = 1, n, получаемых последовательным сравнением всех пар объектов, эксперту
n (n -1)
приходится отвечать на —--- вопросов о величинах v j .
В связи с этим метод Т. Саати оправдан только для небольшого числа критериев и объектов.
На конкретном примере показано, что метод аналитической иерархии является приближенным, при этом оценена его погрешность относительно оптимального решения, полученного методом аппроксимационной матрицы формирования весов объектов в многокритериальных задачах.
Учитывая математическое обоснование, а также простоту вычислений формирования весов объектов, представленный метод может быть рекомендован вместо метода Т. Саати при решении прикладных задач.
Литература
1. Ашихмин И. В., Ройзензон Г. В. Выбор лучшего объекта на основе парных сравнений на подмножествах критериев // Методы поддержки принятия решений : сб. тр. Ин-та систем. анализа Рос. акад. наук / под ред. О. И. Ларичева. М. : Эдиториал УРСС. 2001. С. 51-71.
2. Ногин В. Д. Принятие решений в многокритериальной среде: количественный подход. М. : ФИЗМАТЛИТ, 2004. 176 с.
3. Multiple criteria decision analysis: state of the art surveys multiple criteria decision analysis: state of the art surveys / Ed. by Figueira J., Greco S., Ehrgott M. Springer, 2005. 1048 p.
4. Корнеенко В. П. Оптимизационный метод выбора результирующего ранжирования объектов, представленных в ранговой шкале измерения // Управление большими системами. Вып. 82. 2019. С. 44-60.
5. Sigford S. V., Parvin R. H. Project PATTERN a Methodology for Delernining Relevance in Complex Decision-Making // IEEE Trans. 1965. Vol. 12, No. 1. P. 9-13.
6. Saaty T. L. Axiomatic Foundation of the Analytic Hierarchy Process // Management Science. 1986. Vol. 32, Iss. 7. P. 841-855.
7. Саати Т. Принятие решений при зависимостях и обратных связей: аналитические сети. М. : Изд-во ЛКИ, 2008. 360 с.
8. Хорн Р., Джонсон Ч. Матричный анализ. М. : Мир, 1989. 655 с.
9. Ногин В. Д. Упрощенный вариант метода анализа иерархий на основе нелинейной свертки критериев // Журн. вычислит. математики и матем. физики. 2004. Т. 44, № 7. С.1259-1268.
10. Миркин Б. Г. Проблема группового выбора. М. : Наука, 1974. 256 с.
11. Expert Choice. URL: https://www.expertchoice.com/2021 (дата обращения: 05.02.2021).