clc
clear
close?all
%load?shuju
%load?xcal
%load?xpre
load?ycal
load?ypre
load?cab_1152_jm_GF
load?cab_1152_yz_GaussianFilter
%Xcal=xcal;
Xcal=cab_1152_jm_Sheet3_GaussianFilter;
%Xval=xpre;
Xval=cab_1152_yz_Sheet3_GaussianFilter;
yval=ypre;
m_min=1;
m_max=20;

N?=?size（Xcal1）;?
K?=?size（Xcal2）;?


%?Phase?1:?Projection?operations?for?the?selection?of?candidate?subsets
????
normalization_factor?=?std（Xcal）;
????
for?k?=?1:K
????x?=?Xcal（:k）;
????Xcaln（:k）?=?（x?-?mean（x））?/?normalization_factor（k）;
end


SEL?=?zeros（m_maxK）;
for?k?=?1:K
????SEL（:k）?=?projection（Xcalnkm_max）;
end

????
%?Phase?2:?Evaluation?of?the?candidate?subsets?according?to?the?PRESS?criterion

PRESS?=?Inf*ones（m_maxK）;

for?k?=?1:K
????for?m?=?m_min:m_max
????????var_s