matlab_流水车间_能耗调度_流水强开_能耗_能耗调度

function [BestGene,aa]=GA(MaxGeneration,GeneSize,GeneNum,pcross,pmute) Parent = Init(GeneSize,GeneNum); [BestGene,Parent] = KeepBest(Parent); aa = []; for i = 1:MaxGeneration [i E(BestGene)]; Child = chose(Parent); Child = cross(Child,pcross); Child = mute(Child,pmute); [BestGene,Parent] = KeepBest(Child); aa = [aa;E(BestGene)]; end end function GeneInit = Init(GeneSize,GeneNum) %GeneInit = [2 1 4 3;1 3 4 2;3 1 2 4;4 3 2 1]; GeneInit=[]; i=0; while(i~=GeneSize) temp=randperm(GeneNum); if ~ismember(temp,GeneInit,'rows') i=i+1; GeneInit=[GeneInit;temp]; end end end function Child = chose(Parent) GeneSize = size(Parent,1); for i = 1:GeneSize x = Parent(i,:); val(i) = 1/E(x); end ValSum = sum(val); val = val / ValSum; for i = 2:GeneSize val(i) = val(i) + val(i-1); end for i = 1:GeneSize randval = rand; if randval <= val(1) Child(i,:) = Parent(1,:); end for j = 2:GeneSize if randval > val(j-1) && randval <= val(j) Child(i,:) = Parent(j,:); break; end end end Child(end,:) = Parent(end,:); end function Child = cross(Parent,pcross) [GeneSize,GeneNum] = size(Parent); GeneSize = GeneSize - 1; Child = Parent; for i = 1:GeneSize/2 if rand < pcross flag = 0; while( flag==0 ) randval1 = floor((GeneNum-1)*rand) + 1; randval2 = floor((GeneNum-1)*rand) + 1; if randval1 ~= randval2 flag = 1; end end temp = Child(2*i-1,randval1:randval2); Child(2*i-1,randval1:randval2) = Child(2*i,randval1:randval2); Child(2*i,randval1:randval2) = temp; [Child(2*i-1,:),Child(2*i,:)]=checkgene(Child(2*i-1,:),Child(2*i,:)); end end end function Child = mute(Parent,pmute) [GeneSize,GeneNum] = size(Parent); GeneSize = GeneSize - 1; Child = Parent; for i = 1:GeneSize if rand < pmute t1=randi(GeneNum); t2=randi(t1); t3=randi([t1 GeneNum]); t4=Child(i,t2); Child(i,t2)=Child(i,t3); Child(i,t3)=t4; end end end function [BestGene,Parent] = KeepBest(Child) [GeneSize,GeneNum] = size(Child); for i = 1:GeneSize val(i) = E(Child(i,:)); end [~ , flag]=min(val); BestGene = Child(flag,:); Parent = Child; Parent(1,:) = BestGene; Parent(end,:) = BestGene; end function [V1, V2]=checkgene(X1, X2) %1 find the repetation of X1 t=X1; t=sort(t); rX1=t(diff(t)==0); t=X2; t=sort(t); rX2=t(diff(t)==0); %2 find the pos. of rX [m,n]=size(rX1); if n==0 V1=X1; V2=X2; return; end for i=1:n t1=find(X1==rX1(i)); pos1(i)=t1(1); t2=find(X2==rX2(i)); pos2(i)=t2(1); end rX1=X1(sort(pos1)); rX2=X2(sort(pos2)); %2 exchange for i=1:n t1=find(X1==rX1(i)); X1(t1(2))=rX2(i); t2=find(X2==rX2(i)); X2(t2(2))=rX1(i); end V1=X1; V2=X2; end