压缩感知 莱斯大学的一些代码

莱斯大学关于压缩感知的一些代码，供大家参考

% CSSA1.m:  CONDENSING SORT AND SELECT ALGORITHM （1-d）
%
% B:     input data
% L:     number of levels in wavelet transform
% vol:   stopping volume of mask （default is full volume）
%        NOTE: does NOT include scaling coeff volume!
%
% mask:  final supernode configurations ordered 
%
% Assumes that the underlying tree is the 1-d dyadic wavelet tree
%
% RGB INI September 1998

function mask = cssa1（BLvol）;

%---------------------------------------------------------------------------%
% INITIALIZATIONS

% make input data a column vector
B    = B（:）;
N    = length（B）;
cvol = 0;
sno  = 2;
% scaling coefficient with highest index in wavelet array
parroot = N/2^L;

if nargin < 3
  vol = （N-N/2^L）;
end
if vol > （N-N/2^L）
  disp（‘ERROR: volume parameter cannot be > number of wavelet coeffs‘）
  mask = [];
  return
end

% set up each node as a supernode.  use an Nx5 matrix:
%    col 1: utp  （uptree pointer）
%    col 2: rp   （supernode root pointer）
%    col 3: snv  （supernode value）
%    col 4: num  （number of internal nodes）
%    col 5: mask （mask/kernel value）

sn（1:N 1） = （1:N）‘;      % point uptree top pointers to self
sn（1:N 2） = -ones（N1）;  % each node is a sn root => -1 flag
sn（1:N 3） = B;           % each snv is the data itself
sn（1:N 4） = ones（N1）;   % num = 1 for each sn   
sn（1:N 5） = zeros（N1）;  % output mask is zero at start
sn（1:N/2^L 5） = ones（N/2^L1）;  % except for parents of roots （“zero“ nodes）

% set up SNV queue  （include scaling coeff values just so that Q never empties
snvQ（1:N 1） = （1:N）‘;    % each node is its own sn topnode
snvQ（1:N 2） = B;         % each snv is the wavelet data itself
snvQ（1:N/2^L2） = NaN*snvQ（1:N/2^L2）;


%---------------------------------------------------------------------------%
% MAIN LOOP （follow pseudo code on p. 138 and steps on p. 142）

% notation:   MAX refers to the current sn we take from snv queue
%             CHECK refers to the parent sn it abutts

while （cvol < vol）

  % STEP 1
  % find supernode with largest SNV （sn MAX）
  % delete it from the snvQ as well
  topmax = Qfind（snvQ）;
  snvQ = Qdel（snvQtopmax）;
  
  % STEP 2
  % find the root node of sn MAX and then its parent
  % this latter node is a leaf of the CHECK sn above ours
  if sn（topmax2） == -1
    % then topmax is the sn root
    leafcheck = ceil（topmax/2）;
  else
    % go to the sn root
    leafcheck = ceil（sn（topmax2）/2）;
  end
  
  
  % check whether CHECK sn has mask=10 （first find top node of CHECK）
%  topcheck = Tfind（snleafcheck）;
  [topchecksn] = Tfind（snleafcheck）;       % need to pass sn back
                                             % if we do path comp step
    
  if sn（topcheck5） > 0 
    % we are in the money - CHECK sn has mask > 0.  so take sn MAX
    % STEP 3 （usual case - not the first sn taken）: 
    sn = umask（sntopmaxsnoN）;
    cvol = cvol + sn（topmax4）; 
    sno = sno + 1;
 
  else
    %

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

     文件       6148  2009-08-05 19:08  RICE UNIVERSITY代码\modelcs_v1.1\1D trees\.DS_Store

     文件       3184  2008-07-15 12:35  RICE UNIVERSITY代码\modelcs_v1.1\1D trees\cssa1.m

     文件       1960  2008-07-15 12:35  RICE UNIVERSITY代码\modelcs_v1.1\1D trees\greedy1.m

     文件       2090  2009-08-05 10:40  RICE UNIVERSITY代码\modelcs_v1.1\1D trees\treemp.m

     文件       1747  2009-08-04 16:48  RICE UNIVERSITY代码\modelcs_v1.1\1D trees\treemp_example.m

     文件       2022  2009-08-04 14:58  RICE UNIVERSITY代码\modelcs_v1.1\1D trees\treemp_fun.m

     文件       6148  2009-08-05 19:08  RICE UNIVERSITY代码\modelcs_v1.1\2D trees\.DS_Store

     文件       1565  2008-08-13 19:10  RICE UNIVERSITY代码\modelcs_v1.1\2D trees\cosamp_2D_nf.m

     文件       3938  2008-07-30 15:14  RICE UNIVERSITY代码\modelcs_v1.1\2D trees\cssa2.m

     文件       2451  2008-07-15 12:35  RICE UNIVERSITY代码\modelcs_v1.1\2D trees\greedy2.m

     文件     538749  2008-08-13 17:10  RICE UNIVERSITY代码\modelcs_v1.1\2D trees\peppers.png

     文件       1974  2009-08-05 10:34  RICE UNIVERSITY代码\modelcs_v1.1\2D trees\treemp_2D_example.m

     文件       2439  2009-08-05 10:42  RICE UNIVERSITY代码\modelcs_v1.1\2D trees\treemp_greedy_2D.m

     文件       2611  2009-08-05 10:26  RICE UNIVERSITY代码\modelcs_v1.1\2D trees\treemp_nf_2D.m

     文件       6148  2009-08-05 19:08  RICE UNIVERSITY代码\modelcs_v1.1\Block sparsity\.DS_Store

     文件       2346  2009-08-05 11:06  RICE UNIVERSITY代码\modelcs_v1.1\Block sparsity\jsmp.m

     文件       1301  2009-08-05 11:04  RICE UNIVERSITY代码\modelcs_v1.1\Block sparsity\jsmp_example.m

     文件       2454  2009-08-05 11:06  RICE UNIVERSITY代码\modelcs_v1.1\Block sparsity\jsmp_fun.m

     文件       6148  2009-08-05 19:09  RICE UNIVERSITY代码\modelcs_v1.1\Clustered sparsity\.DS_Store

     文件       1873  2009-08-05 13:40  RICE UNIVERSITY代码\modelcs_v1.1\Clustered sparsity\kc.m

     文件        632  2009-08-05 13:37  RICE UNIVERSITY代码\modelcs_v1.1\Clustered sparsity\kc_example.m

     文件       6148  2009-08-05 19:09  RICE UNIVERSITY代码\modelcs_v1.1\Delta\.DS_Store

     文件        659  2009-08-05 14:28  RICE UNIVERSITY代码\modelcs_v1.1\Delta\bestdelta.m

     文件       2468  2009-08-05 14:44  RICE UNIVERSITY代码\modelcs_v1.1\Delta\deltarec.m

     文件        760  2009-08-05 15:00  RICE UNIVERSITY代码\modelcs_v1.1\Delta\delta_example.m

     文件       2678  2012-02-15 09:19  RICE UNIVERSITY代码\modelcs_v1.1\README.txt

     文件       6148  2009-08-05 19:09  RICE UNIVERSITY代码\modelcs_v1.1\Sparsity\.DS_Store

     文件          0  2010-10-01 16:58  RICE UNIVERSITY代码\modelcs_v1.1\Sparsity\;.m

     文件       1743  2009-07-31 15:14  RICE UNIVERSITY代码\modelcs_v1.1\Sparsity\cosamp.m

     文件       1127  2009-08-04 11:20  RICE UNIVERSITY代码\modelcs_v1.1\Sparsity\cosamp_example.m

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