main_1d_other_run.m

% This code takes a 1D signal f of specified type, applies a sampling 
% or approximation strategy, and computes the energy spectrum. The goal
% is to estimate the energy spectrum well on a logarithmic scale.
% Depending on DECODER, the following methods are used. 
% (32)      -a fixed rate sampling and fourier interpolation
% (31,312)  -best M/2 or M-term wavelet approximation by thresholding
%           -compressive sampling/compressed sensing with linear sampling y=A*f, 
%            and reconstruction by
% (1)          -Lumped Orthogonal Maching Pursuit (LOMP), or
% (82)         -Spectrum Estimation by Sparse Optimization (SpESO) using
%               Quasi-Oracle Multilevel OMP (QOMOMP)
% The results can be saved into .mat files, and plotted with the plot scripts.
% 
% Copyright (C) 2014  Gudmundur Adalsteinsson
% See file LICENCE for licence and warranty details
% See file README.md for more information
%

timerID=tic;

SAVEDLOC='./';
path(path,genpath('functions'))

%LOADVAR=0;	% load saved variables?
SAVEVAR=1;	% save variables?

DWTTYPE=1;	% choose wavelet package
			% 1:Wavelab, 2:MATLAB's Wavelet toolbox, 0:cdf9/7

JH=0;
RHHDA=0;
SYNTH=0;

% =============data type===============
% JH data
%JH=1;  mstr='jh';
%N=1024;

% RHHDA hot-wire data
% RHHDA=1;  mstr='rhh';
% N = 1024*32;

% SYNTH data
SYNTH=1;  mstr='synth';
N = 1024*32;
synthtype=1;	 % 1:FFT,     2:DWT
syntslope=1;	 % 1:(3,5/3), 2:(5/3,3)
% =====================================

% undersampling ratio
do1vec=8;		% for SpESO (DECODER=82) use half:  82:(32,16,8), 1:(16,8,4)
SIMN=length(do1vec);
SIMN0=4;		% number of experiments for each case

EXTRA='';		% extra text for special runs, '' for nothing
K=284;			% filter length

Mvec=ceil((N+K-3)./do1vec);   % any M opt db
Mmin=Mvec(end);

% [1,31,312,32] , 82
for MATRIX=1      % 1=filter,  not used: 2=conv., 3=full
	for DECODER=[1,31,312,32]  %  1=LOMP, 31=M/2Best, 312=MBest, 32=shannon, 82=SpESO/QOMOMP
		
		FIGQUAL=sprintf('_%d_%d_%d_%s_%d_%d_%d_%s',N,Mmin,SIMN0,mstr,K,MATRIX,DECODER,EXTRA);
		SAVEFILE = [SAVEDLOC,'saved_1d',FIGQUAL,'.mat']; % for variables
		
		disp(' ')
		disp(FIGQUAL)
		
		spectra_error_1d_calc;  % Calculation
		
	end
end


toc(timerID)

disp('done')




% Copyright (C) 2014  Gudmundur Adalsteinsson
% See file LICENCE for licence and warranty details