稀疏表示SRC算法

经典的稀疏表示SRC算法，适合广大人脸识别的同学使用参考-The classic sparse representation SRC algorithm

function [xx_debiasobjectivetimesdebias_startmses]= ...
    GPSR_Basic（yAtauvarargin）
%  
% GPSR_Basic version 5.0 December 4 2007
%
% This function solves the convex problem 
% arg min_x = 0.5*|| y - A x ||_2^2 + tau || x ||_1
% using the algorithm GPSR-Basic described in the following paper
%
% “Gradient Projection for Sparse Reconstruction: Application
% to Compressed Sensing and Other Inverse Problems“
% by Mario A. T. Figueiredo Robert D. Nowak Stephen J. Wright
% Journal of Selected Topics on Signal Processing December 2007
% （to appear）.
%
% -----------------------------------------------------------------------
% Copyright （2007）: Mario Figueiredo Robert Nowak Stephen Wright
% 
% GPSR is distributed under the terms
% of the GNU General Public License 2.0.
% 
% Permission to use copy modify and distribute this software for
% any purpose without fee is hereby granted provided that this entire
% notice is included in all copies of any software which is or includes
% a copy or modification of this software and in all copies of the
% supporting documentation for such software.
% This software is being provided “as is“ without any express or
% implied warranty.  In particular the authors do not make any
% representation or warranty of any kind concerning the merchantability
% of this software or its fitness for any particular purpose.“
% ----------------------------------------------------------------------
% 
% Please check for the latest version of the code and paper at
% www.lx.it.pt/~mtf/GPSR
%
%  ===== Required inputs =============
%
%  y: 1D vector or 2D array （image） of observations
%     
%  A: if y and x are both 1D vectors A can be a 
%     k*n （where k is the size of y and n the size of x）
%     matrix or a handle to a function that computes
%     products of the form A*v for some vector v.
%     In any other case （if y and/or x are 2D arrays） 
%     A has to be passed as a handle to a function which computes 
%     products of the form A*x; another handle to a function 
%     AT which computes products of the form A‘*x is also required 
%     in this case. The size of x is determined as the size
%     of the result of applying AT.
%
%  tau: usually a non-negative real parameter of the objective 
%       function （see above）. It can also be an array the same 
%       size as x with non-negative entries; in this  case
%       the objective function weights differently each element 
%       of x that is it becomes
%       0.5*|| y - A x ||_2^2 + tau^T * abs（x）
%
%  ===== Optional inputs =============
%
%  
%  ‘AT‘    = function handle for the function that implements
%            the multiplication by the conjugate of A when A
%            is a function handle. If A is an array AT is ignored.
%
%  ‘StopCriterion‘ = type of stopping criterion to use
%                    0 = algorithm stops when the relative 
%                        change in the number of non-zero 
%                        components o

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

     文件      21922  2008-01-15 16:22  SRC\GPSR_Basic.m

     文件       4555  2010-07-14 14:01  SRC\SRC.m

     目录          0  2014-09-27 11:39  SRC

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

                26477                    3