Andrew Ng机器学习笔记+Weka相关算法实现(五)SVM最优间隔和核方法
这一章主要解说Ng的机器学习中SVM的兴许内容。主要包括最优间隔分类器求解。核方法。
- 最优间隔分类器的求解
利用以一篇讲过的的原始对偶问题求解的思路,我们能够将相似思路运用到SVM的求解上来。
详细的分析例如以下:
对于SVM求解的问题:
我们把约束条件略微变形一下:
仅仅有函数间隔是1的点才干使上式取等号,也就是有意义的。例如以下图:
叉叉和圈圈分别代表正反例,能够看出,仅仅有落在边缘的点的α≠0,这些点才是支持向量。其它的点α=0,对切割超平面没有意义。上图的支持向量一共同拥有3个。
写出拉格朗日表达式:
依据对偶问题求解思路:
解对偶问题时,首先要固定α,以 w,b 为变量,最小化 L;
最小化 L 时,求解 L 对 w和 b 的偏导。并将导数设为 0。能够得到:
解出来:
代入原始的拉格朗日表达式:
部分推导步骤例如以下(感谢大神分享):
最后一项是0,原式化简为:
注意到:是向量的内积,我们能够进一步简化书写。
此时的拉格朗日函数仅仅包括了变量α,我们求出了α就能够得到b和w。
接下来极大化最小值:
这里的目标函数和线性约束都是凸函数,不存在等式约束 h。存在 w使得对于全部的 i。gi(W)< 0。因此,一定存在w∗, α∗使
得w∗是原问题的解, α∗是对偶问题的解。接下来就是求得α。
然后就能够得到:
也能够得到W:
即离超平面近期的正的函数间隔要等于离超平面近期的负的函数间隔。这个不难理解,就是仅仅有是平分线时候才离两側都最远。
求得W后,我们得到:
对于怎样求得α,在下一篇博客中介绍。
- 核方法
对于開始时提出来的房屋价格与面积的关系问题,设XY各自是房屋面积和价格,假设我们知道x和y满足三次函数关系,那么我们就能够使用三次多项式拟合x和y的关系。可是我们首先要将特征扩展到(x,x^2,x^3),这个过程叫做特征映射。我们表示为φ,在这个问题中:
假设我们希望将映射后的特征用于SVM的分类,就要将内积
映射为:
假设例子可能存在线性不可分的情况,而将特征映射到高维空间后。往往就可分了。
我们定义核函数:
我们仅仅须要做计算,就能够完毕映射,可是这样做复杂度非常高。比方最初的特征是 n 维的,我们将其映射到n^2维,然后再计算,这样须要O(n^2)的时间。我们考虑简化方法:
对于一个特例:假设 x 和 z 都是 n 维的,
展开后:
不难看出。这个时候发现我们能够仅仅计算原始特征 x 和 z 内积的平方,就等价与计算映射后特征的内积。此时:
对于还有一个相似的核函数:
我们能够表示为:
更一般的表示。对于一个形如:
的核函数。其映射后的特征是 维度。 向量中每一个元素都是最高为 d 阶的变量的组合。
假设直接计算核函数的值。其复杂度均为 O(n),这就避免了直接计算内积结果的复杂度。
假设ϕ(x)和ϕ(z)在其相应的维度空间中位置接近,那么内积值 K 会非常大,
反之则内积会小。这有点相似两个向量的乘积。这意味着核函数 K 是一个度量向量 x 和向量 z接近程度的指标 。
假设 x 和 z 相差非常大( ||x − z|| ≫ 0)。那么核函数值约等于 0。
由于这个函数相似于高斯分布,因此称为高斯核
函数,也叫做径向基函数(Radial Basis Function 简称 RBF)。
它能够把原始特征映射到无穷维(泰勒展开式有无穷项多项式)。
使用了核函数后,原始的 变成了
,我们不必求解核函数表达式,仅仅须要将内积替换为
就可以。
那么什么样的函数才干作为核函数呢?
给定一个函数K,是否存在φ,使得 对全部的x和z都成立?也就是说,假设核函数合法,那么必定能写成两个映射函数乘积的形式。
首先定义核矩阵。对于一个数据集, 定义一个
m*m 的矩阵 K( K 既代表核函数也代表核矩阵)。 K 中的每一个元素定义例如以下:
这样能够计算出 m*m 的核函数矩阵( Kernel Matrix)。
依据核函数定义,有下式子:
也就是说,矩阵K是对称的。我们令 代表矩阵K的第k维,那么对于随意的向量Z,我们得到:
从这个公式我们能够看出,假设 K 是个有效的核函数。那么。在训练集上得到的核函数矩阵 K 应该是半正定的( K ≥ 0)。结合曾经的结论,K矩阵应该是对称半正定的。
这个条件也是充分的,由 Mercer 定理来表达:
核函数不仅仅用在 SVM 上,仅仅要是计算中有内积出现。就能够替换为核函数。
- 不可分情况的处理
对于低维不可分问题,我们能够通过核函数映射到高维,一般就能够线性可分,可是仍然不能保证一定能够线性可分,那么我们应该採用零哟中思路找到一个最合适的超平面进行切割。例如以下图所看到的:
第二图由于一个噪声点可能会造成切割平面移动甚至线性不可分。我们应该同意一些点游离在我们定义的规则之外,此时的条件例如以下:
原有的问题增加了非负參数ξ。可是放松限制条件后。我们需
要又一次调整目标函数,以对离群点进行处罚,也就是增加了惩处项C,C越大。代表越不希望出现离群点。
能够得到修改后的拉格朗日表达式:
依照之前的求解思路:
能够看到公式并没有什么变化,仅仅是多了一个限制条件。可是事实上b的求解也会发生变化,详细到SMO算法的时候再解说。
KKT条件也会发生变化:
Andrew Ng机器学习笔记+Weka相关算法实现(五)SVM最优间隔和核方法的更多相关文章
- Andrew Ng机器学习笔记+Weka相关算法实现(四)SVM和原始对偶问题
这篇博客主要解说了Ng的课第六.七个视频,涉及到的内容包含,函数间隔和几何间隔.最优间隔分类器 ( Optimal Margin Classifier).原始/对偶问题 ( Primal/Dual P ...
- 斯坦福大学Andrew Ng - 机器学习笔记(5) -- 支持向量机(SVM)
大概用了一个月,Andrew Ng老师的机器学习视频断断续续看完了,以下是个人学习笔记,入门级别,权当总结.笔记难免有遗漏和误解,欢迎讨论. 鸣谢:中国海洋大学黄海广博士提供课程视频和个人笔记,在此深 ...
- 斯坦福大学Andrew Ng - 机器学习笔记(4) -- 机器学习算法的选择与评估
大概用了一个月,Andrew Ng老师的机器学习视频断断续续看完了,以下是个人学习笔记,入门级别,权当总结.笔记难免有遗漏和误解,欢迎讨论. 鸣谢:中国海洋大学黄海广博士提供课程视频和个人笔记,在此深 ...
- 吴恩达(Andrew Ng)——机器学习笔记1
之前经学长推荐,开始在B站上看Andrew Ng的机器学习课程.其实已经看了1/3了吧,今天把学习笔记补上吧. 吴恩达老师的Machine learning课程共有113节(B站上的版本https:/ ...
- 斯坦福大学Andrew Ng - 机器学习笔记(8) -- 推荐系统 & 大规模机器学习 & 图片文字识别
大概用了一个月,Andrew Ng老师的机器学习视频断断续续看完了,以下是个人学习笔记,入门级别,权当总结.笔记难免有遗漏和误解,欢迎讨论. 鸣谢:中国海洋大学黄海广博士提供课程视频和个人笔记,在此深 ...
- 斯坦福大学Andrew Ng - 机器学习笔记(7) -- 异常检测
大概用了一个月,Andrew Ng老师的机器学习视频断断续续看完了,以下是个人学习笔记,入门级别,权当总结.笔记难免有遗漏和误解,欢迎讨论. 鸣谢:中国海洋大学黄海广博士提供课程视频和个人笔记,在此深 ...
- 斯坦福大学Andrew Ng - 机器学习笔记(6) -- 聚类 & 降维
大概用了一个月,Andrew Ng老师的机器学习视频断断续续看完了,以下是个人学习笔记,入门级别,权当总结.笔记难免有遗漏和误解,欢迎讨论. 鸣谢:中国海洋大学黄海广博士提供课程视频和个人笔记,在此深 ...
- 斯坦福大学Andrew Ng - 机器学习笔记(3) -- 神经网络模型
大概用了一个月,Andrew Ng老师的机器学习视频断断续续看完了,以下是个人学习笔记,入门级别,权当总结.笔记难免有遗漏和误解,欢迎讨论. 鸣谢:中国海洋大学黄海广博士提供课程视频和个人笔记,在此深 ...
- 斯坦福大学Andrew Ng - 机器学习笔记(2) -- 逻辑回归 & 正则化
大概用了一个月,Andrew Ng老师的机器学习视频断断续续看完了,以下是个人学习笔记,入门级别,权当总结.笔记难免有遗漏和误解,欢迎讨论. 鸣谢:中国海洋大学黄海广博士提供课程视频和个人笔记,在此深 ...
随机推荐
- servlet多线程同步问题
Servlet/JSP技术和ASP.PHP等相比,由于其多线程运行而具有很高的执行效率.•由于Servlet/JSP默认是以多线程模式执行的,所以,在编写代码时需要非常细致地考虑多线程的同步问题.•如 ...
- poj2728(最小比率生成树)
poj2728 题意 给出 n 个点的坐标和它的高度,求一颗生成树使得树上所连边的两点高度差之和除以距离之和最小. 分析 01分数规划+最小生成树. 对于所有的边,在求最小生成树过程中有选或不选的问题 ...
- Java多线程中的异常处理
在java多线程程序中,所有线程都不允许抛出未捕获的checked exception,也就是说各个线程需要自己把自己的checked exception处理掉.这一点是通过java.lang.Run ...
- Count and Say (Array Length Encoding) -- LeetCode
The count-and-say sequence is the sequence of integers beginning as follows:1, 11, 21, 1211, 111221, ...
- js日常笔记
写在前面: 在工作中,有时候会遇到一些零零碎碎的小知识点,虽然这些网上都可以查询到,但还是想把一些自己不是很熟悉的当做笔记记录下来,方便以后查询. 1.按钮隐藏/显示/可用/不可用 $("# ...
- C++指针、引用知多少?
上午搞了一个小程序,test半天都没有得到想要的结构,原来是递归的时候没有注意的循环的问题,结果直接死循环了.催了...看来当程序出现问题的时候,首先要整理的是算法思路是否有问题,其次是算法的实现,是 ...
- 1.常量&变量
简单体验 // 定义变量 var i = 10 println(i) i = 15 println(i) let j = 20 // 常量一经定义不能自改数值 // j = 25 println(j) ...
- linux的file指令
显示文件的类型,用命令 file 可以使你知道某个文件究竟是ELF格式的可执行文件, 还是shell script文 件或是其他的什么格式 例如:#file startx 语 法:file [-beL ...
- ASP.NET 5基础之中间件
来源https://docs.asp.net/en/latest/fundamentals/middleware.html 一些可以整合进http请求管道的小的应用组件称做中间件.ASP.NET 5集 ...
- JAVA常见算法题(十八)
package com.xiaowu.demo; /** * 两个乒乓球队进行比赛,各出三人.甲队为a,b,c三人,乙队为x,y,z三人,以抽签决定比赛名单. 有人向队员打听比赛的名单:a说他不和x比 ...