机器学习入门(十七):SVM——非线性 SVM 和核函数

非线性分类问题

遇到分类问题的时候,最理想的状态,当然是样本在向量空间中都是线性可分的,我们可以清晰无误地把它们分隔成不同的类别——线性可分 SVM

如果实在不行,我们可以容忍少数不能被正确划分,只要大多数线性可分就好——线性 SVM

可是,如果我们面对的分类问题,根本就是非线性的呢?比如像下面这样:

image

图中红色的点是正类样本,蓝色的点是负类样本。通过我们的观察可知,它们之间的界限是很分明的,用图中绿色的圈本来可以把它们完全分开。

很可惜,“圆圈”在二维空间里无法用线性函数表示,也就是说这些样本在二维空间里根本线性不可分。所以,无论是线性可分 SVM 还是线性 SVM,都无法在这些样本上良好工作。

这可怎么办呢?难道,这种情况我们就处理不了了?

并不是!

我们可以想个办法,让这些在二维空间中线性不可分的样本,在更高维度的空间里线性可分

比如说,如果我们能把上图中那些正负类的样本映射到三维空间中,并且依据不同的类别给它们赋予不一样的高度值——z 轴取值(就像下图这样),那么不就线性可分了嘛。

enter image description here

如此一来,在二维空间团团转的正负例,在三维空间中分为两层,中间用一个超平面,就可以完美分隔了。

非线性 SVM

非线性 SVM 分隔超平面

对于在有限维度向量空间中线性不可分的样本,我们将其映射到更高维度的向量空间里,再通过间隔最大化的方式,学习得到支持向量机,就是非线性 SVM

我们将样本映射到的这个更高维度的新空间叫做特征空间

注意:如果是理想状态,样本从原始空间映射到特征空间后直接就成为线性可分的,那么接下来的学习是可以通过硬间隔最大化的方式来学的。

不过,一般的情况总没有那么理想,因此,通常情况下,我们还是按照软间隔最大化,在特征空间中学习 SVM。

简单理解就是:非线性 SVM = 核技巧 + 线性 SVM

我们用向量 x 表示位于原始空间中的样本,ϕ(x) 表示 x 映射到特征空间之后的新向量。

非线性 SVM对应的分隔超平面为:f(x)=wϕ(x)+b。

非线性 SVM 的对偶问题

套用上一篇线性 SVM 的对偶问题,此处非线性 SVM 的对偶问题就变成了:

min(12∑mi=1∑mj=1αiαjyiyjϕ(xi)⋅ϕ(xj)−∑mi=1αi)

s.t.∑mi=1αiyi=0

0⩽αi⩽C,i=1,2,...,m

大家可以看到,和线性 SVM 唯一的不同就是:之前的 xi 与 xj 的内积(点乘) 变成了 ϕ(xi) 与 ϕ(xj) 的内积

核函数

对于有限维的原始空间,一定存在更高维度的空间,使得前者中的样本映射到新空间后可分。但是新空间(特征空间)的维度也许很大,甚至可能是无限维的。这样的话,直接计算 ϕ(xi)⋅ϕ(xj) 就会很困难。

为了避免计算 ϕ(xi) 和 ϕ(xj) 的内积,我们需要设置一个新的函数——k(xi,xj):

k(xi,xj)=ϕ(xi)⋅ϕ(xj)

原始空间中的两个样本 xi 和 xj 经过 k(⋅,⋅) 函数计算所得出的结果,是它们在特征空间中映射成的新向量的内积。

如此一来,我们就不必真的计算出 ϕ(xi) 点乘 ϕ(xj) 的结果,而可以直接用 k(⋅,⋅) 函数代替它们。

我们把这个 k(⋅,⋅) 函数叫做核函数。现在我们给出它的正式定义:

设 X 为原始空间(又称输入空间),H 为特征空间(特征空间是一个带有内积的完备向量空间,又称完备内积空间或希尔伯特空间)。

如果存在一个映射: X×X ,使得对所有 xi,xj∈X,函数 k(xi,xj) 满足条件:k(xi,xj)=ϕ(xi)⋅ϕ(xj),即 k(⋅,⋅) 函数为输入空间中任意两个向量映射到特征空间后的内积。则称 k(⋅,⋅) 为核函数,ϕ(⋅) 为映射函数。

运用核技巧求解非线性 SVM 的对偶问题

有了核函数,我们就可以将非线性 SVM 的对偶问题写成:

min(12∑mi=1∑mj=1αiαjyiyjk(xi,xj)−∑mi=1αi)

s.t.∑mi=1αiyi=0

0⩽αi⩽C,i=1,2,...,m

之后的求解过程与线性 SVM 一致:先根据对偶问题求解 α,再根据 α 的结果计算 w,然后根据支持向量求解 b,在此就不赘述了。

相应地,最终求出的非线性 SVM 在特征空间的最大分隔超平面也就成了:

f(x)=wϕ(x)+b=∑mi=1αiyiϕ(xi)⋅ϕ(x)+b=∑mi=1αiyik(x,xi)+b

上述运用核函数求解的过程,称为核技巧。

核函数的性质

如果我们已经知道了映射函数是什么,当然可以通过两个向量映射后的内积直接求得核函数。

但是,在我们还不知道映射函数本身是什么的时候,有没有可能直接判断一个函数是不是核函数呢?

换句话说,一个函数需要具备怎样的性质,才是一个核函数

以下是函数 k(⋅,⋅) 可以作为核函数的充分必要条件

设 X 是输入空间,k(⋅,⋅) 是定义在 X×X 上的对称函数,当且仅当任意 x1,x2,...,xm∈X, 核矩阵 K(如下所示)总是半正定时,k(⋅,⋅) 就可以作为核函数使用。

核矩阵:

image

注意:下面是几个线性代数的基础概念。

对称函数:函数的输出值不随输入变量的顺序改变而改变的函数叫做对称函数。

因为输入空间中的 xi和xj 是向量,特征空间中的 ϕ(xi)和ϕ(xj) 也是向量,k(xi,xj) 表示特征空间向量的内积,而两个向量的内积并不因为其顺序变化而变化。因此,k(xi,xj)=k(xj,xi),即 k(⋅,⋅) 为对称函数。

相应地,核矩阵 K 为对称矩阵。

半正定矩阵:一个 n×n 的实对称矩阵 M 为半正定,当且仅当对于所有非零实系数向量 z,均有: zTMz⩾0。

其中“非零实系数向量”的含义是:一个向量中所有元素均为实数,且其中至少有一个元素值非零。

核函数的种类

要知道,非线性 SVM 的关键在于将输入空间中线性不可分的样本映射到线性可分的特征空间中去。特征空间的好坏直接影响到了 SVM 的效果。

****如何选择核函数也就成了一个关键性的问题。

虽然我们已经学习了核函数的定义和性质,但让我们凭空去自己构建一个核函数出来,还是一个非常困难的任务。

即使真得构造出来了,这个核函数在具体问题上的效果如何,还需要通过大量测试来证明,很可能费了好大劲,最后效果还不理想。

好在前人已经给我们留下来一些经历了长久磨练的常用核函数,让我们可以直接拿来用。接下来,我们分别看下这些核函数。

线性核(Linear Kernel)

k(xi,xj)=xTixj

使用时无须指定参数(Parameter),它直接计算两个输入向量的内积。经过线性核函数转换的样本,特征空间与输入空间重合,相当于并没有将样本映射到更高维度的空间里去。

很显然这是最简单的核函数,实际训练、使用 SVM 的时候,在不知道用什么核的情况下,可以先试试线性核的效果。

多项式核(Polynomial Kernel)

k(xi,xj)=(γxTixj+r)d,γ>0,d⩾1

需要指定三个参数:γ、r 和 d。

这是一个不平稳的核,适用于数据做了归一化(Normalization,参见本文最后一节)的情况。

RBF 核(Radial Basis Function Kernel)

k(xi,xj)=exp(−γ||xi−xj||2),γ>0

RBF 核又名高斯核 (Gaussian Kernel),是一个核函数家族。它会将输入空间的样本以非线性的方式映射到更高维度的空间(特征空间)里去,因此它可以处理类标签和样本属性之间是非线性关系的状况。

它有一个参数:γ,这个参数的设置非常关键!

如果设置过大,则整个 RBF 核会向线性核方向退化,向更高维度非线性投影的能力就会减弱;但如果设置过小,则会使得样本中噪音的影响加大,从而干扰最终 SVM 的有效性。

不过相对于多项式核的3个参数,RBF 核只有一个参数需要调,还是相对简单的。

当线性核效果不是很好时,可以用 RBF 试试。或者,很多情况下可以直接使用 RBF。

著名的 LIBSVM(台湾大学林智仁/Lin Chih-Jen 教授等设计开发的一款简单、易用、快速有效的 SVM/SVR 支持包)的默认核函数,就是 RBF 核。

Sigmoid 核(Sigmoid Kernel)

k(xi,xj)=tanh(γxTixj+r)

有两个参数,γ 和 r,在某些参数设置之下,Sigmoid 核矩阵可能不是半正定的,此时 Sigmoid 核也就不是有效的核函数了。因此参数设置要非常小心。

整体而言,Sigmoid 核并不比线性核或者 RBF 核更好。但是,当参数设置适宜时,它会有不俗的表现。

在具体应用核函数时,最好针对具体问题参照前人的经验。

构建自己的核函数

除了常见的核函数,我们还可以根据以下规律进行核函数的组合。

  1. 与正数相乘: k(⋅,⋅) 是核函数,对于任意正数(标量)α⩾0,αk(⋅,⋅) 也是核函数。

  2. 与正数相加: k(⋅,⋅) 是核函数,对于任意正数(标量)α⩾0,α+k(⋅,⋅) 也是核函数。

  3. 线性组合: k(⋅,⋅) 是核函数,则如下线性组合也是核函数:∑ni=1αiki(⋅,⋅),αi⩾0。

  4. 乘积: k1(⋅,⋅) 和 k2(⋅,⋅) 都是核函数,则它们的乘积——k1(⋅,⋅)k2(⋅,⋅)——也是核函数。

  5. 正系数多项式函数: 设 P 为实数域内的正系数多项式函数,k(⋅,⋅) 是核函数,则 P(k(⋅,⋅)) 也是核函数。

  6. 指数函数: k(⋅,⋅) 是核函数,则 exp(k(⋅,⋅)) 也是核函数。

掌握了这些规律,我们也可以尝试根据需要构建自己的核函数。

数据归一化(Data Normalization)

数据归一化是一种数据处理方法,具体所做的就是对取值范畴不同的数据进行归一化处理,使它们处在同一数量级上。

最常见的,就是把各种数据都变成 (0,1) 之间的小数。

enter image description here

上图是一个归一化过程在二维中的直观显示。

大家可以看到, “扁长”的数据分布,经过归一化处理之后,变成了一个“正圆”。

常用的归一化算法有:

(1)线性转换

x′=(x−min)(max−min)

(2)标准分

x′=(x–μ)γ

原本的 x 符合正态分布,μ 为其分布的均值,而 γ 为分布的方差。

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 158,736评论 4 362
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 67,167评论 1 291
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 108,442评论 0 243
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 43,902评论 0 204
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 52,302评论 3 287
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 40,573评论 1 216
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 31,847评论 2 312
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 30,562评论 0 197
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 34,260评论 1 241
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 30,531评论 2 245
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 32,021评论 1 258
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 28,367评论 2 253
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 33,016评论 3 235
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 26,068评论 0 8
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 26,827评论 0 194
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 35,610评论 2 274
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 35,514评论 2 269

推荐阅读更多精彩内容