%% Parameter setup
% Simulation parameters
ResInt=2; % Resolution factor (integer) for dence discrete time-grid; recommened greater than 2
% System parameters
N=127; % Length of pseudo-random patterns
SeqType=1; % Pseudo-random patterns are (0: Samples of Bernoulli random vector) or (1: Irregular cyclic shifts of the maximum length sequence).
% Signal parameters
fmax=10e9; % Maximum frequency of x(t)
B=5e6; % Maximum bandwidth of narrow bands in x(t)
W=15; % Length of compressed samples
%% Simulation script
NTrials=300; % Number of Monte Carlo trials
SymmetricSupports=false; % if the spectral supports of x(t) are symmetric, OMP can find two symmetric supports per iteration.
KB=10; % Number of narrow bands in x(t). It must be an even number due to real signal x(t).
SNR_vars=[10^4]; % in dB
M_vars=[3]; %Number of channels
p_vars=1:4; % Aliasing parameter
Data1=cell(length(p_vars),length(M_vars)); % Memory for saving the supports recovery results
pindex=0;
for p=p_vars
if p==1, NonIdealFilter=false; else NonIdealFilter=true;end; % 0: ideal LPF, 1: random LPF
p
pindex=pindex+1;
mindex=0;
for M=M_vars
mindex=mindex+1;
M
q_vars=(p+1):(13*p+1); % Sampling rate=f_p*q/p;
q_vars=q_vars(mod(q_vars,2)==1);
q_vars=q_vars(gcd(q_vars,p)==1);
data1_vector=zeros(1, length(q_vars));
data2_vector=data1_vector;