Clustering 聚类
谱聚类
上文我们引入了是聚类,并介绍了第一种聚类算法K-means。今天,我们来介绍一种流行的聚类算法——谱聚类(Spectral Clustering),它的实现简单,而且效果往往好于传统的聚类算法,如k-means,但是其背后的原理涉及了很多重要而复杂的知识,如图论,矩阵分析等。别担心,今天小编就带你一举攻克这些难关,拿下谱聚类算法。
Q:什么是谱聚类?
A: 谱聚类是最流行的聚类算法之一,它的实现简单,而且效果往往胜过传统的聚类算法,如K-means。它的主要思想是把所有数据看作空间中的点,这些点之间用带权重的边相连,距离较远的点之间的边权重较低,距离较近的点之间边权重较高,通过对所有数据点和边组成的图进行切图,让切图后不同子图间边权重和尽可能低,而子图内边权重和尽可能高来达到聚类的目的。
接下来我们先介绍图的基础知识以及图拉普拉斯,然后引入谱聚类算法,最后从图分割的角度来解释谱聚类算法。
图的概念
对于一个图 G ,我们通常用 G(V,E) 来表示它,其中V代表数据集合中的点{v1,v2,...vn},E代表边集(可以有边相连,也可以没有)。在接下来的内容中我们用到的都是带权重图(即两个顶点vi和vj之间的连边带有非负权重wij>0).由此可以得到一个图的 加权邻接矩阵W =(Wij)i,j=1,...n。如果Wij=0代表两个顶点之间无边相连。无向图G中wij=wji,所以权重矩阵是对称的。
对于图中的任意一个点vi,定义它的度为与它相连的所有边的权重之和,即
![3492d3a8d20e080ea737bfcac5e8d185.png](https://i-blog.csdnimg.cn/blog_migrate/2c45de87953a5d185497f569a85838e2.jpeg)
。 度矩阵 就是对角元素分别为每个点的度,非对角元素为0的矩阵。
给定点的一个子集A属于V,定义A的补集为
![47005124d58a4545a10965ee02096e0f.png](https://i-blog.csdnimg.cn/blog_migrate/38d793cbbda50cc2680ac40dac916bb2.jpeg)
A的指示向量为
![34190ce9c800d6e6a7992a373949de1d.png](https://i-blog.csdnimg.cn/blog_migrate/b3edbc9ff735407cca109a1b50e5b18d.jpeg)
如果fi=1,则顶点vi在子集A中,否则为0。为了方便,在下文中,我们使用
![884035adf2968f270d019f104058fceb.png](https://i-blog.csdnimg.cn/blog_migrate/9efbe9c2dd34fdae3603069b78398b45.jpeg)
来表示顶点i在集合A中。对于两个不相交子集A,B属于V,我们定义
![dfd9ba821783aa6e394dd347e9d3f00b.png](https://i-blog.csdnimg.cn/blog_migrate/eb5230efba4a8b62643c0fb3622e3b52.jpeg)
表示两个子集之间的权重和。
有两种度量V中子集A大小的方法:|A|表示A中顶点的个数,vol(A)表示A中顶点的度的和。
相似图
思考我们构建相似性图的目的是什么?是为了对点之间的局部邻域关系建模。那么根据我们所关注的邻域关系,相似性图的定义也可以不同,以下提到的图都经常在谱聚类中使用:
ε邻域图 :将所有距离小于ε的点相连,由于有边相连的点之间都差不多,加权边不会包含关于图中数据点更多的关系,因此,这种图常常是无权重的;
k近邻图 :将vi与它前k近的顶点相连,要注意,这种定义的图是有向图,因为k近邻关系并不是对称的(A是B的k近邻,而B不一定是A的k近邻)。我们可以通过两种方式将它变成无向图:一种是只要A是B的k近邻,AB之间就会连一条边;另一种是必须同时满足A,B是彼此的k近邻,才能在AB之间连一条边。在构建好图之后,在根据点的相似性给边赋权重;
全连接图 :任意两点之间都连一条边,并根据两点之间的相似性给边赋权重。由于图必须能代表局部邻居关系,所以使用的相似性度量方法必须能对这种关系建模。比如,高斯相似性方程,其中参数σ控制邻域的宽度(作用类似ε邻域图中的ε参数)。
图拉普拉斯矩阵及其基本性质
图拉普拉斯矩阵 是谱聚类的主要工具,对这些矩阵的研究有一个专门的领域,叫做谱图理论,感兴趣的同学可以去了解下。在这章中我们定义不同的图拉普拉斯并介绍它们最重要的性质。
在接下来的介绍中,我们定义的图 G 是无向的,带权重的(权重矩阵为 W ),我们假定所有特征向量都是标准化的(如常向量C和aC是一样的),特征值是有序且可重复的,如果提到了前k个特征向量,就意味着这k个特征向量对应k个最小的特征值。D代表上文提到的度矩阵,W代表权重矩阵。
非规范化图拉普拉斯矩阵:
非规范化图拉普拉斯矩阵定义为 L:D-W ,它具有如下重要特性:
1. 对于任意特征向量
![8ac651d1eb98070566da03118b34e278.png](https://i-blog.csdnimg.cn/blog_migrate/f5b2ffd2d7c4f3082bdff38a7e477678.jpeg)
,有:
![0f2c6156d4e755b2ce512053130bbf8f.png](https://i-blog.csdnimg.cn/blog_migrate/13ab190e56cf4a0150e6c46bf50786d6.jpeg)
;
证明:
![08aea49463d8c0657d24fd405d3418f7.png](https://i-blog.csdnimg.cn/blog_migrate/e2092526cfb6ded1b296a9d3e9336c56.jpeg)
2. L是对称半正定的;
证明:D,W都是对称的,且特性1的证明表明L是半正定的;
3. L的最小特征值是0,对应的特征向量是单位向量;
4. L有n个非负的实值特征向量 :
![dc8259852515f1820d7bde85923631c2.png](https://i-blog.csdnimg.cn/blog_migrate/3e5640c1fc43b6225fd1a3862670c3ef.jpeg)
证明:由前三个特性可以直接得出;
注意非规范化图拉普拉斯矩阵以及它的特征向量,特征值可以被用在描述图的很多特性,其中在谱聚类中的一个重要性质如下:
连通分量数与L的普 :G是一个非负权重无向图,那么L的特征值0的重数k就等于图中联通分量A1,...Ak的个数,特征值0的特征空间由这些联通分量的指示向量
![a4dea6ce51a2fe0e3c549ef8649c2e42.png](https://i-blog.csdnimg.cn/blog_migrate/6fd8e638a028b476b339cafadc091598.jpeg)
表示;
证明:以k=1为例,代表这是一个连通图,假设特征值0对应的特征向量为f,则
![095ae66c446749e21444b165054c2702.png](https://i-blog.csdnimg.cn/blog_migrate/c001ebcf30105fb58cd3a32adb3dac9c.jpeg)
,由于权重w是非负的,要使和为0当且仅当每一项都为0。因此,如果两个顶点vi和vj有边相连(wij>0),fi必须等于fj,从这里可以看出如果图中的顶点可以有一条路径相连,那么f必须是常向量。此外,由于无向图中连通分量的所有顶点都可以通过一条路径相连,所以f在整个连通分量上必须是常数。因此在只有一个连通分量的图中,只有一个常向量作为0的特征向量(作为唯一一个联通部分的指示向量)。
现在考虑k大于1的场景,不失一般性,我们假设顶点是按照它们所属的连通分量排序的,在这种情形下,邻接权重矩阵W可以写成分块对角矩阵的形式,同样,L可以写作
![b68b7a03049ad595e2645632201ddcf7.png](https://i-blog.csdnimg.cn/blog_migrate/0580aa7edc6408e089bfab2ae0159af7.jpeg)
注意,每一个块Li 也是一个图拉普拉斯矩阵,分别对应图的第i个连通分量,L的普就由Li的普组合而成,而且L对应的特征向量就是Li的特征向量。由于Li一个连通图的图拉普拉斯矩阵,而每一个Li有重数为1的特征值0,也就是对应第i个联通子图的常向量。
规范化图拉普拉斯矩阵
规范化图拉普拉斯矩阵有如下两种定义:
![d084eeeb367de8922d84fd9c1a06dc94.png](https://i-blog.csdnimg.cn/blog_migrate/f9fda62f47d0f58a13008e58e6338c45.jpeg)
我们先来介绍它们的性质:
1. 对任意特征向量
![04ac3801d47c50608aa64ba1b52fbcf9.png](https://i-blog.csdnimg.cn/blog_migrate/ed00cf91166d39ea2916cfc864c6f450.jpeg)
,我们有:
![e682b730fa913eacc128037e48ff1144.png](https://i-blog.csdnimg.cn/blog_migrate/0c3c19ea99c1cf79d0cb05663a27afdd.jpeg)
;
2. λ是Lrw的特征向量u的特征值 当且仅当 λ是Lsym的特征向量
![81dd6e90f0550578ba11ae0021cea15e.png](https://i-blog.csdnimg.cn/blog_migrate/9a557c0c78ff60987fb2fe6b905f6c0e.jpeg)
的特征值;
3. λ是Lrw的特征向量u的特征值当且仅当λ和u满足Lu=λDu;
4. 0是Lrw的常数特征向量
![6f997724b2b30db84494983743f68ae8.png](https://i-blog.csdnimg.cn/blog_migrate/7fee260892a2794f555c84b8e7e1397d.jpeg)
的特征值,那么0也是Lsym的特征向量
![b456ba4554c4764275cba90aba591870.png](https://i-blog.csdnimg.cn/blog_migrate/b260647c8708c9320bbf42c23d61ad6c.jpeg)
的特征值。;
5. Lsym和Lrw是半正定的且有n个非负的实特征值;
与非规范化图拉普拉斯矩阵相同,Lsym和Lrw的0特征值的重数k就等于图中联通分量的个数,证明与上面类似,不再赘述。
谱聚类算法
现在我们来介绍最常见的谱聚类算法。
![abe6c8f4753e2099236fd6b220f9d6db.png](https://i-blog.csdnimg.cn/blog_migrate/02356c922d51c8a25bb0a93ec7007ef8.jpeg)
解释:首先构建相似性图,并用W表示权重邻接矩阵,构建度矩阵D,求出拉普拉斯矩阵L;计算L的前m个特征向量,以这m个特征向量作为列组成n*m的矩阵U,其中每一行作为一个m维的样本,共n个,对这n个行向量进行kmeans聚类。
对于规范化的谱聚类有两种不同的版本,依赖于两种标准化图拉普拉斯矩阵。
第一种是:
![7fc53fddb9b3785cb4bb57a14da09bf3.png](https://i-blog.csdnimg.cn/blog_migrate/0df656a145e6ffbaa12d18d1b5eb9322.jpeg)
注意这里使用的是generalized eigenvectors(广义特征向量),即矩阵Lrw对应的特征向量,所以这一算法针对的是标准化的拉普拉斯矩阵Lrw。
下一个算法也是标准化的谱聚类,不过用的是Lsym,该算法介绍了一种额外的行标准化步骤:
![67008c11ee1a309efcb14340f2eb5cdb.png](https://i-blog.csdnimg.cn/blog_migrate/516a936a69ab54446217675443f45a06.jpeg)
上面三种算法的步骤其实都差不多,只不过是使用的拉普拉斯矩阵有差别。在算法中,主要的工作就是将x抽象表征为k维的数据点(降维),接下来我们将揭晓为什么这样做能提高聚类的效果。我们将从图分割的角度来介绍谱聚类的工作原理。
图切分的观点
聚类算法的思想就是根据数据点之间的 相似性 将它们 划分 到不同组。当我们的数据以 相似性图 的形式给出时,聚类问题又可以这样解释:我们想要找到图的一种划分,不同组点的边之间有很低的权重而同一个组中点之间的边有较高的权重。在本节中,我们将介绍如何将谱聚类问题近似于图的划分问题。
给定一个相似性图,邻接权重矩阵W,最简单直接的方法来创建一个图的分割就是解决最小割问题。前面我们已经提到过
![14440b842560d09566753b7585734fa6.png](https://i-blog.csdnimg.cn/blog_migrate/380a9c13ba9da0594bd49091abf34c02.jpeg)
和A的补集的概念。给定子集个数k,最小割方法就是选择一种划分A1,...,Ak使得
![08c1a4afd21179911a1bb21079097401.png](https://i-blog.csdnimg.cn/blog_migrate/fa47cb4ea9100285d241c9e6dda5372c.jpeg)
最小。特别是对k=2,最小割问题相对是一个简单的问题,可以被有效地解决。然而在实际中,最小割的解并不能很好地将图划分,因为它只是简单地将一个独立的点从剩余的图中划分出来,这显然不符合聚类的需求。解决这种问题的一种方法是必须保证子集A1,...,Ak足够大,实现这一限制的最常见的目标函数有两种:RatioCut和Ncut。RatioCut用子集A中点的个数表示A的大小,Ncut用子集A中边的权重和来度量A的大小,它们的定义分别如下:
![ed1caa255ad273b2485c84b6a38e9bda.png](https://i-blog.csdnimg.cn/blog_migrate/bfdada7dfe816d1782d59ee6ac94461f.jpeg)
假如簇Ai不是特别小,两个目标函数的值都会有较小的值,因此这两种函数都试图达到聚类的平衡,但是,引入这一“平衡”条件使得之前简单的最小割问题变成了np难问题。谱聚类就是解决这一问题的一种松弛版本,我们将看到松弛Ncut将会导致规范化谱聚类,而松弛RatioCut将导致非标准化谱聚类。
从RatioCut切图的角度解释谱聚类算法
我们先讨论k=2的情形,我们的目标函数是
![a2a40cc36e9e6f27b05e9112f0d0a1e1.png](https://i-blog.csdnimg.cn/blog_migrate/bb53460a01e7fea08a920c51bb0d0064.jpeg)
,首先重写这个问题。给定子集 A ,我们定义向量
![55af2fa4542bb69ccd070d7e00468984.png](https://i-blog.csdnimg.cn/blog_migrate/94fdce01728ba775cc700cf0940daa7b.jpeg)
,其中
![c947d16f4be38fb91900ac77659043da.png](https://i-blog.csdnimg.cn/blog_migrate/5bd57f4bd681880ba0941a81a379cde7.jpeg)
现在RatioCut目标函数可以写作非规范化图拉普拉斯矩阵的形式:
![aafefcb1b909958381da5c52f31fcf8b.png](https://i-blog.csdnimg.cn/blog_migrate/fd44a9d51ce39ce83debbfca7e8c6e90.jpeg)
而,
![19a469392050d5df21d4dd2e50ff881b.png](https://i-blog.csdnimg.cn/blog_migrate/a672cc244b26c9e0a9d6fe33d5200a47.jpeg)
最后,由于f满足
![b816a9c25da1c4347af7e391fbbbd319.png](https://i-blog.csdnimg.cn/blog_migrate/0f0d97c9013e2ca2e37ebdae1c6f811b.jpeg)
所以,优化问题重写为 :
![414840f8f91e220a9eeeeebe8cf52192.png](https://i-blog.csdnimg.cn/blog_migrate/5152cf4365b5ad0e89049f3800afc188.jpeg)
![12bb5a197cf92d2014d4c4d1d4e50d22.png](https://i-blog.csdnimg.cn/blog_migrate/28a40548e44d3aec1c8e4398917292c7.jpeg)
这是一个离散优化问题因为仅允许解向量f中的每一个值只能是两个特殊值,仍是Np难问题。最常用的松弛就是将离散条件改为fi可以取任意实数,松弛后的目标函数为
![a60062e45819bd6ed353853f68a3d9fd.png](https://i-blog.csdnimg.cn/blog_migrate/ee90ab5d7bfd4dad7b296f2b8bcac76f.jpeg)
根据Rayleigh-Ritz理论,问题的解f为拉普拉斯矩阵L的第二小的特征值所对应的特征向量(最小的特征值为0,对应常向量)。所以我们可以通过求L的第二特征向量来解决RatioCut问题。然而我们得到的f是实值向量,要把它转化为离散的指示向量,最简单的方式就是使用sign,
![dff39a4c4369d7e10b0970536196d623.png](https://i-blog.csdnimg.cn/blog_migrate/c2e838f604a47145794c1d67f24d2352.jpeg)
推广到k>2的场景,将V划分成k个子集A1,...Ak,我们定义k个指示向量
![fd40e5c24b76d8c34957c2c99cc9c468.png](https://i-blog.csdnimg.cn/blog_migrate/2b4bd5281812af4b98f349be1e7f0f88.jpeg)
,其中
![90c9b52a299af91629c710c940215b47.png](https://i-blog.csdnimg.cn/blog_migrate/0756dfd7840221fd637038957b419de5.jpeg)
代表每个点的划分情况。定义矩阵
![ec1f7c1b2e226d24bcfcb12b125d98a4.png](https://i-blog.csdnimg.cn/blog_migrate/26a26bbde04520b3c164adcc0be5b476.jpeg)
作为以k个指示向量为列向量的矩阵,显然H中的列相互正交
![4bced000c8789f22633c68882bb39641.png](https://i-blog.csdnimg.cn/blog_migrate/3ca591646f0aa784c050b6326f19c848.jpeg)
,于上面的计算类似,我们可以得到
![7134b38192fd65e43339a5f338278986.png](https://i-blog.csdnimg.cn/blog_migrate/1820e4498faf67978062119cde117309.jpeg)
,进而
![d0e4f163d17a2ef7bc1cfbe15333ad0c.png](https://i-blog.csdnimg.cn/blog_migrate/5d09d97f0c097d5d4e9a199847b2220c.jpeg)
,然后得到RatioCut问题的定义:
![43e80fb8f86385adcde19f0caf4c0252.png](https://i-blog.csdnimg.cn/blog_migrate/0c32012d0e46522a7fde3b3100e6f82c.jpeg)
Tr表示矩阵的迹,所以最小化RatioCut问题可以写作
![817b142403a34e99aeb0a0b62f45b50a.png](https://i-blog.csdnimg.cn/blog_migrate/5121533b789bb9924cb76dd231ad4b07.jpeg)
同样的,我们运行H中的取值为任意实数,松弛后的问题为:
![d4fe521079db26917a7d8384b872b4e1.png](https://i-blog.csdnimg.cn/blog_migrate/f7bc8941603f6a3d812f0d09178aee0b.jpeg)
这是迹极小化问题的标准形式,同样根据Rayleigh-Ritz理论,问题的解就是选择包含L的前k个特征向量作为列而组成的H矩阵。其实矩阵H就是上面非标准化谱聚类算法中提到的矩阵U。然后将实值矩阵转化成离散的形式。最后用kmeans对U的每一行进行聚类。
从NCut的角度解释谱聚类算法:
与在RatioCut中用到的相似的技术同样可以用到规范化谱聚类作为最小化Ncut的松弛。当k=2时我们定义聚类的指示向量f中的每个元素为
![ec6a3e24002c45e14a525699d384557e.png](https://i-blog.csdnimg.cn/blog_migrate/8fc7ad21450487f07240f7f7a1df7d01.jpeg)
和上面类似我们可以检查 ,
![1b81d5fcd726ce746155cf7330f08af3.png](https://i-blog.csdnimg.cn/blog_migrate/27cd1b51844c36d038ae80c5eea66070.jpeg)
,以及
![f46bffd1296d92681d5f023fec4b0ed1.png](https://i-blog.csdnimg.cn/blog_migrate/5a6dfe4b318023400773bfa148d938b9.jpeg)
因此,可以重写最小化Ncut问题为
![1a707ea5631ca7f9ce0dac08ab7bea26.png](https://i-blog.csdnimg.cn/blog_migrate/81a4074dd814be84da2db42fbf243928.jpeg)
然后对问题进行松弛:
![69871d0d89d493cad49cac8908853b4c.png](https://i-blog.csdnimg.cn/blog_migrate/66539513d7f440c937b58945aa5350d4.jpeg)
定义
![3ca5d745bfbee1c24ac868714d240534.png](https://i-blog.csdnimg.cn/blog_migrate/6f4b7b0f40789165fb2cabe2a0e4e5a6.jpeg)
,有
![4a71794e0192ad3d8e1678e469210730.png](https://i-blog.csdnimg.cn/blog_migrate/dc6e30145857b5bab3290df7c95c171b.jpeg)
其中
![59014d3451fc100eee8ce383d7b7d13d.png](https://i-blog.csdnimg.cn/blog_migrate/1f782501a0da9b63b946f980fda1bcb9.jpeg)
,
![0b4fe99e176ce7e30109a21dd2474078.png](https://i-blog.csdnimg.cn/blog_migrate/74ef0cded5e2e97f88f47e11d766751e.jpeg)
是Lsym的第一特征向量,vol(V)是常数。因此,上述目标函数符合Rayleigh-Ritz理论,最优解g就是Lrw的第二特征向量。
推广到k>2的场景,我们定义指示向量
![f9c7bb769c352df191365cf0b29781c6.png](https://i-blog.csdnimg.cn/blog_migrate/ac5a44e6579f7b02842e98b87976d725.jpeg)
,且
![c2d7d007e2f931cf14d928a84b82964f.png](https://i-blog.csdnimg.cn/blog_migrate/248781b62a3aac0f373fa912202e9898.jpeg)
定义H为包含k个指示向量(作为列)的矩阵。显然
![1bd8d8fe53d5edc49a16f8c46fda4549.png](https://i-blog.csdnimg.cn/blog_migrate/71ce22cb77fa56d1b469b6d40359224d.jpeg)
,
![2c1c0d70eae583f85bae7c74938e994c.png](https://i-blog.csdnimg.cn/blog_migrate/c954d14ae7f3ac03bd97947b2ca53a3c.jpeg)
那么最小化Ncut问题可以写作
![53eda68b9e4360431beed3c93e1651c9.png](https://i-blog.csdnimg.cn/blog_migrate/84926fffcb78478353725e48ac9900c9.jpeg)
,松弛后为
![8cb50fd8657449cbd8e6ee5f51690982.png](https://i-blog.csdnimg.cn/blog_migrate/0d3535f21a2f6c1799715e04534b1ff8.jpeg)
这又是一个标准的迹最小化问题,问题的解T是包含Lsym的前k个特征向量为列向量的矩阵。再将
![8e77e3c5f07020bcef23afa3d4cd6355.png](https://i-blog.csdnimg.cn/blog_migrate/210576ec40fbda4109ba9187c09d9148.jpeg)
代入,可以发现H包含矩阵Lrw的前k个特征向量。
谱聚类算法总结
谱聚类的优点 :
1. 对于处理稀疏数据的聚类效果很有效;
2. 使用了降维,在处理高维数据聚类时比传统聚类好;
3. 当聚类的类别个数较小的时候,谱聚类的效果会很好;
4. 谱聚类算法建立在谱图理论上,与传统的聚类算法相比,具有能在任意形状的样本空间上聚类且收敛于全局最优解。
谱聚类的缺点 :
1. 谱聚类对相似图和聚类参数的选择非常敏感;
2. 谱聚类适用于均衡分类问题,即簇之间点的个数差别不大,对于簇之间点的个数相差悬殊的问题不适用。