age = [47 58 82 97 133 134];
moy = [mean(j_47) mean(j_58) mean(m_82) mean(m_87) mean(v_133) mean(v_134)];
va = [var(j_47) var(j_58) var(m_82) var(m_87) var(v_133) var(v_134)];
va_moy = va./moy;

fm = zeros(6,1);
rap1 = fm;
rap2 = fm;
vad = fm;
mamp = fm;
x = j_47;
n=1;
[t,y] = reech_stride(x,4);
[es,ei,ampx,xmd] = enveloppes(y,45);
vad(n) = var(xmd);
mamp(n) = mean(ampx);
[fmoy,rapport,rapport2] = estim_DSP_3(xmd,4,1);
pause
fm(n) = fmoy;
rap1(n) = rapport;
rap2(n) = rapport2;
x = j_58;
n=2;
[t,y] = reech_stride(x,4);
[es,ei,ampx,xmd] = enveloppes(y,45);
vad(n) = var(xmd);
mamp(n) = mean(ampx);
[fmoy,rapport,rapport2] = estim_DSP_3(xmd,4,1);
 pause
fm(n) = fmoy;
rap1(n) = rapport;
rap2(n) = rapport2;
x = m_82;
n=3;
[t,y] = reech_stride(x,4);
[es,ei,ampx,xmd] = enveloppes(y,45); 
vad(n) = var(xmd);
mamp(n) = mean(ampx);
[fmoy,rapport,rapport2] = estim_DSP_3(xmd,4,1);
pause
fm(n) = fmoy;
rap1(n) = rapport;
rap2(n) = rapport2;
x = m_87;
n=4;
[t,y] = reech_stride(x,4);
[es,ei,ampx,xmd] = enveloppes(y,45);
vad(n) = var(xmd);
mamp(n) = mean(ampx);
[fmoy,rapport,rapport2] = estim_DSP_3(xmd,4,1);
 pause
fm(n) = fmoy;
rap1(n) = rapport;
rap2(n) = rapport2;
x = v_133;
n=5;
[t,y] = reech_stride(x,4);
[es,ei,ampx,xmd] = enveloppes(y,45);
vad(n) = var(xmd);
mamp(n) = mean(ampx);
[fmoy,rapport,rapport2] = estim_DSP_3(xmd,4,1);
 pause
fm(n) = fmoy;
rap1(n) = rapport;
rap2(n) = rapport2;
x = v_134;
n=6;
[t,y] = reech_stride(x,4);
[es,ei,ampx,xmd] = enveloppes(y,45);
vad(n) = var(xmd);
mamp(n) = mean(ampx);
[fmoy,rapport,rapport2] = estim_DSP_3(xmd,4,1);
 pause
fm(n) = fmoy;
rap1(n) = rapport;
rap2(n) = rapport2;
fmoy = fm;
