CN107747930B - 一种基于万有引力加速布谷鸟算法的圆度误差评定方法 - Google Patents

一种基于万有引力加速布谷鸟算法的圆度误差评定方法 Download PDF

Info

Publication number
CN107747930B
CN107747930B CN201710873935.4A CN201710873935A CN107747930B CN 107747930 B CN107747930 B CN 107747930B CN 201710873935 A CN201710873935 A CN 201710873935A CN 107747930 B CN107747930 B CN 107747930B
Authority
CN
China
Prior art keywords
algorithm
individual
evolution
value
substance
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201710873935.4A
Other languages
English (en)
Other versions
CN107747930A (zh
Inventor
傅文渊
凌朝东
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huaqiao University
Original Assignee
Huaqiao University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Huaqiao University filed Critical Huaqiao University
Priority to CN201710873935.4A priority Critical patent/CN107747930B/zh
Publication of CN107747930A publication Critical patent/CN107747930A/zh
Application granted granted Critical
Publication of CN107747930B publication Critical patent/CN107747930B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
    • G01B21/20—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring contours or curvatures, e.g. determining profile
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/004—Artificial life, i.e. computing arrangements simulating life
    • G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Computation (AREA)
  • General Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Health & Medical Sciences (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

一种基于万有引力加速布谷鸟算法的圆度误差评定方法，基于万有引力搜索无需学习外部环境因素的变化亦能感知全局最优信息的特点，将布谷鸟巢穴赋予不同的个体质量，因此其在优化过程中遵循万有引力定律；对布谷鸟算法中的Levy飞行随机游动和偏好随机游动同时利用优化个体间存在的万有引力进行加速搜索，获得相应的布谷鸟寄生巢穴个体更新位置；在万有引力的作用下，有效的平衡了布谷鸟算法的全局搜索能力和局部寻优能力，避免算法执行末期陷入局部极值点而出现的迟滞现象，提高算法的全局搜索效率和收敛精度。本发明待优化的圆度误差E能快速趋于稳定的最优值，并求解出两个同心圆的理想圆心和从而使得两同心圆之间的区域为最小区域。

Description

一种基于万有引力加速布谷鸟算法的圆度误差评定方法
技术领域
本发明涉及圆周误差评定技术领域，更具体地说，涉及一种基于万有引力加速布谷鸟算法的圆度误差评定方法。
背景技术
圆度误差评定是回转体零件加工过程重要的评价指标，其精度的高低直接影响加工工件的质量及使用寿命。评定圆度误差经典的方法有最小二乘法、最小外接圆法、最大内切圆法和最小区域圆法。最小二乘法理论相对成熟，但是涉及非线性运算，算法复杂度高。同时，其评定结果误差较大，不能满足最小条件的圆度误差评定结果。最小外接圆法是把实际被测圆的最小外接圆作为外包圆，以最小外接圆的圆心为中心作实际圆的内包圆，因此，圆度误差为此两圆的半径差。最大内切圆法与最小外接圆法评定误差原理类似，都存在工件评定区域过大，造成评定过程累积的误差增大，导致评定的圆度误差增大。最小区域圆法是国际通用的误差评定方法。相比于其他3种方法，最小区域圆法具有更高的测量精度，并且是唯一的。它能够真实的反映加工工件的圆度误差数值，因此本发明中圆度误差评定采用最小区域圆法。
布谷鸟算法(Cuckoo Search,CS)作为一种新型优化算法，与其他智能启发式算法类似，也存在搜索效率不高的问题，尤其是易于陷入局部极值点，算法执行后期出现迟滞现象。为此国内外学者对该算法进行相应的改进以提高算法的收敛性能。归纳起来，目前学者研究的改进算法主要包括两方面：布谷鸟算法控制参数的改进以及与其它算法融合的混合策略。文献“Ong P.Adaptive cuckoo search algorithm for unconstrainedoptimization[J].The Scientific World Journal,2014,14(9):1-8.”对CS算法的搜索步长设置为自适应变化，提高了算法的收敛性能。文献“Naik M,Nath M R,Wunnava A,etal.A new adaptive Cuckoo search algorithm[C].IEEE International Conference onRecent Trends in Information Systems.2015:1-5.”在文献文献“Ong P.Adaptivecuckoo search algorithm for unconstrained optimization[J].The ScientificWorld Journal,2014,14(9):1-8.”的基础上进一步研究，将搜索步长的更新调整为由函数适应度值自适应的确定，而不依赖Levy分布的步长因子。文献“Wang Lijin,Yin Yilong,etal.Cuckoo search with varied scaling factor[J].Frontiers of Computer Science,2015,9(4):623-635.”设置Levy分布的步长因子为均匀分布的随机数，极大提高了CS算法的全局搜索性能。但该方法没有考虑局部搜索的性能，导致算法执行后期易发生迟滞现象。文献“王李进,尹义龙,钟一文.逐维改进的布谷鸟搜索算法[J].软件学报,2013,24(11):2687-2698.”提出一种逐维更新的函数评价策略，将各维逐一更新，同时采用贪婪更新方式接受改善解。该算法具有一定的竞争力，但是采用逐维更新的评价策略导致函数评价次数急剧增加，因此算法执行效率较低。文献“Walton S,Hassan O,Morgan K,et al.Modifiedcuckoo search:A new gradient free optimisation algorithm[J].Chaos Solitons&Fractals,2011,44(9):710-718.”将发现概率设置为动态变化数值，增强了算法的收敛性能。但是由于改进的算法搜索步长可容许的范围过小，因此算法全局收敛性能较差。另一方面，针对布谷鸟算法与其它算法融合的混合策略，不少学者取得了较好的研究成果。文献“Wang G G,Gandomi A H,Chu H C,et al.Hybridizing harmony search algorithm withcuckoo search for global numerical optimization[J].Soft Computing,2016,20(1):1-13.”将和声优化算法与布谷鸟算法结合，在算法执行过程中增加和声算法中的变异操作判定结果，并把该结果反馈到搜索过程，加快了算法的收敛速度。文献“Guo J,Sun Z,TangH,et al.Hybrid optimization algorithm of particle swarm optimization andcuckoo search for preventive maintenance period optimization[J].DiscreteDynamics in Nature&Society,2016,1(1):1-12.”借鉴粒子群算法良好的局部优化性能，在CS算法中引入粒子群组件，提高了CS算法的收敛精度。文献“Li XT,Yin MH.Parameterestimation for chaotic systems using the cuckoo search algorithm with anorthogonal learning method[J].Chinese Physics B,2012,21(5):113-118.”和“Li XT,Wang J,Yin MH.Enhancing the performance of cuckoo search algorithm usingorthogonal learning method[J].Neural Computing&Applications,2014,24(6):1233-1247.”在Levy飞行之后的偏好随机扰动引入一种正交学习机制，以提高CS算法的整体搜索性能。综合起来，此类改进算法虽然提高了搜索过程中的综合性能，例如种群的多样性，收敛精度等，但代价是增加了函数评价次数及算法的复杂度。当处理高维度优化问题时，特别是多峰函数、脊峰函数和奇异函数时，算法的自适应调节能力较差，导致寻优效果不理想。同时，现有的此类改进算法通常是将两种算法机械的结合在一起，没有深入挖掘种群寻优过程的内部机理，算法执行效率较低。因此，研究新的改进算法具有积极的意义。
发明内容
本发明的目的是在现有的方法上做改进，提供一种基于万有引力加速布谷鸟算法的圆度误差评定方法，基于万有引力搜索无需学习外部环境因素的变化亦能感知全局最优信息的特点，将布谷鸟巢穴赋予不同的个体质量，因此其在优化过程中遵循万有引力定律；对布谷鸟算法中的Levy飞行随机游动和偏好随机游动同时利用优化个体间存在的万有引力进行加速搜索，获得相应的布谷鸟寄生巢穴个体更新位置；在万有引力的作用下，有效的平衡了布谷鸟算法的全局搜索能力和局部寻优能力，避免算法执行末期陷入局部极值点而出现的迟滞现象，提高算法的全局搜索效率和收敛精度。本发明将万有引力加速布谷鸟算法应用于圆度误差评定，具有较好的鲁棒性，待优化的圆度误差E能快速趋于稳定的最优值，并求解出两个同心圆的理想圆心和从而使得两同心圆之间的区域为最小区域。
本发明采用如下技术方案：
一种基于万有引力加速布谷鸟算法的圆度误差评定方法，包括如下步骤：
(1)优化的圆度误差函数为E(X(1),X(2))＝min(Ry-Rx)，其中,(X(1),X(2))为待优化的圆心，E为圆度误差，Rx和Ry分别为嵌套在工件实际轮廓的两个同心圆的半径，(ui,vi)为工件轮廓测量点的坐标值，i∈[1,s]，s为测量点数；设置算法的最大进化代数为W、发现概率为Pa、初始时刻万有引力系数为G0、算法控制参数为θ、γ0、p、宿主巢穴的初始移动速度为V0、初始宿主巢穴位置X0,m,m∈[1,N]、N为种群个数、优化的圆度误差函数空间维度为D＝2；
(2)根据E(X(1),X(2))＝min(Ry-Rx)计算初宿主巢穴位置X0,m对应的适应度函数值
(3)计算第r次进化时物质个体的万有引力常数Gr；同时计算当前优化函数的适应度最优值Er,best和最差值Er,worst，以及对应的最优解Xr,gb，其中r为进化代数；
(4)计算出作用物质个体m在第r次进化时的惯性质量Mr,pm和作用物质个体k在第r次进化时的惯性质量Mr,ak；
(5)根据Gr、Mr,pm、Mr,ak、第r代第m个种群的候选解Xr,m的第j维度上的数值及第r代第k个种群的候选解Xr,k的第j维度上的数值计算当前进化代数的宿主巢穴所受到的万有引力合力和加速度其中，j表示优化的圆度误差函数的第j维度，j＝1,2；
(6)计算得到Levy飞行随机游动方式的宿主巢穴Xr+1,m，同时，按照发现概率Pa舍弃一部分第r+1代宿主巢穴位置Xr+1,m；其中，表示点对点乘法，L(β)服从Levy概率分布；
(7)计算Xr+1,m＝Xr,gb+p·(Xr,k-Xr,z-ar,m)，得到偏好随机游动方式产生的宿主巢穴Xr+1,m并替换步骤(6)中被舍弃的相同部分的第r+1代宿主巢穴位置Xr+1,m，其中k,z∈[1,N]且k,z均为随机整数；
(8)计算步骤(7)中种群产生的宿主巢穴位置Xr+1,m对应的适应度函数值同时更新当前最优值Er+1,best和最差值Er+1,worst，以及对应的最优解Xr+1,gb；
(9)若满足算法进化代数为W，则输出当前进化的最优解和并停止算法，否则返回步骤(3)重复执行算法，和即为两个同心圆的理想圆心。
优选的，适应度最优值Er,best和最差值Er,worst的求解过程如下：
其中，Er,m是圆度误差的适应度值，表示进化至第r代的第m个物质个体的适应度值。
优选的，惯性质量Mr,pm和惯性质量Mr,ak的求解过程如下：
其中，Er,k是圆度误差的适应度值，表示进化至第r代的第k个物质个体的适应度值。
优选的，万有引力合力和加速度的求解过程如下：
其中，dj为区间(0,1)内均匀分布的一个随机数，b为物质个体质量按降序排列后的前排个体数目；表示在维度j上，物质个体m对物质个体k的万有引力；
由牛顿第二定律，物质个体m在维度j上第r次进化时的加速度定义如下：
其中，Mr,mm为进化至第r代时作用于物质个体m本身的惯性质量。
优选的，万有引力的求解过程如下：
其中，ε为一个无穷小的常数，Mr,pm表示被作用物质个体m在第r次进化时的惯性质量，Mr,ak表示为作用物质个体m在第r次进化时的惯性质量，Rr,mk为物质个体m和物质个体k的欧氏空间距离，即Rr,mk＝||Xr,m,Xr,k||2，Gr表示第r次进化时物质个体的万有引力常数。
优选的，万有引力常数Gr的求解过程如下：
Gr＝G(G0,r)＝G0·e-θ·r/W
其中，G0表示进化初始时刻物质个体的万有引力系数，θ为算法控制参数。
优选的，控制参数θ＝20。
优选的，L(β)的求解过程如下：
L(β)～u＝t-1-β,0<β≤2
搜索路径与时间t具有幂次方概率密度函数的特性，通过数学代换得到，
其中β＝1.5,u和v服从标准高斯分布，即u,v∈N(0,1)。
提供一种基于万有引力加速布谷鸟算法(Gravitational Acceleration Searchbased Cuckoo Search,GASCS)的圆度误差评定方法，将布谷鸟算法(Cuckoo Search,CS)和万有引力搜索算法(Gravitational search algorithm,GSA)相结合应用到圆度误差评定上，在万有引力的作用下，有效的平衡了布谷鸟算法的全局搜索能力和局部寻优能力，避免算法执行末期陷入局部极值点而出现的迟滞现象，提高算法的全局搜索效率和收敛精度，待优化的圆度误差E能快速趋于稳定的最优值，并求解出两个同心圆的理想圆心和从而使得两同心圆之间的区域为最小区域。
附图说明
图1是本发明的GASCS算法与GSA及CS算法的圆度误差随进化代数变化图；
图2是本发明的GASCS算法与GSA及CS算法运行次数与圆度误差随变化图。
具体实施方式
以下通过具体实施方式对本发明作进一步的描述。
圆度误差属于形位误差，由国家标准规定需要符合理想位置的最小化条件。圆度误差的最小区域圆法是用两个同心圆覆盖评定的加工工件实际轮廓，两同心圆之间的区域为最小区域，则圆度误差为这两同心圆的半径差。假设(ui,vi),(i＝1,2,…,s)为工件轮廓测量点的坐标值，s为测量点的数目。同时，假设理想圆的圆心为待优化的圆心为(X(1),X(2))，圆度误差为E。
根据最小区域圆的误差评定可得两同心圆之间的区域为最小区域。因此需要满足以下数学关系式：
其中，Rx，Ry分别为嵌套在工件实际轮廓的两个同心圆的半径。由如上两式可以得到，因此Ry≥Rx。则两个同心圆的半径差为
r＝Ry-Rx
因此，算法搜索的目的是在满足工件轮廓测量点坐标值对应的及条件下，尽可能使r＝Ry-Rx最小化，也即是求解E(X(1),X(2))。因此圆度误差评定问题转化为求解目标函数极小值问题。
基于上述要解决的问题，本发明一种基于万有引力加速布谷鸟算法的圆度误差评定方法，包括如下步骤：
(1)优化的圆度误差函数为E(X(1),X(2))＝min(Ry-Rx)，其中,(X(1),X(2))为待优化的圆心，E为圆度误差，Rx和Ry分别为嵌套在工件实际轮廓的两个同心圆的半径，(ui,vi)为工件轮廓测量点的坐标值，i∈[1,s]，s为测量点数；设置算法的最大进化代数为W、发现概率为Pa、初始时刻万有引力系数为G0、算法控制参数为θ、γ0、p、宿主巢穴的初始移动速度为V0、初始宿主巢穴位置X0,m,m∈[1,N]、N为种群个数、优化的圆度误差函数空间维度为D＝2；
(2)根据E(X(1),X(2))＝min(Ry-Rx)计算初宿主巢穴位置X0,m对应的适应度函数值
(3)计算第r次进化时物质个体的万有引力常数Gr；同时计算当前优化函数的适应度最优值Er,best和最差值Er,worst，以及对应的最优解Xr,gb，其中r为进化代数；
(4)计算出作用物质个体m在第r次进化时的惯性质量Mr,pm和作用物质个体k在第r次进化时的惯性质量Mr,ak；
(5)根据Gr、Mr,pm、Mr,ak、第r代第m个种群的候选解Xr,m的第j维度上的数值及第r代第k个种群的候选解Xr,k的第j维度上的数值计算当前进化代数的宿主巢穴所受到的万有引力合力和加速度其中，j表示优化的圆度误差函数的第j维度，j＝1,2；
(6)计算得到Levy飞行随机游动方式的宿主巢穴Xr+1,m，同时，按照发现概率Pa舍弃一部分第r+1代宿主巢穴位置Xr+1,m；其中，表示点对点乘法，L(β)服从Levy概率分布；
(7)计算Xr+1,m＝Xr,gb+p·(Xr,k-Xr,z-ar,m)，得到偏好随机游动方式产生的宿主巢穴Xr+1,m并替换步骤(6)中被舍弃的相同部分的第r+1代宿主巢穴位置Xr+1,m，其中k,z∈[1,N]且k,z均为随机整数；
(8)计算步骤(7)中种群产生的宿主巢穴位置Xr+1,m对应的适应度函数值同时更新当前最优值Er+1,best和最差值Er+1,worst，以及对应的最优解Xr+1,gb；
(9)若满足算法进化代数为W，则输出当前进化的最优解和并停止算法，否则返回步骤(3)重复执行算法，和即为两个同心圆的理想圆心。
优选的，适应度最优值Er,best和最差值Er,worst的求解过程如下：
其中，Er,m是圆度误差的适应度值，表示进化至第r代的第m个物质个体的适应度值。
优选的，惯性质量Mr,pm和惯性质量Mr,ak的求解过程如下：
其中，Er,k是圆度误差的适应度值，表示进化至第r代的第k个物质个体的适应度值。
优选的，万有引力合力和加速度的求解过程如下：
其中，dj为区间(0,1)内均匀分布的一个随机数，b为物质个体质量按降序排列后的前排个体数目；表示在维度j上，物质个体m对物质个体k的万有引力；
由牛顿第二定律，物质个体m在维度j上第r次进化时的加速度定义如下：
其中，Mr,mm为进化至第r代时作用于物质个体m本身的惯性质量。
优选的，万有引力的求解过程如下：
其中，ε为一个无穷小的常数，Mr,pm表示被作用物质个体m在第r次进化时的惯性质量，Mr,ak表示为作用物质个体m在第r次进化时的惯性质量，Rr,mk为物质个体m和物质个体k的欧氏空间距离，即Rr,mk＝||Xr,m,Xr,k||2，Gr表示第r次进化时物质个体的万有引力常数。
优选的，万有引力常数Gr的求解过程如下：
Gr＝G(G0,r)＝G0·e-θ·r/W
其中，G0表示进化初始时刻物质个体的万有引力系数，θ为算法控制参数。
优选的，控制参数θ＝20。
优选的，L(β)的求解过程如下：
L(β)～u＝t-1-β,0<β≤2
搜索路径与时间t具有幂次方概率密度函数的特性，通过数学代换得到，
其中β＝1.5,u和v服从标准高斯分布，即u,v∈N(0,1)。
为了验证本发明提出的基于GASCS算法的圆度误差评定，将该算法同GSA算法和CS算法同时应用于圆度误差评定。算法初始化设置：种群为N＝20个宿主巢穴，发现概率Pa＝0.2，控制参数θ＝20，γ0＝0.1和p＝1，s＝100，算法最大进化代数W＝100。GASCS算法、GSA算法和CS算法的圆度误差随进化代数变化关系如图1所示。同时，对这3种算法分别独立运行30次，每次优化结果与运行次数的变化关系如图2所示。
由图1和图2可知，GASCS算法较CS算法和GSA算法具有明显的寻优优势，圆度误差评定结果最小。GSA算法在圆度误差评定中的算法搜索性能与CS算法相比，具有明显的差距，搜索性能不佳。并且在30次算法独立运行中，每次圆度误差评定的结果偏差过大，即算法的鲁棒性较差，使得寻优结果不够稳定。对加工工件的圆度误差评定的结果表明，该算法是一种高效的圆度误差评定算法。
上述实施例仅仅用以说明本发明，而并非用作对本发明的限定。只要是依据本发明的技术实质，对上述实施例进行变化、变型等都将落在本发明的权利要求范围内。

Claims (1)

1.一种基于万有引力加速布谷鸟算法的圆度误差评定方法，其特征在于，包括如下步骤：
(1)优化的圆度误差函数为E(X(1)，X(2))＝min(Ry-Rx)，其中，(X(1)，X(2))为待优化的圆心，E为圆度误差，Rx和Ry分别为嵌套在工件实际轮廓的两个同心圆的半径，(ui，vi)为工件轮廓测量点的坐标值，i∈[1，s]，s为测量点数；设置算法的最大进化代数为W、发现概率为Pa、初始时刻万有引力系数为G0、算法控制参数为θ、γ0、p、宿主巢穴的初始移动速度为V0、初始宿主巢穴位置X0，m，m∈[1，N]、N为种群个数、优化的圆度误差函数空间维度为D＝2；
(2)根据E(X(1)，X(2))＝min(Ry-Rx)计算初宿主巢穴位置X0，m对应的适应度函数值
(3)计算第r次进化时物质个体的万有引力常数Gr；同时计算当前优化函数的适应度最优值Er，best和最差值Er，worst，以及对应的最优解Xr，gb，其中r为进化代数；
(4)计算出作用物质个体m在第r次进化时的惯性质量Mr，pm和作用物质个体k在第r次进化时的惯性质量Mr，ak；
(5)根据Gr、Mr，pm、Mr，ak、第r代第m个种群的候选解Xr，m的第j维度上的数值及第r代第k个种群的候选解Xr，k的第j维度上的数值计算当前进化代数的宿主巢穴所受到的万有引力合力和加速度其中，j表示优化的圆度误差函数的第j维度，j＝1，2；
(6)计算得到Levy飞行随机游动方式的宿主巢穴Xr+1，m，同时，按照发现概率Pa舍弃一部分第r+1代宿主巢穴位置Xr+1，m；其中，表示点对点乘法，L(β)服从Levy概率分布；
(7)计算Xr+1，m＝Xr，gb+p·(Xr，k-Xr，z-ar，m)，得到偏好随机游动方式产生的宿主巢穴Xr+1，m并替换步骤(6)中被舍弃的相同部分的第r+1代宿主巢穴位置Xr+1，m，其中k，z∈[1，N]且k，z均为随机整数；
(8)计算步骤(7)中种群产生的宿主巢穴位置Xr+1，m对应的适应度函数值同时更新当前最优值Er+1，best和最差值Er+1，worst，以及对应的最优解Xr+1，gb；
(9)若满足算法进化代数为W，则输出当前进化的最优解和并停止算法，否则返回步骤(3)重复执行算法，和即为两个同心圆的理想圆心；
适应度最优值Er，best和最差值Er，worst的求解过程如下：
其中，Er，m是圆度误差的适应度值，表示进化至第r代的第m个物质个体的适应度值；
惯性质量Mr，pm和惯性质量Mr，ak的求解过程如下：
其中，Er，k是圆度误差的适应度值，表示进化至第r代的第k个物质个体的适应度值；
万有引力合力和加速度的求解过程如下：
其中，dj为区间(0，1)内均匀分布的一个随机数，b为物质个体质量按降序排列后的前排个体数目；表示在维度j上，物质个体m对物质个体k的万有引力；
由牛顿第二定律，物质个体m在维度j上第r次进化时的加速度定义如下：
其中，Mr，mm为进化至第r代时作用于物质个体m本身的惯性质量；
万有引力的求解过程如下：
其中，ε为一个无穷小的常数，Mr，pm表示被作用物质个体m在第r次进化时的惯性质量，Mr，ak表示为作用物质个体m在第r次进化时的惯性质量，Rr，mk为物质个体m和物质个体k的欧氏空间距离，即Rr，mk＝||Xr，m，Xr，k||2，Gr表示第r次进化时物质个体的万有引力常数；
万有引力常数Gr的求解过程如下：
Gr＝G(G0，r)＝G0·e-θ·r/W
其中，G0表示进化初始时刻物质个体的万有引力系数，θ为算法控制参数；
控制参数θ＝20；
L(β)的求解过程如下：
L(β)～u＝t-1-β，0＜β≤2
搜索路径与时间t具有幂次方概率密度函数的特件，通过数学代换得到，
其中β＝1.5，u和v服从标准高斯分布，即u，v∈N(0，1)。
CN201710873935.4A 2017-09-25 2017-09-25 一种基于万有引力加速布谷鸟算法的圆度误差评定方法 Active CN107747930B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710873935.4A CN107747930B (zh) 2017-09-25 2017-09-25 一种基于万有引力加速布谷鸟算法的圆度误差评定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710873935.4A CN107747930B (zh) 2017-09-25 2017-09-25 一种基于万有引力加速布谷鸟算法的圆度误差评定方法

Publications (2)

Publication Number Publication Date
CN107747930A CN107747930A (zh) 2018-03-02
CN107747930B true CN107747930B (zh) 2019-12-31

Family

ID=61255338

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710873935.4A Active CN107747930B (zh) 2017-09-25 2017-09-25 一种基于万有引力加速布谷鸟算法的圆度误差评定方法

Country Status (1)

Country Link
CN (1) CN107747930B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108804384B (zh) * 2018-06-21 2020-02-07 哈尔滨工业大学 一种用于最小区域球度评定的最优导向自适应搜索方法
CN109115161B (zh) * 2018-06-21 2020-02-07 哈尔滨工业大学 一种基于空间定向收缩及改进布谷鸟搜索算法的球度评定方法
CN108804390B (zh) * 2018-06-21 2020-02-07 哈尔滨工业大学 一种基于改进布谷鸟搜索策略的最小区域球度评定方法
CN109323677B (zh) * 2018-08-21 2020-08-11 上海隧道工程有限公司 改进布谷鸟搜索算法的圆度误差评定算法
CN109931903B (zh) * 2019-02-26 2021-07-06 上海大学 一种基于改进鲸鱼优化算法的圆柱度评定方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102506863A (zh) * 2011-11-07 2012-06-20 北京航空航天大学 一种基于万有引力搜索的无人机航路规划方法
CN103488825A (zh) * 2013-09-10 2014-01-01 太原理工大学 一种用于齿轮齿廓曲线重构的数据拟合方法
CN103646178A (zh) * 2013-12-18 2014-03-19 中国石油大学(华东) 一种基于改进引力搜索算法的多目标优化方法
CN103760907A (zh) * 2013-12-30 2014-04-30 哈尔滨工程大学 一种基于布谷鸟搜索算法的水下潜器三维路径规划方法
CN107169557A (zh) * 2017-05-12 2017-09-15 淮阴师范学院 一种对布谷鸟优化算法进行改进的方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060173874A1 (en) * 2005-02-03 2006-08-03 Yunqiang Chen Method and system for interactive parameter optimization using multi-dimensional scaling

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102506863A (zh) * 2011-11-07 2012-06-20 北京航空航天大学 一种基于万有引力搜索的无人机航路规划方法
CN103488825A (zh) * 2013-09-10 2014-01-01 太原理工大学 一种用于齿轮齿廓曲线重构的数据拟合方法
CN103646178A (zh) * 2013-12-18 2014-03-19 中国石油大学(华东) 一种基于改进引力搜索算法的多目标优化方法
CN103760907A (zh) * 2013-12-30 2014-04-30 哈尔滨工程大学 一种基于布谷鸟搜索算法的水下潜器三维路径规划方法
CN107169557A (zh) * 2017-05-12 2017-09-15 淮阴师范学院 一种对布谷鸟优化算法进行改进的方法

Also Published As

Publication number Publication date
CN107747930A (zh) 2018-03-02

Similar Documents

Publication Publication Date Title
CN107747930B (zh) 一种基于万有引力加速布谷鸟算法的圆度误差评定方法
CN106447098B (zh) 一种光伏超短期功率预测方法及装置
CN110942205B (zh) 一种基于himvo-svm的短期光伏发电功率预测方法
CN105843189B (zh) 一种用于半导体生产线基于简化仿真模型的高效调度规则选择方法
CN106067034B (zh) 一种基于高维矩阵特征根的配电网负荷曲线聚类方法
CN108053077A (zh) 一种基于区间二型t-s模糊模型的短期风速预测方法与系统
Shen et al. 3D printing build orientation optimization for flexible support platform
CN109670687B (zh) 一种基于粒子群优化支持向量机的质量分析方法
CN108363622B (zh) 一种基于多核dsp运算的无源时差定位算法的实现方法
Liao et al. Short term load forecasting and early warning of charging station based on PSO-SVM
CN110363214A (zh) 一种基于gwa-svm的机器人装配的接触状态识别方法
CN115841176A (zh) 充电桩变频负荷的安全态势感知方法及系统
CN109460875B (zh) 基于果蝇优化算法的电力系统经济调度方法
Wang et al. Otsu multi-threshold image segmentation algorithm based on improved particle swarm optimization
CN114783215B (zh) 无人机分簇方法、装置和电子设备
CN110059871B (zh) 光伏发电功率预测方法
CN107145911A (zh) 一种基于reABC‑SVM的葡萄酒质量鉴别方法
CN114462684B (zh) 耦合数值天气预报与实测数据的风速多点同步预测方法
CN116542139A (zh) 一种液体喷射抛光表面粗糙度的预测方法及装置
CN106960260B (zh) 一种便于电力调度的风电功率预测系统
CN109683558A (zh) 基于粒子群优化的生产过程系统多目标优化模型控制方法
Liu et al. Chaos adaptive improved particle swarm algorithm for solving multi-objective optimization
CN110991743B (zh) 一种基于聚类分析和优化神经网络的风电功率短期组合预测方法
Ji et al. Regional power load forecasting based on PSOSVM
Zheng et al. Support vector machine model based on glowworm swarm optimization in application of vibrant fault diagnosis for hydro-turbine generating unit

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant