【预测模型-ELAMN预测】基于哈里斯鹰算法优化ELMAN神经网络实现数据回归预测matlab代码

1 简介风能,作为一种重要,有潜力,无污染,可再生、可持续的能源,已经成为全球发电最为迅速的能源之一,越来越受到世界各国的青睐。近年来,为缓解能源短缺问题,改善环境,实现经济乃至人类的可持续发展,世界各国纷纷大力发展风能资源。然而,在实际操作中,风能固有的波动性和间歇性通常会增加风能入网的难度,致使风电场输出功率不稳定,严重影响输出的电能质量,大大增加了电网安全性和稳定性运行风险,因此,风速预测在风电相关工作中至关重要,准确、可靠的预测结果不仅有利于调度人员事先掌握风电场的功率变化情况,及时制定调度运行计划,提高能量转换效率、降低风险,增加发电量等,同时也有利于风电并网稳定运行和有效消纳,对可能存在影响电网安全稳定运行的风险及时预警,从而避免风电功率随机波动造成电能损耗甚至电网崩溃。近年来,国内外学者进行了大量风速预测相关研究,风速预测水平得到一定程度的提升。单一预测模型虽简单容易实施,但其预测精度往往较低,难以满足风电场发电的需求。相比之下,基于优化算法和数据分解法的混合预测模型对风速预测性能有很大的提升。因此,这些模型往往会降低预测结果的稳定性,对风电场发电以及风电并网的安全性和稳定性带来巨大的挑战。针对如上问题,本文构建出一个基于哈里斯鹰算法和Elman神经网络的混合预测模型本文所提出的混合预测模型不仅可以降低预测误差,提高风速预测精度,增强预测结果的稳定性,同时对风电场发电的完善,发展以及风电并网安全管理也大有裨益。Elman 神经网络是一种局部反馈递归神经网络包括输入层、隐藏层、承接层和输出层由于其良好的动态记忆和时变能力自1990年Elman首次提出以来Elman神经网络在时间序列预测中得到了广泛的应用但由于随机选择初始值和阈值且采用梯度下降法寻优其网络学习速度较慢且预测的精度比较低。本文利用 HHO优化 Elman神经网络的初始权值及阈值提高了 Elman 神经网络的预测精度。第一步 初始化 Elman 网络所需参数并新建 Elman 网 络netnewelmminmaxinput111tansigpurelintraingdx其中 input 是训练集数据tansig 是隐藏层的激活函数purelin是输出层激活函数traingdx是梯度下降函数第二步 初始化HHO各参数包括鹰的数量m 、最大迭代次数Maxiter参数的上界ub及下界lb第三步 不断迭代哈里斯鹰位置优化网络直到迭代结束形成最优个体其中最优个体包括权值w1、w2及w3阈值b1和b2第四步 Elman 神经网络解析最优个体利用解析后的权值及阈值训练Elman神经网络并预测风速数据。2 部分代码function [sFeat,Sf,Nf,curve] jBHHO(feat,label,N,max_Iter,HO)beta 1.5;ub 1;lb 0;fun jFitnessFunction;dim size(feat,2);X zeros(N,dim);for i 1:Nfor d 1:dimif rand() 0.5X(i,d) 1;endendendfitR inf;fit zeros(1,N);Y zeros(1,dim);Z zeros(1,dim);curve inf;t 1;%---Iteration start-------------------------------------------------while t max_Iterfor i 1:Nfit(i) fun(feat,label,X(i,:),HO);if fit(i) fitRfitR fit(i);Xrb X(i,:);endendX_mu mean(X,1);for i 1:NE0 -1 2 * rand();E 2 * E0 * (1 - (t / max_Iter));if abs(E) 1q rand();if q 0.5k randi([1,N]);r1 rand();r2 rand();for d 1:dimXn X(k,d) - r1 * abs(X(k,d) - 2 * r2 * X(i,d));S 1 / (1 exp(-Xn));if rand() SX(i,d) 1;elseX(i,d) 0;endendelseif q 0.5r3 rand();r4 rand();for d 1:dimXn (Xrb(d) - X_mu(d)) - r3 * (lb r4 * (ub - lb));S 1 / (1 exp(-Xn));if rand() SX(i,d) 1;elseX(i,d) 0;endendendelseif abs(E) 1J 2 * (1 - rand());r rand();if r 0.5 abs(E) 0.5for d 1:dimDX Xrb(d) - X(i,d);Xn DX - E * abs(J * Xrb(d) - X(i,d));S 1 / (1 exp(-Xn));if rand() SX(i,d) 1;elseX(i,d) 0;endendelseif r 0.5 abs(E) 0.5for d 1:dimDX Xrb(d) - X(i,d);Xn Xrb(d) - E * abs(DX);S 1 / (1 exp(-Xn));if rand() SX(i,d) 1;elseX(i,d) 0;endendelseif r 0.5 abs(E) 0.5LF jLevyDistribution(beta,dim);for d 1:dimYn Xrb(d) - E * abs(J * Xrb(d) - X(i,d));S 1 / (1 exp(-Yn));if rand() SY(d) 1;elseY(d) 0;endZn Y(d) rand() * LF(d);S 1 / (1 exp(-Zn));if rand() SZ(d) 1;elseZ(d) 0;endendfitY fun(feat,label,Y,HO);fitZ fun(feat,label,Z,HO);if fitY fit(i)fit(i) fitY;X(i,:) Y;endif fitZ fit(i)fit(i) fitZ;X(i,:) Z;endelseif r 0.5 abs(E) 0.5LF jLevyDistribution(beta,dim);for d 1:dimYn Xrb(d) - E * abs(J * Xrb(d) - X_mu(d));S 1 / (1 exp(-Yn));if rand() SY(d) 1;elseY(d) 0;endZn Y(d) rand() * LF(d);S 1 / (1 exp(-Zn));if rand() SZ(d) 1;elseZ(d) 0;endendfitY fun(feat,label,Y,HO);fitZ fun(feat,label,Z,HO);if fitY fit(i)fit(i) fitY;X(i,:) Y;endif fitZ fit(i)fit(i) fitZ;X(i,:) Z;endendendendcurve(t) fitR;fprintf(\nIteration %d Best (BHHO) %f,t,curve(t))t t 1;endPos 1:dim;Sf Pos(Xrb 1);Nf length(Sf);sFeat feat(:,Sf);endfunction LF jLevyDistribution(beta,dim)nume gamma(1 beta) * sin(pi * beta / 2);deno gamma((1 beta) / 2) * beta * 2 ^ ((beta - 1) / 2);sigma (nume / deno) ^ (1 / beta);u randn(1,dim) * sigma;v randn(1,dim);step u ./ abs(v) .^ (1 / beta);LF 0.01 * step;end3 仿真结果4 参考文献[1]朱昶胜, 康亮河, and 冯文芳. 基于自适应鲸鱼优化算法结合Elman神经网络的股市收盘价预测算法. 计算机应用 40.5(2020):9.[2]郭剑鹰, 夏李玮, and 刘华军. 基于Elman神经网络预测的机动目标跟踪滤波算法. 航天电子对抗 34.3(2018):6.[3]郭雨鑫, 刘升, 高文欣, 张磊. (2021). 多策略改进哈里斯鹰优化算法. 微电子学与计算机, 38(7), 7.​