电力系统经济调度-蜜蜂算法

这几天在做一些关于经济调度的问题，有一个很不错的启发式算法叫做蜂群经济调度算法，分享一下！
通过设置蜂群数量，上下界，定义好了损失函数等参数，可以很快达到稳定。
想要直接用的可以去我资源界面下载，与博客同名。

具体代码：

%% Problem 
model=MakeModel();
CostFunction=@(x) Cost(x,model);        % Cost Function
nVar=model.nPlant;                      % Number of Decision Variables
VarSize=[1 nVar];                       % Size of Decision Variables Matrix
VarMin=0;                               % Lower Bound of Variables
VarMax=1;                               % Upper Bound of Variables

%% Bees Algorithm Parameters
MaxIt = 50;                               % Maximum Number of Iterations
nScoutBee = 5;                            % Number of Scout Bees
nSelectedSite = round(0.5*nScoutBee);     % Number of Selected Sites
nEliteSite = round(0.4*nSelectedSite);    % Number of Selected Elite Sites
nSelectedSiteBee = round(0.5*nScoutBee);  % Number of Recruited Bees for Selected Sites
nEliteSiteBee = 2*nSelectedSiteBee;       % Number of Recruited Bees for Elite Sites
r = 0.1*(VarMax-VarMin);	              % Neighborhood Radius
rdamp = 0.95;                             % Neighborhood Radius Damp Rate

%% Start
% Empty Bee Structure
empty_bee.Position = [];
empty_bee.Cost = [];
empty_bee.Sol = [];
% Initialize Bees Array
bee = repmat(empty_bee, nScoutBee, 1);
% Create New Solutions
for i = 1:nScoutBee
bee(i).Position = unifrnd(VarMin, VarMax, VarSize);
[bee(i).Cost bee(i).Sol] = CostFunction(bee(i).Position);
end
% Sort
[~, SortOrder] = sort([bee.Cost]);
bee = bee(SortOrder);
% Update Best Solution Ever Found
BestSol = bee(1);
% Array to Hold Best Cost Values
BestCost = zeros(MaxIt, 1);

%% Bees Algorithm Body
for it = 1:MaxIt
    % Elite Sites
    for i = 1:nEliteSite
        bestnewbee.Cost = inf;
        for j = 1:nEliteSiteBee
            newbee.Position = BeeDance(bee(i).Position, r);
            [newbee.Cost newbee.Sol] = CostFunction(newbee.Position);
            if newbee.Cost<bestnewbee.Cost
                bestnewbee = newbee;
            end
        end
        if bestnewbee.Cost<bee(i).Cost
            bee(i) = bestnewbee;
        end
    end
    % Selected Non-Elite Sites
    for i = nEliteSite+1:nSelectedSite
        bestnewbee.Cost = inf;
        for j = 1:nSelectedSiteBee
            newbee.Position = BeeDance(bee(i).Position, r);
            [newbee.Cost newbee.Sol] = CostFunction(newbee.Position);
            if newbee.Cost<bestnewbee.Cost
                bestnewbee = newbee;
            end
        end
        if bestnewbee.Cost<bee(i).Cost
            bee(i) = bestnewbee;
        end
    end
    % Non-Selected Sites
    for i = nSelectedSite+1:nScoutBee
        bee(i).Position = unifrnd(VarMin, VarMax, VarSize);
        [bee(i).Cost bee(i).Sol] = CostFunction(bee(i).Position);
    end
    % Sort
    [~, SortOrder] = sort([bee.Cost]);
    bee = bee(SortOrder);
    % Update Best Solution Ever Found
    BestSol = bee(1);
    % Store Best Cost Ever Found
    BestCost(it) = BestSol.Cost;
    % Display Iteration Information
    disp(['In Iteration No ' num2str(it) ': Bees Cost is = ' num2str(BestCost(it))]);
    % Damp Neighborhood Radius
    r = r*rdamp;
end

% Plot ITR
plot(BestCost,'-or', 'LineWidth', 1);
xlabel('Bees Iteration');
ylabel('Best Cost Value');
ax = gca; 
ax.FontSize = 12; 
set(gca,'Color','y')
legend({'Bees ED'},'FontSize',12,'FontWeight','bold','TextColor','r');
%% Results
BestSol
%
BestSol.Sol
%
Res=BestSol.Sol.PTotal-BestSol.Sol.PL-model.PD

编写不易，求个点赞！！！！！！！
“你是谁？”

“一个看帖子的人。”

“看帖子不点赞啊？”

“你点赞吗？”

“当然点了。”

“我也会点。”

“谁会把经验写在帖子里。”

“写在帖子里的那能叫经验贴？”

“上流！”
cheer！！！