news 2026/7/24 18:09:18

基于遗传算法的最优化求解问题:简单多元函数极值问题的优化与通用性代码

作者头像

张小明

前端开发工程师

1.2k 24
文章封面图
基于遗传算法的最优化求解问题:简单多元函数极值问题的优化与通用性代码

基于遗传算法的最优化求解问题 其通过遗传算法对简单的多元函数求极值问题进行优化求解,得到了最优解和迭代收敛曲线 代码通用性很好 很适合小白入门

今天咱们来聊聊怎么教计算机自己找答案——这事儿听着玄乎,其实用遗传算法就能轻松搞定。举个栗子,假设我们要找函数f(x,y)=x²+y²的最小值,人类肉眼一看就知道在(0,0)处,但计算机可不会代数运算,这时候就该让基因进化大法上场了。

先准备个种群,就像养鱼塘放鱼苗:

import numpy as np def create_population(size, dim, lbound, rbound): return np.random.uniform(lbound, rbound, (size, dim))

这函数生成了20条染色体(解),每条有x,y两个基因,取值范围定在-10到10。好比在鱼塘里随机撒下不同品种的鱼苗,等着看哪些能适应环境。

接下来是适者生存环节:

def fitness(individual): return -sum(individual**2) # 求最小值转换为求负数的最大值 def selection(pop, fitnesses, elite_size=4): elite_indices = np.argsort(fitnesses)[-elone_size:] return pop[elite_indices]

这里有个骚操作:把求最小值转换成求负数的最大值。选择阶段就像选美比赛,表现最好的前4名直接晋级,保留下优质基因。

重点来了——交配和变异:

def crossover(parent1, parent2): cross_point = np.random.randint(1, len(parent1)) return np.hstack((parent1[:cross_point], parent2[cross_point:])) def mutate(individual, mutation_rate=0.1): mask = np.random.rand(len(individual)) < mutation_rate individual[mask] += np.random.normal(0, 0.5, sum(mask)) return individual

交叉操作像父母各剪一段基因拼接成新个体,变异则是随机给某些基因加点噪声。注意变异幅度用正态分布,这样既不会突变太猛也不会完全不变,比固定步长聪明多了。

整套进化流程封装成:

def genetic_algorithm(func, dim=2, iter_num=50): pop = create_population(20, dim, -10, 10) best_values = [] for _ in range(iter_num): fitnesses = np.array([func(ind) for ind in pop]) elites = selection(pop, fitnesses) selected = elites.copy() while len(selected) < 20: parents = elites[np.random.choice(len(elites), 2, replace=False)] child = crossover(*parents) child = mutate(child) selected = np.vstack((selected, child)) pop = selected best_values.append(-np.max(fitnesses)) # 转换回原函数值 return best_values

迭代过程中始终保持种群数量稳定,每一代都保留精英并补充新血。记录最优解的变化趋势,方便后续画收敛曲线。

跑完50代后画个进化轨迹:

import matplotlib.pyplot as plt history = genetic_algorithm(fitness) plt.plot(history, 'g', linewidth=2) plt.title('进化过程曲线') plt.xlabel('迭代次数') plt.ylabel('函数值') plt.show()

你会发现曲线像滑梯一样往下溜,最后稳稳停在0附近。实际跑几次会发现,基本在10代以内就能找到精度1e-4级别的解,比随机搜索高效不止一个量级。

这套代码的妙处在于改个目标函数就能解决其他优化问题。比如把fitness函数改成sum(np.sin(individual)individual*2),立马变成求解震荡函数的最小值。参数调整也简单,种群规模、变异率、选择数量这些都可以根据问题复杂度灵活配置。

刚入门的同学可能会问:为什么不直接用梯度下降?问得好!遗传算法的优势在于:1)不依赖梯度信息,2)能跳出局部最优,3)特别适合离散/非线性问题。当然,计算量确实大些,但对于小白理解优化算法的本质,这可比推导数学公式直观多啦。

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

打开COMSOL看到电磁波模块就手痒?今天拿介质圆柱散射练练手。先搞个半径5μm的氧化铝圆柱(ε_r=9.8),扔到532nm激光里会发生啥?咱们边操作边唠嗑

COMSOL介质圆柱散射效率分析。 也可分析散射截面&#xff0c;消光截面与吸收截面。建模时直接在几何里画个圆&#xff0c;边界条件记得套两层&#xff1a;里面是散射边界&#xff08;别让波反射回来捣乱&#xff09;&#xff0c;外面包个完美匹配层。材料库调出氧化铝参数时注意…

作者头像 李华
网站建设 2026/7/20 12:16:23

当C#遇上工业PLC:手撕多品牌通讯源码实录

C#与三菱&#xff0c;西门子&#xff0c;台达&#xff0c;基恩士&#xff0c;等各品牌plc通讯源码。搞过工控的老铁都知道&#xff0c;PLC通讯就像和不同方言的人聊天——三菱说MC协议&#xff0c;西门子玩S7&#xff0c;台达可能掏出Modbus&#xff0c;基恩士说不定甩个自定义…

作者头像 李华
网站建设 2026/7/20 21:28:50

零基础秒变分析高手!覆盖7大行业的超全模版中心来了

每天打开后台&#xff0c;五花八门的用户数据堆得满满当当&#xff0c;手上的需求还没收尾&#xff0c;新的需求已经排着队找上门&#xff1b;熬夜赶出来的报表写满了&#xff0c;却迟迟找不到问题所在……“数据报表越做越厚&#xff0c;业务洞察却越来越慢”——这或许是许多…

作者头像 李华
网站建设 2026/7/20 17:09:06

适者生存的淘汰逻辑

COA-SVM&#xff1a;土狼优化算法优化SVM&#xff1a;COA-SVM。2018提出一种用于全局优化的自然启发式算法&#xff0c;可以用来写在机器学习炼丹房里&#xff0c;总有些传统算法需要点"野性"才能爆发真正实力。今天咱们聊聊怎么让土狼群撕开SVM参数优化的困局——这…

作者头像 李华
网站建设 2026/7/20 18:20:40

全面讲解如何测试与调试数据库触发器

如何真正掌控数据库触发器&#xff1a;从测试到调试的实战全解在现代系统开发中&#xff0c;有一个“低调却致命”的组件&#xff0c;它不显山露水&#xff0c;却能在关键时刻决定数据是否一致、事务能否提交、甚至整个服务会不会雪崩——那就是数据库触发器。你可能已经用它来…

作者头像 李华