一种基于Sine映射和分段式并行连续禁忌搜索数据处理方法技术

本发明专利技术涉及一种基于Sine映射和分段式并行连续禁忌搜索数据处理方法，包括以下步骤：S1：利用Sline映射生成初始寻优候选成员；S2：采用分段式搜索生成领域及领域候选解；S3：以禁忌准则比较判断每次并行框架下迭代的目标函数适应度值；S4：当寻优迭代次数大于全局控制迭代数时，停止迭代输出CNOP最优解，否则跳至步骤S2继续寻优。与现有技术相比，本发明专利技术基于避免使用伴随模式的考虑，利用Sine映射、分阶段寻优以及并行技术改进了连续禁忌搜索算法，提高CNOP求解的效率，并扩展了其应用领域。

A Data Processing Method Based on Sine Mapping and Piecewise Parallel Continuous Tabu Search

The invention relates to a data processing method based on Sine mapping and piecewise parallel continuous tabu search, which includes the following steps: S1: generating initial candidate members by Sline mapping; S2: generating domain and domain candidate solutions by piecewise search; S3: judging the fitness value of the objective function of each iteration under the parallel framework by Tabu criteria; S4: when the number of iterations is large. When the number of iterations is controlled globally, the output of CNOP optimal solution is stopped, otherwise the optimization can be continued by jumping to 2. Compared with the prior art, the present invention improves the continuous tabu search algorithm by using Sine mapping, stage-by-stage optimization and parallel technology, improves the efficiency of CNOP solution and expands its application field based on avoiding the use of adjoint pattern.

【技术实现步骤摘要】
一种基于Sine映射和分段式并行连续禁忌搜索数据处理方法
本专利技术涉及计算机科学与气象科学的交叉学科研究领域，用于求解条件非线性最优扰动CNOP(ConditionalNonlinearOptimalPerturbation，CNOP)的优化求解，尤其是涉及一种基于Sine映射和分段式并行连续禁忌搜索数据处理方法(ParallelContinuousTabuSearchAlgorithmwithSineMapsandStagedStrategy，CTS-SS)。
技术介绍
CNOP一般应用于气象科学中进行天气或气候的可预报性研究，其求解最常用的方法是谱投影梯度方法。该方法基于气象数值模式的伴随模式提供梯度信息进行寻优，但伴随模式求解梯度的计算量巨大、耗时较长，且谱投影梯度方法无法并行。此外，不是所有的数值模式都开发了相应的伴随模式，且伴随模式的开发工程量巨大，尤其对于很多复杂的大模式来说，该过程可能会消耗数年时间。本专利技术基于避免使用伴随模式的考虑，利用Sine映射、分阶段寻优以及并行技术改进了连续禁忌搜索算法，提高CNOP求解的效率，并扩展了其应用领域。
技术实现思路
本专利技术的目的就是为了克服上述现有技术存在的缺陷而提供一种基于Sine映射和分段式并行连续禁忌搜索数据处理方法，本专利技术的CTS-SS方法主要包括三方面改进：Sine映射优化初始寻优候选解、分阶段寻优、并行化框架。本专利技术的目的可以通过以下技术方案来实现：一种基于Sine映射和分段式并行连续禁忌搜索数据处理方法，包括以下步骤：S1：利用Sline映射生成初始寻优候选成员；S2：采用分段式搜索生成领域及领域候选解；S3：以禁忌准则比较判断每次并行框架下迭代的目标函数适应度值；S4：当寻优迭代次数大于全局控制迭代数时，停止迭代输出CNOP最优解，否则跳至步骤S2继续寻优。进一步地，所述步骤S1包括以下分步骤：S11：在(-1，1)区间内随机生成初始个体X0，随后对其进行迭代K次的混沌计算；S12：通过K次迭代将初始个体X0变成Xk；S13：将迭代后的值赋值给初始候选寻优个体Ci；S14：用Sine映射方法生成了pop个初始候选寻优个体，从中选出一个具有最佳适应度函数值的寻优个体C作为初始寻优候选解。进一步地，所述步骤S1的过程描述公式为：x0,j＝rand(-1,1),j＝1,2,...,nxk+1,j＝sin(πxk,j),xk,j∈(-1,1),k＝0,1,2...KCi.j＝xK,j,i＝1,2,...,pop。C＝min{C1,C2,...,Cpop}i＝1,2,...,pop进一步地，所述步骤S2包括以下分步骤：S21：设定第一领域半径值和第一领域候选解个数进行第一阶段全局搜索；S22：设定第二领域半径值进行第二阶段局部搜索。进一步地，所述步骤S21，其描述公式为：式中，R1为第一阶段领域半径，Si1为第一阶段候选解，j＝1,2,...,n，i1＝1,2,...,m1，m1为第一阶段候选解个数。进一步地，所述步骤S22，其描述公式为：式中，R2为第二阶段领域半径，Si2为第二阶段候选解，j＝1,2,...,n，i2＝1,2,...,m2，m2为第二阶段候选解个数。进一步地，所述步骤S3包括以下分步骤：S31：针对第一搜索阶段，如果Si1不在禁忌表内，判断适应度函数值f(Si1)-best_value>＝ε是否为真，若为真，则将当前最优解C列入禁忌表，同时更新best_value＝f(Si1)，C＝Si1，其中best_value是当前最优解C的适应度函数值，ε是禁忌判定参数，第一阶段共运行maxcycle1次。S32：针对第二搜索阶段，如果Si2不在禁忌表内，判断适应度函数值f(Si2)-best_value>＝ε是否为真，若为真，则将当前最优解C列入禁忌表，同时更新best_value＝f(Si2)，C＝Si2，其中best_value是当前最优解C的适应度函数值，ε是禁忌判定参数，第一阶段共运行maxcycle2次。与现有技术相比，本专利技术具有以下优点：(1)求解效率高，速度快，本专利技术基于避免使用伴随模式的考虑，利用Sine映射、分阶段寻优以及并行技术改进了连续禁忌搜索算法，提高CNOP求解的效率，并扩展了其应用领域。(2)精确度高，稳定性好，本专利技术通过本专利技术在气象数值模式中进行条件非线性最优扰动(CNOP)的优化求解。首先使用Sine映射生成初始寻优候选解，设置两次分阶段寻优过程中的领域半径及其候选解个数，根据当前的硬件资源设置并行计算核心数，再进行主进程和子进程的分配，完成基于Sine映射和分段式策略的并行连续禁忌搜索算法(CTS-SS)高效求解CNOP，精确度高，稳定性好。附图说明图1为本专利技术的方法流程示意图；图2为本专利技术的分段式流程示意图。具体实施方式下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本专利技术的一部分实施例，而不是全部实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都应属于本专利技术保护的范围。实施例如附图1所示，左侧是CTS-SS算法的整体流程，由Master进程(主进程控制)：包括Sine映射生成较优初始成员(初始寻优候选解)、生成邻域及邻域候选解(包括第一阶段和第二阶段)、并行计算各候选解的适应度函数值(目标函数适应度值)、判断是否程序是否终止，若是，输出CNOP；否则更新成员(当前最优解)，再进行阶段性寻优。附图1的右侧为Master-Slaver并行结构计算各候选解适应度值的示意图，两次阶段性寻优的计算过程均如图中所示，区别在于每个Slaver进程所计算的候选解的个数。附图1中Sine映射生成初始寻优候选解的过程如下：x0,j＝rand(-1,1),j＝1,2,...,nxk+1,j＝sin(πxk,j),xk,j∈(-1,1),k＝0,1,2...KCi.j＝xK,j,i＝1,2,...,popC＝min{C1,C2,...,Cpop}i＝1,2,...,pop以上四个式子对应具体步骤过程如下：1：在(-1，1)区间内随机生成初始个体X0，随后对其进行迭代K次的混沌计算；2：通过K次迭代将初始个体X0变成Xk；3：将迭代后的值赋值给初始候选寻优个体Ci；4：用Sine映射方法生成了pop个初始候选寻优个体，从中选出一个具有最佳适应度函数值的寻优个体C作为初始寻优候选解。附图2展示了CTS-SS算法的分段式搜索策略，即将原连续禁忌搜索算法(CTS)的计算过程分为两个阶段：第一阶段选用较大的邻域半径值并设定较小的邻域候选解个数，增大全局搜索力度，此阶段的处理按照下式：式中，R1为第一阶段领域半径，Si1为第一阶段候选解，j＝1,2,...,n，i1＝1,2,...,m1，m1为第一阶段候选解个数。在第一阶段邻域半径R1生成m1候选解个数{S1，S2，…，Si1，…，Sm1}，当满足第一个条件时，说明Si1在半径为R1的领域内，无需调整；当满足第二个条件时，将Si1的每一维调整至半径为R1的领域内。计算Si1的适应度函数值f(Si1)，如果Si1不在禁忌表内，判断f(Si1)-best_v本文档来自技高网...

【技术保护点】
1.一种基于Sine映射和分段式并行连续禁忌搜索数据处理方法，其特征在于，包括以下步骤：S1：利用Sline映射生成初始寻优候选成员；S2：采用分段式搜索生成领域及领域候选解；S3：以禁忌准则比较判断每次并行框架下迭代的目标函数适应度值；S4：当寻优迭代次数大于全局控制迭代数时，停止迭代输出CNOP最优解，否则跳至步骤S2继续寻优。

【技术特征摘要】
1.一种基于Sine映射和分段式并行连续禁忌搜索数据处理方法，其特征在于，包括以下步骤：S1：利用Sline映射生成初始寻优候选成员；S2：采用分段式搜索生成领域及领域候选解；S3：以禁忌准则比较判断每次并行框架下迭代的目标函数适应度值；S4：当寻优迭代次数大于全局控制迭代数时，停止迭代输出CNOP最优解，否则跳至步骤S2继续寻优。2.根据权利要求1所述的一种基于Sine映射和分段式并行连续禁忌搜索数据处理方法，其特征在于，所述步骤S1包括以下分步骤：S11：在(-1，1)区间内随机生成初始个体X0，随后对其进行迭代K次的混沌计算；S12：通过K次迭代将初始个体X0变成Xk；S13：将迭代后的值赋值给初始候选寻优个体Ci；S14：用Sine映射方法生成了pop个初始候选寻优个体，从中选出一个具有最佳适应度函数值的寻优个体C作为初始寻优候选解。3.根据权利要求1所述的一种基于Sine映射和分段式并行连续禁忌搜索数据处理方法，所述步骤S1的过程描述公式为：4.根据权利要求1所述的一种基于Sine映射和分段式并行连续禁忌搜索数据处理方法，所述步骤S2包括以下分步骤：S21：设定第一领域半径值和第一领域候选解个数进行第一阶段全局搜索；S22：设定第二领域半径值进行第二阶段局部搜索。5.根据权利要求1所述的一种基于Sine映射和分段式并行连续禁忌搜索数据处理...

【专利技术属性】
技术研发人员：穆斌，张琳琳，袁时金，钱一闻，
申请(专利权)人：同济大学，
类型：发明
国别省市：上海,31