基于全局和局部代理辅助的粒子群算法的优化方法及应用

本发明涉及工程技术和进化计算，特别涉及一种基于全局和局部代理辅助的粒子群优化算法的优化方法及应用。

背景技术：

随着现代工业和科技的发展，电磁(electromagnetic,em)设备在航天卫星、交通运输和医疗保健等领域中发挥着重要的作用，em设备的设计直接影响到其设备的性能，因此，对em设备进行优化设计具有很大的现实意义与必要性。此外，em设备设计优化问题一直以来都是计算电磁学的研究热点，如何设计出一种有效的方法来优化em设备的效率是值得深思的问题。

目前em设备的优化方法有两种：第一种是通过基于数值计算的em模拟(例如有限元分析)来设计出高性能的em；第二种是通过进化算法对em进行优化设计，由于进化算法能适应不同环境下的问题以及在大多数情况下都能得到令人满意的有效解，其已广泛应用于em设备的优化设计中。然而，为了获得em设备的最佳性能估计，这些优化方法均需要花费大量的模拟仿真，对于设计复杂的em设备，应用是不可能的，因此，在合理的优化时间来优化出em设备的最佳设计是一大挑战。

由于em设备的设计直接影响到其性能，em设备的优化设计显得十分重要。在现实生活中，随着em设备的设计规范与性能要求日益严格，em设备的设计优化也需要大量的仿真模拟才能获得最佳设计，例如一个典型的计算电磁模拟可能需要30分钟，这意味着该优化需要花费大量的计算开销。

一种能够加快算法收敛速度，将搜索与开发达到平衡的优化算法亟待开发。

技术实现要素：

针对现有的优化算法计算开销大的问题，本发明提供一种基于全局和局部代理辅助的粒子群算法的优化方法及应用，能够有效解决上述问题。

为了解决上述技术问题，本发明采用的技术手段是：

一种基于全局和局部代理辅助的粒子群算法的优化方法，包括：

对整个种群进行初始化，得到初始化种群；

使用目标函数对所述初始化种群中每个个体进行真实评估，并将所述初始化种群和所述初始化种群中每个个体的适应度值直接保存到数据库db中；

使用全局搜索与局部搜索交替进行的方式，不断通过更新迭代从而把最优个体保存下来；

迭代结束后，将现存的最优个体来作为最优组合。

所述“对整个种群进行初始化”，包括：

根据设计空间的范围使用latinhypercube采样方法对整个种群进行初始化的步骤。

所述“全局搜索”方式，包括：

通过集成模型对所述新种群中的新样本进行预测，获得每个所述新样本的预测值；

选择预测值最小的个体进行真实适应度评估，并将所述预测值最小的个体以及所述预测值最小的个体对应的适应度值保存到所述数据库db中；

判断所述预测值最小的个体是否比目前的最佳个体更优：

若，更优，则替换掉所述目前的最佳个体，直到达到停止条件；

否则切换到所述局部搜索。

所述“集成模型”，包括：

使用所述数据库db中的全部样本来分别构建：gp模型、多项式回归模型以及rbf模型；

通过加权求和方法来组合所述gp模型、所述多项式回归模型以及所述rbf模型，构成所述集成模型。

所述集成模型，包括：

权重计算公式如下所示：

其中εj为第j个代理模型的预测误差度量，wj是第j个代理模型的权重值；代理模型的预测误差度量值为均方根误差值。

每个代理模型的所述均方根误差值通过下式获得：

其中，均方根误差值用符号ε来表示；则第j个代理模型的均方根误差值值为εj；y(xi)为xi的真实适应度值，为第j个代理模型对xi的预测值；m为代理模型的训练样本个数。

所述“获得每个所述新样本的预测值”，包括：

所述预测值通过下式获得：

其中，是所述集成模型的最终预测值，是第j个代理模型对新样本x^*的预测值，wj是第j个代理模型的权重，k为代理模型的个数。

所述“局部搜索”方式，包括：

使用粒子群优化算法来对整个种群进行进化，得到新种群；

从所述数据库db中选择m个最优的样本作为训练样本集；

使用所述训练样本集来构建gp模型；

使用gp模型来预测出新种群个体的lcb值；

使用粒子群优化算法来对具有最小lcb预测值的个体进行真实适应度评估，并将所述具有最小lcb预测值的个体与所述具有最小lcb预测值的个体对应的适应度值保存到数据库db中；

判断该个体是否比目前的最佳个体更优，若更优，则替换掉最佳个体，直到达到停止条件，否则切换到全局搜索。

所述“使用gp模型来预测出新种群个体的lcb值”，包括：

预测分布对于f(x)的lcb值定义为：

其中ω为常数；s²(x)预测的不确定性；x为新种群中的个体；n为正态分布。

一种如上所述的基于全局和局部代理辅助的粒子群算法的优化方法在em设备进行优化设计方向上的应用。

本发明的有益效果至少包括：

(1)本发明采用具有自适应惯性权重策略的粒子群优化算法作为优化器。由于标准的粒子群优化算法在高维优化问题上存在容易陷入局部最优解的缺点，本发明为了让算法的全局搜索能力和局部搜索能力能达到平衡，因此在标准的粒子群优化算法中引入一种自适应惯性权重策略，将粒子群看做一个整体，从而调整惯性权重。

(2)在本发明的全局和局部代理辅助的粒子群优化算法中，对于em设备设计的整个搜索空间，采用全局搜索的方法，皆在快速定位有希望的区域。在全局搜索中，通过加权求和方法来组合单个代理模型从而形成一个强大的集成模型，这使得在全局搜索中针对不同情景下的em优化设计问题，通过更合理的方式为各代理模型分配合适的权重进行组合，最终作出决策。

(3)在本发明的全局和局部代理辅助的粒子群优化算法中，对于有希望区域进行局部搜索，皆在开发出比当前最佳值更优的解决方案。在局部搜索中，通过从数据库db中选取一定数量的最佳样本作为训练集来训练单个代理模型，进一步使用基于该代理模型的模型管理策略来选择合适的个体进行真实适应度评估，皆在局部区域加快算法的收敛速度，找到更佳的解决方案。

(4)采用一种模型管理策略用于对集成模型进行精准管理。在全局搜索中，先使用代理模型对新种群进行预测并筛选出适应度预测值最佳的个体进行真实适应度评估；而在局部搜索中，通过一定数量的最佳样本来训练代理模型，并使用基于gp模型的lcb准则来筛选出具有最小lcb预测值的个体，从而对其进行真实评估。

本发明所述的方法能够通过全局搜索与局部搜索相结合，有利于加快算法的收敛速度，将搜索与开发达到平衡。

附图说明

图1为本发明所述方法的流程图。

具体实施方式

下面结合附图对本发明进行进一步的说明。

本发明的核心思想是：本发明所述的方法将对em设备的优化设计问题进行全局搜索和局部搜索：在全局搜索中，先使用优化器对种群进行进化，接着通过使用数据库db中的全部样本对集成模型进行训练，然后将集成模型用于预测新种群的适应度值并筛选出适应度预测值最佳的个体进行实际函数评估并将其及对应的真实适应度值加入到数据库db中；在局部搜索中，先通过从数据库db中选定一些当前最佳的解决方案用于构建局部模型，并使用基于gp模型的最优置信下限(lowerconfidencebound,lcb)准则来选取合适个体进行实际函数评估并将其及对应的真实适应度值加入到数据库db中。

值得注意的是，如果某个搜索无法找到比当前全局优化解更优的解决方案，全局搜索和局部搜索将交替进行。通过全局搜索与局部搜索相结合，有利于加快算法的收敛速度，有利于将搜索与开发达到平衡。

具体的，如图1，本发明所述的基于全局和局部代理辅助的粒子群算法的优化方法，包括以下步骤：

s1、根据设计空间的范围使用latinhypercube采样方法对整个种群进行初始化，并使用昂贵的目标函数对种群中每个个体进行真实评估；最后将初始化种群和它们的适应度值直接保存到数据库db中；

s2、若未达到算法停止条件，则进行全局搜索：

s201.使用粒子群优化算法来对整个种群进行进化，得到新的种群；

s202.使用数据库db中的全部样本来分别构建gp模型、多项式回归模型和rbf模型，并通过加权求和方法来组合三个模型，从而构成集成模型；

权重计算公式如下所示：

其中εj为第j个代理模型的预测误差度量，wj是第j个代理模型的权重值。在本文中，用均方根误差值(rmse)来作为代理模型的预测误差度量值：

每个代理模型的rmse可以由以下公式来计算：

其中，rmse用符号ε来表示，则第j个代理模型的rmse值，y(xi)为xi的真实适应度值，为第j个代理模型对xi的预测值，m为代理模型的训练样本个数。

因此，当集成模型对新样本x^*进行预测时，最终的预测值将会由以下公式计算得出：

其中，是集成模型的最终预测值，是第j个代理模型对新样本x^*的预测值，wj是第j个代理模型的权重，k为代理模型的个数。

s203.使用集成模型对整个新种群的候选解进行适应度预测；

s204.选择预测值最小的个体进行真实适应度评估，并将它与其对应的适应度值保存到数据库db中；

s205.判断该个体是否比目前的最佳个体更优，若更优，则替换掉最佳个体，继续执行步骤1,直到达到停止条件，否则切换到局部搜索；

s3、若未达到算法停止条件，则进行局部搜索：

s301.使用粒子群优化算法来对整个种群进行进化，得到新的种群；

s302.从数据库db中选择m个最优的样本作为训练样本集a；

s303.使用训练样本集a来构建gp模型；

s304.使用gp模型来预测出新种群个体的lcb值；给出预测分布对于f(x)的lcb值可以定义为：

其中ω为常数。

s305.使用粒子群优化算法来对具有最小lcb预测值的个体进行真实适应度评估，并将它与其对应的适应度值保存到数据库db中；

s306.判断该个体是否比目前的最佳个体更优，若更优，则替换掉最佳个体，继续执行(1),直到达到停止条件，否则切换到全局搜索；

s4、当迭代结束后，将现存的最优个体来作为em设备的最优变量设计组合。

所述步骤s2和步骤s3中，通过具有自适应惯性权重的粒子群优化算法来作为种群优化器，有利于平衡算法的全局搜索能力和局部搜索能力。

所述步骤s2和步骤s3中，在全局搜索中选取集成模型预测的适应度最佳的个体进行真实适应度评估，从而提高集成模型的精度；在局部搜索中通过基于代理模型的模型管理策略来选择合适的个体进行实际函数评估，进一步加快算法的收敛速度。通过结合全局搜索和局部搜索，从而让算法达到搜索和开发互相平衡。

具体实施方式：

需要说明的是：

1.本文所述的代理模型的选取可以使用不同的机器学习模型来代理，例如人工神经网络和随机森林等等；

2.可以选取其他的智能优化算法作为优化器，例如差分进化算法和遗传算法等等；

3.可以选取其他集成模型组合策略来组合单个代理模型，例如通过stacking的方法来训练出一个人工神经网络，从而对模型进行组合；

4.可以选取其他模型管理策略来提高模型的精度，例如基于委员会查询(qbc)的模型管理策略。

以上，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易变化或替换，都属于本发明的保护范围之内。因此本发明的保护范围所述以权利要求的保护范围为准。