%%?清空環境變量
clc
clear

%%?訓練數據預測數據提取及歸一化
%下載輸入輸出數據
load?date1?input
load?date2?output

%從1到2000間隨機排序
k=rand（1146）;
[mn]=sort（k）;

%找出訓練數據和預測數據
input_train=input（n（1:117）:）‘;?%130組訓練數據，22組預測數據
output_train=output（n（1:117））‘;
input_test=input（n（118:146）:）‘;
output_test=output（n（118:146））‘;

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

%%?BP網絡訓練
%?%初始化網絡結構
net=newff（inputnoutputn10）;

net.trainParam.epochs=10;
net.trainParam.lr=0.1;
net.trainParam.goal=0.00001;

%網絡訓練
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‘-*‘）;
plot

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

?????文件???????1288??2012-09-20?20:08??源碼\BP.asv

?????文件???????1349??2012-09-29?22:55??源碼\BP.m

?????文件???????1425??2010-01-30?20:21??源碼\BP_Hidden.m

?????文件???????1106??2012-09-29?22:54??源碼\date1.mat

?????文件????????820??2012-09-29?22:54??源碼\date2.mat

?????目錄??????????0??2012-09-30?20:31??源碼

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

?????????????????5988????????????????????6