Bacteria Foraging Optimization 多头菌算法

% Bacterial Foraging Algorithm
clear all;
clc;
% parameters initializing
ExtraSpeed=0.1;
Bounds=[-5.125.12;-5.125.12];%boundaries of variable
Dimension=2;% the number of every bacteria features
Number=26;% bacteria number
Iteration=100; %the iteration
MaxMoving=4;%the maximum length of moving in  one direction
C（:1）=0.001*ones（Number1）;%after reverse direction the length of moving forward
NumberReplication=6;
NumberMigration=3;
NumberDivision=Number/2;
ProbabilityMigration=0.25+ExtraSpeed+ExtraSpeed;
DoseAttract=0.05+ExtraSpeed;
SpeedToReleaseDoseAttract=0.05+ExtraSpeed;
DoseRepellant=0.05+ExtraSpeed;
SpeedToReleaseDoseRepellant=0.05+ExtraSpeed;
% initializing bacteria position
for i=1:Number
    P（1i111） = -5.12 + rand*10.24;
    P（2i111） = -5.12 + rand*10.24; 
end
% now for iteration
for l=1:NumberMigration
    for k=1:NumberReplication
        for j=1:Iteration
            for i=1:Number
                %J means fitness
                J（ijkl） = Cost（P（:ijkl））; 
                %Jcc is a akin of pheromone
                Jcc=sum（-DoseAttract*exp（-SpeedToReleaseDoseAttract*（（P（1ijkl）-...
                    P（11:Numberjkl））.^2+（P（2ijkl）-P（21:Numberjkl））.^2））） +...
                    sum（DoseRepellant*exp（-SpeedToReleaseDoseRepellant*（（P（1ijkl）-... 
                    P（11:Numberjkl））.^2+（P（2ijkl）-P（21:Numberjkl））.^2）））; 
                J（ijkl） = J（ijkl） + Jcc;  
                Jlast = J（ijkl）;
                %Delta mean reverse
                Delta（:i） = （2*round（rand（Dimension1））-1）.*rand（Dimension1）; 
                %PHI mean after reversing choose a random direction move forward
                PHI = Delta（:i）/sqrt（Delta（:i）‘*Delta（:i））;  
                P（:ij+1kl） = P（:ijkl） + C（ik）*PHI;
                J（ij+1kl） = Cost（P（:ij+1kl））; 
                %m means roamming
                m =