基于白鲸优化算法BWO优化的VMD-KELM光伏发电短期功率预测MATLAB代码（含详细算法介绍）

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=r*r; 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=r*r; disp('测试集') disp(['平均绝对误差mae为： ',num2str(mae)]) disp(['平均误差me为： ',num2str(me)]) disp(['均方误差根rmse为： ',num2str(rmse)]) disp(['相关系数R2为： ',num2str(R2)]) %% reay = output_train; % 真实数据 prey = train_simu1vmd ; % 仿真数据 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=r*r; disp('训练集') disp('------------------------vmd-kelm----------------------------------') disp(['平均绝对误差mae为： ',num2str(mae)]) disp(['平均误差me为： ',num2str(me)]) disp(['均方误差根rmse为： ',num2str(rmse)]) disp(['相关系数R2为： ',num2str(R2)]) reay = output_test; % 真实数据 prey = test_simu1vmd; % 仿真数据 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=r*r; disp('测试集') disp(['平均绝对误差mae为： ',num2str(mae)]) disp(['平均误差me为： ',num2str(me)]) disp(['均方误差根rmse为： ',num2str(rmse)]) disp(['相关系数R2为： ',num2str(R2)]) %% reay = output_train; % 真实数据 prey = train_simu1bwovmd ; % 仿真数据 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=r*r; disp('训练集') disp('------------------------bwo-vmd-kelm----------------------------------') disp(['平均绝对误差mae为： ',num2str(mae)]) disp(['平均误差me为： ',num2str(me)]) disp(['均方误差根rmse为： ',num2str(rmse)]) disp(['相关系数R2为： ',num2str(R2)]) reay = output_test; % 真实数据 prey = test_simu1bwovmd; % 仿真数据 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=r*r; disp('测试集') disp(['平均绝对误差mae为： ',num2str(mae)]) disp(['平均误差me为： ',num2str(me)]) disp(['均方误差根rmse为： ',num2str(rmse)]) disp(['相关系数R2为： ',num2str(R2)])