aa=load('../Données/doris.data');
t=2004+(31+29+23)/3666+(aa(:,1)-aa(1,1))/365.25;
E=aa(:,8);
N=aa(:,9);
U=aa(:,10);
t0=2004.8;
t1=t(t<t0);
n1=length(t1);
t2=t(t>=t0);
E1=E(t<t0);
E2=E(t>=t0);
N1=N(t<t0);
N2=N(t>=t0);
U1=U(t<t0);
U2=U(t>=t0);
t2=t2(1:n1);
E2=E2(1:n1);
N2=N2(1:n1);
U2=U2(1:n1);
[miE1,maE1]=confintmean(E1,0.05);
[miN1,maN1]=confintmean(N1,0.05);
[miU1,maU1]=confintmean(U1,0.05);
[miE2,maE2]=confintmean(E2,0.05);
[miN2,maN2]=confintmean(N2,0.05);
[miU2,maU2]=confintmean(U2,0.05);
subplot(311)
plot(t1,E1,'kx',t2,E2,'ro',t1,t1*0+maE1,'b',t1,t1*0+miE1,'b'...
    ,t2,t2*0+maE2,'b',t2,t2*0+miE2,'b','linewidth',2);
title('Position de la station de Jiufeng (Chine)');
ylabel('Composante Est (mm)')
subplot(312)
plot(t1,N1,'kx',t2,N2,'ro',t1,t1*0+maN1,'b',t1,t1*0+miN1,'b'...
    ,t2,t2*0+maN2,'b',t2,t2*0+miN2,'b','linewidth',2);
ylabel('Composante Nord (mm)')
subplot(313)
plot(t1,U1,'kx',t2,U2,'ro',t1,t1*0+maU1,'b',t1,t1*0+miU1,'b'...
    ,t2,t2*0+maU2,'b',t2,t2*0+miU2,'b','linewidth',2);
ylabel('Composante Up (mm)')