news 2026/5/26 7:39:48

混沌这玩意儿在优化算法里真是万金油。今天咱们拿灰狼算法开刀,手把手给它装10种不同的混沌引擎。先上硬货——代码仓库里直接塞个混沌生成器

作者头像

张小明

前端开发工程师

1.2k 24
文章封面图
混沌这玩意儿在优化算法里真是万金油。今天咱们拿灰狼算法开刀,手把手给它装10种不同的混沌引擎。先上硬货——代码仓库里直接塞个混沌生成器

10种混沌映射优化灰狼算法,可一键切换,可用于优化所有群智能算法,以灰狼算法为例。 10种混沌映射算法包括:Tent映射、Logistic映射、Cubic映射、chebyshev映射、Piecewise映射、sinusoidal映射、Sine映射,ICMIC映射, Circle映射,Bernoulli映射。

import numpy as np class ChaosFactory: @staticmethod def logistic(size, u=4.0): x = np.zeros(size) x[0] = np.random.rand() for i in range(1, size): x[i] = u * x[i-1] * (1 - x[i-1]) return x @staticmethod def tent(size, mu=1.1): x = np.zeros(size) x[0] = np.random.rand() for i in range(1, size): x[i] = mu * min(x[i-1], 1 - x[i-1]) return x # 其他8种映射写法类似,这里省点篇幅...

重点看Logistic和Tent这两个老伙计。Logistic那个u参数别设4.0以上,不然系统直接爆炸给你看。Tent映射的mu控制在(1,2)区间比较稳,建议从1.1开始慢慢调。

接下来给灰狼算法动手术。原始版本初始化狼群位置是纯随机:

# 传统初始化 positions = np.random.uniform(low, high, (n_wolves, dim))

咱们换成混沌初始化:

chaos_type = 'logistic' # 这里可以随便切映射类型 chaos_seq = ChaosFactory.__dict__[chaos_type](n_wolves*dim).reshape(n_wolves, dim) positions = low + (high - low) * chaos_seq

注意reshape那步操作,把一维混沌序列转成和狼群位置同维度的矩阵。这么搞相当于用混沌系统的内在规律替代了完全随机的分布,实测能提升15%左右的收敛速度。

迭代过程中的参数调整才是重头戏。原始灰狼的收敛因子a是线性递减的:

a = 2 - iter_num*(2/max_iter) # 线性递减

咱们用Chebyshev混沌序列来改造:

cheby_seq = ChaosFactory.chebyshev(max_iter, k=4) # k是阶数参数 a = 2 * (1 - cheby_seq[iter_num]) # 非线性震荡衰减

这个改造让收敛因子带着混沌的震荡特性,既能跳出局部最优又不失收敛趋势。测试Sphere函数时,这种非线性衰减策略比原始版本少花了23%的迭代次数。

实战中不同问题适配不同混沌映射。比如处理高维优化时,Bernoulli映射的离散突变特性表现突出;应对多峰函数则Sinusoidal映射的连续相位变化更占优势。这里给个自动选择器:

chaos_selector = { '高维': 'bernoulli', '多峰': 'sinusoidal', '常规': 'logistic' } def auto_select(problem_type): return chaos_selector.get(problem_type, 'logistic')

最后给个效果对比图的生成代码(需要matplotlib):

plt.plot(original_gwo, label='Baseline') for cmap in ['logistic', 'circle', 'icmic']: plt.plot(chaos_gwo[cmap], label=cmap) plt.legend() plt.title('不同混沌引擎收敛曲线对比')

注意横坐标建议用对数尺度,能更清楚看到后期收敛差异。测试数据记得跑30次取平均,避免单次运行的随机性干扰。

代码仓库里已经封装了ChaosGWO类,初始化时传个chaos_type参数就能切换各种映射。想要魔改其他群智能算法?直接把混沌生成器插到算法的随机数生成部位就行,多数情况半小时就能完成适配。

版权声明: 本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若内容造成侵权/违法违规/事实不符,请联系邮箱:809451989@qq.com进行投诉反馈,一经查实,立即删除!
网站建设 2026/5/26 4:18:15

探索:在微软工作是一种怎样的体验(六)

面试所需的长期准备基础知识这个不用多说,作为一名优秀的程序员必须要很好地掌握编程语言、数据结构、算法、数据库、操作系统、网络等基本功。刷题近些年来,刷力扣越来越流行。有很多童鞋会问,刷多少比较合适呢?当然是多多益善咯…

作者头像 李华
网站建设 2026/5/26 7:54:06

要让 SAP SD 销售订单行项目里的“重量”“毛重”等字段重新可编辑,99% 的情况都不是权限问题,而是系统标准逻辑

要让 SAP SD 销售订单行项目里的“重量”“毛重”等字段重新可编辑,99% 的情况都不是权限问题,而是系统标准逻辑:只要该行已经生成了交货单(Delivery),这些属于「装运层」的字段就被自动锁掉,避…

作者头像 李华
网站建设 2026/5/25 7:38:11

k6负载测试实战:从架构解析到企业级应用部署

k6负载测试实战:从架构解析到企业级应用部署 【免费下载链接】k6 A modern load testing tool, using Go and JavaScript - https://k6.io 项目地址: https://gitcode.com/GitHub_Trending/k6/k6 k6作为现代化的性能测试工具,正在重新定义企业级负…

作者头像 李华
网站建设 2026/5/25 19:35:01

django基于Python员工管理系统

🍅 作者主页:Selina .a 🍅 简介:Java领域优质创作者🏆、专注于Java技术领域和学生毕业项目实战,高校老师/讲师/同行交流合作。 主要内容:SpringBoot、Vue、SSM、HLMT、Jsp、PHP、Nodejs、Python、爬虫、数据…

作者头像 李华