[发明专利]基于量子粒子群优化改进的模糊C-均值聚类方法

权利要求书

1.一种基于量子粒子群优化改进的模糊C-均值聚类方法，其特征是，所述模糊C-均值聚类方法包括如下步骤：

步骤一、聚类参数设定：给定所需的聚类类别c、粒子群规模N、最大迭代次数MAXTER，当前迭代次数t，收缩扩张系数α；

步骤二、种群的初始化：定义距离标准

d(x_k,v_i)＝1-exp(-β||x_k-v_i||²)                            （1）

其中，β=[1n(Σk=1n||xk-x‾||)]-1,]]>且x‾=1n(Σk=1nxk);]]>

利用单次AFCM算法进行快速分类，结合公式（1）得到目标函数

JAFCM(X,U,V)=Σi=1cΣk=1nUikm[1-exp(-β||xk-vi||2)]---(2)]]>

其中，X＝{x₁,x₂,…x_k,…x_n}是给定的数据集，n为数据集中样本的个数，x_k(x_k1,x_k2，…,x_kd)表示具有d个属性的数据对象，聚类个数为c，聚类中心集合V＝{v₁,v₂,…,v_c}，其中v_i(v_i1，v_i2，…,v_id)表示第i个聚类中心，U_ik表示第k个样本属于第i类的模糊隶属度，对于任意k,i和模糊加权指数m,使公式（2）达到最小值的必要条件为：

Uik=1Σj=1c{[1-exp(-β||xk-vi||2)]/[1-exp(-β||xk-vj||2)]}1/m-1,1≤i≤c,1≤k≤n---(3)]]>

中心向量为：

vi=Σk=1nUikm(exp(-β||xk-vi||2))xkΣk=1nUikm(exp(-β||xk-vi||2))---(4)]]>

利用公式（4）得到的聚类中心向量作为种群中的一个初始粒子，并重复进行N次，生成初始种群；

步骤三、粒子的全局搜索：在一个d维的目标搜索空间中，QPSO算法由N个代表潜在问题解的粒子组成群体X(t)＝{X₁(t),X₂(t),…,X_N(t)}，在t时刻，第i个粒子位置为

粒子没有速度向量，个体最优位置Pbest表示为p_i(t)＝(p_i1(t),p_i2(t),…,p_id(t))，全局最优位置Gbest表示为G(t)＝(G₁(t),G₂(t),…,G_j(t)…,G_d(t))，且G(t)＝p_g(t)，其中g为粒子群体处于全局最优位置的下标，g∈{1,2,…,N}；定义所有粒子个体的最优平均位置C(t)，即

C(t)=(C1(t),C2(t),...Cj(t),...Cd(t))]]>

=1NΣi=1Npi(t)=(1NΣi=1Npi1(t),1NΣi=1Npi2(t),...,1NΣi=1Npij(t),...,1NΣi=1Npid(t))---(6)]]>

粒子的个体最优位置Pbest由下述公式确定为

pi(t)=Xi(t)f[Xi(t)]<f[pi(t-1)]pi(t-1)f[Xi(t)]≥f[pi(t-1)]---(7)]]>

其中，f[X]＝J_AFCM(X，U，V)；

整个粒子群的全局最优位置Gbest由下述公式确定为

g=argMin1<i≤N{f[pi(t)]}---(8)]]>

G(t)＝P_g(t)                                                    (9)

更新粒子时，首先得到粒子的个体最优位置Pbest及整个粒子群的全局最优位置Gbest，其次利用公式(2)更新所有粒子个体的适应度值，并将每一个粒子的当前个体位置适应度值与全局最优位置的适应度值比较；若每一个粒子的当前个体位置适应度值优于全局最优位置的适应度值，更新公式(9)中的G(t)；

令

则粒子的进化方程为：

Xi,j(t+1)=Pi,j(t)±α*|Cj(t)-Xi,j(t)|*ln(1ui,j(t)),ui,j(t)~U(0,1)---(11)]]>

其中，u_i，j(t)为符合正态分布的辅助参数；G_j(t)表示粒子群体的全局最优位置；

步骤四、输出聚类结果：重复上述粒子更新操作，如果当前的迭代次数达到预先设定的最大次数MAXTER，则停止迭代输出最佳聚类效果。

2.根据权利要求1所述的基于量子粒子群优化改进的模糊C一均值聚类方法，其特征是：所述收缩扩张系数α为α＝(1.0-0.5)*(MAXITER-T)/MAXITER+0.5。