Машина опорных векторов
Материал из MachineLearning.
м (Разбиение одного блока на несколько блоков для снятия нагрузки с движка, генерирующего изображения (не отображается изображение)) |
|||
(10 промежуточных версий не показаны.) | |||
Строка 2: | Строка 2: | ||
Оптимальная разделяющая гиперплоскость. Понятие зазора между классами (margin). Случай линейной разделимости. Задача квадратичного программирования. Опорные векторы. Случай отсутствия линейной разделимости. Функции ядра (kernel functions), спрямляющее пространство, теорема Мерсера. Способы построения ядер. Примеры ядер. Сопоставление SVM и нейронной RBF-сети. Обучение SVM методом активных ограничений. SVM-регрессия. | Оптимальная разделяющая гиперплоскость. Понятие зазора между классами (margin). Случай линейной разделимости. Задача квадратичного программирования. Опорные векторы. Случай отсутствия линейной разделимости. Функции ядра (kernel functions), спрямляющее пространство, теорема Мерсера. Способы построения ядер. Примеры ядер. Сопоставление SVM и нейронной RBF-сети. Обучение SVM методом активных ограничений. SVM-регрессия. | ||
- | |||
- | |||
== Машина опорных векторов в задачах классификации == | == Машина опорных векторов в задачах классификации == | ||
Строка 11: | Строка 9: | ||
=== Линейно разделимая выборка === | === Линейно разделимая выборка === | ||
- | Рассмотрим задачу нахождения наилучшего в некотором смысле разделения множества векторов на два класса с помощью линейной решающей функции. Пусть имеется множество прецедентов <tex>(\Xi ,Y)</tex>, где <tex>\Xi = \{ {\bf{x}}_1 ,...,{\bf{x}}_N \}</tex> | + | Рассмотрим задачу нахождения наилучшего в некотором смысле разделения множества векторов на два класса с помощью линейной решающей функции. Пусть имеется множество прецедентов <tex>(\Xi ,Y)</tex>, где <tex>\Xi = \{ {\bf{x}}_1 ,...,{\bf{x}}_N \}</tex> — обучающая выборка, а <tex>Y = (y_1 ,...,y_N )</tex> — множество меток двух классов <tex>\omega _1</tex> и <tex>\omega _2</tex>. Требуется по обучающей выборке построить линейную решающую функции, т.е. такую линейную функцию <tex>f({\bf{x}})</tex>, которая удовлетворяла бы условию |
::<tex>f({\bf{x}}_i ) > 0</tex> для всех <tex>{\bf{x}}_i \in \omega _1</tex>, | ::<tex>f({\bf{x}}_i ) > 0</tex> для всех <tex>{\bf{x}}_i \in \omega _1</tex>, | ||
Строка 25: | Строка 23: | ||
::<tex>y_i f({\bf{x}}_i ) > 0 </tex> для всех <tex> {\bf{x}}_i \in \Xi</tex> | ::<tex>y_i f({\bf{x}}_i ) > 0 </tex> для всех <tex> {\bf{x}}_i \in \Xi</tex> | ||
- | Умножая, если нужно функцию <tex>f</tex> на некоторое положительное число нетрудно видеть, что система неравенств {{eqref|1}} равносильна системе | + | Умножая, если нужно, функцию <tex>f</tex> на некоторое положительное число, нетрудно видеть, что система неравенств {{eqref|1}} равносильна системе |
::<tex>y_i f({\bf{x}}_i ) > 1</tex> для всех <tex>{\bf{x}}_i \in \Xi </tex>. | ::<tex>y_i f({\bf{x}}_i ) > 1</tex> для всех <tex>{\bf{x}}_i \in \Xi </tex>. | ||
- | Кроме того, так как <tex>f({\bf{x}})</tex> | + | Кроме того, так как <tex>f({\bf{x}})</tex> — линейная функция, то последняя система неравенств примет вид |
{{eqno|2}} | {{eqno|2}} | ||
::<tex>y_i (({\bf{w}},{\bf{x}}_i ) + b) \ge 1,\quad i = 1,...,N</tex>, | ::<tex>y_i (({\bf{w}},{\bf{x}}_i ) + b) \ge 1,\quad i = 1,...,N</tex>, | ||
- | где <tex>{\bf{w}}</tex> вектор весовых коэффициентов, <tex>b</tex> некоторое число. Тогда разделяющей два класса гиперплоскостью будет <tex>({\bf{w}},{\bf{x}}) + b = 0</tex>. Нетрудно видеть, что и все гиперплоскости вида <tex>({\bf{w}},{\bf{x}}) + b' = 0</tex>, где <tex>b' \in (b - 1,b + 1)</tex> также будут разделяющими (рис.1). Расстояние между граничными гиперплоскостями <tex>({\bf{w}},{\bf{x}}) + b - 1 = 0</tex> и <tex>({\bf{w}},{\bf{x}}) + b + 1 = 0</tex> равно <tex>{{\frac {{2}}{{\left\| {\bf{w}} \right\|}} </tex>. Действительно, <tex>\left( {\frac{{\bf{w}}}{{\left\| {\bf{w}} \right\|}},{\bf{x}}} \right) + \frac{{b - 1}}{{\left\| {\bf{w}} \right\|}} = 0</tex> и <tex>\left( {\frac{{\bf{w}}}{{\left\| {\bf{w}} \right\|}},{\bf{x}}} \right) + \frac{{b + 1}}{{\left\| {\bf{w}} \right\|}} = 0</tex> | + | где <tex>{\bf{w}}</tex> — вектор весовых коэффициентов, <tex>b</tex> — некоторое число. Тогда разделяющей два класса гиперплоскостью будет <tex>({\bf{w}},{\bf{x}}) + b = 0</tex>. Нетрудно видеть, что и все гиперплоскости вида <tex>({\bf{w}},{\bf{x}}) + b' = 0</tex>, где <tex>b' \in (b - 1,b + 1)</tex>, также будут разделяющими (рис.1). Расстояние между граничными гиперплоскостями <tex>({\bf{w}},{\bf{x}}) + b - 1 = 0</tex> и <tex>({\bf{w}},{\bf{x}}) + b + 1 = 0</tex> равно <tex>{{\frac {{2}}{{\left\| {\bf{w}} \right\|}} </tex>. Действительно, <tex>\left( {\frac{{\bf{w}}}{{\left\| {\bf{w}} \right\|}},{\bf{x}}} \right) + \frac{{b - 1}}{{\left\| {\bf{w}} \right\|}} = 0</tex> и <tex>\left( {\frac{{\bf{w}}}{{\left\| {\bf{w}} \right\|}},{\bf{x}}} \right) + \frac{{b + 1}}{{\left\| {\bf{w}} \right\|}} = 0</tex> — нормальные уравнения этих гиперплоскостей. Тогда <tex>p_1 = \frac{{b - 1}}{{\left\| {\bf{w}} \right\|}}</tex> и <tex>p_2 = \frac{{b + 1}}{{\left\| {\bf{w}} \right\|}}</tex> — расстояния от этих гиперплоскостей до начала координат и <tex>{{\frac {{2}}{{\left\| {\bf{w}} \right\|}} </tex> — расстояние между гиперплоскостями. На самих граничных плоскостях может находиться некоторое число обучающих векторов. Эти векторы называются ''опорными''.[[Изображение:pic_l_1.jpg|thumb]] |
- | Для надежного разделения классов необходимо чтобы расстояние между разделяющими гиперплоскостями было как можно большим, т.е. <tex>\left\| {\bf{w}} \right\|</tex> была как можно меньше. Таким образом, ставится задача нахождения минимума квадратичного функционала <tex>0.5({\bf{w}},{\bf{w}})</tex> (коэффициент 0.5 вводится для удобства дифференцирования) в выпуклом многограннике, задаваемым системой неравенств {{eqref|2}}. В выпуклом множестве квадратичный функционал всегда имеет единственный минимум (если это множество не пусто). Из теоремы Куна | + | Для надежного разделения классов необходимо, чтобы расстояние между разделяющими гиперплоскостями было как можно большим, т.е. <tex>\left\| {\bf{w}} \right\|</tex> была как можно меньше. Таким образом, ставится задача нахождения минимума квадратичного функционала <tex>0.5({\bf{w}},{\bf{w}})</tex> (коэффициент 0.5 вводится для удобства дифференцирования) в выпуклом многограннике, задаваемым системой неравенств {{eqref|2}}. В выпуклом множестве квадратичный функционал всегда имеет единственный минимум (если это множество не пусто). Из теоремы Куна — Таккера следует, что решение этой оптимизационной задачи равносильно поиску седловой точки лагранжиана |
::<tex>L({\bf{w}},b,{\bf{\lambda }}) = 0.5({\bf{w}},{\bf{w}}) - \sum\limits_{i = 1}^N {\lambda _i (y_i (({\bf{w}},{\bf{x}}_i ) + b) - 1)} \to \ {\min }\limits_{{\bf{w}},b} \ {\max }\limits_{\bf{\lambda }}</tex> | ::<tex>L({\bf{w}},b,{\bf{\lambda }}) = 0.5({\bf{w}},{\bf{w}}) - \sum\limits_{i = 1}^N {\lambda _i (y_i (({\bf{w}},{\bf{x}}_i ) + b) - 1)} \to \ {\min }\limits_{{\bf{w}},b} \ {\max }\limits_{\bf{\lambda }}</tex> | ||
Строка 112: | Строка 110: | ||
где коэффициенты <tex>\lambda _i </tex> зависят от <tex>y_i </tex> и от значения <tex>\left( {\varphi ({\bf{x}}_i ),\varphi ({\bf{x}}_j )} \right) </tex>. Таким образом, для нахождения решающей функции нужно знать значения скалярных произведений <tex>\left( {\varphi ({\bf{x}}_i ),\varphi ({\bf{x}}_j )} \right) </tex>. Для этого исследуем свойства функции <tex>K({\bf{x}},{\bf{y}}) = \left( {\varphi ({\bf{x}}),\varphi ({\bf{y}})} \right) </tex>, которая называется ядром. Следующая теорема, известная в теории интегральных операторов и доказанная [http://en.wikipedia.org/wiki/James_Mercer_(mathematician) Джеймсом Мерсером] в 1909 году, полностью характеризует ядро. | где коэффициенты <tex>\lambda _i </tex> зависят от <tex>y_i </tex> и от значения <tex>\left( {\varphi ({\bf{x}}_i ),\varphi ({\bf{x}}_j )} \right) </tex>. Таким образом, для нахождения решающей функции нужно знать значения скалярных произведений <tex>\left( {\varphi ({\bf{x}}_i ),\varphi ({\bf{x}}_j )} \right) </tex>. Для этого исследуем свойства функции <tex>K({\bf{x}},{\bf{y}}) = \left( {\varphi ({\bf{x}}),\varphi ({\bf{y}})} \right) </tex>, которая называется ядром. Следующая теорема, известная в теории интегральных операторов и доказанная [http://en.wikipedia.org/wiki/James_Mercer_(mathematician) Джеймсом Мерсером] в 1909 году, полностью характеризует ядро. | ||
- | '''Теорема 1.''' Функция <tex>K({\bf{x}},{\bf{y</tex> является ядром тогда и только тогда, когда она удовлетворяет условиям: | + | '''Теорема 1.''' Функция <tex>K({\bf{x}},{\bf{y}})</tex> является ядром тогда и только тогда, когда она удовлетворяет условиям: |
# <tex>K({\bf{x}},{\bf{y}}) = K({\bf{y}},{\bf{x}})</tex> (симметричность); | # <tex>K({\bf{x}},{\bf{y}}) = K({\bf{y}},{\bf{x}})</tex> (симметричность); | ||
Строка 195: | Строка 193: | ||
== Машина опорных векторов в задачах регрессии == | == Машина опорных векторов в задачах регрессии == | ||
+ | |||
+ | === Постановка задачи === | ||
+ | [[Изображение:Eps-insensitive.jpg|thumb|Синяя - функция потерь <tex>\mathcal{L}(y,g(x))</tex>]] | ||
+ | Ставится задача [[Регрессионный анализ|{{S|регрессии}}]]: | ||
+ | Дана обучающая выборка <tex>X=\{(x_i,y_i)\}_{i=1}^l</tex>, где <tex>x_i</tex>-признаковое описание i-го объекта, <tex>y_i</tex> - характеристика, приписываемая объекту. Целью является поиск такой функции <tex>g(\mathbf{x})</tex>, которая бы аппроксимировала выборку наилучшим образом. Задана функция потерь <tex> \mathcal{L}(y,g(\mathbf{x}))= | ||
+ | \begin{cases} | ||
+ | 0, & \mbox{ } \mid y-g(x)\mid \le \epsilon\\ | ||
+ | \mid y-g(x)\mid - \epsilon, & \mbox{ } \mid y-g(x)\mid > \epsilon | ||
+ | \end{cases} </tex>, | ||
+ | где <tex>\epsilon>0</tex> определяет допустимое отклонение для результата. Считается, что этот параметр задает эксперт. | ||
+ | |||
+ | === Анализ задачи === | ||
+ | Будем искать решение задачи регрессии в линейном случае: | ||
+ | <tex>f(x) = (w,\mathbf{x})-w_0 </tex>. | ||
+ | Функция потерь принимает вид <tex>a(x_i)=\mid (w,x_i)-w_0-y_i \mid_\epsilon</tex> для каждого вектора <tex>(x_i,y_i)</tex>, где <tex>\mid z \mid_\epsilon = max(0,\mid z \mid-\epsilon)</tex>. | ||
+ | |||
+ | В таком случае функционал потерь принимает вид | ||
+ | ::<tex>Q_\epsilon(a,X)=\sum_{i=1}^l \mid (w,x_i)-w_0-y_i \mid_\epsilon + \tau (w,w)^2 \rightarrow \underset{w,w_0}{min}</tex>. | ||
+ | |||
+ | Последнее слагаемое удерживает коэффициенты <tex>w</tex> от бесконечного возрастания. | ||
+ | Как и в задаче классификации, решение зависит от скалярного произведения объектов, а не от самих объектов. Минимизация в данном случае эквивалентна задаче квадаратичного программирования с ограничениями типа неравенств. Покажем это. Положим <tex>C=\frac{1}{2\tau}</tex>, и введем дополнительные переменные <tex>\xi_i^+</tex> и <tex>\xi_i^-</tex>, значения которых равны потере при завышенном и заниженном ответах соответственно | ||
+ | |||
+ | ::<tex>\xi_i^+=(a(x_i)-y_i-\epsilon)_+</tex>, <tex>\xi_i^-=(-a(x_i)+y_i-\epsilon)_-</tex>, <tex>i=1,...,l</tex>. | ||
+ | Теперь мы можем записать задачу минимизации в виде задачи квадратичного программирования. | ||
+ | |||
+ | <tex> | ||
+ | \begin{cases} | ||
+ | \frac{1}{2} (w,w)^2 + C\sum_{i=1}^l(\xi_i^+ + \xi_i^-)\rightarrow \underset{w,w_0,\xi_i^+,\xi_i^-}{min}, \\ | ||
+ | (w,x_i)-w_0 \le y_i + \epsilon + \xi_i^+, & i=1,..,l; \\ | ||
+ | (w,x_i)-w_0 \ge y_i - \epsilon - \xi_i^-, & i=1,..,l; \\ | ||
+ | \xi_i^- \ge 0, \mbox{ } i=1,..,l; \\ | ||
+ | \xi_i^+ \ge 0, \mbox{ } i=1,..,l; \\ | ||
+ | \end{cases} | ||
+ | </tex> | ||
+ | |||
+ | === Опорные векторы и двойственная задача === | ||
+ | Как и в задаче классификации решается двойственная задача. Лагранжиан задается через двойственные переменные <tex> | ||
+ | \lambda_i^+, \lambda_i^-, i=1,...,l | ||
+ | </tex>, | ||
+ | а скалярные произведения можно заменить ядром <tex>K(x_i,x_j)</tex> | ||
+ | |||
+ | Получим результат: | ||
+ | |||
+ | <tex> | ||
+ | \begin{cases} | ||
+ | \mathcal{L}(\lambda^+, \lambda^-)=-\epsilon \sum_{i=1}^l (\lambda_i^+ + \lambda_i^-)+\sum_{i=1}^l (\lambda_i^- - \lambda_i^+)y_i - \frac{1}{2}\sum_{i,j=1}^l (\lambda_i^- - \lambda_i^+)(\lambda_j^- - \lambda_j^+)K(x_i,x_j) \rightarrow \underset{\lambda^+,\lambda^-}{max}, \\ | ||
+ | 0 \le \lambda_i^+ \le C,\mbox{ } i=1,...,l \\ | ||
+ | 0 \le \lambda_i^- \le C,\mbox{ } i=1,...,l \\ | ||
+ | \sum_{i=1}^l (\lambda_i^+ + \lambda_i^-) = 0;\\ | ||
+ | \end{cases} | ||
+ | </tex> | ||
+ | |||
+ | Таким образом, все объекты <tex>x_i, i=1,...,l</tex> делятся на 5 типов: | ||
+ | |||
+ | <tex>1. \mid a(x_i)-y_i \mid < \epsilon;\mbox{ } \lambda_i^+=\lambda_i^-=\xi_i^+=\xi_i^-=0;</tex><br> | ||
+ | <tex>2. a(x_i)=y_i + \epsilon;\mbox{ }0 < \lambda_i^+ < C,\mbox{ } \lambda_i^-=\xi_i^+=\xi_i^-=0;\\</tex><br> | ||
+ | <tex>3. a(x_i)=y_i + \epsilon;\mbox{ }0 < \lambda_i^- < C,\mbox{ } \lambda_i^+=\xi_i^+=\xi_i^-=0;\\</tex><br> | ||
+ | <tex>4. a(x_i)>y_i + \epsilon;\mbox{ }\lambda_i^+ = C,\mbox{ } \lambda_i^-=0,\mbox{ }\xi_i^+=a(x_i) - y_i - \epsilon,\mbox{ } \xi_i^-=0;\\</tex><br> | ||
+ | <tex>5. a(x_i)<y_i - \epsilon;\mbox{ }\lambda_i^+ = 0,\mbox{ } \lambda_i^-=C,\mbox{ }\xi_i^+=0,\mbox{ } \xi_i^-=y_i - a(x_i) - \epsilon;\\</tex><br> | ||
+ | |||
+ | |||
+ | Объекты 2-5 называются опорными и учитываются при определении вектора весов. На объектах 1-3 потери равны нулю, а на 4 и 5 потери больше нуля. | ||
+ | Уравнение регрессии выражается через двойственные переменные: | ||
+ | <tex> | ||
+ | a(x)=\sum_{i=1}^l (\lambda_i^- - \lambda_i^+)K(x_i,x)-w0; | ||
+ | </tex>, где параметр <tex>w_0</tex> определяется из уравнений | ||
+ | |||
+ | <tex> | ||
+ | (w,x_i)-w_0 = | ||
+ | \begin{cases} | ||
+ | y+\epsilon, & \mbox{ } \x_i \mbox{ belongs to 2nd type }\\ | ||
+ | y-\epsilon, & \mbox{ } \x_i \mbox{ belongs to 3rd type }\\ | ||
+ | \end{cases} | ||
+ | </tex> | ||
+ | |||
+ | Параметр <tex>w_0</tex> можно определить из любого уравнения, но чтобы избежать вычислительных ошибок берется медиана полученного множества для <tex>w_0</tex>. | ||
+ | |||
+ | Как уже упоминалось, параметр <tex>\epsilon</tex> задает эксперт, а параметр С подбирается по [[Скользящий контроль|{{S|скользящему контролю}}]], что является трудоёмкой процедурой. | ||
+ | |||
+ | == Мультиклассовый метод опорных векторов == | ||
+ | |||
+ | В случае нескольких классов на практике зачастую применяется решающее правило, основанное на разбиении задачи на бинарные по схеме "один против остальных" (One-vs-Rest): | ||
+ | |||
+ | <tex> | ||
+ | f(x) = \arg \max_{y\in Y} (\langle w_y, x \rangle + b_y). | ||
+ | </tex> | ||
+ | |||
+ | Стоит отметить, что дальнейшее обобщение такого решающего правила приводит к методу опорных векторов со структурированным результатом (Structured Output SVM). | ||
+ | |||
+ | Для обучения параметров <tex>w_0, b_0, \dots, w_K, b_K</tex> такой композиции алгоритмов может быть решено <tex>K</tex> стандартных для метода опорных векторов задач оптимизации, однако данный подход имеет свои недостатки. В работе Краммера и Зингера была предложена следующая совместная формулировка задачи оптимизации | ||
+ | |||
+ | <tex> | ||
+ | \min_{w,\xi} \frac{1}{2} ||w_y||^2 + C \sum_{i} \xi_i, | ||
+ | </tex> | ||
+ | |||
+ | при ограничениях <tex>\langle w_{y_i}, x_i \rangle - \langle w_{y}, x_i \rangle \geq 1 - \delta_{y_i, y} - \xi_i</tex>, задающих ответы классификатора на обучающей выборке и ограничениях положительности штрафов <tex>\xi_i \geq 0 </tex>. | ||
+ | |||
+ | Двойственная задача такой формулировки имеет вид | ||
+ | |||
+ | <tex> | ||
+ | \min_{\alpha} \frac{1}{2} \sum_{y\in Y} \sum_{i} \sum_{j} \alpha_i^y \alpha_j^y \langle x_i, x_j \rangle + \sum_{y \in Y} \sum_{i} (1-\delta_{y_i, y} ) \alpha_i^y | ||
+ | </tex> | ||
+ | |||
+ | при неотрицательных коэффициентах <tex>\alpha_i^y \geq 0</tex>. | ||
== Программные реализации == | == Программные реализации == | ||
Строка 200: | Строка 302: | ||
*[http://svmlight.joachims.org/ SVMLight] | *[http://svmlight.joachims.org/ SVMLight] | ||
- | Наиболее развитая и популярная реализация SVM на С++. Библиотека адаптирована для больших выборок и имеет эффективную реализацию скользящего контроля. Включены стандартные | + | Наиболее развитая и популярная реализация SVM на С++. Библиотека адаптирована для больших выборок и имеет эффективную реализацию скользящего контроля. Включены стандартные ядерные функции и допускается использование предварительно вычисленных матриц ядерных функций. |
*[http://www.csie.ntu.edu.tw/~cjlin/libsvm/ LIBSVM] и [http://www.csie.ntu.edu.tw/~cjlin/liblinear/ LIBLINEAR] | *[http://www.csie.ntu.edu.tw/~cjlin/libsvm/ LIBSVM] и [http://www.csie.ntu.edu.tw/~cjlin/liblinear/ LIBLINEAR] | ||
Эффективные и простые в использовании реализации SVM на С++ со схожими интерфейсами. | Эффективные и простые в использовании реализации SVM на С++ со схожими интерфейсами. | ||
- | В | + | В LibLinear реализована только линейная классификация и регрессия. Работают с большими выборками. |
== Литература == | == Литература == | ||
Строка 216: | Строка 318: | ||
* [[Линейный классификатор]] | * [[Линейный классификатор]] | ||
* [[Машинное обучение (курс лекций, К.В.Воронцов)]] | * [[Машинное обучение (курс лекций, К.В.Воронцов)]] | ||
+ | * <i>Alexey Nefedov</i>. [https://svmtutorial.online Support Vector Machines: A Simple Tutorial] | ||
[[Категория:Классификация]] | [[Категория:Классификация]] | ||
[[Категория:Энциклопедия анализа данных]] | [[Категория:Энциклопедия анализа данных]] | ||
[[Категория:Линейные классификаторы]] | [[Категория:Линейные классификаторы]] |
Текущая версия
Машина опорных векторов — является одной из наиболее популярных методологий обучения по прецедентам, предложенной В. Н. Вапником и известной в англоязычной литературе под названием SVM (Support Vector Machine).
Оптимальная разделяющая гиперплоскость. Понятие зазора между классами (margin). Случай линейной разделимости. Задача квадратичного программирования. Опорные векторы. Случай отсутствия линейной разделимости. Функции ядра (kernel functions), спрямляющее пространство, теорема Мерсера. Способы построения ядер. Примеры ядер. Сопоставление SVM и нейронной RBF-сети. Обучение SVM методом активных ограничений. SVM-регрессия.
Содержание |
Машина опорных векторов в задачах классификации
Понятие оптимальной разделяющей гиперплоскости
Линейно разделимая выборка
Рассмотрим задачу нахождения наилучшего в некотором смысле разделения множества векторов на два класса с помощью линейной решающей функции. Пусть имеется множество прецедентов , где
— обучающая выборка, а
— множество меток двух классов
и
. Требуется по обучающей выборке построить линейную решающую функции, т.е. такую линейную функцию
, которая удовлетворяла бы условию
для всех
,
для всех
.
Без ограничения общности можно считать, что метки классов равны
Тогда поставленную выше задачу можно переформулировать следующим образом. Требуется найти линейную решающую функцию , которая бы удовлетворяла условию
для всех
Умножая, если нужно, функцию на некоторое положительное число, нетрудно видеть, что система неравенств (1) равносильна системе
для всех
.
Кроме того, так как — линейная функция, то последняя система неравенств примет вид
,
Для надежного разделения классов необходимо, чтобы расстояние между разделяющими гиперплоскостями было как можно большим, т.е. была как можно меньше. Таким образом, ставится задача нахождения минимума квадратичного функционала
(коэффициент 0.5 вводится для удобства дифференцирования) в выпуклом многограннике, задаваемым системой неравенств (2). В выпуклом множестве квадратичный функционал всегда имеет единственный минимум (если это множество не пусто). Из теоремы Куна — Таккера следует, что решение этой оптимизационной задачи равносильно поиску седловой точки лагранжиана
в ортанте по множителям Лагранжа , при условии, что
.
Последнее условие равносильно тому, что
или
Из необходимых условий существования седловой точки (полагая ) имеем
Откуда следует, что вектор следует искать в виде
,
причем
.
В силу (3) в сумму (4) с ненулевыми коэффициентами входят только те векторы, для которых
. Такие векторы называют опорными, так как это именно те векторы, через которые будут проходить граничные гиперплоскости, разделяющие классы. Для найденного весового вектора
смещение
можно вычислить как
для любого опорного вектора
.
Найдем значения множителей Лагранжа, как критических точек лагранжиана. Для этого подставим (4) и (5) в лагранжиан, получим
.
Таким образом, задача сводится к нахождению критических точек функции
.
Так как эта функция представляет собой разность линейной и квадратичной функций, причем квадратичная функция отрицательно определена, то требуется найти наибольшее значение функции при условии
в области
. Существует много алгоритмов (в теории оптимизации) решения этой задачи (например, градиентные методы, метод покоординатного спуска и т.д.).
- Замечания.
- Суммирования в (6) осуществляются не по всем векторам, а только по опорным, которых может быть гораздо меньше, чем обучающих.
- Линейная решающая функция в результате имеет вид
, где
зависят только
и от значений скалярного произведения
, причем суммирования осуществляются только по опорным векторам.
- После того, как решающая функция
вычислена, вектор
следует относить классу
, если
и классу
, если
. Вероятность неправильной классификации можно оценить с помощью некоторой непрерывно убывающей функции
, удовлетворяющей условиям:
,
при
. Тогда вероятность
неправильной классификации вектора
будет равна
, если
(
), где
:
,
. То есть
, если
.
- В такой постановке алгоритм линейный классификации был разработан В. Вапником в 1963 году.
- Пример. Методом опорных векторов разделите классы
и
, если
,
,
.
- Решение. Положим
,
,
. Тогда функция
будет иметь вид
,
причем . Тогда
. Составим и решим нормальную систему для функции
:
Если же , то
и
,
при
. Следовательно,
,
и
.
Таким образом, наибольшее значение функции в области
при условии
достигается в точке
. В этом случае,
Таким образом, и
. Ширина разделяющей полосы будет равна
, а прямые
и
будут ее границами (см. рис.2).
Линейно неразделимая выборка
В 1992 году в работе Бернарда Бозера (Boser B.), Изабелл Гийон (Guyon I.) и Владимира Вапника был предложен способ адаптации машины опорных векторов для нелинейного разделения классов. В этом случае нужно вложить пространство признаков в пространство
большей размерности с помощью отображения
. Будем считать, что
пространство со скалярным произведением. Тогда, рассматривая алгоритм опорных векторов для образов
обучающей выборки, сведем решение задачи к линейно разделимому случаю, т.е. разделяющую функцию будем искать в виде
,
,
где коэффициенты зависят от
и от значения
. Таким образом, для нахождения решающей функции нужно знать значения скалярных произведений
. Для этого исследуем свойства функции
, которая называется ядром. Следующая теорема, известная в теории интегральных операторов и доказанная Джеймсом Мерсером в 1909 году, полностью характеризует ядро.
Теорема 1. Функция является ядром тогда и только тогда, когда она удовлетворяет условиям:
-
(симметричность);
-
неотрицательно определена, т.е. матрица
,
является неотрицательно определенной для любых векторов
.
Упражнение. Докажите, что следующие функции являются ядрами:
-
;
-
;
-
.
Теорема 2. Справедливы следующие свойства ядер:
- сумма ядер – ядро;
- произведение ядер – ядро;
- сумма равномерно сходящегося ряда ядер – ядро;
- композиция ядра и любого отображения (т.е.
) – ядро.
Следствие.
- многочлен с положительными коэффициентами от ядра – ядро;
- экспонента от ядра – ядро;
- функция
ядро.
Доказательство. Утверждения 1 и 2 следуют из пунктов 1, 2 и 3 теоремы. Справедливость утверждения 3 вытекает из того, что , а симметричность и положительная определенность функций
проверяется непосредственно.
Любые векторов могут быть разделены на любые два класса с помощью мономиального отображения степени не больше
. Поэтому, если
,
такое отображение, то ядро, соответствующее этому отображению можно искать в виде
.
Таким образом, это ядро гарантирует разделение любых векторов на любые два класса. В этом случае нахождение разделяющих функций осуществляется следующим образом:
1) найдем наибольшее значение функции
при условии в области
, получим вектор
;
2) разделяющую функцию ищем в виде
.
Преимущества и недостатки SVM:
- это наиболее быстрый метод нахождения решающих функций;
- метод сводится к решению задачи квадратичного программирования в выпуклой области, которая всегда имеет единственное решение;
- метод находит разделяющую полосу максимальной ширины, что позволяет в дальнейшем осуществлять более уверенную классификацию;
- метод чувствителен к шумам и стандартизации данных;
- не существует общего подхода к автоматическому выбору ядра (и построению спрямляющего подпространства в целом) в случае линейной неразделимости классов.
Решение. Так как количество векторов обучающей выборки равно трем, то и в качестве ядра возьмем функцию
. Тогда
.
Составим нормальную систему
Можно показать, что в точке будет достигаться наибольшее значение функции
при условии
в области
. Следовательно,
.
Ядра и спрямляющие пространства
Наиболее распространенные ядра:
- Полиномиальное:
- Полиномиальное:
- Радиальное:
, для
Алгоритмы настройки
Машина опорных векторов в задачах регрессии
Постановка задачи
Ставится задача регрессии:
Дана обучающая выборка , где
-признаковое описание i-го объекта,
- характеристика, приписываемая объекту. Целью является поиск такой функции
, которая бы аппроксимировала выборку наилучшим образом. Задана функция потерь
,
где
определяет допустимое отклонение для результата. Считается, что этот параметр задает эксперт.
Анализ задачи
Будем искать решение задачи регрессии в линейном случае:
.
Функция потерь принимает вид
для каждого вектора
, где
.
В таком случае функционал потерь принимает вид
.
Последнее слагаемое удерживает коэффициенты от бесконечного возрастания.
Как и в задаче классификации, решение зависит от скалярного произведения объектов, а не от самих объектов. Минимизация в данном случае эквивалентна задаче квадаратичного программирования с ограничениями типа неравенств. Покажем это. Положим
, и введем дополнительные переменные
и
, значения которых равны потере при завышенном и заниженном ответах соответственно
,
,
.
Теперь мы можем записать задачу минимизации в виде задачи квадратичного программирования.
Опорные векторы и двойственная задача
Как и в задаче классификации решается двойственная задача. Лагранжиан задается через двойственные переменные ,
а скалярные произведения можно заменить ядром
Получим результат:
Таким образом, все объекты делятся на 5 типов:
Объекты 2-5 называются опорными и учитываются при определении вектора весов. На объектах 1-3 потери равны нулю, а на 4 и 5 потери больше нуля.
Уравнение регрессии выражается через двойственные переменные:
, где параметр
определяется из уравнений
Параметр можно определить из любого уравнения, но чтобы избежать вычислительных ошибок берется медиана полученного множества для
.
Как уже упоминалось, параметр задает эксперт, а параметр С подбирается по скользящему контролю, что является трудоёмкой процедурой.
Мультиклассовый метод опорных векторов
В случае нескольких классов на практике зачастую применяется решающее правило, основанное на разбиении задачи на бинарные по схеме "один против остальных" (One-vs-Rest):
Стоит отметить, что дальнейшее обобщение такого решающего правила приводит к методу опорных векторов со структурированным результатом (Structured Output SVM).
Для обучения параметров такой композиции алгоритмов может быть решено
стандартных для метода опорных векторов задач оптимизации, однако данный подход имеет свои недостатки. В работе Краммера и Зингера была предложена следующая совместная формулировка задачи оптимизации
при ограничениях , задающих ответы классификатора на обучающей выборке и ограничениях положительности штрафов
.
Двойственная задача такой формулировки имеет вид
при неотрицательных коэффициентах .
Программные реализации
Наиболее развитая и популярная реализация SVM на С++. Библиотека адаптирована для больших выборок и имеет эффективную реализацию скользящего контроля. Включены стандартные ядерные функции и допускается использование предварительно вычисленных матриц ядерных функций.
Эффективные и простые в использовании реализации SVM на С++ со схожими интерфейсами. В LibLinear реализована только линейная классификация и регрессия. Работают с большими выборками.
Литература
- Вапник В. Н. Восстановление зависимостей по эмпирическим данным. — М.: Наука, 1979. — 448 с. (подробнее)
- Hastie, T., Tibshirani, R., Friedman, J. The Elements of Statistical Learning, 2nd edition. — Springer, 2009. — 533 p. (подробнее)
Ссылки
- SMO — Sequential Minimal Optimization — эффективный алгоритм настройки SVM.
- Джон Платт — сайт автора SMO: статьи, реализации и тестовые примеры.
- Линейный классификатор
- Машинное обучение (курс лекций, К.В.Воронцов)
- Alexey Nefedov. Support Vector Machines: A Simple Tutorial