💥 💥 💞 💞 欢迎来到本博客 ❤️ ❤️ 💥 💥
🏆 博主优势: 🌞 🌞 🌞博客内容尽量做到思维缜密,逻辑清晰,为了方便读者。
⛳ 座右铭:行百里者,半于九十。
📋 📋 📋 本文目录如下: 🎁 🎁 🎁
目录
💥1 概述
📚2 运行结果
🎉3 参考文献
🌈4 Python代码实现
💥1 概述
文献来源:
摘要:柔性作业车间调度问题(Flexible Job-shop Scheduling Problem,FJSP)是经典作业车间调度问题的一个扩展,前者更接近于实际生产。以最小化最大完工时间为目标,提出了一种改进的离散粒子群优化算法。传统粒子群优化算法一般适用于优化连续模型问题,FJSP作为复杂度比较高的组合优化问题,是一种典型的离散模型。提出的算法采用机器负荷平衡机制初始化粒子种群,在粒子的更新过程中引入了3个操作算子来更新粒子的工序排序部分和机器分配部分,这3个算子分别为基于工序排序或机器分配的变异、与个体最优位置之间进行工序先后顺序保留的交叉(POX)操作、与全局最优位置进行随机点保存的交叉(RPX)操作。先后执行以上3个算子以完成粒子的一次更新。这种操作能够使种群较快地收敛于最优解。对标准测试案例进行实验的结果表明,所提算法对解决FJSP具有有效性,并且能够快速地搜索到近似最优解;与其他同类算法相比,所提算法在求解效果和收敛速度上均具有优越性。
关键词:
作业车间调度;离散优化问题;柔性;粒子群优化;
📚2 运行结果
#总共15个Brandimarte文件 for i in range(15): #每个数据用例都测试10次,取最好的一次结果,如果为了效率可以每个数据都测试1次 results = [solve_FJSP(i,j) for j in range(10)] Pg_list = [result[0] for result in results ] fitness_list = [result[1] for result in results ] job_op_num = results[0][2] p_table = results[0][3] best_fitness_index = np.argmax(np.array(fitness_list)) best_fitness = fitness_list[best_fitness_index] best_Pg = Pg_list[best_fitness_index] #画图,写入.txt文档 path= './BestFitness/BrandimarteMk'+str(i+1)+'/' Decode.decode(best_Pg,job_op_num,p_table,'save',path) print(best_Pg,best_fitness) with open(path+'best_schedule.txt', 'w') as f: f.write(str(best_Pg)+'\n'+str(best_fitness))
部分代码:
# 生成初始种群 # 种群大小,可以根据m和n的值来调整大小,如C*m*n c为一个常系数 # Popsize = 5*p_table.shape[1]*len(job_op_num) Popsize = 200 encode = Encode(Popsize, p_table, job_op_num) # 全局选择的染色体 global_chrs = encode.global_selection() # #局部选择的染色体 local_chrs = encode.local_selection() # #随机选择的染色体 random_chrs = encode.random_selection() # 合并三者,得到初始的种群 chrs = np.vstack((global_chrs, local_chrs, random_chrs)) # 以下是关于操作染色体的代码 # 初始的超参数赋值 o_mega = 0.15 c1 = 0.5 c2 = 0.7 pf_max = 0.8 pf_min = 0.2 # 迭代次数,也可以根据m和n的值来调整大小, # Iter = 5*p_table.shape[1]*len(job_op_num) Iter = 200 # 得到初始的个体最优位置 P = copy.deepcopy(chrs) # 得到初始的全局最优位置 # Decode.decode(chr,job_op_num,p_table,'decode'),其中的‘decode’表示不画图,只是计算适应度 fitness_list = [Decode.decode(chr, job_op_num, p_table, 'decode',None) for chr in P] Pg = P[np.argmin(fitness_list)] for iter in range(Iter): # 计算pf pf = pf_max - (pf_max - pf_min) / Iter * iter # 更新种群中所有的染色体 copy_chrs = copy.deepcopy(chrs) chrs = [pso.f_operator(job_op_num, p_table, chr, P[index], Pg, pf, o_mega, c1, c2) for index, chr in enumerate(copy_chrs)] # 更新个体最优位置 P = np.array([chr1 if Decode.decode(chr1, job_op_num, p_table, 'decode',None) <= Decode.decode(chr2, job_op_num, p_table, 'decode',None) else chr2 for chr1, chr2 in zip(P, chrs)]) # 更新全局最优位置 fitness_list = [Decode.decode(chr, job_op_num, p_table, 'decode',None) for chr in P] Pg = P[np.argmin(fitness_list)] # for chr in chrs: # print(Decode.decode(chr, job_op_num, p_table, 'decode',None)) # print("第" + str(iter + 1) + '次循环的最优fitness:', Decode.decode(Pg, job_op_num, p_table, 'decode',None)) print("第"+str(file_num+1)+'个数据集,第'+str(run_times+1)+'次运行'+'迭代:'+str(iter+1)+'/'+str(Iter)) fitness = Decode.decode(Pg, job_op_num, p_table, 'decode',None)
🎉3 参考文献
部分理论来源于网络,如有侵权请联系删除。
[1]丁舒阳,黎冰,侍洪波.基于改进的离散PSO算法的FJSP的研究[J].计算机科学,2018,45(04):233-239+256.