clear all
aa=load('..\Données\GLB.TS+dSST.txt');
for i=1:12
    T(i:12:126*12)=aa(:,i+1);
    t(i:12:126*12)=aa(:,1)+i/12;
end
for k=1900:2000
   tt=zeros(size(t));
   tt(t>=k)=t(t>=k)-k;
   plot(t,T)
   A=[0*t+1; t; tt; sin(2*pi*t); cos(2*pi*t)];
   y=A'\T';
   modele=A'*y; 
   z(k-1899)=var(T-modele');
end
plot(1900:2000,z)