simulation.rar_militaryfqq_optimization_改进的粒子群_约束型求解器_遗传粒子群

function T_temp = f2_single(Params,pop_temp,show_res) N_temp = size(pop_temp,1); M_temp = zeros(1,N_temp); T_temp = []; for i = 1:N_temp temp_C = Params.C; %% 考虑加工工序,如果工序不合格，则temp_fyz = 1，重新生成 temp_fyz = 0; for j = 1:Params.N %对于所有工序 if j == 1 %第一个工序肯定满足时间需求，且不用对比工作环境，直接分配工作站 gzgx = [pop_temp(i,j)]; %记录工作站内的工序 temp_C = temp_C - Params.ZYSJ(pop_temp(i,j)); %计算工作站剩余工作时间 elseif temp_C >= Params.ZYSJ(1,pop_temp(i,j)) &&... Params.HJ(pop_temp(i,j))==Params.HJ(pop_temp(i,j-1)) %工作站还有工作时间时，且与前一个工序工作站的环境不同时 temp_C = temp_C - Params.ZYSJ(pop_temp(i,j)); %计算工作站剩余工作时间 if j == Params.N %如果是最后一道工序 M_temp(i) = M_temp(i) + 1; %该工作站已完成工作 T_temp = [T_temp,Params.C - temp_C]; gzgx = [gzgx,pop_temp(i,j)]; %添加该工序到工作站内的工序记录 if show_res disp(strcat('第',num2str(M_temp(i)),'个工作站内的工序为：',num2str(gzgx))); end else gzgx = [gzgx,pop_temp(i,j)]; %添加该工序到工作站内的工序记录 end else %如果工作站没有剩余工作时间或者工作环境需要变化时 T_temp = [T_temp,Params.C - temp_C]; %记录工作站花费时间 M_temp(i) = M_temp(i) + 1; %开启新的工作站 if show_res disp(strcat('第',num2str(M_temp(i)),'个工作站内的工序为：',num2str(gzgx))); end gzgx = [pop_temp(i,j)]; temp_C = Params.C - Params.ZYSJ(pop_temp(i,j)); %由于原工作站没有剩余工作时间，因此重置新工作站的时间、特殊环境标记 end end %% 先后序判定 for l = 1:Params.N %任意工序l for m = 1:1:Params.N %任意工序m if Params.QX(m,l) == 1 && ~any(pop_temp(i,[1:find(pop_temp(i,:)==l)])==m) %如果工序m是工序l的前序,且m不在工序l之前 temp_fyz = 1; end end end %% 罚函数 if temp_fyz == 1 M_temp(i) = Params.N + M_temp(i); %如果temp_fyz为1，也就是不符合前序时，进行惩罚，即工作站比工序还多 end end if show_res disp(strcat('系统最优工作站数目为：',num2str(M_temp))); end end