连续投影算法SAP-matlab程序

连续投影算法，matlab程序可用，直接将数据放入excel就可运行。

function chain = projection（XkM）

% Projections routine for the Successive Projections Algorithm using the
% built-in QR function of Matlab
%
% chain = projections（XkM）
%
% X --> Matrix of predictor variables （# objects N x # variables K）
% k --> Index of the initial column for the projection operations
% M --> Number of variables to include in the chain
%
% chain --> Index set of the variables resulting from the projection operations

X_projected = X;

norms = sum（X_projected.^2）;    % Square norm of each column vector
norm_max = max（norms）; % Norm of the “largest“ column vector

X_projected（:k） = X_projected（:k）*2*norm_max/norms（k）; % Scales the kth column so that it becomes the “largest“ column

[dummy1dummy2order] = qr（X_projected0）; 
chain = order（1:M）‘;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        193  2017-10-31 20:05  SPA\boduan_shang.mat

     文件        799  2013-05-26 18:55  SPA\projection.m

     文件       2745  2019-01-02 19:01  SPA\SPA.m

     文件       1133  2007-02-18 13:40  SPA\validation.m

     目录          0  2019-01-02 19:01  SPA

----------- ---------  ---------- -----  ----

                 4870                    5