load?data
ldata?=?data.SYSLoad;
LData?=?ldata（:）‘;
[LdataPS]?=?mapminmax（LData01）;%歸一化

ori_X?=?zeros（8000048）;
for?j?=?1:80000
????for?i?=?1:48
???????ori_X（ji）?=?Ldata（j+i-1）;
????end
end

ori_Y?=?zeros（800001）;
for?i?=?1:80000
????ori_Y（i）?=?Ldata（48+6+i）;
end


index?=?randperm（80000）;
train_index?=?index（1:2）;
X?=?ori_X（train_index:）;
Y?=?ori_Y（train_index）;


test_X?=?ori_X（test_index:）;
test_nor_Y?=?ori_Y（test_index）;%歸一化
test_Y?=?mapminmax（‘reverse‘test_nor_YPS）;


%train?dbn
dbn.sizes?=?[24?12];
opts.numepochs?=???50;
opts.batchsize?=?100;
opts.momentum??=???0;
opts.alpha?????=???1;
dbn?=?dbnsetup_bp（dbn?ori_X?opts）;
dbn?=?dbntrain（dbn?ori_X?opts）;

%unfold?dbn?to?nn
nn?=?dbnunfoldtonn（dbn?1）;
nn.activation_function?=?‘sigm‘;


%train?nn
opts.numepochs?=??5000;
opts.batchsize?=?100;
nn?=?nntrain（nn?X?Y?opts）;

[er?bad?ybar]?=?nntest（nn?test_X?test_nor_Y）;


reverseForecastResult?=?mapminmax（‘reverse‘ybarPS）;


%未歸一化數值計算mape
mape?=?sum（abs（test_Y-reverseForecastResult）./test_Y）/2000;
%歸一化數值計算rmse
rmse?=?sqrt（sum（（ybar?-?test_nor_Y）.^2）?/?2000）;?
disp（‘dbn‘）


figure;?%創建一個繪圖窗口
plot（test_nor_Y）;
hold?on
plot（ybar）;
legend（‘Actual?wind?speed‘‘Forecast?wind?speed‘）;
ylabel（{‘Wind?Speed?（m/s）‘}）;
xlabel（{‘Number?of?day‘}）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件??????886166??2016-10-18?17:58??classification.pdf
?????文件????????1377??2019-11-12?15:20??ddbn.m
?????文件??????191390??2016-10-21?11:47??ranzato-06.pdf