精华内容
下载资源
问答
  • 4.7 向量优化 广义和凸的向量优化问题 最优解与值 Pareto最优解与值 标量化 多准则优化 例子 广义和凸的向量优化问题 广义向量优化问题 其中是优化变量,是正常锥,是目标函数,是不等式约束函数,是等式...

    4.7 向量优化

    1. 广义和凸的向量优化问题
    2. 最优解与值
    3. Pareto最优解与值
    4. 标量化
    5. 多准则优化
    6. 例子

    广义和凸的向量优化问题

    广义向量优化问题

    minimize(K)\,\, f_0(x)\\ subject \, \, to \, \,\begin{matrix} f_i(x)\leq 0,i=1,2,\cdots m\\ h_i(x)=0,i=1,2\cdots p \end{matrix}

    其中x\in R^n是优化变量,K\subseteq R^q是正常锥,f_0:R^n\rightarrow R^q是目标函数,f_i:R^n\rightarrow R是不等式约束函数,h_i:R^n\rightarrow R是等式约束函数。其余标准优化问题的区别在于目标函数不是标量而是向量,且比较目标函数值的方法不再是简单的数的大小,而是有正常锥K定义的广义不等式。故也称标准优化问题为标量优化问题。

    凸向量优化问题

    minimize(K)\,\, f_0(x)\\ subject \, \, to \, \,\begin{matrix} f_i(x)\leq 0,i=1,2,\cdots m\\ Ax=b \end{matrix}

    A\in R^{p\times n}

    即如果目标函数是K-凸的,其不等式约束函数f_i,i=1,2\cdots p是凸的,等式约束函数是仿射的,则称问题是凸向量优化问题。

     

    对于向量优化问题,我们是要寻找一个目标函数值在K锥上最优,即最小。但是因为向量优化问题的目标函数是向量,向量并不一定是可比较的,所以就要分情况考虑,类似第二章的最小元和极小元,在向量优化问题中,有两种情况,即最优解和Pareto最优解。

    最优解与值

    考虑可行点的目标值的集合:

    O=\left \{ f_0(x)|\exists x \in D,f_i(x)\leq 0,i=1,2\cdots m,h_i(x)=0,i=1,2,\cdots p \right \}\sqsubseteq R^q

    称为可达目标值集合,显然如果这个集合中有最小元,那么问题就有最优解,最优解对应的目标函数值为最优值。

    回顾最小元的定义:

    一个集合C有最小元x当且仅当\forall y \in C,x\preceq _K y,当且仅当C\subseteq x+K(x+K表示可以与x相比并且大于等于x的所有元素。)即C中所有元素都可以与x比较,且小于等于x。

    考虑向量优化问题有最优解的情况,即可达目标集有最小元的情况,因此向量问题含有最优解x^*,即O\subseteq f_0(x^*)+K,如下图:

    Pareto最优解与值

    现考虑可达目标集不含最小元的情况,此时问题不含最优解和最优值,只能在科大目标集中寻找极小元,我们称可达目标集中的极小元为Pareto最优,其对应的目标函数值为Pareto最优值。

    回顾集合极小元的定义:

    一个集合C有最小元x当且仅当\exists y \in C,x\preceq _K y\Rightarrow y=x,当且仅当C\cap (x-K)=\left \{ x \right \}(x-K表示可以与x相比并且小于等于x的所有元素。)即C中某些元素可以与x比较,且所有可以与x比较的元素都小于等于x。

    所以点x是Pareto最优解,当且仅当它是可行的,而且(f_0(x)-K)\cap O=\left \{ f_0(x) \right \},如下图:

    一个向量优化问题可以有很多Pareto最优解,Pareto最优值的集合记为P,它满足P\subseteq O\cap bd(O),即每一个Pareto最优值都位于可达目标集合变价上的可达标目标值。

    标量化

    基于对偶广义不等式的极小和最小点的特征,选择任意的\lambda \succeq _{K^*} 0,即在任意对偶广义不等式中为正的向量,考虑标量优化问题

    minimize \, \,\lambda ^Tf_0(x) \\ subject \, \, to \, \,\begin{matrix} f_i(x)\leq 0 & i=1,2\cdots m\\ h_i(x)=0 & i=1,2\cdots p \end{matrix}

    标量优化问题的最优解是对应的向量优化问题的Pareto最优解。向量\lambda有时称为权向量。

    上图显示了不同的\lambda值对应不同的Pareto最优解,注:某些Pareto最优点不能有任何权向量\lambda \succeq _{K^*} 0找到。

    从几何角度来看,点x是标量化问题的最优解,即在可行集上极小化了\lambda ^Tf_0,当且仅当对于所有可行的y,有\lambda ^T(f_0(y)-f_0(x))\geq 0相当于\left \{ u|-\lambda ^T(u-f_0(x))=0 \right \}是可达目标值集合O在f_0(x)处的支撑超平面。

    对于凸向量优化问题,当权向量遍历K^*-非负的非零向量,标量化方法可以得到所有Pareto最优解。

    多准则优化

    当向量优化函数关于锥K=R_+^q时,他称为多准则或多目标优化问题。

    f_0(x)=(F_1(x),F_2(x),\cdots F_q(x)),q个不同的目标函数,简单地说就是同时极小化这q个目标函数的值。

    如果f_1\cdots f_q,h_1\cdots h_p,F_1\cdots F_q是凸的,则该多准则优化问题是凸问题。

    因为是多个目标函数同时极小化,一种好的情况是可以得到一个最优解x同时极小化这q个目标函数,但很多情况下往往找不到一个同时极小化这q个目标函数的最优解,很可能只能极小化其中的一部分,·而不能极大化其他的目标函数,这时候就需要作出一些取舍和权衡,选择极小化更重要的目标函数。

    标量化多准则问题

    跟标量化类似,通过加权和目标函数来标量化多准则问题

    \lambda ^Tf_0(x)=\sum _{i=1}^q\lambda _iF_i(x)

    所以通过设置权向量就可以有倾向地选择目标函数,比如可以为更重要的目标函数分量设置更大的权值。

    例子

    正则化最小二乘

    minimize \, \, f_0(x)=(\begin{Vmatrix} Ax-b \end{Vmatrix}_2^2,\begin{Vmatrix} x \end{Vmatrix}_2^2)

    \lambda ^Tf_0(x)=\lambda _1F_1(x)+\lambda _2F_2(x)=\lambda _1\begin{Vmatrix}Ax-b \end{Vmatrix}_2^2+\lambda _2\begin{Vmatrix} x\end{Vmatrix}_2^2

    =\lambda _1(x^TA^TAx-2b^TAx+b^Tb)+\lambda _2x^Tx

    =x^T(\lambda _1A^TA+\lambda _2I)x-2\lambda _1b^TAx+\lambda _1b^Tb

    对x求导:

    \frac{\partial \lambda ^Tf_0(x)}{\partial x}=2(\lambda _1A^TA+\lambda _2I)x-2\lambda _1A^Tb

    令其为0,得到x=(\lambda _1A^TA+\lambda _2I)^{-1}\lambda _1A^Tb,整理得到x=(\lambda _1A^TA+\lambda _2I)^{-1}\lambda _1A^Tb =(A^TA+\frac{\lambda _2}{\lambda _1}I)^{-1}\frac{\lambda _1}{\lambda _1}A^Tb=(A^TA+\mu I)^{-1}A^Tb

    其中\mu =\lambda _2 /\lambda _1\mu趋于无穷大时,即\lambda _1=0时,最优解就是下图左侧的点,\mu趋于0时,即\lambda _2=0时,最优解就是下图右侧的点。

     

    来源:https://blog.csdn.net/wangchy29/article/details/86666185

    展开全文
  • 4.7向量优化 广义和凸的向量优化问题 最优解与值 Pareto最优解与值 标量化 多准则优化 例子 广义和凸的向量优化问题 广义向量优化问题 其中是优化变量,是正常锥,是目标函数,是不等式约束函数,是等式...

    4.7向量优化

    1. 广义和凸的向量优化问题
    2. 最优解与值
    3. Pareto最优解与值
    4. 标量化
    5. 多准则优化
    6. 例子

    广义和凸的向量优化问题

    广义向量优化问题

    minimize(K)\,\, f_0(x)\\ subject \, \, to \, \,\begin{matrix} f_i(x)\leq 0,i=1,2,\cdots m\\ h_i(x)=0,i=1,2\cdots p \end{matrix}

    其中x\in R^n是优化变量,K\subseteq R^q是正常锥,f_0:R^n\rightarrow R^q是目标函数,f_i:R^n\rightarrow R是不等式约束函数,h_i:R^n\rightarrow R是等式约束函数。其余标准优化问题的区别在于目标函数不是标量而是向量,且比较目标函数值的方法不再是简单的数的大小,而是有正常锥K定义的广义不等式。故也称标准优化问题为标量优化问题。

    凸向量优化问题

    minimize(K)\,\, f_0(x)\\ subject \, \, to \, \,\begin{matrix} f_i(x)\leq 0,i=1,2,\cdots m\\ Ax=b \end{matrix}

    A\in R^{p\times n}

    即如果目标函数是K-凸的,其不等式约束函数f_i,i=1,2\cdots p是凸的,等式约束函数是仿射的,则称问题是凸向量优化问题。

     

    对于向量优化问题,我们是要寻找一个目标函数值在K锥上最优,即最小。但是因为向量优化问题的目标函数是向量,向量并不一定是可比较的,所以就要分情况考虑,类似第二章的最小元和极小元,在向量优化问题中,有两种情况,即最优解和Pareto最优解。

    最优解与值

    考虑可行点的目标值的集合:

    O=\left \{ f_0(x)|\exists x \in D,f_i(x)\leq 0,i=1,2\cdots m,h_i(x)=0,i=1,2,\cdots p \right \}\sqsubseteq R^q

    称为可达目标值集合,显然如果这个集合中有最小元,那么问题就有最优解,最优解对应的目标函数值为最优值。

    回顾最小元的定义:

    一个集合C有最小元x当且仅当\forall y \in C,x\preceq _K y,当且仅当C\subseteq x+K(x+K表示可以与x相比并且大于等于x的所有元素。)即C中所有元素都可以与x比较,且小于等于x。

    考虑向量优化问题有最优解的情况,即可达目标集有最小元的情况,因此向量问题含有最优解x^*,即O\subseteq f_0(x^*)+K,如下图:

    Pareto最优解与值

    现考虑可达目标集不含最小元的情况,此时问题不含最优解和最优值,只能在科大目标集中寻找极小元,我们称可达目标集中的极小元为Pareto最优,其对应的目标函数值为Pareto最优值。

    回顾集合极小元的定义:

    一个集合C有最小元x当且仅当\exists y \in C,x\preceq _K y\Rightarrow y=x,当且仅当C\cap (x-K)=\left \{ x \right \}(x-K表示可以与x相比并且小于等于x的所有元素。)即C中某些元素可以与x比较,且所有可以与x比较的元素都小于等于x。

    所以点x是Pareto最优解,当且仅当它是可行的,而且(f_0(x)-K)\cap O=\left \{ f_0(x) \right \},如下图:

    一个向量优化问题可以有很多Pareto最优解,Pareto最优值的集合记为P,它满足P\subseteq O\cap bd(O),即每一个Pareto最优值都位于可达目标集合变价上的可达标目标值。

    标量化

    基于对偶广义不等式的极小和最小点的特征,选择任意的\lambda \succeq _{K^*} 0,即在任意对偶广义不等式中为正的向量,考虑标量优化问题

    minimize \, \,\lambda ^Tf_0(x) \\ subject \, \, to \, \,\begin{matrix} f_i(x)\leq 0 & i=1,2\cdots m\\ h_i(x)=0 & i=1,2\cdots p \end{matrix}

    标量优化问题的最优解是对应的向量优化问题的Pareto最优解。向量\lambda有时称为权向量。

    上图显示了不同的\lambda值对应不同的Pareto最优解,注:某些Pareto最优点不能有任何权向量\lambda \succeq _{K^*} 0找到。

    从几何角度来看,点x是标量化问题的最优解,即在可行集上极小化了\lambda ^Tf_0,当且仅当对于所有可行的y,有\lambda ^T(f_0(y)-f_0(x))\geq 0相当于\left \{ u|-\lambda ^T(u-f_0(x))=0 \right \}是可达目标值集合O在f_0(x)处的支撑超平面。

    对于凸向量优化问题,当权向量遍历K^*-非负的非零向量,标量化方法可以得到所有Pareto最优解。

    多准则优化

    当向量优化函数关于锥K=R_+^q时,他称为多准则或多目标优化问题。

    f_0(x)=(F_1(x),F_2(x),\cdots F_q(x)),q个不同的目标函数,简单地说就是同时极小化这q个目标函数的值。

    如果f_1\cdots f_q,h_1\cdots h_p,F_1\cdots F_q是凸的,则该多准则优化问题是凸问题。

    因为是多个目标函数同时极小化,一种好的情况是可以得到一个最优解x同时极小化这q个目标函数,但很多情况下往往找不到一个同时极小化这q个目标函数的最优解,很可能只能极小化其中的一部分,·而不能极大化其他的目标函数,这时候就需要作出一些取舍和权衡,选择极小化更重要的目标函数。

    标量化多准则问题

    跟标量化类似,通过加权和目标函数来标量化多准则问题

    \lambda ^Tf_0(x)=\sum _{i=1}^q\lambda _iF_i(x)

    所以通过设置权向量就可以有倾向地选择目标函数,比如可以为更重要的目标函数分量设置更大的权值。

    例子

    正则化最小二乘

    minimize \, \, f_0(x)=(\begin{Vmatrix} Ax-b \end{Vmatrix}_2^2,\begin{Vmatrix} x \end{Vmatrix}_2^2)

    \lambda ^Tf_0(x)=\lambda _1F_1(x)+\lambda _2F_2(x)=\lambda _1\begin{Vmatrix}Ax-b \end{Vmatrix}_2^2+\lambda _2\begin{Vmatrix} x\end{Vmatrix}_2^2

    =\lambda _1(x^TA^TAx-2b^TAx+b^Tb)+\lambda _2x^Tx

    =x^T(\lambda _1A^TA+\lambda _2I)x-2\lambda _1b^TAx+\lambda _1b^Tb

    对x求导:

    \frac{\partial \lambda ^Tf_0(x)}{\partial x}=2(\lambda _1A^TA+\lambda _2I)x-2\lambda _1A^Tb

    令其为0,得到x=(\lambda _1A^TA+\lambda _2I)^{-1}\lambda _1A^Tb,整理得到x=(\lambda _1A^TA+\lambda _2I)^{-1}\lambda _1A^Tb =(A^TA+\frac{\lambda _2}{\lambda _1}I)^{-1}\frac{\lambda _1}{\lambda _1}A^Tb=(A^TA+\mu I)^{-1}A^Tb

    其中\mu =\lambda _2 /\lambda _1\mu趋于无穷大时,即\lambda _1=0时,最优解就是下图左侧的点,\mu趋于0时,即\lambda _2=0时,最优解就是下图右侧的点。

    ​​​​​​​

     

     

     

    展开全文
  • 非负矩阵分解NMF

    万次阅读 多人点赞 2016-08-03 12:37:26
    http://blog.csdn.net/pipisorry/article/details/52098864非负矩阵分解(NMF,Non-negative matrix factorization)NMF的发展及原理 著名的科学杂志《Nature》于...该文提出了一种新的矩阵分解思想——非负矩阵分解(Non

    http://blog.csdn.net/pipisorry/article/details/52098864

    非负矩阵分解(NMF,Non-negative matrix factorization)

    NMF的发展及原理

      著名的科学杂志《Nature》于1999年刊登了两位科学家D.D.Lee和H.S.Seung对数学中非负矩阵研究的突出成果。该文提出了一种新的矩阵分解思想——非负矩阵分解(Non-negative Matrix Factorization,NMF)算法,即NMF是在矩阵中所有元素均为非负数约束条件之下的矩阵分解方法。该论文的发表迅速引起了各个领域中的科学研究人员的重视:一方面,科学研究中的很多大规模数据的分析方法需要通过矩阵形式进行有效处理,而NMF思想则为人类处理大规模数据提供了一种新的途径;另一方面,NMF分解算法相较于传统的一些算法而言,具有实现上的简便性、分解形式和分解结果上的可解释性,以及占用存储空间少等诸多优点。为高效处理这些通过矩阵存放的数据,一个关键的必要步骤便是对矩阵进行分解操作。通过矩阵分解,一方面将描述问题的矩阵的维数进行削减,另一方面也可以对大量的数据进行压缩和概括。
      利用矩阵分解来解决实际问题的分析方法很多,如PCA(主成分分析)、ICA(独立成分分析)、SVD(奇异值分解)、VQ(矢量量化)等。在所有这些方法中,原始的大矩阵V被近似分解为低秩的V=WH形式。这些方法的共同特点是,因子W和H中的元素可为正或负,即使输入的初始矩阵元素是全正的,传统的秩削减算法也不能保证原始数据的非负性。在数学上,从计算的观点看,分解结果中存在负值是正确的,但负值元素在实际问题中往往是没有意义的。例如图像数据中不可能有负值的像素点;在文档统计中,负值也是无法解释的。

    NMF的基本思想

      NMF的基本思想可以简单描述为:对于任意给定的一个非负矩阵A,NMF算法能够寻找到一个非负矩阵U和一个非负矩阵V,使得满足 ,从而将一个非负的矩阵分解为左右两个非负矩阵的乘积。

    NMF问题描述

       传统的NMF问题可以描述如下

        给定矩阵,寻找非负矩阵和非负矩阵,使得


        分解前后可理解为:原始矩阵的列向量是对左矩阵中所有列向量的加权和,而权重系数就是右矩阵对应列向量的元素,故称为基矩阵,为系数矩阵。一般情况下的选择要比小,即满足,这时用系数矩阵代替原始矩阵,就可以实现对原始矩阵进行降维,得到数据特征的降维矩阵,从而减少存储空间,减少计算机资源。

        由于分解前后的矩阵中仅包含非负的元素,因此,原矩阵A中的一列向量可以解释为对左矩阵U中所有列向量(称为基向量)的加权和,而权重系数为右矩阵V中对应列向量中的元素。这种基于基向量组合的表示形式具有很直观的语义解释,它反映了人类思维中“局部构成整体”的概念。研究指出,非负矩阵分解是个NP问题,可以划为优化问题用迭代方法交替求解U和V。NMF算法提供了基于简单迭代的求解U,V的方法,求解方法具有收敛速度快、左右非负矩阵存储空间小的特点,它能将高维的数据矩阵降维处理,适合处理大规模数据。利用NMF进行文本、图像大规模数据的分析方法,较传统的处理算法速度更快、更便捷。

    由于NMF不允许基图像或中间的权重矩阵中出现负值,因此只有相加组合得到的正确基图像才允许,最后通过处理后的重构图像效果是比较满意的(对矩阵非负的限制使得这种分解能够达到用部分表达整体的效果,简单地说就是,整体由部分的叠加而没有了正负抵消 )。[Learning the parts of objects by non-negative matrix factorization]

    非负矩阵分解NMF的一个示例解释

      在VQ分解中,每一列的被约束为一个一元矢量。其中只有一个元素为1,其他元素为0。若的第一列中,第r1个元素为1,那么中第一列的脸,就完全由基图像中的第r1列数据表示。此时得到的基图像称为原型基图像,这些原型图像表示一张原型脸。

           在PCA分解中,的各列之间相互正交,各行之间相互正交。这个约束比VQ的松弛很多,也就是,中的元素可为正也可为负。中每一张脸的每一个像素点都是中各列对应的像素点的一个加权和。由于权重矩阵中元素符号的任意性,所以基矩阵表示出来并不像VQ中原型脸那样的直观可解释。此时将W的列数据画出来并不一定能直接看到一张“脸”。但是在统计上可以解释为最大方差方向,我们把这些“脸”称为“特征脸”。

           在NMF中,由于加了非负约束。与VQ的单一元素不为0不同,NMF允许基图像H间的加权结合来表示脸部图像V;与PCA不同,NMF的加权系数H中的元素都为非负的。前两者得到的都是一个完整的脸部特征基图像,而NMF得到的是脸部子特征。通俗点说,VQ是用一张完整的图像直接代表源脸部图像;PCA是将几个完整人脸加减压成一张脸;而NMF是取甲的眼睛,乙的鼻子,丙的嘴巴直接拼成一张脸。这样解释虽然细节上略有不妥,但不失其概念上的意义。

    通过图1中的面部特征提取例子可领略NMF处理数据的方式。最左边的大矩阵由一系列的小图组成,这些小图是分析数据库中包含的2429个脸部图像的结果,每幅图像由19×19个像素组成。传统方法中这样的小图是一幅完整的人脸图像,但是在NMF方法中,每个小图是通过一组基图像乘以一个权重矩阵而产生的面部特征图,经过这样处理的每幅小图像恰好表示了诸如“鼻子”、“嘴巴”、“眼睛”等人脸局部概念特征,这便大大压缩了存放的图像数据量。左边的大矩阵由每幅小图像的19列一起组成矩阵的一列,那样它就是19×19=361行,2429列。这个例子中,NMF方法用基图像来代表眼、眉毛、鼻子、嘴、耳朵、胡子等,它们一起组成了数据库中的脸。这样给人最先的直觉就是它很好地压缩了数据。事实上Lee和Seung在他们的论文中更深入地指出,与人类识别事物的过程相似,NMF也是一种优化的机制,近似于我们的脑分析和存储人脸数据的过程。这个例子中,原图像表示这些局部特征的加权组合,这与人类思维中“局部构成整体”的概念是相吻合的。因此,NMF算法似乎体现了一种智能行为。


                     图1  NMF提取面部特征的实例

    上图的第一个方块为矩阵W,组成的图像。其中每一个小格为W的一列的19*19个元素重构而成的19*19的矩阵图像。第二个方块为H矩阵,其中红色表示负数,灰黑表示正数, 颜色程度表示大小。右边的图像只是V矩阵的一列的19*19个元素组成的一张原始脸。

    NMF限制V,W,H都是非负的,而VQ限制H必须必须由单位向量(只有一个元素为1,其余为0)组成,PCA限制W的列是规范正交的,H的行是正交的,这导致了有图中三种不同的分解结果。
    从结果可以看出:NMF分解结果与我们部分形成整体的直觉是相一致的。NMF学习出一种基于部分的表达。

    皮皮blog


    非负矩阵分解NMF的应用

      事实上,在Lee和Seung发表他们的研究成果之前,针对非负矩阵的研究早在20世纪70年代已经有数学家做了一些相关的工作,但是没有引起过多的关注。20世纪90年代早期,科学家开始将数学上非负矩阵的研究成果用于环境处理和卫星遥控的应用,但是对于非负矩阵的应用意义和价值的理解仍只局限于少数科学家中,人们还没有广泛重视这种方法。直到1999年Lee和Seung的非负矩阵研究成果发表在《Nature》杂志之后,这一切才得以改变。尽管同年有另两位科学家也发表了与Lee和Seung相近的研究结果,但由于论文刊登在并非如《Nature》那样具有极高声誉的学术杂志上,因此其工作并没有得到如Lee和Seung同样的关注,这也从一个侧面折射了高水平学术杂志对研究工作的推动作用。
      二、应用领域
      NMF是一个很有效的算法,它力图在大规模的矩阵数据中发现具有解释功能的关系,相比当前文献中公布的其他方法来说,使用NMF的算法也是非常精确和快速的。NMF算法思想能为世界上权威的学术刊物所接受并非偶然,因为该理论本身蕴涵了巨大的潜能,这种潜在的力量将通过各种具体的应用来得以体现。

    在众多应用中,NMF能被用于发现数据库中的图像特征,便于快速自动识别应用;能够发现文档的语义相关度,用于信息自动索引和提取;能够在DNA阵列分析中识别基因等等。NMF能用于发现数据库中图像的特征,便于快速识别应用,比如实现录入恐怖分子的照片,然后在安检口对可疑人员进行盘查。在文档方面,NMF能够发现文档的语义相关度,用于信息的自动索引和提取。在生物学中,在DNA阵列分析中识别基因等。在语音识别系统中NMF也能发挥重要作用。有人已经将NMF应用于盲信号分离BSS领域,取得一定的效果,不过其研究的BSS有非欠定的要求(即源数不大于传感器数),且混合矩阵等都是正的。

      (1) 图像分析
      NMF最成功的一类应用是在图像的分析和处理领域。图像本身包含大量的数据,计算机一般将图像的信息按照矩阵的形式进行存放,针对图像的识别、分析和处理也是在矩阵的基础上进行的。这些特点使得NMF方法能很好地与图像分析处理相结合。人们已经利用NMF算法,对卫星发回的图像进行处理,以自动辨别太空中的垃圾碎片;使用NMF算法对天文望远镜拍摄到的图像进行分析,有助于天文学家识别星体;美国还尝试在机场安装由NMF算法驱动的识别系统,根据事先输入计算机的恐怖分子的特征图像库来自动识别进出机场的可疑恐怖分子。
      (2) 文本聚类/数据挖掘
      文本在人类日常接触的信息中占有很大分量,为了更快更精确地从大量的文本数据中取得所需要的信息,针对文本信息处理的研究一直没有停止过。文本数据不光信息量大,而且一般是无结构的。此外,典型的文本数据通常以矩阵的形式被计算机处理,此时的数据矩阵具有高维稀疏的特征,因此,对大规模文本信息进行处理分析的另一个障碍便是如何削减原始数据的维数。NMF算法正是解决这方面难题的一种新手段。NMF在挖掘用户所需数据和进行文本聚类研究中都有着成功的应用例子。由于NMF算法在处理文本数据方面的高效性,著名的商业数据库软件Oracle在其第10版中专门利用NMF算法来进行文本特征的提取和分类。为什么NMF对于文本信息提取得很好呢?原因在于智能文本处理的核心问题是以一种能捕获语义或相关信息的方式来表示文本,但是传统的常用分析方法仅仅是对词进行统计,而不考虑其他的信息。而NMF不同,它往往能达到表示信息的局部之间相关关系的效果,从而获得更好的处理结果。
      (3) 语音处理
      语音的自动识别一直是计算机科学家努力的方向,也是未来智能应用实现的基础技术。语音同样包含大量的数据信息,识别语音的过程也是对这些信息处理的过程。NMF算法在这方面也为我们提供了一种新方法,在已有的应用中,NMF算法成功实现了有效的语音特征提取,并且由于NMF算法的快速性,对实现机器的实时语音识别有着促进意义。也有使用NMF方法进行音乐分析的应用。复调音乐的识别是个很困难的问题,三菱研究所和MIT(麻省理工学院)的科学家合作,利用NMF从演奏中的复调音乐中识别出各个调子,并将它们分别记录下来。实验结果表明,这种采用NMF算法的方法不光简单,而且无须基于知识库。
      (4) 机器人控制
      如何快速准确地让机器人识别周围的物体对于机器人研究具有重要的意义,因为这是机器人能迅速作出相应反应和动作的基础。机器人通过传感器获得周围环境的图像信息,这些图像信息也是以矩阵的形式存储的。已经有研究人员采用NMF算法实现了机器人对周围对象的快速识别,根据现有的研究资料显示,识别的准确率达到了80%以上。
      (5) 生物医学工程和化学工程
      生物医学和化学研究中,也常常需要借助计算机来分析处理试验的数据,往往一些烦杂的数据会耗费研究人员的过多精力。NMF算法也为这些数据的处理提供了一种新的高效快速的途径。科学家将NMF方法用于处理核医学中的电子发射过程的动态连续图像,有效地从这些动态图像中提取所需要的特征。NMF还可以应用到遗传学和药物发现中。因为NMF的分解不出现负值,因此采用NMF分析基因DNA的分子序列可使分析结果更加可靠。同样,用NMF来选择药物成分还可以获得最有效的且负作用最小的新药物。[非负矩阵分解在分析癌症突变异质性中的作用]

      此外,NMF算法在环境数据处理、信号分析与复杂对象的识别方面都有着很好的应用。近年来采用NMF思想的应用才刚展开,相信以后会有更多的成功应用。这些成功的应用反过来也将促进NMF的进一步研究。
      三、结束语
      数学如同计算机的灵魂,NMF通过计算机与各个领域结合后的应用取得了令人叹服的成效。NMF的故事还在继续,NMF的应用领域还有待进一步的发掘,针对NMF的进一步研究也没有停止过,其中诸如分解的存在性、惟一性和收敛性以及收敛的速度等问题的深入探讨必将使该思想能更好地服务于人类。

    [非负矩阵分解:数学的奇妙力量]

    NMF的应用实例

    示例1

    无论是在需要将新闻分门别类的应用中还是在信息检索中判断两个网页文本内容的相似度(从而判断是否同类)的应用中,我们都需要从文档中学习到一些主要的特征词以代表这篇文档并据此聚类。

    有许多应用都需要将模式(文档,图像,声音)进行压缩或者概括,或者说,我们需要从这些模式中提取出少量的重要的特征,原来的模式可以用这些特征进行还原而不出现大的误差。将矩阵分解为两个低秩矩阵相乘提供一种解决此问题的思路,两个矩阵相乘的结果是原矩阵低维的近似。

    在图像领域,需要将图像压缩表达为一些特征的组合。

    实战中可能多次实验结果不同,需要强调这是局部最优,结果受初始值选取的影响。

    示例2

    请设计合适的主成分分析模型, 并提供相应算法, 使得当输⼊为下图左⼿边信息时, 可以返回右⼿边的主成分元素。

    左图是通过右图8个部分通过一定的随机稀疏组合而成的8×8=64幅飞机图,我们要做的就是给你左图,求出右图的成分。

    这时使用NMF就可以轻松求解:(求解出的components向量可能是归一化的,而原图是01组成的,我们可以通过阈值来将浮点数变成01)

    estimators = decomposition.NMF(n_components=n_components, init='nndsvd', )
    estimator.fit(data)
    components_ = estimator.components_[:n_components]  # numpy.ndarray
    plot_gallery(components_)

    components_ = components_[sort_index]
    plot_gallery('criterias', criterias, image_shape=image_shape)
    components_[components_ > 0.1] = 1
    components_[components_ <= 0.1] = 0
    plot_gallery('%s - Train time %.1fs' % (name, train_time), components_, image_shape=image_shape)
    error = spatial.distance.euclidean(criterias.reshape(-1), components_.reshape(-1))

    不过如果使用下面高级主题中提到的Structured PCA,可能效果更好。

    [LiuXin: 大数据优化算法]

    皮皮blog



    非负矩阵分解的算法和实现

    NMF实现原理:NMF若干更新法则

    NMF求解问题实际上是一个最优化问题,利用乘性迭代的方法求解,非负矩阵分解是一个NP问题。NMF问题的目标函数有很多种,应用最广泛的就是欧几里得距离和KL散度。

    关于更新法则,Daniel D. Lee和H. Sebastian Seung的文章《Algorithms for Non-negative Matrix Factorization》有详细的公式推导证明。由于W与H的乘积是对V的近似估计,所以评价分解好坏的标准便是两者之间的差异。文中在不同的代价函数下提出了不同的更新法则,包括乘性更新法则与加性更新法则。文中还通过构造辅助函数对迭代算法的收敛性进行了证明。

    写法1

       在NMF的分解问题中,假设噪声矩阵为,那么有

                    

       现在要找出合适的使得最小。

    假设噪声服从不同的概率分布,通过最大似然函数会得到不同类型的目标函数。接下来会分别以噪声服从高斯分布和泊松分布来说明。

        (1)噪声服从高斯分布

           假设噪声服从高斯分布,那么得到最大似然函数为

          

            取对数后,得到对数似然函数为

           

            假设各数据点噪声的方差一样,那么接下来要使得对数似然函数取值最大,只需要下面目标函数值最小。

           

           该损失函数为2范数损失函数,它是基于欧几里得距离的度量。又因为

           

           那么得到

           

           同理有

           

            接下来就可以使用梯度下降法进行迭代了。如下

           

           如果选取

           

            那么最终得到迭代式为

           

            可看出这是乘性迭代规则,每一步都保证了结果为正数,一直迭代下去就会收敛,当然收敛性的证明省略。

       (2)噪声服从泊松分布

           若噪声为泊松噪声,那么得到损失函数为

           

           同样经过推到得到

           

    [非负矩阵分解(NMF) ]

    写法2


    •  Cost function 1,欧几里得距离( Euclidean distance):

                                                                                                    (1)

               在此代价函数下的乘性更新法则有:

           及                                       (2)
     

    • cost function 2,分离度(divergence):

                                                             (3)

                  在此代价函数下的乘性更新法则为:

             以及                                         (4)


           另外有加性更新法则如下:

           欧几里得距离( Euclidean distance)下的加性更新为:

                                       (5)

     此时只要迭代步长的所有元素设为足够小的正数那么欧式距离会随着迭代而减小。而当迭代步长满足:,可以得到此时的加性更新等价于式子(2)中的乘性更新。

            分离度(divergence)下的加性更新为:

                                                    (6)

    同样迭代步长若为足够下的正数的话,分离度代价函数的值会随着更新而减小,而当迭代步长满足:,可同样得到此时的加性更新等价于式子(4)中的乘性更新法则。

    [ 非负矩阵分解(NMF,Nonnegtive Matrix Factorization) ]

    收敛性证明

    如上所述,《Algorithms for Non-negative Matrix Factorization》文中还通过构造辅助函数对迭代算法的收敛性进行了证明。小小整理如下:


    [凯:非负矩阵分解Algorithms for Non-negative Matrix Factorization]

    NMF的相关算法

    Local NMF(LNMF),Discriminant NMF(DNMF),non-negative sparse coding(NNSC),non-smooth NMF(nsNMF),projective NMF,temporal NMF with spatial decorrelation constraints,shifted NMF,incremental NMF,sparse higher order NMF and polynomial NMF 等等,对于NMF分解的唯一性探讨也是NMF最近比较受关注的方面。

    非负矩阵分解已经有了很多算法,例如Multiplicative   Update,Projected Gradient Method,Gradient Descent。交替最小二乘法比较符合上面的解释,具有更好的统计意义,而且收敛性质很好,计算速度快。已有的交替最小二乘法都使用冷冰冰的二次规划算法求解非负回归,在bignmf中使用的最小角回归的思路解非负回归,更有统计的感觉,而且速度更快。

     非负矩阵分解的具体算法

    输入参数:X,R,MAXITER,其中X为被分解的矩阵,R为降阶后B的秩,MAXITER为迭代次数
    输出参数:B,H
    1):初始化矩阵B,H为非负数,同时对B的每一列数据归一化
    2):for i=1:MAXITER
        a:更新H矩阵一行元素:H(i,j)=H(i,j)*(B'*X)(i,j)/(B'*B*H)(i,j)
        b:更新B的一列元素:B(k,j)=B(k,j)*(X*H')(k,j)/(B*H*H')(k,j);
        c: 重新对B进行列归一化
    3)end

    非负矩阵分解的实现代码

    matlab源程序如下:

    dim=size(X);                                    %计算x的规格
    X=double(X);
    B=10*rand(dim(1),r);                            %初始化BH,为非负数
    B=B./(ones(dim(1),1)*sum(B));                   %归一化B的每一列

    H=10*rand(r,dim(2));

    maxiter=100;                                    %最大迭代次数
    for iter=1:maxiter
        H=H.*(B'*(X./(B*H)));
        B=B.*((X./(B*H))*H');
        B=B./(ones(dim(1),1)*sum(B));
    end

    [matlab中的NMF函数及其使用说明NMFfuction- Mathworks]

    [Chih-Jen Lin关于NMF的一些源代码 NMF工具

    MATLAB Scripts download nmf.m.

    Python: nmf.py by Anthony Di Franco (numpy needed; see an example.py here).

    Go: nmf.go by Dan Kortschak and examples.

    C: NMF.c by Dong Li]

    [NMF算法简介及python实现]

    [NMF的R实现:bignmf:https://github.com/panlanfeng/bignmf]

    [NMF的Julia实现:https://github.com/JuliaStats/NMF.jl]

    皮皮blog



    非负矩阵分解相关的其它高级主题

    [NMF与pLSA的对比《Non-negative Matrix Factorization and Probabilistic Latent Semantic Analysis》]

    [除了NMF分解以外,还有一种非负分解也用的比较多,而且很有效,叫做非负张量分解法:http://www.doc88.com/p-8942237517189.html]

    [浅析基于张量分解的隐变量模型参数学习]

    [Structured PCA带结构的主成分分析:主成分分析是找到给定数据的结构信息的主要手段之一,然而对于某些特殊的输入数据传统的主成分分析, 并不能真正得到问题的主成分, 这时我们在建模时就需要考虑到这些特殊的结构。

    NMF分解和SSPAC分解在主成分较小的数据上,表现similarity,在主成分较多的数据上,SSPAC表现较好。(Jenatton R, Obozinski G, Bach F R. Structured Sparse Principal Component Analysis[C]//AISTATS. 2010: 366-373)

    Structured Variable Selection with Sparsity-Inducing Norms

    Structured  Principal Component Analysis.tar.gz ]

    from: http://blog.csdn.net/pipisorry/article/details/52098864

    ref: [wikipedia: Non-negative matrix factorization]

    [1] Daniel D lee and H.SebastianSeung(1999).”Learning the parts of objects by non-negative matrix factorization”. Nature
    [2] Daniel D.Leeand H.SebastianSeung(2001). Algorithms for Non negative Matrix Factorization.Advancesin Neural Information
    [3] Lee D D, Seung H S. Learning the parts of objects by non-negative matrix factorization[J]. Nature, 1999, 401(6755): 788-791.
    [4] Bach F, Jenatton R, Obozinski G. Structured sparse principal component analysis[J]. J. Mach. Learn. Res, 366-373.
    [5] Jenatton R, Obozinski G, Bach F R. Structured Sparse Principal Component Analysis[C]//AISTATS. 2010: 366-373.
    [6] python sklearn http://scikitlearn.org/stable/modules/decomposition.html#pca

    展开全文
  • 支持向量机通俗导论(理解SVM的三层境界)

    万次阅读 多人点赞 2012-06-01 22:48:43
    支持向量机通俗导论(理解SVM的三层境界) 作者:July 。致谢:pluskid、白石、JerryLead。 说明:本文最初写于2012年6月,而后不断反反复复修改&amp;优化,修改次数达上百次,最后修改于2016年11月。 声明:...

                支持向量机通俗导论(理解SVM的三层境界)

    作者:July 。致谢:pluskid、白石、JerryLead。
    说明:本文最初写于2012年6月,而后不断反反复复修改&优化,修改次数达上百次,最后修改于2016年11月。
    声明:本文于2012年便早已附上所有参考链接,并注明是篇“学习笔记”,且写明具体参考了pluskid等人的文章。文末2013年的PDF是为证。另,如果本文的图片/公式无法正常显示,请点击本文的七月在线题库版

    前言

        动笔写这个支持向量机(support vector machine)是费了不少劲和困难的,原因很简单,一者这个东西本身就并不好懂,要深入学习和研究下去需花费不少时间和精力,二者这个东西也不好讲清楚,尽管网上已经有朋友写得不错了(见文末参考链接),但在描述数学公式的时候还是显得不够。得益于同学白石的数学证明,我还是想尝试写一下,希望本文在兼顾通俗易懂的基础上,真真正正能足以成为一篇完整概括和介绍支持向量机的导论性的文章。

        本文在写的过程中,参考了不少资料,包括《支持向量机导论》、《统计学习方法》及网友pluskid的支持向量机系列等等,于此,还是一篇学习笔记,只是加入了自己的理解和总结,有任何不妥之处,还望海涵。全文宏观上整体认识支持向量机的概念和用处,微观上深究部分定理的来龙去脉,证明及原理细节,力保逻辑清晰 & 通俗易懂。

        同时,阅读本文时建议大家尽量使用chrome等浏览器,如此公式才能更好的显示,再者,阅读时可拿张纸和笔出来,把本文所有定理.公式都亲自推导一遍或者直接打印下来(可直接打印网页版或本文文末附的PDF)在文稿上演算,从而享受随时随地思考、演算的极致快感。

        OK,还是那句话,有任何问题,欢迎任何人随时不吝指正 & 赐教,感谢。

    第一层、了解SVM

        支持向量机,因其英文名为support vector machine,故一般简称SVM,通俗来讲,它是一种二类分类模型,其基本模型定义为特征空间上的间隔最大的线性分类器,其学习策略便是间隔最大化,最终可转化为一个凸二次规划问题的求解。

    1.1、分类标准的起源:Logistic回归

        理解SVM,咱们必须先弄清楚一个概念:线性分类器。

        给定一些数据点,它们分别属于两个不同的类,现在要找到一个线性分类器把这些数据分成两类。如果用x表示数据点,用y表示类别(y可以取1或者-1,分别代表两个不同的类),一个线性分类器的学习目标便是要在n维的数据空间中找到一个超平面(hyper plane),这个超平面的方程可以表示为( wT中的T代表转置):

                                                               

        可能有读者对类别取1或-1有疑问,事实上,这个1或-1的分类标准起源于logistic回归

        Logistic回归目的是从特征学习出一个0/1分类模型,而这个模型是将特性的线性组合作为自变量,由于自变量的取值范围是负无穷到正无穷。因此,使用logistic函数(或称作sigmoid函数)将自变量映射到(0,1)上,映射后的值被认为是属于y=1的概率。

        假设函数

        其中x是n维特征向量,函数g就是logistic函数。

        而的图像是

        可以看到,将无穷映射到了(0,1)。

        而假设函数就是特征属于y=1的概率。

        从而,当我们要判别一个新来的特征属于哪个类时,只需求即可,若大于0.5就是y=1的类,反之属于y=0类。

        此外,只和有关,>0,那么,而g(z)只是用来映射,真实的类别决定权还是在于。再者,当时,=1,反之=0。如果我们只从出发,希望模型达到的目标就是让训练数据中y=1的特征,而是y=0的特征。Logistic回归就是要学习得到,使得正例的特征远大于0,负例的特征远小于0,而且要在全部训练实例上达到这个目标。

        接下来,尝试把logistic回归做个变形。首先,将使用的结果标签y = 0和y = 1替换为y = -1,y = 1,然后将)中的替换为b,最后将后面的替换为(即)。如此,则有了。也就是说除了y由y=0变为y=-1外,线性分类函数跟logistic回归的形式化表示没区别。

        进一步,可以将假设函数中的g(z)做一个简化,将其简单映射到y=-1和y=1上。映射关系如下:

    1.2、线性分类的一个例子

        下面举个简单的例子。如下图所示,现在有一个二维平面,平面上有两种不同的数据,分别用圈和叉表示。由于这些数据是线性可分的,所以可以用一条直线将这两类数据分开,这条直线就相当于一个超平面,超平面一边的数据点所对应的y全是-1 ,另一边所对应的y全是1。

        这个超平面可以用分类函数表示,当f(x) 等于0的时候,x便是位于超平面上的点,而f(x)大于0的点对应 y=1 的数据点,f(x)小于0的点对应y=-1的点,如下图所示:

        注:有的资料上定义特征到结果的输出函数与这里定义的实质是一样的。为什么?因为无论是,还是,不影响最终优化结果。下文你将看到,当我们转化到优化的时候,为了求解方便,会把yf(x)令为1,即yf(x)是y(w^x + b),还是y(w^x - b),对我们要优化的式子max1/||w||已无影响。

        (有一朋友飞狗来自Mare_Desiderii,看了上面的定义之后,问道:请教一下SVM functional margin 为=y(wTx+b)=yf(x)中的Y是只取1和-1 吗?y的唯一作用就是确保functional margin的非负性?真是这样的么?当然不是,详情请见本文评论下第43楼

        换言之,在进行分类的时候,遇到一个新的数据点x,将x代入f(x) 中,如果f(x)小于0则将x的类别赋为-1,如果f(x)大于0则将x的类别赋为1。

        接下来的问题是,如何确定这个超平面呢?从直观上而言,这个超平面应该是最适合分开两类数据的直线。而判定“最适合”的标准就是这条直线离直线两边的数据的间隔最大。所以,得寻找有着最大间隔的超平面。

    1.3、函数间隔Functional margin与几何间隔Geometrical margin 

        在超平面w*x+b=0确定的情况下,|w*x+b|能够表示点x到距离超平面的远近,而通过观察w*x+b的符号与类标记y的符号是否一致可判断分类是否正确,所以,可以用(y*(w*x+b))的正负性来判定或表示分类的正确性。于此,我们便引出了函数间隔(functional margin)的概念。

        定义函数间隔(用表示)为:

        而超平面(w,b)关于T中所有样本点(xi,yi)的函数间隔最小值(其中,x是特征,y是结果标签,i表示第i个样本),便为超平面(w, b)关于训练数据集T的函数间隔:

         mini  (i=1,...n)

        但这样定义的函数间隔有问题,即如果成比例的改变w和b(如将它们改成2w和2b),则函数间隔的值f(x)却变成了原来的2倍(虽然此时超平面没有改变),所以只有函数间隔还远远不够。

        事实上,我们可以对法向量w加些约束条件,从而引出真正定义点到超平面的距离--几何间隔(geometrical margin)的概念。

        假定对于一个点 x ,令其垂直投影到超平面上的对应点为 x0 ,w 是垂直于超平面的一个向量,为样本x到超平面的距离,如下图所示:

    根据平面几何知识,有

        其中||w||为w的二阶范数(范数是一个类似于模的表示长度的概念),是单位向量(一个向量除以它的模称之为单位向量)。

        又由于x0 是超平面上的点,满足 f(x0)=0,代入超平面的方程,可得,即

        随即让此式的两边同时乘以,再根据,即可算出γ: 

      为了得到的绝对值,令乘上对应的类别 y,即可得出几何间隔(用表示)的定义:

        从上述函数间隔和几何间隔的定义可以看出:几何间隔就是函数间隔除以||w||,而且函数间隔y*(wx+b) = y*f(x)实际上就是|f(x)|,只是人为定义的一个间隔度量,而几何间隔|f(x)|/||w||才是直观上的点到超平面的距离。

    1.4、最大间隔分类器Maximum Margin Classifier的定义

        对一个数据点进行分类,当超平面离数据点的“间隔”越大,分类的确信度(confidence)也越大。所以,为了使得分类的确信度尽量高,需要让所选择的超平面能够最大化这个“间隔”值。这个间隔就是下图中的Gap的一半。

        通过由前面的分析可知:函数间隔不适合用来最大化间隔值,因为在超平面固定以后,可以等比例地缩放w的长度和b的值,这样可以使得的值任意大,亦即函数间隔可以在超平面保持不变的情况下被取得任意大。但几何间隔因为除上了,使得在缩放w和b的时候几何间隔的值是不会改变的,它只随着超平面的变动而变动,因此,这是更加合适的一个间隔。换言之,这里要找的最大间隔分类超平面中的“间隔”指的是几何间隔。

       于是最大间隔分类器(maximum margin classifier的目标函数可以定义为:

        同时需满足一些条件,根据间隔的定义,有

        其中,s.t.,即subject to的意思,它导出的是约束条件。

        回顾下几何间隔的定义,可知:如果令函数间隔等于1(之所以令等于1,是为了方便推导和优化,且这样做对目标函数的优化没有影响,至于为什么,请见本文评论下第42楼回复),则有 = 1 / ||w||且,从而上述目标函数转化成了

        相当于在相应的约束条件下,最大化这个1/||w||值,而1/||w||便是几何间隔。   

        如下图所示,中间的实线便是寻找到的最优超平面(Optimal Hyper Plane),其到两条虚线边界的距离相等,这个距离便是几何间隔,两条虚线间隔边界之间的距离等于2,而虚线间隔边界上的点则是支持向量。由于这些支持向量刚好在虚线间隔边界上,所以它们满足(还记得我们把 functional margin 定为 1 了吗?上节中:处于方便推导和优化的目的,我们可以令=1),而对于所有不是支持向量的点,则显然有

        OK,到此为止,算是了解到了SVM的第一层,对于那些只关心怎么用SVM的朋友便已足够,不必再更进一层深究其更深的原理。

    第二层、深入SVM

    2.1、从线性可分到线性不可分

    2.1.1、从原始问题到对偶问题的求解

        接着考虑之前得到的目标函数:

         由于求的最大值相当于求的最小值,所以上述目标函数等价于(w由分母变成分子,从而也有原来的max问题变为min问题,很明显,两者问题等价):

        因为现在的目标函数是二次的,约束条件是线性的,所以它是一个凸二次规划问题。这个问题可以用现成的QP (Quadratic Programming) 优化包进行求解。一言以蔽之:在一定的约束条件下,目标最优,损失最小。

        此外,由于这个问题的特殊结构,还可以通过拉格朗日对偶性(Lagrange Duality)变换到对偶变量 (dual variable) 的优化问题,即通过求解与原问题等价的对偶问题(dual problem)得到原始问题的最优解,这就是线性可分条件下支持向量机的对偶算法,这样做的优点在于:一者对偶问题往往更容易求解;二者可以自然的引入核函数,进而推广到非线性分类问题。

         那什么是拉格朗日对偶性呢?简单来讲,通过给每一个约束条件加上一个拉格朗日乘子(Lagrange multiplier),定义拉格朗日函数(通过拉格朗日函数将约束条件融合到目标函数里去,从而只用一个函数表达式便能清楚的表达出我们的问题):

        然后令

        容易验证,当某个约束条件不满足时,例如,那么显然有(只要令即可)。而当所有约束条件都满足时,则最优值为,亦即最初要最小化的量。

        因此,在要求约束条件得到满足的情况下最小化,实际上等价于直接最小化(当然,这里也有约束条件,就是≥0,i=1,…,n)   ,因为如果约束条件没有得到满足,会等于无穷大,自然不会是我们所要求的最小值。

        具体写出来,目标函数变成了:

        这里用表示这个问题的最优值,且和最初的问题是等价的。如果直接求解,那么一上来便得面对w和b两个参数,而又是不等式约束,这个求解过程不好做。不妨把最小和最大的位置交换一下,变成:

        交换以后的新问题是原始问题的对偶问题,这个新问题的最优值用来表示。而且有,在满足某些条件的情况下,这两者相等,这个时候就可以通过求解对偶问题来间接地求解原始问题。

        换言之,之所以从minmax的原始问题,转化为maxmin的对偶问题,一者因为的近似解,二者,转化为对偶问题后,更容易求解。

        下面可以先求L 对w、b的极小,再求L 对的极大。

    2.1.2、KKT条件

        上文中提到“在满足某些条件的情况下,两者等价”,这所谓的“满足某些条件”就是要满足KKT条件。

      勘误:经读者qq_28543029指出,这里的条件不应该是KKT条件,要让两者等价需满足strong duality (强对偶),而后有学者在强对偶下提出了KKT条件,且KKT条件的成立要满足constraint qualifications,而constraint qualifications之一就是Slater条件。所谓Slater 条件,即指:凸优化问题,如果存在一个点x,使得所有等式约束都成立,并且所有不等式约束都严格成立(即取严格不等号,而非等号),则满足Slater 条件。对于此处,Slater 条件成立,所以d*≤p*可以取等号。

        一般地,一个最优化数学模型能够表示成下列标准形式:

        其中,f(x)是需要最小化的函数,h(x)是等式约束,g(x)是不等式约束,p和q分别为等式约束和不等式约束的数量。

        同时,得明白以下两点:

    • 凸优化的概念:\mathcal{X} \subset \mathbb{R}^n 为一凸集, f:\mathcal{X}\to \mathbb{R} 为一凸函数。凸优化就是要找出一点,使得每一 x \in \mathcal{X} 满足 。
    • KKT条件的意义:它是一个非线性规划(Nonlinear Programming)问题能有最优化解法的必要和充分条件。

        而KKT条件就是指上面最优化数学模型的标准形式中的最小点 x* 必须满足下面的条件:

        经过论证,我们这里的问题是满足 KKT 条件的(首先已经满足Slater条件,再者f和gi也都是可微的,即L对w和b都可导),因此现在我们便转化为求解第二个问题。

        也就是说,原始问题通过满足KKT条件,已经转化成了对偶问题。而求解这个对偶学习问题,分为3个步骤:首先要让L(w,b,a) 关于 w 和 b 最小化,然后求对的极大,最后利用SMO算法求解对偶问题中的拉格朗日乘子。

    2.1.3、对偶问题求解的3个步骤

        (1)、首先固定要让 L 关于 w 和 b 最小化,我们分别对w,b求偏导数,即令 ∂L/∂w 和 ∂L/∂b 等于零(对w求导结果的解释请看本文评论下第45楼回复):

        将以上结果代入之前的L: 

        得到:

        提醒:有读者可能会问上述推导过程如何而来?说实话,其具体推导过程是比较复杂的,如下图所示:

          最后,得到:

        如 jerrylead所说:“倒数第4步”推导到“倒数第3步”使用了线性代数的转置运算,由于ai和yi都是实数,因此转置后与自身一样。“倒数第3步”推导到“倒数第2步”使用了(a+b+c+…)(a+b+c+…)=aa+ab+ac+ba+bb+bc+…的乘法运算法则。最后一步是上一步的顺序调整。

        从上面的最后一个式子,我们可以看出,此时的拉格朗日函数只包含了一个变量,那就是(求出了便能求出w,和b,由此可见,上文第1.2节提出来的核心问题:分类函数也就可以轻而易举的求出来了)。

        (2)、求对的极大,即是关于对偶问题的最优化问题。经过上面第一个步骤的求w和b,得到的拉格朗日函数式子已经没有了变量w,b,只有。从上面的式子得到:

        这样,求出了,根据,即可求出w,然后通过,即可求出b,最终得出分离超平面和分类决策函数。

        (3)在求得L(w, b, a) 关于 w 和 b 最小化,以及对的极大之后,最后一步则可以利用SMO算法求解对偶问题中的拉格朗日乘子

           上述式子要解决的是在参数上求最大值W的问题,至于都是已知数要了解这个SMO算法是如何推导的,请跳到下文第3.5节、SMO算法。

        到目前为止,我们的 SVM 还比较弱,只能处理线性的情况,下面我们将引入核函数,进而推广到非线性分类问题。

    2.1.4、线性不可分的情况

        OK,为过渡到下节2.2节所介绍的核函数,让我们再来看看上述推导过程中得到的一些有趣的形式。首先就是关于我们的 hyper plane ,对于一个数据点 x 进行分类,实际上是通过把 x 带入到算出结果然后根据其正负号来进行类别划分的。而前面的推导中我们得到 

        因此分类函数为:

        这里的形式的有趣之处在于,对于新点 x的预测,只需要计算它与训练数据点的内积即可(表示向量内积),这一点至关重要,是之后使用 Kernel 进行非线性推广的基本前提。此外,所谓 Supporting Vector 也在这里显示出来——事实上,所有非Supporting Vector 所对应的系数都是等于零的,因此对于新点的内积计算实际上只要针对少量的“支持向量”而不是所有的训练数据即可。

        为什么非支持向量对应的等于零呢?直观上来理解的话,就是这些“后方”的点——正如我们之前分析过的一样,对超平面是没有影响的,由于分类完全有超平面决定,所以这些无关的点并不会参与分类问题的计算,因而也就不会产生任何影响了。

        回忆一下我们2.1.1节中通过 Lagrange multiplier得到的目标函数:

         注意到如果 xi 是支持向量的话,上式中红颜色的部分是等于 0 的(因为支持向量的 functional margin 等于 1 ),而对于非支持向量来说,functional margin 会大于 1 ,因此红颜色部分是大于零的,而又是非负的,为了满足最大化,必须等于 0 。这也就是这些非Supporting Vector 的点的局限性。 

        至此,我们便得到了一个maximum margin hyper plane classifier,这就是所谓的支持向量机(Support Vector Machine)。当然,到目前为止,我们的 SVM 还比较弱,只能处理线性的情况,不过,在得到了对偶dual 形式之后,通过 Kernel 推广到非线性的情况就变成了一件非常容易的事情了(相信,你还记得本节开头所说的:“通过求解对偶问题得到最优解,这就是线性可分条件下支持向量机的对偶算法,这样做的优点在于:一者对偶问题往往更容易求解;二者可以自然的引入核函数,进而推广到非线性分类问题”)。

    2.2、核函数Kernel

    2.2.1、特征空间的隐式映射:核函数

        事实上,大部分时候数据并不是线性可分的,这个时候满足这样条件的超平面就根本不存在。在上文中,我们已经了解到了SVM处理线性可分的情况,那对于非线性的数据SVM咋处理呢?对于非线性的情况,SVM 的处理方法是选择一个核函数 κ(⋅,⋅) ,通过将数据映射到高维空间,来解决在原始空间中线性不可分的问题。

        具体来说,在线性不可分的情况下,支持向量机首先在低维空间中完成计算,然后通过核函数将输入空间映射到高维特征空间,最终在高维特征空间中构造出最优分离超平面,从而把平面上本身不好分的非线性数据分开。如图所示,一堆数据在二维空间无法划分,从而映射到三维空间里划分:

        而在我们遇到核函数之前,如果用原始的方法,那么在用线性学习器学习一个非线性关系,需要选择一个非线性特征集,并且将数据写成新的表达形式,这等价于应用一个固定的非线性映射,将数据映射到特征空间,在特征空间中使用线性学习器,因此,考虑的假设集是这种类型的函数:

        这里ϕ:X->F是从输入空间到某个特征空间的映射,这意味着建立非线性学习器分为两步:

    1. 首先使用一个非线性映射将数据变换到一个特征空间F,
    2. 然后在特征空间使用线性学习器分类。

        而由于对偶形式就是线性学习器的一个重要性质,这意味着假设可以表达为训练点的线性组合,因此决策规则可以用测试点和训练点的内积来表示:

        如果有一种方式可以在特征空间中直接计算内积〈φ(xi · φ(x),就像在原始输入点的函数中一样,就有可能将两个步骤融合到一起建立一个非线性的学习器,这样直接计算法的方法称为核函数方法:

        核是一个函数K,对所有x,z(-X,满足,这里φ是从X到内积特征空间F的映射。

    2.2.2、核函数:如何处理非线性数据

        来看个核函数的例子。如下图所示的两类数据,分别分布为两个圆圈的形状,这样的数据本身就是线性不可分的,此时咱们该如何把这两类数据分开呢(下文将会有一个相应的三维空间图)?

        事实上,上图所述的这个数据集,是用两个半径不同的圆圈加上了少量的噪音生成得到的,所以,一个理想的分界应该是一个“圆圈”而不是一条线(超平面)。如果用来表示这个二维平面的两个坐标的话,我们知道一条二次曲线(圆圈是二次曲线的一种特殊情况)的方程可以写作这样的形式:

        注意上面的形式,如果我们构造另外一个五维的空间,其中五个坐标的值分别为, , , ,那么显然,上面的方程在新的坐标系下可以写作:

        关于新的坐标,这正是一个 hyper plane 的方程!也就是说,如果我们做一个映射,将 按照上面的规则映射为,那么在新的空间中原来的数据将变成线性可分的,从而使用之前我们推导的线性分类算法就可以进行处理了。这正是 Kernel 方法处理非线性问题的基本思想。

        再进一步描述 Kernel 的细节之前,不妨再来看看上述例子在映射过后的直观形态。当然,你我可能无法把 5 维空间画出来,不过由于我这里生成数据的时候用了特殊的情形,所以这里的超平面实际的方程是这个样子的(圆心在轴上的一个正圆):

        因此我只需要把它映射到这样一个三维空间中即可,下图即是映射之后的结果,将坐标轴经过适当的旋转,就可以很明显地看出,数据是可以通过一个平面来分开的(pluskid:下面的gif 动画,先用 Matlab 画出一张张图片,再用 Imagemagick 拼贴成):

        核函数相当于把原来的分类函数:

        映射成:

        而其中的可以通过求解如下 dual 问题而得到的:

        这样一来问题就解决了吗?似乎是的:拿到非线性数据,就找一个映射,然后一股脑把原来的数据映射到新空间中,再做线性 SVM 即可。不过事实上好像并没有这么简单。

        细想一下,刚才的方法是不是有问题?

    • 在最初的例子里,我们对一个二维空间做映射,选择的新空间是原始空间的所有一阶和二阶的组合,得到了五个维度;
    • 如果原始空间是三维(一阶、二阶和三阶的组合),那么我们会得到:3(一次)+3(二次交叉)+3(平方)+3(立方)+1(x1*x2*x3)+2*3(交叉,一个一次一个二次,类似x1*x2^2) = 19维的新空间,这个数目是呈指数级爆炸性增长的,从而势必这给的计算带来非常大的困难,而且如果遇到无穷维的情况,就根本无从计算了。

        这个时候,可能就需要 Kernel 出马了。

        不妨还是从最开始的简单例子出发,设两个向量,而即是到前面说的五维空间的映射,因此映射过后的内积为:

            (公式说明:上面的这两个推导过程中,所说的前面的五维空间的映射,这里说的前面便是文中2.2.1节的所述的映射方式,回顾下之前的映射规则,再看那第一个推导,其实就是计算x1,x2各自的内积,然后相乘相加即可,第二个推导则是直接平方,去掉括号,也很容易推出来)

        另外,我们又注意到:

         二者有很多相似的地方,实际上,我们只要把某几个维度线性缩放一下,然后再加上一个常数维度,具体来说,上面这个式子的计算结果实际上和映射

         之后的内积的结果是相等的,那么区别在于什么地方呢?

    1. 一个是映射到高维空间中,然后再根据内积的公式进行计算;
    2. 而另一个则直接在原来的低维空间中进行计算,而不需要显式地写出映射后的结果

        (公式说明:上面之中,最后的两个式子,第一个算式,是带内积的完全平方式,可以拆开,然后,通过凑一个得到,第二个算式,也是根据第一个算式凑出来的)

        回忆刚才提到的映射的维度爆炸,在前一种方法已经无法计算的情况下,后一种方法却依旧能从容处理,甚至是无穷维度的情况也没有问题。

        我们把这里的计算两个向量在隐式映射过后的空间中的内积的函数叫做核函数 (Kernel Function) ,例如,在刚才的例子中,我们的核函数为:

        核函数能简化映射空间中的内积运算——刚好“碰巧”的是,在我们的 SVM 里需要计算的地方数据向量总是以内积的形式出现的。对比刚才我们上面写出来的式子,现在我们的分类函数为:

        其中 由如下 dual 问题计算而得:

        这样一来计算的问题就算解决了,避开了直接在高维空间中进行计算,而结果却是等价的!当然,因为我们这里的例子非常简单,所以我可以手工构造出对应于的核函数出来,如果对于任意一个映射,想要构造出对应的核函数就很困难了。

    2.2.3、几个核函数

        通常人们会从一些常用的核函数中选择(根据问题和数据的不同,选择不同的参数,实际上就是得到了不同的核函数),例如:

    • 多项式核,显然刚才我们举的例子是这里多项式核的一个特例(R = 1,d = 2)。虽然比较麻烦,而且没有必要,不过这个核所对应的映射实际上是可以写出来的,该空间的维度是,其中 m 是原始空间的维度。
    • 高斯核,这个核就是最开始提到过的会将原始空间映射为无穷维空间的那个家伙。不过,如果选得很大的话,高次特征上的权重实际上衰减得非常快,所以实际上(数值上近似一下)相当于一个低维的子空间;反过来,如果选得很小,则可以将任意的数据映射为线性可分——当然,这并不一定是好事,因为随之而来的可能是非常严重的过拟合问题。不过,总的来说,通过调控参数,高斯核实际上具有相当高的灵活性,也是使用最广泛的核函数之一。下图所示的例子便是把低维线性不可分的数据通过高斯核函数映射到了高维空间:
    • 线性核,这实际上就是原始空间中的内积。这个核存在的主要目的是使得“映射后空间中的问题”和“映射前空间中的问题”两者在形式上统一起来了(意思是说,咱们有的时候,写代码,或写公式的时候,只要写个模板或通用表达式,然后再代入不同的核,便可以了,于此,便在形式上统一了起来,不用再分别写一个线性的,和一个非线性的)。

    2.2.4、核函数的本质

            上面说了这么一大堆,读者可能还是没明白核函数到底是个什么东西?我再简要概括下,即以下三点:

    1. 实际中,我们会经常遇到线性不可分的样例,此时,我们的常用做法是把样例特征映射到高维空间中去(如上文2.2节最开始的那幅图所示,映射到高维空间后,相关特征便被分开了,也就达到了分类的目的);
    2. 但进一步,如果凡是遇到线性不可分的样例,一律映射到高维空间,那么这个维度大小是会高到可怕的(如上文中19维乃至无穷维的例子)。那咋办呢?
    3. 此时,核函数就隆重登场了,核函数的价值在于它虽然也是将特征进行从低维到高维的转换,但核函数绝就绝在它事先在低维上进行计算,而将实质上的分类效果表现在了高维上,也就如上文所说的避免了直接在高维空间中的复杂计算。

        最后引用这里的一个例子举例说明下核函数解决非线性问题的直观效果。

        假设现在你是一个农场主,圈养了一批牛群,但为预防狼群袭击牛群,你需要搭建一个篱笆来把牛群围起来。但是篱笆应该建在哪里呢?你很可能需要依据牛群和狼群的位置建立一个“分类器”,比较下图这几种不同的分类器,我们可以看到SVM完成了一个很完美的解决方案。

        这个例子从侧面简单说明了SVM使用非线性分类器的优势,而逻辑模式以及决策树模式都是使用了直线方法。

        OK,不再做过多介绍了,对核函数有进一步兴趣的,还可以看看此文

    2.3、使用松弛变量处理 outliers 方法

        在本文第一节最开始讨论支持向量机的时候,我们就假定,数据是线性可分的,亦即我们可以找到一个可行的超平面将数据完全分开。后来为了处理非线性数据,在上文2.2节使用 Kernel 方法对原来的线性 SVM 进行了推广,使得非线性的的情况也能处理。虽然通过映射将原始数据映射到高维空间之后,能够线性分隔的概率大大增加,但是对于某些情况还是很难处理。

        例如可能并不是因为数据本身是非线性结构的,而只是因为数据有噪音。对于这种偏离正常位置很远的数据点,我们称之为 outlier ,在我们原来的 SVM 模型里,outlier 的存在有可能造成很大的影响,因为超平面本身就是只有少数几个 support vector 组成的,如果这些 support vector 里又存在 outlier 的话,其影响就很大了。例如下图:

        用黑圈圈起来的那个蓝点是一个 outlier ,它偏离了自己原本所应该在的那个半空间,如果直接忽略掉它的话,原来的分隔超平面还是挺好的,但是由于这个 outlier 的出现,导致分隔超平面不得不被挤歪了,变成途中黑色虚线所示(这只是一个示意图,并没有严格计算精确坐标),同时 margin 也相应变小了。当然,更严重的情况是,如果这个 outlier 再往右上移动一些距离的话,我们将无法构造出能将数据分开的超平面来。

        为了处理这种情况,SVM 允许数据点在一定程度上偏离一下超平面。例如上图中,黑色实线所对应的距离,就是该 outlier 偏离的距离,如果把它移动回来,就刚好落在原来的 超平面 蓝色间隔边界上,而不会使得超平面发生变形了。

        插播下一位读者@Copper_PKU的理解:换言之,在有松弛的情况下outline点也属于支持向量SV,同时,对于不同的支持向量,拉格朗日参数的值也不同,如此篇论文《Large Scale Machine Learning》中的下图所示:

        对于远离分类平面的点值为0;对于边缘上的点值在[0, 1/L]之间,其中,L为训练数据集个数,即数据集大小;对于outline数据和内部的数据值为1/L。更多请参看本文文末参考条目第51条。

        OK,继续回到咱们的问题。我们,原来的约束条件为:

        现在考虑到outlier问题,约束条件变成了:

        其中称为松弛变量 (slack variable) ,对应数据点允许偏离的 functional margin 的量。当然,如果我们运行任意大的话,那任意的超平面都是符合条件的了。所以,我们在原来的目标函数后面加上一项,使得这些的总和也要最小:

        其中 是一个参数,用于控制目标函数中两项(“寻找 margin 最大的超平面”和“保证数据点偏差量最小”)之间的权重。注意,其中 是需要优化的变量(之一),而 是一个事先确定好的常量。完整地写出来是这个样子:

        用之前的方法将限制或约束条件加入到目标函数中,得到新的拉格朗日函数,如下所示:

         分析方法和前面一样,转换为另一个问题之后,我们先让针对最小化:

         将  带回 并化简,得到和原来一样的目标函数:

         不过,由于我们得到而又有(作为 Lagrange multiplier 的条件),因此有,所以整个 dual 问题现在写作:

        把前后的结果对比一下(错误修正:图中的Dual formulation中的Minimize应为maxmize):

        可以看到唯一的区别就是现在 dual variable  多了一个上限 。而 Kernel 化的非线性形式也是一样的,只要把换成即可。这样一来,一个完整的,可以处理线性和非线性并能容忍噪音和 outliers 的支持向量机才终于介绍完毕了。

        行文至此,可以做个小结,不准确的说,SVM它本质上即是一个分类方法,用w^T+b定义分类函数,于是求w、b,为寻最大间隔,引出1/2||w||^2,继而引入拉格朗日因子,化为对拉格朗日乘子a的求解(求解过程中会涉及到一系列最优化或凸二次规划等问题),如此,求w.b与求a等价,而a的求解可以用一种快速学习算法SMO,至于核函数,是为处理非线性情况,若直接映射到高维计算恐维度爆炸,故在低维计算,等效高维表现。

        OK,理解到这第二层,已经能满足绝大部分人一窥SVM原理的好奇心,然对于那些想在证明层面理解SVM的则还很不够,但进入第三层理解境界之前,你必须要有比较好的数理基础和逻辑证明能力,不然你会跟我一样,吃不少苦头的。

    第三层、证明SVM

        说实话,凡是涉及到要证明的东西.理论,便一般不是怎么好惹的东西。绝大部分时候,看懂一个东西不难,但证明一个东西则需要点数学功底,进一步,证明一个东西也不是特别难,难的是从零开始发明创造这个东西的时候,则显艰难(因为任何时代,大部分人的研究所得都不过是基于前人的研究成果,前人所做的是开创性工作,而这往往是最艰难最有价值的,他们被称为真正的先驱。牛顿也曾说过,他不过是站在巨人的肩上。你,我则更是如此)。

        正如陈希孺院士在他的著作《数理统计学简史》的第4章、最小二乘法中所讲:在科研上诸多观念的革新和突破是有着很多的不易的,或许某个定理在某个时期由某个人点破了,现在的我们看来一切都是理所当然,但在一切没有发现之前,可能许许多多的顶级学者毕其功于一役,耗尽一生,努力了几十年最终也是无功而返。

        话休絮烦,要证明一个东西先要弄清楚它的根基在哪,即构成它的基础是哪些理论。OK,以下内容基本是上文中未讲到的一些定理的证明,包括其背后的逻辑、来源背景等东西,还是读书笔记。

    本部分导述

    • 3.1节线性学习器中,主要阐述感知机算法;
    • 3.2节非线性学习器中,主要阐述mercer定理;
    • 3.3节、损失函数;
    • 3.4节、最小二乘法;
    • 3.5节、SMO算法;
    • 3.6节、简略谈谈SVM的应用;

    3.1、线性学习器

    3.1.1、感知机算法

        这个感知机算法是1956年提出的,年代久远,依然影响着当今,当然,可以肯定的是,此算法亦非最优,后续会有更详尽阐述。不过,有一点,你必须清楚,这个算法是为了干嘛的:不断的训练试错以期寻找一个合适的超平面(是的,就这么简单)。

        下面,举个例子。如下图所示,凭我们的直觉可以看出,图中的红线是最优超平面,蓝线则是根据感知机算法在不断的训练中,最终,若蓝线能通过不断的训练移动到红线位置上,则代表训练成功。

        既然需要通过不断的训练以让蓝线最终成为最优分类超平面,那么,到底需要训练多少次呢?Novikoff定理告诉我们当间隔是正的时候感知机算法会在有限次数的迭代中收敛,也就是说Novikoff定理证明了感知机算法的收敛性,即能得到一个界,不至于无穷循环下去。

    •  Novikoff定理:如果分类超平面存在, 仅需在序列S上迭代几次,在界为的错误次数下就可以找到分类超平面,算法停止。

        这里为扩充间隔。根据误分次数公式可知, 迭代次数与对应于扩充(包括偏置)权重的训练集的间隔有关。

        顺便再解释下这个所谓的扩充间隔即为样本到分类间隔的距离,即从引出的最大分类间隔。OK,还记得上文第1.3节开头的内容么?如下:

        在给出几何间隔的定义之前,咱们首先来看下,如上图所示,对于一个点x,令其垂直投影到超平面上的对应的为x0,由于w是垂直于超平面的一个向量,为样本x到分类间隔的距离,我们有

        然后后续怎么推导出最大分类间隔请回到本文第一、二部分,此处不重复板书。

        同时有一点得注意:感知机算法虽然可以通过简单迭代对线性可分数据生成正确分类的超平面,但不是最优效果,那怎样才能得到最优效果呢,就是上文中第一部分所讲的寻找最大分类间隔超平面。此外,Novikoff定理的证明请见这里

    3.2、非线性学习器

    3.2.1、Mercer定理

        Mercer定理 :如果函数K是上的映射(也就是从两个n维向量映射到实数域)。那么如果K是一个有效核函数(也称为Mercer核函数),那么当且仅当对于训练样例,其相应的核函数矩阵是对称半正定的。 

        要理解这个Mercer定理,先要了解什么是半正定矩阵,要了解什么是半正定矩阵,先得知道什么是正定矩阵(矩阵理论博大精深,关于矩阵推荐我正在看的一本《矩阵分析与应用》)。然后这里有一个此定理的证明,可以看下。

        正如@Copper_PKU所说:核函数在SVM的分类效果中起了重要的作用,最后这里有个tutorial可以看看。

    3.3、损失函数

        在本文1.0节有这么一句话“支持向量机(SVM)是90年代中期发展起来的基于统计学习理论的一种机器学习方法,通过寻求结构化风险最小来提高学习机泛化能力,实现经验风险和置信范围的最小化,从而达到在统计样本量较少的情况下,亦能获得良好统计规律的目的。”但初次看到的读者可能并不了解什么是结构化风险,什么又是经验风险。要了解这两个所谓的“风险”,还得又从监督学习说起。

        监督学习实际上就是一个经验风险或者结构风险函数的最优化问题。风险函数度量平均意义下模型预测的好坏,模型每一次预测的好坏用损失函数来度量。它从假设空间F中选择模型f作为决策函数,对于给定的输入X,由f(X)给出相应的输出Y,这个输出的预测值f(X)与真实值Y可能一致也可能不一致,用一个损失函数来度量预测错误的程度。损失函数记为L(Y, f(X))。

        常用的损失函数有以下几种(基本引用自《统计学习方法》):

         

        如此,SVM有第二种理解,即最优化+损失最小,或如@夏粉_百度所说“可从损失函数和优化算法角度看SVM,boosting,LR等算法,可能会有不同收获”。

        OK,关于更多统计学习方法的问题,请参看此文

        关于损失函数,如下文读者评论中所述:可以看看张潼的这篇《Statistical behavior and consistency of classification methods based on convex risk minimization》。各种算法中常用的损失函数基本都具有fisher一致性,优化这些损失函数得到的分类器可以看作是后验概率的“代理”。此外,张潼还有另外一篇论文《Statistical analysis of some multi-category large margin classification methods》,在多分类情况下margin loss的分析,这两篇对Boosting和SVM使用的损失函数分析的很透彻。

    3.4、最小二乘法

    3.4.1、什么是最小二乘法?

        既然本节开始之前提到了最小二乘法,那么下面引用《正态分布的前世今生》里的内容稍微简单阐述下。

        我们口头中经常说:一般来说,平均来说。如平均来说,不吸烟的健康优于吸烟者,之所以要加“平均”二字,是因为凡事皆有例外,总存在某个特别的人他吸烟但由于经常锻炼所以他的健康状况可能会优于他身边不吸烟的朋友。而最小二乘法的一个最简单的例子便是算术平均。

        最小二乘法(又称最小平方法)是一种数学优化技术。它通过最小化误差的平方和寻找数据的最佳函数匹配。利用最小二乘法可以简便地求得未知的数据,并使得这些求得的数据与实际数据之间误差的平方和为最小。用函数表示为:

      使误差「所谓误差,当然是观察值与实际真实值的差量」平方和达到最小以寻求估计值的方法,就叫做最小二乘法,用最小二乘法得到的估计,叫做最小二乘估计。当然,取平方和作为目标函数只是众多可取的方法之一。

       最小二乘法的一般形式可表示为:

        有效的最小二乘法是勒让德在 1805 年发表的,基本思想就是认为测量中有误差,所以所有方程的累积误差为

        我们求解出导致累积误差最小的参数即可:

        勒让德在论文中对最小二乘法的优良性做了几点说明:

    •  最小二乘使得误差平方和最小,并在各个方程的误差之间建立了一种平衡,从而防止某一个极端误差取得支配地位
    •  计算中只要求偏导后求解线性方程组,计算过程明确便捷
    • 最小二乘可以导出算术平均值作为估计值

        对于最后一点,从统计学的角度来看是很重要的一个性质。推理如下:假设真值为,x1, ... , xn为n次测量值, 每次测量的误差为ei = xi - ,按最小二乘法,误差累积为

        求解 使达到最小,正好是算术平均

        由于算术平均是一个历经考验的方法,而以上的推理说明,算术平均是最小二乘的一个特例,所以从另一个角度说明了最小二乘方法的优良性,使我们对最小二乘法更加有信心。

        最小二乘法发表之后很快得到了大家的认可接受,并迅速的在数据分析实践中被广泛使用。不过历史上又有人把最小二乘法的发明归功于高斯,这又是怎么一回事呢。高斯在1809年也发表了最小二乘法,并且声称自己已经使用这个方法多年。高斯发明了小行星定位的数学方法,并在数据分析中使用最小二乘方法进行计算,准确的预测了谷神星的位置。

        说了这么多,貌似跟本文的主题SVM没啥关系呀,别急,请让我继续阐述。本质上说,最小二乘法即是一种参数估计方法,说到参数估计,咱们得从一元线性模型说起。

    3.4.2、最小二乘法的解法

        什么是一元线性模型呢? 请允许我引用这里的内容,先来梳理下几个基本概念:

    • 监督学习中,如果预测的变量是离散的,我们称其为分类(如决策树,支持向量机等),如果预测的变量是连续的,我们称其为回归。
    • 回归分析中,如果只包括一个自变量和一个因变量,且二者的关系可用一条直线近似表示,这种回归分析称为一元线性回归分析。
    • 如果回归分析中包括两个或两个以上的自变量,且因变量和自变量之间是线性关系,则称为多元线性回归分析。
    • 对于二维空间线性是一条直线;对于三维空间线性是一个平面,对于多维空间线性是一个超平面...   

        对于一元线性回归模型, 假设从总体中获取了n组观察值(X1,Y1),(X2,Y2), …,(Xn,Yn)。对于平面中的这n个点,可以使用无数条曲线来拟合。要求样本回归函数尽可能好地拟合这组值。综合起来看,这条直线处于样本数据的中心位置最合理。 

        选择最佳拟合曲线的标准可以确定为:使总的拟合误差(即总残差)达到最小。有以下三个标准可以选择:        

    1. 用“残差和最小”确定直线位置是一个途径。但很快发现计算“残差和”存在相互抵消的问题。
    2. 用“残差绝对值和最小”确定直线位置也是一个途径。但绝对值的计算比较麻烦。
    3. 最小二乘法的原则是以“残差平方和最小”确定直线位置。用最小二乘法除了计算比较方便外,得到的估计量还具有优良特性。这种方法对异常值非常敏感。  

         最常用的是普通最小二乘法( Ordinary  Least Square,OLS):所选择的回归模型应该使所有观察值的残差平方和达到最小,即采用平方损失函数。  

         我们定义样本回归模型为:

            

        其中ei为样本(Xi, Yi)的误差。

        接着,定义平方损失函数Q:

        则通过Q最小确定这条直线,即确定为变量,把它们看作是Q的函数,就变成了一个求极值的问题,可以通过求导数得到。

        求Q对两个待估参数的偏导数:

        根据数学知识我们知道,函数的极值点为偏导为0的点。   

        解得:

            

             

        这就是最小二乘法的解法,就是求得平方损失函数的极值点。自此,你看到求解最小二乘法与求解SVM问题何等相似,尤其是定义损失函数,而后通过偏导求得极值。

       OK,更多请参看陈希孺院士的《数理统计学简史》的第4章、最小二乘法。

    3.5、SMO算法

        在上文中,我们提到了求解对偶问题的序列最小最优化SMO算法,但并未提到其具体解法。首先看下最后悬而未决的问题:

        等价于求解:

        1998年,Microsoft Research的John C. Platt在论文《Sequential Minimal Optimization:A Fast Algorithm for Training Support Vector Machines》中提出针对上述问题的解法:SMO算法,它很快便成为最快的二次规划优化算法,特别是在针对线性SVM和数据稀疏时性能更优。

        接下来,咱们便参考John C. Platt的这篇文章来看看SMO的解法是怎样的。

    3.5.1、SMO算法的推导

        咱们首先来定义特征到结果的输出函数:

        注:这个u与我们之前定义的实质是一样的。

        接着,重新定义下咱们原始的优化问题,权当重新回顾,如下:

        求导得到:

        代入中,可得

        通过引入拉格朗日乘子转换为对偶问题后,得:

      s.t:

        注:这里得到的min函数与我们之前的max函数实质也是一样,因为把符号变下,即由min转化为max的问题,且yi也与之前的等价,yj亦如此。

        经过加入松弛变量后,模型修改为:

        从而最终我们的问题变为:

        下面要解决的问题是:在上求上述目标函数的最小值。为了求解这些乘子,每次从中任意抽取两个乘子,然后固定以外的其它乘子,使得目标函数只是关于的函数。这样,不断的从一堆乘子中任意抽取两个求解,不断的迭代求解子问题,最终达到求解原问题的目的。

        而原对偶问题的子问题的目标函数可以表达为:

        其中

        为了解决这个子问题,首要问题便是每次如何选取。实际上,其中一个乘子是违法KKT条件最严重的,另外一个乘子则由另一个约束条件选取。

        根据KKT条件可以得出目标函数中取值的意义:

        这里的还是拉格朗日乘子:

    1. 对于第1种情况,表明是正常分类,在间隔边界内部(我们知道正确分类的点);
    2. 对于第2种情况,表明了是支持向量,在间隔边界上;
    3. 对于第3种情况,表明了是在两条间隔边界之间;

        而最优解需要满足KKT条件,即上述3个条件都得满足,以下几种情况出现将会出现不满足:

    • <=1但是<C则是不满足的,而原本=C
    • >=1但是>0则是不满足的,而原本=0
    • =1但是=0或者=C则表明不满足的,而原本应该是0<<C

        也就是说,如果存在不满足KKT条件的,那么需要更新这些,这是第一个约束条件。此外,更新的同时还要受到第二个约束条件的限制,即

        因此,如果假设选择的两个乘子,它们在更新之前分别是,更新之后分别是,那么更新前后的值需要满足以下等式才能保证和为0的约束:

        其中,是常数。

        两个因子不好同时求解,所以可先求第二个乘子的解(),得到的解()之后,再用的解()表示的解()。

        为了求解,得先确定的取值范围。假设它的上下边界分别为H和L,那么有:

        接下来,综合这两个约束条件,求取的取值范围。

        当y1 != y2时,根据可得,所以有,如下图所示:

        当y1 = y2时,同样根据可得:,所以有,如下图所示:

        如此,根据y1和y2异号或同号,可得出的上下界分别为:

        回顾下第二个约束条件,令上式两边乘以y1,可得

        其中,

        因此可以用表示,,从而把子问题的目标函数转换为只含的问题:

        对求导,可得

        化简下:

        然后将、和代入上式可得:

        令(表示预测值与真实值之差),,然后上式两边同时除以,得到一个关于单变量的解:

        这个解没有考虑其约束条件,即是未经剪辑时的解。

        然后考虑约束可得到经过剪辑后的的解析解为:

        求出了后,便可以求出,得

        那么如何选择乘子呢?

    • 对于,即第一个乘子,可以通过刚刚说的那3种不满足KKT的条件来找;
    • 而对于第二个乘子可以寻找满足条件:的乘子。

        而b在满足下述条件:

        下更新b:

        且每次更新完两个乘子的优化后,都需要再重新计算b,及对应的Ei值。

        最后更新所有,y和b,这样模型就出来了,从而即可求出咱们开头提出的分类函数:

        此外,这里也有一篇类似的文章,大家可以参考下。

    3.5.2、SMO算法的步骤

        综上,总结下SMO的主要步骤,如下:

        意思是,

    1. 第一步选取一对,选取方法使用启发式方法;
    2. 第二步,固定除之外的其他参数,确定W极值条件下的表示。

        假定在某一次迭代中,需要更新对应的拉格朗日乘子,那么这个小规模的二次规划问题写为:

        那么在每次迭代中,如何更新乘子呢?引用这里的两张PPT说明下:

        知道了如何更新乘子,那么选取哪些乘子进行更新呢?具体选择方法有以下两个步骤:

    1. 步骤1:先“扫描”所有乘子,把第一个违反KKT条件的作为更新对象,令为a1;
    2. 步骤2:在所有不违反KKT条件的乘子中,选择使|E1 −E2|最大的a2进行更新,使得能最大限度增大目标函数的值(类似于梯度下降. 此外,而,求出来的E代表函数ui对输入xi的预测值与真实输出类标记yi之差)。

        最后,每次更新完两个乘子的优化后,都需要再重新计算b,及对应的Ei值。

        综上,SMO算法的基本思想是将Vapnik在1982年提出的Chunking方法推到极致,SMO算法每次迭代只选出两个分量ai和aj进行调整,其它分量则保持固定不变,在得到解ai和aj之后,再用ai和aj改进其它分量。与通常的分解算法比较,尽管它可能需要更多的迭代次数,但每次迭代的计算量比较小,所以该算法表现出较好的快速收敛性,且不需要存储核矩阵,也没有矩阵运算。

    3.5.3、SMO算法的实现

        行文至此,我相信,SVM理解到了一定程度后,是的确能在脑海里从头至尾推导出相关公式的,最初分类函数,最大化分类间隔,max1/||w||,min1/2||w||^2,凸二次规划,拉格朗日函数,转化为对偶问题,SMO算法,都为寻找一个最优解,一个最优分类平面。一步步梳理下来,为什么这样那样,太多东西可以追究,最后实现。如下图所示:

        至于下文中将阐述的核函数则为是为了更好的处理非线性可分的情况,而松弛变量则是为了纠正或约束少量“不安分”或脱离集体不好归类的因子。

        台湾的林智仁教授写了一个封装SVM算法的libsvm库,大家可以看看,此外这里还有一份libsvm的注释文档。

        除了在这篇论文《fast training of support vector machines using sequential minimal optimization》中platt给出了SMO算法的逻辑代码之外,这里也有一份SMO的实现代码,大家可以看下。

    3.6、SVM的应用

        或许我们已经听到过,SVM在很多诸如文本分类,图像分类,生物序列分析和生物数据挖掘,手写字符识别等领域有很多的应用,但或许你并没强烈的意识到,SVM可以成功应用的领域远远超出现在已经在开发应用了的领域。

    3.6.1、文本分类

        一个文本分类系统不仅是一个自然语言处理系统,也是一个典型的模式识别系统,系统的输入是需要进行分类处理的文本,系统的输出则是与文本关联的类别。由于篇幅所限,其它更具体内容本文将不再详述。

        OK,本节虽取标题为证明SVM,但聪明的读者们想必早已看出,其实本部分并无多少证明部分(特此致歉),怎么办呢?可以参阅《支持向量机导论》一书,此书精简而有趣。本节完。

    读者评论

       本文发表后,微博上的很多朋友给了不少意见,以下是节选的一些精彩评论:

    1. “压力”陡增的评论→//@藏了个锋:我是看着July大神的博文长大的啊//@zlkysl:就是看了最后那一篇才决定自己的研究方向为SVM的。--Sina Visitor System
    2. @张金辉:SVM的三重境界,不得不转的一篇。其实Coursera的课堂上Andrew Ng讲过支持向量机,但显然他没有把这作为重点,加上Ng讲支持向量机的方法我一时半会难以完全消化,所以听的也是一知半解。真正开始了解支持向量机就是看的这篇“三重境界”,之后才对这个算法有了大概的概念,以至如何去使用,再到其中的原理为何,再到支持向量机的证明等。总之,这篇文章开启了我长达数月的研究支持向量机阶段,直到今日--http://zhan.renren.com/profile/249335584?from=template#!//tag/三重境界
    3. @孤独之守望者:"最后,推出svm的cost function 是hinge loss,然后对比其他的方法的cost function,说明其实他们的目标函数很像,那么问题是svm为什么这么popular呢?您可以再加些VC dimension跟一些error bound的数学,点一下,提供一个思路和方向"。--Sina Visitor System
    4. @夏粉_百度:在面试时,考察SVM可考察机器学习各方面能力:目标函数,优化过程,并行方法,算法收敛性,样本复杂度,适用场景,调参经验,不过个人认为考察boosting和LR也还不错啊。此外,随着统计机器学习不断进步,SVM只被当成使用了一个替代01损失hinge研究,更通用的方法被提出,损失函数研究替代损失与贝叶斯损失关系,算法稳定性研究替代损失与推广性能关系,凸优化研究如何求解凸目标函数,SVM,boosting等算法只是这些通用方法的一个具体组建而已。
    5. @居里猴姐:关于SVM损失函数的问题,可以看看张潼老师的这篇《Statistical behavior and consistency of classification methods based on convex risk minimization》。各种算法中常用的损失函数基本都具有fisher一致性,优化这些损失函数得到的分类器可以看作是后验概率的“代理”。此外,张潼老师还有另外一篇论文《Statistical analysis of some multi-category large margin classification methods》,在多分类情况下margin loss的分析,这两篇对Boosting和SVM使用的损失函数分析的很透彻。
    6. @夏粉_百度:SVM用了hinge损失,hinge损失不可导,不如其它替代损失方便优化并且转换概率麻烦。核函数也不太用,现在是大数据时代,样本非常大,无法想象一个n^2的核矩阵如何存储和计算。 而且,现在现在非线性一般靠深度学习了。//@Copper_PKU:请教svm在工业界的应用典型的有哪些?工业界如何选取核函数,经验的方法?svm的训练过程如何优化?
    7. @Copper_PKU:July的svm tutorial 我个人觉得还可以加入和修改如下部分:(1) 对于支持向量解释,可以结合图和拉格朗日参数来表达,松弛中sv没有写出来. (2) SMO算法部分,加入Joachims论文中提到的算法,以及SMO算法选取workset的方法,包括SMO算法的收敛判断,还有之前共轭梯度求解方法,虽然是较早的算法,但是对于理解SMO算法有很好的效果。模型的优化和求解都是迭代的过程,加入历史算法增强立体感。--  Sina Visitor System
    8. //@廖临川: 之所以sgd对大训练集的效果更好,1.因为SGD优化每次迭代使用样本子集,比使用训练全集(尤其是百万数量级)要快得多;2.如果目标函数是凸的或者伪凸的,SGD几乎必然可以收敛到全局最优;否则,则收敛到局部最优;3.SGD一般不需要收敛到全局最优,只要得到足够好的解,就可以立即停止。//@Copper_PKU:sgd的核心思想:是迭代训练,每拿到一个样本就算出基于当前w(t) 的loss function,t代表训练第t次,然后进行下一w(t+1)的更新,w(t+1)=w(t)-(learning rate) * loss function的梯度,这个类比神经网络中bp中的参数训练方法。 sample by sample就是每次仅处理一个样本 而不是一个batch。
    9. //@Copper_PKU:从损失函数角度说:primal问题可以理解为正则化项+lossfunction,求解目标是在两个中间取平衡 如果强调loss function最小则会overfitting,所以有C参数。 //@研究者July:SVM还真就是在一定限定条件下,即约束条件下求目标函数的最优值问题,同时,为减少误判率,尽量让损失最小。
    10. ...

        非常享受这种全民大讨论的年代,没有谁一定就对或一定就错,而是各自发表各自的理解见解,真棒!

    参考文献及推荐阅读

    1. 支持向量机导论》,[英] Nello Cristianini / John Shawe-Taylor 著
    2. 支持向量机导论一书的支持网站:http://www.support-vector.net/
    3. 《数据挖掘导论》,[美] Pang-Ning Tan / Michael Steinbach / Vipin Kumar 著
    4. 《数据挖掘:概念与技术》,(加)Jiawei Han;Micheline Kamber 著;
    5. 《数据挖掘中的新方法:支持向量机》,邓乃扬 田英杰 著;
    6. 支持向量机--理论、算法和扩展》,邓乃扬 田英杰 著;
    7. 支持向量机系列pluskid支持向量机系列 « Free Mind
    8. 支持向量机及其应用
    9. 数据挖掘十大经典算法初探
    10. 《模式识别支持向量机指南》,C.J.C Burges 著;
    11. 统计学习方法》,李航著;
    12. 《统计自然语言处理》,宗成庆编著,第十二章、文本分类;
    13. SVM入门系列,Jasper:文本分类技术 - Jasper's Java Jacal - BlogJava
    14. 最近邻决策和SVM数字识别的实现和比较,作者不详;
    15. 纯白板手推SVM:视频播放
    16. 斯坦福大学机器学习课程原始讲义:PDF版学习笔记 - JerryLead - 博客园
    17. 斯坦福机器学习课程笔记Machine Learning - 标签 - JerryLead - 博客园
    18. 支持向量机SVM(一) - JerryLead - 博客园
    19. SMO算法的数学推导:支持向量机(五)SMO算法 - JerryLead - 博客园
    20. 数据挖掘掘中所需的概率论与数理统计知识、上
    21. 关于机器学习方面的文章,可以读读:机器学习 - 随笔分类 - Leo Zhang - 博客园
    22. 数学系教材推荐:数学系教材推荐_Flatlander_新浪博客
    23. 《神经网络与机器学习(原书第三版)》,[加] Simon Haykin 著;
    24. 正态分布的前世今生:http://t.cn/zlH3Ygc
    25. 数理统计学简史》,陈希孺院士著;
    26. 《最优化理论与算法(第2版)》,陈宝林编著;
    27. A Gentle Introduction to Support Vector Machines in Biomedicinehttp://www.nyuinformatics.org/downloads/supplements/SVM_Tutorial_2010/Final_WB.pdf,此PPT很赞,除了对引入拉格朗日对偶变量后的凸二次规划问题的深入度不够之外,其它都挺好,配图很精彩,本文有几张图便引自此PPT中;
    28. 来自卡内基梅隆大学carnegie mellon university(CMU)的讲解SVM的PPT:http://www.autonlab.org/tutorials/svm15.pdf
    29. 发明libsvm的台湾林智仁教授06年的机器学习讲义SVM:http://wenku.baidu.com/link?url=PWTGMYNb4HGUrUQUZwTH2B4r8pIMgLMiWIK1ymVORrds_11VOkHwp-JWab7IALDiors64JW_6mD93dtuWHwFWxsAk6p0rzchR8Qh5_4jWHC
    30. http://staff.ustc.edu.cn/~ketang/PPT/PRLec5.pdf
    31. Introduction to Support Vector Machines (SVM),By Debprakash Patnai M.E (SSA),https://www.google.com.hk/url?sa=t&rct=j&q=&esrc=s&source=web&cd=1&ved=0CCwQFjAA&url=http%3a%2f%2fwww%2epws%2estu%2eedu%2etw%2fccfang%2findex%2efiles%2fAI%2fAI%26ML-Support%2520Vector%2520Machine-1%2eppt&ei=JRR6UqT5C-iyiQfWyIDgCg&usg=AFQjCNGw1fTbpH4ltQjjmx1d25ZqbCN9nA
    32. 多人推荐过的libsvmLIBSVM -- A Library for Support Vector Machines
    33. 《machine learning in action》,中文版为《机器学习实战》;
    34. SMO算法的提出:Sequential Minimal Optimization A Fast Algorithm for Training Support Vector MachinesMicrosoft researchers and engineers working around the world
    35. 《统计学习理论的本质》,[美] Vladimir N. Vapnik著,非常晦涩,不做过多推荐;
    36. 张兆翔,机器学习第五讲之支持向量机http://irip.buaa.edu.cn/~zxzhang/courses/MachineLearning/5.pdf
    37. VC维的理论解释:http://www.svms.org/vc-dimension/,中文VC维解释http://xiaoxia001.iteye.com/blog/1163338
    38. 来自NEC Labs America的Jason Weston关于SVM的讲义http://www.cs.columbia.edu/~kathy/cs4701/documents/jason_svm_tutorial.pdf
    39. 来自MIT的SVM讲义:http://www.mit.edu/~9.520/spring11/slides/class06-svm.pdf
    40. PAC问题:http://www.cs.huji.ac.il/~shashua/papers/class11-PAC2.pdf
    41. 百度张潼老师的两篇论文:《Statistical behavior and consistency of classification methods based on convex risk minimization》http://home.olemiss.edu/~xdang/676/Consistency_of_Classification_Convex_Risk_Minimization.pdf,《Statistical analysis of some multi-category large margin classification methods》;
    42. 没有找到站点
    43. 《矩阵分析与应用》,清华张贤达著;
    44. SMO算法的实现:SVM算法实现(一)_techq's blog-CSDN博客_svm算法
    45. 常见面试之机器学习算法思想简单梳理:机器学习&数据挖掘笔记_16(常见面试之机器学习算法思想简单梳理) - tornadomeet - 博客园
    46. 矩阵的wikipedia页面:http://zh.wikipedia.org/wiki/%E7%9F%A9%E9%98%B5
    47. 最小二乘法及其实现:一元线性回归模型与最小二乘法及其C++实现_Liam Q的专栏-CSDN博客
    48. 统计学习方法概论:统计学习笔记(1)——统计学习方法概论_Liam Q的专栏-CSDN博客_统计学习方法学习笔记
    49. 程序员训练机器学习 SVM算法分享-CSDN.NET
    50. A Tutorial on Support Vector Regression:http://alex.smola.org/papers/2003/SmoSch03b.pdf;SVR简明版:http://www.cmlab.csie.ntu.edu.tw/~cyy/learning/tutorials/SVR.pdf
    51. SVM Org:Account Suspended
    52. R. Collobert. Large Scale Machine Learning. Université Paris VI phd thesis. 2004:http://ronan.collobert.com/pub/matos/2004_phdthesis_lip6.pdf
    53. Making Large-Scale SVM Learning Practical:http://www.cs.cornell.edu/people/tj/publications/joachims_99a.pdf
    54. 文本分类与SVM:文本分类与SVM_张知临的专栏-CSDN博客_文本分类svm
    55. Working Set Selection Using Second Order Information
      for Training Support Vector Machines:http://www.csie.ntu.edu.tw/~cjlin/papers/quadworkset.pdf
    56. SVM Optimization: Inverse Dependence on Training Set Size:http://icml2008.cs.helsinki.fi/papers/266.pdf
    57. Large-Scale Support Vector Machines: Algorithms and Theory:http://cseweb.ucsd.edu/~akmenon/ResearchExam.pdf
    58. 凸优化的概念:http://cs229.stanford.edu/section/cs229-cvxopt.pdf
    59. 《凸优化》,作者: Stephen Boyd / Lieven Vandenberghe,原作名: Convex Optimization;
    60. Large-scale Non-linear Classification: Algorithms and Evaluations,Zhuang Wang,讲了很多SVM算法的新进展:http://ijcai13.org/files/tutorial_slides/te2.pdf
    61. 基于SMO算法实现SVM:http://www.cs.iastate.edu/~honavar/smo-svm.pdf
    62. copper推荐的一些SVM相关的论文(当然,其中不少论文在上面的条目中都已经提到):http://c.blog.sina.com.cn/profile.php?blogid=68d0b92d89000h35
    63. 在线编辑Latex 公式:http://private.codecogs.com/latex/eqneditor.php?lang=zh-cn

    后记

        OK,此文从最初2012年5月开始动笔,到后续不断的修改,创造了三个之最,即所写时间最长,所花心血最大,所改次数最多,因为我的目标是让没有任何机器学习基础的都能看懂此文,所以总是不停的改,不停的改,不想放过任何一个小的细节。再者,引用侯捷的一句话是:天下大作,必作于细。

        最后,非常感谢pluskid及诸多朋友们的文章及著作,让我有机会在其基础上总结、深入。有任何问题,敬请广大读者随时不吝批评指正,感谢。


    本文PDF版

        本文会一直不断翻新,再者,上述 4 个PDF的阅读体验也还不是最好的,如果有朋友制作了更好的PDF,欢迎分享给我:Sina Visitor System,谢谢。

        July、二零一六年一月十七日第N次修改(N > 100)。

    展开全文
  • 它是一种二类分类模型,其基本模型定义为特征空间上的间隔最大的线性分类器,即支持向量机的学习策略便是间隔最大化,最终可转化为一个凸二次规划问题的求解。 2、求解过程 1、数据分类—SVM引入 假设在一个二维...
  • 优化理论是针对最优化问题而言的,通常情况下,我们的最优化问题涉及到以下三种情况:1)无约束优化问题形如:min f(x);对于此类问题,我们采用fermat定理,即对函数求导得到导数为0的点,该点也就是极值。将候选点...
  • 非负矩阵分解

    2015-06-04 11:15:50
    由于分解前后的矩阵中仅包含非负的元素,因此,原矩阵A中的一列向量可以解释为对左矩阵U中所有列向量(称为基向量)的加权和,而权重系数为右矩阵V中对应列向量中的元素。这种基于基向量组合的表示形式具有很直观的...
  • 非负矩阵分解(NMF)

    千次阅读 2017-11-28 20:31:07
    NMF的基本思想可以简单描述为:对于任意给定的一个非负矩阵A,NMF算法能够寻找到一个非负矩阵U和一个非负矩阵V,使得满足 ,从而将一个非负的矩阵分解为左右两个非负矩阵的乘积。
  • 上一次介绍支持向量机,结果说到 Maximum Margin Classifier ,到最后都没有说“支持向量”到底是什么东西。不妨回忆一下上次最后一张图: 可以看到两个支撑着中间的 gap 的超平面,它们到中间的 separating ...
  • paper4中介绍了支持向量机,结果说到 Maximum Margin Classifier ,到最后都没有说“支持向量”到底是什么东西。不妨回忆一下上次最后一张图: 可以看到两个支撑着中间的 gap 的超平面,它们到中间的 separating...
  • 支持向量机SVM分析

    万次阅读 2013-07-06 18:27:08
    1995年Vapnik等人[2]提出一种机器学习的新方法支持向量机(Support Vector Machine,SVM)之后,支持向量机成为继人工神经网络之后又一研究热点。SVM 算法的核心在于最优分类超平面的确定,即通过训练样本确定分类器...
  • 基于支持向量机的图像分类(上篇)

    万次阅读 多人点赞 2018-03-31 21:25:39
    摘要:本文通过图文详细介绍如何利用支持向量机对图像进行分类。这篇文章从什么是图像分类任务开始一步步详细介绍支持向量机原理,以及如何用它解决图像多分类任务。将这部分内容分为上下两篇:上篇重点详细介绍实现...
  • 优化

    2018-08-21 11:47:39
    后者连局部极值都不是,在鞍点处Hessian矩阵不定,即既非正定,也非负定。   凸优化通过对目标函数,优化变量的可行域进行限定,可以保证不会遇到上面两个问题。   凸优化是一类特殊的优化问题,它要求: ...
  • NMF非负矩阵分解

    2018-01-26 10:33:00
    该文提出了一种新的矩阵分解思想――非负矩阵分解(Non-negative Matrix Factorization,NMF)算法,即NMF是在矩阵中所有元素均为非负数约束条件之下的矩阵分解方法。该论文的发表迅速引起了各个领域中的科学研究人员...
  • 非负矩阵分解 即将一个非负的大矩阵(V)分解成两个小的非负矩阵(W和H)。 WH≈V WH≈V WH≈V Vn*m ,Wn*k ,Hk*m 矩阵 ,其中(n+m)r &lt;&lt;nm。 可以参考: Learning the parts of objects by non-negative...
  • 非负矩阵分解(NMF) NMF的基本思想 为什么分解的矩阵式非负? 为什么要运用非负矩阵分解? NMF的基本思想: 对于任意给定的一个非负矩阵V,NMF算法能够通过计算,从原矩阵提取权重矩阵和特征矩阵两个不同的...
  • NMF(非负矩阵分解)算法

    千次阅读 2016-01-04 11:12:24
    由于分解前后的矩阵中仅包含非负的元素,因此,原矩阵A中的一列向量可以解释为对左矩阵U中所有列向量(称为基向量)的加权和,而权重系数为右矩阵V中对应列向量中的元素。这种基于基向量组合的表示形式具有很直观的...
  • NMF 非负矩阵分解

    2014-07-21 16:46:39
    该文提出了一种新的矩阵分解思想――非负矩阵分解(Non-negative Matrix Factorization,NMF)算法,即NMF是在矩阵中所有元素均为非负数约束条件之下的矩阵分解方法。该论文的发表迅速引起了各个领域中的科学研究人员...
  • 支持向量

    2017-10-17 16:54:25
    支持向量机通俗导论(理解SVM的三层境界) 作者:July、pluskid ;致谢:白石、JerryLead 出处:结构之法算法之道blog。 前言  动笔写这个支持向量机(support vector machine)是费了不少劲...
  • C++ 向量 vector

    2020-03-07 09:59:22
    并且每一个元素均可由(非负)编号唯一指代,并且**可以直接访问**,A[i]的物理地址是A+i*s(s是单个元素暂用的空间量),也称作线性数组。 向量由数组抽象来,由一组元素按照线性次序封装,各个元素[0,n)内的秩(int ,c++...
  • 深入理解支持向量

    2020-06-01 22:06:01
    文章目录从感知机说起三大核心支持向量KKT条件核方法总结Preference 从感知机说起 三大核心 支持向量 KKT条件 核方法 总结 Preference
  • 为了寻找到一组最优的分解 向量,可以将NMF优化问题转变为如下非凸问题的求解公式1.1: 其中, 代表着矩阵 和 中的元素都是非负的。NMF已经广泛应用在矩阵填充、文本聚类和语音信号处理等领域。 Goodfellow在其标志...

空空如也

空空如也

1 2 3 4 5 ... 20
收藏数 7,212
精华内容 2,884
关键字:

向量非负优化