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

%%?參數初始化

nc?=?4;?????????????????%?最大連續前進步數

%%?環境參數初始化
his?=?zeros（1MaxIter）;?????????%?全局最優點
his_idx?=?1;%全局最優點出現的代數

%%?定義控制變量的上下界
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）;????%?每個細菌控制變量的下界
UpperBound?=?Boundu.*ones（NDimPopSize）;????%?每個細菌控制變量的上界

%%?初始化轉向角
angle?=?pi/4.*ones（NDim-1PopSize）;?????%?每個細菌的轉向角
a?=?round（（（NDim+1）^.5））;???????????????%?根據控制變量個數確定的細菌的最大轉向角

%%?初始化細菌群體
population?=??rand（NDim?PopSize）.*（UpperBound-LowerBound）?+?LowerBound;?????%?隨機初始化每個細菌的位置
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）;????%?定義菌群中每個個體前進的步長
direction?=?zeros（NDimPopSize）;????????????????%?定義細菌轉向的方向向量direction

%%?評價初始種群
fvalue?=?zeros（1PopSize）;??????????????%?初始化每個細菌個體的適應值
for?i?=?1:PopSize???????????????????????%?評價種群的中每個細菌個體
????fvalue（i）?=?FUN（population（:i））;
end
[fbestvalindex]?=?min（fvalue）;?????????%記錄群體中最小的適應值及其編號????
bestmember?=?population（:index）;???????%?記錄群體中最好的個體所在的位置

%%?尋優主循環
iteration?=?0;
while?iteration????iteration?=?iteration?+1;???????????
????
????for?j?=?1:PopSize???????????????????
????????R1?=?randn（1）;??????????????????%?產生步長公式中的隨機數
????????R2?=?rand（NDim-11）;????????????%?隨機轉向角公式中的隨機數矩陣
????????R3?=?rand（NDim?1）;?????????????%?密度感應中吸引表達式中的隨機數
????????
???????%%?確定轉向角
????????Tumbleangle?=?-（pi/（a^2））.*R2+angle（:j）;???%?隨機選擇轉向角
????????distance（:j）?=?r*R1;???????????????????????%?隨機選擇前進步長
????????direction（1j）?=?prod（cos（Tumbleangle（1:NDim-1）））;??????%?將轉向角轉換為笛卡爾坐標系中的轉向向量
????????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）;????????????%?轉向運動后細菌的新位置
????????Outflag?=?（NewPosition>UpperBound（:j）?|?NewPosition????????NewPosition?=?NewPosition-Outflag.*distance（:j）.*direction（:j）;???????%?將變量控制在要求范圍之內
????????Newfvalue?=?FUN（NewPosition）;???%計算新位置的適應值

?????%%?前進運動
????????searchcounter?=?0;??????????????????????????%?前進步數計數器
????????while?searchcounter?????????????if?Newfvalue?????????????????population（:j）?=?NewPosition;??????%?新位置若更優秀的話細菌前進到新位置
????????????????angle（:j）?=?Tumbleangle;???????????%?細菌的轉向角更新為當前的轉向角
????????????????NewPosition?=?population（:j）+distance（:j）.*direction（:j）;????????????????%?更新下一步的位置
????????????????Outflag?=?（NewPosition>UpperBound（:j）?|?NewPosition