机器学习——PLA

对研究生要跟的导师还不确定,暑假打算学习点专业课与英语,也迟迟没有做好。我想不管哪个实验室应该都不会离开机器学习吧。因此开始看这方面的东西。

今天接触了一个算法叫PLA(Percetron Learning Algorithm),用来做线性分类的算法。它应用的前提是样本集是线性可分的。用二维特征值的样本来举例子:

linear separable

我们要的就是通过算法来找到这条线。对于二维的样本(特征分别为x1x_1,x2x_2),则分类结果为-1,+1,分别表示为negative,positive(这与另一种常用的分类算法不一致)。思想就是每个特征对应一个参数,乘积之后如果大于某个阈值,则设定为分类结果为+1,小于则为-1。若为特征添加一个特征$ x_0 $恒等于一,则用向量化可以将分类过程写成如下形式: result(X)=sign(WXT) result(X) = sign(W{X^T}) 其中W为特征的参数向量,我们使用yy来表示样本真实类别。

整个算法的思想其实很简单,刚开始画出一条线,如果分错了,则往正确的方向旋转这条线。但是如何旋转这条线,旋转多少角度,还是很有意思的。一般的想法都是从图直观上来看,利用代价函数(cost function)来解决,不过对于这个简单的线性分类,用代价函数进行梯度下降还是相对来说计算量还是比较大的。这个算法让我觉得厉害的地方在于它向量化不光是为了提高计算效率,而是从向量的角度来考虑一步步向结果逼近的:遇到的分类错误有两种情况,如果yy应该是正,但是WXTW{X^T}得到的是负的,从向量角度来说,WWXX的夹角太大,因此内积为负,我们要减小两个向量的角度,可以令W=W+XW=W+X,相反,yy应该是负,结果却为正,那么WWXX的夹角太小,可以令W=WXW=W-X,这样就增大了夹角,从下图可以很直观的看出来这个道理:

向量加减

合并两种情况,则每次遇到被错误分类的样本XnX_n,我们对参数向量WtW_t进行如下更新(其中tt表示修正的次数):Wt+1=Wt+ynXn W_ {t+1} = W_t + y_nX_n 下面主要来证明一下,对于线性可分的样本来说,这个算法一定会停止,找到那条符合的参数向量WfW_f。我们可以知道:

ynWfTXn>=min[ynXnWfT]>0 y_n{W_f}^TX_n >= _ {min}[ y_nX_nW_f^T]>0

Wt+1WfT=WtWft+ynXnWfT>=WtWfT+min[ynXnWfT] W_ {t+1}W_f^T = W_tW_f^t + y_nX_nW_f^T >= W_tW_f^T + _ {min}[ y_nX_nW_f^T]

因此我们可以得到:

Wt+1WfT>=WtWfT>=W0WfT+(t+1)*min[ynXnWfT] W_ {t+1}W_f^T >= W_tW_f^T >= W_0W_f^T +(t+1) *_ {min}[ y_nX_nW_f^T]

因为我们每次找到一个出错点才进行修正,出错点为(X_n,y_nX\_n,y\_n),可以知道$ y_nW_tX_n^T<0 $,则:

Wt+12=Wt2+ynXn2+2ynWtXnT>=Wt2+max[Xn2] \Vert W_ {t+1} \Vert ^2 = \Vert W_t \Vert ^2 + \Vert y_nX_n \Vert ^2+2y_nW_tX_n^T>= \Vert W_t \Vert ^2+_ {max}[ \Vert X_n \Vert ^2] 因此我们可以得到:

Wt+12=Wt2+ynXn2+2×ynWtXnT>=W02+(t+1)×max[Xn2] \Vert W_ {t+1} \Vert ^2 = \Vert W_t \Vert ^2 + \Vert y_nX_n \Vert ^2+2 \times y_nW_tX_n^T>= \Vert W_0 \Vert ^2+ (t+1) \times _ {max}[ \Vert X_n \Vert ^2]

最后,假设最开始W0=0W_0=0,通过正规化,我们可以算出WtW_tWfW_f之间的角度是不断逼近的(角度重合时候也就得到了正确的参数,实际上角度相同后我们不在乎向量的模的大小,因为它不会影响WtTXnW_t^TX_n的符号):

WtWfTWtWf>=t×min[ynWfXn]t×maxXnWf=tqR \frac{W_tW_f^T}{ \Vert W_t \Vert \Vert W_f \Vert } >= \frac{t \times _ {min}[y_nW_fX_n]}{\sqrt{t} \times _ {max} \Vert X_n \Vert \Vert W_f \Vert } = \sqrt{t} \frac q R

其中:q=min[ynWfXn],R2=(maxXn)2q = _ {min}[y_nW_fX_n],R^2 = ( _ {max} \Vert X_n \Vert )^2,而正规化后内积是小于等于1的(此时方向一致),可以得到$ t<= {q^2} ,可以得到,最多经过 {q^2} $次修正即可得到正确结果。

到这里就证明结束了,以后会上传相关的代码。

可以看到虽然算法实现很简单,但其中的推理还是不容易的。