一种可控加速性能的改进粒子群优化技术的制作方法
【专利摘要】本发明公开了一种可控加速性能的改进粒子群优化技术，其特征在于，该方法包括如下步骤：1)群粒子搜索速度更新扩展建模；2)三种速度更新模式的定义；3)群粒子变速寻优轨迹建立；4)完全可控粒子群优化的代数学原理分析；5)完全可控粒子群优化的动力学原理分析。本发明的有益效果是，方法设计合理，使用简明，准确地求解多峰值函数的优化问题，运算效率高。
【专利说明】
-种可控加速性能的改进粒子群优化技术
技术领域
[0001] 本发明基于优化理论，在基本粒子群优化的基础上提出一种具有可控速度因子的 改进粒子群算法。在不同速度调控策略下优选具有最佳变速特性的寻优轨迹。通过采用可 控的寻优速度因子，粒子群优化性能得到明显改进，特别是采用完全控制策略，不仅可获得 较高的优化精度，而且收敛速度更快，表现出更好的综合性能。该计算方法为解决许多领域 中设及的多变量优化问题奠定了基础。
【背景技术】
[0002] 基本粒子群优化(BPSO)算法是一种全局优化算法，源于对鸟类群体捕食行为的模 拟EU，在许多学科领域中获得应用，用于解决其中面临的大量非线性、不可微和多峰值的多 变量优化问题
[0003] 然而，对于解决多峰值函数及多变量优化问题，BPSO算法在进化后期收敛速度较 慢，且有时会出现捜索粒子在全局最优解附近"振荡"的现象。为此，SHI等W提出更新粒子 速度时速度加权因子W沿线性下降的办法，结果显示算法性能未得到明显改进，计算代价却 增加较多。此外还有通过采用模糊规则动态修改W，使算法自适应地调整全局系数，试图兼 顾捜索效率和精度W。但对许多复杂的非线性优化问题，试图通过调整一个全局系数来提 高捜索精度是不够的，有必要在速度调控策略、寻优轨迹及寻优效率等方面进一步提高粒 子群算法的收敛精度W及运算效率。
【发明内容】
[0004] 本发明的目的是为了解决上述问题，开发了一种可控加速性能的改进粒子群优化 技术。
[0005] 实现上述目的本发明的技术方案为，一种可控加速性能的改进粒子群优化技术， 其特征在于，该方法包括如下步骤：
[0006] 1)群粒子捜索速度更新扩展建模；
[0007] 2)=种速度更新模式的定义；
[000引 3)群粒子变速寻优轨迹建立；
[0009] 4)完全可控粒子群优化的代数学原理分析；
[0010] 5)完全可控粒子群优化的动力学原理分析。
[0011] 所述群粒子捜索速度更新扩展建模计算式为：
[001 2] Vi - WiXvi-1+Cl X 日1 X (JDbest, i-广 pipresent, i-1) +C2 X Pi X ( gbest, i-广 pipresent, i-1),
[0013] i = l，2，l^，N.
[0014] 式中Pbest为个体极值，即单个粒子所捕捉到的最优解。gbest为全局极值，即整个粒 子群所捜索到的最优解。Ppresent为粒子的捜索速度与当前位置向玲。Wi为第i次学习的速度 加权，且有Wi G [Wmin ,Wmax]。V为粒子的捜索速度向量。Cl与C2均为方向系数，Qi与Pi为单个与 全局变化的权因子，两个因子均为学习迭代数i的单调递减函数。N为整个进化次数。
[0015] 所述=种速度更新模式的定义计算式为：
[0016]
[0017]式中CRM、PRM和CCM分别为完全随机、部分随机和完全可控速度更新模式，Wi为捜 索速度惯性加权因子，Wmax和Wmin分别为Wi的上、下界，Qi和扣分别为个体和全局变异因子， rand为[0 1]之间服从均匀分布的任一随机数，f(i)为速度变轨函数且满足f(i) G [Wmin Wmax] O
[001引所述群粒子变速寻优轨迹建立计算式为：
[0019]
[0020]
[0021]
[0022] 式中0?1'、0?1'、61^1'、阳1'和5口61'分别为完全随机、受限随机、高斯径向基、幕指数和 二次幕指数速度寻优轨迹函数，A W = Wmax-Wmin加权因子变化量，〇2、C和n均为函数响应f(i) 对变量i的散布参数且。1，N为群粒子数目。
[0023] 所述完全可控粒子群优化的代数学原理分析计算式为：
[0024]
[0025] 式中为加权因子，w、a和0分别为搜索速度惯性、个体和全局变异的加权因子，V为 捜索速度，X为粒子位置，P为全局变异向量。
[0026] 所述完全可控粒子群优化的动力学原理分析计算式为：
[0027]
[002引式中y = p-x为粒子群的全局更新位置。
【附图说明】
[0029] 图1是本发明所述一种可控加速性能的改进粒子群优化技术的流程示意图；
[0030] 图2是几种典型的寻优变速轨迹图；
[0031 ]图3简化PSO系统中一个捜索粒子的运动轨迹-从代数观点；
[0032] 图4 CCM速度更新模式下一个捜索粒子的速度与位置变化轨迹-从动力学观点；
[0033] 图5 CCM速度更新模式下简化PSO系统的固有频率与阻尼比；
[0034] 图6几个典型的PSO优化器的在线与离线性能-实验一；
[0035] 图7几个典型的PSO优化器的在线与离线性能-实验二；
[0036] 图8改进的CCM-PET基PSO优化器400个捜索粒子的两聚类分布。
【具体实施方式】
[0037] 下面结合附图对本发明进行具体描述，如图I是本发明所述一种可控加速性能的 改进粒子群优化技术的流程示意图，采用高效的捜索粒子速度调控策略和寻优轨迹，对基 本粒子群算法的寻优精度W及寻优效率进行改进。
[0038] 本技术方案从代数学和动力学角度对只有一个捜索粒子的简化粒子群系统进行 了仿真分析，并W两个有代表性的优化例子阐述改进粒子群优化方法的工作过程，其基本 原理为:通过调整寻优轨迹函数f(i)散布参数〇2、C和n可W很方便地使f(i)的响应满足其 上、下界的界值要求，如图2。即
[0039]
[0040] 实施例1
[0041 ]完全可控粒子群优化的代数学原理
[0042] 使速度更新加权因子W随着变速轨迹PET在区间[0.2 0.引内单调递减变化，基于 此对简化的BPSO系统进行分析，分析结果如图3所示。捜索粒子速度V的模与相角分别描述 于图3(a)和图3(c)，捜索粒子位置y的模与相角分别描述于图3(b)和图3(d)
[0043] 对于每个时间间隔t，w的解为：
[0044]
[0045] 从而，CCM速度更新模式对单捜索粒子的速度和位置轨迹的影响可W表达：
[0046]
[0047] 显然，BPSO中捜索粒子的速度与位置轨迹总是在一个狭窄的区域不断重复且完全 随机变化，大量的时间都浪费在非监督与低效的捜索过程之中。与BPSO不同，CCM-PET优化 器使得捜索粒子能W平滑的阻尼振荡方式在一个更宽广的解空间中运动。而且，CCM速度更 新模式使捜索粒子速度与位置随着权因子W的降低而同步成比例地改变，从而可W通过选 择合适的变速轨迹来调控粒子群优化过程。
[004引实施例2
[0049] 完全可控粒子群优化的动力学原理
[0050] 事实上，捜索速度V并非常数，而是一个不断更新的变量，运使得PSO系统的动力学 特性变得更加复杂。在CCM速度更新模式下捜索粒子的速度和位置的模轨迹如图4所示，固 有频率和阻尼比的模轨迹如图5所示。
[0051] 只考虑t = 2的情况，二阶运动微分方程为：
[0化2]
[0053] 运动方程的响应解为：
[0化4]
[0化5]
[0化6]
[0化7]
[0化引
[0059] 由图4和图5可见，变速轨迹的瞬时范围越大，捜索粒子的振荡范围就越宽，振荡持 续时间就越长。CCM基粒子群优化器所具有的运种动力学特性，对于快速侦测并俘获全局最 优解十分重要。
[0060] 实施例3
[0061] 基于改进的PSO优化器的多变量优化
[0062] 实验一中几个典型的PSO优化器的性能如图6所示。
[0063] 实验一优化问题为：
[0064]
[0065] 整个进化迭代次数N为500次，Wmin和Wmax分别为0.2和0.8。捜索粒子群为只有30个 粒子的小样本群，初始的最优历史解与最佳适配度均设置为无穷，全局最优解的初始值设 为[0 0]，单个与全局方向值均设为2,最大速度为0.3。
[0066] 在小样本情况（只有30个捜索粒子），所有10次测试中BPSO从未捜索到全局最优 解，平均运算代价为274次迭代。正如所期望的那样，几乎所有的CCM基改进PSO优化器(除了 CCM-SPET)性能都超越了 BPSO，无论是在优化精度还是运算效率方面。如表1.
[0067] 表1几种典型的PSO基优化器的优化结果一一实验一
[006引
[0069] 注：'S'一一成功；'F'一一失败；'S~'一一近似成功(只要进一步提高进化次数N)
[0070] 实验二优化问题为：
[0071]
[0072] 几种典型的PSO基优化器的优化性能如表2所示。
[0073] 表1几种典型的PSO基优化器的优化结果一一实验二
[0074]
[0075] 实验二中几个典型的PSO优化器的性能如图7所示。
[0076] 在两次优化实验中，可W注意到一个十分有趣的现象，即每次当CCM基优化器成功 俘获全局最优解后，捜索粒子总是收敛为两个明显的聚类，即成功者（与全局最优解重合） 和失败者（明显偏离全局最优解），如图8所示。
[0077]
[0078] 由此，可W为进一步改进优化方法提供有益启示：
[0079] (a)在进化期间，及时停止对那些注定失败的捜索粒子的速度和位置的更新，从而 节约运算代价。
[0080] (b)对于优化失败的情况，通过改变失败粒子的速度更新模式和变速轨迹进行二 次优化捜索。
[0081 ] 参考文献
[0082] [1]KE順抓Y J,EBERHART R C.Particle swarm optimization.Proc.of the IE邸 Conf.on 化ural 化tworks,IV.F*e;rth:IE邸 Press,1995:l 942-1 948.
[0083] [2]ANT0NI0 ISMAEL de F V,EDITE MANUELA da G P. Optimization of nonlinear constrained particle swarm. UKI 0TechnoIoginis ir Ekonominis Vystmas,2006,12(1)：30-36.
[0084] [ 3 ]高飞，童恒庆.基于改进粒子群优化算法的混浊系统参数估计方法.物理学报， 2006,55(2):577-582.
[0085] [4]SHI Y，EBERHART R C.A modified swarm optimizer.Proc.Of the IE趾 Int'I Conf.of Evolutionary Computation.Piscataway: IEEE Press,1998：69-73.
[0086] [ 5 ] SHI Y , EBERHART R C.Fuzzy adaptive particle swarm optimization.Proc.of the Congress on Evolutionary Computation.Piscataway: IE 邸，2001:101-106.
[0087] 上述技术方案仅体现了本发明技术方案的优选技术方案，本技术领域的技术人员 对其中某些部分所可能做出的一些变动均体现了本发明的原理，属于本发明的保护范围之 内。
【主权项】
1. 一种可控加速性能的改进粒子群优化技术，其特征在于，该方法包括如下步骤： 1) 群粒子搜索速度更新扩展建模； 2) 三种速度更新模式的定义； 3) 群粒子变速寻优轨迹建立； 4) 完全可控粒子群优化的代数学原理分析； 5) 完全可控粒子群优化的动力学原理分析。2. 根据权利要求1所述的一种可控加速性能的改进粒子群优化技术，其特征在于，所述 群粒子搜索速度更新扩展建模计算式为： Vi - WiXVi-l + ClXQiX ( pbest, I-I-Ppresent, 1-1)+02 Χβ? X ( gbest, I-I-Ppresent, i-1 ) , ? - I , 2 , L , Ν. 式中Pbest为个体极值，即单个粒子所捕捉到的最优解。gbest为全局极值，即整个粒子群 所搜索到的最优解。Ppr^nt为粒子的搜索速度与当前位置向玲。Wi为第i次学习的速度加权， 且有WiG [Wmin,Wmax] ^为粒子的搜索速度向量。Cl与C2均为方向系数，CIi与为单个与全局变 化的权因子，两个因子均为学习迭代数i的单调递减函数。N为整个进化次数。3. 根据权利要求1所述的一种可控加速性能的改进粒子群优化技术，其特征在于，所述 三种速度更新模式的定义计算式为：式中CRM、PRM和CCM分别为完全随机、部分随机和完全可控速度更新模式，Wi为搜索速度 惯性加权因子，Wm^Pwmin分别为Wi的上、下界，CXi和01分别为个体和全局变异因子，rand为[0 1]之间服从均匀分布的任一随机数，f(i)为速度变轨函数且满足f(i) e [wmin wmax]。4. 根据权利要求1所述的一种可控加速性能的改进粒子群优化技术，其特征在于，所述 群粒子变速寻优轨迹建立计算式为：式中CRT、LRT、GRBT、PET和SPET分别为完全随机、受限随机、高斯径向基、幂指数和二次 幂指数速度寻优轨迹函数，Aw = Wmax-Wmin加权因子变化量，σ2、ξ和η均为函数响应f(i)对变 量i的散布参数且ξ>1，N为群粒子数目。5. 根据权利要求1所述的一种可控加速性能的改进粒子群优化技术，其特征在于，所述 完全可控粒子群优化的代数学原理分析计算式为：式中为加权因于，w、a和β分别为搜累速度惯性、个体和全局变异的加权因子，ν为搜索 速度，X为粒子位置，P为全局变异向量。6.根据权利要求1所述的一种可控加速性能的改进粒子群优化技术，其特征在于，所述 完全可控粒子群优化的动力学原理分析计算式为：式中y = P_x为粒子群的全局更新位置。
【文档编号】G06N3/00GK106022462SQ201610301518
【公开日】2016年10月12日
【申请日】2016年5月6日
【发明人】焦卫东, 杨志强
【申请人】浙江师范大学