今天为大家带来一期基于DBO-TCN的电力负荷预测。不得不说，TCN在时序数据中的预测效果在一定程度上要优于LSTM。
原理详解
文章对TCN网络的五个参数进行优化，分别是：
正则化参数，学习率，滤波器个数，滤波器大小，区块数。这五个参数属于TCN网络中比较重要的几个参数，参数的选择会影响最后的精度。
采用蜣螂优化TCN，参数设置范围分别是：
|-------|------------------|
| 正则化参数 | [0.001 0.1] |
| 学习率 | [ 0.0001 0.01] |
| 滤波器个数 | [20 100] |
| 滤波器大小 | [2 10] |
| 区块数 | [2 10] |
将DBO种群数设置为：30，迭代次数设置为30。
数据准备
对电力负荷数据进行处理
本次数据包含负荷值，温度，湿度，风速，压强，降水量，能见度，水汽压和体感温度，部分数据截图如下：
选取1200个样本作为训练集，每个样本组成为：当天24个小时的全部数据，因此训练集的输入数据大小为1200*216，其中216=24*9，24代表24个小时，9代表9个特征。训练集的输出数据大小为：1200*1。1代表未来一小时的负荷值。
选取300个样本作为测试集，同理，测试集的输入数据大小为300*216，训练集的输出数据大小为：300*1。
结果展示
采用TCN网络对电力负荷数据进行训练和预测
TCN的预测结果如下：
可以看到，未优化的TCN网络预测效果不是很理想，但是相比于上一期给出的LSTM的预测效果已经好很多了。
DBO-TCN模型的预测效果如下：
误差对比图如下：
可以看到DBO-TCN预测效果有了明显提升，DBO-TCN的MSE误差为13.821，相比于未优化TCN的23.238有了很大提升！
DBN-TCN的回归拟合图：
误差直方图：
部分代码
cs
% 网络训练
net0 = trainNetwork(inputn,outputn,lgraph,options0 );
an0 = net0.predict(inputn_test);
test_simu0=mapminmax('reverse',an0,outputps); %把仿真得到的数据还原为原始的数量级
%误差指标
error0 = output_test - test_simu0;
mse0=mse(output_test,test_simu0)
%% 标准TCN神经网络作图
figure
plot(output_test,'b-','markerfacecolor',[0.5,0.5,0.9],'MarkerSize',6)
hold on
plot(test_simu0,'r--','MarkerSize',6)
title(['TCN的mse误差：',num2str(mse0)])
legend('真实y','预测的y')
xlabel('样本数')
ylabel('负荷值')
box off
set(gcf,'color','w')
%% 调用DBO优化TCN
disp('调用DBO优化TCN......')
% DBO优化参数设置
SearchAgents = 20; % 种群数量 30
Max_iterations = 30; % 迭代次数 30
lowerbound = [0.001 0.0001 20 2 2]; %五个参数的下限分别是正则化参数，学习率，滤波器个数:numFilters，滤波器大小:filterSize，区块数:numBlocks
upperbound = [0.1 0.01 100 10 10]; %五个参数的上限
dimension = length(lowerbound);%数量，即要优化的TCN参数个数
[fMin,Best_pos,Convergence_curve,bestnet] = DBOforTCN(SearchAgents,Max_iterations,lowerbound,upperbound,dimension,inputn,outputn,inputn_test,outputps,output_test,numFeatures,outputSize);
L2Regularization = Best_pos(1,1); % 最佳L2正则化系数
InitialLearnRate = Best_pos(1,2) ;% 最佳初始学习率
numFilters = fix(Best_pos(1,3));
filterSize = fix(Best_pos(1,4));
numBlocks = fix(Best_pos(1,5));
disp(['最优参数：',num2str(Best_pos)])
disp('优化结束，将最佳net导出并用于测试......')
setdemorandstream(pi);
%% 对测试集的测试
an = bestnet.predict(inputn_test);
test_simu = mapminmax('reverse',an,outputps);
error = test_simu-output_test;
msee = mse(output_test,test_simu);
figure
hold on
plot(test_simu,'g')
plot(output_test,'b')
legend('DBO-TCN预测值','实际值')
ylabel('预测结果')
xlabel('预测样本')
title(['DBO-TCN的mse误差：',num2str(msee)])
box off
set(gcf,'color','w')
figure
plot(abs(error),'g-*')
hold on
plot(abs(error0),'-or')
title('预测误差对比图','fontsize',12)
xlabel('预测样本','fontsize',12)
ylabel('误差绝对值','fontsize',12)
legend('DBO-TCN','TCN')
%% 回归图与误差直方图
figure;
plotregression(test_simu,output_test,['优化后回归图']);
set(gcf,'color','w')
figure;
ploterrhist(test_simu-output_test,['误差直方图']);
set(gcf,'color','w')
%% 打印出评价指标
% 预测结果评价
ae= abs(test_simu-output_test);
rmse = (mean(ae.^2)).^0.5;
mse = mean(ae.^2);
mae = mean(ae);
mape = mean(ae./test_simu);
R = corr(output_test,test_simu);
R2 = 1 - norm(output_test - test_simu)^2 / norm(output_test-mean(output_test ))^2;
disp('预测结果评价指标：')
disp(['RMSE = ', num2str(rmse)])
disp(['MSE = ', num2str(mse)])
disp(['MAE = ', num2str(mae)])
disp(['MAPE = ', num2str(mape)])
disp(['相关系数R = ', num2str(R)])
disp(['决定系数R^2为： ',num2str(R2)])代码获取
完整代码获取，后台回复关键词：
DBOTCN