PL - 14-15 - Partiel + Correction

Télécharger au format pdf ou txt
Télécharger au format pdf ou txt
Vous êtes sur la page 1sur 4

FACULTE DES SCIENCES DE TUNIS 04/11 /2014

Dépt. des Sciences de l’Informatique

DEVOIR SURVEILLE

Section : INFO3
Matière : Programmation Linéaire
Durée : 1H30mns
Documents Non Autorisés
------------------------------------------------------------------------------------------------------------------------------------
Le barème est à titre indicatif

Exercice 1 : (5 pts)

Considérons le programme linéaire suivant :

AX = b
(P) X ≥0
CX = Z(Max)

Où A est une matrice ayant m lignes et n colonnes.

Montrer que si X1, X2, . . ., Xk sont solutions de (P) (Xi ∈ Rn ) alors

X = α1X1 + α2X2 + . . . +αkXk où αi ∈ [0 ; 1] et Σαi = 1

est solution de (P).

Problème : (15 pts) (Toutes les réponses doivent être justifiées)

On considère le programme linéaire suivant :

-X1 + 4X2 ≤ 12
(P) X1 + X2 ≤ 8 Xi ≥ 0
-X1 + 3X2 ≥ 0
2X1 + X2 = Z(Max)

1) Résoudre (P) par l’algorithme du simplexe.

2) Déduire à partir du dernier tableau du simplexe, l’inverse de la matrice : A = ቀ1 1



1 −3
Et exprimer A-1 sous forme de produit de matrices de pivotage.

3) Ecrire le dual (D) de (P) et déduire du dernier tableau du simplexe, une solution optimale de (D).
4) Ecrire (D) sous forme standard. Soit (DS) le programme linéaire obtenu. Donner alors la base de (DS) relative à la
solution trouvée à la question 3).

Bon travail
F. Ben Charrada
1
FACULTE DES SCIENCES DE TUNIS Tunis le, 04/11/2014
Dépt. des Sciences de l’Informatique

Devoir Surveillé

Correction

Exercice :
Réalisabilité :

AX = A(α1X1 + α2X2 + . . . +αkXk) = ΣαiAXi = Σαib = b*(Σαi) = b


X est donc solution réalisable de (P).

Optimalité :

Posons Zo, la valeur de Z à l’optimum. On a : CXi = Zo ( i = 1,…,k)


D’autre part CX = C*Σ αiXi = ΣαiCXi = ΣαiZo = Zo* Σαi = Zo
X est donc solution optimale de (P).

Problème:
1) Résolution de (P) par l’algorithme du simplexe :

En multipliant la 3ème contrainte par (-1) et en ajoutant des variables d’écart, (P) s’écrit sous forme standard :

-X1 + 4X2 + X3 = 12
(PS) X1 + X2 + X4 = 8 Xi ≥ 0
X1 - 3X2 + X5 = 0
2X1 + X2 = Z(Max)

2
X1 X2 X3 X4 X5 b

-1 4 1 0 0 12

1 1 0 1 0 8

1 -3 0 0 1 0

2* 1 0 0 0 0 =-Z

0 1 1 0 1 12

0 4 0 1 -1 8

1 -3 0 0 1 0

0 7* 0 0 -2 0 =-Z

0 0 1 -1/4 5/4 10

0 1 0 1/4 -1/4 2

1 0 0 3/4 1/4 6

0 0 0 -7/4 -1/4 -14 =-Z

Jo = {1, 2, 3} est une base optimale. La solution de base associée est :

X1 = 6 ; X2 = 2 ; X3 = 10 et X4 = X5 = 0 ; Z = 14

2) La matrice A n’est autre que la matrice des coefficients de X1 et X2 dans la 2ème et la 3ème ligne du 1er
tableau du simplexe (cf. tableau du simplex). En ignorant la 1ère ligne, les deux matrices de pivotage
associées à la 1ère et 2ème itération sont respectivement :

1 −1 1/4 0
ቁ ൰
0 1 3/4 1
P1=ቀ et P2=൬

Or au tableau I, la matrice unité correspond aux coefficients de X4 et X5. On est passé donc du tableau I au tableau III, en
multipliant le 1er tableau par la matrice B correspondant aux coefficients de X4 et X5 dans le tableau III. C’est-à-dire que :

1/4 −1/4
൰ = P2*P1
3/4 1/4
B =൬

0 1
Par conséquent A*B = ቀ ቁ matrice des coefficients de X1 et X2 dans le tableau III. (cf. tableau du simplexe).
1 0
D’après la propriété des produits de matrices, A-1 est obtenue à partir de B en inversant ses lignes et par
conséquent on a : A-1= P’2*P1 où P’2 est obtenue en inversant les lignes de P2. C’est-à-dire :

3
3/4 1/4 3/4 1 1 −1
A-1= ൬ ൰= ൬ ൰*ቀ ቁ
1/4 −1/4 1/4 0 0 1

3) Le dual (D) de (P) est :

-Y1 + Y2 - Y3 ≥ 2
(D) 4Y1 + Y2 + 3Y3 ≥ 1 Y1, Y2 ≥ 0 et Y3 ≤ 0
12Y1 + 8Y2 = W(min)

D’après le cours, le vecteur multiplicateur π est solution de (D). On a :


1
c’3 = 0 =c3 - πA = 0 - (π1, π2, π3)*൭0൱ = - π1
3
π1 = 0
0
0
c’4 = -7/4 = c4 - πA = 0 - (π1, π2, π3)*൭1൱ = - π2
4
π2 = 7/4
0
0
c’5 = -1/4 = c5 - πA5 = 0 - (π1, π2, π3)*൭0൱ = - π3 π3 = 1/4
1
Or, pour pouvoir appliquer le simplexe, on a multiplié la 3ème contrainte de (P) par (-1) pour avoir une
inégalité de la forme « ≤ ». La solution de (D) est donc : Y = (0, 7/4, -1/4)

4) (D) s’écrit sous forme standard :

-Y1 + Y2 - Y3 - Y4 = 2
(DS 4Y1 + Y2 + 3Y3 - Y5 = 1 Y1, Y2 , Y4 , Y5 ≥ 0 et Y3 ≤ 0
12Y1 + 8Y2 = W(min)

La base contient 2 éléments, et puisque Y2 et Y3 sont non nuls alors JO = {2,3}.

Vous aimerez peut-être aussi