news 2026/4/3 3:56:51

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

作者头像

张小明

前端开发工程师

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/3/27 19:19:08

握着电烙铁调试C51板子的深夜,突然想起当年课程设计被密码锁支配的恐惧。今天就带大家用Proteus和Keil复刻这个经典项目,重点聊聊那些仿真时会遇到的坑

C51电子密码锁Proteus仿真 以AT89C51为主控芯片&#xff0c;LCD1602为显示屏设计密码为全数字的电子号码锁。 功能&#xff1a; (1)硬件部分&#xff1a;设计单片机外围电路、LCD显示电路、键盘电路等。 (2)软件部分&#xff1a;开锁、上锁、重新输入密码、设置新密码、保存新密…

作者头像 李华
网站建设 2026/3/28 18:56:11

x64dbg调试多线程程序注意事项

用x64dbg调试多线程程序&#xff1f;别让线程“乱跑”毁了你的分析你有没有遇到过这种情况&#xff1a;在x64dbg里设了个断点&#xff0c;结果一运行&#xff0c;程序频繁中断——不是你想调试的那个线程触发的&#xff0c;而是某个后台心跳线程、日志刷新线程或者GUI重绘线程不…

作者头像 李华
网站建设 2026/3/26 23:54:05

12bit SAR ADC电路:可直接仿真与频谱分析的电子工程应用

12bit sar adc电路&#xff0c;可直接仿真&#xff0c;逻辑模块也是实际电路&#xff0c;可利用cadence或者matlab进行频谱分析延申科普&#xff1a;ADC&#xff08;Analog-to-Digital Converter&#xff09;是一种电子设备&#xff0c;用于将连续的模拟信号转换为离散的数字信…

作者头像 李华
网站建设 2026/3/31 21:09:07

直流电机驱动电路在智能小车原理图中的应用实战案例

智能小车电机驱动设计实战&#xff1a;从原理到PCB落地的硬核经验你有没有遇到过这样的情况&#xff1f;代码写得没问题&#xff0c;传感器也正常读数&#xff0c;可一给电机通电——主控突然复位、轮子抖动不止、声音嗡嗡作响……最后排查半天&#xff0c;发现“罪魁祸首”竟是…

作者头像 李华
网站建设 2026/3/15 10:18:32

手把手玩转昆仑通泰触摸屏与V20变频器USS通讯

昆仑通泰触摸屏与v20变频器uss通讯&#xff0c;控制变频器 1&#xff0c;通过触摸屏与变频器uss通讯 2&#xff0c;通过触摸屏读取变频器电压&#xff0c;电流&#xff0c;频率 3&#xff0c;通过触摸屏设定变频器正反转&#xff0c;设定频率&#xff0c;加速&#xff0c;减速时…

作者头像 李华