一种簇特征加权的模糊紧致散布聚类方法

摘要

针对现有的WFCM算法在聚类时没有考虑样本硬划分实际情况且对样本分布不均衡情况数据聚类效果不好，FCS算法没有考虑硬划分边界点的情况以及忽略样本特征参数对各类聚类影响的问题，本发明公开了一种簇特征加权的模糊紧致散布聚类方法,通过对样本隶属度、特征权重进行调整，遵循了样本硬划分的实际情况，并充分考虑样本特征参数对各类聚类的影响，尽可能使得样本类内紧致、类间分散，解决了位于硬划分边界的样本隶属度问题，在样本分布不均衡情况下对于噪声数据和异常数据实现了更有效的划分。本方法聚类性能良好，收敛速度快、迭代效率高。聚类性能良好，迭代效率高，适于应用在工业控制中样本分布不均衡、实时性、准确率要求高的场合。

主权项

一种簇特征加权的模糊紧致散布聚类方法，其特征在于，包括如下步骤：步骤一：设置隶属度指数m、特征加权指数α∈[‑10，‑1]∪(1,10]、β∈{0.005,0.05,0.5,1}，初始迭代次数p＝0以及迭代误差ε＞0，随机生成初始聚类中心ai，(s为特征参数个数)；步骤二：根据下式计算系数η_i：${\eta }_i=\frac{\beta }{4}\frac{{\min }_{i\ne i^{\prime }}{\left|\right|a_i-a_{i^{\prime }}\left|\right|}^2}{{\max }_t{\left|\right|a_t-\bar{X}\left|\right|}^2}$其中，为样本均值；步骤三：根据下式更新样本隶属度μ_ij：${\mu }_{\mathrm{ij}}=\frac{{\left(\underset{k=1}{\overset{s}{\Sigma }}{\omega }_{\mathrm{ik}}^{\alpha }({\left|\right|x_{\mathrm{jk}}-a_{\mathrm{ik}}\left|\right|}^2-{\eta }_i{\left|\right|a_{\mathrm{ik}}-\overline{X_k}\left|\right|}^2)\right)}^{\frac{1}{1-m}}}{\underset{t=1}{\overset{c}{\Sigma }}{\left(\underset{k=1}{\overset{s}{\Sigma }}{\omega }_{\mathrm{tk}}^{\alpha }({\left|\right|x_{\mathrm{jk}}-a_{\mathrm{tk}}\left|\right|}^2-{\eta }_t{\left|\right|a_{\mathrm{tk}}-\overline{X_k}\left|\right|}^2)\right)}^{\frac{1}{1-m}}}$记${\Delta }_{\mathrm{ij}}=\underset{k=1}{\overset{s}{\Sigma }}{\omega }_{\mathrm{ik}}^{\alpha }({\left|\right|x_{\mathrm{jk}}-a_{\mathrm{ik}}\left|\right|}^2-{\eta }_i{\left|\right|a_{\mathrm{ik}}-\overline{X_k}\left|\right|}^2)$当样本点x_j存在落在硬划分边界上时，此时Δ_ij＝0，在保证各样本点相对于第i类的距离尺度不变的前提下，对所有Δ_ij≥0的利用P(Δ_ij)进行调整：${\Delta }_{\mathrm{ij}}=P({\Delta }_{\mathrm{ij}}\ge 0)={\Delta }_{\mathrm{ij}}+\mathrm{rand}*\underset{j}{\min }({\Delta }_{\mathrm{ij}}>0)(j=1,...,n)$调整后利用下式计算新的μ_ij：${\mu }_{\mathrm{ij}}=\frac{{{\Delta }_{\mathrm{ij}}}^{\frac{1}{1-m}}}{\underset{t=1}{\overset{c}{\Sigma }}{{\Delta }_{\mathrm{tj}}}^{\frac{1}{1-m}}}$因为有样本点x_j落在第i类硬划分区域内，所以会有μ_ij＜0，因此对μ_ij进行硬划分调整：$\left\{\begin{array}{cc}{\mu }_{\mathrm{ij}}=1,& {\Delta }_{\mathrm{ij}}<0\\ {\mu }_{i^{\prime }j}=0,& i^{\prime }\ne i\end{array}\right.$步骤四：根据下式计算特征权重ω_ik：${\omega }_{\mathrm{ik}}=\frac{{\left(\underset{j=1}{\overset{n}{\Sigma }}{\mu }_{\mathrm{ij}}^m({\left|\right|x_{\mathrm{jk}}-a_{\mathrm{ik}}\left|\right|}^2-{\eta }_i{\left|\right|a_{\mathrm{ik}}-\overline{X_k}\left|\right|}^2)\right)}^{\frac{1}{1-\alpha }}}{\underset{t=1}{\overset{s}{\Sigma }}{\left(\underset{j=1}{\overset{n}{\Sigma }}{\mu }_{\mathrm{tj}}^m({\left|\right|x_{\mathrm{jt}}-a_{\mathrm{it}}\left|\right|}^2-{\eta }_i{\left|\right|a_{\mathrm{it}}-\overline{X_t}\left|\right|}^2)\right)}^{\frac{1}{1-\alpha }}}$记${\Delta }_{\mathrm{ik}}=\underset{j=1}{\overset{n}{\Sigma }}{\mu }_{\mathrm{ij}}^m({\left|\right|x_{\mathrm{jk}}-a_{\mathrm{ik}}\left|\right|}^2-{\eta }_i{\left|\right|a_{\mathrm{ik}}-\overline{X_k}\left|\right|}^2)$若Δ_ik＜0，因为ω_ik∈[0,1]，所以需将Δ_ik投影到大于0的区间且保证各样本的第k个特征参数与第i类的硬划分区的距离尺度不变，于是利用下式调整Δ_k：${\Delta }_{\mathrm{ik}}={\Delta }_{\mathrm{ik}}-\underset{k}{\min }\left({\Delta }_{\mathrm{ik}}\right)+\underset{k}{\min }({\Delta }_{\mathrm{ik}}>0)$调整后利用特征权重公式计算新的ω_ik；步骤五：根据下式计算聚类中心a_ik：$a_{\mathrm{ik}}=\frac{\underset{j=1}{\overset{n}{\Sigma }}{\mu }_{\mathrm{ij}}^m(x_{\mathrm{jk}}-{\eta }_i\overline{X_k})}{\underset{j=1}{\overset{n}{\Sigma }}{\mu }_{\mathrm{ij}}^m(1-{\eta }_i)}$步骤六：令迭代次数p＝p+1，直到否则转到步骤二；步骤七：将第p次迭代得到的μ_ij输出，根据即第j个样本属于第i类。