细菌觅食算法求解经济调度问题

function BFO（）
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%*********************细菌觅食算法**********************
%%%%%%%%%%%%%%%%%%%-----BFA算法-----%%%%%%%%%%%%%%%%%%%
    clear;
    clc;
    %-----（1）初始化参数-----
    %机组出力上下限
    pmin=[40  60  80 24 26   68 110 135 135 130    94  94 125 125 125   125 125 220 220 242   242 254 254 254 254   254 254  10  10  10   20 20 20 20 18   18 20 25 25 25 ];
    pmax=[80 120 190 42 42  140 300 300 300 300   375 375 500 500 500   500 500 500 500 550   550 550 550 550 550   550 550 150 150 150   70 70 70 70 60   60 60 60 60 60];
    
    p = 40;    % 搜索范围的维度
    s = 50;   % 细菌的个数
    Nc = 60;  % 趋化的次数
    Ns = 8;   % 趋化操作中单向运动的最大步数
    Nre = 6;    % 复制操作步骤数
    Ned = 6;    % 驱散（迁移）操作数
    Sr = s/2;   % 每代复制（分裂）数
    Ped = 0.25; %