【光伏发电功率预测】蜣螂算法优化变分模态分解结合核极限学习机DBO-VMD-KELM光伏发电功率预测【含Matlab源码 3424期】

✅博主简介：热爱科研的Matlab仿真开发者，修心和技术同步精进，Matlab项目合作可私信。
🍎个人主页：海神之光
🏆代码获取方式：
海神之光Matlab王者学习之路—代码获取方式
⛳️座右铭：行百里者，半于九十。

更多Matlab仿真内容点击👇
Matlab图像处理（进阶版）
路径规划（Matlab）
神经网络预测与分类（Matlab）
优化求解（Matlab）
语音处理（Matlab）
信号处理（Matlab）
车间调度（Matlab）

⛄一、蜣螂算法优化变分模态分解结合核极限学习机简介

摘要：蜣螂优化算法( Dung beetle optimizer, DBO), 是由 Jiankai Xue 等于2022 年提出的一种群体智能优化算法。其灵感来源于蜣螂的生物行为过程，具有寻优能力强，收敛速度快的特点。
1 蜣螂优化算法
众所周知，蜣螂是自然界中一种常见的昆虫，动物的粪便为食，在全世界内分布广泛，扮演着自然界中分解者的角色，对生态系统平衡起着至关重要的作用。蜣螂有一个有趣的习惯，它们会把粪便捏成球，然后把它滚出来，目的是能够尽可能快速、有效地移动粪球，防止被其他蜣螂抢夺。蜣螂的可以利用天体线索(特别是太阳、月亮和偏振光)来导航，让粪球沿着直线滚动，如果完全没有光源(也就是在完全黑暗的环境中)，蜣螂的就不再走直线，而是弯曲的，有时甚至略圆，有很多因素(如风、地面不平)都会导致蜣螂偏离原来的方向，蜣螂在滚粪球的过程如遇到障碍物而无法前进时，通常会爬到粪球上面"跳舞"(包括一系列的旋转和停顿)，决定它们的运动方向。

从蜣螂的习性中观察发现，其获取粪球主要有以下两个目的：①用来产卵和养育下一代；②作为食物。蜣螂会把粪球埋起来，雌性蜣螂会在粪球里产卵，粪球不仅是蜣螂幼虫的发育场所，也是必需的食物。所以，粪球对蜣螂的生存起着不可替代的作用。

本位介绍了一种新的群体智能优化算法------DBO（Dung beetle optimizer）技术，其灵感主要来源于蜣螂的滚球、跳舞、觅食、偷窃、和繁殖等行为。

1.1 结构和算法
根据上面的讨论，蜣螂在滚动过程中需要通过天体线索导航，以保持粪球在直线路径上滚动。为了模拟滚球行为，要求蜣螂在整个搜索空间中沿着给定的方向移动。蜣螂的运动轨迹如图1所示。在图1中，蜣螂利用太阳来导航，其中红色箭头表示的是滚动的方向,同时，我们假设光源的强度也会影响蜣螂的路径。

1.2 计算步骤
DBO 算法作为一种新颖的基于 SI 的优化技术, 主要有六个步骤:
(1) 初始化蜣螂群和 DBO 算法的参数;
(2) 根据目标函数计算出所有目标代理的适应度值;
(3) 更新所有蛲螂的位置;
(4) 判断每个目标代理是否超出边界;
(5) 更新当前最优解及其适应度值;
(6) 重复上述步骤, 直到 t 满足终止准则, 输出全局最优解及其适应度值。

2 蜣螂算法优化变分模态分解结合核极限学习机
蜣螂算法优化变分模态分解结合核极限学习机是一种用于光伏发电功率预测的算法流程。该算法结合了蜣螂算法、变分模态分解和核极限学习机三种算法，通过蜣螂算法优化变分模态分解的过程，将光伏发电功率信号分解成多个本征模态函数，然后利用核极限学习机对这些本征模态函数进行非线性建模，从而实现对光伏发电功率的预测。

具体流程如下：
（1）采集光伏发电功率信号，并将其进行预处理。
（2）利用蜣螂算法优化变分模态分解，将光伏发电功率信号分解成多个本征模态函数。
（3）利用核极限学习机对这些本征模态函数进行非线性建模，得到预测模型。
（4）利用预测模型对未来的光伏发电功率进行预测。

⛄二、部分源代码

clear
clc
close all

load(‘kelm.mat’)
train_simu1kelm=train_simu1;
test_simu1kelm=test_simu1;
load(‘vmd-kelm.mat’)
train_simu1vmd=train_simu1;
test_simu1vmd=test_simu1;
load(‘vmd-bwo-kelm.mat’)
train_simu1bwovmd=train_simu1;
test_simu1bwovmd=test_simu1;

figure
plot(output_train)
hold on
plot(train_simu1kelm)
plot(train_simu1vmd)
plot(train_simu1bwovmd)
legend(‘真实值’,‘预测值’)
title(‘训练集’)
legend(‘真实值’,‘kelm’,‘vmd-kelm’,‘vmd-bwo-kelm’)
xlabel(‘时间/5min’)
ylabel(‘光伏发电功率/kW’)

figure
plot(output_test )
hold on
plot(test_simu1kelm)
plot(test_simu1vmd)
plot(test_simu1bwovmd)
legend(‘真实值’,‘预测值’)
title(‘测试集’)
legend(‘真实值’,‘kelm’,‘vmd-kelm’,‘vmd-bwo-kelm’)
xlabel(‘时间/5min’)
ylabel(‘光伏发电功率/kW’)
%%
reay = output_train; % 真实数据
prey = train_simu1kelm ; % 仿真数据
reay = reshape(reay,1,size(reay,1)*size(reay,2)); % 真实数据
prey = reshape(prey,1,size(prey,1)*size(prey,2)); % 仿真数据
num=length(reay);%统计样本总数
error=prey-reay; %计算误差
mae=sum(abs(error))/num; %计算平均绝对误差
me=sum((error))/num; %计算平均绝对误差
mse=sum(error.error)/num; %计算均方误差
rmse=sqrt(mse); %计算均方误差根
r=min(min(corrcoef(prey,reay)));
R2=rr;
disp(‘训练集’)
disp(‘------------------------kelm----------------------------------’)
disp(['平均绝对误差mae为： ',num2str(mae)])
disp(['平均误差me为： ',num2str(me)])
disp(['均方误差根rmse为： ',num2str(rmse)])
disp(['相关系数R2为： ',num2str(R2)])
reay = output_test; % 真实数据
prey = test_simu1kelm; % 仿真数据
reay = reshape(reay,1,size(reay,1)*size(reay,2)); % 真实数据
prey = reshape(prey,1,size(prey,1)*size(prey,2)); % 仿真数据
num=length(reay);%统计样本总数
error=prey-reay; %计算误差
mae=sum(abs(error))/num; %计算平均绝对误差
me=sum((error))/num; %计算平均绝对误差
mse=sum(error.error)/num; %计算均方误差
rmse=sqrt(mse); %计算均方误差根
r=min(min(corrcoef(prey,reay)));
R2=rr;
disp(‘测试集’)
disp(['平均绝对误差mae为： ',num2str(mae)])
disp(['平均误差me为： ',num2str(me)])
disp(['均方误差根rmse为： ',num2str(rmse)])
disp(['相关系数R2为： ',num2str(R2)])

⛄三、运行结果

⛄四、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1]张炜, 崔娟. 基于麻雀算法优化的核极限学习机在风电功率预测中的应用[J]. 电力科学与工程, 2018, 34(2): 1-7.

[2]师晓龙, 王楠, 李颖, 等. 基于改进的麻雀算法和核极限学习机的风电功率预测[J]. 电力系统保护与控制, 2016, 44(16): 39-45.

[3]殷定庆, 程晓杰, 程肖飞, 等. 基于麻雀算法优化的核极限学习机在风电功率预测中的应用[J]. 计算机与数字工程, 2014, 42(12): 2858-2863

3 备注
简介此部分摘自互联网，仅供参考，若侵权，联系删除

🍅 仿真咨询
1 各类智能优化算法改进及应用
生产调度、经济调度、装配线调度、充电优化、车间调度、发车优化、水库调度、三维装箱、物流选址、货位优化、公交排班优化、充电桩布局优化、车间布局优化、集装箱船配载优化、水泵组合优化、解医疗资源分配优化、设施布局优化、可视域基站和无人机选址优化

2 机器学习和深度学习方面
卷积神经网络（CNN）、LSTM、支持向量机（SVM）、最小二乘支持向量机（LSSVM）、极限学习机（ELM）、核极限学习机（KELM）、BP、RBF、宽度学习、DBN、RF、RBF、DELM、XGBOOST、TCN实现风电预测、光伏预测、电池寿命预测、辐射源识别、交通流预测、负荷预测、股价预测、PM2.5浓度预测、电池健康状态预测、水体光学参数反演、NLOS信号识别、地铁停车精准预测、变压器故障诊断

3 图像处理方面
图像识别、图像分割、图像检测、图像隐藏、图像配准、图像拼接、图像融合、图像增强、图像压缩感知

4 路径规划方面
旅行商问题（TSP）、车辆路径问题（VRP、MVRP、CVRP、VRPTW等）、无人机三维路径规划、无人机协同、无人机编队、机器人路径规划、栅格地图路径规划、多式联运运输问题、车辆协同无人机路径规划、天线线性阵列分布优化、车间布局优化

5 无人机应用方面
无人机路径规划、无人机控制、无人机编队、无人机协同、无人机任务分配

6 无线传感器定位及布局方面
传感器部署优化、通信协议优化、路由优化、目标定位优化、Dv-Hop定位优化、Leach协议优化、WSN覆盖优化、组播优化、RSSI定位优化

7 信号处理方面
信号识别、信号加密、信号去噪、信号增强、雷达信号处理、信号水印嵌入提取、肌电信号、脑电信号、信号配时优化

8 电力系统方面
微电网优化、无功优化、配电网重构、储能配置

9 元胞自动机方面
交通流 人群疏散 病毒扩散 晶体生长

10 雷达方面
卡尔曼滤波跟踪、航迹关联、航迹融合