支持向量机SVM算法学习(二)
2016-04-07 22:32
218 查看
【转载请注明出处】http://www.cnblogs.com/jerrylead
6 拉格朗日对偶(Lagrange duality)
先抛开上面的二次规划问题,先来看看存在等式约束的极值问题求法,比如下面的最优化问题:
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131234477453.png)
目标函数是f(w),下面是等式约束。通常解法是引入拉格朗日算子,这里使用
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131234489306.png)
来表示算子,得到拉格朗日公式为
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131234483700.png)
L是等式约束的个数。
然后分别对w和
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131234492729.png)
求偏导,使得偏导数等于0,然后解出w和
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131234502073.png)
。至于为什么引入拉格朗日算子可以求出极值,原因是f(w)的dw变化方向受其他不等式的约束,dw的变化方向与f(w)的梯度垂直时才能获得极值,而且在极值处,f(w)的梯度与其他等式梯度的线性组合平行,因此他们之间存在线性关系。(参考《最优化与KKT条件》)
然后我们探讨有不等式约束的极值问题求法,问题如下:
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131234511451.png)
我们定义一般化的拉格朗日公式
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131234524417.png)
这里的
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131234538254.png)
和
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131234538155.png)
都是拉格朗日算子。如果按这个公式求解,会出现问题,因为我们求解的是最小值,而这里的
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131234547217.png)
已经不是0了,我们可以将
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131234556071.png)
调整成很大的正值,来使最后的函数结果是负无穷。因此我们需要排除这种情况,我们定义下面的函数:
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131234552417.png)
这里的P代表primal。假设
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131234563464.png)
或者
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/20110313123456366.png)
,那么我们总是可以调整
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131234579395.png)
和
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131234586821.png)
来使得
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131234582643.png)
有最大值为正无穷。而只有g和h满足约束时,
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131234593690.png)
为f(w)。这个函数的精妙之处在于
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235003.png)
,而且求极大值。
因此我们可以写作
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235002969.png)
这样我们原来要求的min f(w)可以转换成求
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/20110313123501427.png)
了。
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235029805.png)
我们使用
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235035246.png)
来表示
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235036293.png)
。如果直接求解,首先面对的是两个参数,而
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235045637.png)
也是不等式约束,然后再在w上求最小值。这个过程不容易做,那么怎么办呢?
我们先考虑另外一个问题
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235057523.png)
D的意思是对偶,
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235053553.png)
将问题转化为先求拉格朗日关于w的最小值,将
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235061818.png)
和
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235072799.png)
看作是固定值。之后在
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235087193.png)
求最大值的话:
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/20110313123508158.png)
这个问题是原问题的对偶问题,相对于原问题只是更换了min和max的顺序,而一般更换顺序的结果是Max Min(X) <= MinMax(X)。然而在这里两者相等。用
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235105043.png)
来表示对偶问题如下:
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235108008.png)
下面解释在什么条件下两者会等价。假设f和g都是凸函数,h是仿射的(affine,
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235118814.png)
)。并且存在w使得对于所有的i,
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235113449.png)
。在这种假设下,一定存在
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235129272.png)
使得
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235124189.png)
是原问题的解,
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/20110313123514710.png)
是对偶问题的解。还有
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/20110313123514644.png)
另外,
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235158418.png)
满足库恩-塔克条件(Karush-Kuhn-Tucker,
KKT condition),该条件如下:
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235189201.png)
所以如果
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235188611.png)
满足了库恩-塔克条件,那么他们就是原问题和对偶问题的解。让我们再次审视公式(5),这个条件称作是KKT
dual complementarity条件。这个条件隐含了如果
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235197956.png)
,那么
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235199003.png)
。也就是说,
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/2011031312352016.png)
时,w处于可行域的边界上,这时才是起作用的约束。而其他位于可行域内部(
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235215839.png)
的)点都是不起作用的约束,其
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235215216.png)
。这个KKT双重补足条件会用来解释支持向量和SMO的收敛测试。
这部分内容思路比较凌乱,还需要先研究下《非线性规划》中的约束极值问题,再回头看看。KKT的总体思想是将极值会在可行域边界上取得,也就是不等式为0或等式约束里取得,而最优下降方向一般是这些等式的线性组合,其中每个元素要么是不等式为0的约束,要么是等式约束。对于在可行域边界内的点,对最优解不起作用,因此前面的系数为0。
7 最优间隔分类器(optimal margin classifier)
重新回到SVM的优化问题:
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235229818.png)
我们将约束条件改写为:
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235232260.png)
从KKT条件得知只有函数间隔是1(离超平面最近的点)的线性约束式前面的系数
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235239686.png)
,也就是说这些约束式
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235243557.png)
,对于其他的不在线上的点(
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235241016.png)
),极值不会在他们所在的范围内取得,因此前面的系数
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235254820.png)
.注意每一个约束式实际就是一个训练样本。
看下面的图:
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235267786.png)
实线是最大间隔超平面,假设×号的是正例,圆圈的是负例。在虚线上的点就是函数间隔是1的点,那么他们前面的系数
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235273227.png)
,其他点都是
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235282571.png)
。这三个点称作支持向量。构造拉格朗日函数如下:
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235281633.png)
注意到这里只有
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235299898.png)
没有
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235307847.png)
是因为原问题中没有等式约束,只有不等式约束。
下面我们按照对偶问题的求解步骤来一步步进行,
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235316385.png)
首先求解
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235313288.png)
的最小值,对于固定的
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235321585.png)
,
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235323712.png)
的最小值只与w和b有关。对w和b分别求偏导数。
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235339186.png)
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235346088.png)
并得到
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235346579.png)
将上式带回到拉格朗日函数中得到,此时得到的是该函数的最小值(目标函数是凸函数)
代入后,化简过程如下:
![](http://pic002.cnblogs.com/images/2011/279228/2011051016205860.png)
![](http://pic002.cnblogs.com/images/2011/279228/2011051016211679.png)
最后得到
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235354876.png)
由于最后一项是0,因此简化为
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235352302.png)
这里我们将向量内积
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235369204.png)
表示为
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235371647.png)
此时的拉格朗日函数只包含了变量
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235371232.png)
。然而我们求出了
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235383957.png)
才能得到w和b。
接着是极大化的过程
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235394447.png)
,
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235401316.png)
前面提到过对偶问题和原问题满足的几个条件,首先由于目标函数和线性约束都是凸函数,而且这里不存在等式约束h。存在w使得对于所有的i,
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235409855.png)
。因此,一定存在
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235414216.png)
使得
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235424117.png)
是原问题的解,
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235437606.png)
是对偶问题的解。在这里,求
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235445871.png)
就是求
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235451868.png)
了。
如果求出了
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235467865.png)
,根据
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235463131.jpg)
即可求出w(也是
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235479128.png)
,原问题的解)。然后
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235489062.png)
即可求出b。即离超平面最近的正的函数间隔要等于离超平面最近的负的函数间隔。
关于上面的对偶问题如何求解,将留给下一篇中的SMO算法来阐明。
这里考虑另外一个问题,由于前面求解中得到
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/20110313123548981.png)
我们通篇考虑问题的出发点是
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235482900.png)
,根据求解得到的
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235492485.png)
,我们代入前式得到
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235502419.png)
也就是说,以前新来的要分类的样本首先根据w和b做一次线性运算,然后看求的结果是大于0还是小于0,来判断正例还是负例。现在有了
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235513192.png)
,我们不需要求出w,只需将新来的样本和训练数据中的所有样本做内积和即可。那有人会说,与前面所有的样本都做运算是不是太耗时了?其实不然,我们从KKT条件中得到,只有支持向量的
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235529190.png)
,其他情况
![](http://images.cnblogs.com/cnblogs_com/jerrylead/201103/201103131235531042.png)
。因此,我们只需求新来的样本和支持向量的内积,然后运算即可。这种写法为下面要提到的核函数(kernel)做了很好的铺垫。这是上篇,先写这么多了。
相关文章推荐
- Asp.Net常用方法
- 分布式设计很简单--guzz分布式切表功能正式发布
- Bitmap
- 价值观作业
- POJ-3411 Paid Roads
- [Unity] unity3d中脚本生命周期(MonoBehaviour lifecycle)
- selenium操作隐藏的元素
- 理解逻辑斯蒂回归模型
- centos 6.5 用户管理
- 支持向量机SVM算法学习(一)
- IOS中的数据存储方式,特点,使用情况
- 问卷调查
- Linux系统的下载及VMware的下载
- 并查集 hdu3038 How Many Answers Are Wrong
- XML解析之pull解析
- IOS中的数据存储方式,特点,使用情况
- Android实现图片滚动控件,含页签功能,让你的应用像淘宝一样炫起来
- Linux内核【链表】整理笔记(1)
- angular UI-Router
- Delphi在创建和使用DLL的时候如果使用到string,请引入ShareMem单元