function?[hisfbestvalbestmemberfbestfmean]?=?BFA（FUNNDimPopSizeMaxIter）

%%?參數(shù)初始化

nc?=?4;?????????????????%?最大連續(xù)前進(jìn)步數(shù)

%%?環(huán)境參數(shù)初始化
his?=?zeros（1MaxIter）;?????????%?全局最優(yōu)點(diǎn)
his_idx?=?1;%全局最優(yōu)點(diǎn)出現(xiàn)的代數(shù)

%%?定義控制變量的上下界
Bound?=?FUN（）;??????????
Boundl=ones（NDimPopSize）;
Boundu=ones（NDimPopSize）;
for?i=1:1:PopSize
????Boundl（:i）=Bound（:1）;
????Boundu（:i）=Bound（:2）;
end
LowerBound?=?Boundl.*ones（NDimPopSize）;????%?每個(gè)細(xì)菌控制變量的下界
UpperBound?=?Boundu.*ones（NDimPopSize）;????%?每個(gè)細(xì)菌控制變量的上界

%%?初始化轉(zhuǎn)向角
angle?=?pi/4.*ones（NDim-1PopSize）;?????%?每個(gè)細(xì)菌的轉(zhuǎn)向角
a?=?round（（（NDim+1）^.5））;???????????????%?根據(jù)控制變量個(gè)數(shù)確定的細(xì)菌的最大轉(zhuǎn)向角

%%?初始化細(xì)菌群體
population?=??rand（NDim?PopSize）.*（UpperBound-LowerBound）?+?LowerBound;?????%?隨機(jī)初始化每個(gè)細(xì)菌的位置
vmax?=?ones（NDimPopSize）;??????????????%?初始化菌群的速度向量
for?i?=?1:NDim?????????????????????????
????vmax（i:）?=?（UpperBound（i:）-LowerBound（i:））;?%?最大的速度向量是變量的范圍
end
r?=?norm（vmax（:1））;????????????????????%?r?是最大速度向量的模
distance?=?r*repmat（ones（1PopSize）NDim1）;????%?定義菌群中每個(gè)個(gè)體前進(jìn)的步長(zhǎng)
direction?=?zeros（NDimPopSize）;????????????????%?定義細(xì)菌轉(zhuǎn)向的方向向量direction

%%?評(píng)價(jià)初始種群
fvalue?=?zeros（1PopSize）;??????????????%?初始化每個(gè)細(xì)菌個(gè)體的適應(yīng)值
for?i?=?1:PopSize???????????????????????%?評(píng)價(jià)種群的中每個(gè)細(xì)菌個(gè)體
????fvalue（i）?=?FUN（population（:i））;
end
[fbestvalindex]?=?min（fvalue）;?????????%記錄群體中最小的適應(yīng)值及其編號(hào)????
bestmember?=?population（:index）;???????%?記錄群體中最好的個(gè)體所在的位置

%%?尋優(yōu)主循環(huán)
iteration?=?0;
while?iteration????iteration?=?iteration?+1;???????????
????
????for?j?=?1:PopSize???????????????????
????????R1?=?randn（1）;??????????????????%?產(chǎn)生步長(zhǎng)公式中的隨機(jī)數(shù)
????????R2?=?rand（NDim-11）;????????????%?隨機(jī)轉(zhuǎn)向角公式中的隨機(jī)數(shù)矩陣
????????R3?=?rand（NDim?1）;?????????????%?密度感應(yīng)中吸引表達(dá)式中的隨機(jī)數(shù)
????????
???????%%?確定轉(zhuǎn)向角
????????Tumbleangle?=?-（pi/（a^2））.*R2+angle（:j）;???%?隨機(jī)選擇轉(zhuǎn)向角
????????distance（:j）?=?r*R1;???????????????????????%?隨機(jī)選擇前進(jìn)步長(zhǎng)
????????direction（1j）?=?prod（cos（Tumbleangle（1:NDim-1）））;??????%?將轉(zhuǎn)向角轉(zhuǎn)換為笛卡爾坐標(biāo)系中的轉(zhuǎn)向向量
????????for?i?=?2:NDim-1
????????????direction（ij）?=?sin（Tumbleangle（i））.*prod（cos（Tumbleangle（i:NDim-1）））;
????????end
????????direction（NDimj）?=?sin（Tumbleangle（NDim-1））;
????????NewPosition?=?population（:j）+distance（:j）.*direction（:j）;????????????%?轉(zhuǎn)向運(yùn)動(dòng)后細(xì)菌的新位置
????????Outflag?=?（NewPosition>UpperBound（:j）?|?NewPosition????????NewPosition?=?NewPosition-Outflag.*distance（:j）.*direction（:j）;???????%?將變量控制在要求范圍之內(nèi)
????????Newfvalue?=?FUN（NewPosition）;???%計(jì)算新位置的適應(yīng)值

?????%%?前進(jìn)運(yùn)動(dòng)
????????searchcounter?=?0;??????????????????????????%?前進(jìn)步數(shù)計(jì)數(shù)器
????????while?searchcounter?????????????if?Newfvalue?????????????????population（:j）?=?NewPosition;??????%?新位置若更優(yōu)秀的話細(xì)菌前進(jìn)到新位置
????????????????angle（:j）?=?Tumbleangle;???????????%?細(xì)菌的轉(zhuǎn)向角更新為當(dāng)前的轉(zhuǎn)向角
????????????????NewPosition?=?population（:j）+distance（:j）.*direction（:j）;????????????????%?更新下一步的位置
????????????????Outflag?=?（NewPosition>UpperBound（:j）?|?NewPosition