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