GPSR 梯度投影法

压缩感知中的梯度投影法，可以更好地进行图像重构与压缩。

close all
% n is the original signal length
n = 2^12;

% k is number of observations to make
k = 2^10;

% number of spikes to put down
% n_spikes = floor（.01*n）;
n_spikes = 160;


% random +/- 1 signal
f = zeros（n1）;
q = randperm（n）;
f（q（1:n_spikes）） = sign（randn（n_spikes1））;
%f（q（1:n_spikes）） = randn（n_spikes1）;

% measurement matrix
disp（‘Creating measurement matrix...‘）;
R = randn（kn）;

% orthonormalize rows
R = orth（R‘）‘;

if n == 8192  
   % in this case we load a precomputed
   % matrix to save some time
   load Rmatrix_2048_8192.mat
end
%
disp（‘Finished creating matrix‘）;

hR = @（x） R*x;
hRt = @（x） R‘*x;

% noisy observations
sigma = 0; 
y = hR（f） + sigma*randn（k1）;

% regularization parameter
tau = 0.005*max（abs（R‘*y））;

first_tau_factor = 0.8*（max（abs（R‘*y））/tau）;
steps = 5;

debias = 0;

stopCri=3;
tolA=1.e-5;

disp（‘Starting GPSR BB monotonic‘）
[x_BB_monox_debias_BB_monoobj_BB_mono...
    times_BB_monodebias_start_BB_monomse_BB_mono]= ...
         GPSR_BB（yRtau...
         ‘Debias‘debias...
         ‘Monotone‘1...
         ‘Initialization‘0...
         ‘MaxiterA‘10000...
         ‘True_x‘f...
         ‘StopCriterion‘stopCri...
       	 ‘ToleranceA‘tolA...
         ‘Verbose‘0）;
t_BB_mono = times_BB_mono（end）;


disp（‘Starting GPSR BB monotonic with continuation‘）
[x_BB_mono_contx_debias_BB_mono_contobj_BB_mono_cont...
    times_BB_mono_contdebias_start_BB_monomse_BB_mono_cont]= ...
         GPSR_BB（yRtau...
         ‘Debias‘debias...
         ‘Continuation‘1...
         ‘ContinuationSteps‘steps...
         ‘FirstTauFactor‘first_tau_factor...
         ‘Monotone‘1...
         ‘Initialization‘0...
         ‘True_x‘f...
         ‘StopCriterion‘stopCri...
       	 ‘ToleranceA‘tolA...
         ‘Verbose‘0）;
t_BB_mono_cont = times_BB_mono_cont（end）;

disp（‘Starting GPSR Basic‘）
[x_Basicx_debias_Basicobj_Basic...
    times_Basicdebias_start_Basicmse_Basic]= ...
         GPSR_Basic（yRtau...
         ‘Debias‘debias...
         ‘Initialization‘0...
         ‘MaxiterA‘10000...
         ‘True_x‘f...
         ‘StopCriterion‘stopCri...
       	 ‘ToleranceA‘tolA...
         ‘Verbose‘0）;
t_Basic = times_Basic（end）;

disp（‘Starting GPSR Basic with continuation‘）
[x_Basic_contx_debias_Basic_contobj_Basic_cont...
    times_Basic_contdebias_start_Basicmse_Basic_cont]= ...
         GPSR_Basic（yRtau...
         ‘Debias‘debias...
         ‘Continuation‘1...
         ‘ContinuationSteps‘steps...
         ‘FirstTauFactor‘first_tau_factor...
         ‘Initialization‘0...
         ‘True_x‘f...
         ‘StopCriterion‘stopCri...
       	 ‘ToleranceA‘tolA...
         ‘Verbose‘0）;
t_Basic_cont = times_Basic_cont（end）;



fprintf（1‘\n\n-------------------------------------------------\n‘）   
fprintf（1‘-------------------------------------------------\n‘）   
fprintf（1‘Problem: n = %g  k = %g number of spikes = %g\n‘nkn_spikes）
fprintf（1‘Parameters: sigma = %g 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        5276  2007-12-04 11:46  figure_5.m
     文件        8489  2007-10-19 16:40  figures_1_2_3.m
     文件        7046  2007-09-23 20:49  demo_continuation.m
     文件        5225  2007-09-12 18:21  figure_4.m
     文件       12630  2007-09-06 15:10  IST.m
     文件        5849  2007-07-13 16:40  demo_image_deblur.m
     文件        2377  2007-07-12 18:41  figure_6.m
     文件        8569  2007-06-16 15:41  l1_ls.m
     文件       65126  2007-01-25 13:54  Camera.tif
     文件       23882  2009-01-16 13:25  GPSR_BB.m
     文件     1296396  2008-04-30 17:42  immi1.mat
     文件       21922  2008-01-15 16:22  GPSR_Basic.m
     文件        3178  2007-12-05 19:06  README.txt
     文件        6148  2007-12-04 11:46  .DS_Store