replace.m

模糊聚类算法

function r1=replace(holdm,v1,hist,c,population,f)
%%子优代父，基于浓度的群体更新
tempP=rand(1);%生成0-1间的随机数
C=0.5;%抗体浓度
a=0.3;%参数1
b=0.5;%参数2

    father=f;
    for i=1:population
        [V(i,:),f(i)]=fcm(hist,c,v1(i,:));
    end
    child=f;
    Ps=a*C*abs(1-f./max(f))+b*f./max(f);%基于浓度更新的公式
    hold=v1;
    for i=1:population
       if child(i)>father(i) & tempP<Ps
          hold(i)=v1(i);
       else hold(i)=holdm(i);
      end
   end
   r1=hold;