////////////////////////////////////////
// l : longueur de la chaine
// e : proba de passer de 1 à 2 pour la chaine cachée
function esti(l,e)
  // Paramètres
r=[1-e e;2*e 1-2*e]
pi1=[.3 .3 .3 .1;.3 .3 .1 .3;.3 .1 .3 .3;.1 .3 .3 .3]
pi2=[.5 .3 .1 .1;.4 .4 .1 .1;.4 .1 .4 .1;.5 .3 .1 .1]
// Trajectoire cachée
   // A faire
// Trajectoire observée
   // A faire
// Initialisations
P=zeros(2,l);P(1,1)=1;
F=zeros(2,l);
L=zeros(2,l);
// Prédiction-Filtrage
   // A faire
// Lissage
   // A faire
// Affichage
xbasc();
plot2d([1,l],[.5,.5],1)
plot2d([1:l],L(1,:),2)
plot2d([1:l],L(2,:),5)
plot2d2([1:l],U(2:l+1),3)
legends(['Trajectoire','Proba état 1','Proba état 2'],[3,2,5],1)
// Evaluation 
  // A faire
endfunction