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