Метод опорных векторов (Support vector machine)

Линейно разделимый случай SVM

Пусть дана выборка $(x_{i},y_{i}),~i=1\ldots N,~x_{i}\in R^{D},~y_{i}\in \{-1;1\}$

Линейный классификатор:

$a(x)=sign(w^{T}x+w_{0}),~w\in R^{D},~w_{0}\in R$

Какой геометрический смысл у линейного классификатора? Он строит гиперплоскость, которая задается уравнением:

$w^{T}x+w_{0}=0$

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

$w^{T}x+w_{0}=b,~~~b>0$

Найдем расстояние от него до гиперплоскости. Известно, что расстояние может быть посчитано по формуле: ${\frac {|w^{T}x+w_{0}|}{\|w\|}}$ В нашем случае расстояние до гиперплоскости равно ${\frac {b}{\|w\|}}$.

Таким образом ширина разделяющей полосы будет равна ${\frac {2b}{\|w\|}}$ и наша задача заключается в нахождении вектора $w$. Он должен быть таким, чтобы гиперплоскость линейно разделяла нашу выборку и ширина зазора между классами была максимальной. Зачем это нужно? Мы знаем, что чем больше значение $\langle w,x\rangle$тем классификатор увереннее в ответе. Поэтому мы и хотим увеличить ширину зазора между классами, чтобы повысить качество классификации.

Итак, перейдем к математической постановке задачи:

${\begin{cases}{\frac {2b}{\|w\|}}\rightarrow {\underset {w,w_{0}}{\max }}\\\forall x_{i},~y=1~~w^{T}x_{i}+w_{0}\geq b\\\forall x_{i},~y=-1~~w^{T}x_{i}+w_{0}\leq -b\end{cases}}$

Заметим, что второе и третье уравнение можно переписать в виде одного:

$y_{i}\left(w^{T}x_{i}+w_{0}\right)\geq b$

Если точка $(w,w_{0},b)$ - решение задачи, то для $\forall \alpha >0~~(\alpha w,~\alpha w_{0},~\alpha b)$ - тоже решение. Положим $b=1$, т.к. мы можем масштабировать $(w,w_{0},b)$.

Итого получим следующую задачу:

${\begin{cases}\|w\|\rightarrow {\underset {w,w_{0}}{\min }}\\y_{i}\left(w^{T}x_{i}+w_{0}\right)\geq 1~~~\forall i\in 1,\dots ,N\\\end{cases}}$

Выделяют два типа объектов:

• Опорные вектора: $y_{i}\left(w^{T}x_{i}+w_{0}\right)=1$
• Неинформативные векторы: $y_{i}\left(w^{T}x_{i}+w_{0}\right)>1$

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

Линейно неразделимый случай SVM

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

${\begin{cases}\|w\|\rightarrow {\underset {w,w_{0}}{\min }}\\y_{i}\left(w^{T}x_{i}+w_{0}\right)\geq 1-\xi _{i}~~~\forall i\in 1,\dots ,N\\\xi _{i}\geq 0,~~i=1,\dots ,N\end{cases}}$

$\xi$ надо как-то штрафовать, чтобы они не были слишком большими. В итоге получим следующую задачу:

${\begin{cases}{\frac {1}{2}}\|w\|^{2}+C\sum _{i=1}^{N}\xi _{i}\rightarrow {\underset {w,w_{0},\xi }{\min }}\\y_{i}\left(w^{T}x_{i}+w_{0}\right)\geq 1-\xi _{i}~~~\forall i\in 1,\dots ,N\\\xi _{i}\geq 0,~~i=1,\dots ,N\end{cases}}$

Параметр $C$ определяет цену ошибки классификации. Заметим, что в качестве штрафа мы также могли взять $C\sum _{i}^{N}\xi _{i}^{2}$.

Классификация типов объектов

Выделяют несколько типов объектов:

• Неинформативные: $y_{i}\left(w^{T}x_{i}+w_{0}\right)>1$
• Опорные векторы: $y_{i}\left(w^{T}x_{i}+w_{0}\right)\leq 1$
  • Граничные опорные векторы: $y_{i}\left(w^{T}x_{i}+w_{0}\right)=1$
  • Опорные вектора-нарушители:
    • $y_{i}\left(w^{T}x_{i}+w_{0}\right)>0$ - векторы, лежащие внутри разделяющей полосы, но в своем классе
    • $y_{i}\left(w^{T}x_{i}+w_{0}\right)<0$ - векторы, которые классифицируются неверно

Какой функции потерь и регуляризатору соответствует SVM?

Рассмотрим еще раз линейно неразделимый случай SVM.

${\begin{cases}{\frac {1}{2}}\|w\|^{2}+C\sum _{i=1}^{N}\xi _{i}\rightarrow {\underset {w,w_{0},\xi }{\min }}\\y_{i}\left(w^{T}x_{i}+w_{0}\right)\geq 1-\xi _{i}~~~\forall i\in 1,\dots ,N\\\xi _{i}\geq 0,~~i=1,\dots ,N\end{cases}}$

Постараемся избавиться от $\xi$. Ясно, что из последних двух неравенства можем получить:${\begin{cases}\xi _{i}\geq 1-y_{i}\left(w^{T}x_{i}+w_{0}\right)\\\xi _{i}\geq 0,~~i=1,\dots ,N\end{cases}}$

Или иначе это можно записать: $\xi _{i}=\max \left(1-y_{i}\left(w^{T}x_{i}+w_{0}\right),0\right)\geq 0$

И тогда нашу задачу можно переписать в виде:

${\frac {1}{2}}\|w\|^{2}+C\sum _{i=1}^{N}\max \left(1-y_{i}\left(w^{T}x_{i}+w_{0}\right),0\right)\rightarrow {\underset {w,w_{0}}{\min }}$

Это и есть функция потерь, которую мы пытаемся минимизировать. Видно, что в данном случае используется $L_{2}$ регуляризация, а функция потерь является кусочно - линейной.

Вывод решения

Запишем еще раз условия задачи:

${\begin{cases}{\frac {1}{2}}\|w\|^{2}+C\sum _{i=1}^{N}\xi _{i}\rightarrow {\underset {w,w_{0},\xi }{\min }}\\y_{i}\left(w^{T}x_{i}+w_{0}\right)\geq 1-\xi _{i}~~~\forall i\in 1,\dots ,N\\\xi _{i}\geq 0,~~i=1,\dots ,N\end{cases}}$

Запишем Лагранжиан:

$L(w,b,\xi ,\lambda ,\mu )={\frac {1}{2}}\|w\|^{2}+C\sum _{i=1}^{N}\xi _{i}-\sum _{i=1}^{N}\lambda _{i}[y_{i}(w^{T}x_{i}+b)-1+\xi _{i}]-\sum _{i=1}^{N}\mu _{i}\xi _{i}$

Выпишем Условия Куна-Таккера:

${\begin{cases}\bigtriangledown _{w}L=w-\sum _{i=1}^{l}\lambda _{i}y_{i}x_{i}=0\Longrightarrow w=\sum _{i=1}^{l}\lambda _{i}y_{i}x_{i}\\\bigtriangledown _{b}L=-\sum _{i=1}^{l}\lambda _{i}y_{i}=0\Longrightarrow \sum _{i=1}^{l}\lambda _{i}y_{i}=0\\\bigtriangledown _{\xi _{i}}L=C-\lambda _{i}-\mu _{i}\Longrightarrow \mu _{i}+\lambda _{i}=C\\\lambda _{i}[y_{i}(w^{T}x_{i}+b)-1+\xi _{i}]=0\Longrightarrow (\lambda _{i}=0)\ or\ (y_{i}(w^{T}x_{i}+b)=1-\xi _{i})\\\mu _{i}\xi _{i}=0\Longrightarrow (\mu _{i}=0)\ or\ (\xi _{i}=0)\\\xi _{i}\geq 0,\ \lambda _{i}\geq 0,\ \mu _{i}\geq 0\end{cases}}$
Воспользуемся этими условиями и перепишем Лагранжиан:

$L={\frac {1}{2}}\|\sum _{i=1}^{l}\lambda _{i}y_{i}x_{i}\|^{2}-\sum _{i,j=1}^{l}\lambda _{i}\lambda _{j}y_{i}y_{j}x_{i}^{T}x_{j}-b\sum _{i=1}^{l}\lambda _{i}y_{i}+\sum _{i=1}^{l}\lambda _{i}+\sum _{i=1}^{l}\xi _{i}(C-\lambda _{i}-\mu _{i})=$ $=\sum _{i=1}^{l}\lambda _{i}-{\frac {1}{2}}\sum _{i,j=1}^{l}\lambda _{i}\lambda _{j}y_{i}y_{j}x_{i}^{T}x_{j}$

Таким образом мы перешли к двойственной задаче:

${\begin{cases}\sum _{i=1}^{l}\lambda _{i}-{\frac {1}{2}}\sum _{i,j=1}^{l}\lambda _{i}\lambda _{j}y_{i}y_{j}x_{i}^{T}x_{j}\rightarrow {\underset {\lambda }{max}}\\0\leq \lambda _{i}\leq \ C,\ i=1,2,...,l\\\sum _{i=1}^{l}\lambda _{i}y_{i}=0\\\end{cases}}$

Обобщение через ядра

Вместо скалярного произведения можно использовать ядра: ${\begin{cases}\sum _{i=1}^{l}\lambda _{i}-{\frac {1}{2}}\sum _{i,j=1}^{l}\lambda _{i}\lambda _{j}y_{i}y_{j}K(x_{i},x_{j})\rightarrow {\underset {\lambda }{max}}\\0\leq \lambda _{i}\leq \ C,\ i=1,2,...,l\\\sum _{i=1}^{l}\lambda _{i}y_{i}=0\\\end{cases}}$

В этом случае классификатор будет иметь вид:

$a(x)=sign(\sum _{i=1}^{l}\lambda _{i}y_{i}K(x_{i},x)+b)$