Génération d’une marche aléatoire sur un graphe d’état : randmarkov

8.4.57 Génération d’une marche aléatoire sur un graphe d’état : randmarkov

randmarkov(M,i0,n) génère une suite de n états (chaine de Markov) partant de i0 dont les probabilites de transitions sont données par M (matrice stochastique) ou
randmarkov(v,i0) génère une matrice stochastique ayant p boucles récurrentes v=[n1,...,np] et i0 états transients

randmarkov(M,i0,n) renvoie la liste des états qui a n+1 éléments.
On tape :

randmarkov([[0,0,1/2,0,1/2],[0,0,1,0,0],[1/4,1/4,0,1/4,1/4], [0,0,1/2,0,1/2],[0,0,0,0,1]],2,20)

On obtient par exemple :

[2,3,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4]

On tape :

randmarkov([1,2,1,3],4)

On obtient :

une matrice carrée de dimension [11,11]

en effet 1+2+1+3+4=11