模式识别-1-非监督学习方法：聚类分析 边肇祺第二版

1、第一章 非监督学习方法：聚类分析根本概念相似性测度与聚类准那么基于试探的聚类搜索算法系统聚类分解聚类动态聚类1.1 根本概念分类与聚类的区别分类：用类别的样本训练集来设计分类器监督学习聚类集群：用事先不知样本的类别，而利用样本的先验知识来构造分类器无监督学习举例：小孩区分桔子和苹果相似性与距离聚类相似性：模式之间具有一定的相似性，这既表现在实物的显著特征上，也表现在经过抽象以后特征空间内的特征向量的分布状态上。聚类分析定义：对一批没有标出类别的模式样本集，按照样本之间的相似程度分类，相似的归为一类，不相似的归为另一类，这种分类称为聚类分析，也称为无监督分类。分类依据：一个样本的特征向量相当于特

2、征空间中的一点，整个模式样本集合的特征向量可以看成特征空间的一些点，点之间的距离函数可以作为模式相似性的度量，并以此作为模式的分类依据。聚类分析是按不同对象之间的差异，根据距离函数的规律进行模式分类的。距离函数的定义特征向量的特性聚类分析的有效性：聚类分析方法是否有效，与模式特征向量的分布形式有很大关系。假设向量点的分布是一群一群的，同一群样本密集距离很近，不同群样本距离很远，那么很容易聚类；假设样本集的向量分布聚成一团，不同群的样本混在一起，那么很难分类；对具体对象做聚类分析的关键是选取适宜的特征。特征选取得好，向量分布容易区分，选取得不好，向量分布很难分开。特征空间维数特征信息的冗余性：在

3、对象分析和特征提取中，往往会提取一些多余的特征，以期增加对象识别的信息量。高维特征空间分析的复杂性：特征空间维数越高，聚类分析的复杂性就越高高维特征空间降维降维方法：相关分析：特征向量的相关矩阵R，分析相关性主成分分析：以正交变换为理论根底独立成分分析：以独立性为根底特征的表示数值表示：对于实际问题，为了便于计算机分析和计算，特征必须进行量化。对不同的分析对象，量化方法是不一样的。连续量的量化：用连续量来度量的特征，只需取其量化值，如长度、重量等。分级量的量化：度量分析对象等级的量，用有序的离散数字进行量化，比方学生成绩的优，良，中，差可用1，2，3，4等量化表示。定性量的量化：定性指标，没有

4、数量关系，也没有次序要求。比方，性别特征：男和女，可用0和1来进行表示。两类模式分类的实例 区分一摊黑白围棋子选颜色作为特征进行分类，用“1代表白，“0代表黑，那么很容易分类；选大小作为特征进行分类，那么白子和黑子的特征相同，不能分类。1.2 相似性测度和聚类准那么 一、相似性的测度欧氏距离: 表征两个模式样本在特征空间中的Euclid距离，模式X和Z间的距离愈小，那么愈相似注意：X和Z的量纲必须一致 消除量纲不一致对聚类的影响：特征数据的正那么化也称标准化、归一化，使特征变量与量纲无关。马氏距离：表征模式向量X与其均值向量m之间的距离平方，C是模式总体的协方差矩阵，引入协方差矩阵，排除了样本

5、之间的相关性。 欧式距离中，如果特征向量中某一分量的值非常大，那么就会掩盖值小的项所起到的作用，这是欧式距离的缺乏；当采用马氏距离，就可以屏蔽这一点。因为相关性强的一个分量，对应于协方差矩阵C中对角线上的那一项的值就会大一些。再将这一项取倒数，减小该影响。当协方差为对角矩阵时，各特征分量相互独立；当协方差为单位矩阵时，马氏距离和欧氏距离相同。其中 分别是样本向量的第k个分量；当m2时，明氏距离就是欧氏距离；当m1时，就是街坊city block)距离： 一般化的明氏距离角度相似性函数：表征了模式向量x和z之间夹角的余弦，反映了几何上的相似性，当坐标系旋转或者尺度变换，夹角余弦测度均保持不变对位

6、移和线性变换不成立如果x和z的分量用二值来表示，0表示不具有某种特征，1表示具有某种特征，那么夹角余弦测度表示x和z具有共有特征数目的相似性测度。二、聚类准那么确实定试探法凭直观和经验，针对实际问题选择相似性测度并确定此相似性测度的阈值，然后选择一定的训练样本来检验测度和阈值的可靠程度，最后按最近邻规那么指定某些模式样本属于某一个聚类类别。举例：对于欧氏距离，它反映了样本间的近邻性，但将一个样本分到不同类别时，还必须规定一距离测度的阈值准那么作为聚类的判别准那么聚类准那么函数法 聚类就是将样本进行组合分类以使类别可分性为最大，因此聚类准那么应是反映类别间相似性或可分性的函数；同时，类别又由一个

7、个样本组成，因此类别的可分性与样本间的差异性直接相关。基于此，聚类准那么函数J，应是模式样本集x和模式类别Sj, j=1,2,c的函数，即J代表了分属于c个聚类类别的全部模式样本与其对应类别模式均值之间的误差平方和；对于不同的聚类形式， J值是不同的，聚类的目的是：使J值到达极小；由此可见：聚类分析转化为寻找准那么函数极值的最优化问题；此种聚类方法通常称为最小方差划分，适用于各类样本密集且数目相差不多，而不同类间的样本又明显分开的情况图例解释把握类内距离与类间距离的问题；聚类准那么函数有许多其他形式。1.3 基于试探的聚类搜索算法一、按最邻近规那么的简单试探法 给N个待分类的模式样本 ，要求按

8、距离阈值T分类到聚类中心算法过程：Step 1：取任意的样本xi作为一聚类中的初始值，如令z1=x1，计算假设D21T，确定一新的聚类中心z2=x2否那么x2以z1为中心的聚类；Step 2：假设已有聚类中心z1和z2，计算 假设D31T和D32T ，那么确定一新的聚类中心z3=x3；Step i： 讨论这种方法的优点：计算简单，假设模式样本的集合分布的先验知识，那么可获得较好的聚类结果。在实际中，对于高维模式样本很难获得准确的先验知识，因此只能选用不同的阈值和起始点来试探，并对结果进行验证。这种方法在很大程度上依赖于以下因素：第一个聚类中心的位置初始化问题待分类模式样本排列次序聚类样本的选择

9、问题距离阈值T的大小判决准那么问题样本分布的几何性质样本的固有特性问题二、最大最小距离算法根本思想：根据实际问题选择距离函数，以试探类间距离为最大作为预选出聚类中心的条件。核心就是：最大类间距离，最小类内距离。算法过程描述：先按照距离最大最小的方法预选出聚类中心，在按照按最邻近规那么将模式分类到聚类中心。对于N个待分类的模式样本 ，要求按最大最小距离法分类到聚类中心 。Step 1：选任意一模式样本xi作为第一聚类中心z1Step 2：选离z1最远距离的样本xj作为第二聚类中心z2Step 3：逐个计算各模式样本 与 之间的距离，并选出其中的最小距离。Step 4：在所有模式样本的最小值中选出

10、最大距离，假设该最大值到达 的一定分数比值以上，那么将相应的样本取为第三聚类中心。Step i： 算法性能分析：算法复杂度增加，在选聚类中心过程中消耗较大的资源。1.4 系统聚类系统聚类：先把每个样本作为一类，然后根据它们间的相似性或相邻性聚合，类别由多到少，直到获得适宜的分类要求为止；相似性、相邻性用距离表示。聚合的关键就是每次迭代中形成的聚类之间以及它们和样本之间距离的计算，不同的距离函数会得到不同结果。两类间距离计算准那么：(注意理解1. 最短距离：两类中相距最近的两样本间的距离 2. 最长距离 ：两类中相距最远的两个样本间的距离。 3. 中间距离：最短距离和最长距离都有片面性，因此有时

11、用中间距离。设1类和23类间的最短距离为d12，最长距离为d13， 23类的长度为d23，那么中间距离为：上式推广为一般情况：4. 重心距离：均值间的距离5. 类平均距离：两类中各个元素两两之间的距离平方相加后取平均值 6. 离差平方和：设N个样本原分q类，那么定义第i类的离差平方和为：离差平方和增量：设样本已分成p,q两类，假设把p,q合为r类，那么定义离差平方增量：算法过程描述：Step1：初始距离矩阵的计算D(0)说明：1距离矩阵元素的值是类与类之间的距离，距离的定义有多种。2距离矩阵，是对称矩阵。对角上线的元值表示同类之间的距离，即为0。Step2：对于第n次迭代的距离矩阵D(n)进行

12、聚合 说明：距离矩阵中选择距离最小的，如果有相同的可以任选其中一个，要忽略对角线上的元素；也可以把相同的全部聚合。Step3：根据第n次聚合结果，计算合并后的新类别之间的距离矩阵D(n+1) 说明：合并类的距离计算应该符合距离的运算规那么。假设距离反映的是两类的重心距离，那么合并后，应该仍然反映的重心的距离。Step4：收敛性判决距离阈值D的设定 说明：算法的收敛条件判断准那么确实定。例1：如以下图所示简单的一维情况 1、设全部样本分为6类，2、计算距离矩阵D(0)1234561029031160449166405254364066425811903、求最小元素：4、把1, 3合并7=(1,3

13、)4, 6合并8=(4,6)5、作距离矩阵D(1)，按最小距离准那么7285702908491605254406、假设合并的类数没有到达要求，转3。否那么停止。3、求最小元素：4、8, 5, 2合并, 9=2,5,4,61.5 分解聚类分解聚类：把全部样本作为一类，然后根据相似性、相邻性分解。目标函数为：两类均值方差 N：总样本数， ：1类样本数 ：2类样本数，分解聚类框图初始分类调整分类方案最终结果目标函数达到最优？例2：21个样本，每个样本取二个特征，原始资料矩阵如下表：样本号12345678910 x10022445667x26553431210111213141516171819202

14、1-4-2-3-3-5100-1-1-3322021-1-2-1-3-5目标函数解：第一次分类时计算所有样本，分别划到 时的E值，找出最大E值对应的样本。 1、开始时， 2、分别计算当 划入 时的E值把 划入时有 然后再计算把 划入 时对应的E值，找出一个最大的E值。 一直计算下去 把 划为 的E值最大。 再继续进行第二，第三次迭代计算出 E(2) , E(3) , 次数 E值 第10次迭代 划入 时，E最大。于是分成以下两类：每次分类后要重新计算 的值。可用以下递推公式：1.6 动态聚类兼顾系统聚类和分解聚类一、动态聚类方法概要 先选定某种距离作为样本间的相似性度量; 确定评价聚类结果的准那

15、么函数; 给出某种初始分类，用迭代法找出使准那么函数取极值的最好的聚类结果。选代表点初始分类分类合理否最终分类修改分类YN动态聚类框图 二、代表点种子点的选取方法：代表点就是初始分类的聚类中心数K 凭经验选代表点，根据问题的性质、数据分布，从直观上看来较合理的代表点K; 将全部样本随机分成K类，计算每类重心，把这些重心作为每类的代表点; 用前K个样本点作为代表点。 按密度大小选代表点： 以每个样本作为球心，以d为半径做球形；落在球内的样本数称为该点的密度，并按密度大小排序。首先选密度最大的作为第一个代表点，即第一个聚类中心。再考虑第二大密度点，假设第二大密度点距第一代表点的距离大于d1人为规定

16、的正数那么把第二大密度点作为第二代表点，否那么不能作为代表点，这样按密度大小考察下去，所选代表点间的距离都大于d1。d1太小，代表点太多，d1太大，代表点太小，一般选d12d。对代表点内的密度一般要求大于T。T0为规定的一个正数。三、初始分类和调整 选一批代表点后，代表点就是聚类中心，计算其它样本到聚类中心的距离，把所有样本归于最近的聚类中心点，形成初始分类，再重新计算各聚类中心，称为成批处理法。 选一批代表点后，依次计算其它样本的归类，当计算完第一个样本时，把它归于最近的一类，形成新的分类。再计算新的聚类中心，再计算第二个样本到新的聚类中心的距离，对第二个样本归类。即每个样本的归类都改变一次

17、聚类中心。此法称为逐个处理法。 直接用样本进行初始分类，先规定距离d，把第一个样品作为第一类的聚类中心，考察第二个样本，假设第二个样本距第一个聚类中心距离小于d，就把第二个样本归于第一类，否那么第二个样本就成为第二类的聚类中心，再考虑其它样本，根据样本到聚类中心距离大于还是小于d，决定分裂还是合并。最正确初始分类：如下图，随着初始分类K的增大，准那么函数下降很快，经过拐点A后，下降速度减慢。拐点A就是最正确初始分类。四、K-均值算法：成批处理法例：有20个样本，每个样本有2个特征，数据分布如以下图第一步：令K=2，选初始聚类中心为样本序号x1x2x3x4x5x6x7x8x9x10特征x1010

18、1212367特征x20011122266x11x12x13x14x15x16x17x18x19x2086789789896777788899因为所以，同理同样，把其他所有样本 与两个聚类中心的距离计算出来，可判断得因此，可分类两大类第三步：根据新分成的两类建立新的聚类中心第四步：转第二步。第二步：重新计算 到Z1(2) , Z2(2) 的距离，把它们归为最近聚类中心，重新分为两类，第三步，更新聚类中心第四步，第二步，第三步，更新聚类中心说明： (1)K是指需要分成K类，均值是指每类的中心,就是该类所有样本的平均值，不一定就有某个样本在这个位置上。 (2)算法的收敛性判别：前后两次迭代的结果，也就是每迭代分类后，分类都是一样的，此时停止。 (3)K值和初始聚类中心对分类的结果影响很大。通常需要其它的算法来确定这两个的选取