function [BESTX,BESTY,ALLX,ALLY]=gafcm(K,N,Pm,LB,UB,D,c,m)
%% ´Ëº¯ÊýʵÏÖÒÅ´«Ëã·¨£¬ÓÃÓÚÄ£ºýC£­¾ùÖµ¾ÛÀà
% GreenSimÍŶÓÔ­´´×÷Æ·£¬×ªÔØÇë×¢Ã÷
% »¶Ó­·ÃÎÊGreenSim¡ª¡ªËã·¨·ÂÕæÍŶӡú[url]http://blog.sina.com.cn/greensim[/url]
% ¡¾×ªÔØ×Ô¡¿[url]http://wenku.baidu.com/view/12c5dcee6294dd88d0d26b1e.html[/url]
%% ÊäÈë²ÎÊýÁбí
% K µü´ú´ÎÊý
% N ÖÖȺ¹æÄ££¬ÒªÇóÊÇżÊý
% Pm ±äÒì¸ÅÂÊ
% LB ¾ö²ß±äÁ¿µÄϽ磬M¡Á1µÄÏòÁ¿
% UB ¾ö²ß±äÁ¿µÄÉϽ磬M¡Á1µÄÏòÁ¿
% D ԭʼÑù±¾Êý¾Ý£¬n¡ÁpµÄ¾ØÕó
% c ·ÖÀà¸öÊý
% m Ä£ºýC¾ùÖµ¾ÛÀàÊýѧģÐÍÖеÄÖ¸Êý
%% Êä³ö²ÎÊýÁбí
% BESTX K¡Á1ϸ°û½á¹¹£¬Ã¿Ò»¸öÔªËØÊÇM¡Á1ÏòÁ¿£¬¼Ç¼ÿһ´úµÄ×îÓŸöÌå
% BESTY K¡Á1¾ØÕ󣬼Ǽÿһ´úµÄ×îÓŸöÌåµÄÆÀ¼Ûº¯ÊýÖµ
% ALLX K¡Á1ϸ°û½á¹¹£¬Ã¿Ò»¸öÔªËØÊÇM¡ÁN¾ØÕ󣬼Ǽȫ²¿¸öÌå
% ALLY K¡ÁN¾ØÕ󣬼Ǽȫ²¿¸öÌåµÄÆÀ¼Ûº¯ÊýÖµ
%% µÚÒ»²½£º
M=length(LB);%¾ö²ß±äÁ¿µÄ¸öÊý
%ÖÖȺ³õʼ»¯£¬Ã¿Ò»ÁÐÊÇÒ»¸öÑù±¾
farm=zeros(M,N);
for i=1:M
x=unifrnd(LB(i),UB(i),1,N);
farm(i,:)=x;
end
%Êä³ö±äÁ¿³õʼ»¯
ALLX=cell(K,1);%ϸ°û½á¹¹£¬Ã¿Ò»¸öÔªËØÊÇM¡ÁN¾ØÕ󣬼Ǽÿһ´úµÄ¸öÌå
ALLY=zeros(K,N);%K¡ÁN¾ØÕ󣬼Ǽÿһ´úÆÀ¼Ûº¯ÊýÖµ
BESTX=cell(K,1);%ϸ°û½á¹¹£¬Ã¿Ò»¸öÔªËØÊÇM¡Á1ÏòÁ¿£¬¼Ç¼ÿһ´úµÄ×îÓŸöÌå
BESTY=zeros(K,1);%K¡Á1¾ØÕ󣬼Ǽÿһ´úµÄ×îÓŸöÌåµÄÆÀ¼Ûº¯ÊýÖµ
k=1;%µü´ú¼ÆÊýÆ÷³õʼ»¯
%% µÚ¶þ²½£ºµü´ú¹ý³Ì
while k<=K
%% ÒÔÏÂÊǽ»²æ¹ý³Ì
newfarm=zeros(M,2*N);
Ser=randperm(N);%Á½Á½Ëæ»úÅä¶ÔµÄÅä¶Ô±í
A=farm(:,Ser(1));
B=farm(:,Ser(2));
P0=unidrnd(M-1);
a=[A(1:P0,:);B((P0+1):end,:)];%²úÉú×Ó´úa
b=[B(1:P0,:);A((P0+1):end,:)];%²úÉú×Ó´úb
newfarm(:,2*N-1)=a;%¼ÓÈë×Ó´úÖÖȺ
newfarm(:,2*N)=b;
for i=1:(N-1)
A=farm(:,Ser(i));
B=farm(:,Ser(i+1));
P0=unidrnd(M-1);
a=[A(1:P0,:);B((P0+1):end,:)];
b=[B(1:P0,:);A((P0+1):end,:)];
newfarm(:,2*i-1)=a;
newfarm(:,2*i)=b;
end
FARM=[farm,newfarm];
%% Ñ¡Ôñ¸´ÖÆ
SER=randperm(3*N);
FITNESS=zeros(1,3*N);
fitness=zeros(1,N);
for i=1:(3*N)
Beta=FARM(:,i);
FITNESS(i)=FIT(Beta,D,c,m);
end
for i=1:N
f1=FITNESS(SER(3*i-2));
f2=FITNESS(SER(3*i-1));
f3=FITNESS(SER(3*i));
if f1<=f2&&f1<=f3
farm(:,i)=FARM(:,SER(3*i-2));
fitness(:,i)=FITNESS(:,SER(3*i-2));
elseif f2<=f1&&f2<=f3
farm(:,i)=FARM(:,SER(3*i-1));
fitness(:,i)=FITNESS(:,SER(3*i-1));
else
farm(:,i)=FARM(:,SER(3*i));
fitness(:,i)=FITNESS(:,SER(3*i));
end
end
%% ¼Ç¼×î¼Ñ¸öÌåºÍÊÕÁ²ÇúÏß
X=farm;
Y=fitness;
ALLX{k}=X;
ALLY(k,:)=Y;
minY=min(Y);
pos=find(Y==minY);
BESTX{k}=X(:,pos(1));
BESTY(k)=minY;
%% ±äÒì
for i=1:N
if Pm>rand&&pos(1)~=i
AA=farm(:,i);
BB=GaussMutation(AA,LB,UB);
farm(:,i)=BB;
end
end
disp(k);
k=k+1;
end
%% »æͼ
BESTY2=BESTY;
BESTX2=BESTX;
for k=1:K
TempY=BESTY(1:k);
minTempY=min(TempY);
posY=find(TempY==minTempY);
BESTY2(k)=minTempY;
BESTX2{k}=BESTX{posY(1)};
end
BESTY=BESTY2;
BESTX=BESTX2;
plot(BESTY,'-ko','MarkerEdgeColor','k','MarkerFaceColor','k','MarkerSize',2)
ylabel('º¯ÊýÖµ')
xlabel('µü´ú´ÎÊý')
grid on