SPA连续投影算法

SPA连续投影算法

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）‘;

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

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

     文件     119789  2015-12-31 16:44  shuju.mat

     文件       2742  2016-01-12 17:39  SSPPAA.m

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

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

               124463                    4