Method of indirect preferences for forming criterion weights with a mul-ti-level structure

Cover Page

Cite item

Full Text

Abstract

The method of forming criterion weights in tasks of multi-criteria evaluation of objects with a multi-level structure of indicators presented in the form of a hierarchical tree, is shown. The weights values are calculated on the basis of indirect measurement of preferences of adjacent pairs of local weights of criteria included in the vertices of a higher level of the hierarchical tree, in the form of expert assessments in the quantitative scale of relations. The calculation of quantitative weights first comes down to lexicographic sorting by descending order of criteria importance at each level of the hierarchy that are included in the vertices of a higher level, as a result of which the number of expert comparisons of adjacent pairs in the scale of relations is reduced. The formation of local coefficients of the criteria importance is mathematically justified and is based on a matrix with special properties. A comparative analysis of the proposed method for forming quantitative weights with the method of hierarchy analysis, the least squares method and the approximation matrix method based on the matrix of paired comparisons, is given. An example is given of solving the problem of multi-criteria evaluation of combat aircraft that participated in the tender, according to the tactical and technical characteristics presented in various measurement scales.

Full Text

Введение

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

Для объектов с многоуровневой структурой характерно представление системы критериев в виде иерархического дерева [1–3]. При этом возникает проблема оценки количественной важности локальных критериев, входящих в вершину более высокого уровня.

Обзору и классификации методов формирования весовых коэффициентов посвящено большое число работ [4–9]. Известно более двадцати различных методов определения весов, часть из которых подробно рассмотрена в статье [8].

Настоящая статья посвящена развитию методов косвенного измерения предпочтений критериев в шкале отношений [10].

1. Определение важности критериев

Одна из первых методик решения многокритериальных задач с многоуровневой структурой была разработана для оценки эффективности вариантов прогнозирования и планирования при разработке сложных научных и научно-технических программ (ПАТТЕРН, PATTERN1, США) [2]. Методика включает построение дерева целей, которое служит для оценки относительной важности критериев, и применение прямого метода экспертного оценивания количественных весов критериев.

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

В соответствии с теорией измерений [11] веса важности критериев могут быть измерены в градациях различных шкал, например: в номинальных шкалах с двумя градациями «равноважны» и «неравноважны»; в количественных (разности, отношений) с различными весами важности (предпочтительности) по числу критериев.

Понятие важности критериев c точки зрения теории измерений можно рассмотреть следующим образом. Пусть объекты множества A=aq:q=1,nA¯ оцениваются по fjF=f1,,fm критериям.

Определение 1. В количественной шкале отношений критерий fi, i=1,m¯, считается более важным, чем критерий fj (обозначение  fifj), если количественный вес wi= wfi критерия fi будет превосходить количественный вес wj=wfj критерия fj в wiwj>1 раз, т.е.

fifjwfi>wfj

и наоборот, менее важен, если wiwj<1, а при wiwj=1, т.е. равенстве весов, критерии равноважны

fifjwfi=wfj.

Известно следующее определение однородных критериев по важности [12, 13].

Определение 2. Утверждение «критерий Ki важнее критерия Kj» (обозначение ij) означает, что каждая векторная оценка x в которой xi>xj, предподчительнее, чем xij.

Здесь xij – обозначение векторной оценки, которая получена из векторной x=(x1,,xm) оценки, например, если x=5, 4, 3, 4, то x23=5, 3, 4, 4.

Замечание 1. Сравнение определений 1 и 2 приводит к выводу о том, что определение 1 не привязано к векторным оценкам, в отличие от определения 2. Можно показать, что утверждение «каждая векторная оценка x в которой xi>xj, предподчительнее, чем xij» следует из предпочтительности критерия fi относительно fi и условия xi>xj. 

Действительно, пусть критерий fi важнее критерия fi, т.е. для количественных весов справедливо неравенство: wi>wj. Очевидно, векторы x и xij несравнимы по предпочтительности, однако в силу однородности шкалы измерения в задачах векторной оптимизации они сравнимы по аддитивной свёртке [5]:

Fx,w=s=1mwsxs,

где ws=wfs – количественный вес критерия fs, s=1,  m¯; x=x1,,xs,,xm.

Агрегированные оценки векторов x и xij представляются в виде:

Fx,w=wixi+wjxj+s=1,sj,simwsxs. Fxij,w=wixj+wjxi+s=1,sj,simwsxs.

Разность Fx,wFxij,w агрегированных оценок x и xij имеет вид:

Fx,wFxij,w =(wixi+wjxj)wixj+wjxi=wiwjxixj>0,

т.е. векторная оценка x, в которой xi>xj, предподчительнее, чем xij в силу предпочтительности критерия fi над критерием fj.

В работе [14] в качестве количественных весов предлагается использовать числовую последовательность П. Фишберна [15]:

wi=2m+1imm+1, (1)

где i – номер критерия в порядке убывания его важности, при этом i=1mwi=1.

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

ans=1n2  при s=1, Cs+2n1n2n1 при 1<sn, (2)

где s - группа важности, ans – универсальные коэффициенты важности в политике выбора для n критериев, в которой каждая группа важности включает ровно один критерий. Задающий параметр Cs для n=2, 3, 4,, 10 принимает значения 0, 6, 18,…, 216.

В силу того, что нормированные веса частных критериев характеризуют долю вклада в обобщённый критерий и формируются обычно методами экспертного оценивания [8], то вследствие этого формулы (1) и (2) не могут быть применимы в качестве весов критериев, связанных с различными предметными областями решаемых многокритериальных задач.

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

2. Алгоритм метода

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

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

Для деревьев в основном применяются два способа перечисления: «по ветвям», когда индекс вершины указывает путь к этой вершине; «по уровням», когда по очереди рассматриваются все уровни сверху вниз, а вершины одного уровня нумеруются подряд слева направо. Способом перечисления «по ветвям» дерево задаётся в виде множества упорядоченных вершин [18]:

ID=F,D,

где F=F0, Fj1,,Fj1jk jk= 1,nj1jk1¯;k=1,nF¯ – множество вершин (критериев), в которых индекс j1jk вершины Fj1jk указывает путь к этой вершине от корневой вершины F0 (k = 0); 0, 1, ,k,,nF – номера уровней в дереве;

 D=Fj1jk1,Fj1jk jk= 1,nj1jk1¯;k=1,nF¯ – множество дуг, в которых множество вершин {Fj1jk}, инцидентно вершине Fj1jk1;

F0 – глобальный (обобщённый) критерий верхнего (нулевого) уровня иерархии;

Fj1 – групповые критерии 1-го уровня иерархии, являющиеся концевыми вершинами множества дуг {(F0,Fj1) j1=1,n0¯ }; n0 – число дуг, инцидентных вершине F0;

Fj1jk – групповые критерии k-го уровня, являющиеся концевыми вершинами дуг Fj1jk1,Fj1jkjk= 1,nj1jk1¯ ; nj1jk1 – число дуг, инцидентных вершине Fj1jk1.

Концевые вершины nF-го нижнего уровня условно обозначены строчными буквами fj1jn. Каждому Эs эксперту, где s=1, nЭ¯, предъявляется множество подкритериев

FξFj1jk1ξ,

на k-м уровне иерархии упорядоченных по убыванию важности и входящих в вершину Fjk1 Fj1jk1:

F1F2Fξ Fnjk1, ξ=1,njk1¯, njk1nj1jk1, (3)

где ≽ – обозначение нестрогого предпочтения, которое означает для пары критериев строгого предпочтения (обозначение ≻) или отношение равнозначности (обозначение ).

Fj1jk11,,Fj1jk1nj1jk1 – множество подкритериев (k-1)-го уровня, входящих в вершину Fj1jk k-го уровня; jk= 1,nj1jk1¯;k=1,nF¯. Подкритерии, входящие в вершины более высокого уровня, упорядочиваются по важности:

Fj1jk11Fj1jk1nj1jk1  jk= 1,nj1jk1¯;k=1,nF¯.

Шаг 2. Экспертное оценивание степеней превосходства смежных критериев на k-м уровне иерархического дерева в шкале отношений. Эксперту предъявляется множество подкритериев FξFj1jk1ξ , на k-м уровне иерархии упорядоченных по убыванию важности и входящих в вершину (см. выражение (3)).

Результаты сравнения wξ,ξ+1wlξwlξ+1 смежных пар подкритериев Fξ,Fξ+1 на k-м уровне иерархии, входящих в вершину Fjk1Fj1jk1, локальные веса wlξ=wlFξ которых заранее неизвестны, представляются в виде (njk11) наддиагональных элементов

w12*, w23*,,wξ, ξ+1*,,wnjk11,njk1* (4)

квадратной матрицы попарного сравнения смежных подкритериев Fξ,Fξ+1 WFjk1wξ, ξ+1 :

WFjk1=11w1,2*0w1,2*100  w2,3*101    1wnjk11,njk*00wnjk11,njk*1. (5)

При этом количественная оценка  wξ,ξ+1 означает, что эксперт считает критерий Fξ важнее критерия Fξ+1FξFξ+1 в wξ,ξ+1wlξ/wlξ+1 раз, если wξ,ξ+1>1, а при wξ,ξ+1=1 критерии равноважны (FξFξ+1).

Исходя из свойства обратной симметрии, равенство wξ+1,ξ=1 wξ,ξ+1<1, означает, что количественно критерий Fξ+1 в wξ+1,ξ раз менее важен критерия Fξ.

При таком подходе достаточно всего лишь (njk11)-го экспертного сравнения вместо njk12njk12 при попарном сравнении как, например, в методе анализа иерархий (МАИ) Т. Саати [1].

Шаг 3. Формирование элементов столбцов матрицы коэффициентов важности критериев. По экспертным оценкам (4) предпочтительности смежных критериев, входящих в вершину критерия более высокого уровня иерархического дерева, формируются элементы столбцов матрицы WFjk1 (5):

 wξ=w1ξ,w2ξ,,wηξ,,wnjk1ξ    ξ=1,njk1¯. (6)

Любой элемент произвольного столбца матрицы можно вычислить по экспертным оценкам относительно критерия, принимаемого за базовый, которому соответствует единичный элемент wξξ =1  на главной диагонали матрицы.

Теорема 1. Пусть даны экспертные оценки wξ, ξ+1*,  ξ=1,njk11¯, предпочтительности смежных подкритериев Fξ на k-м уровне, входящие в вершину Fjk-1 на k-1-м уровне иерархии. Тогда любой элемент столбца матрицы WFjk1 (4) может быть вычислен по элементам оценок, расположенным над главной диагональю матрицы по формуле:

wηξ =t=ηξ1wt, t+1*, η<v,  η=1,ξ1¯,                1,                  η=ξ,                           t=ξη1wt, t+1*1, η>ξ,  η=ξ+1,njk1¯. (7)

Доказательство. Действительно, для номеров строк η<ξ можно записать:

wηξ =t=ηξ1wt, t+1*=wlηwlη+1×wlη+1wlη+2××wlξ2wlξ1×wlξ1wlξ=wlηwlξ ,  

аналогично для номеров строк η>ξ можно записать:

wηξ =t=ξη1wt, t+1*1=wlξwlξ+1×wlξ+1wlξ+2××wlη2wlη1×wlη1wlη1=wlξwlη1=wlηwlξ,

что и требовалось доказать.

Пример 1. Пусть для четырёх упорядоченных по важности f1f2f3f4 критериев известны количественные веса wlf1=4, wlf2=3, wlf3=2, wlf4=1.

Их нормированные веса имеют значения 0,4; 0,3; 0,2 и 0,1. Пусть локальные веса заранее неизвестны, но эксперт результаты сравнения смежных пар представил в виде:

w1243; w2332 ; w3421,

которые приняты за наддиагональные элементы матрицы. С учётом обратной симметрии матрицу можно представить в виде

W=14/3*003/4*13/2*002/3*12*001/2*1.

Остальные элементы столбцов wj=(w1j,   w2j, w3j,, w4j)T, j=1, 4¯, матрицы вычисляются по формулам, в которые обязательно входят наддиагональные элементы:

w1: w31=w32w21=23*34*=12;

w41=w43w31=12*12=14;

w2: w42=w43w32=12*23*=13;

w3: w31=w23w12=32*43*=2;

w4: w14=w12*w23*w34*=43*32*2*=4;

w24=w23*w34*=32*2*=3.

В результате получается матрица:

 (8)

Замечание 2. Можно убедиться, что ранг матрицы W (8) равен 1, и она является мультипликативной, например, для элемента w18 матрицы (8) справедливо w14=w13w34=w12w24.

Матрица WFjk1 (4) с вычисленными элементами по формуле wηξ (7) называется полной матрицей.

Шаг 4. Формирование локальных коэффициентов важности критериев. За локальные коэффициенты принимаются нормированные элементы любого столбца матрицы WFjk1 вычисленные по формуле (7). Можно доказать, что нормированные элементы полной матрицы совпадают между собой и равны нормированным локальным весам критериев.

Теорема 2 (о взаимосвязи между нормированными элементами полной матрицы и нормированными весами критериев). Нормированные элементы столбцов wξ 6,ξ=1,njk1¯, вычисленные по наддиагональным элементам по формуле wηξ (6) матрицы   WFjk1 (4), совпадают между собой и равны нормированным весам критериев, т.е.

w~ξ=w~1ξ w~njk1ξ =wl~1wl~njk1  ξ=1,njk1¯,

где w~ηξ=wηξη-1nwηξ – нормированный элемент ξ-го столбца wξ;

wl~η=wlηη=1njk1wlη – нормированный вес Fη локального критерия.

Доказательство. Из формулы (6) следует, что для полной матрицы для любого столбца wξ 5,ξ=1,njk1¯ справедливо wηξ =wlηwlξ, η=1,njk1¯. Откуда

w~ηξ=wlη/wlξη=1nwlη/wlξ=wlηη=1nwlη=wl~η η, ξ=1,njk1¯,

что и требовалось доказать.

Если в примере 1 пронормировать столбцы матрицы W (8), то они совпадут между собой и будут равны нормированным весам критериев:

W=0,40,40,40,40,30,30,30,30,20,20,20,20,10,10,10,1.

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

При решении прикладных задач проще использовать элементы последнего столбца матрицы, которые можно вычислить через произведения wξ, ξ+1*,  ξ= 1,njk11¯ экспертных оценок по формуле:

wη,njk1 =t=ηnjk11wt, t+1*, η<njk1,            1,          η=njk1,                      

а затем нормировать и принять за коэффициенты важности критериев.

Из примера 1 видно, что по формуле wη,njk1  легко вычисляются элементы столбца w4.

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

Пусть веса wlj1j2jk нормированы, т.е. их сумма по всем подвершинам произвольной вершины Fj1jk1 равна единице:

jk=1nj1jk1wlFj1j2jk=1  k=1,n¯.

Произведение весов вершин, лежащих на пути от корневой вершины F0 к произвольной концевой вершине fj1j2jn, будет представлять собой интегральный (глобальный) вес wgfj1j2jn данной концевой вершины:

null

где jn= 1,nj1jn1¯; jk= 1,nj1jk1¯;,j1= 1,n0¯. 

Тогда постановку многокритериальной задачи оценки объектов с многоуровневой структурой можно представить в виде аддитивной свёртки:

rΣl=j1=1n0j2=1nj1jn=1nj1j2jnwgj1j2jnrj1j2jnl, (9)

где rj1j2jnl=fj1j2jnal – оценка alA объекта в концевой вершине fj1j2jn дерева ID в результирующей однородной шкале.

Решение задачи можно представить в виде упорядочения объектов в соответствии с значениями агрегированных оценок:

aq1aqnArΣq1rΣqnA.

В качестве примера на рисунке 1 представлено трёхуровневое иерархическое дерево упорядоченных по важности критериев, где веса обозначены в виде:

wlj1=wlFj1,  j1=1,m¯;  wlj1j2= wlfj1j2,  j2=1,nj1¯.

 

Рисунок 1 – Трёхуровневое иерархическое дерево критериев оценки объектов

 

Замечание 3. Глобальные коэффициен-ты, вычисляемые на шаге 5 данного метода, аналогично вычисляются и в методе ПАТТЕРН [2] и в МАИ [1]. При этом глобальный вес концевого критерия в иерархическом дереве характеризует (относительную) долю, которую он вносит в глобальный критерий верхнего (нулевого) уровня иерархии.

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

Для сравнения методов используются данные из работы [19], в которой оцениваются сравнительные расстояния между городами Каир, Токио, Чикаго, Сан-Франциско, Лондон, Монреаль. Матрица парных сравнений получена методом экспертных оценок [19]:

A=11/383373193391/81/911/61/521/31/3611/361/31/353161/71/91/21/61/61.

Упорядочив города в виде Токио, Каир, Лондон, Сан-Франциско, Чикаго, Монреаль по убыванию расстояния по методу данной работы, можно сравнить расстояния между смежными парами городов токио-Каир, Каир-Лондон, Лондон-Сан-Франциско, Сан-Франциско-Чикаго, Чикаго-Монреаль, результаты которого можно представить наддиагональными элементами мультипликативной матрицы и вычисленными элементами правого столбца:

WМКП=11,200010,02011,6008,350011,505,2200012,93,48000011,2000001,

нормированные элементы которого принять за количественные веса расстояний и представить в виде вектора: wМКП=0,342;0,285;0,178;0,119;0,041;0,034.

Поскольку по нормированным элементам весов критериев восстанавливаются все элементы wij=wiwj мультипликативной матрицы, то за критерий сравнения принимается матричный критерий близости восстановленной мультипликативной матрицы W к исходной матрице A в виде матричной l1-норма [20]:

dA,W=AW1 =i=1nj=1naijwiwj. (10)

Мультипликативная матрица обладает особым свойством – любой элемент можно представить через произведение пары других

wij=wikwkj, i,j,k=1,n¯, (11)

где wji=1wij; wii=1, i,j=1,n¯.

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

Теорема 3. Между любыми элементами wij мультипликативной матрицы W=wij  i,j=1,n¯ и любой парой компонент wi, wj вектора важности критериев w= w1,,wj,wnT, j= 1,n¯ справедливо биективное отображение wiwjwij, ставящее каждому отношению wiwj в однозначное соответствие элемент wij матрицы W, и обратно wijwiwj, при этом справедливо равенство:

wij=wiwj, i,j=1,n¯. (12)

Доказательство. Так как матрица W=wij  мультипликативная, то если для любой пары чисел (wi,  wj) из w1,,wn, представленных в виде отношения wiwj, wj>0,  выполняется условие мультипликативности (11), то между любыми элементами матрицы W=wij,   wij=wiwj, существует взаимно однозначное отображение между элементами wij и wiwj. Действительно, если wik=wiwk и wkj=wkwj, то согласно (12): wikwkj=wiwk×wkwj=wiwj=wij, т.е. для всех i,j,k=1,n¯ выполняется условие мультипликативности (11), что и требовалось доказать.

Представленный МКП нахождения весов по парным смежным (n-1) сравнениям можно сравнить с МАИ [1], методом наименьших квадратов (МНК) [21] и методом аппроксимационной матрицы (МАМ) [22].

По МАИ решением служит нормированный собственной вектор матрицы , соответствующий максимальному собственному значению. По МНК в качестве искомого вектора принимается решение оптимизационной задачи [21]:

i=1nj=1naijwjwi2min

при ограничениях i=1nwi=1 и wi>0 для всех i=1,n¯.

Эти методы трудоёмки, и для построения матрицы A необходимо проведение n2n/2 попарных сравнений. В МАИ для числа объектов не менее пяти процедура нахождения собственных значений матрицы степеней превосходства важностей критериев или предпочтений альтернатив осуществляется с применением численных методов нахождения корней полинома, реализованных в [23]. По МАМ в качестве искомого вектора w принимается решение оптимизационной задачи [22]:

i=1nj=1na~ijw~i2minw~1,,w~n, (13)

где a~ij=aiji=1naij – нормированные элементы матрицы суждений A=aij,

 – нормированный элемент -го столбца мультипликативной матрицы, совпадающий с любым элементом нормированного столбца (см. теорему 2).

Исходные данные для сравнения методов представлены в таблице 1.

 

Таблица 1 – Нормализованные относительные расстояния

Город

1

wИР

wМАИ

wМНК

wМАМ

wМКП

1

2

3

4

5

6

токио

0,361

0,397

0,459

0,389

0,342

Каир

0,278

0,263

0,221

0,255

0,285

Лондон

0,177

0,164

0,141

0,167

0,178

Сан-франциско

0,132

0,116

0,107

0,126

0,119

Чикаго

0,032

0,033

0,036

0,036

0,041

монреаль

0,019

0,027

0,036

0,028

0,034

 

В столбцах 2 – 6 таблицы 1 представлены нормализованные относительные расстояния, где в столбце 2 – истинные нормализованные относительные расстояния, wИР; 3 – расстояния по МАИ, wМАИ; 4 – расстояния по МНК, wМНК; 5 – расстояния по МАМ, wМАМ; 6 – расстояния по МКП, wМКП. Данные в столбцах 2, 3 и 4 взяты из работы [19, с.188].

Результаты сравнений по критерию (10) близости к исходной матрице парных сравнений представлены в таблице 2.

 

Таблица 2 – Результаты сравнений по критерию

Метод

ИР

МАИ

МНК

МАМ

МКП

dA, W

36,41

24,62

29,22

23,10

20,26

 

Эффективность методов можно оценить по критерию близости между векторами:

dwi,wj=k=1nwikwjk.

Результаты сравнения представлены в таблице 3.

 

Таблица 3 – Результаты сравнений по критерию dwi,wj 

Метод

МАИ

МНК

МАМ

МКП

dwi,wj

0,089

0,237

0,080

0,064

 

Таким образом, методы, базирующиеся на матрице парных сравнений, можно ранжировать по убыванию эффективности: МКПМАММАИМНК.

Замечание 4. Преимущество МКП относительно методов, базирующихся на матрице парных сравнений, в том, что этап построения иерархического дерева упорядоченных критериев по убыванию важности позволяет снизить число экспертных смежных пар сравнений критериев до n-1 по сравнению с числом при парном сравнении n2n/2.

4. Пример решения задачи многокритериального оценивания

Эффективность МКП можно рассмотреть на примере решения многокритериальной задачи сравнения боевых самолётов по тактико-техническим характеристикам (ТТХ), принимавших участие в индийском тендере MMRCA [24]. Исходные данные по девяти показателям в виде множества F={f1,,f9} и оценки xjl=fjCl шести самолётов в виде множества C={C1,,C6} представлены в таблице 4, где приняты обозначения самолётов: C1 – Dassault Rafale; C2 – Eurofighter Typhoon; C3 – F-16IN Super Viper; C4 – F/A-18E/F Super Hornet; C5 – JAS 39 NG (IN); C6 – МиГ-35. Эталонные самолёты обозначены: C* – наихудший по показателям ТТХ и стоимости (эталон – худший вариант); C* – наилучший по показателям ТТХ и стоимости (эталон – лучший вариант).

 

Таблица 4 – ТТХ самолётов, участвовавших в тендере

Показатели ТТХ

Оценки самолётов в

исходных шкалах,

Эталоны

C1

C2

C3

C4

C5

C6

C*

C*

f1 - боевая нагрузка, т

9,5

7,5

7,8

8,05

5,3

7,0

5,3

9,5

f2 - управляемый вектор тяги

нет

нет

нет

нет

нет

есть

нет

есть

f3 - скороподъёмность, м/с

305

315

254

228

255

330

228

330

f4 - максимальная взлётная масса, т

24,5

23,5

21,8

29,9

14,3

23,5

14,3

29,9

f5- максимальное число Маха на высоте

1,80

2,25

2,00

1,80

2,00

2,25

1,80

2,25

f6 - практический потолок, км

15,24

19,81

18,00

15,00

15,24

17,50

15,00

19,81

f7 - стоимость, млн. $ (2011 г.)

124

120

50

55

48

45

124

45

f8 - тяговооружённость

1,03

1,18

1,10

0,93

1,18

1,10

0,93

1,18

f9 - масса топлива, т

4,70

5,00

3,37

6,78

3,36

4,80

3,36

6,78

 

В соответствии с шагом 1 алгоритма представленные в таблице 4 показатели используются как критерии. Пусть эксперты разбили критерии на три группы по предпочтительности: F1 – целевая эффективность; F2 – техническая эффективность и F3 – технико-экономическая эффективность:

F0: F1F2F3.

В каждой группе подкритерии также упорядочены по убыванию важности:

F1: f11f12f13; F2: f21f22f23; F3: f31f32f33.

На рисунке 2 представлено трёхуровневое иерархическое дерево.

 

Рисунок 2 – Иерархическое дерево критериев оценки эффективности самолётов

 

В соответствии с шагами 2 и 3 представлены матрицы попарных смежных сравнений и вычисленные правые столбцы:

F0: W0=15/35/2013/2001,

F1:W1=14/34/3011001,

F2:W2=18/78/6017/6001,

F3:W3=19/89/8011001.

В соответствии с шагом 4 сформированы локальные коэффициенты важности критериев:

F0: wl1=0,5; wl2=0,3; wl3=0,2;

F1: wl11=0,4; wl12=0,3; wl13=0,3;

F2: wl21=0,38; wl22=0,33; wl23=0,29.

F3: wl31=0,36; wl32=0,32; wl33=0,32.

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

F1: wg11=0,5×0,40=0,20; wg12=0,5×0,3=0,15; wg13=0,5×0,3=0,15;

F2: wg21=0,3×0,38=0,114; wg22=0,3×0,33=0,099; wg23=0,3×0,29=0,087.

F3: wg31=0,2×0,36=0,072; wg32=0,2×0,32=0,064; wg33=0,2×0,32=0,064.

Поскольку исходные данные представлены в различных шкалах измерения, то для применения аддитивного механизма агрегирования необходимо перейти к результирующей однородной шкале [19]. Правило перехода πi от исходных оценок xil самолётов в количественной шкале по fi критерию к ril оценкам в результирующей балльной шкале можно представить в виде множественно-точечного отображения:

πifi*+r1hi,fi*+rhir,

где him=fi*fi* m – шаг дискретизации шкалы fi критерия, fi* - минимальное значение, fi* – максимальное значение критерия;

Xi, rx r1x r=fi*+r1hi,fi*+rhi – отрезки разбиения исходной шкалы;

rR=1, 2,..., b – градации результирующей шкалы;

b – число шкальных градаций в порядковой (балльной) шкале

При этом каждой оценке, попадающей в класс Xi,r разбиения ставится в соответствие оценка ril в результирующей шкале по правилу:  xilπiril. В случае попадания оценки объекта на смежные классы ей присваивается связанный ранг, равный среднему арифметическому значению смежных рангов [25]:  xilπir¯il=r+r+12,  если  xilXi,rXi,r+1.

В качестве механизма агрегирования оценок самолётов, представленных в результирующей шкале, применён интегральный метод с учётом весов и без учёта в виде аддитивной свёртки оценок aq, q=1÷6, по концевым критериям с весами важности в виде [26]:

Fal;wg1,,wgm=j=1mwgfiril,  (14)

где ril=fial – оценка самолёта alA=all=1,nA¯, nA=6 в результирующей шкале по критерию fi,i=1,m¯, m=9; wgi=wgfi – глобальный вес fi критерия.

В таблице 5 представлены ТТХ самолётов в 100-балльной результирующей шкале, где rl=r1l,r2l,,ril,,r9l – профиль Cl самолёта, l=1÷9, а также результат агрегирования без учёта весов критериев (в последней строке таблицы).

 

Таблица 5 - Оценки самолётов в однородной 100-балльной шкале без учёта весов критериев

Показатели

Оценки самолётов в однородной 100-балльной шкале, ril=πixil

Эталоны

r1

r2

r3

r4

r5

r6

r*

r*

f1f11

100

53

60

66

1

41

1

100

f2f12

1

1

1

1

1

100

1

100

f3f13

76

86

26

1

27

100

1

100

f4f21

66

59

49

100

1

59

1

100

f5f22

1

100

45

1

45

100

1

100

f6f23

5

100

63

1

5

52

1

100

f7f31

1

6

94

88

97

100

1

100

f8f32

40

100

68,5

1

100

68,5

1

100

f9f33

40

48

1,0

100

1

43

1

100

rΣl=i=19ril

330,0

553,0

407,5

359,0

278,0

663,5

9

900

 

Ранжирование самолётов по обобщённым оценкам без учёта весов критериев можно представить в виде:

МиГ35EurofighterTyphoonF16 INSuper ViperF/A 18E/FSuper HornetDassaultRafale

Результаты агрегирования с учётом глобальных весов по формуле (14) представлены в таблице 6, а также на рисунке 3.

 

Таблица 6 - Показатели самолётов в результирующей шкале с учётом весов критериев

Показатели

Глобальные веса критриев

Оценки самолётов в однородной 100-балльной шкале с учётом весов

Эталоны

C1

C2

C3

C4

C5

C6

C*

C*

f1f11

0,200

20,00

10,60

12,00

13,20

0,20

8,20

0,20

20,00

f2f12

0,150

0,15

0,15

0,15

0,15

0,15

15,00

0,15

15,00

f3f13

0,150

11,40

12,90

3,90

0,15

4,05

15,00

0,15

15,00

f4f21

0,114

7,52

6,73

5,59

11,40

0,11

6,73

0,11

11,40

f5f22

0,099

0,10

9,90

4,46

0,10

4,46

9,90

0,10

9,90

f6f23

0,087

0,44

8,70

5,48

0,09

0,44

4,52

0,09

8,70

f7f31

0,072

0,07

0,43

6,77

6,34

6,98

7,20

0,07

7,20

f8f32

0,064

2,56

6,40

4,38

0,06

6,40

4,38

0,06

6,40

f9f33

0,064

2,56

3,07

0,06

6,40

0,06

2,75

0,06

6,40

i=19wgiril

44,80

58,88

42,79

37,89

22,85

73,69

1,00

100,00

 

Рисунок 3 – Результат оценивания самолётов в 100-балльной шкале с учётом весов критериев

 

В этом случае ранжирование самолётов по обобщённым оценкам с учётом весов критериев можно представить в виде:

МиГ35EurofighterTyphoonDassaultRafaleF16 INSuper ViperF/A 18E/FSuper HornetJAS 39NG IN.

На первом месте оказался самолёт МиГ-35 так как по трём показателям у него максимальные баллы. Кроме того, преимущество самолёта связано с наличием управляемого вектора тяги, который отсутствует у других самолётов.

По результатам агрегирования рейтинг боевых самолётов по эффективности представлен в таблице 7. Таким образом, самолёт МиГ-35 оказался лучшим по ТТХ и по стоимости.

 

Таблица 7 - Рейтинг самолётов по эффективности

Самолёт

Балльная оценка

Место в рейтинге

Без учёта весов

С учётом весов

Без учёта весов

С учётом весов

1

Dassault Rafale

37

45

5

3

2

Eurofighter Typhoon

60

58

2

2

3

F-16 IN Super Viper

45

43

3

4

4

F/18E/FSuper Hornet

40

38

4

5

5

JAS 39 NG(IN)

31

23

6

6

6

МиГ-35

63

59

1

1

Заключение

Основное преимущество предложенного МКП - в его простоте относительно МНК и МАИ, в отсутствии выполнения всех попарных сравнений, число которых нелинейно зависит от исходного числа объектов и равно (n2-n)/2. Снижение числа сравнений смежных пар критериев до значения (n-1) достигнуто за счёт предварительного упорядочения критериев по важности (предпочтительности) и особых свойств матрицы с наддиагональными экспертными оценками критериев, измеренных в шкале отношений.

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

 

1 PATTERN (от англ. Planning Assistance Through Technical Evaluation from Relevans Number) поддержка планирования посредством относительных показателей технической оценки.

×

About the authors

Viktor P. Korneenko

RUDN University; V.A. Trapeznikov Institute of Management Problems of the Russian Academy of Sciences

Author for correspondence.
Email: vkorn@ipu.ru
ORCID iD: 0000-0002-3643-1609
Scopus Author ID: 57222359318
ResearcherId: JGD-2807-2023

(b. 1953) graduated from the Faculty of Radio Electronics of the Leningrad Military Engineering Academy named after A.F. Mozhaisky in 1975. In 1986 he graduated from the Faculty of Computational Mathematics and Cybernetics of Lomonosov Moscow State University. Senior Researcher at the V.A. Trapeznikov Institute of Management Problems of the Russian Academy of Sciences, Associate Professor of the Department of Mathematical Methods and Information Technologies of the Higher School of Industrial Policy and Entrepreneurship of the RUDN named after Patrice Lumumba. Candidate of Technical Sciences, Associate Professor. The list of scientific papers includes more than 100 works in the field of combinatorial optimization, mathematical methods in analytical activity, and multi-criteria evaluation of objects with a multi-level structure

Russian Federation, Moscow; Moscow

References

  1. Saati T. Decision-making with dependencies and feedbacks: Analytical networks [In Russian]. M.: LKI, 2008. 360 p.
  2. Lopukhin MM. "PATTERN", a method of planning and forecasting scientific papers [In Russian]. M.: Soviet radio, 1971. 160 p.
  3. Keeney RL., Raiffa H. Decisions with Multiple Objectives: Preferences and Value Trade-Offs. New York: Wiley, 1976. 569 p.
  4. Belton V., Stewart TJ. Multiple criteria decision analysis. An integrated approach. Boston: Cluwer, 2003. 374 p.
  5. Bormotov AN. Substantiation of the method of forming the weighting coefficients of the criterion of practical opti-maleness based on the results of mathematical modeling of composites [In Russian]. Technical Sciences. 2016; 8: 14-18.
  6. Polishchuk LI. On generalized criteria with importance coefficients in vector optimization problems [In Russian]. Automaton. and telemech. 1982; 2: 55-60.
  7. Zotyev DB. On the problem of determining weight coefficients based on expert assessments [In Russian]. Factory laboratory. Diagnostics of materials. 2011; 77(1): 75-78.
  8. Anokhin AM. Glotov VA., Paveliev VV., Cherkashin AM. Methods for determining the coefficients of the importance of criteria [In Russian]. Automaton and telemech. 1997; 8: 3-35.
  9. Dmitriev MG., Lomazov VA. Evaluation of the sensitivity of linear convolution of particular criteria in the expert determination of weighting coefficients [In Russian]. Artificial intelligence and decision-making. 2014; 1: 52-56.
  10. Korneenko VP. Method of indirect preferences [In Russian]. Review of applied and industrial mathematics. 2008; 15(5): 890-891.
  11. Pfanzsagl J. Theory of measurement. Berlin, Heidelberg: Spriger-Verlag, 1971. 235 p.
  12. Podinovsky VV. Quantitative importance of criteria [In Russian]. AiT. 2000; 5: 110-123.
  13. Podinovsky VV. Ideas and methods of the theory of the importance of criteria in multi-criteria decision-making tasks [In Russian]. Moscow: Nauka, 2019. 103 p.
  14. Remesnik E.S.Application of Fishburne sequences in models with quantitative factors // Theory and Practice of Economics and Entrepreneurship / Proceedings of the XVI All-Russian Scientific and Practical conference with International participation "Theory and Practice of Economics and Entrepreneurship" (Simferopol–Gurzuf, 2019). Simferopol: IP Zueva T.V., 2019. P.210-212.
  15. Fishburn PC. Utility Theory for Decision Making. New York: Wiley, 1970. 234 p.
  16. Piyavsky SA. Forms for calculation of universal coefficients when adopting multiple critical decisions [In Russian]. Ontology of designing. 2019; 9(2): 282-298. doi: 10.18287/2223-9537-2019-9-2-282-298.
  17. Korneenko VP. Optimization method for selecting the resulting ranking of objects presented in the rank scale of measurement [In Russian]. Management of large systems. 2019; 82: 44-60.
  18. Korneenko VP. Methods of multi-criteria evaluation of objects with a multi-level structure of performance indicators [In Russian]. Moscow: MAKS Press, 2018. 292 p.
  19. Yushmanov SV. A method for finding weights that does not require a complete matrix of pairwise comparisons [In Russian]. AiT. 1990; 2: 186-189.
  20. Horn R., Johnson Ch. Matrix analysis. New York: Cambridge University, 1990. 561 p.
  21. Chu A., Kalaba RE., Springarn K. A Comparison of two methods for determining the weights of Belonging to fuzzy sets. Journal of Optimization Theory and Applications. 1979; 27: 531–538.
  22. Korneenko VP. Method of approximation matrix of object weights formation in multi-criteria selection problems [In Russian]. Bulletin of Cybernetics. 2021; 1: 51-62.
  23. Expert Choice. https://www.expertchoice.com/2021 (accessed: 04/25/2023).
  24. Comparative technical characteristics of aircraft participating in the Indian MMRCA tender [In Russian]. https://ru.wikipedia.org/wiki/Шаблон:Сравнительные_ТТХ_самолётов_принимавших_участие_в_индийском_тендере_MMRCA.
  25. Kendall MG. Rank correlation methods. New York: Oxford University, 1990. 260 p.
  26. Korneenko VP. Method for local aggregation of data of objects with multi-level structure in sequential scales [In Russian]. Proceedings of the 14th International Conference "Management of Large-scale systems Development" (MLSD-2021). Moscow: IPU RAS, 2021. P.485-493. https://mlsd2021.ipu.ru/proceedings/485-493.pdf.

Supplementary files

Supplementary Files
Action
1. JATS XML
2. Figure 1 – Three-level hierarchical tree of object evaluation criteria

Download (216KB)
3. Figure 2 – Hierarchical tree of aircraft efficiency evaluation criteria

Download (268KB)
4. Figure 3 – Result of aircraft evaluation on a 100-point scale taking into account the weights of the criteria

Download (189KB)

Copyright (c) 2024 Korneenko V.P.

Creative Commons License
This work is licensed under a Creative Commons Attribution 4.0 International License.

СМИ зарегистрировано Федеральной службой по надзору в сфере связи, информационных технологий и массовых коммуникаций (Роскомнадзор).
Регистрационный номер и дата принятия решения о регистрации СМИ: серия ФС 77 - 70157 от 16.06.2017.

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies