[发明专利]基于局部密度估计和近邻关系传播的多路谱聚类方法

权利要求书

1.一种基于局部密度估计和近邻关系传播的多路谱聚类方法，包括如下步骤：

(1)输入数据集X＝{x₁，x₂，...，x_n}∈R^d，x_n表示数据集中的第n个样本，n为样本个数，d为样本维数；

(2)估计样本的局部密度：

(2a)寻找样本x的K个最近邻样本，并构成集合N(x)＝{y₁，y₂，...，y_k}∈R^d，其中，y_K表示样本x的第K个最近邻样本；

(2b)计算x的距离集合：D(x)=∪i=1Kd(yi,Nearest(yi)),]]>其中，Nearest(y_i)表示样本y_i的最近邻样本，d(y_i，y_j)表示y_i与y_j的欧氏距离：

d(yi,yj)=(yi1-yj1)2+(yi2-yj2)2+...+(yid-yjd)2]]>

其中，y_id和y_jd分别表示第i个样本和第j个样本的第d维属性值；

(2c)定义样本集Near(x)：先判断D(x)中元素是否可以分成两类，若可分，设元素个数较多的一类样本表示为y′₁，y′₂，...，y′_m，其中，m表示样本个数，则认为这m个样本对x的密度估计有较大影响，定义Near(x)＝{y′_i|i＝1：m}，否则，K个样本整体对x的密度估计都有较大影响，定义Near(x)＝N(x)；

(2d)计算x的判别密度估计

f^P(x)=Σxi∈Near(x)e-d(x,xi)22σ2]]>

其中，Near(x)根据步骤(2c)可得，d(x，x_i)为样本x与x_i的欧氏距离，σ²为窗宽；

(2e)计算x的单独密度估计

f^S(x)=e-d(x,y3)22σ2]]>

其中，y₃为x的第3个近邻样本；

(2f)计算x的局部密度估计定义局部密度估计为判别密度估计与单独密度估计之和，即f^L(x)=f^P(x)+f^S(x);]]>

(3)样本密度升维：将局部密度作为第d+1维扩展到样本x中，构成新的数据集X*={x1*,x2*,...,xn*}∈Rd+1,]]>其中，xn*[xnf^L(x)]T;]]>

(4)计算X^*样本之间的欧氏距离，得到距离矩阵B，即B＝[b_ij]_n×n；

(5)计算X^*样本之间的相似度，得到相似度矩阵W＝[w_ij]_n×n，w_ij＝exp(-b_ij²/2σ²)，其中，w_ij表示第i个样本和第j个样本之间的相似性；

(6)根据距离矩阵B计算阈值ε，并确定B中元素对应样本之间的近邻关系，获得初始的近邻关系矩阵T，其中，

(7)根据近邻关系传播原则分别更新T和W，得出亲合矩阵A；

(8)构建度矩阵D和拉普拉斯矩阵Lsym，其中，D为对角矩阵，对角元素表示第i个样本x_i的度，Lsym＝D^-1/2AD^-1/2；

(9)计算Lsym的前k个最大特征值所对应的特征向量，并构成矩阵V，然后进一步单位化为矩阵Y＝[y_ij]_n×k，其中，

(10)将Y的每一行作为k维空间中的一个样本点，通过K-means算法将这些样本点聚成k类；当且仅当Y的第i行被分配为第j类时，将样本x_i分配为第j类。