基于马氏距离的模糊聚类优化算法——KM-FCM

聚类是数据处理的重要方法。模糊聚类建立了数据样本对于类别的不确定性的描述，表达了样本类属的模糊性，能够更客观地反映现实世界，具有较强的聚类效果与数据表达能力。

关于模糊聚类的研究应用最广的是基于目标函数的模糊聚类方法，此方法把聚类问题描述为一个带约束的优化问题，通过求解优化问题的解来确定数据集的模糊划分和聚类结果。FCM算法是最经典的基于目标函数的模糊聚类算法，但也存在一些问题，因此基于FCM算法的改进算法和应用被研究者所关注。蔡静颖[1]为改进FCM算法不能处理非球形簇且未考虑样本矢量各特征重要程度、处理高相关数据集时错误率增加的缺点，使用马氏距离得到FCM-M算法及MF-FCM算法；蔡威[2]提出自适应距离度量的GK算法。NATACHA等[3]研究了用马氏距离和闵可夫斯基距离来取代欧氏距离的模糊聚类的方法，以提高聚类检测能力，并对聚类结果进行了可视化分析。在应用方面，张敏等[4]将马氏距离和模糊c-均值聚类结合，研究抠图算法；康韦晓[5]将基于马氏距离的PFCM算法应用于非线性系统故障诊断；赵泉华等[6]将马氏距离的模糊聚类算法用于遥感图像分割。另外，通过对FCM的初始化方法的改进研究逐渐形成了山峰函数法或势函数法[7]、减法聚类[8]等方法，进一步实现了快速减法聚类的模糊聚类算法，及基于密度函数的近似初始化方法[9]。上述文献中基于马氏距离的模糊聚类算法的研究多是对协方差估算方面的改进，并且已有的初始化方法对基于马氏距离的模糊聚类算法并不适用。为此，本研究通过对经典聚类算法和马氏距离特性的分析，提出具有新初始化方法的基于马氏距离模糊聚类的优化算法。

1 经典模糊聚类算法与马氏距离

1.1 经典模糊聚类算法(FCM)

在基于目标函数的模糊聚类算法中，模糊c-均值算法(fuzzy c-means，FCM)的理论最为完善，应用最为广泛。FCM类型的算法最早是由“硬”聚类算法HCM导出的[10]，DUNN[11]把它的目标函数扩展到隶属度属于模糊情形的类内加权平均误差和函数后来BEZDEK[12]又引入了一个参数m，把J2推广到一个目标函数的无限簇，并给出了交替优化算法，即形成了经典的模糊c-均值算法。

FCM算法的核心思想如下：设X={x1,x2,…,xi,…,xn}为n元数据集合。FCM聚类方法就是把X划分为c个子集S1，S2,…,Si,…,Sc,若用V={v1,v2,…,vi,…,vc}表示这c个子集的聚类中心，uij表示元素xj对Si的隶属度，则FCM算法的优化目标函数为

 

(1)

uij满足如下约束条件：

 

(2)

这里U={uij}为c×n矩阵，V={v1,v2,…,vi,…,vc}为s×c矩阵，dij为xj与vj的距离，经典的FCM算法里使用的是欧氏距离，推荐使用m=2。通过如上设定，最佳聚类结果可使得目标函数取得极小值。

FCM算法的具体步骤如下：

步骤1 设定聚类个数c(1<c<n)和模糊指数m(1≤m<+∞)；初始化各类中心V0；设置收敛的精度ε>0；令迭代次数L=0，给出迭代数目的最大值Lmax。

步骤2 根据来计算U(L+1)。

4）A.He thinks it would be really fun. B.He will definitely join a walking tour.

(3) 可扩充性原则.可扩充性是模块设计的重要优势.特别是在航天发射场,模块化设计单元应考虑未来保障新型号发射任务的可能性,应预留足够的接口和空间,以便进一步发展.

步骤4 判断2次聚类中心矩阵的欧氏距离与给定阈值的大小，如果满足终止条件：‖V(L)-V(L-1)‖≤ε,L≥1，或者迭代次数不小于给定的最大迭代次数，则迭代停止，否则，重复步骤2和步骤3。

在对比实验中，FCM聚类的运行相对于HCM算法速度较慢，但聚错样本数明显减少，并且该算法的收敛性已经得以证明。但是由于经典模糊聚类算法就是反复修改聚类中心矩阵和隶属度矩阵的分类过程，并且FCM对聚类中心的初始化依赖，使得经典模糊聚类算法不能确保得到全局最优解。另一方面，与HCM算法一样，需要预先指定聚类数目c，而实际中聚类数目通常都是未知的，并且使用欧氏距离只适于发现球状类型等非凸面形状的簇，不能处理椭球形结构簇，在处理高维数据时效果欠佳。

1.2 马氏距离

3) 用式(3)计算隶属矩阵U；

设样本集合为X={x1,x2,…,xi,…,xm}，共有m个样本，xi是n维特征矢量，i∈{1,2,…,m},令X代表m×n的输入矩阵，每行为一个样本，则样本的均值、自相关矩阵和协方差矩阵可用矩阵表示为

步骤3 用计算聚类中心矩阵V(L+1)，令L=L+1。

 

其中代表元素均为的m维列矢量。

样本Xi到样本总体X的马氏距离定义为

由于对于聚类样本的总体分布通常是未知的，而样本协方差是总体协方差的无偏估计，常用样本协方差矩阵代替总体协方差：若协方差矩阵是奇异的，将导致无法直接求马氏距离，故大量文献对基于马氏距离的研究多是对计算其数据协方差的逆不存在的情况提出解决方案，王振丽[17]为改进马氏距离使用加权MP马氏距离进行研究；吴香华等[18]对马氏距离聚类分析中协方差矩阵的估算进行改进；赵小强等[19]通过改进协方差矩阵的估计来提高马氏距离聚类分析效率。

解读这则史料的关键在“庶人工商遂，人臣隶圉免”这句话。“庶人工商遂”，杜注：“得遂进仕”。“人臣隶圉免”，杨伯峻先生参武亿《义证》注云：

基于马氏距离的聚类算法有如下特点。

1)协方差矩阵本身的意义是在多维向量之间找出一个自适应的权重，即马氏距离的计算是建立在总体样本基础上的，有利于加强聚类的准确性。

2)在计算马氏距离过程中，要求总体样本个数大于样本的维数，否则得到的总体样本协方差逆矩阵不存在。

再比如说擦地板。我家两间卧室铺的木地板，过去一直是妻子拿抹布擦，像个日本女人一样，跪在地板上一块一块地擦。不是每天擦一遍，最起码三天就得擦一遍。妻子说上拖把拖，拐拐角角的哪能拖干净？现在妻子不拿抹布擦地板，就算把拖地的任务都交给了我。

第1步：求出样本集合O中两两样本之间的距离，找出距离最小的一对较小样本x1,x2；

3)在实际应用中，“总体样本个数大于样本的维数”这个条件很容易满足，即在绝大多数的情况下，马氏距离是可以顺利计算的，从而有关基于马氏距离的模糊聚类的优化可以先忽略协方差问题，针对高维大样本的数据做出优化的改进算法。但是马氏距离的计算是不稳定的，不稳定来源于协方差矩阵，这也是马氏距离与欧氏距离的最大差异之处，故在研究中应注意结合用加强稳定性的优化方案来弥补马氏距离的不足。

今人对黄庭坚词集的整理颇为重视，版本众多。然而，每个版本所收黄庭坚词的数量却多寡不同。不计断句，仅就篇幅完整的词作进行统计，则《全宋词》和《全宋词补辑》二者合计189首；龙榆生校本《豫章黄先生词》收词181首；严寿澂校点《山谷词》收词187首；刘琳等《黄庭坚全集》收词184首，又在补遗中收录11首；朱德才《增订注释秦观黄庭坚词》收录黄词191首；马兴荣、祝振玉校《山谷词校注》收词179首，另有补遗9首，存疑18首。这种状况源于对部分词作著作权的确认辨别有异。这说明，对部分词作著作权的归属问题进行辨识仍有待于学界努力。以下笔者尝试对此作简要考辨。

其次，在交替优化方面，结合新的聚类初始化方法，并引用文献[9]中的利用粒度分析原理的GD有效性准则函数在算法循环中合并聚类方法，得到基于马氏距离的模糊聚类优化算法——KM-FCM。

2 模糊聚类优化算法——KM-FCM

2.1 新的聚类初始化方法

由于基于山峰函数和减法聚类的初始化方法的实现性并不令人满意，近几年对初始化方法的研究出现了结合进化及生物仿生算法进行的初始化方法，李静[20]结合粒子群算法初始聚类中心，使得聚类结果以较快收敛速度接近最优解，并且无需给定聚类数目，但其缺点是引入了其他参数；NAIK等[21]使用基于教学学习(TLBO)的优化算法，聚类结果以较快收敛速度接近最优解，但是局限于欧氏距离且需要给定聚类数目。

宫颈癌 出现同房后出血，停经多年又突然阴道出血，白带增多，呈血性或洗肉水样，下腹部及腰部疼痛，要警惕。

“在舞者前，我放置了一个有多只闪光灯的导光栅格条形柔光箱。舞者的右后方有一束在低位定点标记的持续光源。我们要求她在室内做了多次同样的动作，并用同步后视图捕获影像，如此一来，快门释放后曝光的瞬间就能闪光。我们选择了几张图片，并用PS进行了合成。”

考虑到基于马氏距离的模糊聚类算法比经典的模糊聚类算法运行速度慢，为减少算法速度负担，在此提出的初始化方法也是将k-means聚类算法结果作为模糊聚类的初始聚类中心。由于k-means对聚类中心初始化敏感，易陷入局部最优且需给定聚类数，故给出基于启发式搜索算法与k-means算法结合的聚类初始化方法，对基于马氏距离的模糊算法进行优化。

新的初始化方法主要思想如下。为实现全局优化并且避免重复初始聚类中心，首先考虑在一定的聚类个数范围内用启发式搜索聚类中心，然后利用k-means聚类算法得出初始聚类中心。最佳聚类数通常满足为样本数据个数，本研究也作这样的设定。

新的初始化方法的具体步骤如下。

订书器，为我们提供了把许多页纸装订在一起加以保存的可能。一般订书器都需要用订书钉，而这款无针订书器就显得与众不同了。你只需轻轻一压，文件就被钉在一起了，纸上只留下一小排美丽的波浪形装订线。它的高科技之处在于独特的按压锁扣技术。人手在向下按时，其金属齿能瞬间产生200千克的力，把局部的纸张压成波浪形，如此一来，纸张就被牢牢地固定在一起了。更加巧妙的是，只要用订书器的底部把装订线磨平，就又能把订好的纸张分开了。

三十三天为额尔和木·哈日求情之后，又给他出谋划策，告诉他如何战胜哈冉惠，而哈冉惠放过额尔和木·哈日之后，说出了这样的话：

第2步：把x1,x2放到新的集合A1中，并将它们从样本集合O中删除；

还有一种聚类初始化方法是将已有的复杂度低的聚类算法结果作为模糊聚类的初始聚类中心，如直接运用k-means算法选出初始聚类中心[9]；结合概率抽样方法的k-means++算法[22-23]是比前者初始效果更好的聚类初始化方法，但引入了一个缺乏理论支撑的参数p。

第3步：计算A1中样本的均值a1，求出a1与O中每个样本的距离，找出O中与a1最近的较小样本x3，将x3并入集合A1，且从O中删除，如此重复，直到A1中样本个数达到的向下取整值；

(1)教学经费及训练办公场所.据了解，目前民族武术社同样没有固定的办公地点和训练场所，包括开会地点也是由社团负责人自行安排，若负责人找不到合适的地点开会，就由社团成员一起想办法解决.而训练场所也是由传承人临时找地方，如公园、广场、街道等地方都是最常用的训练地点.同时，通背拳研究会的教学情况也并不乐观，研究会的多数成员私下都会组织授课，但都是以个人名义进行，真正以研究会名义打开局面进行传承的暂时没有.

第4步：再把O中现有样本中距离最小的一对较小样本放入新的集合A1中，并将它们从O中删除，重复第3步的过程，直到形成m个集合A1,A2,…,Am,m预设成比大的值，n为数据集样本的个数，这里设定成

第5步：分别求出A1,A2,…,An的均值对数据集运用k-means算法，得到的个聚类中心作为用于模糊聚类算法的初始聚类中心。

针对每个聚类数目进行模糊聚类时，都要重新初始化所造成的聚类数目不稳定的情况，经过算法迭代后，采用合并聚类中心的方式。这使得聚类个数自适应，无需给定聚类数目c。

2.2 KM-FCM算法

首先，在聚类相似性准则方面，在经典FCM算法的目标函数中用马氏距离替代欧氏距离，并且在目标函数上引进一个协方差调节因子：得到的优化目标函数为

 

该优化问题的拉格朗日乘子式为

 

最小化J，对ci,uij，Σi求偏导，并令其结果等于零，得:

 

(3)

 

(4)

 

针对以上问题，在基于马氏距离的模糊聚类算法的研究基础上，本研究并不侧重对马氏距离协方差问题的优化，而是另辟蹊径，提出了新的初始化优化方法，形成一种具有全局优化性能、利于处理高维数据的、新的基于马氏距离的模糊聚类算法。

1) 取定最大初始聚类个数为样本数据个数，模糊加权指数m=2，迭代停止阈值Lmin=1×10-5，权重因子α=0.6，最大迭代次数Lmax=100，并令迭代计数器L=0；

2) 由设定的初始聚类个数cmax，运用新的初始化方法，得到初始聚类中心结果矩阵C并输出；

欧氏距离只适用于样本的属性在相互独立的条件下同等对待每个属性对聚类影响的情形。当属性之间相关时，对相关的属性计算欧氏距离时将产生重复数据，影响了聚类效果和最佳聚类数的确定。同时，欧氏距离受属性量纲的影响，对多维数据的处理是不利的。针对样本向量中各维特征对模式分类的不同影响，李洁等[13]提出了基于特征加权的模糊聚类新算法，但收敛速度有所下降。用马氏距离来取代欧式距离，可以有效解决以上困扰。马氏距离是一种有效计算2个未知样本集相似度的方法，与欧氏距离不同的是它考虑到各种特性之间的联系，并且是与尺度无关的。在文献[14]中说明了当聚类算法在用于入侵检测[15]时，马氏距离比欧式距离具有明显的优势，并且马氏距离的模糊聚类在图像分割上也比欧氏距离的效果更优[16]。

4) 用式(4)更新聚类中心矩阵C；

因此，学生选择的几个关键词，最好能拥有一个共同的特征，或者说是有意识地选择那些能够体现一个特点的关键词。我们教师指导学生写作的时候，首先要想想主题词的特征是什么，哪些关键词能展现这个特征，然后再去思考这些关键词有没有内容可写，这样作文自然有话可讲了。

5) 若聚类中心达到迭代停止阈值，则输出模糊分类隶属矩阵U和聚类中心C，若未达到迭代停止阈值，则令L= L+1，转向步骤3)；

6)计算有效性函数

其中:

 

将GD值保存起来；

7)将类间两两之间的距离最小的两类合并为一类，得到c-1个聚类中心；

8)令c=c-1，若c<2，则转向步骤9)，否则转向步骤3)，直到达到最大迭代次数Lmax=100；

9)选择最小值GD对应的聚类结果，即为最佳聚类结果，算法结束。

3 结果与分析

为了证明应用本研究提出的初始化优化后的KM-FCM算法比FCM算法及未优化的基于马氏距离的模糊聚类算法FCM-M[1]，M-FCM[19]在多属性的样本聚类实验中具有优势，笔者将FCM，FCM-M，M-FCM，KM-FCM算法应用于来自UCI数据库的Iris，Wine和Pima等3个标准数据集。Iris数据集由4维空间的150个样本组成，分为3个类别，第1类与其他2类完全分离，而第2类与第3类之间有交叉，有很多学者认为Iris数据也可以分为两类[24-25]；Wine数据由13维空间的178个样本组成，分为3个完全分离的类别；Pima数据集由8维空间的768个样本组成，分为2个相互交叠的类别。

分别使用FCM，FCM-M，M-FCM聚类算法与KM-FCM聚类算法对Iris数据各进行20次实验，取20次实验的平均值。其中，FCM，FCM-M，M-FCM聚类算法设定聚类类别数，4种算法实验中相同参数设置如下：隶属度矩阵的指数m=2，最大迭代次数Lmax=100，迭代停止准则为Lmin=1×10-5。对3种数据集的4种聚类算法比较分别如表1，表2和表3所示。

由表1可知，KM-FCM算法对低维属性的Iris数据集进行聚类是有效的；由表2可知，在高维属性的模糊聚类中，FCM的聚类效果明显下降，而其他基于马氏距离的模糊聚类算法聚类精度较高，KM-FCM算法聚错个数明显减少；由表3可知，KM-FCM算法在有交叠的、类边界并不清晰的较大样本Pima数据集上聚类效果依然良好。对于3个不同的数据集，从时间上来看，由于文献[1]中的FCM-M算法并未使用优化技巧，运行速度大于文献[19]中用到数据标准化处理、总体协方差估计优化方法的M-FCM算法以及本研究提出的KM-FCM算法。但是，经过初始化优化的KM-FCM算法聚类精度比其他3种聚类算法的聚类精度都要好，误分个数少，具有全局优化效果，并且无需设定聚类个数。在20次实验中，由于马氏距离的不稳定性，并且未进行协方差奇异问题处理，KM-FCM算法在有交叠样本数据中的聚类稳定性有待进一步研究。

沥青混凝土路面施工是连续作业，各工序联系紧密。同时在施工现场存在着各种不可预见的因素和风险，所以在人员配备上需要合理安排，管理人员和技术人员要均衡分组。通常采取双班或三班制，在关键工序上要配备责任心强、技术水平高的员工。

 

表1 对于Iris数据集的20次实验4种聚类算法比较Fig.1 Comparison of the four clustering algorithms for 20 experiments on Iris datasets

  

算法误分个数平均耗时/s聚类精度/%FCM160．163789．33FCM⁃M80．210594．67M⁃FCM140．231491．23KM⁃FCM70．220895．33

 

表2 对于Wine数据集的20次实验4种聚类算法比较Fig.2 Comparison of the four clustering algorithms for 20 experiments on Wine datasets

  

算法误分个数平均耗时/s聚类精度/%FCM920．825848．31FCM⁃M480．913173．03M⁃FCM490．924372．03KM⁃FCM410．938276．97

 

表3 对于Pima数据集的20次实验4种聚类算法比较Fig.3 Comparison of the four clustering algorithms for 20 experiments on Pima datasets

  

算法误分个数平均耗时/s聚类精度/%FCM3160．902158．89FCM⁃M2770．975663．93M⁃FCM2380．984268．90KM⁃FCM2330．996369．66

4 结 语

通过对经典聚类算法和马氏距离特性的研究，以及观察关于初始化方法改进的最新动态，本研究给出了适用于马氏距离模糊聚类算法的初始化方法，提出基于马氏距离模糊聚类的优化算法KM-FCM。用马氏距离替换经典的模糊聚类算法中的欧氏距离，采用启发式搜索与k-means算法结合的初始化方法，与FCM，FCM-M，M-FCM等3种算法在3个标准数据集上进行仿真对比实验，验证了新的算法的有效性和全局优化作用。研究结果可为基于马氏距离的模糊聚类算法的优化提供参考。

参考文献/References：

[1] 蔡静颖.模糊聚类算法及应用[M].北京：冶金工业出版社，2015.

[2] 蔡威. 模糊聚类算法在数据挖掘中的应用研究[D]. 兰州：兰州交通大学,2012.

CAI Wei. Research on the Application of Fuzzy Clustering Algorithm in Data Mining[D]. Lanzhou：Lanzhou Jiaotong University,2012.

[3] NATACHA G, IREN V, GEORGE G, et al. Fuzzy C-means clustering with mahalanobis and minkowski distance metrics[J]. Procedia Computer Science, 2017,114: 224-233.

[4] 张敏,闵乐泉,张群,等. 基于马氏距离和模糊C均值聚类的抠图算法与应用[J].北京科技大学学报,2014,36(5):688-694.

ZHANG Min, MIN Lequan, ZHANG Qun, et al. Matting algorithm and application based on Mahalanobis distance and the fuzzy C-means clustering algorithm[J]. Journal of University of Science and Technology Beijing, 2014,36(5):688-694.

[5] 康韦晓. 基于马氏距离的PFCM算法的非线性系统故障诊断方法[D]. 哈尔滨：哈尔滨工业大学,2016.

KANG Weixiao.Fault Diagnosis Method for Nonlinear System based on PECM Algorithm with Mahalanobis Distance[D]. Harbin：Harbin Institute of Technology, 2016.

[6] 赵泉华,李晓丽,赵雪梅,等. 结合马氏距离的区域化模糊聚类遥感图像分割[J]. 中国矿业大学学报,2017,46(1):222-228.

ZHAO Quanhua, LI Xiaoli, ZHAO Xuemei, et al. Remote sensing image segmentation algorithm with regional fuzzy cluster and Mahalanobis distance[J]. Journal of China University of Mining & Technology, 2017, 46(1): 222-228.

[7] YAGER R R, FILEV D P.Approximate clustering via the mountain method[J]. IEEE Transaction on Systems Man & Cybemetics, 2002, 24(8): 1279-1284.

[8] CHIU S L. Fuzzy model identification based on cluster estimation[J]. Journal of Intelligent and Fuzzy System, 1994,2(3):267-278.

[9] 陈东辉. 基于目标函数的模糊聚类算法关键技术研究[D]. 西安：西安电子科技大学,2012.

CHEN Donghui. Research of Key Techniques in Fuzzy Clustering Based on Objective Function[D]. Xi’an：Xidian University, 2012.

[10] 陈新泉.聚类算法中的优化方法应用[M].成都:电子科技大学出版社，2014.

[11] DUNN J C.Well-separated clusters and optimal fuzzy partitions[J].Journal of Cybernetics,1974,4(1):95-104.

[12] BEZDEK J C.Pattern Recognition with Fuzzy Objective Function Algorithms[M].New York:Plenum Press, 1981.

[13] 李洁, 高新波, 焦李成.基于特征加权的模糊聚类新算法[J].电子学报, 2006, 34(1): 89-92.

LI Jie, GAO Xinbo, JIAO Licheng. A new feature weighted fuzzy clustering algorithm[J]. Acta Electronica Sinica, 2006,34(1): 89-92.

[14] 易倩,滕少华,张巍. 基于马氏距离的K均值聚类算法的入侵检测[J]. 江西师范大学学报(自然科学版),2012,36(3):284-287.

YI Qian, TENG Shaohua, ZHANG Wei. Mahalanobis distance-based on K-means clustering algorithm for intrusion detection [J]. Journal of Jiangxi Normal University (Natural Science), 2012, 36(3): 284-287.

[15] LIU J, CHEN H, ZHONG Z, et al. Intrusion detection algorithm for the wormhole attack in Ad Hoc network[C]// Proceedings of International Conference on Computer Science and Information Technology Advances in Intelligent Systems and Computing[S.l.]:[s.n.].2014:147-157.

[16] LIU H C, JENG B C, YlH J M, et al. Fuzzy C-means algorithm based on standard Mahalanobis distances[C]. Proc of the 2009 international Symposium on Information Processing. Huangshan:[s.n.], 2009:422-427.

[17] 王振丽. 基于加权MP马氏距离的GS方法研究[D]. 南京：南京理工大学,2016.

WANG Zhenli. Research about GS Method Based on the Weighted MP Mahalanobis Distance[D]. Nanjing：Nanjing University of Science & Technology, 2016.

[18] 吴香华,牛生杰,吴诚鸥,等. 马氏距离聚类分析中协方差矩阵估算的改进[J]. 数理统计与管理,2011,30(2):240-245.

WU Xianghua, NIU Shengjie, WU Chengou, et al. An improvement on estimating covariance matrix during clusteranalysis using mahalanobis distance[J]. Journal of Applied Statistics and Management, 2011, 30(2): 240-245.

[19] 赵小强,李雄伟. 基于改进马氏距离的模糊C聚类研究[J]. 中南大学学报(自然科学版),2013,44(sup2):195-198.

ZHAO Xiaoqiang, LI Xiongwei. A fuzzy C-means clustering algorithm based on improved Mahalanobis distance[J]. Journal of Central South University (Science and Technology), 2013, 44(sup2): 195-198.

[20] 李静. 模糊聚类算法的研究及应用[D]. 无锡：江南大学,2014.

LI Jing. Research and Application of Fuzzy Clustering Algorithm[D].Wuxi：Jiangnan University, 2014.

[21] NAIK A，SATAPATHY S C，PARVATH K. Improvement of initial cluster center of C-means using teaching learning based optimization[J]. Procedia Technology, 2012,6(4)：428-435.

[22] CELEBI M E, KINGRAVI H A,VELA P A. A comparative study of efficient initialization methods for the k-means clustering algorithm[J]. Expert Systems with Applications,2013, 40(1): 200-210.

[23] STETCO A, ZENG Xiaojun, KEANE J. Fuzzy C-means++: Fuzzy C-means with effective seeding initialization[J]. Expert Systems with Applications, 2015，42(21)：7541-7548.

[24] 汪西莉, 焦李成.一种基于马氏距离的支持向量快速提取算法[J].西安电子科技大学学报(自然科学版)，2004, 31(4):639-643.

WANG Xili, JIAO Licheng. A fast algorithm for extracting the support vector on the Mahalanobis distance[J]. Journal of Xidian University, 2004, 31(4): 639-643.

[25] RUIZ A, LPEZDETERUEL P E. Nonlinear kernel-based statistical pattern analysis[J]. IEEE Trans Neural Netw, 2001,12 (1):16-32.