 !c:
 lat=dynre2he(:,5);long=dynre2he(:,6);
 vx=dynre2ve;vy=dynre2vn;
 he=dynre2he;
 temps=he(:,2)*60+he(:,3);
 plot(temps)
 %test des donnees lat et long

 plot(long,lat)
 %ok sans erreur

 vnref=mean(vy(:,25:30)');
 veref=mean(vx(:,25:30)');
 [n,z]=min(temps)

size(temps)


 temps(690:753)=temps(690:753)+1440*(ones(64,1));
 subplot(211);plot(temps,veref)
 xlabel('temps en minutes du premier jour')
 title('veref')
 subplot(212);plot(temps,vnref)
 xlabel('temps en minutes du premier jour')
 title('vnref')
 print
 for j=1:24
verel(:,j)=vx(:,j)-veref';
vnrel(:,j)=vy(:,j)-vnref';
end
 plot(temps',verel(:,1))
 plot(temps',vnrel(:,1))
 hold on
 [b,a]=butter(5,1/10);
 for j=1:24
verelf(:,j)=filtfilt(b,a,verel(:,j));
vnrelf(:,j)=filtfilt(b,a,vnrel(:,j));
end
 !cd \matlab\bin
 bintra=[24 20 15 10 7 4 3 2 1];

 

 subplot(211);veccant(verelf,vnrelf,300,temps,bintra,5, 748,5)
 title('Dynaproc reseau 2   Vrel /vref(bin25-30)')
 subplot(212);plovec(43.3,43.6,7.65,10,lat,long,verelf(:,1),vnrelf(:,1),2)
 title('vitesse bin1/bin25-30')
 xlabel('longitude')
 ylabel('latitude')
 print
 return
