Chap 3 Dualité
Chap 3 Dualité
Chap 3 Dualité
Primal Dual
Max Z = 8 x1 + 6 x2
sujet : Min W = 3000 y1 + 2400 y2 + 1800 y3
5 x1 + 3 x2 3000 sujet :
5 y1 + 2 y2 + y3 8
2 x1 + 3 x2 2400 3 y + 3 y + 3 y 6
x1 + 3 x2 1800 1 2 3
y1 , y2 , y3 0
x1 , x2 0
Avec :
X = ( x j ) : vecteur des quantits fabriques.
(
= CT Y T A X . )
Si > 0 , on gagne plus en travaillant, donc il doit produire.
Si < 0 , on gagne plus en revendant, donc il doit revendre et ne pas produire.
Exemple :
Un industriel possde, dans lun de ses ateliers de production, trois machines
M1, M2 et M3 qui sont partiellement inutilises :
- M1 est libre 18 heures par semaine.
- M2 est libre 8 heures par semaine.
- M3 est libre 14 heures par semaine.
Pour mieux rentabiliser ces machines on envisage de fabriquer deux produits
P1 et P2 dont lcoulement sur le march ne poserait pas de problme.
- Chaque unit de P1 rapporterait 1 dinars et consommerait : 1 heure de M1,
1 heure de M2 et 2 heures de M3.
- Chaque unit de P2 rapporterait 2 dinar et consommerait : 3 heures de M1,
1 heure de M2 et 1 heure de M3.
Le (PL) permettant de maximiser ses bnfices est :
Max Z = x1 + 2 x2
sujet :
x1 + 3 x2 18
.
x1 + x2 8
2 x1 + x2 14
x1 , x2 0
1 1 1
On lui propose de louer ses trois machines aux prix Y T = , , lheure de
2 4 2
fonctionnement. Est-ce que la proposition intresse le producteur ou non ?
( *** )
Calculons ?
( )
= CT Y T A X .
1 3
1 1 1 x
C = (1 2 ) ; A = 1 1 ; Y T = , , et X = 1
T
2 1 2 4 2 x2
1 3
1 1 1 x1 7 9 x1
= (1 2 ) 1 1 x = (1 2 )
2 4 2 2 4 4 x2
2 1
3 1 x
= 1 < 0
4 4 x2
Quelles sont les quantits des bruts B1 et B2 quon doit utiliser pour
minimiser les dpenses.
Soient :
x1 : quantit de B1 utilise.
x2 : quantit de B2 utilise.
La dpense minimiser est : W = 20 x1 + 25 x2 .
Les contraintes sont :
0,1x1 + 0,5 x2 0,5 x1 + 5 x2 5
0,1x + 0, 2 x 0, 4 x + 2x 4
1
1 2
2
0,3 x1 + 0, 2 x2 0, 6 3 x1 + 2 x2 6
x1 , x2 0 x1 , x2 0
Avec :
X = ( x j ) : vecteur des quantits de bruts de produit de base consomm.
= ( C T Y T A) X .
m
Si yi est libre alors aij xi = c j
i =1
Exemple1 :
Primal Dual
Max Z = x1 + 2 x2
sujet : Min W = 18 y1 + 8 y2 + 14 y3
x1 + 3 x2 18 sujet :
y1 + y2 + 2 y3 1
x1 + x2 8 3 y + y + y 2
2 x1 + x2 14 1 2 3
y1 , y2 , y3 0
x1 , x2 0
Exemple2 :
Primal Dual
Min Z = 2 x1 + 2,5 x2
sujet : Max W = 5 y1 + 4 y2 + 6 y3
x1 + 5 x2 5 sujet :
y1 + y2 + 3 y3 2
x1 + 2 x2 4 5 y + 2 y + 2 y 2,5
3 x1 + 2 x2 6 1 2 3
y1 , y2 , y3 0
x1 , x2 0
Exemple3 : (Forme non standard)
Min Z = x1 + x2 + x3
sujet :
x1 3 x2 + 4 x3 = 5
Primal : Dual ?
x1 2 x2 3
2 x2 x3 4
x1 , x2 0 ; x3 est libre
2 x2 x3 4
x1 , x2 0 ; x3 est libre
2 x2 x3+ + x3 4
+
x1 , x2 , x3 , x3 0
Primal Dual
? X 0 ? Y 0
T
AX B A Y C Y T A CT
C T X = Z max BT Y = W min
Si on multiplie les contraintes du primal par le vecteur Y T et celles du dual
par le vecteur X , on peut crire :
Y T AX Y T B et Y T AX C T X
Ceci implique que : C T X Y T AX Y T B CT X Y T B .
loptimal : C T X = Y T B Z max = W min
Ceci relve une proprit importante des couples de problmes duals :
Si le primal ou le dual a une solution optimale finie, lautre problme a aussi
une solution optimale finie et les valeurs optimales de leurs fonctions
conomiques sont gales.
3.6 Thorme de complmentarit des carts
Soit le problme primal suivant :
a11 x1 + a12 x2 + ... + a1n xn xn +1 = b1
a x + a x + ... + a2 n xn xn + 2 = b2
21 1 22 2
am1 x1 + am 2 x2 + ... + amn xn xn + m = bm
x j 0
c1 x1 + c2 x2 + ... + cn xn = W min
Le problme dual correspondant est :
a11 y1 + a21 y2 + ... + am1 ym + ym +1 = c1
a y + a y + ... + am 2 ym + ym + 2 = c2
12 1 22 2
a1n y1 + am 2 y2 + ... + amn ym + yn + m = cn
y j 0
b1 y1 + b2 y2 + ... + bm ym = Z max
Si on multiplie chaque quation i, du primal, par yi et soustraire les
diffrentes lignes multiplies par yi de la ligne de la fonction objectif on obtient :
m m m m m
c1 x1 ai1 x1 yi + c2 x2 ai 2 x2 yi + + cn xn aim xn yi + xn + i yi = W min bi yi
i =1 i =1 i =1 i =1 i =1
m
m
m
m m
c1 ai1 yi x1 + c2 ai 2 yi x2 + + cn ain yi xn + xn +i yi = W min bi yi
i =1 i =1 i =1 i =1 i =1
ym+1 Z max
m
ym +1 x1 + ym + 2 x2 + + ym + n xn + xn +i yi = W min Z max
i =1
n m
loptimal W Z = 0 ym+ j x j + xn+i yi = 0
j =1 i =1
xn +i yi = 0
x j ym + j = 0
y1 , y2 , y3 0
x1 , x2 0
Si on les crit sous la forme standard on obtient :
Primal Dual
Min Z = 2 x1 + 3x2
sujet : Max W = 3 y1 + 5 y2 + 4 y3
2 x1 + x2 x3 = 3 sujet :
2 y1 + 2 y2 + y3 + y4 = 2
2 x1 x2 x4 = 5 y y + 4y + y = 3
x1 + 4 x2 x5 = 4 1 2 3 5
yi 0, i = 1,...,5
x j 0, j = 1,...,5
x = 3
La solution 1 est elle optimale ?
x2 = 1
La complmentarit des carts permet dcrire :
y1 ( 2 x1 + x2 3) = 0
0 y1 = 0
y2 ( 2 x1 x2 5 ) = 0 y2 0
=0 y = 0
3
y
3 1( x + 4 x2 4 ) = 0
0
x1 ( 2 2 y1 2 y2 y3 ) = 0 2 0 2 y2 0 = 0
et
x2 ( 3 y1 + y2 4 y3 ) = 0 3 0 + y2 0 = 0
y =1
Do : 2
y2 = 3
La solution est non optimale.
26
x1 = 9
Essayons maintenant la solution ?
x = 7
2 9
La complmentarit des carts permet dcrire :
y1 ( 2 x1 + x2 3) = 0
0 y1 = 0
y2 ( 2 x1 x2 5 ) = 0 y2 0
=0 y 0
3
y3 ( x1 + 4 x2 4 ) = 0
0
x1 ( 2 2 y1 2 y2 y3 ) = 0
et
x2 ( 3 y1 + y2 4 y3 ) = 0
5
2 0 2 y2 y3 = 0 y2 = 9
3 0 + y2 4 y3 = 0 y = 8
3 9
La solution est optimale.
Cette proprit peut tre nonce pour nimporte quelle paire de problmes
duals. Elle est trs importante car elle permet dobtenir la solution optimale dun
des programmes partir du tableau final de lautre programme.
Une fois crit sous forme standard cette paire de programmes duals possde
le mme nombre de variables dans le primal et dans le dual. chaque contrainte
du primal sont associes une variable dcart primale et une variable de dcision
duale. De mme, chaque contrainte du dual sont associes une variable dcart
duale et une variable de dcision primale. Il existe donc une correspondance
biunivoque entre les variables de dcisions primales et les variables dcart
duales dune part, et entre les variables de dcisions duales et les variables
dcart primales dautre part.
Dans cet exemple on a les correspondances suivantes :
x1 y4 ; x2 y5 et y1 x3 ; y2 x4 ; y3 x5 .
Nous avons dit que nous pouvons obtenir la solution optimale de lun des
programmes partir du tableau final de lautre programme. Examinons ceci pour
la paire de programmes duals suivante :
Primal Dual
Max Z = 2 x1 + 3 x2
sujet : Min W = 18 y1 + 8 y2 + 14 y3
x1 + 3 x2 18 sujet :
y1 + y2 + 2 y3 2
x1 + x2 8 3 y + y + y 3
2 x1 + x2 14 1 2 3
y1 , y2 , y3 0
x1 , x2 0
Si on les crit sous la forme standard on obtient :
Primal Dual
Max Z = 2 x1 + 3x2
sujet : Min W = 18 y1 + 8 y2 + 14 y3
x1 + 3 x2 + x3 = 18 sujet :
y1 + y2 + 2 y3 y4 = 2
x1 + x2 + x4 = 8 3 y + y + y y = 3
2 x1 + x2 + x5 = 14 1 2 3 5
yi 0, i = 1,...,5
x j 0, j = 1,...,5
Le tableau final du primal est :
x1 x2 x3 x4 x5 b
x2 0 1 1/2 -1/2 0 5
x1 1 0 -1/2 3/2 0 3 Solution optimale
x5 0 0 1/2 -5/2 1 3
Z 0 0 1/2 3/2 0 21
x1 y4 ; x2 y5 et y1 x3 ; y2 x4 ; y3 x5 .
x1 x2 x3 x4 x5 x6 b
x4 -1/4 -1/2 -1/4 1 0 0 -3
x5 -2/5 -3/10 -3/10 0 1 0 -2 ngatifs
x6 -2/5 -2/5 -1/5 0 0 1 -2,5
-Z 15 20 12 0 0 0 0
non ngatifs
Dans ce tableau les cots rduits de dpart sont non ngatifs et le cot droit
est ngatif, ce qui implique que la solution de base duale correspondante est
ralisable et que la solution primale nest pas ralisable. Par consquent il
continuer la rsolution jusqu ce que le cot droit devient non ngatif, c'est--dire
jusqu ce que la solution primale devienne ralisable, ce qui revient rsoudre le
problme dual.
2me tape :(choix de la variable sortante)
Identifier la variable de base qui est la plus ngative et la choisir comme
variable sortante de la base.
3me tape :(recherche de limitation)
Si la ligne s est celle de la variable sortante, choisir comme nouvelle variable
de base celle qui correspond :
cj c
xr = min asj < 0 = r
j a asr
sj
Si tous les asj* ne sont pas ngatifs, alors le programme dual est sans borne.
Dans cet exemple la variable la plus ngative est x4 (variable sortante). La
variable rentrante, qui correspond la plus petite valeur de la limitation
20
( = 40 ), est x2.
( 1 2)
4me tape : (transformation du tableau)
Changer de base, rsoudre pour trouver la valeur des nouvelles variables de
base en transformant le tableau comme dans la mthode primale du simplexe.
Retourner la deuxime tape.
Ce qui donne le tableau suivant :
x1 x2 x3 x4 x5 x6 b
x2 1/2 1 1/2 -2 0 0 6
x5 -1/4 0 -3/20 -3/5 1 0 -1/5 t
x6 -1/5 0 0 -4/5 0 1 -1/10
-Z 5 0 2 40 0 0 -120
Limitation 20 - 40/3 200/3 - -
u
x1 x2 x3 x4 x5 x6 b
x2 0 1 0 -8/3 10/3 -5/3 5,5
x3 0 0 1 -8/3 -20/3 25/3 1/2 Solution optimale
x1 1 0 0 4 0 -5 1/2
-Z 0 0 0 76/3 40/3 25/3 -123,5