автореферат диссертации по информатике, вычислительной технике и управлению, 05.13.17, диссертация на тему:Метод опорных объектов для обучения распознаванию образов в произвольных метрических пространствах

Абрамов Вадим Игоревич

На правах рукописи

4г

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

Специальность 05.13.17 — Теоретические основы информатики

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико-математических наук

2 7 НОЯ 2014

Москва, 2014

005556076

005556076

Работа выполнена на кафедре «Интеллектуальные системы» Федерального государственного автономного образовательного учреждения высшего профессионального образования «Московский физико-технический институт (государственный университет)».

Научный руководитель:

Официальные оппоненты:

Ведущая организация:

доктор технических наук, профессор Моттль Вадим Вячеславович

доктор технических наук, профессор

Дорофеюк Александр Александрович, Институт проблем управления РАН, заведующий лабораторией №55.

доктор физико-математических наук, профессор Двоенко Сергей Диниловнч, кафедра информационной безопасности, Государственное образовательное учреждение высшего профессионального образования «Тульский государственный университет».

Федеральное государственное унитарное предприятие «Государственный научно-исследовательский институт авиационных систем».

Защита состоится «18» декабря 2014 г. в 15:00 часов на заседании диссертационного совета Д 002.017.02 при Федеральном государственном бюджетном учреждении науки «Вычислительный центр им. A.A. Дородницына Российской академии наук», расположенном по адресу 119991, г. Москва, ул. Вавилова, 40.

С диссертацией можно ознакомиться в библиотеке ВЦ РАН.

Автореферат разослан «-//>> ноября 2014 г.

Ученый секретарь диссертационного совета

доктор физико-математических наук ^ Рязанов В.В

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность. Проблема обучения распознаванию образов составляет важнейший аспект современной теории машинного обучения (Machine Learning)1. Пусть в пределах генеральной совокупности объектов реального мира coefi всякий объект характеризуется скрытой от наблюдателя принадлежностью к одному из конечного множества классов = Природа выбирает некоторый объект и требует, чтобы наблюдатель

«угадал» класс объекта, всякий раз «наказывая» за ошибку j)(co)^ j(co). В данной работе мы ограничимся рассмотрением задачи распознавания образов для двух классов ¥ = {-1,1}.

Теория машинного обучения основана на понятии конечной обучающей совокупности, в которой известны истинные значения индексов классов объектов

и которая является единственным источником информации о скрытых свойствах природы, доступным наблюдателю. Задача обучения заключается в поиске такой аппроксимации у(со) функции v(co), известной лишь в пределах обучающей совокупности, которая позволила бы продолжить эту функцию на все множество объектов .у(со): П—>¥ = {-1,1}.

В данной диссертации гипотеза компактности, сформулированная Э.М. Браверманом2, рассматривается как основополагающий принцип машинного обучения. В самой общей формулировке для задач обучения распознаванию образов гипотеза компактности заключается в предположении, что объекты, отнесенные природой к одному и тому же классу, имеют, как правило, большее сходство между собой, чем объекты разных классов. Способ измерения сходства либо несходства объектов должен выработать наблюдатель.

В качестве математически корректной модели интуитивного суждения о несходстве объектов реального мира со s Q в данной диссертации рассматривается некоторая метрика, выбранная наблюдателем:

р(ю', со") : П х fi -> R, р(со', со") = р(ю", со') > 0, р(со', со") > 0, если со' # со", р(ш', со") + р(со", со") > р(со', со"). ' '

Выбор метрики удачен, если выполняется гипотеза компактности: когда р(со',со")—»0, то в большинстве случаев классы объектов совпадают _у(со')=_у(со").

В диссертации показано, что именно гипотеза компактности, выражаемая некоторой метрикой (2), лежит в основе чрезвычайно популярного метода машинного обучения, получившего в мировой литературе название кернельного метода (Kernel Based Approach)3, являющегося развитием известного метода потенциальных функций4. Под кернелом (потенциальной функцией) понимается всякая симметричная числовая функция, определенная на множестве объектов /T(co',co"):QxQ —»R, образующая неотрицательно определенную матрицу значений ^(co,.,coy).'J=l.—.w] Ддя всякой конечной совокупности объектов {co1,...,cov}cQ. Известно, что всякий кернел определяет погружение исходного множества объектов П в некоторое гипотетическое линейное пространство, быть может, бесконечномерное, на которое эта функция естественным образом продолжается, и в котором она играет роль скалярного произведения, определяя тем самым евклидову метрику.

Наиболее популярный метод обучения распознаванию двух классов объектов, получивший название метода опорных векторов (Support Vector Machiné) основан на идее поиска такой гиперплоскости в этом воображаемом линейном пространстве (дискриминант-ной гиперплоскости), которая обеспечила бы наиболее правильную классификацию объектов обучающей совокупности в смысле максимизации их евклидовых расстояний до гиперплоскости с нужной стороны от нее.

1 Вапник В.Н. Восстановление зависимостей по эмпирическим данным. М., Наука, 1979,449 с.

2 Браверман Э.М. Опыты по обучению машины распознаванию зрительных образов. Автоматика и телемеханика, 1962, т. 23, № 3, с. 349-364.

3 Vapnik V. Statistical Learning Theory. John-Wiley & Sons, Inc., 1998, 736 p.

4 Айзерман M.A., Браверман Э.М., Розоноэр Л.И. Метод потенциальных функций в теории обучения машин. М.: Наука, 1970, 386 с.

(1)

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

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

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

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

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

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

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

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

Задачи исследования. Для достижения поставленной цели в диссертации сформулированы и решены следующие задачи:

1. Разработка математического аппарата погружение метрического пространства с произвольной метрикой в псевдоевклидово линейное пространство.

2. Разработка параметрического семейства решающих прав™ различения точек двух классов в псевдоевклидовом линейном пространстве.

3. Разработка критерия обучения распознаванию объектов двух классов в произвольном метрическом пространстве как обобщение классического метода опорных векторов.

4. Разработка численных методов и алгоритмов обучения распознаванию объектов двух классов в произвольном метрическом пространстве.

5. Количественное экспериментальное исследование разработанных методов и алгоритмов.

Методы исследования. Теоретическое исследование базируется на общих принципах теории метрических пространств, аналитической геометрия и линейной алгебры, методе опорных векторов, методах выпуклой и невыпуклой оптимизации. Экспериментальное исследование проводилось с использованием программно-алгоритмического комплекса, разработанного автором.

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

1. Математический аппарат погружение метрического пространства с произвольной метрикой в псевдоевклидово линейное пространство.

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

3. Критерий обучения распознаванию объектов двух классов в произвольном метрическом пространстве, являющийся обобщением классического критерия опорных векторов.

4. Численные методы и алгоритмы обучения распознаванию объектов двух классов в произвольном метрическом пространстве.

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

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

Связь с плановыми научными исследованиями. Работа выполнена при поддержке грантов Российского фонда фундаментальных исследований №№ 11-07-00409-а, 14-07-00661-а, 13-07-13132 офимРЖД.

Апробация работы. Основные положения и результаты диссертации докладывались на конференциях «Интеллектуализация обработки информации ИОИ-2010» (Республика Кипр, г. Пафос, 2010 г.), «Интеллектуализация обработки информации ИОИ -2012» (Черногория, г. Будва, 2012 г.), «Интеллектуализация обработки информации ИОИ -2014» (Греция, о. Крит, 2014 г.), «Математические методы распознавания образов ММРО - 2013» (г. Казань, 2013 г.).

Публикации. По тематике работы опубликовано 4 статьи, в том числе 2 статьи в журналах, рекомендованных ВАК.

Структура и объем работы. Диссертация состоит из введения, 4 глав основного содержания, заключения и библиографии. Работа содержит 82 страницы основного текста.

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

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

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

Самым распространённым математическим инструментом измерения несходства между объектами реального мира сое£1 является представление объектов в компьютере векторами их числовых признаков х(ш) = (*,(<»),...,д:„(со))б М" с последующим вычислением евклидовой метрики между векторами признаков р(ю',ю") = р(х(со'),х(со"))=||х'-х"|[=

[(х'-х")г(х'-х")]^2. Пожалуй, наиболее популярным в мировой литературе методом обучения распознаванию образов с двумя классами объектов в пространстве векторов действительных признаков является метод опорных векторов. Его особое удобство при решении широкого класса прикладных задач объясняется двумя обстоятельствами.

Во-первых, получаемое решающе правило распознавания выражено в терминах не самих векторов признаков, а лишь их попарных скалярных произведений (х(со'))'х(со"): fi х Q -> М. Линейное дихотомическое решающее правило с/(х) (от английского слова decision) имеет вид дискриминаптной гиперплоскости в пространстве признаков

требуя вычисления скалярного произведения вектора признаков нового объекта х(ш) с векторами признаков объектов обучающей совокупности ху = х(со^) (1).

Во-вторых, как правило, большинство неотрицательных коэффициентов Xj, определяющих направляющий вектор дискриминантной гиперплоскости а = в ока~ зываются равными нулю, и сохранять в памяти достаточно лишь небольшое число векторов признаков остальных объектов обучающей совокупности J = {у : X,. > 0} с J = {1,..., N}, называемых опорными, которые дали название методу опорных векторов.

Это обстоятельство лежит в основе общего подхода к распознаванию образов, в котором предполагается, что объекты представлены лишь некоторой двухместной функцией их парного сравнения ^(co',cd"):QxQ—>R вместо индивидуальных векторов признаков x(co)eR". Такой способ представления объектов особенно адекватен широкому классу приложений, в

которых трудно выбрать числовые признаки отдельных объектов, но достаточно легко вычислить некоторую числовую характеристику отношения между любыми двумя объектами. Для того, чтобы процесс построения дискриминантной гиперплоскости (3) сохранил все преимущества исходного метода опорных векторов, традиционно принято считать, что функция ^(со',ю"):ПхП->К должна быть кернелом, т.е. образовывать неотрицательно определенные матрицы

Кл=[К(Ю;,ш,),у\/ = 1,...,ЛГ], crKvc>0, ceR* (4)

для любой конечной совокупности объектов, в частности, удовлетворять условию К(ы, га) > 0 при N = \. Всякий кернел погружает множество объектов реального мира юеГ2 в большее гильбертово линейное пространство ПэО, в котором он играет роль скалярного произведения. Решающее правило распознавания имеет в этом случае вид, аналогичный (3), с тем лишь отличием, что кернел K(ü>j , со) используется вместо скалярного произведения векторов признаков:

= ^>0, yei. (5)

Здесь уже нет векторов признаков в явном виде, поэтому метод (5) уместно называть методом опорных объектов.

В то же время, требование неотрицательной определенности для функции парного сравнения объектов оказывается слишком обременительным для многих прикладных задач анализа данных. Альтернативный подход был предложен Р. Дьюиным и его коллегами5,6 под названием реляционного дискриминантного анализа (Relational Discriminant Analysis) и независимо в работах ' . {беспризнаковое распознавание образов). Идея заключается в том, чтобы интерпретировать значения произвольной функции парного сравнения между всяким объектом coeQ и всеми объектами обучающей совокупности {co,,...,a>N} как вектор вторичных признаков этого объекта х(со) = (хДсо) = ^(со,,со),i = 1 ,...,N), и применить затем обычный метод опорных векторов в R,v (3):

Именно такой подход развивается в 9.

С формальной точки зрения решающее правило (6) остается в классе правил опорных объектов, поскольку каждая сумма ^"^(со^сoJ)K(a>i,co) выражает сравнение нового

объекта со е П только с опорными объектами обучающей совокупности {со;: j е J}. Однако применение правила (6) требует вычисления его вторичных признаков относительно всех объектов обучающей совокупности хДсо) = К(сй1.,со), г е J = {1,...,7V}, которые все равно надо держать в памяти. Последнее обстоятельство разрушает основное преимущество метода опорных объектов (5) — возможность не запоминать обучающие объекты, не являющиеся опорными. В работе 10 показано, что решающее правило (6) эквивалентно правилу

Duin R.P.VV, De Ridder D., Tax D.M.J. Featureless classification. Proceedings of the Workshop on Statistical Pattern Recognition, Prague, June 1997, pp.37-42.

6 Duin R., Pekalska E., De Ridder D. Relational discriminant analysis. Pattern Recognition Letters, Vol. 20, 1999, pp. 1175- -

7 1181.

Mottl, V., Dvoenko, S., Seredin, O., Kulikowski, C., & Muchnik, I. Featureless pattern recognition in an imaginary Hilbert space and its application to protein fold classification. Machine Learning and Data Mining in Pattern Recognition, Springer Berlin Heidelberg, 2001, pp. 322-336.

Середин О.С. Методы и ачгоритмы беспризнакового распознавания образов. Дисс. к.ф.-м.н. М, ВЦ РАН, 2001. Середин О.С. Линейные методы распознавания образов на множествах объектов произ-вольной природы, представленных попарными сравнениями. Общий случай. Известия ТулГУ, Естественные науки. Тула: Изд-во ТулГУ, Вып. 1,2012, с. 141-152.

Середин О.С. Потенциальная функция на множестве объектов распознавания как инстру-мент их попарного сравнительного представления Известия ТулГУ, Естественные науки. Тула: Изд-во ТулГУ, Вып. 1, 2013, с. 177-189.

опорных объектов (5) в том и только том случае, когда функция парного сравнения объектов К(а>, со") обладает свойствами кернела (4).

В первой главе разработан основной понятийный аппарат, лежащий в основе данного диссертационного исследования, и выдвинуты две основных концепции, развиваемые в диссертации.

Существование континуума кернелов, метрически эквивалентных всякому заданному кернелу. В первой главе доказано, что если на множестве объектов реального мира соеП определен кернел (потенциальная функция) .K(co',co"):f2xi2—>К, то существует континуум других кернелов, полностью эквивалентных ему в смысле решающего правила распознавания вида (5), причем с тем же множеством опорных объектов J с: J = {1,..., А^} для всякой обучающей совокупности (1), которые отличаются друг от друга только значениями коэффициентов Xj > 0, j е J.

Все эти эквивалентные кернелы определяют одну и ту же метрику на множестве объектов р(ю',со"): Q х Q —> R, удовлетворяющую, в дополнение к неравенству треугольника

р(со', со") + р(<в", со") > р(со', со"), (7)

еще и требованию условной неотрицательной определенности матриц

Pw=[-p2(my.m,),y,/ = l,...,JV], crPvc>0, lrc = X"lCy=0, (8)

для всех конечных совокупностей объектов. Метрики такого вида названы в настоящей диссертации пред-евклидовыми метриками (proto-Euclidean metric) 11 в отличие от стандартного понятия евклидовых метрик, определяемых в некотором линейном пространстве выбором скалярного произведения.

Недостаточность понятия пред-евклидовых метрик и, следовательно, кернела. В первой главе проведен анализ типичных классов прикладных задач беспризнакового анализа данных, которые приводят к необходимости использования метрик более широкого типа, нежели пред-евклидовые.

В частности, в задачах классификации биологических полимеров, аминокислотных цепей белков либо нуклеотидных цепей ДНК, общепринятым является способ оптимизационного выравнивания длин двух сравниваемых последовательностей (pair-wise alignment)12' 3. Аналогичный оптимизационный метод парного сравнения сигналов разной длины, например в задачах компьютерного анализа речи, получил в англоязычной литературе название dynamic time warping14' 5. Ещё одним примером задачи, в которой естественным образом порождается функция парного отношения между объектами анализа - это распознавание форм бинарных растровых изображений на основе сравнения их скелетных представлений' . Все эти методы приводят к некоторым метрикам на соответствующем множестве последовательностей либо сигналов (7), но эти метрики принципиально не являются пред-евклидовыми (8). Не обладает свойством пред-евклидовости и метрика на множестве динамических подписейГ8.

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

и Термин предложен К.В. Воронцовым.

12 Needleman S.B., Wunsch, C.D. A general method applicable to the search for similarities in the amino acid sequence of two

proteins. Journal of Molecular Biology, Vol. 48, Issue 3, March 1970, pp. 443-453.

13 Smith T.F., Waterman M.S. Identification of Common Molecular Subsequences. Journal of Molecular Biology, Vol. 147,

1981, pp. 195-197.

M Salvador S., Chan Ph. Toward accurate dynamic*time warping in linear time and space. KDD Workshop on Mining Temporal and Sequential Data, 2004, pp. 70-80.

15 Xiaoyue W., et al. Experimental comparison of representation methods and distance measures for time series data. Data

Mining and Knowledge Discovery, March 2013, Volume 26, Issue 2, pp. 275-309.

16 Местецкий Л.М. Непрерывная морфология бинарных изображений. Фигуры, скелеты, циркуляры. М.: Фшматлит,

2009,288 с.

17 Kushnir О., Seredin О. Parametric Description of Skeleton Radial Function by Legendre Polynomials for Binary Images

Comparison. ICISP 2014, Lecture notes in Computer Science 8509, A. Elmoataz et al. (Eds.),, pp. 520-530. Springer International Publishing Switzerland, (2014. '* Mottl V., Lange M., Sulimova V., Ermakov A. Signature verification based on lusion of on-line and off-line kernels. Proceedings of the 19th International Conference on Pattern Recognition. Tampa, USA, December 8-11, 2008, pp. 1-4.

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

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

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

В второй главе изложена предлагаемая методология погружение метрического пространства с произвольной метрикой в псевдоевклидово линейное пространство.

Пусть О - некоторое множество (генеральную совокупность) объектов реального мира с заданной на нем произвольной метрикой р^'.ю") (7). Выберем некоторый элемент ф е П в качестве «центра» метрического пространства и образуем двухместную функцию (ш\ со") = (1/2) [ р2 (со', ф) + р2 (со", ф) - р2 (со', со")] • (9)

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

Из определения общности (9) следует, что общность элемента метрического пространства сое £2 с самим собой равна квадрату его расстояния до центра Агф(оз,со) = р2(со,ф). Следовательно, расстояние между любыми двумя элементами однозначно определяется их общностями друг с другом и с самим собой независимо от выбора центра:

р2 (со', со") = А"ф (со', со') + А"ф (со", со") - 2 (со', со"). (10)

Теорема 1. Если ЛГф(со',со") - общность двух элементов метрического пространства относительно центра ф, то их общность относительно другого центра ф равна

К-(ы',(,Г) = К( (со', со") - АГф (со', ф) - АГф (со", ф) + (ф, ф). (11)

Теорема 2. Если'метрика на множестве объектов является пред-евклидовой {&), то для любой конечной совокупности объектов при любом выборе центра симметрическая матрица неотрицательно определена (4).

Из утверждения теоремы 2 следует, что всякая пред-евклидова метрика определяет множество кернелов на множестве объектов реального мира О.

Предположим для простоты, что множество П конечно |П| = А/, О = {ш],...,сом}, причем порядок нумерации элементов в дальнейшем не будет играть роли. Составим симметрическую матрицу значений общности относительно выбранного центра Кф=[^ф(со,,соу),

г,} =1 Число строк и столбцов этой матрицы может быть сколь угодно велико, но

обращаться к полному множеству ее элементов не придется ни в каких алгоритмах.

В силу симметричности матрицы Кф (М х М) все её собственные числа и собственные векторы действительны, причем собственные векторы попарно ортонормированы. Разумеется, собственные числа и векторы этой матрицы зависят от выбора центра ф е П:

.....!= 61

.....Чг

О, у.

(12)

Матрица Кф всегда может быть представлена в виде Кф = • Для произволь-:

ной метрики матрица К+, вообще говоря, не является положительно определенной, поэтому среди ее собственных чисел могут быть как положительные, так и отрицательные числа. Упорядочим собственные числа в порядке убывания, полагая, что 2,ф1 >0,...,с,фр >0, и

^Ф/ч+1 м Пару целых чисел принято называть сигнатурой матри-

цы Кф 19'20, которая в нашем случае является матрицей общности заданной метрики.

Теорема З2'. Сигнатура матрицы Кф не зависит от выбора центра метрического пространства феГ2.

Теорема 3 показывает, что сигнатура является характеристикой заданной метрики, а не выбора центра. Значения собственных чисел и, вообще говоря, зависят от выбора

центра ф е Г2, но после расположения их в порядке убывания число неотрицательных р и отрицательных q чисел остается неизменным в силу теоремы 3. Условимся все числа ¡;ф1 считать неотрицательными и введем обозначения:

=ЧФ.М> 0.

(13)

Составим матрицу (МхМ) из собственных векторов матрицы Кф как из столбцов гф=(гф1 •••гф м), и свяжем с каждым элементом со, конечного множества П — {со, ,...,юд(} действительный вектор х£ш=хф,=(хфЬ.,...,хфМ(.)еЕм, представляющий собой г'-ю строку матрицы Zф, элементы которой у = 1,...,р умножены на один и тот же коэффициент ( , а

дальнейшие элементы /' = р + \,...,М - на коэффициент . Введем в рассмотрение квадратную матрицу (Мх М), которую будем называть единичной матрицей сигнатуры р:

рхр рх(М-р)

И-р)хр ~1(М-рММ-р)

(14)

В этих обозначениях представление Кф = как разности Кф=^'']^ф|гф,гф/-

1 ^ ,2ф , позволяет записать матрицу общности элементов конечного метрического пространства в виде

кф=

Ф.1 р

X I У

\Т Т Г

р ЬМ у

\Х,миФ.1

иФ,миФ,.« )

(

•ХФ Г

и.

(15)

где векторы и(;бГ и уф|«К ' являются частями векторов хф|..

Таким образом, мы связали элементы произвольного конечного метрического пространства £1 = {й,,...,сйм}, в котором выбран центральный элемент феО, с М -мерными векторами действительных признаков элементов хфшеКм,...,хф еКм, определяемыми М собственными векторами матрицы значений общности этого конечного множества и М = р + (М - р) собственными числами >0, £ф/,., <0 (12), причем

центральному элементу соответствует нулевой вектор хф=0еК" .

19 Беклемишев Д.В. Курс аналитической геометрии и линейной алгебры. Учебное пособие для вузов. М., Высшая школа, 320 с.

* Умнов А.Е. Анаштическая геометрия и линейная алгебра: учебное пособие. М., МФТИ, 2011, 544 с. Доказательство теоремы получено Е.О. Черноусовой.

Для произвольных пар векторов (х',х") линейного пространства общее выражение для элемента матрицы значений общности (15) определяет двухместную числовую функцию ЛГ(х',х"): К'ихКА'->К согласно (14)

К(х',х") = (х'-хф)г .ТДх' - хф) = х'гЛрх". (16)

Здесь нижний индекс сознательно опущен К(х',\"), а не К^(х',х"), поскольку в силу теоремы 3 сигнатура не зависит от выбора центрального элемента ф е П.

Тем самым мы погрузили конечное метрическое пространство О = {со,,...,со,,} с центральным элементом ф е П. в М -мерное линейное пространство К'", в котором согласно (15) определена двухместная числовая функция (16).

Будем говорить, что такое линейное пространство К" натянуто на конечное метрическое пространство П = {со,,...,«)„}.

Рассмотрим свойства двухместной функции К(х', х") (16).

1)Симметричность К(х',х") = К(х",х'). (17)

2)Билинейность К(с'х' + с"х",х'") = с'К(х',х"') + с"К(х",х"'). (18) Эти два свойства совпадают со свойствами скалярного произведения в линейном пространстве, которое должно обладать еще и третьим свойством:

3) Неотрицательность при совпадающих значениях аргументов К(х,х)>0. (19)

Однако последним свойством двухместная функция К(х', х") (16) не обладает, поскольку для некоторых векторов хеК" может иметь место неравенство £(х,х)<0. Действительно, согласно (17) для произвольного вектора согласно (16)

^(х,х) = хг^х = (игуг)Г откуда следует, что

3') К(х, х) < 0 для векторов х е КЛ', таких что иги < \ту. (20)

Двухместную числовую функцию К(х', х") (16) принято называть индефинитным скалярным произведением с сигнатурой р + д = М, а само линейное пространство с индефинитным скалярным произведением - псевдоевклидовым пространством22. В частности, в случае сигнатуры р = М, q = 0, получается обычное евклидово пространство с обычным скалярным произведением, для которого К(х,х) >0.

Как следствие, в псевдоевклидовом пространстве 1ХЛ', натянутом на метрическое пространство П = {ю,,...,сом}, для всех пар векторов ((х',х"),х'т = (и'г,у'г), х"г = (и"г, у"г)) определено значение

/•2(х,,х") = (х'-х")г1р(х'-х") = (и'-и")г(и'-и")-(у'-у")г(у,-у"), (21)

которое может быть положительным г2(х',х")>0, равным нулю г2(х',х") = 0 даже при х' фх" , либо отрицательным г2 (х', х") < 0. Метрика в обычном понимании этого термина

г(х',х") = ^(и' - и")г(и' - и") - (V' - у"У (у' - у') (22)

определена только для тех пар векторов, для которых (и'-и")г(и'-и")-(у'-у")г(у'-у")>0. Если же (и'-и")г(и'-и")-(у'-у")г(у'-у")<0, то значение метрики оказывается мнимым как квадратный корень из отрицательного числа. . -

Таким образом, построенное псевдоевклидово пространство Е" не является, вообще говоря, метрическим, поскольку в нем не для всех пар точек х',х" е МЛ' определено значение метрики г(х',х"), но значение квадрата г2(х',х") (21) определено для всех пар векторов, являясь отрицательным для некоторых из них.

2* Азизов Т.Я., Копачевский Н.Д. Введение в теорию пространств Крейна. Специальный курс лекций. Симферополь, ООО Форма, 2010, 112 с.

Изометрический образ метрического пространства в псевдоевклидовом линейном пространстве. В псевдоевклидовом линейном пространстве, натянутом на конечное метрическое пространство П = {со, ,...,со„}, индефинитное скалярное произведение любых двух векторов х'еК" и х'ей" отличается от обычного скалярного произведения х'Тх" только единичной матрицей сигнатуры /> + <? = М (16). В то же время собственные числа и собственные векторы матрицы значений общности Кф для конечного метрического пространства с выбранным центральным элементом феП (12) определяют координаты конечного множества векторов в линейном пространстве {хфш, >—>хф <„„ со е П => хфш е Км, зависящих от выбора центрального элемента, в которые отображаются элементы самого метрического пространства и которые играют роль базиса в этом линейном пространстве. Условимся обозначать это конечное подмножество в Л" символом

= >■■•,*♦,,„}<= К" (23)

и называть образом метрического пространства в псевдоевклидовом линейном пространстве. Разумеется, образ метрического пространства (23), существенно зависит от центра через собственные числа и собственные векторы матрицы общности (12).

Как говорилось выше (22), метрика в обычном смысле определена в псевдоевклидовом пространстве Ем только для тех пар векторов (х',х"), для которых (и'-и")'(и' - и") -(у'-у")т(\'~ у")>0. Однако векторы из конечного образа метрического пространства ®п ф (23) являются именно такими. Действительно,

Г (^,0)4 >Хф,ю, ) = (хф.ш, — Хф,о>,) — хф.(0, ) = Хф,(04 ■^рХф.ш, + Хф,ш,'"рХф,ы, — ^Хф.ш, >

т.е. согласно (15) и, далее, (10)

г1 (хф,м, • хф,.., > = к.*' мф.* >'+ ^ф к,', ) - 2ЛГф (юфд, шф,) = р2 (ю4, ) > 0, со,, со, е П. (24)

Обратим внимание на тот факт, что хотя линейное пространство Е " , в которое погружено исходное конечное метрическое пространство П = {а)1,...,<оЛ/}, не зависит от выбора центрального элемента феГ2, погружение (23) существенно зависит от ф, однако изомет-ричность погружения опять же никак не зависит от этого выбора.

Частный случай: Погружение метрического пространства с пред-евклидовой метрикой в евклидово линейное пространство. Предположим теперь, что исходная метрика удовлетворяет, кроме неравенству треугольника (7), еще и требованию условной неотрицательной определенности матриц Р„ =^-р2(со>,со(),7,/ = 1,...,Аг] для всякой конечной совокупности элементов (8), т.е. сгРЛ,с>0, 1 гс = с1 = 0. Такая метрика названа выше пред-евклидовой.

Теорема 4. В случае пред-евклидовой метрики матрица Кф=£лГф(со,.,о)Д г\/ = 1,..., А/^

значений общности элементов конечного метрического пространства (9) неотрицательно определена для любого центра ф е П.

В силу утверждения этой теоремы, все собственные числа матрицы Кф неотрицательны ¡;ф| >0,...,^фм >0, и сигнатура этой матрицы равна р = М, (? = 0. Следовательно, матрица Кф может быть представлена как матрица обычных скалярных произведений

( Т т \

кф= ; ; .

т т

X X X X

Итак, как и в общем случае, элементы произвольного конечного метрического пространства П = {со|,...,со„}, в котором выбран центральный элемент феО, оказались связаны с М-мерными векторами действительных признаков элементов хф еМ",..., хф<0 еКЛ', определяемыми М собственными векторами матрицы значений общности этого конечного множества и М собственными числами ¡;ф1 >0,...,%фм >0, причем центральному элементу соответствует нулевой вектор хф= ОеКм .

Отличие частного случая евклидовой метрики от общего случая произвольной метрики заключается в том, что теперь метрическое пространство погружено в обычное евклидово линейное пространство, а не в псевдоевклидово, как прежде, т.е. в линейное пространство с обычным скалярным произведением, обладающим всеми тремя свойствами (17), (18) и (19). Соответственно, метрика (21)-(22) определена для всех пар точек х',х" е М" .

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

Рассмотрим произвольную конечную неупорядоченную совокупность {х,,...,хп} элементов псевдоевклидова пространства х^ 6 К" (15), п<М. Пусть с = (с|-"с<)£3" - вектор коэффициентов при элементах совокупности, в сумме составляющих единицу У^= 1гс = 1, где 1 еМ" — вектор, составленный из единиц. Линейная комбинация

(25)

называется аффинной комбинацией элементов {х,,...,х„} коэффициентами с = (с, ■■■сп) .

Теорема 5. Квадрат расстояния любого элемента псевдоевклидова пространства х е К4' до аффинной комбинации хс е Мм (25) согласно (21) определяется равенством

'•2(хс^) = Е1.,с,'-2(х4,х)-(1/2)Х;.1Е,11с,с/'-2К.х/)- (26)

Если совокупность {х,,...,хп} составлена из произвольных элементов псевдоевклидова пространства, и элемент х также выбран произвольно, то среди квадратов расстояний /-2(х,,х,) и г2(хл,х) могут быть, вообще говоря, отрицательные, и квадрат расстояния т-2(хс,х) может также оказаться отрицательным.

Но если в качестве х и \к выступают векторы, соответствующие объектам исходного метрического пространства х = хф ш, х^ = хф ш , со, <в4 е П,

г2(хс,ХФ, „) = х 1.1с^2(ХФ.»«'х*.»> -0/2)Е1=, ^ (*Ф,И<'»♦.«,) • (27)

то расстояния КхфЮ1,хфш) = р(со*,сй) и г(хф(1)1,хфл) = р(со4,сй,) являются вещественными, ^ • х».„) = Р2 К ,со) > 0 и г2 (хф_М1,хф ) = р2(ш,,со,) > 0.

Тем не менее, метрические расстояния г(хс,хфш) существуют только в том случае, если г2 (хс,хф и) > 0 для любых с. Для произвольной метрики на множестве П последнее условие, вообще говоря, не выполняется, и существуют совокупности элементов исходного метрического пространства со, со,,..., о)л еП, такие, что г2(хс,хфю)<0. Более того, для всякой такой совокупности может существовать множество аффинных комбинаций хс е К" (25).

Это означает, что вектор хс е КА/, являющийся аффинной комбинацией векторов хфм ,...,хфш е К'" с некоторыми коэффициентами с (25), может не иметь прообраза ни в каком расширении метрического пространства О, т.е. может оказаться невозможным даже мысленно добавить соответствующий элемент шс в Г2, поскольку для него определены

квадраты расстояний до всех других элементов ?"2(хс,хф0)), но не могут быть определены метрические расстояния г(хг,хфш) = ^г2(хс,хфю), если г2(хс,хфщ)<0.

В частности, если п = 2, то коэффициенты аффинной комбинации векторов х' б Нд' и х'еК" определяются одним действительным числом с, например, с' = 1-с и с" = с. Тогда согласно (25)

хс =(1-с)х' + сх"еЕл'. (28)

Такой вектор будем называть соосным упорядоченной паре векторов <х',х">. Согласно (26) для соосного вектора г2(хс,х') = с2г2(х',х"), 7-2(хс,х") = (1-с)2г2(х',х"). Если г2 (х',х") > 0, то определены метрические расстояния

r(xc,x')=|c|Kx',x'), Кхс>0=|1-ск(х',х"). (29)

Для любого вектора х е Ш"

г2 (Хе, X) = (1 - с)г2 (х\ X) + сг2 (х", х) - с(1 - с)г2 (х', х"), (30)

а квадраты расстояний соосного вектора до исходных векторов определяются равенствами r2(xc,x') = cV(x',x"), r2(xf,x") = (l-c)V(x',x"). (31)

Вообще говоря, квадрат расстояния г7(х',х") может быть отрицательным, как и квадраты расстояний г2(хс,х'), г2(хс,х") и г2(хс,х).

Но если в качестве х' и х" выступают векторы, соответствующие объектам исходного метрического пространства х' = хф01. и х" = хфш., со', со"еQ, то расстояние '"(хфш.,хф(0.) =

р(ш',со") является вещественным, г2(хфш.,хфт.) > 0, тогда согласно (31) г2(хс,хфм,)>0 и г2 (хс, хф 0). ) > 0. Для векторов хф 0), соответствующих другим элементам метрического пространства coeQ,

г2 (хс - хф,ш ) = 0 - СУ (x+v. хф,<„ ) + с2 (ХФ,Ш-. хф,ш ) - с(1 - с)г2 (хф , хф ) =

(1 - с)р2 (со', со) + ср2 (со", со) - с( 1 - с)р2 (ш', со"). Однако для произвольной метрики на множестве Ù существуют тройки элементов исходного метрического пространства со,ев',со"еQ, такие, что г2(хс,хф())<0, и невозможно даже мысленно добавить соосный элемент <вс в П, поскольку для него не определены метрические расстояния до других элементов г(хс, хф ш ) = ^/-2(хс,хфо>) .

Аффинное псевдоевклидово пространство, натянутое на изометрический образ метрического пространства. Аффинная комбинация (25) определена для любой конечной совокупности {х,,...,хп} векторов псевдоевклидова пространства R", в частности, для образа метрического пространства Кпф={хфШ1 ,---,хф Ши }cl," (23), т.е. для совокупности всех M векторов, соответствующих элементам исходного метрического пространства Q = {со, ,...,сом} с выбранным центральным элементом феП. Всякая совокупность действительных коэффициентов (с, ■ • • си ), в сумме составляющих единицу = 1, определяет вектор x^^^q.x^eR*', зависящий от выбора центрального элемента, но находящийся от других

векторов х е RM на некоторых расстояниях (27), полностью определяемых коэффициентами (с, -"Cj,). Хотя вектор аффинной комбинации хс и зависит от центра, квадрат его расстояния г2 (хс, х) до всякого вектора xeR:U, быть может, отрицательный, от центра не зависит.

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

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

Пусть по-прежнему Q есть множество объектов реального мира с заданной на нем метрикой (7). Обучающая совокупность (1) представлена в компьютере как матрица попарных расстояний между объектами с вместе с известными индексами их принадлежности к одному из двух классов

{р(со;,со/),у,/ = 1,...,Л', yl=y(a,j) = ± 1,у = 1.....JV}. (32)

Целью наблюдателя является построение решающего правила распознавания классов новых объектов соеП, не представленных в обучающей совокупности, причем единственным свойством каждого нового объекта, доступным наблюдателю, является совокупность его расстояний до объектов обучающей выборки p(coj ;co), j - 1,..., N.

Пусть R" — псевдоевклидово линейное пространство, натянутое на конечное множество объектов реального мира |Q| = M. Это псевдоевклидово пространство однозначно

«привязано» к метрическому пространству О., поскольку в силу теоремы 3 его сигнатура фиксирована и определяется только метрикой р(ю',со"):ОхГ2—»К. Всякий выбор некоторого элемента ф е П в качестве центрального ставит в соответствие каждому элементу со е О соответствующий ему вектор хф ш е , т.е. определяет изометрический образ М^СК" (23) метрического пространства в ЕЛ', т.е. г2(Хф М1,хф(|) )=р5(со)1,ш;) (24).

Дискршшнантпьш диполем будем называть упорядоченную пару векторов <а_,,а, >еКА/хК", а сами векторы а_,,а, еКм - полюсами диполя23. Рассмотрим множество всех векторов хс еКм, соосных паре <а_,,а, > в смысле (28) со всеми действительными коэффициентами хг =(1-с)а_, +са, еМм. Центральной точкой диполя будем называть вектор х1/2 = (1/2)(а_, - а,), с = 1/2. Условимся рассматривать только такие диполи, квадрат расстояния между полюсами которых является положительным и равен единице согласно (21) г2(а_,,а,) = г(а.,,а1) = 1.

Пусть х е КА' - произвольный вектор в псевдоевклидовом пространстве, например, вектор хф т, соответствующий некоторому объекту реального мира со е П. Тогда формула (30) определяет квадрат расстояния между хс. и х:

г2 (хс, х) = (1 - с)г2 (а_,, х) + сг2 (а,, х) - с( 1 - с)г2 (а,,, а,). (33)

В силу предположения г2(а_,,а,) = 1 эта функция является квадратичной и строго выпуклой функцией действительного коэффициента сеК., но может принимать, вообще говоря, и отрицательные значения. В частности, отрицательным может быть ее минимальное значение штге,г!(1сд). Тем не менее, пусть с1=аг£тшсбЯг2(хс,х) - точка минимума, тогда вектор х; =(1-с,)а_| +с5а, еМ'" естественно называть проекцией точки хеК" на луч, образованный диполем < а_,, а, >, вектор х; характеризуется метрическим расстоянием до центра диполя

»■(*;. >х,/2 I ) = -1/2|г{а_,,а,) = |<3, -1/2|. (34)

Центральная идея методологии обучения распознаванию образов в произвольных метрических пространствах, предлагаемая в данной диссертации, заключается в использовании расстояния (34) с учетом знака (с1 —1/2) как параметрического семейства дискрими-нантных функций, каждая из которых задает некоторое разбиение псевдоевклидова пространства КЛ' натри части ¿/(х |а_,,а,) = (сх - 1/2)г(а_1,а,) = 0, определяемое выбором

диполя < а_,,а, >еК'и х»1'.

Теорема 6. В метрических терминах дискриминантная функция ¿/(х|а_|,а,) = (¿\ -1/2)г(а_,,а,) имеет вид

1 / - Г> 0 => положительная часть,

с/(х | а_,, а,) = —( г2 (ач, х) - г2 (а,, х Ш = 0 => нейтральная часть, (35)

2 [< 0 => отрицательная часть.

Будем называть нейтральное множество, определяемое согласно (35), дискриминант-ной гиперплоскостью в псевдоевклидовом пространстве К" :

Ща_, ,а|) = {хеКм:г2(а„х) = г2 (а,, х)} = {х е Км: с!(х | а_,, а,) = о} . (36)

В частности, при конкретном выборе центрального элемента ф е С1 каждому реальному объекту со е Г2 соответствует вектор хф 0) е КЛ', поэтому дискриминантная функция в псевдоевклидовом пространстве (35) фактически определяет дискриминантную функцию на множестве объектов реального мира:

2 [> 0 => класс 1,

¿/(со|а_,,а,) = — (г2(а_1,хфш)-г2(а,,хфо))н=0=> нейтральное решение, (37)

2 [< 0 => класс — 1.

*3 Терминология предложена К.В. Воронцовым.

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

Представляется естественным искать дискриминантную функцию, наилучшим образом разделяющую обучающую совокупность (32) в смысле ¿/(со. |а_,,а,)^0, выражая узлы дискриминантного диполя как неизвестные аффинные комбинации векторов {хф ,...,хфю }, в которые отображаются объекты самой обучающей совокупности:

Теорема 7 .Для всякой точки и диполя .а^М1' (38) выполняется равенство

В равенстве (39) только первая сумма в правой части зависит от предъявленного объекта шбП, причем в качестве коэффициентов выступают разности с_, ] - с, /, сумма которых для любого диполя а_,,а, е должна равняться нулю согласно (38):

а1 = с-и _си> Х>у=°- (4°)

Теорема 8. Требование единичной длины диполя выражается равенством

(У2)Е,1,ХГ-,(-Р2(шу,т,))*/», = 1 ■ (41)

Значения коэффициентов (а1,...,аг/) (40) определяют ориентацию диполя в псевдоевклидовом пространстве относительно образов объектов обучающей совокупности {хфЮ1 ,...,хф(0 }, оставляя свободными как «параллельный перенос» диполя, так и его «сдвиг» вдоль своей оси. Именно этот «сдвиг» и характеризует вторая двойная сумма в правой части (39), которая является константой по отношению к предъявленному объекту ю е . Обозначим ее символом

Ь = -с,Л,)р2(ш.,со,). (42)

Мы получили эквивалентное выражение для дискримирантной функции, которая, как оказалось, полностью определяется N + 1 действительными числами (о,,...,аЛ.,6) при двух ограничениях типа равенств:

J у=1 2

Метод опорных объектов для обучения распознаванию образов в произвольном метрическом пространстве. Практически буквальным выражением принципа оптимальной дискриминантной гиперплоскости, лежащего в основе метода опорных векторов, является критерий обучения в евклидовом метрическом пространстве, требующий максимизации зазора (марджина в терминологии В.Н. Вапника) между объектами двух классов в обучающей совокупности у )>£ > 0, я —> гпах . Мы будем использовать здесь модификацию этого принципа, предложенную в кандидатской диссертации А.И. Татарчука24, заключающуюся в том, что всякое нарушение этого неравенства «наказывается» функцией потерь

В качестве критерия обучения естественно искать такую гиперплоскость, которая разделяла бы обучающую выборку на два класса, с одной стороны, с как можно большей величиной зазора 1/е2—»гшп, а с другой, с как можно меньшей величиной суммарного штрафа

24 Татарчук А.И. Байесовские методы опорных векторов для обучения распознаванию образов с управляемой селективностью отбора признаков. Дисс. к.ф.-м.н., ВЦ РАН, 2014.

для ошибочно классифицированных объектов обучающей выборки . Баланс таких

требований к процессу обучения выражается оптимизационным критерием \/г~г+ С^^=15у->тт, где С - структурный параметр, определяющий соотношение требований

максимизации зазора и минимизации суммы потерь. С учетом ограничения (43) такой критерий обучения относительно искомых узлов диполя естественно записать как следующую задачу оптимизации с ограничениями в исходном метрическом пространстве:

' 0/2)уу + фе(1-5у), (45)

Очевидно, что задача обучения инвариантна к изменению масштаба приятой метрики у р2 (со', со"), если пропорционально изменить также порог функции потерь у е:

(№[ЕГЛ(-Р2(юясо,))а, 6(1-5,). (46)

Замена переменных а, = а,/2ъ,Ь = Ь/2е и С = С/2 позволяет избежать невыпуклой области варьирования переменных и записать критерий как задачу квадратичного программирования, близкую к общепринятой традиции метода опорных векторов:

Здесь не учтены ограничения на коэффициенты при объектах обучающей совокупности

у

(45), однако решение задачи будет автоматически удовлетворять этим условиям.

В отличие от традиционного метода опорных векторов, это невыпуклая задача квадратичного программирования, поскольку квадратичная форма в целевой функции не является условно неотрицательно определенной при Задача становится выпуклой,

только если заранее известно, что метрика является пред-евклидовой (8).

В силу невыпуклости задачи обучения можно говорить лишь о поиске локального минимума критерия с линейными ограничениями-неравенствами (47). Такая задача эквивалентна поиску седловой точки функции Лагранжа, аргументами которой, наряду с целевыми переменными (я1,...,ал,,6,51,...,5л), являются также неотрицательные множители Лагранжа, которые обозначим как (А.,,...ДЛ,) для первой группы ограничений и (ц,,...,цЛ,) для второй.

Теорема 9 .Двойственная форма задачи обучения (47) имеет вид:

-ОАОЕ^Е?., зол

Ее решение (X.,,... ДЛ,) полностью определяет параметры решающего правила распознавания, удовлетворяющие необходимому условию минимума критерия (47):

с1 /2>Е;„<,.<лЕл,,>„р2(®я«>^А,"

= о/2^А> =°'Ь =-"-^-Г-"-. (49)

и значение максимального зазора в (45) 1/е2 = (1/2)^ Хм^'Л (~Р2 • (50)

Однако в случае произвольной метрики решение «наивной» двойственной задачи (48) может привести к отрицательному значению квадрата максимального зазора е2<0 в (50). Требование максимизации зазора 1/е2—>тт в исходном критерии (45) равносильно требо-

ванию минимизации квадратичной формы (50), которая в общем случае не является условно неотрицательно определенной, и может принимать сколь угодно большие по модулю отрицательные значения. В то же время по своей сути зазор является действительным числом, т.е. е2>0, поэтому требование >т5п должно быть дополнено в задаче обучения (47) ограничением Если метрика

является евклидовой, то это ограничение никогда не будет нарушаться и никак не исказит идею обучения.

Учет этого ограничения равносилен введению в двойственную задачу (48) дополнительного ограничения (-р2(<ау,со/))Х.;Х;>0. Вместо строгого ограничения в диссертации используется штрафная функция, «почти точно» выражающая его требование:

/(Х1,...ДЛ,) = ДехР{-£>Х^,Е;>уЛ(-р2(соу,со,))^,}, (51)

где О - достаточно большое число, значение которого само должно ассоциироваться с «очень большим» штрафом, например, О = 10'°.

Действительно, такая функция принимает большое /(*■......>•»)

значение И при нулевом значении аргумента (квадратичной формы), быстро стремится к нулю при даже небольших положительных значениях, и уходит в бесконечность при отрицательных значениях аргумента: £

->оо.

7=1 Í=1

=я>>0, ¿ ¿ (и/))Л А =0'

j=i ы

М N

Е Ё y¡y< (V (°v А

Рис. 1. Штрафная функция, не допускающая отрицательных значений J-l '-1 квадратичной формы.

С учетом штрафной функции предлагается вместо «наивной» двойственной задачи (48) решать задачу

\ Л. Л. . . Г N N

4

X h-7ХХ УуУ/ (-Р2 К' ю,)) \Л ~D ехР \ -DX%yJy> ("р2 (со;' ))'XJX' j->maxO.. - А* )>

7=1

*7=1 (=1

7=1 /=1

(52)

J=|

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

Обучение можно считать удачным, если в результате решения двойственной задачи (52) удалось найти решающее правило (43), правильно классифицирующее объекты обучающей совокупности с достаточно большим значением зазора между классами е (50), т.е.

с малым значением Это значение будет гарантирован-

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

Класс метрических дискриминантных решающих правил возрастающей сложности. Очевидно, что разделяющая способность класса решающих правил í/(co|a,,...,ov,¿)^0 (43), выражаемая максимально достижимым зазором е на фиксированной обучающей совокупно-

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

p(m',(0'|ß) = -^[l-exp(-ßp2(c0>"))]'/2, ß>0. (53)

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

Теорема 10. Если р(со',со") есть метрика на ninim

П, то р(со',со"|Р) тоже метрика при любом «..................

ß > 0, причем р(со',со" | ß) —> p(co',co")«/)i/ß —> 0.

Поскольку (53) есть метрика, то все сказанное в предыдущих разделах остается справедли- _!_,?, вым по отношению к ней. Непосредственная Ja подстановка (53) в решающее правило классификации объектов метрического пространства (43), в исходную «наивную» задачу обучения (45) и далее во все выражения вплоть до двойственной : задачи (52) приводит к простой замене

(-р2К,со,)) на (l/Vß)[-l + exp(-ßp2(co>"))],

причем в силу замечания (46) коэффициент не с --

имеет значения, так что подставлять достаточно

[-l+exP(-ßPV,0)] • РИС' 2" Пре°бра30ВаНИС МС1РИКИ-

Кроме того, учет условий и' соответственно, в (48) делает

подстановку еще более простой - достаточно заменить (-р2(со;,со;)) на exp(-ßp2(co',co")), и использовать класс решающих правил

Таким образом, вместо одного параметрического семейства решающих правил (43) при ß —> 0 и одной задачи обучения (47), мы получили класс параметрических семейств решающих правил и соответствующих им задач обучения, определяемых выбором еще одного дополнительного параметра ß > 0.

Обучение во вложенных семействах дискриминантных решающих правил возрастающей сложности. С ростом параметра ß увеличивается «свобода» выбора решающего правила классификации в процессе обучения. В силу последнего обстоятельства параметр ß уместно назвать структурным параметром критерия обучения. Все сказанное выше есть частный случай при ß —> 0.

В третьей главе доказано, что для любой метрики р(со', со"): О х Г> —> R существуют такие обучающие совокупности (32), которые не могут быть правильно классифицированы ни при каких значениях параметров решающего правила (ait...,aN,ti) (43). В то же время, при достаточно большом значении структурного параметра ß—>га существует решающее правило (54), правильно классифицирующее любую обучающую совокупность. Выбор подходящего значения ß по единственной обучающей совокупности возможен только на основе некоторого внешнего метода контроля обобщающей способности обучения, например, путем применения той либо иной процедуры кросс-валидации.

В четвертой главе рассмотрены пути численной реализации двойственной задачи обучения распознаванию образов в множестве объектов с произвольной метрикой (52), (49)-(50), и приводятся результаты экспериментального исследования предложенных в диссертации методов обучения.

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

n/r, R1

«злостной», во всех экспериментах глобальный минимум критерия достигался в результате применения обычных итерационных методов, разработанных для решения задач выпуклого программирования. В данной работе использовался пакет программ, созданный Н.А. Разиным2 и основанный на стандартном методе внутренней точки.

Экспериментальное исследование разработанной методологии обучения распознаванию образов в массивах объектов, рассматриваемых как метрические пространства с произвольной метрикой, проведено на примере практической задачи верификация личности по динамической подписи (on-line signature). Задача заключается в проверке нулевой гипотезы о том, что рассматриваемая подпись действительно принадлежит заявленному автору (genuine signature), против альтернативной гипотезы, что подпись является сознательной подделкой (skilled forgery).

Рассматриваются динамические подписи из базы данных SVC 200426, каждая из которых вводится в компьютер непосредственно в процессе написания (online), и представлена многокомпонентным дискретным сигналом со индивидуальной длины, отражающим ее геометрические и динамические особенности. Степень несходства подписей р(со',ш"), играющая роль метрики, вычисляется на основе парного выравнивания соответствующих сигналов разной длины27. Такая метрика, вообще говоря, не является евклидовой в том смысле, что условие (8) может не выполняться для некоторых совокупностей объектов. Проявление этого факта мы увидим в результатах эксперимента.

Массив содержит динамические подписи 40 лиц z е Z по 20 для каждого из них, причем 10 подписей являются настоящими и еще 10 подделками. Таким образом, обучающая совокупность для каждой персоны QI=|û)y,y=l,...,A^J состоит из N = 20 подписей, представленных матрицей (20x20) попарных расстояний р(юу,ю,) и снабженных индексом класса 1 (настоящая) либо у =-\ (подделка).

Общий эксперимент заключался в обучении верификации истинности подписи каждого лица zeZ = {l,...,40} согласно (52), (49)-(50) и далее (54) с выбором структурного параметра преобразования метрики (3 по (53) методу скользящего контроля. В 37 частных экспериментах скользящий контроль дал значение структурного параметра, близкое к нулю, что говорит о хорошей разделимости обучающей совокупности с исходной метрикой. В остальных 3-х экспериментах исходная метрика не обеспечила разделимости обучающей совокупности, и по скользящему контролю были выбраны значения р = 0,3 ч- 0,6.

Ошибка скользящего контроля составила в среднем 3,65% и колебалась от 0% (21 случай) до 15% (2 случая).

В заключении перечислены следующие основные результаты:

1. Разработан математического аппарата погружение метрического пространства с произвольной метрикой в псевдоевклидово линейное пространство.

2. Разработано параметрического семейства решающих правил различения точек двух классов в псевдоевклидовом линейном пространстве.

3. Разработан критерия обучения распознаванию объектов двух классов в произвольном метрическом пространстве как обобщение классического метода опорных векторов.

4. Разработаны численные методы и алгоритмы обучения распознаванию объектов двух классов в произвольном метрическом пространстве.

5. Проведено количественное экспериментальное исследование разработанных методов и алгоритмов на реальных данных.

Разин Н.А. Выпуклые критерии и параллелизуемые алгоритмы селективного комбинирования разнородных представлений объектов в задачах восстановления зависимостей по эмпирическим данным. Дисс. к.ф.-м.н. ВЦ РАН, 2013. SVC 2004. First International Signature Verification Competition. [Электронный ресурс] URL: http://www.cs.ust.hk/svc2004/index.html

27 Mottl V., Lange M., Sulimova V., Ermakov A. Signature verification based on fusion of on-line and off-line kernels. Proceedings of the 19th International Conference on Pattern Recognition. Tampa, USA, December 8-11, 2008, pp. 1-4.

Список публикаций

[1] Абрамов В.И., Середин О.С., Сулимова В.В., Моттль В.В. Эквивалентность потенциальных функций и линейных пространств представления объектов произвольной природы. Доклады 8-й Международной конференции «Интеллектуализация обработки информации», Кипр, Пафос, 17-24 октября 2010 г., с. 9-12. ООО «МАКС Пресс».

[2] Абрамов В.И., Середин О.С., Сулимова В.В., Моттль В.В. Метод опорных объектов для обучения распознаванию образов в евклидовых метрических пространствах. Доклады 9-й международной конференции «Интеллектуализация обработки информации ИОИ-2012», Будва, Черногория, 16-22 сентября 2012 г. М.: Торус Пресс, 2012, с. 5-8.

[3] Абрамов В.И., Середин О.С., Моттль В.В. Обучеиие распознаванию образов в евклидовых метрических пространствах по методу опорных объектов. Известия ТулГУ, Естественные науки. Тула: Изд-во ТулГУ, Вып. 2, 2013, с. 119-136.

[4] Середин О.С., Абрамов В.И., Моттль В.В. Аффинные операции в псевдоевклидовом линейном пространстве. Известия ТулГУ, Естественные науки. Тула: Изд-во ТулГУ, Вып. 3, 2014, с. 178-196.

Жирным шрифтом выделены статьи в журналах, рекомендованы ВАК. Личный вклад соискателя в работах с соавторами заключается в следующем:

[1] Доказано существование континуума линейных пространств со скалярным произведением, порождающих одну и ту же пред-евклидову метрику.

[2,3] Разработан метод обучения распознаванию образов на множествах объектов, в которых напрямую определена лишь пред-евклидова метрика вместо потенциальной функции, существенно более сложной по своей математической структуре.

[4] Разработан математический аппарат обобщения методологии обучения в метрических пространствах с пред-евклидовой метрикой на случай произвольной метрики.

Абрамов Вадим Игоревич

Метод опорных объектов для обучения распознаванию образов в произвольных метрических пространствах

АВТОРЕФЕРАТ

Подписано в печать 12.11.2014. Формат 60 х 84 '/16. Усл. печ. л. 1,0. Тираж 100 экз. Заказ № 433.

Федеральное государственное автономное образовательное учреждение высшего профессионального образования «Московский физико-технический институт (государственный университет)» Отдел оперативной полиграфии «Физтех-полиграф» 141700, Московская обл., г. Долгопрудный, Институтский пер., 9