稀疏保持投影matlab代码

运用稀疏保持投影进行特征提取（已应用于人脸识别领域）

function [disc_setdisc_valueMean_Image]=Eigenface_f（Train_SETEigen_NUM）

% the magnitude of eigenvalues of this function is corrected right !!!!!!!!!
% Centralized PCA
[NNTrain_NUM]=size（Train_SET）;

if NN<=Train_NUM % for small sample size case
    
   Mean_Image=mean（Train_SET2）;  
   Train_SET=Train_SET-Mean_Image*ones（1Train_NUM）;
   R=Train_SET*Train_SET‘/（Train_NUM-1）;
   
   [VS]=Find_K_Max_Eigen（REigen_NUM）;
   disc_value=S;
   disc_set=V;

else % for small sample size case
    
   Mean_Image=mean（Train_SET2）;  
   Train_SET=Train_SET-Mean_Image*ones（1Train_NUM）;

  R=Train_SET‘*Train_SET/（Train_NUM-1）;
  
  [VS]=Find_K_Max_Eigen（REigen_NUM）;
  disc_value=S;
  disc_set=zeros（NNEigen_NUM）;
  
  Train_SET=Train_SET/sqrt（Train_NUM-1）;
  for k=1:Eigen_NUM
    disc_set（:k）=（1/sqrt（disc_value（k）））*Train_SET*V（:k）;
  end

end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function [Eigen_VectorEigen_Value]=Find_K_Max_Eigen（MatrixEigen_NUM）

[NNNN]=size（Matrix）;
[VS]=eig（Matrix）; %Note this is equivalent to; [VS]=eig（StSL）; also equivalent to [VS]=eig（SnSt）; %

S=diag（S）;
[Sindex]=sort（S）;

Eigen_Vector=zeros（NNEigen_NUM）;
Eigen_Value=zeros（1Eigen_NUM）;

p=NN;
for t=1:Eigen_NUM
    Eigen_Vector（:t）=V（:index（p））;
    Eigen_Value（t）=S（p）;
    p=p-1;
end

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

     文件       1521  2014-08-10 08:52  Orl_SPP_SRC 代码\Eigenface_f.m

     文件       1657  2014-08-10 08:52  Orl_SPP_SRC 代码\eigen_reconstruction.m

     文件        437  2014-08-10 08:52  Orl_SPP_SRC 代码\Find_K_Max_Gen_Eigen.m

     文件       6148  2014-08-10 08:52  Orl_SPP_SRC 代码\inexact_alm_rpca\inexact_alm_rpca\.DS_Store

     文件        560  2014-08-10 08:52  Orl_SPP_SRC 代码\inexact_alm_rpca\inexact_alm_rpca\choosvd.m

     文件       2885  2014-08-10 08:52  Orl_SPP_SRC 代码\inexact_alm_rpca\inexact_alm_rpca\inexact_alm_rpca.m

     文件        202  2014-08-10 08:52  Orl_SPP_SRC 代码\inexact_alm_rpca\inexact_alm_rpca\PROPACK\Afunc.m

     文件        224  2014-08-10 08:52  Orl_SPP_SRC 代码\inexact_alm_rpca\inexact_alm_rpca\PROPACK\AtAfunc.m

     文件        232  2014-08-10 08:52  Orl_SPP_SRC 代码\inexact_alm_rpca\inexact_alm_rpca\PROPACK\Atransfunc.m

     文件        987  2014-08-10 08:52  Orl_SPP_SRC 代码\inexact_alm_rpca\inexact_alm_rpca\PROPACK\bdsqr.m

     文件      74546  2014-08-10 08:52  Orl_SPP_SRC 代码\inexact_alm_rpca\inexact_alm_rpca\PROPACK\bdsqr.mexglx

     文件      25208  2014-08-10 08:52  Orl_SPP_SRC 代码\inexact_alm_rpca\inexact_alm_rpca\PROPACK\bdsqr.mexsg

     文件      26109  2014-08-10 08:52  Orl_SPP_SRC 代码\inexact_alm_rpca\inexact_alm_rpca\PROPACK\bdsqr.mexsg64

     文件      95852  2014-08-10 08:52  Orl_SPP_SRC 代码\inexact_alm_rpca\inexact_alm_rpca\PROPACK\bdsqr.mexsol

     文件      55808  2014-08-10 08:52  Orl_SPP_SRC 代码\inexact_alm_rpca\inexact_alm_rpca\PROPACK\bdsqr.mexw32

     文件       2102  2014-08-10 08:52  Orl_SPP_SRC 代码\inexact_alm_rpca\inexact_alm_rpca\PROPACK\bdsqr_mex.c

     文件        293  2014-08-10 08:52  Orl_SPP_SRC 代码\inexact_alm_rpca\inexact_alm_rpca\PROPACK\Cfunc.m

     文件       1504  2014-08-10 08:52  Orl_SPP_SRC 代码\inexact_alm_rpca\inexact_alm_rpca\PROPACK\compute_int.m

     文件        445  2014-08-10 08:52  Orl_SPP_SRC 代码\inexact_alm_rpca\inexact_alm_rpca\PROPACK\dbdqr.f

     文件     256353  2014-08-10 08:52  Orl_SPP_SRC 代码\inexact_alm_rpca\inexact_alm_rpca\PROPACK\helio.mat

     文件       3544  2014-08-10 08:52  Orl_SPP_SRC 代码\inexact_alm_rpca\inexact_alm_rpca\PROPACK\lanbpro.doc

     文件      19514  2014-08-10 08:52  Orl_SPP_SRC 代码\inexact_alm_rpca\inexact_alm_rpca\PROPACK\lanbpro.m

     文件       3544  2014-08-10 08:52  Orl_SPP_SRC 代码\inexact_alm_rpca\inexact_alm_rpca\PROPACK\lanbpro.txt

     文件       2522  2014-08-10 08:52  Orl_SPP_SRC 代码\inexact_alm_rpca\inexact_alm_rpca\PROPACK\laneig.doc

     文件       9695  2014-08-10 08:52  Orl_SPP_SRC 代码\inexact_alm_rpca\inexact_alm_rpca\PROPACK\laneig.m

     文件       2522  2014-08-10 08:52  Orl_SPP_SRC 代码\inexact_alm_rpca\inexact_alm_rpca\PROPACK\laneig.txt

     文件       3336  2014-08-10 08:52  Orl_SPP_SRC 代码\inexact_alm_rpca\inexact_alm_rpca\PROPACK\lanpro.doc

     文件      14762  2014-08-10 08:52  Orl_SPP_SRC 代码\inexact_alm_rpca\inexact_alm_rpca\PROPACK\lanpro.m

     文件       3336  2014-08-10 08:52  Orl_SPP_SRC 代码\inexact_alm_rpca\inexact_alm_rpca\PROPACK\lanpro.txt

     文件       2386  2014-08-10 08:52  Orl_SPP_SRC 代码\inexact_alm_rpca\inexact_alm_rpca\PROPACK\lansvd.doc

............此处省略38个文件信息