tion
wav_param = 500; % corresponding parameter for spectrum and AM estimations

NbScales = 106;
scalesAM = 2.^(linspace(-1,3.5,NbScales));
subrate = 7; % subsampling step for the scales to ensure the covariance invertibility
scalesWP = scalesAM(1:subrate:end);

stopWP = 2e-2; % minimal gap between two steps in the gradient
itWP = 6; % number of gradient iterations

r = 1e-3; % regularization parameter

Nf = 2500; % number of frequencies for spectrum estimation
NbScalesS = 212;
scalesS = 2.^(linspace(-1,3.5