应用高斯和声搜索的振动筛优化设计方法
【技术领域】
[0001] 本发明涉及振动筛优化设计领域，尤其是涉及一种应用高斯和声搜索的振动筛优 化设计方法。
【背景技术】
[0002] 振动筛是工业加工和制造中的一种常用机械设备，它在工业生产中具有非常广泛 的应用。振动筛优化设计的一般过程是首先根据待优化设计的振动筛的物理机械特性建立 起优化目标函数，然后确定其所需优化设计的变量，并确定振动筛优化设计所需满足的约 束条件。通常，振动筛的优化设计结果是否合理常常会影响振动筛的最终筛分效果。因此， 振动筛的优化设计方法对于振动筛的最终筛分效果是至关重要的。
[0003] 传统的振动筛优化设计方法往往都是利用优化目标函数的数学特性来求解的， 例如梯度下降法、共辄梯度法等。传统的振动筛优化设计方法往往要求优化目标函数需 要满足连续、可导等特性。然而，在工业生产中实际所需的振动筛是多种多样，形式各异 的，它们的优化目标函数往往不一定能满足连续、可导等特性。针对这种情形，研究人员 将智能优化算法应用到振动筛优化设计中。例如，李智提出了一种基于粒子群算法的往 复振动筛运动参数优化设计方法（李智.基于粒子群算法的往复振动筛运动参数优化 设计[J].粮食与饲料工业，2004, 04:12-13);高瞩提出了一种基于遗传算法的圆形振动 筛减振弹簧的优化设计方法（高瞩.基于遗传算法的圆形振动筛减振弹簧的优化设计 [J].轻工机械，2005, 04:68-71.);王苗等提出了一种基于粒子群优化算法的圆振动筛 的优化设计（王苗，段志善.基于粒子群优化算法的圆振动筛的优化设计[J].煤矿机 械，2011，12:16-18)。
[0004] 和声搜索算法是一种近几年提出来的现代智能优化算法，它已经成功地解决了许 多工程优化问题。但是传统和声搜索算法应用于振动筛优化设计时往往存在着收敛速度 慢，易于陷入局部最优的缺点。
【发明内容】
[0005] 本发明主要是解决现有技术所存在的技术问题；针对传统和声搜索算法应用于振 动筛优化设计时往往存在着收敛速度慢，易于陷入局部最优的缺点，提出一种应用高斯和 声搜索的振动筛优化设计方法。本发明能够加快振动筛优化设计的收敛速度，提高振动筛 优化设计的精度。
[0006] 本发明的技术方案：一种应用高斯和声搜索的振动筛优化设计方法，包括以下步 骤：
[0007] 步骤1，根据待优化设计的振动筛的物理机械特性建立以下形式的最小化优化目 标的数学模型：
[0008] 最小化振动筛的优化设计目标函数f(X)，并满足优化设计约束条件：gk(X) < 0，k =1，2,. . .，M，其中X =[X1 x2?..xD]为振动筛的优化设计变量组成的向量；D为振动筛 待优化设计的变量个数；gk (X) < O为第k个优化设计约束条件，M为优化设计约束条件个 数；
[0009] 步骤2,用户初始化参数，所述初始化参数包括振动筛待优化设计的变量个数D， 和声库大小Popsize，记忆库学习率HMCR，音调调整率PAR，最大评价次数MAX_FEs ;
[0010] 步骤3,当前演化代数t = 0,当前评价次数FEs = 0 ;
[0011] 步骤4,随机产生初始和声库/> k .皮，....,其中：个体下标i = 1，. . .，Popsize，并且K/K;… ]为和声库Pt中的第i个个体，其随机初始 化公式为：
[0012] Bl1： =Loj +rand(0J)x(Up；~Lo.f)
[0013] 其中下标j = 1，. . .，D，并且D表示振动筛有多少个要优化设计的变量；名丨为在和 声库Pt中的第i个个体，存储了 D个优化设计变量的值，rand (0, 1)为在[0, 1]之间服从均 匀分布的随机实数产生函数，Lo jP Up ,分别为第j个优化设计变量的取值范围的下界和上 界；
[0014] 步骤5,采用演化优化领域中常用的惩罚函数法计算和声库Pt中每个个体爽的适 应值其中个体下标i = 1，. . .，Popsize，个体#丨_的适应值朽7(5丨）的计算公式如 下：
[0015]
[0016] 其中pow为幂函数，abs为取绝对值函数，exp为指数函数，max为取最大值函数；
[0017] 步骤6,当前评价次数FEs = FEs+Popsize ;
[0018] 步骤7,保存和声库Pt中适应值最小的个体为最优个体Best N
[0019] 步骤8,结合高斯音调调整策略产生一个试验个体lf，并计算试验个体Ut的适应值 Fit (Ut)，具体步骤如下：
[0020] 步骤8. 1，令计数器j = 1 ;
[0021] 步骤8. 2,在[0, 1]之间随机产生一个实数rl，如果rl小于记忆库学习率HMCR，则 转到步骤8. 3,否则转到步骤8. 11 ;
[0022] 步骤8.3,在[1，Popsize]之间随机产生一个正整数RIl ;
[0023] 步骤 8.4,令 W =爲1>/:;
[0024] 步骤8. 5,在[0, 1]之间随机产生一个实数r2,如果r2小于音调调整率PAR，则转 到步骤8. 6,否则转到步骤8. 12 ;
[0025] 步骤 8. 6,令优秀个体数量 pBestN = max (2, PopsizeXr3)，实数 r3 为[0? 01，0? 2] 之间随机产生的一个实数；
[0026] 步骤8. 7,从当前和声库的前pBestN个优秀个体中随机选择出一个个体;
[0027] 步骤 8. 8,令随机权值 RW = rand (0, 1);
[0028] 步骤8. 9，令均值滅細x+ (1.-犮:,并令方差 SLW =abs(^;l; - ^o,,. ,)x rand((Xl)；
[0029] 步骤 8. 10, 1? =勸麵lRa_(臟anK 〇，其中 NormalRand (MeanV, SDV)表示 以MeanV为均值，SDV为方差的高斯随机数产生函数；然后转到步骤8. 12 ;
[0030] 步骤 8. 11，令以::=乙十rand(),l)x(〇^ -Lo,);
[0031] 步骤8. 12,令计数器j = j+1，如果计数器j小于或等于D，则转到步骤8. 2,否则 转到步骤8. 13 ;
[0032]步骤8. 13,计算试验个体Ut的适应值Fit (U )，转到步骤9 ;
[0033] 步骤9,令当前评价次数FEs = FEs+1 ;
[0034] 步骤10,找到当前和声库中适应值最差个体的下标WorstI ;
[0035] 步骤11，按以下公式在当前和声库中的最差个体与个体Ut两者之间选择出 更优个体进入下一代和声库：
[0036]
[0037] 步骤12,随机选择出一个个体，然后对选择出来的个体执行差分变异局部搜索操 作得到个体Dlf，具体步骤如下：
[0038] 步骤12. 1，在[1，Popsize]之间随机产生一个正整数RI2 ;
[0039] 步骤12. 2,令计数器i = 1，并令差分变异局部搜索最大次数DSN = max (5, D/5)， 其中max为取最大值函数；
[0040] 步骤12. 3,如果计数器i小于或等于DSN，则转到步骤12. 4,否则转到步骤13 ;
[0041] 步骤12. 4,在[1，Popsize]之间随机产生两个互不相等，并且也不等于RI2的正 整数RI3和RI4 ;
[0042] 步骤12. 5,令计数器j = 1 ;
[0043] 步骤12. 6,如果计数器j小于或等于D，则转到步骤12. 7,否则转到步骤12. 9 ;
[0044] 步骤 12. 7,令汉巧=仏+raiid(〇,l)x(S;7;,-义/4 ;.);
[0045] 步骤12. 8,令计数器j = j+1，转到步骤12. 6 ;
[0046] 步骤12. 9,计算差分变异局部搜索操作得到的个体DUt的适应值，并令当前评价次 数 FEs = FEs+1 ;
[0047] 步骤12. 10,如果个体DUt的适应值比个体的适应值更优，则在当前和声库中 用个体DUt替换个体，并转到步骤13,否则转到步骤12. 11 ;
[0048] 步骤12. 11，令计数器i = i+1，转到步骤12. 3 ;
[0049] 步骤13,保存和声库Pt中适应值最小的个体为最优个体Best N
[0050] 步骤14,当前演化代数t = t+1 ;
[0051 ] 步骤15,重复步骤8至步骤14直至当前评价次数FEs达到MAX_FEs后结束，执行 过程中得到的最优个体BestIP为振动筛优化设计的结果。
[0052] 本发明在和声搜索算法中设计了高斯音调调整策略，利用当前和声库中的优秀个 体和随机个体的搜索信息，并结合高斯随机数产生函数来产生新个体，在保持和声库的多 样性的同时尽可能加快收敛速度。此外，利用来差分变异策略来进行局部搜索，增强算法的 局部搜索能力，加快振动筛优化设计的收敛速度；与同类方法相比，本发明能够加快振动筛 优化设计的收敛速度，提高振动筛优化设计的精度。
【附图说明】
[0053] 图1为本发明的流程图。
【具体实施方式】
[0054] 下面通过实施例，并结合附图，对本发明的技术方案作进一步具体的说明。
[0055]实施例：
[0056] 本实施例基于文献（王苗，段志善.基于粒子群优化算法的圆振动筛的优化设计 [J].煤矿机械，2011，32(12)，16-18.)中的圆振动筛的优化设计问题为例，本发明的具体 实施步骤如下：
[0057] 步骤1，根据待优化设计的圆振动筛的物理机械特性建立以下形式的最小化优化 目标的数学模型：
[0058] 最小化圆振动筛的优化设计目标函数：
[0062] 其中X = [X1 X2 ... xD]为圆振动筛的优化设计变量组成的向量，X1为圆振动筛 的振幅，X2为圆振动筛的抛射强度，X 3为圆振动筛的筛面倾角；D = 3为圆振动筛的优化设 计变量个数；gk (X)彡0, k = 1，2, ...，M，为第k个优化设计约束条件，M = 6为优化设计约 束条件个数；
[0063] 步骤2,用户初始化参数，所述初始化参数包括圆振动筛的优化设计变量个数D = 3,和声库大小Popsize = 50,记忆库学习率HMCR = 0.95,音调调整率PAR = 0. 3,最大评价 次数 MAX_FEs = 100000 ;
[0064] 步骤3,当前演化代数t = 0,当前评价次数FEs = 0 ;
[0065] 步骤4，随机产生初始和声库巧={5丨,5丨，...,5丨，.,.，<_^，其中：个体下标1 = 1，. . .，Popsize，并且蹲=[氣!,:1 ……]为和声库Pt中的第i个个体，其随机初始 化公式为：
[0066] = Loj +rand(0,l)x(Vp；-Lo^)
[0067] 其中下标j = 1，. . .，D，并且D表示圆振动筛有多少个要优化设计的变量；S丨为在 和声库Pt中的第i个个体，存储了 D个优化设计变量的值，rand (0, 1)为在[0, 1]之间服从 均匀分布的随机实