r: Adrien MEYNARD
% Email: adrien.meynard@univ-amu.fr
% Created: 2017-12-19

%% Load signal
clear all; close all; clc;

warning off;
addpath('cwt');
addpath('deform_estimation');

load('signals/toy_sig_paper');
T = length(y);
%% Joint estimation

Dt = 200; % temporal subsampling for the deformation estimation
dgamma0 = ones(1,T); % gamma'(t) initialization
a0 = ones(1,T); % a(t) initialization

wav_typ = 'sharp'; % wavelet type (cf. cwt_JEFAS.m)
wav_paramWP = 25; % corresponding parameter for warping estima