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

УДК 519.81

ВЫЧИСЛИТЕЛЬНЫЕ ЗАДАЧИ СИНТЕЗА МОДЕЛИ ВЫБОРА МЕТОДОМ МАКСИМАЛЬНОГО ПРАВДОПОДОБИЯ

Ю.В. Бугаев, Б.Е. Никитин, А.С. Чайковский

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

Ключевые слова: альтернатива, функция обобщенного критерия, метод максимального правдоподобия, плотность распределения

Введение

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

Одним из способов решения этой проблемы является метод экстраполяции экспертных оценок (МЭЭО). Он был предложен в [1].

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

Следовательно, данный подход представляется перспективным.

Обоснование метода простейших вариантов группового МЭЭО было дано в [3]. В этой работе рассматривались два способа проведения

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

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

допущений не тривиальна, и предполагает решение

Бугаев Юрий Владимирович - ВГТА, д-р техн. наук, профессор, тел. (4732) 55 - 25 - 50

Никитин Борис Егорович - ВГТА, канд. техн. наук, доцент, тел. (4732) 55 - 25 - 50

Чайковский Андрей Сергеевич - ВГТА, аспирант, Е- mail: [email protected]

ряда вспомогательных задач, решению которых и посвящена настоящая статья.

1. ММП в экстраполяции экспертных оценок

Исходной предпосылкой для МЭЭО является существование функции обобщенного критерия

р(хЪ) =^Ьи/и (х) = bтf{x), (1)

где ^ (х) - известные функции; Ъи - неизвестные

параметры (веса). На основе экспертного сравнения альтернатив из обучающей выборки строится система равенств и неравенств, описывающая

область возможных значений вектора Ъ

коэффициентов функции (1). Более подробное изложение ММЭО можно найти в работе [3].

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

С(г) м> > 0 , (2)

где ^ - вектор полезностей альтернатив; С(г), г = 1, к, N - матрица, отображающая структуру упорядочения г -го эксперта; N - число экспертов. Способ построения матриц С(г) покажем на

следующем примере.

Пример 1. Пусть имеем выборку из 4 альтернатив А1, к, А4. Предположим, что эксперт

проранжировал выборку на порядковой шкале следующим образом: А2 ^А1 ^А3 ^А4. В соответствие с данным упорядочением имеем систему неравенств: ^2 - w1 > 0; w1 - w3 > 0; ^3 - ^4 > 0;. Следова-

тельно, матрица С(1) , отображающая структуру упорядочения эксперта будет иметь вид:

C(1) =

- i i O O

i O - i O

O O i -i

Суть применения ММП основывается на предположении о независимости и нормальности распределения результатов экспертного оценивания. При этом значение экспертной оценки полезности I -й альтернативы wi интерпретируется как некая случайная величина с математическим ожиданием а1. Будем также считать выборку однородной, в силу чего оцененные полезности имеют одинаковую

дисперсию а , которую можно интерпретировать, как величину разногласия мнений экспертов. Каждый предложенный вариант ранжирования (2) интерпретируется, как некое случайное событие, вероятность которого вычисляется по формуле

Рг (в) = ^g(х, а,а)4х, Бг = {х е Ет | С(г)х > о] (3)

вг

где g (х, а, а) - плотность т - мерного нормального

распределения.

На основе этих данных строится функция правдоподобия (ФП)

L{kl, к 2,...,

N1

р1к1(в) ... р;* (в),

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

Подобным же способом можно найти ММП-оценки коэффициентов функции ОК, ранжирующего альтернативы адекватно экспертизе. Пусть вектор х1 соответствует векторной критериальной оценке 1 -й альтернативы и известна с точностью до коэффициентов, функция ОК (1). Тогда наложив ограничения Е(х', Ь) = а1, и определив соответствующую точку условного максимума ФП, можно получить ММП-оценки коэффициентов Ьи .

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

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

2. Вычисление конечных пределов и подбор процедуры численного интегрирования. Очевидно, области Бг, определяемые выражением (3), являются полубесконечными, поэтому необходимо свести пределы интегрирования к конечным, чтобы не использовать итерационные методы.

2. Замена переменных

Пусть ранжирование альтернатив по полезности г -м экспертом выражается неравенством (2). Сделаем замену переменных

1

^ = Bw = — С(г) w > 0. а

(4)

В соответствии с правилами преобразования нормального распределения [4], данная замена, в случае невырожденности матрицы ССт, преобразует плотность нормального распределения к виду

" 2 -1 0 . . 0"

-1 2 -1 . . 0

0 -1 2 . . 0

0 0 0 . . 2

§ Су, а ):

(аеі(ССт))

I (2п)п

1 --( у-а )т (ССт )-1( у -а )

-е 2 , (5)

где ё = — Са (верхний индекс матрицы С(г) опу-а

щен). Область интегрирования Бг, преобразуется при этом в неотрицательный ортант в п -мерном пространстве.

Проблема данной замены переменных состоит в необходимости обеспечения невырожденности

матрицы ССТ. Выполнение этого условия определяется видом экспертной упорядоченности. Рассмотрим несколько случаев упорядочения.

1. Полное упорядочивание по убыванию полезности. Система (2) примет вид:

V > ^2 > к > Wm ,

номера альтернатив для простоты соответствуют их месту в упорядочении. Этот случай соответствует первому варианту ранжирования примера 1, матрица С преобразования (4) имеет т -1 строк и т столбцов. Тогда

ССТ =

Следовательно, при полном упорядочении альтернатив получаем невырожденную матрицу ССТ и преобразование (4) допустимо.

Согласно известному свойству [5], ранг матриц С и ССТ совпадает. Следовательно, при приведенном упорядочении матрица С имеет полный ранг, равный числу ее строк.

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

1 четыре альтернативы могут быть проранжированы следующим образом: А2 уА3 уА4 ; А1 ^А4. В данном случае отношение предпочтения в парах (А1, А2) и (А1, А3) не установлено.

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

построенная на его основе система неравенств (2) непротиворечива и множество ее решений имеет непустую внутренность. При таком упорядочении возможны случаи вырожденности матрицы ССт .

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

Пусть имеем систему линейных однородных неравенств вида

С,1 Х1 + С,2 Х 2 + к + СгшХш > ^ 1 1 ., П . (6)

Обозначим X - область её решений. В общем случае часть неравенств является следствием ос-

тальных и их можно исключить из системы, область решений при этом не изменится.

Определение. Назовем существенными те неравенства системы (6), исключение которых приведет к изменению области решений. Соответственно, несущественными назовем неравенства, исключение которых не повлечет изменения области решений.

Построим правило обнаружения несущественных неравенств.

Теорема 1. Для того, чтобы к -е неравенство системы (6) было существенным, необходимо и достаточно, чтобы существовал ненулевой вектор и, удовлетворяющий всем неравенствам системы, кроме к -го.

Доказательство. Необходимость. Пусть к -е неравенство является существенным, а, значит, его исключение повлечет за собой изменение области X, причем, очевидно, новая область решений будет шире исходной. Это означает, что найдется вектор и , который удовлетворяет новой системе неравенств и не удовлетворяет старой. Но эти системы отличаются только к -м неравенством. Значит, и удовлетворяет всем неравенствам системы (6), кроме к -го.

Также очевидно, что тривиальное решение х = 0 удовлетворяет любому неравенству системы (6). Поэтому, чтобы вектор и обладал указанным свойством, он должен быть отличным от 0.

Достаточность. Пусть нашелся вектор и Ф 0, удовлетворяющий всем неравенствам системы, кроме к -го. Предположим противное, что к -е неравенство несущественное. Тогда, как было сказано ранее, любое решение системы (6), из которой исключено к -е неравенство, автоматически должно удовлетворять и к -му. Т. е. получаем противоречие.

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

С11и1 + С12и2 + ••• + c1mum c21u1 + С22и 2 + ••• + С2тит

Ск1и1 + Ск2и2 + • + Сктит <#

Сп1и1 + Сп2и 2 + • + Сптит

£ ^ max , -1 < и j < 1, j = 1,...,т ,

при к = 1,..., п . Последнее ограничение необходимо, чтобы избежать бесконечно больших значений и

J ■

Если в результате решения £

> O , то к -е

неравенство является существенным, если же

Г

O , то несущественным.

3. Вычисление интегралов

Рассмотрим интеграл (3), с плотностью распределения вида (5).

Тогда вероятность ранжирования примет вид

“+ад +ад

Po = JJ ^ J dx

(2n)

m J m-i •

0 0

+ад i t t i г --(y-d)T(CC )-i(y-d)

dxi

Как было сказано выше, для вычисления интеграла необходимо воспользоваться

безытерационной процедурой численного интегрирования. Одной из таких является

квадратурная формула Гаусса, но для того чтобы ею воспользоваться необходимо взять конечный

верхний предел, так чтобы остаток был < е. В одномерном случае имеем соотношение вида

ад р ад

^ / (х)ёх = ^ / (х)ёх + ^ / (х)ёх = 11 +12,

а а р

для которого можно определить значение р такое что |/2| < е . В многомерном случае точное значение остатка (12) рассчитать сложнее, поэтому будем определять его с запасом. Рассмотрим трехкратный интеграл

ад ад ад

J dxi J dx2 J f (Xi, x2, X3 )dx3

используя свойство аддитивности, представим его как сумму двух интегралов

ад ад

ад ад

^ йхх ^ ёх 2 | /ёх3 + ^ йхх ^ ёх2 | /йхъ .

а Ь ё А Ь ё

Проведя серию аналогичных преобразований, и в силу того, что /(х1, х2, х3) > 0 имеем

А В В А В ад

+

J dxi J dx 2 J fdx3 + J dxi J dx2 J fdx

a b d a b D

ад ад ад A ад ад

J dxi J dx 2 J fdx3 + J dxi J dx 2 J fdx

A b d a B d

A B D ад ад ад

J dxi J dx 2 J fdx3 + J dxi J dx 2 J fdx

a b d A -ад -ад

ад ад ад ад ад ад

J dxi J dx 2 J fdx3 + J dxi J dx2 J fdx

з -

3 -

-ад -ад D

А В В

< ^ ёх1 | ёх 2 | /йхъ +е а Ь ё

В этих интегралах порядок номеров переменных интегрирования, считая, от внутреннего интеграла меняется посредством циклического сдвига: (3, 2, 1), (1, 3, 2), (2, 1, 3). Таким образом, оценку каждого интеграла можно проводить по единой формуле, с учетом порядка номеров переменных интегрирования.

е

ад

+

+

-ад B

ад

Очевидно, в общем случае имеем оценку вида

Ai A2 Am

ад ад ад

Jdxi Jdx2 ... J fdXm < Jdxi Jdx2 ... J fdxn

ад ад ад

ад ад ад

J dxi J dx 2 ... J fdxm + к + J dXm J dxi к J fdxm-i

Ai -ад -ад

Am -ад -ад

Подберем пределы Ai из условия Ii < — , где

т

Ii - i -й по порядку интеграл. Имеем

/(y,d)

(det(CCT ))-1 -^У-^ )т (ccT )-i(y-d)

-------------Р 2

(2п)

Сделаем замену переменных х = y - d и

(CCT) 1 = Z , в результате получим

f ( xb---, xm ) =

(det(CC )) І (2п)”

T ))-1 -ЛxTZx

X1

; z

(7)

- матрица обратная ковариационной.

Представим матрицу 2 в виде

V Т . Т . . Т

2 = и^и^ + и^и2 + • к + и^ит

где вектор и' имеет (1 -1) последних координат

Г 5, ^

0

равных 0, а um =

V 0 У

где S1 - некоторый скаляр.

Подобное разложение соответствует представленной квадратичной форме хТ 2х в виде

Т2 ( Т )2 + ( Т )2 + + ( Т )2

х = (х и1) + (х и2) + • + (х ит) ,

где 1-ое слагаемое содержит все т переменных, 2ое - не содержит хт и т.д. (хТит )2 = х12512. Отсюда

ад ад ад

= J dxi J к J fdx2 ...

Ai -ад -ад

(det(CCT ))-1 f - і S? x2

(П

2 dxi

ад 1 т т -—x um_ium-ix е 2 dx 2

J

ад 1 т т С — x uiui x

Ге 2 dx„

(S)

Покажем, что 10 = Ф(-А151), где Ф(х) - интеграл вероятности. Согласно известной теореме

адад

теории вероятности

J к J fdx 2 к dxm = fi(xi)

-ад -ад

индивидуальная плотность распределения Х1 и

fi(xi)=■

i

i xi

a14ln следовательно

е у , но в силу (S) 1/ У1 = S1

ад 1 Л Л - Лл 1 Л Л

S1 f -Vxi2 S1 f -^S2Xi2

I1 = _____ і е 2 dx1 = ______ і е 2 dx1 =

» л/2п J

л/2П

Л1

i

= J е 2 dt = ФМ^). j

л/2п

Тогда, из условия /1 =— имеем А1 =—— Ф—

т 51 ^ т

Аналогично получаются все пределы А'. Для этого надо использовать разложение (7) с разным порядком следования исключаемых переменных интегрирования.

Литература

1. Пустыльник Е.И., Сысоев В.В., Чирко М.С. Использование линейной модели для экстраполяции экспертных оценок // Автоматизация проектирования. М.: МДНТП, 1981. С. 46 - 50.

2. Десятов Д.Б, Сысоев В.В., Чирко М.С. Принятие решений на основе экспертных оценок с использованием метода максимального правдоподобия // Автоматизация проектирования производственных систем. Воронеж: ВПИ, 1984. С. 32 - 36.

3. Бугаев Ю.В. Экстраполяция экспертных оценок в оптимизации технологических систем // Изв. АН. Теория и системы управления. 2003. № 3. С. 90 -96.

4. Андерсон Т. Введение в многомерный статистический анализ / Пер. с англ. М.: Физматгиз, 1963. 500 с.

5. Демиденко Е. З. Линейная и нелинейная регрессии. М.: Финансы и статистика, 1981. 302 с.

ai a2

a

m

+

ад

2

е

m

I

0

ад

-ад

Воронежская государственная технологическая академия

THE CALCULATING PROBLEMS OF THE CHOICE MODEL SYNTHESIS BY METHOD OF THE MAXIMUM CREDIBILITY

Yu.V. Bugaev, B.E. Nikitin, A.S. Chaikovsky

In article computing problems of synthesis of model of a choice by a method on the basis of extrapolation of expert estimations are considered. Definitions of insignificant inequalities are entered. The theorem of an exception of insignificant inequalities is formulated and proved. The method of calculation of multidimensional integrals is presented by a method of Gaussa

Key words: alternative; function of the generalised criterion; the method of maximum credibility; distribution density