连续投影法降维spa代码

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