Pages de MNM1 - Exos - Doc1
Pages de MNM1 - Exos - Doc1
Pages de MNM1 - Exos - Doc1
Méthodes Numériques
1 Programmation linéaire 1
Méthode du simplexe . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
Raffinerie de pétrole . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
Méthode des variables ajoutées . . . . . . . . . . . . . . . . . . . . . . . . 6
Indices d’octane . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
Fabrique de pièces détachées . . . . . . . . . . . . . . . . . . . . . . . . . 13
Plan de production de moteurs . . . . . . . . . . . . . . . . . . . . . . . . 15
Excavation et matériaux de carrière . . . . . . . . . . . . . . . . . . . . . . 17
2 Dualité 19
Main d’oeuvre et équipements . . . . . . . . . . . . . . . . . . . . . . . . 19
Trois techniques de production . . . . . . . . . . . . . . . . . . . . . . . . 21
Production en heures-machines . . . . . . . . . . . . . . . . . . . . . . . . 22
1 Programmation linéaire
Programme 1
Max (x1 + 2x2 )
x1 + 3x2 ≤ 21
−x1 + 3x2 ≤ 18
x1 − x2 ≤ 5
x et x ≥ 0
1 2
On introduit des variables d’écart, ce qui conduit aux équations suivantes pour les
contraintes du problème :
x1 + 3x2 + x3
= 21
−x1 + 3x2 + x4 = 18
x1 − x2 + x5 =5
1
x1 x2 x3 x4 x5
1 3 1 0 0 21 x3
-1 3 0 1 0 18 x4
1 -1 0 0 1 5 x5
-1 -2 0 0 0 0
2
Ce tableau correspond à l’optimum car il n’y a plus de termes négatifs dans la
dernière ligne. On obtient donc comme solution :
∗
x 1 = 9
x∗2 = 4
x∗3 = 0
x∗4 = 15
x ∗ = 0
5
Programme 2
Min (x1 − 3x2 )
3x1 − 2x2 ≤ 7
−x1 + 4x2 ≤ 9
−2x1 + 3x2 ≤ 6
x et x ≥ 0
1 2
3
x1 x2 x3 x4 x5
5/3 0 1 0 2/3 11 x3
5/3 0 0 1 -4/3 1 x4
-2/3 1 0 0 1/3 2 x2
-1 0 0 0 1 6
Il n’y a plus de terme négatif dans la dernière ligne et on est donc à l’optimum. La
solution est :
x∗1 = 3/5
∗
x2 = 12/5
x∗3 = 10
x∗4 = 0
x ∗ = 0
5
4
Si on note x3 , x4 , x5 les variables d’écart, les contraintes deviennent :
5x1 + 7x2 + x3 = 16500
x1 + x2 + x4 = 2500
9x1 + 7x2 + x5 = 21300
x2 entre et x3 sort.
Tableau 2
x1 x2 x3 x4 x5
5/7 1 1/7 0 0 16500/7 x2
2/7 0 -1/7 1 0 1000/7 x4
4 0 -1 0 1 4800 x5
-1/7 0 4/7 0 0 66000/7
x1 entre et x4 sort.
Tableau 3
x1 x2 x3 x4 x5
0 1 1/2 -5/2 0 2000 x2
1 0 -1/2 7/2 0 500 x1
0 0 1 -14 1 2800 x5
0 0 1/2 1/2 0 9500
Il n’y a plus de terme négatif dans la dernière ligne et on est donc à l’optimum. La
solution est :
x∗1 = 500
∗
x2 = 2000
x∗3 = 0
x∗4 = 0
x∗ = 2800
5
5
Corrigé ex. 3 : Méthode des variables ajoutées
Les deux programmes d’optimisation de cet exercice présentent une difficulté sup-
plémentaire pour appliquer la méthode du simplexe : on ne peut pas démarrer le sim-
plexe à partir de l’origine (c’est-à-dire à partir du point de coordonnées nulles) car ce
point ne vérifie pas les contraintes. L’origine ne fait pas partie du domaine réalisable.
Il faut donc trouver un point de départ dans le domaine réalisable, autrement dit
trouver un point à coordonnées positives qui vérifie les équations des contraintes. On
utilise pour cela la méthode des variables ajoutées. Elle consiste à introduire des va-
riables supplémentaires x1,a , x2,a , . . . dans les contraintes et à chercher à les annuler.
Comme ce sont des variables positives, il suffit d’annuler leur somme et on en fait un
problème d’optimisation en fixant comme objectif de minimiser cette somme :
X
Min xj,a
j
Programme 1
Max (x1 − x2 + x3 )
−3x1 + 2x2 + x3 = 1
x1 − x2 − x3 + x4 = 3
x1 + 4x2 + 2x3 − 2x4 = 1
x1 , x2 , x3 et x4 ≥ 0
On introduit 3 variables positives x1,a , x2,a , x3,a dans les contraintes et on cherche
à minimiser la fonction objectif x1,a + x2,a + x3,a . On se ramène à un problème de
maximisation en changeant le signe de cette fonction objectif. Le problème s’écrit donc
sous la forme suivante
Max (−x1,a − x2,a − x3,a )
−3x1 +2x2 +x3 +x1,a = 1
x1 −x2 −x3 +x4 +x2,a = 3
x1 +4x2 +2x3 −2x4 +x3,a = 1
Très important : il faut veiller à ce que la fonction objectif (−x1,a − x2,a − x3,a )
soit exprimée en fonction des variables hors-base. C’est une règle qui doit toujours être
vérifiée :
6
À tous les stades de la méthode du simplexe, la fonction objectif et les variables
de base doivent être exprimées en fonction des variables hors-base.
On doit donc, avant de commencer, extraire x1,a , x2,a , x3,a en fonction de x1 , x2 , x3 , x4
et les remplacer dans la fonction objectif. On a :
x1,a = 1 + 3x1 − 2x2 − x3
x2,a = 3 − x1 + x2 + x3 − x4
x3,a = 1 − x1 − 4x2 − 2x3 + 2x4
D’où
−x1,a − x2,a − x3,a = −5 − x1 + 5x2 + 2x3 − x4
À partir de là, la méthode du simplexe s’applique sans problèmes.
Tableau 1
x1 x2 x3 x4 x1,a x1,a x3,a
-3 2 1 0 1 0 0 1 x1,a
1 -1 -1 1 0 1 0 3 x2,a
1 4 2 -2 0 0 1 1 x3,a
1 -5 -2 1 0 0 0 -5
x2 entre et x3,a sort.
Tableau 2
x1 x2 x3 x4 x1,a x1,a x3,a
-7/2 0 0 1 1 0 -1/2 1/2 x1,a
5/4 0 -1/2 1/2 0 1 1/4 13/4 x2,a
1/4 1 1/2 -1/2 0 0 1/4 1/4 x2
9/4 0 1/2 -3/2 0 0 5/4 -15/4
x4 entre et x1,a sort.
Tableau 3
x1 x2 x3 x4 x1,a x1,a x3,a
-7/2 0 0 1 1 0 -1/2 1/2 x4
3 0 -1/2 0 -1/2 1 1/2 3 x2,a
-3/2 1 1/2 0 1/2 0 0 1/2 x2
-3 0 1/2 0 3/2 0 1/2 -3
x1 entre et x2,a sort.
Tableau 4
x1 x2 x3 x4 x1,a x1,a x3,a
0 0 -7/12 1 5/12 7/6 1/12 4 x4
1 0 -1/6 0 -1/6 1/3 1/6 1 x1
0 1 1/4 0 1/4 1/2 1/4 2 x2
0 0 0 0 1 1 1 0
Dans le dernier tableau, les trois variables ajoutées sont sorties de la base. Elles
sont donc nulles, ce qui était l’objectif. Cela signifie que les variables qui sont main-
tenant dans la base constituent une solution à coordonnées positives pour le système
7
des contraintes. On a donc trouvé un point de départ pour résoudre le problème de
l’exercice. C’est le point de coordonnées (ici x3 est nulle car elle est hors-base ) :
x1 = 1
x = 2
2
x 3 =0
x4 = 4
On peut donc maintenant traiter le problème posé à partir du point trouvé. On com-
mence par supprimer, dans le dernier tableau calculé, les colonnes des variables ajou-
tées :
x1 x2 x3 x4
0 0 -7/12 1 4 x4
1 0 -1/6 0 1 x1
0 1 1/4 0 2 x2
Dans ce tableau, on voit que les variables x1 , x2 et x4 sont dans la base et que
la variable x3 est hors-base. On peut l’interpréter comme le système de contraintes
suivant :
− 7 x + x4 = 4
12 3
x1 − 16 x3 = 1
x2 + 14 x3 = 2
8
Programme 2
Max (x1 + 2x2 + 3x3 )
x1 + x2 ≤ 5
2x1 + 2x2 − x3 = 6
12x1 + 8x2 − 5x3 = 32
x1 , x2 et x3 ≥ 0
Dans ce problème, la première contrainte est une inégalité, donc il faut commencer par
introduire une variable d’écart x4 . On introduit ensuite les variables ajoutées comme
dans l’exercice précédent. Le problème s’écrit sous la forme
Très important : il faut veiller à ce que la fonction objectif (−x1,a − x2,a − x3,a )
soit exprimée en fonction des variables hors-base. On doit donc, avant de commencer,
extraire x1,a , x2,a , x3,a en fonction de x1 , x2 , x3 , x4 et les remplacer dans la fonction
objectif. On a :
x1,a = 5 − x1 − x2 − x4
x2,a = 6 − 2x1 − 2x2 + x3
x3,a = 32 − 12x1 − 8x2 + 5x3
D’où
−x1,a − x2,a − x3,a = −43 + 15x1 + 11x2 − 6x3 + x4
À partir de là, la méthode du simplexe s’applique sans problèmes :
Tableau 1
x1 x2 x3 x4 x1,a x1,a x3,a
1 1 0 1 1 0 0 5 x1,a
2 2 -1 0 0 1 0 6 x2,a
12 8 -5 0 0 0 1 32 x3,a
-15 -11 6 -1 0 0 0 -43
x1 entre et x3,a sort.
Tableau 2
9
x1 x2 x3 x4 x1,a x1,a x3,a
0 1/3 5/12 1 1 0 -1/12 7/3 x1,a
0 2/3 -1/6 0 0 1 -1/6 2/3 x2,a
1 2/3 -5/12 0 0 0 1/12 8/3 x1
0 -1 -1/4 -1 0 0 5/4 -3
x2 entre et x2,a sort.
Tableau 3
x1 x2 x3 x4 x1,a x1,a x3,a
0 0 1/2 1 1 -1/2 0 2 x1,a
0 1 -1/4 0 0 3/2 -1/4 1 x2
1 0 -1/4 0 0 -1 1/4 2 x1
0 0 -1/2 -1 0 3/2 1 -2
x4 entre et x1,a sort.
Tableau 4
x1 x2 x3 x4 x1,a x1,a x3,a
0 0 1/2 1 1 -1/2 0 2 x4
0 1 -1/4 0 0 3/2 -1/4 1 x2
1 0 -1/4 0 0 -1 1/4 2 x1
0 0 0 0 1 1 1 0
Dans le dernier tableau, les trois variables ajoutées sont sorties de la base. Elles
sont donc nulles, ce qui était l’objectif. Cela signifie que les variables qui sont main-
tenant dans la base constituent une solution à coordonnées positives pour le système
des contraintes. On a donc trouvé un point de départ pour résoudre le problème de
l’exercice. C’est le point de coordonnées (ici x3 est nulle car elle est hors-base ) :
x1 = 1
x = 2
2
x 3 =0
x4 = 2
On peut donc maintenant traiter le problème posé à partir du point trouvé. On com-
mence par supprimer, dans le dernier tableau calculé, les colonnes des variables ajou-
tées :
x1 x2 x3 x4
0 0 1/2 1 2 x4
0 1 -1/4 0 1 x2
1 0 -1/4 0 2 x1
Dans ce tableau, on voit que les variables x1 , x2 et x4 sont dans la base et que
la variable x3 est hors-base. On peut l’interpréter comme le système de contraintes
suivant :
1
x3 + x4 = 2
2
1
x2 − x3 = 1
4
1
x1 − x3 = 2
4
10
La dernière ligne doit contenir la fonction objectif initiale x1 + 2x2 + 3x3 mais
celle-ci doit être exprimée, comme toujours, en fonction de la ou des variable(s) hors-
base uniquement. Le système précédent permet facilement de tout exprimer en fonction
de x3 qui est ici l’unique variable hors-base. On trouve :
15
x1 + 2x2 + 3x3 = 4 + x3
4
Le tableau du simplexe s’écrit donc
x1 x2 x3 x4
0 0 1/2 1 2 x4
0 1 -1/4 0 1 x2
1 0 -1/4 0 2 x1
0 0 -15/4 0 4
La variable x3 entre et x4 sort. Par pivot, on obtient le tableau suivant :
x1 x2 x3 x4
0 0 1 2 4 x3
0 1 0 1/2 2 x2
1 0 0 1/2 3 x1
0 0 0 15/2 19
On est maintenant à l’optimum et la solution du problème est :
∗
x1 = 3
x ∗ = 2
2
∗
x 3 =4
∗
x4 = 0
11
Les containtes de disponibilité des ressources P1 et P2 s’écrivent comme ceci :
(
x1A + x1B ≤ 3900
x2A − x2B ≤ 5000
Notons x01 , x02 , x03 , x04 les variables d’écart. Le tableau de départ de la méthode du
simplexe s’écrit :
Tableau 1
x1A x2A x1B x2B x01 x02 x03 x04
25 -3 0 0 1 0 0 0 0 x01
0 0 1 -1 0 1 0 0 0 x02
1 0 1 0 0 0 1 0 3900 x03
0 1 0 1 0 0 0 1 5000 x04
-2,5 -1,5 -1,5 -0,5 0 0 0 0 16125
12
x2B entre et x03 sort.
Tableau 5
x1A x2A x1B x2B x01 x02 x03 x04
1 0 0 0 1/22 3/22 -3/22 3/22 150 x1A
0 0 1 0 -1/22 -3/22 25/22 -3/22 3750 x1B
0 0 0 1 -1/22 -25/22 25/22 -3/22 3750 x2B
0 1 0 0 1/22 25/22 -25/22 25/22 1250 x2A
0 0 0 0 1/11 14/11 2,5/11 19,5/11 25875
La marge sur coût variable unitaire réalisée pour les lots de type A, compte-tenu du
nombre d’unités d’oeuvre requis et de leur coût de fabrication, est :
13
La marge totale est c1 x1 + c2 x2 . On cherche donc à maximiser la marge :
14
Corrigé ex. 6 : Plan de production de moteurs
Les prix de vente sont fixés à 215 e pour le modèle A et 150 e pour le modèle B.
On désigne par x1 et x2 les quantités de moteurs des deux types A et B qui vont
être produites.
Calculons les marges bénéficiaires résultant de ces fabrications. Pour le modèle
A, le prix de vente unitaire est de 215 et on doit retirer les coûts de fabrication qui
dépendent du temps passé dans les trois ateliers (attention les temps sont en minutes et
les coûts sont exprimés à l’heure). On trouve donc :
Max(50 x1 + 25 x2 )
Il y a d’autre part une contrainte de marché qui impose un quota maximal sur le
nombres de moteurs de type A :
x1 ≤ 1800
On introduit des variables d’écart x3 , x4 , x5 , x6 dans les quatre contraintes :
50x1 +40x2 +x3 = 2500 × 60
30x1 +20x2 +x4 = 1000 × 60
20x 1 +10x 2 +x 5 = 800 × 60
x1 +x6 = 1800
15
x1 x2 x3 x4 x5 x6
60 40 1 0 0 0 150000 x3
30 20 0 1 0 0 60000 x4
20 10 0 0 1 0 48000 x5
1 0 0 0 0 1 1800 x6
-50 -25 0 0 0 0 0
Tableau 2
x1 x2 x3 x4 x5 x6
0 40 1 0 0 -60 42000 x3
0 20 0 1 0 -30 6000 x4
0 10 0 0 1 -20 12000 x5
1 0 0 0 0 1 1800 x1
0 -25 0 0 0 50 90000
Tableau 3
x1 x2 x3 x4 x5 x6
0 0 1 -2 0 0 30000 x3
0 1 0 0,05 0 -1,5 300 x2
0 0 0 -0,5 1 -5 9000 x5
1 0 0 0 0 1 1800 x1
0 0 0 1,25 0 12,5 97500
La deuxième et la quatrième contrainte sont saturées : les valeurs 1,25 et 12,5 dans
la dernière ligne du tableau sont les prix duaux π4 et π6 associés. On fabrique le maxi-
mum envisagé de moteurs de type A et l’atelier de soudure fonctionne à plein. Le
premier atelier (emboutissage) est sous-utilisé : il reste 500 (=30000/60) heures dispo-
nibles. De même, dans le troisième atelier, il reste 9000/60=150 heures disponibles.
16
Corrigé ex. 7 : Excavation et matériaux de carrière
On désigne par x1 et x2 les quantités qui seront extraites des deux carrières. La
redevance à acquitter est de 19, 40 × 103 x1 + 20 × 103 x2 et on cherche à la minimiser.
Pour simplifier les calculs, on divise les coefficients par 103 , d’où le problème :
Les rendements liés au concassage des matériaux conduisent aux inéquations sui-
vantes qui expriment que les quantités obtenues doivent pouvoir couvrir les besoins
imposés par le contrat :
3 3
0, 36 × 10 x1 + 0, 45 × 10 x2 ≥ 13500
0, 40 × 103 x1 + 0, 20 × 103 x2 ≥ 11200
0, 16 × 103 x1 + 0, 10 × 103 x2 ≥ 5000
Il existe des méthodes pour trouver un point de démarrage pour le simplexe : cela
revient, une fois qu’on a introduit les variables d’écart, à résoudre un système d’équa-
tions linéaires en coordonnées positives. Une des méthodes possibles est la méthode
des valeurs ajoutées (mais on ne l’utilisera pas ici).
Dans le cas particulier de cet exercice, on peut se contenter plus simplement de
déterminer un point valide sur l’un des axes. Par exemple, si on regarde les intersec-
tions des contraintes avec l’axe vertical (x1 = 0), on trouve les valeurs 30, 50, 56.
Comme le domaine se trouve au-dessus de ces points, on choisit le point de coordon-
nées (0, 56) comme point de départ. On peut vérifier qu’il satisfait effectivement les
trois contraintes : il est sur la droite de la deuxième contrainte, donc x4 = 0.
On part donc de la situation suivante :
x1 = 0 x4 = 0
x2 = 56 x3 = 130 x5 = 30
Les variables x1 et x4 sont des variables hors-base. On doit donc exprimer les autres
variables en fonction de celles-ci. On obtient :
x2 = 56 − 2x1 + x4
x3 = 130 − 6x1 + 5x4
x5 = 30 − 2x1 + 5x4
17
et donc
2x1 +x2 −x4 = 56
6x1 +x3 −5x4 = 130
2x1 −5x4 +x5 = 30
18
2
S
1
isoquante
O 1
2 Dualité
Max {3/2 x1 + x2 }
2x + x ≤ 2
1 2
x 2 ≤ 1
x1 et x2 ≥ 0
8-1)
En introduisant des variables d’écart, les contraintes s’écrivent :
2x1 + x2 +x3 =2
x2 +x4 =1
avec
c = (3/2 1 0 0)
2 1 1 0
B=
0 1 0 1
19
Dans cette notation, le vecteur des prix duaux s’écrit
et par conséquent
−1
2 1 1 1 −1
π(I) = (3/2 1) = (3/2 1) = (3/4 1/4)
0 1 2 0 2
f ∗ = π 1 b1 + π 2 b2
∆f = π1 ∆b1 + π2 ∆b2
Cette condition est valable quelle que soit l’échelle des extensions, car tant que le
rapport entre b1 et b2 reste inchangé, on maintient la relation b2 ≤ b1 trouvée à la
question 2 qui est la condition pour que les prix duaux trouvés restent valides.
20
Corrigé ex. 9 : Trois techniques de production
9-1)
En appelant x1 , x2 et x3 les quantités de bien fabriquées selon l’une des trois tech-
niques, le programme d’optimisation s’écrit de la manière suivante :
Max {3x1 + 4x2 + 5x3 }
0, 5 x + 1, 5 x + 2 x ≤ 12
1 2 3
2 x1 + 1, 5 x2 + 0, 5 x3 ≤ 15
x1 , x 2 , x 3 ≥ 0
9-2)
On introduit des variables d’écart x4 et x5 :
0, 5 x1 +1, 5 x2 +2 x3 +x4 = 12
2 x1 +1, 5 x2 +0, 5 x3 +x5 = 15
21
C’est l’optimum :
x∗1 = 6, 4 = 96/15
∗
x 2 =0
x∗3 = 4, 4 = 66/15
x∗4
=0
x ∗
=0
5
10-1) Calcul des marges sur coût variable unitaires pour les deux produits P1 et
P2 .
On évalue les données en fonction des résultats du dernier exercice mensuel. La
fabrication du produit P1 a requis 650 heures-machines dans l’atelier d’usinage alors
que celle du produit P2 en a requis 350. Le produit P1 utilise 650/(650+350)=65% du
22
temps de l’atelier d’usinage. Les coûts variables dans cet atelier ont été de 80000, donc
le coût correspondant au produit P1 est de 80000 × 0.65. Comme il a été fabriqué 5000
pièces de P1 , on a finalement un coût unitaire de
80000
× 0.65 = 10.4.
5000
De la même manière, on calcule le coût unitaire pour le produit P2 dans l’atelier d’usi-
nage comme ceci :
80000
× 0.35 = 4.
7000
On calcule de même les coûts unitaires dans l’atelier de finition. La fabrication du
produit P1 a requis 150 heures-machines dans l’atelier de finition alors que celle du
produit P2 en a requis 350. Le produit P1 utilise 150/(150+350)=30% du temps de
l’atelier de finition. Les coûts variables dans cet atelier ont été de 50000, donc le coût
correspondant au produit P1 est de 50000 × 0.30. Comme il a été fabriqué 5000 pièces
de P1 , on a finalement un coût unitaire de
50000
× 0.30 = 3.
5000
De la même manière, on calcule le coût unitaire pour le produit P2 dans l’atelier de
finition comme ceci :
50000
× 0.70 = 5.
7000
Compte-tenu des autres coûts imputés directement aux produits, la marge sur coût
variable du produit P1 est finalement :
m1 = 50 − 12 − 10.4 − 3 = 24.6
m2 = 30 − 5.8 − 4 − 5 = 15.2
f = 24.6x1 + 15.2x2
23
Le programme d’optimisation est finalement :
Max (24.6 x1 + 15.2 x2 )
0.13 x1 + 0.05 x2 ≤ 1100
0.03 x1 + 0.05 x2 ≤ 550
On trouve (
x∗1 = 5500
x∗2 = 7700
La valeur de la fonction objectif en ce point optimal est f = 252340.
10-2) Les conditions d’optimalité s’écrivent :
10-3) La formule exprimant la marge totale en fonction des prix duaux est π1 b1 +
π2 b2 . On trouve ici :
24
pour un nombre d’heures-machines de fonctionnement de 650 + 350 = 1000. Le prix
d’usage sera donc de r1 = 85.
Pour l’atelier de finition, on trouve de même 40000 + 25000 = 65000 pour un
nombre d’heures-machines de fonctionnement de 150 + 350 = 500. Le prix d’usage
sera donc de r2 = 65000/500 = 130.
0-2 ) En comparant les prix d’usage avec les prix duaux, on se rend compte qu’une
augmentation de capacité dans l’atelier d’usinage ou dans l’atelier de finition sera pro-
fitable car r1 < π1 et r2 < π2 .
25