Mathematics">
Support - RO 2024
Support - RO 2024
Support - RO 2024
FSJES- Agdal
Département des sciences de gestion
S5: Finance comptabilité & AGE
-*-
Recherche opérationnelle
Partie I
1- Concepts fondamentaux
2- Formalisation d’un problème de PL
3- Résolution graphique
4- Résolution par l’algorithme de simplexe
5- Dualité
6- Applications économiques
Partie II- La théorie des graphes
applications économiques.
• Modèle linéaire :
fonction linéaire de plusieurs variables à optimiser (premier
degré)
variables soumises à des contraintes :
» linéaires
» restriction de non négativité
Formalisation Algorithme
Modèle
Problème mathématique Solution
réel
Programmation linéaire
Mise en forme d’un programme linéaire
A 5 3 3 2 1800
B 2 5 4 6 8500
C 1 2 3 4 3500
•Profit par unité 7,30 8,64 8,00 9,00
Modélisation mathématique
But:
Compréhension du problème
Définir les 4 variables de décision :
x1 : la production hebdomadaire du produit I
x2 : la production hebdomadaire du produit II
x3 : la production hebdomadaire du produit III
x4 : la production hebdomadaire du produit IV
Le profit associé à une production (x1, x2, x3, x4) est :
z 7,30 x1 8,64 x2 8,00 x3 9,00 x4
x1 , x2 , x3 , x4 0
Les productions ne sont pas négatives :
Formalisation du problème
Fonction Objectif :
max z 7,30 x1 8,64 x2 8 x3 9 x4
Selon la nature des produits, il faudrait peut être imposer à x1, x2, x3,
x4 de ne prendre que des valeurs entières
Application 1 : Problème d’agriculture
Un agriculteur veut allouer 150 hectares de surface
irrigable entre culture de tomates et celles de piments.
Il dispose de 480 heures de main d’œuvre et de 440 m3
d’eau. Un hectare de tomates demande 1 heure de
Solution optimale
• ax by c 0 (a, b) (0,0)
• u (b, a) est un vecteur directeur de (D).
• y mx p m0
•avec, m a b et p c b
•m est le coefficient directeur de (D)
• * Soit la droite (D) : 2x 3 y 6 0
• u (3,2) est un vecteur directeur de (D).
•* Pour construire cette droite, il suffit de connaître deux points :
• « Par deux points on ne fait passer qu’une droite et une seule »
•* Déterminons, x 0 3
- la droite x 1 y 2 0
• - la droite y 2
• - la droite (D1) passant par A(1,2) et de vecteur directeur v (3,1)
y y x 1
(D)
(D1)
A
u 2 2
1 v
-3 0 3 x 0 1 3 x
y 2
• Soit la droite (D) d’équation :
ax by c 0 a, b 0,0
•Alors, les demi-plans de frontière (D) sont :
•L’ensemble des points de coordonnées (x,y) telles que :
ax by c 0
•L’ensemble des points de coordonnées (x,y) telles que :
ax by c 0
• Exemple : Résolvez graphiquement 2 x 3 y 6 0
y
• Remarquons que pour l’origine :
2(0) 3(0) 6 0
• 2 Donc le demi plan ne contient
pas l’origine
0 3 x
Région admissible (ou faisable
ou possible)
Région limitée par l’ensemble des
équations de contraintes du
Contrainte de non négativité
Point intérieur
Point frontière
Point extrême
x1
2 4 6 8 10 12 14 16 18
Région admissible bornée
x2
x1
Région admissible non bornée
x2
x1
Fonction « Objectif »
Déplacement de la fonction
x2
objectif à l’intérieur de la
région admissible pour
atteindre un extremum
A
x1 x1 x1
Exemple 1
Région admissible bornée :
une solution optimale unique
x2 maximiser z x1 2 x2
A (2,6) est l’unique sujet à
solution optimale
8
2 x1 x2 4
Zmax=14 x
A x2 8
6
1
x1 x2 4
x 5
4 1
x1, x2 0
x1
2 5 8
Exemple 2
Région admissible bornée :
une infinité de solutions optimales
x2 maximiser z 2 x1 2 x2
sujet à
8
2 x1 x2 4
Une infinité de x x2 8
solutions optimales 1
x1 x2 4
x 5
4 1
x1, x2 0
x1
2 5 8
Exemple 3
Région admissible non bornée :
pas de solutions optimales finies
maximiser z x1 x2
sujet à
x2 x1=8
1
2 x1 x2 8
x 8 x2 40
1
x1 8
8
x1 , x2 0
6
x1
2 4 6 8 10 12 14 16 18
Exemple 4
Région admissible non bornée :
une solution optimale unique
minimiser w x1 x2
sujet à
x2 x1=8
B(8,6) est l’unique 1
2 x1 x2 8
solution optimale x 8x 40
1 2
wmin=2 x1 8
8
x1, x2 0
B
6
-x1+8x2=40
4
x1
2 4 6 8 10 12 14 16 18
Exemple 5
Région admissible non bornée :
une infinité de solutions optimales
x2 mi z x1 1 / 2 x2
12 sujet à
x1 x2 8
10 x 8 x 40
1 2
8 2 x1 x2 12
x1, x2 0
6
Une infinité de
-x1+8x2=40
4 solutions optimales
x1
2 4 6 8 10 12 14 16 18
Cas pathologiques
1. Cas de redondance :
une contrainte est redondante lorsqu’elle est inutile et
superflue. La droite représentative de cette contrainte
se trouve entièrement en dehors de la région
admissible.
max z x1 2 x2 (1)
Exemple : (3)
sc
x1 x2 8 (1)
5 x 4 x2 20 ( 2)
1
x1 3 (3) (2)
x2 4 (4)
(4)
x1 , x2 0
La contrainte (1) est une contrainte redondante
car (3)+(4) donne : x1 x2 7; 7 8
Donc la contrainte (1) est automatiquement
vérifiée et donc sera éliminée.
2. Cas de dégénérescence : il y a deux types
Exemple : (4)
(3)
maximiser z 2 x1 2 x2 (2)
sujet à
2 x1 x2 4 (1)
x x2 8 (2)
1
x1 x2 4 (3)
x 5 ( 4)
1
x1, x2 0
Dégénérescence du 2ième type :
une solution optimale est dite dégénérée si
plus de deux contraintes concourent en ce point.
Exemple :
(1)
max z x1 x2 (2)
sc
3x1 2 x2 40 (1)
x 10
1 (2)
x2 5 (3)
x1, x2 0
B=(10,5)
La solution optimale B est dégénérée (3)
Exercice 1
• Un restaurateur dispose de trois types de poissons
• 30 p1, 24 p2, 18 p3 et désire offrir :
•
des plats à 80 UM contenant 5 p1, 2 p2 et 1 p3
Problème :
max z 80 x1 60 x2
sc
5 x1 3 x2 30
2 x 3 x2 24
1
x1 3 x2 18
x1 , x2 0
Résolution graphique
A est un point extrême de la région
admissible dont les coordonnées vérifient max z 80 x1 60 x2
la 1ière et la 3ième contrainte avec égalité
x2
sc
A =(3,5) est une
solution optimale unique 5 x1 3 x2 30
2 x 3 x2 24
et zmax=540 UM 1
10 x1 3 x2 18
8 x1 , x2 0
6
A
0
x1
6 12 18
Interprétation
Compréhension du problème :
• Revenus hebdomadaires
27 x1 21x2
• -
• Coûts hebdomadaires d’achat des matières premières
10 x1 9 x2
•
• -
• Autres coûts variables hebdomadaires
14 x1 10 x2
• d’où la fonction économique à maximiser est :
z 3x1 2 x2
Mise en équations des contraintes économiques : Les variables
x1 et x2 sont limitées par 3 contraintes
max z 3 x1 2 x2
sc
2 x1 x2 100
x x2 80
1
x1 40
x1 , x2 0
Résolution graphique
x2
x1 = 40
max z 3 x1 2 x2
sc
2 x1 x2 100
x x2 80
1
x1 40
x1 , x2 0
x1
0,3x1 0,5 x2 180 10 4
x1 , x2 0
Résolution graphique
x2 106
min w 700 x1 500 x2
sc
0,3x1 0,25 x2 125 104
5
0,4 x1 0,25 x2 135 104
4
A
0,3x1 0,5 x2 180 104
3
x1 , x2 0
x1 106
1 2 3 4 5 6
recette
max z 140 x1 150 x2
Mise en équations des contraintes économiques :
les variables x1 et x2 sont limitées par 3 contraintes de
type contrainte de stockage :
a. Contrainte d’essence
0,2 x1 0,4 x2 1200
b. Contrainte de gazole
0,4 x1 0,2 x2 1200
c. Contrainte de fuel
0,4 x1 0,4 x2 1400
A
2,5
x1 103
1
Interprétation
La solution qui maximise le profit est A(1000,2500), c’est-à-dire,
en traitant 1000t de brut A et 2500t de brut B, le raffineur
réalisera un recette maximale de : 140 1000 150 2500 515 000 UM
A l’optimum, la première et la troisième contrainte sont saturées
tandis que la deuxième contrainte ne l’est pas : pour réaliser la
recette maximale, le raffineur doit vendre tout le stock de l’essence
et du fuel lourd :
1200 0,2 1000 0,4 2500 et 1400 0,4 1000 0,4 2500
toutefois, il lui restera un stock du gazole de 300 t :
A B C D E F Ligne 1 : titre de
l’exercice ou référence
1 Exercice 4
Ligne 3 : noms des
2 variables de décision, ici
3 variables de brut A et brut B
Brut A Brut B
décision
Ligne 4 : coefficients des
4 fonction Profit variables de décision dans
140 150 TOTAL la fonction objectif
objectif unitaire
5 quantité 0 0
Ligne 5 : quantités
6 profit =C4*C5 =D4*D5 =C6+D6 cherchées. Elles seront
remplies par le solveur
7 limites après résolution. Elles
sont nulles avant
8 contraintes essence 0.2 0.4 =C8*C5+D8*D5 1200 résolution
Spécification des
cellules variables :
tapez les références ou
$E$6 les noms des cellules
devant être modifiés par
le Solveur jusqu'à ce que
les contraintes du
problème soient
$C$5:$D$5 respectées et que la
cellules cible atteigne le
résultat recherché.
Spécification des
$E$8<=$F$8 contraintes :A l'aide
$E$9<=$F$9 des boutons ajouter,
Modifier et Supprimer
$E$10<=$F$10 de la boîte de dialogue
Paramètres du Solveur,
établissez votre liste de
contraintes dans la zone
Contraintes.
Exercice 4: Feuille Excel
Les options du Solveur
Résultat et résolution
Une fois tous les paramètres
du problème sont mis en
place, le choix du bouton
Résoudre amorce le
processus de résolution du
problème . Vous obtenez
alors une de ces réponses
Résolution par Excel et
Analyse de la solution optimale
Microsoft Excel 10.0 Answer Report
Worksheet: [exercice 4.xls]probleme primal
Report Created: 11/11/2006 11:13:28 PM
Adjustable Cells
La solution optimale
Cell Name Original Value Final Value
$C$5 quantités Brut A 0 1000
$D$5 quantités Brut B 0 2500
Constraints
Le stock de gazole
Cell Name Cell Value Formula Status Slack
$E$8 essence TOTAL 1200 $E$8<=$F$8 Binding 0
$E$9 gazole TOTAL 900 $E$9<=$F$9 Not Binding 300
$E$10 fuel TOTAL 1400 $E$10<=$F$10 Binding 0
Exercice
Adjustable Cells
Cell Name Original Value Final Value
$C$5 quantités Article A 0 150 La solution optimale
$D$5 quantités Article B 0 100
$E$5 quantités Article C 0 150
Le stock
Constraints
Cell Name Cell Value Formula Status Slack
$F$8 contrinte horaire TOTAL 120 $F$8<=$G$8 Binding 0
$F$9 limite de production de l'ensemble TOTAL 400 $F$9<=$G$9 Not Binding 100
$F$10 limite de production de A TOTAL 150 $F$10<=$G$10 Not Binding 150
$F$11 limite de production de B TOTAL 100 $F$11<=$G$11 Binding 0
$F$12 limite de production de C TOTAL 150 $F$12<=$G$12 Binding 0
Interprétation
La solution qui maximise le bénéfice est (150,100,150), c’est-à-dire, en fabriquant 150
t de papier A, 100 t de papier B et 150 t de papier C, l’entreprise X réalisera un
bénéfice maximal de : 140000 UM.
200 150 500 100 400 150 140 000 UM
A l’optimum, la première, la quatrième et la cinquième contrainte sont saturées
tandis que la deuxième et la troisième contrainte ne le sont pas : pour réaliser le
bénéfice maximal, l’entreprise doit épuiser toutes les disponibilités horaires et vendre
la totalité de la production en papier de type B et C :
120 150
5
100 150
2, 5 3
et 100 t 1100 et 150 t 1150
toutefois, elle lui restera un stock de papier de type A de 150 t :
Recherche opérationnelle
Partie II
Non Oui
Résolution graphique Résolution par la
Forme canonique méthode du simplexe
Forme standard
Fonction objectif
n
max (ou min) z c1 x1 c2 x2 cn xn c j x j
j 1
Contraintes technologiques
n
a x j , , b2
j 1
2j
am 1 x1 am 2 x2 amn xn , , bm
n
xj 0 , j 1,2,, n
xj variables de décision (inconnues)
aij, bi, cj paramètres du programme linéaire
Écriture matricielle de la forme canonique :
recherche d’un maximum
•Considérons toutes les contraintes du type et on suppose que b 0
(hypothèse non générale)
max z t cx
n
max z c j x j
j 1
sc sc
n
aij x j bi , 1 i m Ax b
j 1
•avec, x 0, 1 j n
j x 0
c1 x1
a11 a1n b1
A
b c x
am1 amn bm
cn xn
Vecteur colonne de Vecteur ligne de taille n Vecteur colonne
Matrice (m,n) taille m<n de taille n
des coefficients de la
des coefficients techniques des seconds membres fonction économique des niveaux
des contraintes d’activité
Exercice
• L’entreprise X occupe trois ouvriers à raison de 40 heures par
semaine chacun. Elle fabrique trois types de papier A, B, C.
•Pour des raisons d’approvisionnement en matières premières, la production
totale hebdomadaire de cette entreprise est limitée à 500 tonnes.
•Les rendements horaires pour la fabrication des papiers A, B, C sont
respectivement de 5 tonnes, 2.5 et 3 tonnes. Les bénéfices sont
respectivement de 200, 500 et 400 UM par tonne de A, B et C.
•L’expérience a montré que, si la production était limitée respectivement à
300, 100 et 150 tonnes par semaine, elle aurait été entièrement absorbée
par le marché.
•Objectif :
• On veut déterminer les quantités hebdomadaires à produire pour
que le bénéfice soit maximum et toute la production soit vendue.
Exemple
La solution optimale
Adjustable Cells
Cell Name Original Value Final Value
$C$5 quantités Article A 0 150
$D$5 quantités Article B 0 100
Le stock
$E$5 quantités Article C 0 150
Constraints
Cell Name Cell Value Formula Status Slack
$F$8 contrinte horaire TOTAL 120 $F$8<=$G$8 Binding 0
$F$9 limite de production de l'ensemble TOTAL 400 $F$9<=$G$9 Not Binding 100
$F$10 limite de production de A TOTAL 150 $F$10<=$G$10 Not Binding 150
$F$11 limite de production de B TOTAL 100 $F$11<=$G$11 Binding 0
$F$12 limite de production de C TOTAL 150 $F$12<=$G$12 Binding 0
Interprétation
La solution qui maximise le bénéfice est (150,100,150), c’est-à-dire, en fabriquant 150
t de papier A, 100 t de papier B et 150 t de papier C, l’entreprise X réalisera un
bénéfice maximal de : 140000 UM.
15 1
2,5
1
3 120
1 1 500
1 200 x1
A 1 0 0 b 300 c 500 x x2
0 1 0 100 400 x3
0 0 1 150
Écriture matricielle de la forme
canonique : recherche d’un minimum
•Considérons toutes les contraintes du type et on suppose queb 0
(hypothèse non générale)
min w t cx
n
min w c j x j c1 x1
j 1
a11 a1n b1
sc sc
A , b , c , x
n
aij x j bi , 1 i m Ax b
j 1 am1 amn bm
x 0, 1 j n
j
x 0 cn xn
i 1,2,, m a x
j 1
ij j bi , bi 0 a x
j 1
ij j ei bi , ei 0
Fonction objectif
max z c1 x1 c2 x2 cn xn 0e1 0e2 0em
A x + Ime = b
Contraintes
a11x1 a12 x2 a1n xn e1 0 0 b1
technologiques a21x1 a22 x2 a 2 n xn 0 e2 b2
Ax 0
am1 x1 am 2 x2 amn xn 0 0 em bm
Contraintes de non négativité Im e
x j 0 j 1,2,, n
ei 0 i 1,2,, m
Écriture matricielle de la forme standard
t
c x max z c~~
t
x
0
max z e
~~
Ax b
x
sc ~
A I m
eb x 0
sc
x 0
c x
A A | Im
e 0 ~ ~ ~
où c x 0
0 e
~
• On a m équations. Généralement les m vecteurs colonnes de A sont
linéairement indépendants
Variables d’écart
xN
~~
Ax b N B
x
b
B
•On pose :x N 0 et en multipliant les deux termes par B 1 xB B 1b
•Comme Im est inversible, on a une solution de base admissible initiale (x =0, l’origine)
x
•Pour une base donnée : ( A | I m ) b xN x 0 et xB e b
e
max z cx 0e max z c N x N cB xB
max z c N x N cB xB
A | I b x
x
sc sc N | I N
b sc Nx N xB b
e xB
•D’où, xB b NxN xB i bi Ni xN
Solution de base : xB vecteur déterminé en :
Solution optimale :
solution de base admissible qui maximise la fonction objectif ;
B
admissible initiale :
e1 3 4 1 0 160
On attribue une valeur nulle à la
e2 6 3 0 1 180 marge globale : les variables Hors
Base (HB) sont nulles x1=0 et x2 =0,
c’est la solution de base admissible
-z 1200 1000 0 0 0
de départ
coefficients Coefficients de la valeur de la on fabrique 0 pièces P1 et 0
fonction économique
techniques fonction économique pièces P2 : aucun intérêt pratique
Les valeurs des variables dans la
La solution de base admissible initiale est : Base (B) sont données par :
(x1, x2 ,e1 ,e2) = (0, 0, 160, 180)
e1 160 3x1 4 x2 160
La dernière cellule donne la valeur de –z. Donc la
e2 180 6 x1 3x2 180
marge z =0. en effet : z 1200 x1 1000 x2 0
Elles se lisent dans la colonne b
La solution n’est pas optimale. On recherche donc
une solution de base meilleure : autre itération. e1 = 160 et e2 = 180
La dernière ligne donne les valeurs marginales (ou
il reste 160 h d'utilisation
taux de substitution) des variables hors base possible de A1 et 180 h de A2
x1 x2 e1 e2 b R
Déterminons une autre e1 3 4 1 0 160 160/3
solution de base e2 6 3 0 1 180 180/6
admissible
-z 1200 1000 0 0 0
La solution de base admissible initiale est (x1,x2,e1,e2)=(0,0,160,180) avec z =
0
La dernière ligne donne les valeurs marginales ou taux de substitution :
•Si x1 = 1, x2 = 0, e1 = 160, e2 =180 alors, z = 1200 dhs
•une augmentation de 1 unité de x1 ferait croître la fonction objectif de 1200
dhs,
•Si x1 = 0, x2 = 1, e1 = 160, e2 =180 alors, z = 1000 dhs
•une augmentation de 1 unité de x2 ferait croître la fonction objectif de 1000
dhs.
On a intérêt à augmenter la valeur de la fonction objectif le plus rapidement possible
donc à augmenter la variable ayant le plus grand coefficient strictement positif (cas
de maximisation) de la dernière ligne : x1 variable entrante dans la base
Supposons x1 augmente et x2 = 0, alors 6 x1 3x2 e2 180 x1 30
z 1200 x1 1000 x2 1200 30 36000
On a intérêt à prendre le maximum de x1 (les variables hors base restant nulles) :
x1=30. La deuxième contrainte sera saturée e2=0 : e2 variable sortante de (B).
z = 36000 dhs, en effet :
On augmente la fonction objective en faisant entrer une variable dans la base
prenant la place d'une variable qui va sortir de la base.
Critère de sélection
Variable entrante dans (B) Variable sortante de (B)
on exprime la fonction objectif en fonction des On effectue le rapport des seconds membres
seules variables hors base et on choisit la variable des contraintes aux coefficients strictement
pondérée par le cœfficient le plus élevé strictement positifs correspondants de la variable entrante
positif : on sélectionne la variable HB ayant le plus : on sélectionne la variable de la base ayant le
grand coefficient strictement positif dans la plus petit rapport positif dans la colonne R
dernière ligne
Le rapport des coefficients de la colonne C par les
Règle
coefficients strictement positifs de la colonne de la variable
Le pivot doit être égal à 1
entrante dans la base (x1)
Les coefficients de la ligne du
pivot sont divisés par le pivot
Le pivot est égal à 6
Ligne de pivot
e2 sortant dans la base
Les coefficients de la colonne
du pivot (sauf le pivot) sont nuls
xx11 HB
B x2 e1 e2 b R Les autres coefficients sont
e1 30 5/2
4 1 -1/2
0 70
160 160/3 obtenus par la règle du rectangle
e2 161 1/2
3 00 1/6
1 30
180 180/6
pivot Lpivot
-z 0
1200 1000
400 00 0
-200 0
-36000
x1 entrant de la base
Colonne de pivot Cpivot Av
Nv : nouvelle valeur Av : ancienne valeur
Nv Av (Cpivot Pivot ) Lpivot
Cpivot : colonne pivot Lpivot : ligne pivot
HB x1 x2 e1 e2
B
b R Les variables hors base sont nulles :
e1 3 4 1 0 160 160/3 x2 = 0, e2 = 0
e2 6 3 0 1 180 30
on fabrique 0 pièces P2 et il reste
0 h d'utilisation disponible à A2
-z 1200 1000 0 0 0
La contrainte associée à e2 est
Tableau 1
saturée.
B HB x1 x2 e1 e2 b
La nouvelle solution de base réalisable
e1 0 5/2 1 -1/2 70
x1 1 1/2 0 1/6 30 est : (x1, x2, e1, e2) = (30, 0, 70, 0)
-z 0 400 0 -200 -36000 on fabrique 30 P1 et il reste 70h
d'utilisation disponibles à A1
La dernière ligne donne les valeurs marginales :
180 e2 3x2
z 1200 1000 x2 36000 200e2 400 x2 la marge est égale à z =36 000 dhs
6
•Une augmentation de 1 unité de x2 ,ici on augmente la production de 1 pièce de P2, ferait
croître la fonction objectif de 400 dhs, et une augmentation de 1 unité de la variable d’écart
e2 (c- à –d, diminution du second membre de l'équation correspondante de 1 unité, ici on
diminue la disponibilité de 1 h à A2) ferait diminuer la fonction objectif de 200 dhs
On a intérêt a augmenter x2 : x2 variable entrante dans la base
e1 70 5 2 x2 0 x 28
2 x2 28
x1 30 1 2 x2 0 x2 60
e1 est la variable sortante de la base
x1 x2 e1 e2 b R
e1 0 5/2 1 -1/2 70 28 Les variables hors base sont nulles
x1 1 1/2 0 1/6 30 60 (Les contraintes associées sont
-z 0 400 0 -200 -36000
saturées) : e1=0 et e2=0
il reste 0 h d'utilisation
•Tableau final disponible aux ateliers A1 et A2
x1 x2 e1 e2 b
x2 0 1 2/5 -1/5 28 La nouvelle solution de base réalisable
x1 1 0 -1/5 4/15 16
est : (x1, x2, e1, e2) = (16, 28, 0, 0)
-z 0 0 -160 -120 -47200
on fabrique 16 pièces P1 et 28
pièces de P2
Pour augmenter e1 d’une unité, il faut : Et la marge est égale à 47 200 dhs
• diminuer x2 de 2/5
• diminuer x1 de (-1/5)
La variation correspondante de la fonction économique : -(-1/5) 1200-(2/5) 1000 = -160
D’où, augmenter e1 d’une unité diminuerait la fonction objectif de 160 dhs
Pour augmenter e2 d’une unité, il faut :
• diminuer x2 de (-1/5)
• diminuer x1 de 4/15
La variation correspondante de la fonction économique :
-(4/15) 1200-(-1/5) 1000 = -120
D’où, augmenter e2 d’une unité diminuerait la fonction objectif de 120 dhs
On retrouve ainsi la solution optimale de la résolution graphique.
Interprétation graphique de la
méthode du simplexe
A = (0,0) et z = 0 dhs
50
C = (16,28) et z = 47 200 dhs
40
30 C
20
Donc C est l’unique solution optimale et
A 10 20 30 40 50 60
Critères d'arrêt des itérations de la méthode du simplexe lors
de la résolution d’un problème de maximisation
Si tous les coefficients de la dernière ligne, relatifs aux variables HB, sont négatifs ou nuls,
l’algorithme s’arrête et la solution trouvée est optimale.
x1 x2 e1 e2 b
x2 0 1 2 -1 20
x1 1 0 -1 2,5 6
-z 0 0 -160 -80 -40
S'il existe une variable HB (non artificielle) ayant un coefficient positif dans la dernière ligne et
telle que tous les coefficients correspondants dans le tableau soient nuls ou négatifs, alors la
solution est infinie
x1 x2 e1 e2 b
e1 0 -2 1 -1 20
x1 1 0 0 2,5 3
-z 0 16 0 -3 -32
•Si, à la fin des itérations, une variable est HB avec un coefficient nul dans la dernière ligne, alors
on a une arête (plan …) optimale. Les autres sommets solutions sont obtenus en faisant rentrer
cette variable dans la base. x x e e b
1 2 1 2
x2 0 1 2/5 -1/5 28
x1 1 0 -1/5 4/15 16
-z 0 0 0 -120 -47200
Exercice 7
déjà traité par la méthode graphique et par le solveur d’Excel
Brut A B
Essence 0,2 0,4
Gazole 0,4 0,2
Fuel lourd 0,4 0,4
Le traitement d’une tonne de brut A procure une recette de 140 UM Le
traitement d’une tonne de brut B procure une recette de 150 UM
Du fait des contraintes de stockage, la fabrication de chaque produit est
limitée de la manière suivante : Essence : 1 200 t, Gazole : 1 200 t, Fuel
lourd : 1 400 t.
Quelles quantités de ces pétroles bruts devra t-on traiter pour maximiser
la recette ?
Modélisation mathématique
Variables économiques ou de décision :
x1 = quantité de brut A à traiter
x2 = quantité de brut B à traiter
La forme standard la forme canonique
x1 x22 e1 e2 e3 b
e1 0,2
0,5 0,4
1 1
2,5 0
0 0 1200
3000
e2 0,4
0,3 0,2
0 0
-0,5 1 0 1200
600
e3 0,4
0,2 0,4
0 0
-1 0 1 1400
200
-z 65
140 150
0 0
-375 0 0 0
-450000
x1 x2 e1 e2 e3 b
Tableau final
x2 0 1 5 0 -2,5 2500
e2 0 0 1 1 -1,5 300
x1 1 0 -5 0 5 1000
astucieuse.
intéressantes
Passage du programme primal au programme
dual et interprétation économique du dual
•
• A tout programme linéaire primal correspond un programme linéaire dual
max z c1 x1 c2 x2 cn xn
a11x1 a12 x2 a1n xn b1 min w b1 y1 b2 y2 bm ym
a x a x a x b
21 1 22 2 2n n 2
a11 y1 a21 y2 am1 ym c1
sc a y a y a y c
a x a x a x b
12 1 22 2 m2 m 2
m1 1 m2 2 mn n m sc
1 2
x , x , , x n 0 a y a y a y c
1n 1 2n 2 mn m n
y1 , y2 , , ym 0
La variable de décision xj en base (xj > 0) La variable d’écart tj hors base (tj = 0)
La variable d’écart ei en base (ei > 0) La variable de décision yi hors base (yi = 0)
La variable d’écart ei hors base (ei = 0) La variable de décision yi en base (yi > 0)
Signe
Ligne xj en base Colonne tj hors base
opposé
Signe
Ligne ei en base Colonne yi hors base
opposé
Signe
Colonne xj hors base Ligne tj en base
opposé
Signe
Colonne ei hors base Ligne yi en base
opposé
Signe
Colonne second membre Taux marginaux de substitution
opposé
Taux marginaux de substitution des Signe Colonne second membre des variables de
variables hors base opposé base associées
Exercice 6: déjà traité
Problème de maximisation - Problème Primal
•Une usine fabrique 2 pièces P1 et P2 usinées dans deux ateliers A1 et A2 . Les temps
d'usinage sont :
pour P1 : de 3 heures dans l'atelier A1 et de 6 heures dans A2
y1 y2 max
x1 3 6 1200
x2 4 3 1000
Dual : écriture en ligne
min 160 180
max z 1200 x1 1000 x2 0e1 0e2 min w 160 y1 180 y2
3x1 Résolution
4 x2 e1 0 160 3 y1 6 y2 1200
sc 6 x1 3 x2 0 e2
par 180 sc 4 y1 3 y2 1000
x , y ,
1 x2 , e1 ,
simplexee2 0 1 y2 0
•Tableau final du programme primal
HB
x1 x2 e1 e2 b
B
Interprétation : A l’optimum,
x2 0 1 2/5 -1/5 28
La variable de base x1 = 16 donc le taux marginal de
x1 1 0 -1/5 4/15 16
substitution de la variable d’écart t1 est égal à -16 et la
-z 0 0 -160 -120 -47200
valeur de t1 est nulle (t1 hors base)
La variable de base x2 = 28 donc le taux de substitution
•Tableau final du programme dual de la variable d’écart t2 est égal à -28 et la valeur de t2
est nulle (t2 hors base)
HB
y1 y2 t1 t2 b
B Le taux marginal de substitution de la variable d’écart
- hors base e1 est égal à -160 donc la variable de base
y2 0 1 1/5 120
4/15
y1 = 160
y1 1 0 1/5 -2/5 160
Le taux marginal de substitution de la variable hors base
-w 0 0 -16 -28 -47200
e2 est égal à -120 donc la variable de base y2 = 120
Les deux contraintes du dual sont saturées, donc les
•A l’optimum, on a : deux valeurs optimales du dual sont non nulles, ainsi, les
z = w = 47200 dhs contraintes du primal correspondantes sont saturées.
(x1 , x2) = (16 , 28)
(y1 , y2) = (160 , 120)
Exercice 4: déjà traité
Problème de maximisation - Problème Primal
Quelles quantités des pétroles Bruts devra t-on traiter pour maximiser la recette ?
x1 : quantité de brut A à traiter
x2 : quantité de brut B à traiter
y1 y2 y3 max
Dual : écriture en ligne
max z 2 x1 3 x2
x1 x2 1
2 x 3 x2
1 4
PLG sc
x1 2 x2 3
x1 , x2 0
sa forme équivalente :
max z 2 x1 3 x2
x1 x2 e1 1
2 x 3 x2 e2 4
PL sc 1
x1 2 x2 3
x1 , x2 , e1 , e2 0
Ce problème sera équivaut à (PL =) que si tous les ai sont nuls (ici, a2=a3 =0)
Soit donc le modèle (PLa) où la valeur de a2 a3 doit être minimisée
Ainsi pour la résolution de (PLG) on doit passer par deux étapes :
• Phase 1 min za a2 a3 Recherche d’une solution de base admissible initiale.
•.Phase 2 max z 2 x1 3x2 Recherche de la solution optimale.
Phase I : Min
0 0 0 0 1 1
cj
min za a2 a3 Base
Valeur Limite
x1 x2 e1 e2 a2 a3
x1 x2 e1 1
Coeff Var
2 x 3x e a 4
0 e1 1 -1 1 0 0 0 1 --------
1 2 2 2 1 a2 2 3 0 -1 1 0 4 4/ 3
( PLa )
x1 2 x2 a3 3 1 a3 1 2 0 0 0 1 3 3/ 2
zj 3 5 0 -1 1 1
x1 , x2 , e1 , e2 , a2 , a3 0 cj - zj -3 -5 0 1 0 0
7
Min cj 0 0 0 0 1 1
Forme canonique Base Valeur Limite
x1 x2 e1 e2 a2 a3
min za 3x1 5 x2 e2 7 Coeff Var
0 e1 5/3 0 1 -1/3 1/3 0 7/3 --------
x1 x2 e1 1 0 x2 2/3 1 0 -1/3 1/3 0 4/3 --------
2 x 3x e a 4
1 2 2 2
1 a3 -1/3 0 0 2/3 -2/3 1 1/3 1/2
( PLa ) zj -1/3 0 0 2/3 -2/3 1
x1 2 x2 a3 3 1/3
cj - zj 1/3 0 0 -2/3 5/3 0
x1 , x2 , e1 , e2 , a2 , a3 0 Min cj 0 0 0 0 1 1
Base Valeur
x1 x2 e1 e2 a2 a3
Coeff Var
On obtient ainsi a2 a3 0 0 e1 3/2 0 1 0 0 1/2 5/2
Définitions et Applications
II- Méthodes graphiques d’optimisation
La gestion de projet revêt une importance vitale pour l’entreprise. Celle-ci est
appelée pour survivre, à faire des investissements qui nécessitent une
planification rigoureuse, de l’ordonnancement des activités et des coûts et des
charges,…
A- La planification
La planification permet de définir le plan directeur du projet qui décrit les
objectifs et la nature de l’organisation à mettre en place. Elle détermine avec
précision, tout ce qui doit être fait ( et comment et par qui) en respectant toutes
les contraintes liées au temps, aux coûts et à la qualité. Elle comporte 4
dimensions:
47
1- Méthode de GANTT (1918)
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
A
B
C
D
E
Exemple 2 :
Une entreprise veut fabriquer un prototype de vélo
à moteur. Pour cela plusieurs activités sont
nécessaires selon le bureau méthode.
TAF : tracez le diagramme de GANTT sachant que
la fabrication du prototype doit commencer le 3
octobre.
Tâches Tâches Antérieures Durée en jours
A ---- 2
B ---- 1
C E;B;H 1
D C 2
E A 1
F E;B;H 2
G D;F 1
H ---- 3
3 Oct 4 Oct 5 Oct 6 Oct 7 Oct 8 Oct 9 Oct
A
B
Flottement
C
D
E
F
G
H
5
1. d’établir le graphique Gantt et MPM;
B A
C 4 B
2- de mettre en évidence le chemin D 4 C
critique; E 5 D,F
C 4 B
D 4 C
E 5 D,F
F 16 C
G 8 B
H 4 G
I 6 H,F
J 12 I,E
K 7 J,M
L
15 ---
M 5 L
Définition d’un graphe orienté
• Un graphe orienté G est un couple (X,A)
• où
• X : ensemble de sommets {x1,..., xn}
• A : ensemble de couples orientés (xi,xj) appelés arcs
X3
• X = {x1 , x2 , x3 , x4 , x5 , x6} X5
X4
X1
X5
•x6 •x4, x5 •-
Détermination des niveaux des sommets d’un graphe
sans circuit
Dans un graphe sans circuit, le niveau d'un sommet x est la longueur du plus long
chemin ayant pour extrémité x
• No={sommets de niveau 0}
Sommets x Précédents P(x) • ={sommets n’ayant pas de
précédents}
x1 -
x2 x1, x3 • ={x1}
x3 x1
x4 x2, x3 • Tous les sommets x1 sont barrés
x5 x4
x6 x4, x5 • Les sommets barrés sont considérés
• comme n’existants plus
Sommets x Précédents P(x)
x1 - • N1={sommets de niveau 1}
x2 x1, x3 • ={sommets n’ayant pas de
x3 x1 précédents}
x4 x2, x3 • ={x3}
x5 x4
x6 x4, x5
• Tous les sommets x3 sont barrés
Sommets x Précédents P(x) • N2={sommets de niveau 2}
x1 - • ={sommets n’ayant pas de
précédents}
x2 x1, x3
• ={x2}
x3 x1
x4 x2, x3 • Tous les sommets x2 sont barrés
x5 x4
x6 x4, x5 • Les sommets barrés sont considérés
comme n’existants plus
Sommets x Précédents P(x)
x1 - • N3={sommets de niveau 3}
x2 x1, x3 • ={sommets n’ayant pas de
précédents}
x3 x1
• ={x4}
x4 x2, x3
x5 x4 • Tous les sommets x4 sont barrés
x6 x4, x5
• N4={sommets de niveau 4}
Sommets x Précédents P(x)
• ={sommets n’ayant pas de
x1 - précédents}
x2 x1, x3 • ={x5}
x3 x1
x4 x2, x3 • Tous les sommets x5 sont barrés
x5 x4
x6 x4, x5 • Les sommets barrés sont considérés
comme n’existants plus
Sommets x Précédents P(x)
x1 - • N5={sommets de niveau 5}
x2 x1, x3 • ={sommets n’ayant pas de
précédents}
x3 x1
x4 x2, x3 • ={x6}
x5 x4
x6 x4, x5 • Tous les sommets x6 sont barrés
• N4 = { x5 }
Tous les sommets ayant été barrés
• N5 = { x6 }
• Les niveaux sont :
Sommets x Précédents P(x)
Utilité
• N
de la notion de niveaux
0 = { x1x1}
x2
N1 = { -x3 }
N = {x }
x1, x3
2 2
x3 x1
• N3 = { x4x } N4 = x{ ,xx5 } N5 = { x6 }
4 2 3
x5 x4
x6 x4, x5
construire le graphe ordonné par niveaux
La recherche des chemins optimaux se fait plus facilement sur
un graphe ordonné par niveaux
x1 x3 x4 x5
x6
x2
N0 N1 N2 N3 N4 N5
Définitions
• Un chemin est une suite ordonnée (x1,...,xn) de sommets reliés par des
arcs :
i 1,, n 1 xi 1 S xi
• Un circuit est un chemin (x1,...,xn) tel que x1 = xn
• Dans un graphe sans circuit. A chaque arc (x,y) est associé un nombre
positif V(x,y) appelé la valeur de l'arc
V x , x
i 1
i i 1
Cas particuliers de chemins
• Chemin près-hamiltonien : il passe au moins une fois par chaque
sommet du graphe
• (x1, x2, x4, x3, x2, x6, x8, x7, x5)
x2 x4 x8
x5
x1
x3 x6 x7
Cas particuliers de chemins
• Chemin près-hamiltonien : il passe au moins une fois par chaque
sommet du graphe
• (x1, x2, x4, x3, x2, x6, x8, x7, x5)
• Chemin hamiltonien : il passe une et une seule fois par chaque
sommet du graphe
• (x1, x2, x4, x3, x6, x8, x7, x5)
• Chemin près-eulérien : il passe au moins une fois par chaque arc du
graphe
• (x1, x2, x3, x4, x8, x7, x4, x2, x3, x4, x2, x6, x7, x5)
• Chemin eulérien : il passe une et une seule fois par chaque arc du
graphe
• (x1, x2, x3, x4, x8, x7, x4, x2, x6, x7, x5)
x2 x4 x8
x5
x1
x3 x6 x7
Détermination d'un chemin de valeur optimale
entre les sommets D et F
• L'algorithme de Ford
3. On supprime les sommets et les arcs par lesquels on ne peut pas passer pour
aller de D à F
3 F
7
2
3
4 3
6 6
5
3 6 9
1 9 8
5 1
9 8
2 4
D
4 5
Les problèmes d’ordonnancement
•L’objectif est de :
Une tâche x ne pouvant débuter que lorsque toutes les tâches qui y
aboutissent sont terminées.
On rajoute au graphe un sommet terminal permettant de dater la fin
des travaux
La représentation graphique est ordonnée par niveaux des sommets
(des tâches)
5
T1 T2
4
T3 8 T4
Opération x P(x)
O1 --
O2 O1
O3 --
O4 --
O5 O 2, O 3, O 4
O6 O2, O3
O7 O1
O8 O 5, O 7
O9 O 6, O 8
• N0={O1,O3,O4}
Méthode des Potentiels Metra
Détermination du niveau 1:
Opération x P(x)
O1 --
O2 O1
O3 --
O4 --
O5 O 2, O 3, O 4
O6 O2, O3
O7 O1
O8 O 5, O 7
O9 O 6, O 8
• N1={O2;O7}
Méthode des Potentiels Metra
Détermination du niveau 2:
Opération x P(x)
O1 --
O2 O1
O3 --
O4 --
O5 O 2, O 3, O 4
O6 O2, O3
O7 O1
O8 O 5, O 7
O9 O 6, O 8
• N2={O5, O6}
Méthode des Potentiels Metra
Détermination du niveau 3:
Opération x P(x)
O1 --
O2 O1
O3 --
O4 --
O5 O 2, O 3, O 4
O6 O2, O3
O7 O1
O8 O 5, O 7
O9 O 6, O 8
• N3={O8}
opérations durée opérations
Représentation graphique (tâches) (mois) antérieures
O1 4 -
O2 6 O1
Niveau 0 Niveau 1 Niveau 2 Niveau 3 Niveau 4
{O1,O3,O4} {O2,O7} {O5,O6} {O8} {O9} O3 4 -
O4 12 -
• O5 10 O2,O3,O4
O6 24 O2,O3
O7 7 O1
O8 10 O5,O7
? ? 4 ? ? 6 ? ?
O9 3 O6,O8
o1 o2 o6
4 4 24
6
? ? ? ? 7 ? ? ? ? 3 ? ?
o3 o7 o8 10 o9 F
4 10
? ? ? ?
o4 12 o5
Calendrier au plus tôt des tâches
• Tx est la date au plus tôt correspondant à la valeur du chemin de valeur
maximale aboutissant à x (algorithme de Ford)
On commence par les sommets de niveaux les plus faibles jusqu’aux
sommets de niveaux les plus élevés
0 ? 4 4 ? 6 10 ?
o1 o2 o6
4 4 24
6
0 ? 4 ? 7 22 ? 34 ? 3 37 ?
o3 o7 o8 10 o9 F
4 10
0 ? ?
Le chemin de valeur12maximale (durée 37 mois)
o4 aboutissant
12 à F est :o5 (O1, O2, O6, O9)
Calendrier au plus tard des tâches
• T*x est la date au plus tard à laquelle peut commencer une tâche
sans remettre en cause la date de fin des travaux
On commence par les sommets de niveau les plus élevés jusqu’aux
sommets de niveau les plus faibles
T*y1
y1
T*x
x T*y2
y2
Détermination des calendriers au plus tard de la réalisation de
chacune des tâches :
T*F = TF= 37
T*9 = T*F - V(9,F) = 37 - 3 = 34
T*8 = T*9 - V(8,9) = 34 - 10 = 24
T*6 = T*9 - V(6,9) = 34 - 24 = 10
Sur les tâches critiques
•
T*5 = T*8 - V(5,8) = 24 - 10 = 14
on a :
T*7 = T*8 - V(7,8) = 24 - 7 = 17
T*2 = Min [T*5 - V(2,5) ; T*6 - V(2,6)] = Min [14 - 6 ; 10 - 6] = Min [8 ; 4] = 4
T* = Tx
T* = Min [T* - x
1 V(1,2) ; T* - V(1,7)] = Min [4 - 4 ; 17 - 4] = Min [0 ; 13] = 0
2 7
T*3 = Min [T*6 - V(3,6) ; T*5 - V(3,5)] = Min [10 - 4 ; 14 - 4] = Min [6 ; 10] = 6
T*4 = T*5 - V(4,5) = 14 - 12 = 2
0 0 4 4 4 6 10 10
o1 o2 o6
4 24
4 6
0 6 4 17 7 22 24 34 34 3 37 37
o3 o7 o8 10 o9 F
10
4
0 2 12 14
o4 12 o5
Il y‘a deux types de retard relatif à l’exécution
des tâches sans retarder l’achèvement du projet
x
Ty2
y2
Marge totale Marges libres
• mL(x) = min [Ty - Tx - V(x,y)]
• mt(x) = T*x - Tx • Le min étant pris sur les suivants y de x
• mL(O1) = Min[T2-T1-V(1,2);T7-T1-V(1,7)]
• mt(O1) = T* 1 – T1 = 0 - 0 = 0 0• 0 4 4 4 =6 Min[0;0]=0
10 10
24
o1 o2 o6
• mt(O2) = T*2 – T2 = 4 - 4 = 0 • m (O ) = Min[T6-T2-V(2,6);T5-T2-V(2,5)]
4 L 42
• = Min[0;2]=0
6
• mt(O3) = T*3 – T3 = 6 - 0 = 6 0 6 4 17 7
22 24 34 34 3 37 37
•o3 mL(O3)o=7 Min[T6-T3-V(3,6);To58 -T10 o9
3-V(3,5)]
F
• mt(O4) = T* 4 – T4 = 2 - 0 = 2 • = Min[6;8]=6
4
10
• mt(O5) = T*5 – T5 = 14 - 12 = 2 0• 2 mL(O4) = T5 - T12- V(4,5)
4
14 =0
o4 12 o5
• mt(O6) = T* 6 – T6 = 10 – 10 = 0 • mL(O5) = T8 - T5 - V(5,8) = 0
• mL(O9) = TF - T9 - V(9,F) = 0
Exercice 1 :
Une importante société de magasins alimentaires à grande surface
diversifie son activité en créant des commerces dans de petites villes.
La société crée un fonds de commerce qui est ensuite géré de façon
autonome par un commerçant franchisé.
Tout d’abord, la société réalise une étude d’implantation : étude de marché
sur un certain rayon d’action et choix de la localité où sera installé le
commerce.
A partir du jour où l’étude d’implantation est terminée, les tâches
suivantes doivent être exécutées.
Travail à faire :
1. Élaborer la matrice des niveaux
2. Représenter cette succession de tâches par un graphe MPM
3. Déterminer la durée minimale pour que le magasin soit ouvert à la
clientèle
4. Indiquer le chemin critique
5. Préciser à quelles dates au plus tard devront commencer les tâches qui
ne font pas partie du chemin critique.
Nature Durée Antériorité des
Tâche
(jours) tâches
A Recherche d’un local 50 -
B Recherche d’un franchisé 45 -
Constitution du dossier bancaire du
C 15 A, B
franchisé
Constitution du dossier à la chambre de
D 10 A, B
commerce pour les inscriptions obligatoires
E Formation du franchisé 30 B
Aménagement, plâtrerie-peinture du
F 20 A
magasin
G Réfection 8 A
H Équipement de la chambre froide 8 A, F
I Equipement des rayonnages 5 A, F
Implantation du magasin (disposition des
J 6 A,B, E, F, G, H, I
articles)
Tirage en imprimerie des feuillets
K 6 A, B, D
publicitaires
L Distribution de feuillets publicitaires 2 A, B, D, K
Liste et envoi des invitations pour
M 6 A, B, D
l’inauguration
N Inauguration du magasin 1 Toutes les autres
1- La matrice des niveaux :
b - b
c - c
d a d
e
e a
f
f b
g d g
h e,c,f h
2- Le graphe MPM (représentation sagittale):
a 3 d 2 g
0 1 3 5 5 7
2
3
e Fin
3 4 9 9
Début 1
0 0
4
b 2 f 3 h
0 0 2 2 5 5
c 4
0 1
a 1 0
b 0 0
c 1 1
d 2 0
e 1 1
f 0 0
g 2 2
h 0 0
Étude de cas en ordonnancement
Interprétations
Construction du graphe :
un arc correspond à une tâche
la valeur de l'arc représente la durée de la tâche
Un sommet correspond à la fin de certaines taches et au début
d’autres. Il représente un événement, une étape d’avancement
du projet:
•
toutes les tâches qui y arrivent sont terminées