,th,AMsup,'bREPLACE_WITH_DASH_DASH','linewidth',2); grid on;
% title('Amplitude modulation estimation'); 
legend({'JEFAS estimate','$\theta_{n,1}$','$\theta_{n,1} \pm 3\sqrt{\mathrm{CRLB}(\theta_{n,1})}$'},'Interpreter','latex'); axis([0 t(end) 0 3]);
xlabel('Time (s)'); ylabel('Amplitude modulation');

figure;
plot(t,thetaWP_B,'g',t,log2(dgammaML),'r',t,log2(dgamma),'b-.',th,WPinf,'bREPLACE_WITH_DASH_DASH',th,WPsup,'bREPLACE_WITH_DASH_DASH','linewidth',2); grid on; axis tight;
% title('Warping estimation')