基于多机制优化螺旋飞行特征的乌燕鸥算法.pdf

投稿网址:2023 年 第23 卷 第26 期2023,23(26):11299-10科 学 技 术 与 工 程Science Technology and EngineeringISSN 16711815CN 114688/T收稿日期:2022-11-24修订日期:2023-06-21基金项目:国家自然科学基金(62041702);教育部人文社会科学研究项目(20YJA870010);江西省社会科学研究规划项目(19TQ05,21GL12);江西省高校人文社会科学研究项目(TQ20105)第一作者:李光泉(1977),男,汉族,江西鄱阳人,博士,副教授。研究方向:智能算法。E-mail:。通信作者:邵鹏(1983),男,汉族,江西南昌人,博士,讲师。研究方向:智能算法。E-mail:。引用格式:李光泉,刘欣宇,王龙飞,等.基于多机制优化螺旋飞行特征的乌燕鸥算法J.科学技术与工程,2023,23(26):11299-11308.Li Guangquan,Liu Xinyu,Wang Longfei,et al.Sooty tern algorithm based on multi-mechanism optimized spiral flight characteristicsJ.Science Technology and Engineering,2023,23(26):11299-11308.基于多机制优化螺旋飞行特征的乌燕鸥算法李光泉,刘欣宇,王龙飞,邵鹏(江西农业大学计算机与信息工程学院,南昌 330045)摘 要 乌燕鸥优化算法(sooty tern optimization algorithm,STOA)是近年来出现的一种新兴的群智能优化算法,因其独特的螺旋式寻优方式和在实际工程问题中显著的优化效果而备受研究与关注。然而,STOA 算法本身也存在着收敛速度较慢,搜索精度较低,并且容易陷入局部最优等缺点。因此,提出了一种融合混沌映射、自适应惯性权重与高斯变异的多机制乌燕鸥优化算法(multi-mechanism sooty tern optimization algorithm incorporating chaotic mapping,adaptive inertia weight and gaussian mu-tation,GT-STOA),以促进群体多样性并增强算法全局搜索和局部寻优的能力。同时,为验证算法寻优效果的显著性,基于 12个具有不同特征的测试函数,选取 9 种典型的优秀元启发式算法进行对比验证。实验结果表明 GT-STOA 相较于其他 9 种算法具有更高的寻优精度和更快的收敛速度,并且易跳出局部最优找到全局最优解。此外,为了研究 GT-STOA 解决实际问题的能力,对压力容器设计问题进行优化求解,所得实验数据显示 GT-STOA 较传统 STOA 算法在求解精度上提升了 42.54%。关键词 乌燕鸥优化算法;混沌映射;自适应惯性权重;高斯变异中图法分类号 TP301;文献标志码 ASooty Tern Algorithm Based on Multi-mechanismOptimized Spiral Flight CharacteristicsLI Guang-quan,LIU Xin-yu,WANG Long-fei,SHAO Peng(School of Computer and Information Engineering,Jiangxi Agricultural University,Nanchang 330045,China)Abstract The sooty tern optimization algorithm(STOA)is an emerging swarm intelligence optimization algorithm in recent years,which has attracted much research and attention because of its unique spiral optimization search method and the remarkable optimizationeffect in practical engineering problems.However,the STOA has the disadvantages of slower convergence,lower search accuracy,andeasily fall into local optimum.Therefore,a multi-mechanism sooty tern optimization algorithm(GT-STOA),which combines chaoticmapping,adaptive inertia weight and gaussian mutation was proposed to promote population diversity and enhance the global search andlocal optimization capability of the algorithm.Meanwhile,to verify the significant optimization effect of the algorithm,nine typical andexcellent meta-heuristic algorithms were selected for comparative validation on a total of 12 test functions with different features.Theexperimental results show that the GT-STOA has higher optimization accuracy and faster convergence speed than the other nine algo-rithms,and it is easy to jump out of the local optimum to find the global optimum solution.In addition,in order to study the ability ofGT-STOA to solve practical problems,the optimization is conducted on the design of pressure vessel problem,and the experimental da-ta obtained show that the GT-STOA improved the solving accuracy by 42.54%compared with the traditional STOA algorithm.Keywords sooty tern optimization algorithm;chaotic mapping;adaptive inertia weight;gaussian mutation 近些年来,各国学者们由于深受自然界动植物的特征行为而倍受启发,并通过模拟自然界生物的各种觅食行为、攻击行为和繁殖行为等提出了一系列的群体智能优化算法。群智能优化算法种类众多,但因为其结构简单、易于扩充且效果显著而深受广泛学习研究。例如:鲸鱼优化算法1(whaleoptimization algorithm,WOA)、粒子群算法2(parti-cle swarm optimization,PSO)、樽海鞘群体算法3(salp swarm algorithm,SSA)、蝴蝶优化算法4(but-terfly optimization algorithm,BOA)、灰狼优化算法5(grey wolf optimization,GWO)、海鸥优化算法6(seagull optimization algorithm,SOA)、乌燕鸥优化算法7(sooty tern optimization algorithm,STOA)等。其中,乌燕鸥算法就是近年来一种新兴的群智能优投稿网址:化算法,其主要是通过特殊的迁徙行为和攻击行为,实现了算法在全局搜索和局部搜索的主要功能,并因为算法的可操作性强,可实现性高,目前在一些工业问题和分类问题上均取得了良好效果。然而,也因为算法有限的搜索能力和单一的搜索方式,导致了算法较易陷入局部最优,并随着迭代次数的增加会导致种群多样性的不断减小,降低了算法在搜索后期的收敛速度及收敛精度。为此,针对 STOA 存在的不足,中外学者们倾向于从算法混合、单策略优化、多策略优化以及基于实际工程问题优化等多角度提出广泛的改进方法,并在不同的函数集和实际问题上进行仿真实验。贾鹤鸣等8通过将遗传算法和乌燕鸥算法进行混合,并通过平均适应度值进行评估,以此加强乌燕鸥算法的局部搜索能力,证明了该算法在数据挖掘领域具有广阔的工程应用前景;肖永江等9通过将遗传算法的变异思想融入乌燕鸥算法中,提出了一种新颖的改进乌燕鸥算法并进行了 DG 的优化配置,大大提高了算法收敛速度和精度;乔夏君等10针对传统 PID 参数整定存在的问题,通过引入混沌映射提出了一种混沌乌燕鸥算法(chaos sooty ternoptimization algorithm,CSTOA)用以优化发动机参数自整定 PID 控制;Soni 等11通过引入乌燕鸥算法对动态经济排放调度问题进行了优化和解决。同时,孙珂琪等12提出一种混合正余弦算法和 Lvy 飞行的自适应乌燕鸥算法,并将其应用于解决桥式起重机主梁结构优化设计中;雒珊等13为解决算法局部最优能力和寻优能力较低的问题,提出了一种混合Lvy 飞行和热交换混沌乌燕鸥算法,并将算法应用于二级斜齿圆柱齿轮传动机构可靠性轻量化设计,实现了齿轮传动机构轻量化设计的目的;而李家帅等14针对行星齿轮箱振动信号的非线性特征明显、故障特征难以识别的问题,提出了一种基于乌燕鸥算法的故障诊断方法用以提取磨损故障特征频率;李月英15和王国柱等16则在此基础上提出了一种改进的乌燕鸥算法和多策略改进的乌燕鸥算法,有效地解决了移动机器人在复杂环境中的最优路径规划问题和桥式起重机主梁质量减重问题。上述研究不仅在一定程度上平衡、增强了乌燕鸥算法在全局寻优和局部寻优的能力,而且对于算法的收敛速度及收敛精度也有了本质的提升。然而,由于 SOTA 自身独特的寻优方式,导致算法仍然存在着“全局强,局部弱”的显著优缺点。为了能在更大程度上提高算法的局部搜索能力和收敛速度,现提出一种融合混沌映射、自适应惯性权重与高斯变异的多机制乌燕鸥优化算法(GT-STOA)。首先,通过采用 Tent 混沌映射对种群进行初始化;其次,引入自适应惯性权重优化机制,增强算法的局部搜索能力,并在此基础上引入高斯变异机制,此为提升算法跳出局部极值的能力。由此可见,基于多机制协同优化的 GT-STOA 不仅增强了种群的多样性,而且提升了算法在全局搜索和局部搜索的能力,从而使之具备更强的寻优能力和算法鲁棒性。1 乌燕鸥优化算法1.1 乌燕鸥迁徙行为(1)避免碰撞:为了防止乌燕鸥群体间产生相互碰撞,附加变量 SA将会用于计算新的搜索代理位置,以避免相邻乌燕鸥个体之间的碰撞,并计算产生乌燕鸥的新位置:Cst=SAPst(z)(1)式(1)中:Cst为乌燕鸥不与其他个体产生冲突的新位置;Pst为乌燕鸥的当前位置;z 为当前迭代次数;SA为乌燕鸥个体在给定搜索空间中的移动方式。SA=Cf-zCfMaxiterations(2)z=0,1,Maxiterations(3)式中:Cf为一个线性调整 SA的控制变量,它的值从Cf线性减小至 0;z 为当前迭代次数;Maxiterations为最大迭代次数。(2)聚集:乌燕鸥个体在迁徙过程中避免相互碰撞后,群体之间会朝着最佳个体所在的位置方向进行移动。Mst=CBPbst(z)-Pst(z)(4)CB=0.5Rand(5)式中:Mst为乌燕鸥最佳个体所在的位置方向;Pbst为乌燕鸥的最佳位置;CB为一个随机变量,它负责更好地全局探索;Rand为在 0,1 区间的随机数。(3)更新:最终,乌燕鸥可以根据最佳个体所在的位置方向更新其轨迹和最终位置。Dst=Cst+Mst(6)式(6)中:Dst为最佳乌燕鸥个体和当前个体之间的差距。1.2 乌燕鸥攻击行为攻击行为是乌燕鸥种群重要的特征行为之一,在迁徙的途中,它们可以依靠自身重量和羽翼改变它们的速度、高度和攻击角度等。当乌燕鸥群体攻击猎物时,它们便会在空中产生特定的螺旋行为。具体表现形式如下。x=Rsini(7)y=Rcosi(8)z=Ri(9)00311科 学 技 术 与 工 程Science Technology and Engineering2023,23(26)投稿网址:R=uekv(10)式中:R 为乌燕鸥飞行的螺旋半径;i 为在0,2范围内的随机角度;u、v 为形成螺旋飞行的相关常数,设定常数 u、v 值均为 1。由此,根据式(7)式(10)所得乌燕鸥最终位置 Pst更新方式为Pst(z)=Dst(x+y+z)Pbst(z)(11)2 乌燕鸥算法的优化设计2.1 混沌初始化机制乌燕鸥种群的初始位置对于乌燕鸥最优位置的移动搜索会存在较大的影响和一定的约束,在STOA 算法中就初始位置的更新方式而言主要是通过在搜索空间中随机生成位置序列而产生的。随机初始化所得的确定位置,虽然在一定程度上保证了初始位置的随机性,但是就种群的多样性而言同样带来了极大的不确定性,也无法保证个体位置的均匀分布,从而导致种群多样性和寻优速度的降低。通过学者们的不断尝试和复现,并结合混沌映射的遍历性特征和随机性特征,王贺琦等17将 Tent映射用于了优化改进人工蜂群算法,毛清华等18则融合改进 Tent 混沌方法优化了灰狼算法等。以此为鉴,本文研究将引入 Tent 混沌映射机制优化乌燕鸥种群初始化过程,具体映射函数表达式为xt+1=xt,0 xt 0.5(1-xt),0.5 xt 1(12)式(12)中:t 为映射次数,t=0,1,2;xt为第 t 次所得的映射函数值;为混沌参数,在本文中设为定值2。假设在二维搜索空间中,设定种群数量为 200,将随机种群初始化后的位置分布同引入混沌映射后的种群初始化位置分布图进行直观比对。由图 1(a)可以发现,随机初始化种群所得的遍历位置分布散乱,随机性较大。反观图 1(b),在引入 Tent映射优化初始种群位置后所呈现出的分布趋势是较为均匀的,由此直观地反映出 Tent 映射初始化种群具有一定的优化效果。2.2 自适应惯性权重优化机制乌燕鸥优化算法中的惯性权值是衡量算法收敛速度和寻优结果的一个重要指标。当惯性权重取值较大时,算法会因为大权值的存在而表现出显著的全局搜索能力,而当惯性权重取值较小时,算法则更会表现出卓越的局部搜索能力。由此,归于乌燕鸥算法的迭代过程。因为变量 SA在寻优过程中充当着惯性权重的重要作用,所以本文研究将针对变量 SA进行自适应权重优化。通过迭代次数的递增适当地对权值进行自适应缩小,用以增强算法图 1 种群初始化位置分布示意图Fig.1 Schematic diagram of population initializationlocation distribution的局部搜索能力,平衡算法在全局搜索和局部搜索的性能。综上,具体的优化模型如下。SA=Cf-zCfMaxiterations(13)=1,t 0.25 Maxiterations2,0.25 Maxiterations t 0.5Maxiterations4,0.5Maxiterations t 0.75Maxiterations8,0.75Maxiterations t Maxiterations(14)2.3 高斯变异优化机制乌燕鸥优化算法由于小步长搜索特性的存在而表现出了显著的全局搜索能力而被广泛应用于实际问题的优化解决中,但同时也是因为减小了搜索步长给算法的局部搜索带来了一定的负影响。不仅如此,在 STOA 的迭代执行过程中,种群内乌燕鸥位置的多样性会伴随着迭代次数的不断递增而大幅下降,由此削弱了算法精细搜索的能力,使得算法过早的陷入局部最优,降低了算法寻优性能。103112023,23(26)李光泉,等:基于多机制优化螺旋飞行特征的乌燕鸥算法投稿网址:鉴于此,所受李文超等19将高斯变异算子作为自花授粉更新公式,叶坤涛等20通过高斯变异扰动更新种群蜘蛛最优位置为启发,并综于乌燕鸥算法的小步长搜索机理,本文研究在乌燕鸥种群完成迁徙攻击搜索后针对种群内的最优个体位置 Posbest通过高斯变异机制执行优化更新,具体机制模型为Posbest=Posbest0.1+kGaussian(0,1)(15)式(15)中:Posbest为乌燕鸥的最优个体位置;k 为在区间 0,1 中的随机数;Gaussian(0,1)为服从均值为 0、方差为 1 的高斯分布。2.4 GT-STOA 算法流程及描述GT-STOA 算法是一种由多机制优化后的群智能算法,其本质是通过模拟乌燕鸥种群的迁徙行为和攻击行为进行数学模型复现,以达到定位搜索和快速寻优的效果。针对传统乌燕鸥算法存在的不足,通过融合混沌映射,自适应惯性权重和高斯变异扰动展开了多机制辅助优化。主要包括在种群初始化位置的更新,惯性变量权值的自适应性扰动以及在最优个体位置的更新上进行针对性优化,帮助算法丰富了种群多样性,提升了收敛速率和提高了寻优精度等。综上所述,GT-STOA 算法的整体执行过程如下。步骤 1 初始化乌燕鸥种群 xi|i=1,2,N,初始化参数 SA、CB和 Maxiterations,并由式(12)混沌映射函数更新种群初始化位置。步骤 2 计算乌燕鸥个体的适应度值 Pst(x),寻找 适 应 度 值 最 优 所 属 乌 燕 鸥 的 所 处 位 置Pbst(x),记录此时最优解 Xbest。步骤 3 由式(13)和式(14)计算自适应权重值,通过自适应权重的结果更新乌燕鸥个体的位置。步骤 4 由式(1)、式(4)、式(6)计算乌燕鸥个体的运行轨迹和最终位置。步骤 5 由式(7)式(11)计算乌燕鸥最终的螺旋飞行位置 Pst(x+1)。步骤 6 更新计算目标函数适应度值,保留最优乌燕鸥个体位置。步骤 7 由式(15)通过高斯变异对最优乌燕鸥个体位置进行扰动更新,选择保留两者中适应度值更好的解,并更新最优解 Xbest。步骤 8 判断结束条件,如果已满足结束条件,则准备结束程序,输出所得结果,否则跳转步骤 2。步骤 9 程序执行结束,输出最优位置及最优解 Xbest。GT-STOA 的算法伪代码如下。GT-STOA 伪代码Begin初始化 参 数,种 群 规 模 N,测 试 函 数 F,最 大 迭 代 次数 Maxiterations;按式(12)Tent 初始化种群,计算个体适应度值,寻找种群最优位置 Pbst(x);While t Maxiterationsdo for i=1 to N do 按式(13)和式(14)计算自适应权重值,并由此更新种群个体的位置;end for for i=1 to N do 按式(15)对最优个体位置进行扰动更新;end forend WhileEnd2.5 时间复杂度分析算法的时间复杂度对评判算法的性能有着至关重要的作用。为此,计算 STOA 算法在优化后的时间复杂度。当对于种群大小为 N、维度为 D 的启发式算法STOA 而言,其所得的时间复杂度为 O(ND)。而对于多机制优化后的 GT-STOA 而言:Tent 映射初始化参数所需时间的为 ND,自适应惯性权重优化个体位置所需的时间为 ND,高斯变异扰动最优个体位置所需的时间也为 ND。综上,GT-STOA 整体的时间花费为ND+ND+ND=3ND(16)由此可见,本文 GT-STOA 算法的时间复杂度为O(3ND),通过对照算法 STOA 的时间复杂度 O(ND)发现二者的时间复杂度并没有显著性差异,故本文的多策略优化并不会对算法效率造成负面影响。3 实验仿真与结果分析3.1 实验环境及参数设置本次实验的测试环境为 64 位 Windows 10 操作系统,实验软件由 MATLAB R2018b 负责执行算法运算。同时,为了充分证明 GT-STOA 算法的可行性和准确性,也为保证算法执行结果的相对公平性,在本文的仿真实验中将统一算法的测试方法与测试数据。其中,种群数量设定为 30,最大迭代次数为500,并且为了防止算法执行过程的随机性影响,所有实验均独立运行 30 次,记录实验结果的平均值、最优值和标准差。3.2 测试函数为了验证 GT-STOA 的优化效果,选用 12 个基于不同特征的基准测试函数进行仿真实验。其中,F1 F7 是单峰测试函数,F8 F11 是多峰测试函数,F12 则是固定维度的多峰测试函数,测试函数的20311科 学 技 术 与 工 程Science Technology and Engineering2023,23(26)投稿网址:具体信息如下:F1 为 Sphere 函数,F2 为 Schwefel2.22 函数,F3 为 Schwefel 1.2 函数,F4 为 Schwefel2.21 函数,F5 为 Rosenbrock 函数,F6 为 Quartic 函数,F7 为 Sum Squares 函数,F8 为 Rastrigin 函数,F9为 Ackley 函数,F10 为 Griewank 函数,F11 为 Alpine函数,F12 为 Foxholes 函数。3.3 仿真结果分析为验证 GT-STOA 算法的寻优性能,选取 PSO 作为经典智能优化算法的代表。同时,为了衡量 GT-STOA 在新型领域及相关研究领域的寻优效果,选取多策略海鸥优化算法21(multi-mechanism seagulloptimization algorithm,GEN-SOA)、非洲野狗算法22(designoptimizationalgorithm,DOA)、CSTOA、GWO、SSA、BOA、SOA 和 STOA 9 种对照算法进行对比实验。从实验直观性角度出发,将 GT-STOA 同所有对照算法一起,取 F1 F12 测试函数的最优值,平均值和标准差如表 1 中统计所示。同时,根据各算法的寻优结果将绘制部分测试函数的收敛曲线图用于收敛性分析,如图 2 所示。表 1 10 种不同算法的实验对比Table 1 Experimental comparison of 10 different algorithms函数算法最小值平均值标准差函数算法最小值平均值标准差函数算法最小值平均值标准差函数算法最小值平均值标准差STOA3.43 1044.19 1042.90 103STOA4.37 1071.03 1081.61 107SSA4.45 10-082.92 10-74.40 10-7SSA2.73 1013.26 1024.70 102PSO6.23 10-61.39 10-41.39 10-4PSO2.24 1017.73 1013.77 101GWO9.98 10-292.12 10-274.14 10-27GWO2.58 1012.69 1016.16 10-1F1BOA1.05 10-111.29 10-117.97 10-13F5BOA2.88 1012.89 1013.01 10-2SOA1.16 10-11.40 1021.91 102SOA7.09 10-41.99 1057.86 105GEN-SOA0.00 1000.00 1000.00 100GEN-SOA2.35 10-42.55 10-23.77 10-2DOA0.00 1003.51 10-581.89 10-57DOA2.88 1012.89 1014.37 10-2CSTOA6.70 10-201.08 10-152.60 10-15CSTOA2.88 1012.89 1019.19 10-2GT-STOA0.00 1000.00 1000.00 100GT-STOA2.88 1012.89 1014.01 10-2STOA4.08 1016.72 1011.15 101STOA1.02 10-51.02 1003.94 100SSA5.79 10-22.28 1001.47 100SSA5.70 10-21.70 10-18.18 10-2PSO2.72 10-34.07 10-28.86 10-2PSO5.16 10-21.86 10-16.91 10-2GWO2.42 10-178.13 10-174.22 10-17GWO2.97 10-42.04 10-31.18 10-3F2BOA1.44 10-94.60 10-91.21 10-9F6BOA2.91 10-41.46 10-38.23 10-4SOA1.84 10-23.17 1003.78 100SOA2.26 10-31.69 10-13.33 10-1GEN-SOA0.00 1000.00 1000.00 100GEN-SOA1.87 10-51.22 10-31.16 10-3DOA1.79 10-2248.32 10-383.78 10-37DOA1.79 10-53.79 10-43.95 10-4CSTOA3.94 10-112.95 10-95.79 10-9CSTOA1.45 10-75.04 10-63.34 10-6GT-STOA0.00 1000.00 1000.00 100GT-STOA2.13 10-92.11 10-61.93 10-6STOA3.92 1044.94 1044.55 103STOA3.19 1034.77 1036.92 102SSA4.11 1021.59 1031.03 103SSA5.45 10-22.72 1005.73 100PSO3.56 1018.82 1013.44 101PSO1.06 10-43.35 10-37.10 10-3GWO1.19 10-81.91 10-53.97 10-5GWO4.72 10-307.28 10-297.38 10-29F3BOA9.36 10-121.26 10-111.07 10-12F7BOA1.03 10-111.20 10-118.15 10-13SOA4.41 1041.12 1053.98 104SOA5.33 10-31.82 1013.66 101GEN-SOA0.00 1000.00 1000.00 100GEN-SOA0.00 1000.00 1000.00 100DOA0.00 1001.73 10-599.30 10-59DOA0.00 1004.09 10-582.20 10-57CSTOA1.16 10-193.61 10-131.90 10-12CSTOA1.47 10-211.12 10-154.04 10-15GT-STOA0.00 1000.00 1000.00 100GT-STOA0.00 1000.00 1000.00 100STOA6.51 1017.08 1012.59 100STOA4.07 1002.11 1021.04 102SSA4.96 1001.07 1013.65 100SSA2.19 1014.91 1011.62 101PSO6.65 10-11.11 1002.13 10-1PSO3.29 1015.72 1011.64 101GWO1.09 10-76.43 10-75.45 10-7GWO0.00 1003.12 1003.40 101F4BOA5.51 10-96.21 10-94.34 10-10F8BOA0.00 1003.19 1017.09 101SOA9.73 10-16.69 1012.60 101SOA1.11 10-021.41 1012.55 101GEN-SOA0.00 1000.00 1000.00 100GEN-SOA0.00 1000.00 1000.00 100DOA3.09 10-2022.67 10-361.00 10-35DOA0.00 1000.00 1000.00 100CSTOA4.06 10-111.03 10-81.76 10-8CSTOA0.00 1001.89 10-151.02 10-14GT-STOA0.00 1000.00 1000.00 100GT-STOA0.00 1000.00 1000.00 100303112023,23(26)李光泉,等:基于多机制优化螺旋飞行特征的乌燕鸥算法投稿网址:续表函数算法最小值平均值标准差函数算法最小值平均值标准差函数算法最小值平均值标准差函数算法最小值平均值标准差STOA1.93 1011.97 1012.00 10-1STOA2.81 1013.84 1014.51 100SSA1.16 1002.58 1008.25 10-1SSA5.70 10-14.58 1002.05 100PSO2.82 10-32.52 10-15.59 10-1PSO2.89 10-31.49 10-22.34 10-2GWO7.55 10-141.00 10-131.61 10-14GWO5.42 10-175.53 10-47.68 10-4F9BOA4.58 10-95.96 10-94.69 10-10F11BOA2.90 10-101.02 10-99.30 10-10SOA3.29 10-53.32 10-11.30 100SOA9.03 10-56.44 10-11.36 100GEN-SOA8.88 10-168.88 10-160.00 100GEN-SOA0.00 1000.00 1000.00 100DOA8.88 10-168.88 10-160.00 100DOA1.70 10-2381.05 10-375.64 10-37CSTOA6.90 10-119.07 10-92.27 10-8CSTOA1.44 10-122.13 10-103.89 10-10GT-STOA8.88 10-168.88 10-160.00 100GT-STOA0.00 1000.00 1000.00 100STOA3.37 1023.79 1022.34 101STOA9.98 10-11.08 1002.00 10-1SSA6.53 10-42.00 10-21.11 10-2SSA9.98 10-11.10 1002.98 10-1PSO1.27 10-67.24 10-38.31 10-3PSO9.98 10-13.46 1002.73 100GWO0.00 1002.91 10-36.03 10-3GWO9.98 10-14.95 1004.18 100F10BOA1.62 10-125.22 10-122.60 10-12F12BOA9.98 10-11.14 1003.03 10-1SOA7.13 10-42.94 1004.27 100SOA9.98 10-17.56 1008.74 100GEN-SOA0.00 1000.00 1000.00 100GEN-SOA9.98 10-11.60 1009.05 10-1DOA0.00 1000.00 1000.00 100DOA9.98 10-13.30 1002.67 100CSTOA0.00 1001.02 10-144.39 10-14CSTOA9.98 10-11.14 1003.06 10-1GT-STOA0.00 1000.00 1000.00 100GT-STOA9.98 10-11.19 1004.27 10-13.3.1 算法精度分析针对实验的单峰测试函数 F1 F4 而言,GT-STOA 和 GEN-SOA 均在算法执行结果中找到了理论最优值,并且 GT-STOA 相对于其他对照算法而言也实现了精度数值最优的结果,并在 F7、F8、F10 和F12 函数也都寻及理论最优值。对于单峰测试函数F5,SOA 与 GEN-SOA 算法在最终的寻优结果的精度上取得了较好的效果,而 STOA 则相对来说寻优效果较弱。对于测试函数 F6,GT-STOA 通过自适应惯性权重的优化表现出了最优的寻优结果,而 CS-TOA 则仅次于此,在数值精度上平均高出近 2 个数量级左右。然而对于多峰测试函数 F8 F11,就 F8可以发现,GT-STOA 与 5 个算法同样在函数值上找到了理论最优值。对于 F9 而言,GT-STOA 与 GEN-SOA 和 DOA 所能寻及的最优值相同,但同样是所有对照算法中的最优,并且得益于算法的多机制辅助优化,帮助它们所得的标准差均为 0,提升了算法在寻优过程的稳定性,而 GWO,BOA 等则紧随其后。对于测试函数 F10,在最终结果上 GT-STOA、GEN-SOA、GWO、DOA 和 CSTOA 均找到了理论最优值,并且 GT-STOA、GEN-SOA 和 DOA 在平均值和算法稳定性指标上的表现也十分出色。针对测试函数F11 中 GT-STOA 与 GEN-SOA 在结果中找到了理论最优。而在 F12 函数中,虽然 SSA 表现出了较全面的寻优能力,但总体而言,GT-STOA 的寻优能力依旧出色。3.3.2 算法收敛性分析结合收敛图 2 和表 1 可以发现,针对测试函数F1 F12 而言,GT-STOA、GEN-SOA 和 DOA 3 种算法的全局最优值均以相对较高的速率进行收敛下降,并且算法在全局勘测和局部寻优的协调能力中有着显著的优越性,实例表现在测试函数 F5 和 F6等在后期依旧能够通过有力的局部搜索找到最优值。其中,对于测试函数 F1、F3 和 F10,GT-STOA 均在 250 以内便收敛到了全局最优值,并且表现出了收敛速度最快,收敛精度最小的显著效果。而在测试函数 F8 和 F10 中所表现出的收敛性最强,在 50代后便收敛到了算法最优值也是全局最优值。展现了 GT-STOA 算法稳定且显著的全局寻优能力。3.4 高维函数的实验分析综合上述各算法的精度分析及收敛性分析发现,GT-STOA 在低维度寻优过程中有着显著性效果。基于此,为充分验证算法在高维度条件下所显现的寻优能力,对 6 个测试函数在 100 维情况下进行仿真实验,比较 GT-STOA 与各类代表性算法在高维情况下的寻优能力。算法的实验参数与 3.1 节一致,实验结果如表 2 所示。由表 2 的实验结果可知,当测试函数的维度提升至 100 维时,GT-STOA 针对不同类型的单峰函数和多峰函数依旧能在不同的对照算法中表现出更优的寻优能力。在测试函数 F2、F4、F10 和 F11 中GT-STOA 也找到了函数的理论最优值,并且对于实验结果的均值和标准差而言,也展现出了算法在高维情况下显著的稳定性。由上述分析可知,GT-STOA 在处理高维优化问题方面也具有一定的竞争优势。40311科 学 技 术 与 工 程Science Technology and Engineering2023,23(26)投稿网址:图 2 对比算法收敛图Fig.2 Comparison of algorithm convergence graphs3.5 非参数统计检验为避免上述所得的算法结果会出现片面性且无法进行统一定性分析等问题,同时为了更好地验证上述结论,如表3 所示,本文介绍了非参数统计学中的Wilcoxon 检验23和 Friedman 检验24的统计结果。在 Wilcoxon 统计检验中,若显著性水平 P 0.05,则表示算法对于对比算法具有显著性优势。而在 Friedman 统计检验中,将根据排名值进行排序,排名值越低则算法排