news 2026/4/16 21:30:55

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

作者头像

张小明

前端开发工程师

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/4/16 10:39:22

避免常见错误:Python转可执行文件新手教程

打包不翻车&#xff1a;手把手教你把 Python 脚本变成“即点即用”的 .exe 你有没有过这样的经历&#xff1f; 辛辛苦苦写了个数据处理工具&#xff0c;同事双击运行却弹出“找不到 Python”&#xff1b;或者打包完的 .exe 文件大得离谱——一个简单脚本竟有 150MB&#xf…

作者头像 李华
网站建设 2026/4/16 3:33:00

一文说清USB转串口驱动在产线自动化中的核心作用

USB转串口驱动&#xff1a;产线自动化中那根“看不见的神经”在一条高速运转的电子产品生产线里&#xff0c;你或许不会注意到这样一个细节&#xff1a;当一块电路板被送入测试工位&#xff0c;上位机几秒钟内就完成了固件烧录、参数校准和功能自检——整个过程行云流水。但你知…

作者头像 李华
网站建设 2026/4/16 12:28:41

欧姆龙NJ501-1500枕包机:六轴定位控制,初学者的理想之选

枕包机 欧姆龙NJ501-1500,六轴定位控制,凸轮控制&#xff0c; 程序经典&#xff0c;结构清晰&#xff0c;梯形图编程&#xff0c;特别适合初学者入门&#xff1b;最近在调试枕包机时发现欧姆龙NJ501-1500真是块宝藏PLC&#xff0c;尤其是它自带的六轴定位功能&#xff0c;能把送…

作者头像 李华
网站建设 2026/4/16 12:28:10

26. AI-框架工具-LangChain LangGraph

文章目录前言一、LangChain1.什么是LangChain&#xff1f;2. 核心功能与主要特点1. 核心功能2. 主要特点3. 架构组件与工作原理1. 模型I/O&#xff08;Model I/O&#xff09;2. 提示工程&#xff08;Prompt Engineering&#xff09;3. 链&#xff08;Chains&#xff09;4. 记忆…

作者头像 李华
网站建设 2026/4/16 15:37:34

一、需求拆解与参数计算核心

制袋机横切机程序威纶通触摸屏和三菱3UPLC加模拟量4DA模块&#xff0c;控制两台变频器&#xff08;主机和给料辊筒)和一台送料伺服。 通过触摸屏设定一分钟张数&#xff0c;程序自动算出伺服送料速度&#xff0c;和后面给料速度。 主要是参数计算思路。 程序带详细注释&#xf…

作者头像 李华