%%?I.?清空環境變量
clear?all
clc
%%?II.?訓練集/測試集產生
%%
%?1.?導入數據
load?spectra_data.mat
?
%%
%?2.?隨機產生訓練集和測試集
temp?=?randperm（size（NIR1））;
%?訓練集――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.?數據歸一化
[p_train?ps_input]?=?mapminmax（P_train01）;
p_test?=?mapminmax（‘apply‘P_testps_input）;
[t_train?ps_output]?=?mapminmax（T_train01）;
%%?IV.?BP神經網絡創建、訓練及仿真測試
%%
%?1.?創建網絡
net?=?newff（p_traint_train9）;
?
%%
%?2.?設置訓練參數
net.trainParam.epochs?=?1000;
net.trainParam.goal?=?1e-3;
net.trainParam.lr?=?0.01;
?
%%
%?3.?訓練網絡
net?=?train（netp_traint_train）;
?
%%
%?4.?仿真測試
t_sim?=?sim（netp_test）;
?
%%
%?5.?數據反歸一化
T_sim?=?mapminmax（‘reverse‘t_simps_output）;
?
%%?V.?性能評價
%%
%?1.?相對誤差error
error?=?abs（T_sim?-?T_test）./T_test;
?
%%
%?2.?決定系數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.?結果對比
result?=?[T_test‘?T_sim‘?error‘]
%%?VI.?繪圖
figure
plot（1:NT_test‘b:*‘1:NT_sim‘r-o‘）
legend（‘真實值‘‘預測值‘）
xlabel（‘預測樣本‘）
ylabel（‘辛烷值‘）
string?=?{‘測試集辛烷值含量預測結果對比‘;[‘R^2=‘?num2str（R2）]};
title（string）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-10-16?15:40??Prediction?of?gasoline?octane?number\
?????文件????????1432??2018-10-16?14:48??Prediction?of?gasoline?octane?number\ELM.m
?????文件????????3139??2018-10-16?15:01??Prediction?of?gasoline?octane?number\main_2009a.m
?????文件????????3128??2018-10-16?15:01??Prediction?of?gasoline?octane?number\main_2014a.m
?????文件??????171497??2010-10-14?20:24??Prediction?of?gasoline?octane?number\spectra_data.mat
?????文件????????2900??2018-10-16?15:06??Prediction?of?gasoline?octane?number\對汽油辛烷值含量做預測.txt
?????目錄???????????0??2018-10-16?15:59??Prediction?of?gasoline?octane?number\測試\
?????文件?????????362??2018-10-16?15:22??Prediction?of?gasoline?octane?number\測試\result0.txt
?????文件???????54270??2018-10-16?15:21??Prediction?of?gasoline?octane?number\測試\result0汽油辛烷值含量做預測.png
?????文件???????44879??2018-10-16?15:52??Prediction?of?gasoline?octane?number\測試\result1.png
?????文件?????????330??2018-10-16?15:52??Prediction?of?gasoline?octane?number\測試\result1.txt
?????文件???????52817??2018-10-16?15:53??Prediction?of?gasoline?octane?number\測試\result2.png
?????文件?????????330??2018-10-16?15:54??Prediction?of?gasoline?octane?number\測試\result2.txt
?????文件???????26252??2018-10-16?15:55??Prediction?of?gasoline?octane?number\測試\result3.png
?????文件?????????330??2018-10-16?15:55??Prediction?of?gasoline?octane?number\測試\result3.txt
?????文件???????29152??2018-10-16?15:56??Prediction?of?gasoline?octane?number\測試\result4.png
?????文件?????????330??2018-10-16?15:56??Prediction?of?gasoline?octane?number\測試\result4.txt
?????文件???????30145??2018-10-16?15:57??Prediction?of?gasoline?octane?number\測試\result5.png
?????文件?????????332??2018-10-16?15:57??Prediction?of?gasoline?octane?number\測試\result5.txt
?????文件???????25289??2018-10-16?15:58??Prediction?of?gasoline?octane?number\測試\result6.png
?????文件?????????330??2018-10-16?15:58??Prediction?of?gasoline?octane?number\測試\result6.txt
?????文件?????4577709??2018-10-16?15:16??Prediction?of?gasoline?octane?number\測試\Video_2018-10-16_151626.wmv