《一种基于ASU策略的多目标粒子群优化算法.pdf》由会员分享，可在线阅读，更多相关《一种基于ASU策略的多目标粒子群优化算法.pdf（9页完整版）》请在专利查询网上搜索。
1、(10)申请公布号 CN 103854131 A (43)申请公布日 2014.06.11 CN 103854131 A (21)申请号 201410016605.X (22)申请日 2014.01.14 G06Q 10/06(2012.01) G06Q 50/06(2012.01) G06N 3/12(2006.01) (71)申请人 国家电网公司 地址 100033 北京市西城区西长安街 86 号 申请人 国网浙江省电力公司 浙江省电力设计院 (72)发明人 孙可 郑朝明 丁晓宇 陈宏伟 宁康红 徐凯 江全元 邹杨 郑晓 赵萌 (74)专利代理机构 杭州九洲专利事务所有限公 司 33101。
2、 代理人 翁霁明 (54) 发明名称 一种基于 ASU 策略的多目标粒子群优化算法 (57) 摘要 一种基于 ASU 策略的多目标粒子群优化算 法， 该优化算法利用基于 ASU 策略进行精英集缩 减， 能提高最优解分布的均匀性和多样性， 且能搜 索到相对更优的解， 从而提高最终结果的可靠性 ； 它对于一个多目标优化问题， 在精英集缩减过程 中， 将粒子群分成若干个网格， 然后在每个网格 中， 提出利用逐步更新的欧氏拥挤距离去筛选粒 子， 避免一次性剔除一群集中的粒子而破坏解分 布的多样性和均匀性， 且能搜索到相对更优的解。 (51)Int.Cl. 权利要求书 1 页 说明书 4 页 附图 3 。
3、页 (19)中华人民共和国国家知识产权局 (12)发明专利申请 权利要求书1页 说明书4页 附图3页 (10)申请公布号 CN 103854131 A CN 103854131 A 1/1 页 2 1. 一种基于 ASU 策略的多目标粒子群优化算法， 其特征在于该优化算法包括如下步 骤 ： （1） 种群初始化， 设置精英集为空 ； 迭代次数 t=0， 在控制变量变化范围内， 随机初始化 粒子群， 每个粒子的个体极值和全局极值均为初始位置， 设定精英集规模 ； （2） 计算每个粒子的目标函数适应值 ； （3） 判断是否满足结束准则， 达到最大允许迭代次数或最优解对应的目标函数值在给 定的迭代步数。
4、内改变量小于给定值时， 停止优化并输出结果， 形成 Pareto 最优解 ； 否则迭 代次数 t=t+1， 转到步骤 （4） ； （4） 根据 Pareto 支配， 比较各个粒子之间的优劣， 构造粒子群的非支配解集 ； （5） 通过比较非支配解集和精英集中粒子的 Pareto 支配关系， 更新精英集 ； （6） 判断精英集规模是否超过给定规模， 若超过， 则利用 ASU 策略缩减至给定的规模， 并转到步骤 （7） ； 反之， 直接转到步骤 （7） ； （7） 更新粒子的个体极值和全局极值 ； （8） 更新每个粒子的位置和速度， 构造新的粒子， 转到步骤 （2） 。 2. 根据权利要求 1 所述。
5、基于 ASU 策略的多目标粒子群优化算法， 其特征在于所述步骤 5 中， 在精英集的规模超出给定规模时， 利用 ASU 策略进行缩减， 即将粒子群分成若干个网 格， 然后在每个网格中， 利用逐步更新的欧氏拥挤距离去筛选粒子， 即每计算一次欧氏拥挤 距离之后， 剔除欧氏拥挤距离最小的精英解， 然后再计算一次欧氏拥挤距离， 再剔除此时欧 氏拥挤距离最小的精英解， 按此方法， 直至缩减到规定的精英集规模 ； 其中欧氏拥挤距离 d 表示为 x1i选中粒子第 i 维的位置， x2i和 x3i与 x1i最近两个粒子第 i 维的位置， n粒子的维数。 权 利 要 求 书 CN 103854131 A 2 1。
6、/4 页 3 一种基于 ASU 策略的多目标粒子群优化算法 技术领域 0001 本发明涉及的是一种基于 ASU （Asynchronous-Stepwise Updated, 异步 - 逐步更 新） 策略的多目标粒子群优化算法， 属于电力系统规划技术领域。 背景技术 0002 随着我国经济的高速发展和电力系统的不断扩大， 人们对供电的安全性、 可靠性、 质量等方面的要求也日益增长， 势必需要大量的输变电工程投入到电网建设中。为应对电 网建设需求， 每年都会大量的输变电工程上报。在这些工程中， 有一部分是电网建设需要 的， 他们能解决电网的某个缺陷， 但另外有一部分也是没有必要的， 或者上报的方。
7、案并不是 最优的， 需要替换为更好的实施方案。 如何在大量的工程中选择需要的， 并且以一种最好的 方案进行实施， 对于电网建设决策者而言， 至关重要。 0003 输变电工程立项决策方法的基本任务是对上报申请的待选项目， 根据对电网的全 面评价确定其是否有必要立项， 如果立项， 则合理地决策用哪一种方案进行针对性地改造 和建设， 从而使规划方案能满足电网安全可靠运行要求， 适应电力发展要求， 同时使能源资 源得到合理的优化配置。 0004 输变电工程立项决策问题是一个非线性、 多阶段、 多目标的复杂优化问题。目前， 对于输变电工程立项决策问题往往建立成一个多目标优化模型， 利用传统的基于拥挤距离。
8、 策略的方法求解时， 一次性剔除小于某个距离的所有粒子， 导致部分区域粒子缺失， 且不适 用于高维模型 ； 而利用基于网格策略的方法求解时， 会存在两个非常相近的非支配解被划 分到两个相邻格子而一直留在解集中这样的不足。 可见传统的方法求解时不能很好地保持 解的多样性和分布的均匀性， 从而导致最后得到的结果并非最优的。 因此， 急需一种能较好 保持解多样性和分布的均匀性的多目标优化算法求解输变电工程立项决策问题， 从而为工 程实际提供更加可靠的参考结果。 发明内容 0005 本发明的目的是为了在多目标粒子群优化方法中， 提高计算结果的最优性和可靠 性， 为此， 本发明提供了一种基于 ASU 策。
9、略的多目标粒子群优化算法， 该方法利用 ASU 策略 进行精英集缩减， 提高最优解分布的均匀性和多样性， 且能搜索到相对更优的解， 从而提高 最终结果的可靠性。 0006 本发明的目的是通过如下技术方案来实现的， 所述的优化算法包括如下步骤 ： 0007 （1） 种群初始化。设置精英集为空 ； 迭代次数 t=0， 在控制变量变化范围内， 随机初 始化粒子群， 每个粒子的个体极值和全局极值均为初始位置， 设定精英集规模 ； 0008 （2） 计算每个粒子的目标函数的适应值 ； 0009 （3） 判断是否满足结束准则， 达到最大允许迭代次数或最优解对应的目标函数值 在给定的迭代步数内改变量小于给定。
10、值时， 停止优化并输出结果， 形成 Pareto 最优解 ； 否 则迭代次数 t=t+1， 转到步骤 （4） ； 说 明 书 CN 103854131 A 3 2/4 页 4 0010 （4） 根据 Pareto 支配， 比较各个粒子之间的优劣， 构造粒子群的非支配解集 ； 0011 （5） 通过比较非支配解集和精英集中粒子的 Pareto 支配关系， 更新精英集 ； 0012 （6） 判断精英集规模是否超过给定规模， 若超过， 则利用 ASU 策略缩减至给定的规 模， 并转到步骤 （7） ； 反之， 直接转到步骤 （7） ； 0013 （7） 更新粒子的个体极值和全局极值 ； 0014 （8。
11、） 更新每个粒子的位置和速度， 构造新的粒子， 转到步骤 （2） 。 0015 所述步骤 （5） 中， 在精英集的规模超出给定规模时， 利用 ASU 策略进行缩减， 即将 粒子群分成若干个网格， 然后在每个网格中， 利用逐步更新的欧氏拥挤距离去筛选粒子， 即 每计算一次欧氏拥挤距离之后， 剔除欧氏拥挤距离最小的精英解， 然后再计算一次欧氏拥 挤距离， 再剔除此时欧氏拥挤距离最小的精英解， 按此方法， 直至缩减到规定的精英集规 模 ； 其中欧氏拥挤距离 d 表示为 0016 0017 x1i选中粒子第 i 维的位置， 0018 x2i和 x3i与 x1i最近两个粒子第 i 维的位置， 0019 。
12、n粒子的维数。 0020 ASU 策略能避免一次性剔除一群集中的粒子而破坏解分布的多样性和均匀性， 且 能搜索到相对更优的解。 0021 本发明的有益效果是 ： 本发明采用了基于 ASU 策略的多目标粒子群算法， 该方法 利用基于 ASU 策略进行精英集缩减， 能提高最优解分布的均匀性和多样性， 且能搜索到相 对更优的解， 从而提高最终结果的可靠性。 与现有的多目标粒子群算法相比， 本发明提出的 方法主要有以下改进 ： 即 ASU 策略在精英集缩减过程中， 能更好地保持解分布的均匀性和 多样性， 且能搜索到相对更优的解。 附图说明 0022 图 1 是基于 ASU 策略的多目标粒子群优化算法的。
13、流程示意图 ； 0023 图 2 是三种方法对 ZDT1 函数的测试结果图 ； 0024 图 3 是三种方法对 ZDT2 函数的测试结果图 ； 0025 图 4 是三种方法对 ZDT3 函数的测试结果图 ； 具体实施方式 0026 下面将结合附图及实施例对本发明做详细的介绍 ： 结合图 1 所示， 本发明所述的 基于 ASU 策略的多目标粒子群优化算法包括如下步骤 ： 0027 （1） 种群初始化。设置精英集为空 ； 迭代次数 t=0， 在控制变量变化范围内， 随机初 始化粒子群， 每个粒子的个体极值和全局极值均为初始位置， 设定精英集规模 ； 0028 （2） 计算每个粒子的目标函数的适应值。
14、 ； 0029 （3） 判断是否满足结束准则， 达到最大允许迭代次数或最优解对应的目标函数值 在给定的迭代步数内改变量小于给定值时， 停止优化并输出结果， 形成 Pareto 最优解 ； 否 说 明 书 CN 103854131 A 4 3/4 页 5 则迭代次数 t=t+1， 转到步骤 （4） ； 0030 （4） 根据 Pareto 支配， 比较各个粒子之间的优劣， 构造粒子群的非支配解集 ； 0031 （5） 通过比较非支配解集和精英集中粒子的 Pareto 支配关系， 更新精英集 ； 0032 （6） 判断精英集规模是否超过给定规模， 若超过， 则利用 ASU 策略缩减至给定的规 模，。
15、 并转到步骤 （7） ； 反之， 直接转到步骤 （7） ； 0033 （7） 更新粒子的个体极值和全局极值 ； 0034 （8） 更新每个粒子的位置和速度， 构造新的粒子， 转到步骤 （2） 。 0035 ASU 策略是指在精英集缩减时， 将粒子群分成若干个网格， 然后在每个网格中， 利 用逐步更新的欧氏拥挤距离去筛选粒子， 即每计算一次欧氏拥挤距离之后， 剔除欧氏拥挤 距离最小的精英解， 然后再计算一次欧氏拥挤距离， 再剔除此时欧氏拥挤距离最小的精英 解， 按此方法， 直至缩减到规定的精英集规模 ； 0036 实施例 ： 为了体现本发明算法的性能优势， 这里选择 ZDT1 ZDT3 函数来测。
16、试， 这 三个函数分别具有凸特性、 凹特性和离散特性。性能指标选择用世代距离 （generational distance,GD） 体现 Pareto 最优解的最优性， 用解集均匀度 （Spacing,SP） 来体现 Pareto 最 优解分布的均匀性和多样性。 0037 世代距离指算法所得的解与 Pareto 最优前端之间的趋近程度。 0038 0039 di目标空间中解和对应 Pareto 最优解的欧氏距离， n所得解的数量 ； 0040 GD 的值越小， 则算法获得的解越趋近于最优前段 ； 0041 解集均匀度表示算法所得的解在目标空间上分布的均匀程度 ； 0042 0043 di目标空。
17、间中解和对应 Pareto 最优解的欧氏距离， n所得解的数量， d所有 di的平均值 ； SP 的值越小， 则算法获得的解越均匀。 0044 在本发明算法， 以及基于拥挤距离策略的算法和基于网格策略的算法中， 设置种 群规模都为 100， 精英集规模为 100。三种算法对 ZDT1 ZDT3 函数都运行 30 次， 统计 GD 和 SP 的均值和方差， 见表 1 所示 ： 0045 表 1 三种算法的测试结果 说 明 书 CN 103854131 A 5 4/4 页 6 0046 0047 从表中可以看出， 基于 ASU 策略的多目标粒子群优化算法对三个测试函数求得的 GD 的均值和方差、 SP 的均值和方差都小于其他两种算法。由此可见， 本发明算法能够更好 地保持 Pareto 最优解分布的均匀性和多样性， 且能搜索到相对更优的解。 说 明 书 CN 103854131 A 6 1/3 页 7 图 1 说 明 书 附 图 CN 103854131 A 7 2/3 页 8 图 2 图 3 说 明 书 附 图 CN 103854131 A 8 3/3 页 9 图 4 说 明 书 附 图 CN 103854131 A 9 。