function?[GlobalParamsGlobalMinBestGlobalMins]=abc（objfundmax_itersizepopublb）

%%?參數(shù)初始化
NP=sizepop;?%?蜂群規(guī)模
FoodNumber=NP/2;?%?蜜源（解）數(shù)量
limit=100;?%?當(dāng)有蜜源連續(xù)沒(méi)被更新的次數(shù)超過(guò)limit時(shí)，該蜜源將被重新初始化
maxCycle=max_iter;?%?最大迭代次數(shù)
%?待優(yōu)化參數(shù)信息
D=d;?%?待優(yōu)化參數(shù)個(gè)數(shù)，次數(shù)為c和g兩個(gè)
ub=ones（1D）*ub;?%?參數(shù)取值上界
lb=ones（1D）*lb;?%?參數(shù)取值下界
runtime=1;?%?可用于設(shè)置多次運(yùn)行（讓ABC算法運(yùn)行runtime次）以考察程序的穩(wěn)健性

BestGlobalMins=[];?%用于存放全局最優(yōu)解與平均適應(yīng)度值
BestGlobalParams=zeros（runtimeD）;?%?用于存放ABC算法優(yōu)化得到的最優(yōu)參數(shù)
GlobalMin=inf;
GlobalParams=[inf?inf];
%%?初始化蜜源
Range?=?repmat（（ub-lb）[FoodNumber?1]）;
Lower?=?repmat（lb?[FoodNumber?1]）;
Foods?=?rand（FoodNumberD）?.*?Range?+?Lower;

%?計(jì)算每個(gè)蜜源（解）得目標(biāo)函數(shù)值，objfun為計(jì)算SVM預(yù)測(cè)的平均平方誤差（MSE）的函數(shù)，根據(jù)自己的實(shí)際問(wèn)題變異目標(biāo)函數(shù)即可
ObjVal=ones（1FoodNumber）;
for?k?=?1:FoodNumber
????ObjVal（k）?=?objfun（Foods（k:））;
end
Fitness=calculateFitness（ObjVal）;?%?計(jì)算適應(yīng)度函數(shù)值

%?標(biāo)記初始最優(yōu)蜜源（解）
BestInd=find（ObjVal==min（ObjVal））;
BestInd=BestInd（end）;
if?（ObjVal（BestInd）????GlobalMin=ObjVal（BestInd）;
????GlobalParams=Foods（BestInd:）;
end
%%
for?r=1:runtime?%?運(yùn)行ABC算法runtime次
????trial=zeros（1FoodNumber）;?%?用于記錄第i個(gè)蜜源有連續(xù)trail（i）次沒(méi)被更新過(guò)
????iter=1;?%?迭代開(kāi)始
????while?（（iter?<=?maxCycle））?%?循環(huán)條件
????????
????????%%%%%%%%%%%%%%%%%%%%%引領(lǐng)蜂搜索解的過(guò)程%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
????????for?i=1:（FoodNumber）?%?遍歷每個(gè)蜜源（解）
????????????Param2Change=fix（rand*D）+1;?%?隨機(jī)選擇需要變異的參數(shù)
????????????neighbour=fix（rand*（FoodNumber））+1;?%?隨機(jī)選擇相鄰蜜源（解）以準(zhǔn)備變異
????????????%?需要保證選擇的相鄰蜜源不是當(dāng)前蜜源（i）
????????????while（neighbour==i）
????????????????neighbour=fix（rand*（FoodNumber））+1;
????????????end
????????????
????????????sol=Foods（i:）;?%?提取當(dāng)前蜜源（解）對(duì)應(yīng)的的參數(shù)
????????????%?參數(shù)變異得到新的蜜源：v_{ij}=x_{ij}+\phi_{ij}*（x_{kj}-x_{ij}）
????????????sol（Param2Change）=Foods（iParam2Change）+（Foods（iParam2Change）-Foods（neighbourParam2Change））*（rand-0.5）*2;
????????????%?確保參數(shù)取值范圍不越界
????????????ind=find（sol????????????sol（ind）=lb（ind）;
????????????ind=find（sol>ub）;
????????????sol（ind）=ub（ind）;
????????????%?計(jì)算變異后蜜源的目標(biāo)函數(shù)值和適應(yīng)度函數(shù)值
????????????ObjValSol=objfun（sol）;
????????????FitnessSol=calculateFitness（ObjValSol）;
????????????%?更新當(dāng)前蜜源的相關(guān)信息
????????????if?（FitnessSol>Fitness（i））
????????????????Foods（i:）=sol;
????????????????Fitness（i）=FitnessSol;
????????????????ObjVal（i）=ObjValSol;
????????????????trial（i）=0;?%?如果當(dāng)前蜜源被更新了，則對(duì)應(yīng)的trial歸零
????????????else
????????????????trial（i）=trial（i）+1;?%?如果當(dāng)前蜜源沒(méi)有被更新，則trial（i）加1
????????????end
????????end
????????
????????%%%%%%%%%%%%%%%%%%%%%%%%?跟隨蜂搜索解的過(guò)程?%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
????????%?計(jì)算解（蜜源）的概率
????????prob=（0.9.*Fitness./max（Fitness））+0.1;
????????%?循環(huán)初始化
????????i=1;
????????t=0;
????????while（t????????????if（rand????????????????t=t+1;?%?循環(huán)計(jì)數(shù)器加1
????????????????
????????????????Param2Change=fix（rand*D）+1;?%?隨機(jī)確定需要變異的參數(shù)
????????????????neighbour=fix（rand*（FoodNumber））+1;?%?隨機(jī)選擇相鄰蜜源（解）
????????????????%?需要保證選擇的相鄰蜜源不是當(dāng)前蜜源（i）
????????????????while（neighbour==i）
???????????????

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

????.CA....??????5844??2019-05-05?20:37??差分蜂群優(yōu)化算法MATLAB實(shí)現(xiàn)\abc.m

????.CA....??????5714??2019-05-05?20:37??差分蜂群優(yōu)化算法MATLAB實(shí)現(xiàn)\ABC1.m

????.CA....???????182??2008-06-09?19:54??差分蜂群優(yōu)化算法MATLAB實(shí)現(xiàn)\calculateFitness.m

????.CA....???????322??2019-03-30?13:23??差分蜂群優(yōu)化算法MATLAB實(shí)現(xiàn)\Code.m

????.CA....??????1613??2019-03-30?13:32??差分蜂群優(yōu)化算法MATLAB實(shí)現(xiàn)\Cross.m

????.CA....??????2590??2019-05-05?19:44??差分蜂群優(yōu)化算法MATLAB實(shí)現(xiàn)\DE.m

????.CA....??????2496??2020-05-19?13:29??差分蜂群優(yōu)化算法MATLAB實(shí)現(xiàn)\de1.m

????.CA....??????1547??2019-05-05?20:29??差分蜂群優(yōu)化算法MATLAB實(shí)現(xiàn)\DEABC.m

????.CA....????????76??2019-04-29?10:54??差分蜂群優(yōu)化算法MATLAB實(shí)現(xiàn)\f1.m

????.CA....???????186??2019-04-29?10:54??差分蜂群優(yōu)化算法MATLAB實(shí)現(xiàn)\f2.m

????.CA....???????174??2019-04-29?10:51??差分蜂群優(yōu)化算法MATLAB實(shí)現(xiàn)\f3.m

????.CA....???????137??2019-04-29?10:54??差分蜂群優(yōu)化算法MATLAB實(shí)現(xiàn)\f4.m

????.CA....??????2462??2019-05-03?13:14??差分蜂群優(yōu)化算法MATLAB實(shí)現(xiàn)\ga.m

????.CA....??????1239??2020-05-19?13:24??差分蜂群優(yōu)化算法MATLAB實(shí)現(xiàn)\MAIN.m

????.CA....??????1455??2019-03-30?13:54??差分蜂群優(yōu)化算法MATLAB實(shí)現(xiàn)\Mutation.m

????.CA....??????2657??2020-05-19?13:31??差分蜂群優(yōu)化算法MATLAB實(shí)現(xiàn)\pso.m

????.CA....???????806??2019-03-30?13:25??差分蜂群優(yōu)化算法MATLAB實(shí)現(xiàn)\select.m

????.CA....???3456266??2019-03-30?16:28??差分蜂群優(yōu)化算法MATLAB實(shí)現(xiàn)\一種新的混合智能極限學(xué)習(xí)機(jī)_林梅金.caj

????.CA....????????55??2020-05-19?13:32??差分蜂群優(yōu)化算法MATLAB實(shí)現(xiàn)\程序說(shuō)明.txt

????.C.D...?????????0??2020-05-19?13:32??差分蜂群優(yōu)化算法MATLAB實(shí)現(xiàn)

-----------?---------??----------?-----??----

??????????????3485821????????????????????20