差异演化算法求解二次分配问题.pdfVIP

第 11 卷第 34 期 2011 年 12 月 1671--1815(2011)34-制09-04 科学技术与工程 Science Technology and Engineering Vol. 11 No.34 Dec. 2011 @ 2011 Sci. Tech. Engrg. 差异演化算法求解二次分配问题 杨卿誉王志刚拿 (南京师范大学泰州学院数学科学与应用学院，泰州 2253∞) 摘要 二次分配问题是典型的 NP 难题。建立了二次分配问题的数学模型。设计了基于差异演化算法的新方法对其进行 求解。给出了差异演化算法求解该问题的具体方案。对不同的二次分配问题算例进行了仿真实验。结果表明，算法可以有 效、快速地找到二次分配问题的最优解。 关键词差异演化算法 二次分配问题 组合优化 中回法分类号 0232; 文献标志码 A 二次分配问题 (Quadratic Assignment Problem , QAP)是一个典型的组合优化问题，它在选址规划、 集成电路设计、图像处理、键盘布局和作业调度等 众多领域有广泛的应用。因此，对该问题的研究具 有重要的理论与实用价值。目前，许多专家学者都 对其进行了深入的研究[14] ，为该问题的求解提供 了独特的思路。 差异演化算法(Differential Evolution , DE) 是一 种实数编码的基于种群进化的全局优化算法[16] ， 是 Rainer Storn 和 Kenneth Prince 在 1995 年提出的。 差异演化算法具有收敛速度快、操作简单、易编程 实现、极强的稳健性等优点，在很多领域得取到了 良好的应用效果。但这些应用主要集中在连续优 化领域，在离散优化领域的研究和应用相对较少。 本文尝试将差异演化算法应用于求解二次分配问 题，仿真实验的结果说明该方法求解二次分配问题 是有效的和可行的。 1 二次分配问题的描述 二次分配问题的一般提法为:给定 n 个设施和 n个地点，要求给每个设施分配一个位置，定义两个 nXn 矩阵，即 A = (Jij) 和 B = (dij) ，其中马表示设施 2011 年 9 月 13 日收到 ·通信作者简介:王志刚。 E-rnail: wzg19. 配ut@163.com。 i 和设施j之间的流量 ， dij表示地点 i 和地点j 之间 的距离。目标是找到一个排序 P = (P 1 , P 2 , …, Pn )，使其: n n min C(P) = L LhjdL(i)L(j) (1) 式(1)中 L (i) 表示设施 i 被分配的地点 ， P ε II (n) 。 2 差异演化算法 差异演化算法是一种基于实数编码的演化算 法，它的整体结构类似于遗传算法，与遗传算法的 主要区别在变异操作上，差异演化算法的变异操作 是基于染色体的差异向量进行的，其余操作和遗传 算法类似。下面通过求解非线性函数f( 叫，屿，…， Xn ) 的最小值问题 ， Xj 满足 zf 运Xj ~白了， j=LL·-v n ，来介绍差异演化算法的操作过程: 令 Xi ( t) 是第 t 代的第 i 个染色体则 Xi (t) = (Xi1 (t) ,x i2 (t) ,… , X in ( t) ) , i = 1 , 2 ,… , N;t=1 ,2 , …, tmu。其中 ， n 是染色体的长度，即变量的个数 ，N 为群体规模 ， tmax是最大进化代数。 2.1 生成初始种群 在 n 维空间里随机产生满足约束条件的 N个 染色体，实施如下措施: 坷ij ( 0) = randij ( 0 , 1 )( x~ - xt) + xt , i = 1 ,2 , ... , Nj = 1 ,2 ,… ,n (4) 8410 科学技术与工程 11 卷 式(4 )中元;与元;分别是第 j 个变量上界和下界， randij(O ， I) 是[0 ， 1] 之间的随机小数。 2.2 变异操作 从群体中随机选择 3 个染色体与1 ，如，与且 (i笋pl 7p2 7p3) ，则 hij(t + 1) =xp以 t) + λ(zp勾 (t) - xp3j (t) ) (5 ) Xp勾 (t+l) -x向 (t + 1) 为差异化向量，λ 为缩放 因子。 2.3 交叉操作 交叉操作是为了增加群体的多样性，具体操作 如下: rhü(t + 1) ,randiO ,I) ~pc 或j=rand(i) vü(t+l)=~ :. ? 问(t) ,randij(O ,I) pc 且j于是 rand( i) (6) randij(O , 1 )是在 [0 ， 1] 之间的随机小数， pc 为交叉 概率，严 ε[0 ， 1] ， rand(i) 为在 [1 ， n]