10.07.2015 Views

Modèles de Markov triplets en restauration des signaux

Modèles de Markov triplets en restauration des signaux

Modèles de Markov triplets en restauration des signaux

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

CHAPITRE IV. MODÈLE DE MARKOV TRIPLETNous procédons dans la suite à une série d’expéri<strong>en</strong>ces pour étudier le modèle<strong>de</strong> la chaîne <strong>de</strong> <strong>Markov</strong> cachée M-stationnaire et le comparer avec le modèle <strong>de</strong> lachaîne <strong>de</strong> <strong>Markov</strong> cachée.Nous comm<strong>en</strong>çons par simuler une chaîne M-stationnaire à bruit gaussi<strong>en</strong> indép<strong>en</strong>dantà <strong>de</strong>ux classes, <strong>de</strong>ux stationnarités, et <strong>de</strong> taille N= 256×256 (K= 2,M= 2). La loi du couple(X,U) est donnée par (IV.8). Nous comm<strong>en</strong>çons par simulerle processus U <strong>de</strong> loi initiale p(u 1 )= 1 2 et <strong>de</strong> matrice <strong>de</strong> transitionM=(p i,j)avec p i,j =p(u n+1 =i∣u n =j) :M=(0.999 0.0010.001 0.999 ) (IV.11)Conditionnellem<strong>en</strong>t à U les matrices <strong>de</strong> transitions <strong>de</strong> X sont données parM 1 etM 2 qui ont respectivem<strong>en</strong>t comme coeffici<strong>en</strong>ts p(x n+ ∣u n+1 = u 1 ,x n ) et p(x n+ ∣u n+1 =u 2 ,x n ) :M 1 =(0.99 0.010.01 0.99 ) etM 2=(0.7 0.30.3 0.7 ) (IV.12)Dans un premier exemple, conditionnellem<strong>en</strong>t à X, les Y n sont <strong>de</strong>s gaussi<strong>en</strong>nes <strong>de</strong>moy<strong>en</strong>nes m 1 = 0 si x n = ω 1 , m 2 = 2 si x n = ω 2 et <strong>de</strong> même variance σ 2 = 1. Pourvisualiser les processus simulés, prés<strong>en</strong>tés à la figure (IV.2), sous forme d’imagesnous utilisons la transformation d’Hilbert-Peano (I.5).(a) (b) (c)Figure IV.2 – Premier exemple <strong>de</strong> simulation d’une chaîne M-stationnaire à bruitindép<strong>en</strong>dant. (a) : réalisation x 1∶N ; (b) : réalisation u 1∶N (c) : réalisation y 1∶N .Dans le second exemple, on fait augm<strong>en</strong>ter le bruit, au lieu <strong>de</strong> m 2 = 2 pour lesecond classe nous choisissons m 2 =1 et nous gardons les valeurs <strong>de</strong>s autres variablessans modifications. Les réalisations <strong>de</strong>s processus simulés sont représ<strong>en</strong>tées sur lafigure (IV.3).61

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!