改进的鲸鱼算法(EWOA)优化BP神经网络原理及matlab代码

目录

0 引言

1 数学模型

2 模型对比

3 matlab代码

3.1 伪代码

3.2 EWOA-BP

4 视频讲解

0 引言

WOA被广泛应用于解决NPhard问题。然而,WOA种群多样性低和搜索策略差。因此,Mohammad H等人在2022年提出了一种增强的鲸鱼优化算法EWOA,利用池化机制和三种有效的搜索策略优化WOA。

1 数学模型

EWOA是在WOA基础上引入池化机制、迁移、优先选择和丰富的包围猎物等搜索策略,以提高WOA的多样性和全局性

(1)池化机制:该机制是一种交叉策略,利用最好座头鲸和最差鲸鱼进行二进制随机向量和反向向量修正组合形成新的种群,提高模型多样性,其数学模型见下式:

式中B、B(-)为二进制随机向量和反向量,Xbmd最好座头鲸,Xworst为最差的鲸鱼。

(2)迁移搜索策略:使用等式随机分离出座头鲸的一部分覆盖无人区域,改善勘探。增加整个种群的多样性,从而减少局部最佳诱捕。其数学模型见下式:

式(12)为座头鲸迁移策略的位置;式(13)为使用等式计算的搜索空间范围内的一个随机位置,δmin/max为问题上下界;式(14)为最佳座头鲸附近的一个随机位置,δbest-min/max分为该问题对应上下界限。

(3)优先选择搜索策略:取代WOA算法的搜索狩猎行为

式中Prnd1/rnd2是在迭代过程矩阵池中随机选择种群,Ai是经过柯西分布取样系数向量。

(4)丰富的包围猎物搜索策略:取代WOA中搜索猎物阶段

2 模型对比

对比WOA和EWOA优化BP神经网络的模型性能得出对比图,由于BP神经网络受权值阈值随机性影响,模型对比选择固定随机组函数进行,具体代码如下:

rand('seed',333)%333为随机固定组,可修改

模型数据来源UCI回归数据集:

适应度对比:

精度误差对比:

EWOA-BP:

WOA-BP:

3 matlab代码

3.1 伪代码

3.2 EWOA-BP

1)单输出回归预测模型:​​​​​​​单输出回归模型:改进的鲸鱼算法(EWOA)优化BP神经网络模型(EWOA-BP)

2)多输出回归预测模型:多输出回归模型:改进的鲸鱼算法(EWOA)优化BP神经网络模型(EWOA-BP)

3)分类模型:分类模型:改进的鲸鱼算法(EWOA)优化BP神经网络模型(EWOA-BP)

4)随机序列:​​​​​​​时间序列:改进的鲸鱼算法(EWOA)优化BP神经网络模型(EWOA-BP)

4 视频讲解

B站搜索:‘ 不想学习的陈成 ’

  • 22
    点赞
  • 22
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
WOA-BP鲸鱼算法优化BP神经网络一种常见的神经网络优化方法。下面是一些基本概念和实现步骤: 1. WOA-BP鲸鱼算法简介 WOA-BP鲸鱼算法一种基于鲸鱼群智能优化算法BP神经网络优化方法。其基本思想是通过模拟鲸鱼的迁徙和捕食行为来寻找最优解。 2. BP神经网络简介 BP神经网络一种常见的人工神经网络模型,其基本思想是通过反向传播算法来训练神经网络,从而实现对输入输出之间的映射关系进行学习和建模。 3. WOA-BP鲸鱼算法优化BP神经网络步骤 (1) 初始化BP神经网络参数和WOA算法参数; (2) 根据当前种群位置计算适应度函数值; (3) 利用WOA算法更新种群位置; (4) 根据更新后的位置计算新的适应度函数值,并根据新的适应度函数值对种群进行排序; (5) 判断是否满足停止条件,如果满足,则输出最优解,否则返回步骤2。 4. Matlab代码实现 以下是Matlab代码实现的基本框架: ``` % Step 1: 初始化BP神经网络参数和WOA算法参数 % Step 2: 根据当前种群位置计算适应度函数值 % Step 3: 利用WOA算法更新种群位置 % Step 4: 根据更新后的位置计算新的适应度函数值,并根据新的适应度函数值对种群进行排序 % Step 5: 判断是否满足停止条件,如果满足,则输出最优解,否则返回步骤2 % 以下是一个简单的示例代码: % Step 1: 初始化BP神经网络参数和WOA算法参数 pop_size = 10; % 种群大小 max_iter = 100; % 最大迭代次数 dim = 10; % 每个个体的维度 c1 = 2; % 常数c1 c2 = 2; % 常数c2 c3 = 2; % 常数c3 a = 2; % 常数a x_max = 100; % 变量x的上限 x_min = -100; % 变量x的下限 w_max = 1; % 权重w的上限 w_min = -1; % 权重w的下限 pop_position = rand(pop_size,dim); % 随机初始化种群位置 pop_fitness = zeros(1,pop_size); % 初始化种群适应度函数值 % Step 2: 根据当前种群位置计算适应度函数值 for i=1:pop_size pop_fitness(i) = fitness_func(pop_position(i,:)); % 计算适应度函数值 end % Step 3: 利用WOA算法更新种群位置 for t=1:max_iter % 迭代次数循环 for i=1:pop_size % 种群个体循环 r1 = rand(); r2 = rand(); A = 2*a*r1-a; C = 2*r2; b = 1; l = (a-1)*rand()+1; p = rand(); % 随机生成参数p if p<0.5 % 更新个体位置 for j=1:dim if rand()<0.5 D = abs(C*pop_position(i,j)-pop_position(i,j)); pop_position(i,j) = D*exp(b*l)*cos(2*pi*l)+pop_position(i,j); else D = abs(C*pop_position(i,j)-pop_position(i,j)); pop_position(i,j) = D*exp(b*l)*sin(2*pi*l)+pop_position(i,j); end if pop_position(i,j)>x_max % 边界处理 pop_position(i,j) = x_max; elseif pop_position(i,j)<x_min pop_position(i,j) = x_min; end end else % 更新种群位置 for j=1:dim % 根据WOA-BP算法来更新种群中所有个体的位置,并求出每个个体的适应度函数值 r3 = rand(); D = abs(pop_position(i,j)-pop_position(r3,j)); pop_position(i,j) = D*cos(c1*2*pi)*pop_position(r3,j)+D*cos(c2*2*pi)*pop_position(best_index,j)+D*cos(c3*2*pi)*rand(); if pop_position(i,j)>x_max % 边界处理 pop_position(i,j) = x_max; elseif pop_position(i,j)<x_min pop_position(i,j) = x_min; end end for j=1:dim % 根据新位置计算适应度函数值并更新最优解个体序号best_index fitness_val_new = fitness_func(pop_position(i,:)); if fitness_val_new<pop_fitness(i) pop_fitness(i) = fitness_val_new; best_index=i; end if fitness_val_new<pop_fitness(best_index) best_index=i; end end end end end % Step 4: 根据更新后的位置计算新的适应度函数值,并根据新的适应度函数值对种群进行排序 for i=1:pop_size % 根据新位置计算适应度函数值并更新最优解个体序号best_index fitness_val_new = fitness_func(pop_position(i,:)); if fitness_val_new<pop_fitness(i) pop_fitness(i) = fitness_val_new; best_index=i; end if fitness_val_new<pop_fitness(best_index) best_index=i; end end [sorted_fit, sorted_index] = sort(pop_fitness); % 排序 % Step 5: 判断是否满足停止条件,如果满足,则输出最优解,否则返回步骤2 if sorted_fit(1)<min_fitness_val % 达到最小误差则停止迭代,输出最优解 best_solution = pop_position(sorted_index(1),:); fprintf('The best solution is:\n'); disp(best_solution); else % 没有达到最小误差,则继续迭代下去 continue; end % 定义适应度函数fitness_func,根据当前权重计算误差值并返回fitness_val function fitness_val=fitness_func(weights) ... (根据权重weights计算误差并返回fitness_val) end ```
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值