5.4 Adaptation du recuit simulé au
problème
5.4.1 Initialisation
Cette étape consiste à déterminer une
solution réalisable du problème, qui constituera une solution de
départ pour l'algorithme du recuit simulé. l'initialisation est
obtenue par l'heuristique NSCM (Nombre de Stations de Compression en Marche)
décrite précédemment.
5.4.2 Paramètres opérationnelles
· O0 = 500.
· á = 0,5.
· nr = 50.
· Of = 0,0006.
5.4.3 Voisinage
Pour définir une solution
(X',Z(X')) a partir d'une solution
(X,Z(X)), en considérant uniquement la partie 3 (le
vecteur des vitesses), on génère aléatoirement une
position entre 1 et 5, si la station correspondante à la position
générée est en marche alors en génère
aléatoirement (entre 3250 et 6825) une vitesse de rotation des
compresseurs en marche dans cette station.
5.4.4 Principales étapes
(1). Choisir aléatoirement une position k entre
1 et 5 tel que la station k est en marche.
(2). Générer aléatoirement une vitesse de
rotation S[k] entre Smin et
Smax.
(3). Si solution réalisable trouvée alors Fin, si
non aller à (2).
5.4. ADAPTATION DU RECUIT SIMULÉ AU PROBLÈME
?
?????a1 + a2
H [k] :=
q[k] ! q[k] !2
q[k] !3?????? (S [k])2 ;
+ a3 + a4
S [k] S [k] S
[k]
Procédure Voisinage ( E : X : Vecteur, S :
X' : Vecteur)
Début
Pour k := random(1,5);
Si w[k] = 0 alors
Répéter
k := random(1,5);
Jusqu'à w [k] ~ 0;
Fsi;
S[k] := random(3250,6825);
? q[k] !
? q[k] !2 q[k]
H [k] :=
????!3?????? (S [k])2
;
a1 + a2 + a3 + a4
S [k] S [k] S
[k]
q[k] ! q[k] !2
q[k] !3
Eta[k] := b1 + b2 +
b3 + b4 ;
S [k] S [k] S
[k]
Si ((Eta[k] > 1)ou
(Eta[k] < 0)) et ((H [k]
> 45769.670)ou(H [k] <
9092.202)) alors Répéter
S[k] := random(3250,6825);
q[k] ! q[k] !2
q[k] !3
Eta[k] := b1 + b2 +
b3 + b4 ;
S [k] S [k] S
[k]
jusqu'à ((Eta[k] < 1)et
(Eta[k] > 0)) et ((H [k]
< 45769.670)et(H [k] >
9092.202))
Fsi; Fin.
94
5.4. ADAPTATION DU RECUIT SIMULE AU PROBLÈME
5.4.5 Organigramme d'adaptation de recuit simulé au
problème
|
|
|
|
5mlm Smax, a1, a2,
a3, a4 61, 62, b3, b4
0°, Of, a,Tt,
|
|
|
|
|
Solutions initiale X trouvée par l'heuristique, Z(X)
|
|
|
|
nb :=1
|
Choisir X' dans V(X), r :=random(1,5), si w[r] 0 alors S[r]
:=random(3250,6825)
{/ {( `I 2 1j 3 \j 2
H[r].=(al+l1
S[r])+c1(S[r]! +d'(S[r]l 1XS[l]
2 3
eta[r]:=a2+b2(s~r])+(s~r]I
+d2(S[r]~
J
I
AZ := Z (X') -- Z (X)
|
Non
|
Oui
|
|
|
p := randoin[0,1]
|
|
|
|
|
|
|
|
Oui
|
|
|
|
-BE
p < e r
|
Non
|
|
|
|
|
|
|
nb < 20
|
Non
Oui
nb := nb-F1
nb := nb+1
8 := a8
Oui
futioml X. Z(X) ,eta : vecteur
**
Non
4
95
Fin
FIGURE 5.2 -- Organigramme d'adaptation de récuit
simulé au problème
96
| 
