timations</span>

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

stopWP = 2e-2; <span class="comment">% minimal gap between two steps in the gradient</span>
itWP = 6; <span class="comment">% number of gradient iterations</span>

r = 1e-3; <span class="comment">% regularization parameter</span>

Nf = 2500; <span class="comment">% number of frequen