%%?清空環境變量
close?all
clc;?clear


%%?訓練數據預測數據提取及歸一化
%下載輸入輸出數據
load?traindata1011?A?O?
load?goontest?inputtest_may16?outputtest_may16


inputnum=7;hiddennum=6;outputnum=1;?%4-5-1結構?編碼長度31?7-5-1?46

%訓練數據和預測數據
input_train=A（1:360:）‘;
input_test=inputtest_may16（1:24:）‘;
output_train=O（1:360）‘;
output_test=outputtest_may16（1:24）‘;

%選連樣本輸入輸出數據歸一化
[inputninputps]=mapminmax（input_train）;
[outputnoutputps]=mapminmax（output_train）;

%%?BP網絡訓練
%構建網絡
net=newff（inputnoutputnhiddennum{‘tansig‘‘purelin‘}）;

%網絡進化參數
net.trainParam.epochs=100;?%迭代次數，學習率，目標
net.trainParam.lr=0.1;
net.trainParam.goal=0.00001;
net.trainParam.show=100;
net.trainParam.showWindow=0;

%網絡訓練
net=train（netinputnoutputn）;

%%?BP網絡預測
%預測數據歸一化
inputn_test=mapminmax（‘apply‘input_testinputps）;
?
%網絡預測輸出
an=sim（netinputn_test）;
?
%網絡輸出反歸一化
BPoutput=mapminmax（‘reverse‘anoutputps）;

%%?結果分析

figure（1）
plot（BPoutput‘:og‘）
hold?on
plot（output_test‘-*‘）;
legend（‘預測輸出‘‘期望輸出‘）
title（‘BP網絡預測輸出‘‘fontsize‘12）
ylabel（‘函數輸出（W）‘‘fontsize‘12）
xlabel（‘樣本（h）‘‘fontsize‘12）
%預測誤差
error=BPoutput-output_test;
error1=（output_test-BPoutput）./BPoutput;?%相對誤差

%?figure（2）
%?plot（error‘-*‘）
%?title（‘BP網絡預測誤差‘‘fontsize‘12）
%?ylabel（‘誤差‘‘fontsize‘12）
%?xlabel（‘樣本‘‘fontsize‘12）

figure（3）
plot（error1‘-*‘）;
title（‘BP神經網絡預測誤差百分比‘）
ylabel（‘相對誤差‘‘fontsize‘12）
xlabel（‘樣本（h）‘‘fontsize‘12）

figure（4）
hist（error1）;
title（‘神經網絡預測誤差頻率分布直方圖‘）;

MAE=（sum（abs（error1）））/24?%絕對平均誤差?24對誤差
RMSE=sqrt（（sum（error1.^2））/24）%RMSE?均方根誤差公式
%?errorsum=sum（abs（error））

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1813??2014-01-07?21:59??BP?wind?prediction\bp.m

?????文件???????1737??2014-01-07?21:57??BP?wind?prediction\BPNWP.m

?????文件????????874??2013-10-12?14:59??BP?wind?prediction\goontest.mat

?????文件????????974??2013-10-16?15:49??BP?wind?prediction\goontestNWP.mat

?????文件???????6430??2013-10-11?21:12??BP?wind?prediction\traindata1011.mat

?????文件???????7287??2013-10-16?15:38??BP?wind?prediction\traindata1011NWP.mat

?????目錄??????????0??2014-01-14?22:52??BP?wind?prediction

-----------?---------??----------?-----??----

????????????????19115????????????????????7