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

%%?II.?導(dǎo)入數(shù)據(jù)
load?spectra_data;

%%?III.?隨機(jī)劃分訓(xùn)練集與測(cè)試集
temp?=?randperm（size（NIR?1））;
%?temp?=?1:60;
%%
%?1.?訓(xùn)練集——50個(gè)樣本
P_train?=?NIR（temp（1:50）:）;
T_train?=?octane（temp（1:50）:）;
%%
%?2.?測(cè)試集——10個(gè)樣本
P_test?=?NIR（temp（51:end）:）;
T_test?=?octane（temp（51:end）:）;

%%?IV.?主成分分析
%%
%?1.?主成分貢獻(xiàn)率分析
[PCALoadingsPCAScoresPCAVar]?=?princomp（NIR）;
figure
percent_explained?=?100?*?PCAVar?/?sum（PCAVar）;
pareto（percent_explained）
xlabel（‘主成分‘）
ylabel（‘貢獻(xiàn)率（%）‘）
title（‘主成分貢獻(xiàn)率‘）

%%
%?2.?第一主成分vs.第二主成分
[PCALoadingsPCAScoresPCAVar]?=?princomp（P_train）;
figure
plot（PCAScores（:1）PCAScores（:2）‘r+‘）
hold?on
[PCALoadings_testPCAScores_testPCAVar_test]?=?princomp（P_test）;
plot（PCAScores_test（:1）PCAScores_test（:2）‘o‘）
xlabel（‘1st?Principal?Component‘）
ylabel（‘2nd?Principal?Component‘）
legend（‘Training?Set‘‘Testing?Set‘‘location‘‘best‘）

%%?V.?主成分回歸模型
%%
%?1.?創(chuàng)建模型
k?=?4;
betaPCR?=?regress（T_train-mean（T_train）PCAScores（:1:k））;
betaPCR?=?PCALoadings（:1:k）?*?betaPCR;
betaPCR?=?[mean（T_train）-mean（P_train）?*?betaPCR;betaPCR];
%%
%?2.?預(yù)測(cè)擬合
N?=?size（P_test1）;
T_sim?=?[ones（N1）?P_test]?*?betaPCR;

%%?VI.?結(jié)果分析與繪圖
%%
%?1.?相對(duì)誤差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é)果對(duì)比
result?=?[T_test?T_sim?error]

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

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----
?????目錄???????????0??2018-11-27?11:07??Class_13_Code\
?????文件????????1813??2018-11-27?11:08??Class_13_Code\PCR_Demo.m
?????文件????????1302??2015-11-22?20:19??Class_13_Code\PLS_Demo.m
?????目錄???????????0??2015-11-22?20:19??Class_13_Code\html\
?????文件???????10397??2015-11-22?20:05??Class_13_Code\html\PCR_Demo.html
?????文件????????3893??2015-11-22?20:05??Class_13_Code\html\PCR_Demo.png
?????文件????????4596??2015-11-22?20:05??Class_13_Code\html\PCR_Demo_01.png
?????文件????????5188??2015-11-22?20:05??Class_13_Code\html\PCR_Demo_02.png
?????文件????????8450??2015-11-22?20:05??Class_13_Code\html\PCR_Demo_03.png
?????文件????????8883??2015-11-22?20:19??Class_13_Code\html\PLS_Demo.html
?????文件????????3217??2015-11-22?20:19??Class_13_Code\html\PLS_Demo.png
?????文件????????4556??2015-11-22?20:19??Class_13_Code\html\PLS_Demo_01.png
?????文件????????7451??2015-11-22?20:19??Class_13_Code\html\PLS_Demo_02.png
?????文件??????171497??2010-10-14?20:24??Class_13_Code\spectra_data.mat