0 引言

许多现实世界中的优化问题都表现出动态性质，即优化的目标函数或待求解的问题会随时间而发生(随机)变化。例如，目标识别过程中，敏感组件的检测性能会受到周围环境的影响和制约;调度问题中，新工件可能随时到达，机器可能发生随机故障或逐渐磨损，资源量会随时间不断变化，原材料的性能可能会随时间发生变化;金融贸易模型中的市场环境可能会突然改变［1－2］;因此，研究适合于求解这些普遍存在于现实世界的动态优化问题(Dynamic Optimization Problem，DOP)的方法和算法有重要的现实意义。

一般来说优化问题可以分为两类:静态优化问题(Static Optimization Problem，SOP)和动态优化问题(DOP)。其中，静态优化问题有统一的定义，可以认为静态优化问题是只有一个解的单决策问题，然而对DOP目前尚且没有一致的定义。简言之，在动态优化问题中，其最优解的位置会随时间推移而改变。

求解DOP的优化算法要求能够在动态环境中搜索到最优解［3－4］。尽管传统演化算法(Evolutionary Algorithm，EA)如粒子群算法、差分进化算法、人工免疫算法等在求解静态优化问题上取得了很好的效果，但求解DOP需要对传统EA进行相关的改进。

目前，动态优化算法主要可以分为环境变化后增加多样性的方法、运行过程中始终保持多样性的方法、基于记忆机制的方法、多种群方法和基于预测机制方法5类［5］。其中，多种群方法主要应用在一类具有多个局部最优(多峰)函数的动态优化问题上，而多峰函数的优化在现实中较为普遍。

对此国内外许多学者将多种群方法与经典演化算法相结合后应用到对动态优化问题的求解上。如:Yang等［6］使用了分层聚类的方法把种群分成多个子种群，该方法是子种群初始粒子依据在适应值曲面的分布自动产生，而不是依赖于像kmeans聚类方法中的参数k;Brest等［2］提出了自适应参数的多种群差分进化(Self-Adaptive Differential Evolution，jDE)算法，控制参数F和Cr能够有效地自适应;姜立强等［7］提出改进差分进化(Modified Differential Evolution，MDE)算法，该算法将种群分为侦测和搜索两个种群，但该算法过于简单，在求解多峰函数时极易陷入局部最优;陈健等［8］提出多种群骨干粒子群优化 (Multi-swarms Bare Bones Particle Swarm Optimization，MBBPSO)算法，该算法通过设置环境勘探粒子及时检测环境的变化，环境变化后，利用上一个环境搜索的信息初始化新的种群，当种群陷入停滞时，采用新的进化方程以加强粒子的活性和使用多种群策略维持种群的多样性，但该算法在求解大规模动态优化问题上表现不佳;Ali等［9］提出基于成功率的自适应的种群动态衰减的部落差分进化算法(adaptive successbased Tribal DE algorithm with dynamic population Reduction，sTDE-dR)算法，该算法采用自适应的策略来控制F、CR参数，并提出多种群竞争策略。算法中种群的整体规模逐渐减小，各子种群通过比较平均误差来确定各子种群在下一代保存多少，以此来保证自适应的最优策略在最优子种群中引导搜索最优解和更充分地利用有限的评价代价。

虽然国内外学者对改进EA求解DOP作出了很多工作，但是在平衡算法的搜索操作和寻优操作方面，在搜索过程中维持种群的多样性方面，在算法跳出局部最优方面以及在快速追踪移动的最优解方面依然有很多不足。

其三，关于统计与概率的案例研究虽多，但重复性研究较多，活动设计趋于形式化，未能真正落实数据分析观念.基于此，研究者尤其是一线教师联合高校等研究人员应注重积累和开发优秀的教学案例，挖掘教材中蕴涵的观念，充分利用素材体现观念.优秀的案例应以现实生活为背景，使用真实的数据；利用问题驱动统计活动，体会统计思维；使用信息技术，促进学生概念性理解；以学生为主体，给予适当的反馈.

本文在已有研究的基础上，提出一种求解动态优化问题的多种群竞争差分进化算法(Differential Evolution algorithm with Competitive Strategy based on multi-population，DECS)。该算法将一个子种群作为侦测种群，采用新的侦测方法;其余多个子种群作为搜索种群，通过各搜索种群相互竞争来增强种群多样性，且更充分地利用了有限的评价代价;在搜索种群搜索过程中，保证在一个局部最优邻域内有且仅有一个搜索种群进行寻优操作。数值算例表明了该方法的有效性和可行性。

各级各有关部门始终把人民群众的生命安全放在首位，及时发布洪水预警，全面排查危险地区群众，提前转移受洪水威胁人员，有效保障了群众生命安全。灾区各级政府多方筹集资金，加大安置和救济工作力度，保障了受灾群众基本生活。针对房屋损毁多、洪水时间长的实际困难，黑龙江省强化受灾群众集中安置与投亲靠友安置，全力推进房屋修缮重建，确保群众安全温暖过冬，帮助群众尽快返迁。由于“以人为本”始终贯穿于各个工作环节，大江大河洪水未造成一人死亡，灾区社会和谐稳定。

1 相关工作

1．1 动态优化问题的定义和描述

在文献［10－12］中，DOP被定义为在一段时间内的一系列SOP，即每个动态函数由一个基本静态函数和由一组动态规则应用到该基函数所获得的动态函数序列组成。定义如下:

本系统完成了在线考试系统的诸如学生考试，老师出题、出卷、管理考试等基本功能，并加入了考试分析功能、公告功能等拓展功能，大大方便了教学的应用。

定义1 假设ψ为搜索空间，向量x∈ψ，时间t∈T，一个动态函数DFt描述如下:

 

其中:fb0(x)是带有k个可变化特征的基准静态函数;gt是在t－1时刻动态适应度函数的参数函数;cht是在t时刻基函数所有可能的改变集合;st是t时刻特征的改变程度，最后求得t时刻新的评价值。

定义2 函数gt可以被描述为:

 

其中:为在t时刻的m×n种可能的改变方式集合。这个集合中，cic表示哪一个函数特征改变，FLj表示环境函数DFt－1根据icth决定哪一部分的环境改变。

在静态环境下，传统EA中种群个体不断迭代而逐渐集中在搜索空间的一个固定解或者搜索空间的一个有限区域，正常而不是过早地收敛对于传统EA处理静态优化问题是必需的，但是，对于动态优化问题而言，一旦收敛，那么当新的环境到来时，EA就失去了探索新的区域所需要的多样性，从而不能追踪到在新的环境下已经变化了的最优解，这种收敛趋势限制了EA算法的性能;因此在动态优化问题的求解过程中，对不断变化的最优解的快速跟踪甚至比找到最优解本身更有意义［5］，求解DOP的主要挑战在于维持种群多样性和追踪移动的最优解。

而基于群体实参数的差分进化(DE)算法［13－14］是目前最好的随机优化算法之一，该方法原理简单、受控参数少、鲁棒性强，通过启发式搜索策略具有自组织性、自适应性及并行性等特点，并且不要求目标函数连续、可微等信息，具有极好的全局搜索能力［15］。基于此，本文通过改进差分进化算法来求解DOP。

1．2 量子个体

Blackwell等［16］认为所有个体不一定要全部按照相同的规则产生，因此他们提出一种类似于量子机制的量子粒子规则，在种群中最优个体位置的邻域内产生新的个体。

每当环境改变时，记录本次变化期间的离线误差:

 

这种方法在rcloud范围内，会有较高的可能性生成靠近最优值的点。此外，随着维度D的增加，可能性也会增加。

2 多种群竞争的差分进化算法

2．1 侦测环境变化

如何及时监测环境变化是任何进化算法求解DOP必须解决的首要问题，一般是通过设置监测点来实现。目前常用的监测点选取方法是随机选取搜索种群中任意个体或同时选择搜索种群中最优个体和次最优个体［17］，当发现监测点的评价值改变时则判断为环境改变，但这两种方法都不能及时感知环境变化。

DECS专门设定一个种群作为侦测种群，该种群中任何一个个体的评价值发生变化，都将判定为环境发生改变。

步骤7 竞争算法:每迭代m次，计算P1，P2，…，Pn种群中最优个体评价值f(best1)，f(best2)，…，f(bestn);保留f(best1)，f(best2)，…，f(bestn)中最优评价值所对应的种群，其他种群重新初始化;并对保留的唯一种群的下一代执行量子个体生成机制。

2．2 排除规则

使用多个种群来求解动态优化问题需要各个子种群能有效地分开和避免寻找相同的区域，确保没有任意两个种群追踪同一个局部最优位置的邻域(峰)［17］。

医院可以加强对于妇科护理人员的思想教育,让他们能够主动的去与患者进行交流沟通,意识到与患者交流沟通的重要性,让妇科患者了解到妇科护理的重要性和难度所在,积极的配合妇科护理人员完成护理工作,认识到配合护理工作完成的重要性,让患者和婴儿都能够获得更好的护理,避免不必要的安全隐患的发生。

每一次迭代，将各搜索种群中评价值最优个体依次进行比较，如果任意两个搜索种群中最优个体的欧氏距离小于预先定义的阈值rexcl，则有可能两个搜索种群在同一个局部最优位置的邻域寻优。那么，将这两个搜索种群中最优个体的评价值进行比较，将较好评价值个体所在的搜索种群留下，另一个搜索种群重新初始化。

对排斥半径rexcl的取值按照以下的经验规则产生:假设所有的峰均匀分布在整个搜索空间，那么

 

其中:X为搜索范围，peaks为在环境中峰的个数，D为搜索向量的维度。

算法在同一个峰只保留一个最优种群，将其他种群重新初始化，以便充分利用有限的评价代价和提高种群的多样性。

2．3 种群竞争机制

每迭代m次，就将n个搜索种群中最优个体的评价值进行比较，保留较优个体所在的搜索种群，剩余n－1个搜索种群重新初始化，在不断迭代中，多次竞争。

因为对动态优化问题的求解，总的评价次数是预先设定的，而所有搜索种群在独立寻优的过程中，各搜索种群的每一次迭代都需要花费一定的评价代价，如果陷入局部最优的种群继续迭代，会造成评价代价的浪费。但是将竞争失败的种群重新初始化后，就可以将其花费的评价代价用在对环境的搜索上。一旦其找到更优的局部最优区域，那么，其他搜索种群又将重新初始化。

这种竞争机制保持一个搜索种群进行寻优操作，n－1个搜索种群进行环境探索操作，各搜索种群通过竞争决定是寻优操作或探索操作。以此确保在搜索最优解的同时增加了种群的多样性，从而有效地解决了种群容易陷入局部最优的问题。

作文教学是语文教学不可分割的一部分，而写日记又是提高学生写作能力的一条很好的途径。如果让学生在日记中根据内容配画出相应的图画，会保持学生写日记的新鲜感，让他们一直有盎然的兴趣写下去。引导学生绘画、写日记，达到以画促写的效果。

2．4 DECS 描述

基于以上描述，本文提出的多种群竞争差分进化算法的主要步骤如下:

步骤1 初始化。设置算法维度D，设置算法总评价次数，种群初始化范围［Xmin，Xmax］，峰的数量 peaks，初始化侦测种群P0，计算P0的评价值向量f(P0)，初始化各搜索种群P1，P2，…，Pn，算法迭代次数 iter← 0。

步骤2 侦测环境变化算法。

1)若环境发生维度变化，则更新维度参数D，并将所有搜索种群 P1，P2，…，Pn重新初始化，iter←0，跳转步骤8;

2)若侦测种群的评价值向量中任何一个值发生变化，则以变化后的评价值向量更新变化前的评价值向量，并将所有搜索种群 P1，P2，…，Pn重新初始化，iter← 0，跳转步骤 8。

当主程序检测到铆接力和位移大于设定的初始阈值时，则铆接位置识别程序启动。首先采集一幅图像并保存，然后将图像送入在线识别子程序中，输出铆接位置。

步骤3 对n个搜索种群执行变异操作。

步骤4 对n个搜索种群执行交叉操作。

水库始建于1955年，建成于1973年，运行50年来，大坝出现了渗漏、滑坡等安全问题，经2007年大坝安全鉴定，为 “三类坝”，并于2010年年初起对水库进行除险加固，重点对大坝进行了综合防渗处理。工程于2011年10月完工，同年12月通过了验收。实践证明防渗效果显著。

步骤5 对n个搜索种群执行选择操作。

步骤6 排除算法:

1)计算最小排除距离 rexcl的值;计算 P1，P2，…，Pn种群的评价值向量，并找出其中最优的个体besti(i=1，2，…，n);

2)计算种群中最优个体的欧氏距离:rj=d(besti，best((i+1)mod n))，j=1，2，…，n。如果rj＜ rexcel，则计算f(besti)与f(best((i+1)mod n))，保留 besti，best((i+1)mod n)中较优评价值所在搜索种群，另一个搜索种群重新初始化。

而针对环境中维度改变的情况，只要侦测种群发现自身维度与环境维度不一致，则判定为环境发生改变。

接种了一些进口疫苗，比如水痘，如果还是发生了要预防的疾病，可以退还疫苗费用。如果发生不良反应，可以赔偿治疗费用。注意保留接种记录和看病的发票。

  

图1 多种群竞争差分进化算法的流程Fig．1 Flowchart of multi-population-based competitive differential evolution algorithm

3 实验设置与测试函数

仿真实验环境为 Windows 10系统，其中 CPU为 Inter Core i5-4590/3．30 GHz，4核，内存为 8 GB，实验使用 Matlab进行编码。

对于49个动态问题中每一个动态问题，独立运行5次，总评价次数为7．5E06，每评价1．5E05次环境改变，因此单个问题每次运行时环境改变10次。

 

表1 测试基准函数Tab．1 Details of the basic benchmark functions

  

函数名函数范围Sphere f(x)=∑n x2i［－100，100］i=1［－5，5］Weierstrass f(x)=∑n Rastrigin f(x)=∑n(x2i－10 cos(2πxi)+10)i=1 kmax kmax［akcos(πbk)］;a=0．5，b=0．3，kmax=20 i=1(∑k=0［akcos(2πbk(xi+0．5))］) － n∑k=0［－0．5，0．5］Griewank f(x)= 1 4000∑n cos(xi/(i)1/2)+1 ［－100，100］i=1(xi)2－∏n i=1 Ackley f(x)=－20 exp( －0．2 1槡n∑n x2 i=1 i)－exp(1n∑n i=1 cos(2πxi))+20+e ［－32，32］

对GDBG benchmark中每一个测试函数进行7种改变方式，包括小步改变(T1)、大步改变(T2)、随机改变(T3)、混沌改变(T4)、周期性改变(T5)、带噪声的周期性变化(T6)和维度改变(T7)，以此模拟了49个标准动态变化问题。其中仅F1是求最大值问题，且需要分别求解10个峰和50个峰两种情况;F2～F6函数都是求最小值问题。

为了验证算法的性能，这里选取了IEEE CEC 2009提出的标准动态优化问题集 GDBG benchmark进行测试，此benchmark构建了峰的位置、高度、宽度、维度会变化的动态环境，具体函数形式见文献［18］，并将计算结果与复位粒子群优化(restart Particle Swarm Optimization，rPSO)算法［19］(经典算法PSO求解动态环境问题的改进算法，一旦环境变化，种群将重新初始化)、人工免疫算法(Artificial Immune Network for Dynamic optimization，Dopt-aiNet)算法［19］和 MDE 算法［7］进行比较。GDBG中的测试函数见表1。其中F1函数为Rotation peak函数，具体形式见文献［18］，F2函数是Sphere函数的混合函数，F3函数是Rastrigin函数的混合函数，F4函数是Griewank函数的混合函数，F5函数是Ackley函数的混合函数，F6函数是表1中所有函数的混合函数。

在全局最优位置Vg处，以半径rcloud的球体，种群下一代产生规则如下:

 

其中:f(xbest(t))为当前环境中理论最优评估值，f(x*(t))为当前算法求解的实际最优评估值。

当一个动态问题达到总评价次数时，计算算法的平均离线误差期望(Average mean，Avg_mean)、平均离线误差标准差(STandard Deviation，STD)。计算公式见文献［18］。

DECS参数如表2所示，rPSO算法、Dopt-aiNet算法和MDE算法的参数设置分别见文献［7，18－19］。

 

表2 DECS参数Tab．2 Parameter settings for DECS

  

0．5个体维度D 10 排斥距离rexcl 0．1搜索种群数量n 3 量子个体产生半径rcloud 1交叉概率CR 0．5 竞争代数m参数 值 参数 值种群规模NP 10 缩放因子F 200

4 实验结果

表3～9给出了三种算法在GDBG benchmark上的测试结果，包括平均离线误差期望(Avg_mean)、平均离线误差标准差(STD)，表中较优实验结果用加粗字体表示。

比如日本。为了提高全民素质，日本针对不同群体设立了诸多阅读节：把世界阅读日的4月23日定为“儿童读书日”；每年的4月23日到5月12日则为“儿童读书周”；每年10月27日到11月9日为“国民读书周”；日本政府还将2010年定为“国民读书年”，隆重推出一年的国民读书计划。

 

表3 10个峰的F1函数的Average mean和STD性能值Tab．3 Average mean and STD for F1 with 10 peaks

  

改变方式DECS Dopt-aiNet rPSO MDE Avg_mean STD T1 3．724456E －04 1．138583E －03 1．353000E －01 1．006100E+004．879480E+00 6．735230E+00 6．93 Avg_mean STD Avg_mean STD Avg_mean STD 0524E+00 5．216048E+00 T2 6．896028E －01 2．179447E+005．866700E+00 1．027720E+011．065260E+01 1．313120E+011．241597E+01 1．015788E+01 T3 3．863948E+00 6．985084E+004．254500E+00 8．182800E+001．663080E+01 1．775720E+011．131536E+01 9．456498E+00 T4 5．084533E+00 4．643826E+006．356300E+00 8．941400E+002．565840E+01 2．055290E+016．514131E+00 4．908173E+00 T5 9．389123E+00 1．531025E+014．435600E+00 5．554500E+004．596730E+00 5．777510E+002．787631E+01 2．617662E+01 T6 1．536066E+01 1．880491E+019．940700E+00 1．582140E+012．878160E+01 2．594580E+012．815238E+01 1．914847E+01 T7 1．464681E+00 4．614251E+004．211000E+00 8．687300E+001．428170E+01 1．372840E+015．089900E+00 4．393141E+00

 

表4 50个峰的F1函数的Average mean和STD性能值Tab．4 Average mean and STD for F1 with 50 peaks

  

改变方式DECS Dopt-aiNet rPSO MDE Avg_mean STD T1 1．660061E+00 1．905761E+00 3．644000E －01 9．275000E －017．682050E+00 7．217860E+00 1．51 Avg_mean STD Avg_mean STD Avg_mean STD 8407E+01 1．593068E+01 T2 3．580488E －01 1．005898E+004．748500E+00 6．758000E+001．209520E+01 1．064030E+011．241597E+01 1．015788E+01 T3 4．801628E+00 4．844080E+005．253100E+00 6．683000E+001．696330E+01 1．284230E+016．521749E+00 4．148847E+00 T4 6．490611E+00 7．322272E+002．656500E+00 5．977300E+002．096070E+01 2．015650E+016．463105E+00 5．919793E+00 T5 1．319286E+01 1．274812E+012．864100E+00 4．157900E+005．683460E+00 6．568420E+001．427066E+01 1．638993E+01 T6 2．669891E+01 3．177832E+016．833000E+00 1．187900E+013．974220E+01 2．759040E+013．557752E+01 3．175031E+01 T7 2．657925E+00 7．941768E+004．417200E+00 6．452800E+001．369470E+01 1．227190E+012．309890E+01 1．247930E+01

 

表5 F2函数的Average mean和STD性能值Tab．5 Average mean and STD for F2

  

改变方式DECS Dopt-aiNet rPSO MDE Avg_mean STD T1 5．022123E －01 1．112426E+00 9．840000E －02 2．910000E －024．589940E+01 3．322010E+01 2．81 Avg_mean STD Avg_mean STD Avg_mean STD 4148E+01 1．114949E+01 T2 1．789811E+01 2．626851E+018．120900E+00 1．438320E+011．710680E+02 2．153010E+023．929473E+01 2．979214E+01 T3 2．622813E+00 3．038497E+001．799790E+01 6．222590E+011．628950E+02 1．941210E+022．800671E+01 1．016845E+01 T4 2．435926E+02 1．248555E+021．065200E+00 2．826900E+005．364970E+01 6．119330E+012．418529E+02 1．218493E+02 T5 4．103075E+01 3．629524E+011．013840E+02 1．345180E+021．884630E+02 2．021740E+024．109725E+01 3．585491E+01 T6 4．104420E+01 3．648596E+016．519200E+00 1．381720E+017．563170E+01 8．313840E+014．070425E+01 3．577721E+01 T7 2．720657E+01 1．517024E+013．738500E+00 7．954200E+004．843600E+01 7．815160E+011．745018E+02 2．046530E+02

 

表6 F3函数的Average mean和STD性能值Tab．6 Average mean and STD for F3

  

改变方式DECS Dopt-aiNet rPSO MDE Avg_mean STD T1 6．274731E+02 3．683626E+018．108300E+02 6．610850E+016．448880E+02 2．615940E+025．43 Avg_mean STD Avg_mean STD Avg_mean STD 2566E+02 1．658885E+02 T2 5．817382E+02 1．181754E+021．078750E+03 6．412450E+019．365460E+02 1．351750E+025．759969E+02 2．127392E+02 T3 6．013543E+02 1．081089E+011．073430E+03 6．499500E+018．837220E+02 1．359240E+025．161730E+02 2．480069E+02 T4 6．570652E+02 3．828805E+011．031530E+03 2．747490E+027．811840E+02 3．021460E+026．599843E+02 5．219342E+01 T5 5．661777E+02 1．623613E+021．023900E+03 5．787130E+018．952690E+02 1．520220E+025．591761E+02 1．932096E+02 T6 4．104420E+01 3．648596E+011．186900E+03 2．922960E+029．425890E+02 3．254990E+026．235011E+02 1．139852E+02 T7 7．335582E+02 7．933669E+011．061830E+03 1．100980E+028．337610E+02 2．289910E+027．263300E+02 6．702290E+01

 

表7 F4函数的Average mean和STD性能值Tab．7 Average mean and STD for F4

  

改变方式DECS Dopt-aiNet rPSO MDE Avg_mean STD T1 1．376518E+01 9．994010E+00 1．422700E+00 4．545900E+00 5．927570E+01 1．038070E+02 2．81 Avg_mean STD Avg_mean STD Avg_mean STD 4148E+01 1．114949E+01 T2 3．619540E+01 3．671287E+01 1．224410E+02 2．016270E+02 2．822480E+02 2．615470E+02 3．929473E+01 2．979214E+01 T3 2．224852E+00 2．802811E+00 9．866880E+01 1．966950E+02 2．725670E+02 2．583130E+02 2．800671E+01 1．016845E+01 T4 2．881645E+02 1．470517E+02 4．263200E+00 9．725500E+00 1．175910E+02 1．711290E+02 2．871078E+02 1．443914E+02 T5 4．103075E+01 3．629524E+01 3．045660E+02 2．032430E+02 3．634430E+02 2．487290E+02 4．109725E+01 3．585491E+01 T6 4．089044E+01 3．621666E+01 1．263340E+01 5．583860E+01 8．654190E+01 1．517600E+02 4．107515E+01 3．577426E+01 T7 1．449823E+02 2．117304E+02 5．290100E+01 1．305930E+02 1．155390E+02 1．815390E+02 2．700468E+02 2．232451E+02

 

表8 F5函数的Average mean和STD性能值Tab．8 Average mean and STD for F5

  

改变方式DECS Dopt-aiNet rPSO MDE Avg_mean STD T1 1．713926E+01 1．624803E+014．089430E+01 2．212120E+028．216800E+01 1．333800E+021．91 Avg_mean STD Avg_mean STD Avg_mean STD 3104E+02 3．225623E+02 T2 2．173356E+01 3．417936E+013．445310E+01 1．198960E+025．412860E+01 8．843920E+011．566227E+02 2．551123E+02 T3 1．084363E+01 8．155296E+003．494200E+01 1．150250E+024．487450E+01 9．102550E+019．475097E+01 1．441250E+02 T4 1．683532E+01 1．922481E+011．206370E+02 2．935420E+028．421130E+01 1．431610E+023．780549E+02 4．512487E+02 T5 1．119155E+01 1．567978E+019．432230E+02 6．333180E+022．840390E+01 1．992450E+013．761274E+01 2．937091E+01 T6 9．568162E+00 1．446056E+014．803370E+02 6．108020E+029．878240E+01 1．633260E+023．743940E+01 3．191768E+01 T7 2．116921E+01 1．482409E+012．194660E+02 4．278170E+028．463810E+01 1．710430E+023．479443E+02 4．106790E+02

 

表9 F6函数的Average mean和STD性能值Tab．9 Average mean and STD for F6

  

改变方式DECS Dopt-aiNet rPSO MDE Avg_mean STD T1 8．172342E+00 8．277871E+002．044340E+01 7．932300E+015．995650E+01 8．791520E+012．04 Avg_mean STD Avg_mean STD Avg_mean STD 9801E+01 1．894677E+01 T2 3．527640E+01 3．080387E+013．911960E+02 3．954350E+022．356290E+02 3．167010E+024．792410E+01 2．719650E+01 T3 8．376330E －01 1．923465E+004．564410E+02 4．050380E+021．402160E+02 2．509680E+029．782493E+00 9．930505E+00 T4 2．187346E+01 2．280053E+018．396980E+01 2．201770E+028．005960E+01 1．272550E+028．903620E+01 1．760354E+02 T5 2．433638E+01 2．758829E+018．458620E+02 2．512080E+022．778160E+02 3．404960E+023．715577E+01 2．702133E+01 T6 3．374299E+01 3．032667E+014．822070E+02 4．344210E+021．152690E+02 1．970810E+023．310256E+01 2．834234E+01 T7 1．730255E+01 1．460075E+013．724740E+02 3．946680E+021．139600E+02 2．061110E+021．492469E+01 1．664048E+01

  

图2 不同动态优化问题的离线误差收敛趋势Fig．2 Convergence trend of off-line error for different DOP

  

图3 不同动态优化问题的评价值收敛趋势Fig．3 Convergence trend of evaluation value for different DOP

实验中主要衡量的评估量是平均离线误差期望和平均离线误差标准差，平均离线误差期望为离线误差函数Elast(t)在多次独立运行过程中多次动态变化的离线误差期望，能够有效衡量算法对动态环境问题的求解性能。平均误差标准差则反映算法在求解DOP时的鲁棒性。

从表3～表9中可以看出，在10个峰F1函数的T1～T4和T7变化问题，50个峰F1函数的T2、T3、T7变化问题，F2函数的T3、T5变化问题，F3函数的T1～T7变化问题，F4函数的T2、T3、T5变化问题，F5、F6函数的T1～T7变化问题上DECS的平均误差期望小于Dopt-aiNet算法。这些结果统计表明DECS在49个问题中有34个问题的求解结果优于DoptaiNet算法，49个问题中所有问题的求解结果优于rPSO算法和MDE算法。

从以上实验结果看出，DECS在F3函数、F5函数、F6函数上所有变化的求解都要优于Dopt-aiNet算法、rPSO算法、MDE算法。F3函数的基准函数为Rastrigin函数，其有多个局部最优值，是高度多模态的多峰函数，但多个峰的最小值位置是规律分布的。F5函数的基准函数为Ackley函数，其图像特征是外部区域有多个局部最小值邻域，中心区域有一个较大邻域的全局最优位置，且相对全局最优位置的邻域来说，外部区域几乎是平坦的。F6函数的基准函数具有5个函数的众多特征。这说明，DECS能够有效求解环境中特征明显的动态优化问题。

DECS在F2函数和F4函数上所有变化的求解都是只有3个变化超过Dopt-aiNet算法，所有变化超过rPSO算法和MDE算法。F2函数的基准函数为 Sphere函数，其在［－5．12，5．12］内是一个连续的、凸的、单峰的球面函数，尽管本实验中该函数搜索范围［－100，100］较大，但该函数相比本实验中其它函数却不具有较大的复杂性。F4函数的基准函数为Griewank函数，其有许多广泛分布的局部最小值区域，而且是规律分布的，这种复杂性在搜索范围放大的情况下会体现出来，本实验中该函数的搜索范围［－100，100］远大于Rastrigin函数的搜索范围［－5，5］，其峰的数量也远远大于Rastrigin函数峰的数量。实验数据表明，对于F2函数问题，DECS能够找到最优解的邻域，但Dopt-aiNet算法却具有更好的求解性能。而对于F4函数这种具有较大复杂性的环境，DECS性能优于rPSO算法和MDE算法，但和Dopt-aiNet算法性能相近，因此DECS在较大复杂性的环境中不具备明显优势。

步骤8 迭代次数iter←iter+1，如果达到总评价次数，则算法结束;否则，跳转步骤2。

DECS在10峰的F1函数的7个变化中有5种变化超过Dopt-aiNet算法，全部超过rPSO算法和MDE算法，但对于50个峰的F1函数，之所以在T1、T4变化上DECS效果要弱于Dopt-aiNet算法，那是因为相对于50个峰，实验用的搜索种群数量太少。但由于总评价次数是一定的，搜索种群数量太多，会加大对评价代价的消耗，因此搜索种群并不是越多越好。

江苏高考这几年作文题的材料聚焦于一到两个核心语词，行“话题作文”之实，写作者需要围绕核心语词，即“话题”，来构思行文。但也许有人会说，言论型材料基本都含有两三个核心词语或关键概念，构思行文理应围绕它们来写。但我们要注意的是，江苏高考作文题中的核心语词是独立存在型的，而非彼此联系。即使是2014年的“青春”，出题人的立足点也是“青春不朽”这一面，所以，实际上是以“青春不朽”这样的短语形式构成了一个独立的单位，而16年的“个性”或“创新”，两者之间也不一定非要发生关联。

对于49个函数中T1～T7的变化问题，在T3变化中的7个函数问题中，DECS比Dopt-aiNet算法、rPSO算法、MDE算法都要好;因为T3变化是随机变化，GDBG benchmark中T3变化的环境变化偏置量为一个服从正态分布的随机量与一个常数的乘积。这说明DECS能够有效求解环境随机变化的动态优化问题。

图2～3表示在部分问题上，DECS一次独立求解的离线误差收敛图、算法实际最优评价值和理论最优评价值的收敛图。其中:实线表示实际最优评价值，虚线表示理论最优评价值;横轴为评价次数，一次独立运行的评价次数为1．5E06，每评价1．5E05次环境改变。

在T2变化的7个函数问题中，DECS比Dopt-aiNet算法在6个问题中更好，比rPSO算法和MDE算法都要好;T2变化为环境大步的改变，这说明DECS能够有效求解环境变化明显的动态优化问题。

在T5、T7各自的7个函数变化问题中，DECS比DoptaiNet算法都是在5个问题上更好，比rPSO算法和MDE算法都要好。这说明DECS在求解周期的变化问题和维度变化问题上具有较好的性能。

在T1、T4、T6各自的7个函数变化问题中，DECS比DoptaiNet算法分别在3个、3个和4个问题上表现更好，比rPSO算法和MDE算法都要好，实验表明在这类问题上，DECS与Dopt-aiNet算法性能相差不多。这说明DECS在环境小步变化、混沌变化和带噪声的周期变化等问题上不具备明显优势。

在上述问题中，DECS之所以在基函数环境特征较明显的动态优化问题上具有良好性能，是因为其能快速感知环境变化并且在寻优过程中仍有较好的种群多样性。而对环境随机变化、大步变化、周期性变化和维度变化等问题上也具有良好性能，是因为DECS引入竞争机制后，平衡了种群的探索与寻优，使算法具有非常好的环境探索能力。

迟恒记得，办事处印发的综合治理简报上说征地拆迁、安置补偿方面的投入是792万元。网上报道过，国家财政专项拨款工程，地方上虚列开支套空额，有些工程，简直就是圈钱工程。所谓砂库废水、扬尘污染已使周边居民不堪忍受，背景离乡，看来也是夸大其词。

在DECS与Dopt-aiNet算法比较中，DECS在优于DoptaiNet算法的34个问题中有23个问题的平均误差标准差优于Dopt-aiNet算法，这说明在求解49个动态优化问题中，DECS不仅平均性能超过Dopt-aiNet算法，而且具有较好的鲁棒性。DECS与MDE算法比较过程中，DECS在所有49个问题中的平均误差期望均优于MDE算法，而在这49个问题中有41个问题的平均误差标准差要优于MDE算法;在DECS与rPSO算法的比较中，DECS在所有动态问题中的平均误差期望均优于rPSO算法，而在49个问题中有48个问题的平均误差标准差要优于rPSO算法，这说明在求解49个动态优化问题中，DECS不仅平均性能远远超过rPSO算法和MDE算法，而且具有更好的鲁棒性。

从图3看出一旦环境变化，环境中理论最优的评价值也变化，而DECS的实际最优评价值却能够及时监测并且快速追踪。由此看出DECS具备求解动态问题的能力。

5 结语

本文提出一种称为DECS的改进算法来求解DOP，并在统计学上表明是有效的。该算法在传统的DE算法框架内引入了多种群方法和竞争策略，使用多个种群并行搜索来增强搜索过程中的多样性。同时，多个搜索种群之间互相竞争，大大提高了算法的搜索能力，更高效地利用了有限的评估代价。排除规则和量子个体生成机制使搜索种群在搜索空间上具有更好的种群多样性;而独立的侦测种群，能够迅速感知环境变化。

1.3.2.5 心理护理 ①及时解除患者对疾病的焦虑、恐惧等不良心理活动;②加强沟通，使患者建立积极、乐观的健康心理。

实验结果统计表明，DECS与其他演化算法相比表现出更好的优化性能和更好的鲁棒性。未来的工作可以使搜索种群之间互相交流，淘汰掉的个体往往具有更优的多样性;也可以使用交叉概率、放缩因子与变异策略自适应地调整来提高算法的性能。

参考文献(References)

［1］ DAS S， MANDAL A， MUKHERJEE R．An adaptive differential evolution algorithm for global optimization in dynamic environments［J］．IEEE Transactions on Cybernetics， 2014， 44(6):966 －978．

［2］ BREST J， ZAMUDA A， BOSKOVIC B， et al．Dynamic optimization using self-adaptive differential evolution［C］//Proceedings of the 11th Conference on Congress on Evolutionary Computation．Piscataway，NJ:IEEE，2009:415－422．

［3］ JIN Y， BRANKE J．Evolutionary optimization in uncertain environments:a survey［J］．IEEE Transactions on Evolutionary Computation，2005，9(3):303－317．

［4］ TROJANOWSKI K， MICHALEWICZ Z．Evolutionary optimization in non-stationary environments［J］．Journal of Computer Science ＆Technology， 2000(3):1 －32．

［5］ 陈莉，丁立新．动态优化算法综述［J］．武汉大学学报(理学版)，2011，57(3):255 －264．(CHEN L，DING L X．A survey of dynamic optimization algorithms［J］．Journal of Wuhan University:Science Edition，2011，57(3):255－264．

［6］ YANG S， LI C．A clustering particle swarm optimizer for locating and tracking multiple optima in dynamic environments［J］．IEEE Transactions on Evolutionary Computation， 2010， 14(6):959 －974．

［7］ 姜立强，强洪夫，刘光斌．求解动态优化问题的改进差分进化算法［J］．小型微型计算机系统，2013，34(12):2837－2840．(JIANG L Q，QIANG H F，LIU G B．Modified differential evolution for dynamic optimization problems［J］．Journal of Chinese Computer Systems，2013， 34(12):2837 －2840．

［8］ 陈健，申元霞，纪滨．求解动态优化问题的多种群骨干粒子群算法［J］．计算机工程与应用， 2017，53(19):45 －50．(CHEN J，SHEN Y X，JI B．Multi-swarms bare bones particle swarms optimization of solving dynamic optimization problems［J］．Computer Engineering and Applications， 2017， 53(19):45 －50．)

［9］ ALI M Z， AWAD N H， SUGANTHAN P N， et al．An adaptive multipopulation differential evolution with dynamic population reduction［J］．IEEE Transactions on Cybernetics， 2016， 47(9):2768 －2779．

［10］ WEICKER K．An analysis of dynamic severity and population size［C］//Proceedings of the 6th International Conference on Parallel Problem Solving from Nature．London:Springer-Verlag， 2000:159－168．

［11］ ARAGON V S， ESQUIVEL S C．An evolutionary algorithm to track changes of optimum value locations in dynamic environments［J］．Journal of Computer Science＆ Technology，2004，4(3):127－133．

［12］ ROHLFSHAGEN P， YAO X．Dynamic combinatorial optimisation problems:an analysis of the subset sum problem［J］．Soft Computing，2011，15(9):1723－1734．

［13］ STORN R，PRICE K．Differential evolution:a simple and efficient heuristic for global optimization over continuous spaces［J］．Journal of Global Optimization，1997，11(4):341－359．

［14］ DAS S，SUGANTHAN P N．Differential evolution:a survey of the state-of-the-art［J］．IEEE Transactions on Evolutionary Computation，2011，15(1):4－31．

［15］ TASGETIREN M F，SUGANTHAN P N．A multi-populated differential evolution algorithm for solving constrained optimization problem［C］//Proceedings of the 2006 IEEE Congress on Evolutionary Computation．Piscataway， NJ:IEEE，2006:33 －40．

［16］ BLACKWELL T， BRANKE J．Multi-swarm optimization in dynamic environments［C］//Proceedings of the 2004 European Conference on Evolutionary Computation in Combinatorial Optimization．Berlin:Springer， 2004:489 －500．

［17］ MENDES R， MOHAIS A S．DynDE:a differential evolution for dynamic optimization problems［C］//Proceedings of the 2005 IEEE Congress on Evolutionary Computation．Piscataway， NJ:IEEE，2005:2808－2815．

［18］ LI C， YANG S， NGUYEN T T， et al．Benchmark generator for CEC 2009 competition on dynamic optimization［J］．Sensor Letters，2008，11(5):900－909．

［19］ DE F F O， VON Z F J．A dynamic artificial immune algorithm applied to challenging benchmarking problems［C］//Proceedings of the 11th Conference on Congress on Evolutionary Computation．Piscataway， NJ:IEEE，2009:423－430．