%%%%%%  Film pour deux targets avec le programme Mdeuxtarget.c %%%%%%%%
% Attention il faut remettre a jour les parametres taille et n 
% pour que ca marche comme il faut
taille=300;
n=100;
no=2*n;
data=zeros(2*n , taille);
load Mdeuxt300n150p85.txt;
data=Mdeuxt300n150p85( : , : );
observation=data(1:2:no,1:taille)';
exacte=data(2:2:no,1:2);
estimee=data(2:2:no,3:4);
s=0; u=0; v=0;
clear A B XOT YOT XMT YMT XET YET
A(1)=1;  A(2)=n; B(1)=taille+10; B(2)=taille+10;
E(1)=1;  E(2)=n; F(1)=2*taille+30; F(2)=2*taille+30;
ee(1)=1;  ee(2)=n; ff(1)=-2; ff(2)=-2;
c(1)=1;  c(2)=n; d(1)=3*taille+43; d(2)=3*taille+43;
movie_figure1=figure('position',[30 80  900 700]);
M=moviein(n);
for t=1:(n) 
% Vecteur des observations pour les ploter
   for i=1:taille;
     a=observation(i,t);
       if a>0
         s=s+1;
        YOT(s)=i+2*taille+40;
         XOT(s)=t;
     end
    end
% Vecteur de la trajectoire et de son estimée
  YMT(u+1)=exacte(t,1)+taille+20;
  YMT(u+2)=exacte(t,2)+taille+20;
  XMT(u+1)=t;
  XMT(u+2)=t;
  YET(u+1)=estimee(t,1);
  YET(u+2)=estimee(t,2);
  XET(u+1)=t;
  XET(u+2)=t;
  u=u+2;
plot(XOT,YOT,'.b',XET,YET,'.g',XMT,YMT,'+r')
axis off
hold on
plot(A,B,'-k','linewidth',7)
hold on
plot(E,F,'-k','linewidth',7)
hold on
plot(c,d,'-k','linewidth',3)
hold on
plot(ee,ff,'-k','linewidth',3)
%set(gca,'Title',text('String','EN BLEU LES OBSERVATIONS, EN ROUGE LES CIBLES ET EN VERT LES ESTIMATIONS','Color','k'))
axis([ 1 n -3 3*taille+45])
figure(movie_figure1);
M(:,t)=getframe(movie_figure1);
end