【预测模型】基于狼群算法优化BP神经网络实现预测matlab源码

一、简介

BP神经网络具有很好的拟合非线性函数的能力,但传统BP算法收敛慢,容易陷入局部极小值.不少研究人员从BP神经网络的初始权值和阈值着手,期望通过对初始权值和阈值的优化来提高BP神经网络的性能.对此提出一种改进BP神经网络的算法,使用狼群算法优化BP神经网络的初始权值和阈值,用数据训练BP神经网络后预测函数的输出.最后通过一个非线性函数来验证所提出模型的有效性.

二、源代码

```matlab

%% WPA-BP神经网络

%% 清空环境

clear all;

close all;

clc;

data=xlsread('Test.xlsx','Sheet1','A2:G46');%读取数据

[M,NN]=size(data);%M是样本数,N是变量数

inputn=data(:,1:NN-1)';%训练样本输入

outputn=data(:,NN)';%训练样本输出

%节点个数

inputnum=NN-1;%输入个数

hiddennum=12;%隐藏层

outputnum=1;%输出个数

%构建网络

net=newff(inputn,outputn,hiddennum);

%% 狼群算法参数设置

tic

time = toc;

disp(['用时=  ',num2str(time)])

x=zbest;

%% 把最优初始阀值权值赋予网络预测

%% BP网络训练

%网络进化参数

net.trainParam.epochs=100;%bp神经网络的代数

net.trainParam.lr=0.1;

%net.trainParam.goal=0.00001;

%网络训练

[net,per2]=train(net,inputn,outputn);

%% BP网络预测

% load data inputn_test output_test

inputn_test=inputn;%测试数据输入

output_test=outputn;%测试数据输出

an=sim(net,inputn_test);%预测

error=(an-output_test);%误差

error_xiangdui=(an-output_test)./output_test;%相对误差

figure(1)

plot(error_xiangdui,'ko-')

ylabel('相对误差值')

title('狼群优化BP神经网络')

figure(2)

plot(error,'ro-')

ylabel('误差值')

title('狼群优化BP神经网络')

figure(3)

plot(output_test,'r + -');

hold on;

plot(an,'b * -')

hold off

legend('真实值','预测值');

title('狼群优化BP神经网络')

三、运行结果

在这里插入图片描述

点击并拖拽以移动wAAACH5BAEKAAAALAAAAAABAAEAAAICRAEAOw==

四、参考文献

[1]李常洪, and 赵伟栋. "基于狼群算法优化的BP神经网络." 科技创新与生产力 000.001(2016):56-58.

  • 0
    点赞
  • 5
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
以下是MATLAB实现基于BP神经网络灰狼优化算法的多输入寻优的完整代码: ``` function [xopt,fopt]=GWO_BP(Function_name,ub,lb,dim,SearchAgents_no,Max_iter) % 输入: % Function_name: 测试函数名 % ub: 自变量上界 % lb: 自变量下界 % dim: 自变量维度 % SearchAgents_no: 狼群数量 % Max_iter: 最大迭代次数 % 输出: % xopt: 最优解 % fopt: 最优值 % BP神经网络初始化 net=feedforwardnet([10 5 1]); % 三层前馈神经网络 net.trainFcn='trainlm'; % Levenberg-Marquardt算法 net.trainParam.showWindow=0; % 不显示窗口 net.trainParam.showCommandLine=0; % 不显示命令行 net.trainParam.epochs=100; % 训练次数 net=init(net); % 灰狼优化初始化 Alpha_pos=zeros(1,dim); % Alpha狼位置 Alpha_score=inf; % Alpha狼得分 Beta_pos=zeros(1,dim); % Beta狼位置 Beta_score=inf; % Beta狼得分 Delta_pos=zeros(1,dim); % Delta狼位置 Delta_score=inf; % Delta狼得分 Positions=initialization(SearchAgents_no,dim,ub,lb); % 狼群位置 Convergence_curve=zeros(1,Max_iter); % 收敛曲线 % 灰狼优化主循环 for iter=1:Max_iter for i=1:size(Positions,1) % 计算适应度 fitness=calculate_fitness(Positions(i,:),net,Function_name); % 更新Alpha、Beta、Delta狼 if fitness<Alpha_score Alpha_score=fitness; Alpha_pos=Positions(i,:); end if fitness>Alpha_score && fitness<Beta_score Beta_score=fitness; Beta_pos=Positions(i,:); end if fitness>Alpha_score && fitness>Beta_score && fitness<Delta_score Delta_score=fitness; Delta_pos=Positions(i,:); end end % 更新参数a和A a=2-iter*((2)/Max_iter); A=2*a*rand()-a; % 更新狼群位置 for i=1:size(Positions,1) D_alpha=abs(A*Alpha_pos-Positions(i,:)); D_beta=abs(A*Beta_pos-Positions(i,:)); D_delta=abs(A*Delta_pos-Positions(i,:)); X1=Alpha_pos-A*D_alpha; X2=Beta_pos-A*D_beta; X3=Delta_pos-A*D_delta; Positions(i,:)=mean([X1;X2;X3]); end % 更新神经网络参数 for i=1:size(Positions,1) net=configure(net,Positions(i,:)',zeros(dim,1)'); net=train(net,Positions(i,:)',fitness); end % 更新收敛曲线 Convergence_curve(iter)=Alpha_score; end % 输出结果 xopt=Alpha_pos; fopt=Alpha_score; % 绘制收敛曲线 figure; plot(Convergence_curve); xlabel('迭代次数'); ylabel('最优值'); title('收敛曲线'); end % 狼群初始化函数 function Positions=initialization(SearchAgents_no,dim,ub,lb) Positions=rand(SearchAgents_no,dim).*(ub-lb)+lb; end % 计算适应度函数 function fitness=calculate_fitness(Position,net,Function_name) if strcmp(Function_name,'BP') % 神经网络预测 % ... % 计算误差 % ... % 适应度为误差的相反数 fitness=-error; else % 其他测试函数 % ... end end ``` 在上述代码中,我们使用了feedforwardnet函数初始化了一个三层前馈神经网络,并使用trainlm算法训练网络。在calculate_fitness函数中,我们对输入的位置进行神经网络预测,并计算误差,最终将误差的相反数作为适应度。在主循环中,我们按照灰狼优化算法的步骤逐步更新了狼群位置和参数,同时利用每个位置对神经网络进行训练。最终输出最优解和最优值,并绘制了收敛曲线。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

Matlab科研辅导帮

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值