Sans vouloir polémiquer, ca se ressemble tout de même un peu. 
 1: Initialize ? (ai = si, bi = 0.5 × si, d =Random number in (0, 1) )
 
 2: while does not converge do
 3:   for j = 1 : N do
 4:     compute .. based on Equation (11)
 5:   end for
 6: 
 7:   for i = 1 : M do
 8:     compute ... based on Equation (17)
 9:     update ...
10:   end for
11: 
12: end while
13-25: (utiliser le modèle trouvé pour faire des trucs) 
// initial guess for the laws parameters (uniform)
(...)
// performe EM algorithm
for(int loop=0;loop<5000;loop++) {
    for(int i=0;i<N;i++) {
        // ---- E Step ----
    }
    for(int k=0;k<G;k++) {
        // ---- M Step ----
    }
    for(int k=0;k<G;k++) {
        // ---- Update ----
    }
    if( convergence < 1E-10 ) break;
}
// display mixture coefficients and law parameters
(...)
					
Partager