clear?
clc
NI=30;%試驗(yàn)次數(shù)
t=zeros（1NI）;
for?ni=1:NI
????tic;
%人工蜂群算法
%算法參數(shù)
%采蜜蜂個(gè)數(shù)?CN??引領(lǐng)蜂
%蜜源??MN??
%工蜂數(shù)目?GN??跟隨蜂
N=10;
CN=N;
MN=N;
GN=N;
Limit=50;
%問(wèn)題參數(shù)
%問(wèn)題維數(shù)
D=2;
%變量范圍
xmax=550*ones（1D）;
xmin=50*ones（1D）;
%最大迭代次數(shù)
K=1000;
limit=50;
xcn=zeros（ND）;
xgn=zeros（ND）;
xmn=zeros（ND）;
popu=zeros（2*ND）;
fpopu=zeros（12*N）;

%算法初始化
for?i=1:2*N
????popu（i:）=xmin+rand（1D）.*（xmax-xmin）;
end
%適應(yīng)值計(jì)算
for?i=1:2*N
????fpopu（i）=func（popu（i:））;
end
[Bd]=sort（fpopu）;
i=1;
while?i<=N
????xmn（i:）=popu（d（i）:）;
?????fmn（i）=B（i）;
????i=i+1;
end
???
????xcn=xmn;
????fcn=fmn;
????
%迭代開(kāi)始
k=1;
while?k<=K
????g=k;
????%引領(lǐng)蜂探索
????for?i=1:CN
????????j=ceil（rand*D）;
????????d=ceil（N*rand）;
????????xcn（ij）=xmn（ij）+（2*rand-1）*（xmn（ij）-xmn（dj））;
????end
????for?i=1:CN
??????????fcn（i）=func（xcn（i:））;
????????if?fcn（i）???????????xmn（i:）=xcn（i:）;
???????????fmn（i）=fcn（i）;
????????end
????end
????%跟隨蜂跟隨引領(lǐng)蜂
????xgn=xcn;
????fgn=fcn;
????sf=sum（fgn）;%目標(biāo)函數(shù)的總和