%%?I.?清空環(huán)境變量
clear?all
clc

%%?II.?訓(xùn)練集/測試集產(chǎn)生
%%
%?1.?導(dǎo)入數(shù)據(jù)
load?spectra_data.mat

%%
%?2.?隨機產(chǎn)生訓(xùn)練集和測試集
temp?=?randperm（size（NIR1））;??????%打亂60個樣本排序
%?訓(xùn)練集——50個樣本
P_train?=?NIR（temp（1:50）:）‘;??????
T_train?=?octane（temp（1:50）:）‘;
%?測試集——10個樣本
P_test?=?NIR（temp（51:end）:）‘;
T_test?=?octane（temp（51:end）:）‘;
N?=?size（P_test2）;

%%?III.?數(shù)據(jù)歸一化
[p_train?ps_input]?=?mapminmax（P_train01）;
p_test?=?mapminmax（‘a(chǎn)pply‘P_testps_input）;

[t_train?ps_output]?=?mapminmax（T_train01）;??????????????

%%?IV.?BP神經(jīng)網(wǎng)絡(luò)創(chuàng)建、訓(xùn)練及仿真測試
%%
%?1.?創(chuàng)建網(wǎng)絡(luò)
net?=?newff（p_traint_train9）;????%9是隱含層神經(jīng)元的個數(shù)（大家改改測試下結(jié)果影響），連接權(quán)值是3628，講一下怎么計算得到的

%%
%?2.?設(shè)置訓(xùn)練參數(shù)
net.trainParam.epochs?=?1000;???%迭代次數(shù)
net.trainParam.goal?=?1e-3;??????%mse均方根誤差小于這個值訓(xùn)練結(jié)束
net.trainParam.lr?=?0.01;?????????%學(xué)習(xí)率

%%
%?3.?訓(xùn)練網(wǎng)絡(luò)
net?=?train（netp_traint_train）;

%%
%?4.?仿真測試
t_sim?=?sim（netp_test）;?????????%返回10個樣本的預(yù)測值

%%
%?5.?數(shù)據(jù)反歸一化
T_sim?=?mapminmax（‘reverse‘t_simps_output）;???%反歸一化結(jié)果

%%?V.?性能評價
%%
%?1.?相對誤差error
error?=?abs（T_sim?-?T_test）./T_test;

%%
%?2.?決定系數(shù)R^2
R2?=?（N?*?sum（T_sim?.*?T_test）?-?sum（T_sim）?*?sum（T_test））^2?/?（（N?*?sum（（T_sim）.^2）?-?（sum（T_sim））^2）?*?（N?*?sum（（T_test）.^2）?-?（sum（T_test））^2））;?

%%
%?3.?結(jié)果對比
result?=?[T_test‘?T_sim‘?error‘]?????%輸出真實值，預(yù)測值，誤差

%%?VI.?繪圖
figure
plot（1:NT_test‘b:*‘1:NT_sim‘r-o‘）
legend（‘真實值‘‘預(yù)測值‘）
xlabel（‘預(yù)測樣本‘）
ylabel（‘辛烷值‘）
string?=?{‘測試集辛烷值含量預(yù)測結(jié)果對比‘;[‘R^2=‘?num2str（R2）]};
title（string）

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

?????文件???????1706??2017-12-06?19:13??main.m

?????文件?????171497??2017-12-06?19:13??spectra_data.mat

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

???????????????173203????????????????????2