车间调度基于matlab NSGA-2算法求解多目标车间调度问题【含Matlab源码 893期】.zip

clc; clear; close all; %% Problem Definition load CastingData Jm T JmNumber DeliveryTime IntervalTime CostFunction=@(x,Jm ,T ,JmNumber ,DeliveryTime, IntervalTime) MyCost(x,Jm ,T ,JmNumber ,DeliveryTime, IntervalTime); nVar=3; VarSize=[1 nVar]; VarMin=-4; VarMax= 4; pfmax=0.9; pfmin=0.2; VarRange=[VarMin VarMax]; %% NSGA-II Parameters MaxIt=500; nPop=50; pCrossover=0.8; nCrossover=round(pCrossover*nPop/2)*2; pMutation=0.3; nMutation=round(pMutation*nPop); mu=0.3; %% Initialization tic; % PNumber 铸件个数 MNumber 工序个数数组 每个工件对应的工序数量有可能不同 PNumber=size(Jm,1); trace=zeros(2, MaxIt); %寻优结果的初始值 MNumber=[]; for i=1:size(Jm,1) sumTemp=0; for j=1:size(Jm,2) if(length(Jm{i,j}))>0 sumTemp=sumTemp+1; end end MNumber=[MNumber,sumTemp]; end WNumber=sum(MNumber); %工序总个数 %% 初始化 Number=MNumber; D=WNumber*2; %粒子群维度 empty_individual.Position=[]; empty_individual.Cost=[]; empty_individual.Rank=[]; empty_individual.CrowdingDistance=[]; empty_individual.DominatedCount=[]; empty_individual.DominationSet=[]; % 初始化种群 pop=repmat(empty_individual,nPop,1); for i=1:nPop WPNumberTemp=Number; if i<nPop/2 for j=1:WNumber %随机产成工序 val=unidrnd(PNumber); while WPNumberTemp(val)==0 val=unidrnd(PNumber); end %第一层代码表示工序 pop(i).Position(j)=val; %随机初始化位置 WPNumberTemp(val)=WPNumberTemp(val)-1; %第2层代码表示机器 TempT=T{val,MNumber(val)-WPNumberTemp(val)}; % 机器加工时间最少初始化 %[~,minTimeIndex]=min(TempT); % 随机机器初始化 mindex=unidrnd(length(TempT)); %随机产成工序机器 pop(i).Position(j+WNumber)=mindex; end else for j=1:WNumber %随机产成工序 val=unidrnd(PNumber); while WPNumberTemp(val)==0 val=unidrnd(PNumber); end %第一层代码表示工序 pop(i).Position(j)=val; %随机初始化位置 WPNumberTemp(val)=WPNumberTemp(val)-1; %第2层代码表示机器 TempT=T{val,MNumber(val)-WPNumberTemp(val)}; % 机器加工时间最少初始化 [~,minTimeIndex]=min(TempT); % 随机机器初始化 %mindex=unidrnd(length(TempT)); %随机产成工序机器 pop(i).Position(j+WNumber)=minTimeIndex; end end end for i=1:nPop pop(i).Cost=CostFunction(pop(i).Position,Jm ,T ,JmNumber ,DeliveryTime, IntervalTime); end % Non-dominated Sorting [pop ,F]=NonDominatedSorting(pop); % Calculate Crowding Distances pop=CalcCrowdingDistance(pop,F); %% NSGA-II Loop for it=1:MaxIt % Crossover popc=repmat(empty_individual,nCrossover,1); pf=pfmax-(pfmax-pfmin)*it/MaxIt; for k=1:nCrossover i1=BinaryTournamentSelection(pop); i2=BinaryTournamentSelection(pop); % [popc(k,1).Position, popc(k,2).Position]=Crossover(pop(i1).Position,pop(i2).Position,VarRange); popc(k,1).Position= CrossParticle(pop(i1).Position,pop(i2).Position,Jm,pf); popc(k,1).Cost=CostFunction(popc(k,1).Position,Jm ,T ,JmNumber ,DeliveryTime, IntervalTime); end popc=popc(:); % Mutation popm=repmat(empty_individual,nMutation,1); for k=1:nMutation i=BinaryTournamentSelection(pop); if rand()<mu popm(k).Position=Swap(pop(i).Position,Jm); popm(k).Cost=CostFunction(popm(k).Position,Jm ,T ,JmNumber ,DeliveryTime, IntervalTime); else popm(k).Position=pop(i).Position; popm(k).Cost=pop(i).Cost; end end % Merge Pops pop=[pop popc popm]; % Non-dominated Sorting [pop, F]=NonDominatedSorting(pop); % Calculate Crowding Distances pop=CalcCrowdingDistance(pop,F); % Sort Population pop=SortPopulation(pop); % Delete Extra Individuals pop=pop(1:nPop); % Non-dominated Sorting [pop, F]=NonDominatedSorting(pop); % Calculate Crowding Distances pop=CalcCrowdingDistance(pop,F); % Plot F1 PF=pop(F{1}); PFCosts=[PF.Cost]; popCosts=[pop.Cost]; firstObj=popCosts(1,:); secondObj=popCosts(2,:); trace(1, it)=min(firstObj); trace(2, it)=min(secondObj); % 画图 fig=figure(1); set(fig,'NAME','NSGA-MultiObj'); plot(PFCosts(1,:),PFCosts(2,:),'ro'); xlabel('间隔时间拖时'); ylabel('交货延期'); % Show Iteration Information disp(['Iteraion ' num2str(it) ': Number of F1 Members = ' num2str(numel(PF))]); end figure('NAME','NSGA-trace') plot(trace(1,:)); hold on; plot(trace(2,:),'-.');grid; legend('间隔时间拖时','交货延期'); [crowdistances]=cat(1,PF.CrowdingDistance); [~,maxIndex]=max(crowdistances); solutions=cat(1,PF.Position); bestSolution=solutions(maxIndex,:); %根据基因，计算调度工序 P=calP(bestSolution,PNumber); %根据调度工序，计算出调度工序时间 PVal=caltime(bestSolution,P,JmNumber,T,Jm); %% 显示最优解 figure('NAME','NSGA-Res'); MP=bestSolution(1,length(bestSolution)/2+1:length(bestSolution)); for i=1:WNumber val= P(1,i); a=(mod(val,100)); %工序 b=((val-a)/100); %工件 Temp=Jm{b,a}; mText=Temp(MP(1,i)); x1=PVal(1,i); x2=PVal(2,i); y1=mText-1; y2=mText; PlotRec(x1,x2,mText); PlotRec(PVal(1,i),PVal(2,i),mText); hold on; fill([x1,x2,x2,x1],[y1,y1,y2,y2],[1-1/b,1/b,b/PNumber]); text(x1,mText-0.25,num2str(P(i))); end RunTime=toc; %% Results disp(' '); disp(['Total Run Time: ' num2str(RunTime) ' second(o)']); disp(' ');