Pastebin

bouff
def RecuitSimule(etat_initial,T_0, taux_refroidissement, iter_max):
    n = etat_initial.copy()
    f_n = evaluer(n)
    n_best = n.copy()
    f_n_best = f_n
    t= T_0
    N = len(n)
    for i in range(iter_max):
        n_voisin = n.copy()
        index= random.randint(0,N-1)
        n_voisin[index] = random.randint(0,N-1) 
        f_voisin = evaluer(n_voisin)
        if f_voisin < f_n:
            n = n_voisin
        else : 
            delta = abs(f_voisin - f_n)
            proba = 1/(np.exp(delta/t))
            r = random.random()
            if(r < proba):
                n = n_voisin
                f_n = f_voisin
        if f_n < f_n_best:
            n_best = n.copy()
            f_n_best = f_n
        t *= taux_refroidissement

    return n_best,f_n_best
Créé il y a 1 mois.