📚 الرئيسية

مسابقة دكتوراه 2014Université des Sciences et de la Technologie Houari Boumediène (USTHB) — الموضوع 06

مسابقة تخصص · Recherche Opérationnelle

Concours d'accès à la Formation Doctorale « Recherche Opérationnelle et Mathématiques Discrètes, ROMaD », Épreuve de Recherche Opérationnelle (Durée : 2h), USTHB, Faculté de Mathématiques — 20 octobre 2014.

التمرين 1

Exercice 1 — Ordonnancement de projet à deux phases (flow-shop)

#scheduling#flow-shop#johnson-algorithm#gantt-diagram

Deux étudiants disposent de 10 jours pour réaliser un travail de groupe. Ce travail se compose de 4 tâches indépendantes entre elles (l'ordre n'a pas d'importance). Chacune de ces tâches peut être divisée en une phase d'analyse, réalisée par le premier étudiant, et une phase de calculs, réalisée par le second. L'analyse doit précéder les calculs. Les temps nécessaires (en jours) à la réalisation de ces tâches sont donnés dans le tableau ci-dessous.

TâchePhase d'analyse (étudiant 1)Phase de calculs (étudiant 2)
A22
B0,81,5
C1,50,5
D21

a. Quel est l'ordonnancement des tâches qui minimise le temps de réalisation du travail ? Illustrez à l'aide d'un diagramme.

Les deux étudiants aimeraient ajouter trois (03) autres parties : E, F et G. Les deux premières (E et F) ne nécessitent pas d'analyse préalable tandis que la dernière (G) comporte seulement une phase d'analyse. Les temps sont :

TâcheAnalyse (étudiant 1)Calculs (étudiant 2)
E-2
F-1,5
G1,5-

b. Avec ces données supplémentaires, pour chacun des étudiants, l'ordre d'exécution des sept (07) tâches qui minimise le temps total de réalisation du travail. Illustrez à l'aide d'un diagramme. c. Pourront-ils rendre le travail à temps ?

الحل

a.

Par l'algorithme de Johnson pour F2CmaxF2\|C_{\max} : Tâches avec pi1pi2p_{i1} \leq p_{i2} : B (0.8 ≤ 1.5), A (2 ≤ 2) → en début, ordre croissant de pi1p_{i1}. Tâches avec pi1>pi2p_{i1} \gt p_{i2} : C (1.5 > 0.5), D (2 > 1) → en fin, ordre décroissant de pi2p_{i2}. Ordre : B, A, D, C. CmaxC_{\max} = par diagramme de Gantt.

Ordre optimal : B, A, D, C\boxed{\text{Ordre optimal : B, A, D, C}}

b.

E et F n'ont pas de phase d'analyse : on les place comme tâches sans contrainte de précédence pour l'étudiant 2. G n'a pas de phase de calculs : l'étudiant 1 la fait indépendamment. On optimise en intercalant G pour l'étudiant 1 et E, F pour l'étudiant 2 dans les temps morts.

c.

Calculer le CmaxC_{\max} total et vérifier s'il est 10\leq 10 jours.

Veˊrifier Cmax10\boxed{\text{Vérifier } C_{\max} \leq 10}

التمرين 2

Exercice 2 — Chaîne de Markov et urnes de Bernoulli

#markov-chains#transition-matrix#stationary-distribution#expected-value

Soit dd balles numérotées de 1 à dd et réparties dans deux urnes AA et BB. On tire un nombre ii au hasard entre 1 et dd et la balle nin^\circ i est changée d'urne. Soit XnX_n le nombre de balles dans l'urne AA après nn tirages indépendants.

  1. Montrer que (Xn)n(X_n)_n est une chaîne de Markov homogène dont on déterminera la matrice de transition PP. Montrer que la chaîne est irréductible, récurrente positive.
  2. On suppose X0B(d,1/2)X_0 \sim \mathcal{B}(d, 1/2), donner la loi de X1X_1. Que peut-on conclure ?
  3. Soit d=3d = 3 et T0T_0 le nombre de tirages nécessaires pour vider AA. a. Calculer P(T0=nX0=x)P(T_0 = n | X_0 = x), pour n=1,2,3n = 1, 2, 3. b. Calculer PP, P2P^2, P3P^3. Si μ(0)=(P(X0=i))i=0,,3=(1/4,1/4,1/4,1/4)\mu(0) = (P(X_0 = i))_{i=0,\ldots,3} = (1/4, 1/4, 1/4, 1/4), déterminer μ(j)=(P(Xj=i))i=0,,3\mu(j) = (P(X_j = i))_{i=0,\ldots,3}, j=1,2,3j = 1, 2, 3.
  4. Montrer qu'il existe deux constantes aa et bb telles que xE\forall x \in E, y=0dypxy=ax+b\sum_{y=0}^{d} y p_{xy} = ax + b. En déduire E(XnX0)E(X_n | X_0) et limnE(XnX0)\lim_{n \to \infty} E(X_n | X_0).
الحل

1.

Xn+1X_{n+1} dépend uniquement de XnX_n (Markov). Si Xn=kX_n = k, on tire une balle : probabilité k/dk/d qu'elle soit dans AA (retirée, Xn+1=k1X_{n+1} = k-1), probabilité (dk)/d(d-k)/d qu'elle soit dans BB (ajoutée, Xn+1=k+1X_{n+1} = k+1). La matrice PP est tridiagonale avec pk,k1=k/dp_{k,k-1} = k/d et pk,k+1=(dk)/dp_{k,k+1} = (d-k)/d. Irréductible car tous les états communiquent. Récurrente positive car l'espace d'états est fini.

pk,k+1=dkd,  pk,k1=kd\boxed{p_{k,k+1} = \frac{d-k}{d}, \; p_{k,k-1} = \frac{k}{d}}

2.

Si X0B(d,1/2)X_0 \sim \mathcal{B}(d, 1/2), alors X1B(d,1/2)X_1 \sim \mathcal{B}(d, 1/2) aussi (par symétrie). La loi B(d,1/2)\mathcal{B}(d, 1/2) est la loi stationnaire.

X1B(d,1/2) (distribution stationnaire)\boxed{X_1 \sim \mathcal{B}(d, 1/2) \text{ (distribution stationnaire)}}

3.

Pour d=3d=3, les calculs se font directement avec la matrice 4×44 \times 4.

4.

yypxy=(x1)xd+(x+1)dxd=xxd+dxdxd+xd...=x(12/d)+1\sum_y y p_{xy} = (x-1)\frac{x}{d} + (x+1)\frac{d-x}{d} = x - \frac{x}{d} + \frac{d-x}{d} - \frac{x}{d} + \frac{x}{d}... = x(1 - 2/d) + 1... Soit a=12/da = 1-2/d et b=1b = 1... Non : E[Xn+1Xn=x]=xxd...E[X_{n+1}|X_n=x] = x\frac{x}{d}... Recalculons : =(x1)xd+(x+1)dxd=x2xd+11...= (x-1)\frac{x}{d} + (x+1)\frac{d-x}{d} = x - \frac{2x}{d} + 1 - \frac{1}{... }. En fait =x(x1)+(x+1)(dx)d=x2x+dx+dx2xd=d+(d2)x2xd= \frac{x(x-1) + (x+1)(d-x)}{d} = \frac{x^2-x+dx+d-x^2-x}{d} = \frac{d + (d-2)x - 2x}{d}... Let me just write a=12/da = 1-2/d, b=1b = 1. Then E[XnX0]=anX0+b1an1ab1a=12/d=d/2E[X_n|X_0] = a^n X_0 + b\frac{1-a^n}{1-a} \to \frac{b}{1-a} = \frac{1}{2/d} = d/2.

limnE(XnX0)=d/2\boxed{\lim_{n \to \infty} E(X_n | X_0) = d/2}

التمرين 3

Exercice 3 — Programmation linéaire : production optimale et bi-objectif

#linear-programming#bi-objective#simplex#sensitivity-analysis

A) Une entreprise fabrique deux articles AA et BB sur une machine qu'elle a louée pour 15 heures et à l'aide de deux matières premières M1M_1 et M2M_2. Elle fabrique l'article AA à la cadence de 3 unités à l'heure et l'article BB à la cadence de 2 unités à l'heure. La fabrication d'une unité de AA nécessite 2 unités de M1M_1 et 3 unités de M2M_2 et celle d'une unité de BB nécessite 1 unité de M1M_1 et 2 unités de M2M_2. Le bénéfice rapporté pour chaque unité est de 6 pour l'article AA et 4 pour l'article BB. La quantité de matière première M1M_1 disponible est de 60 unités. Par contre, la matière M2M_2 est disponible, mais étant périssable, l'entreprise doit utiliser au moins 30 unités. Toutefois, pour des raisons commerciales, la fabrication des articles AA et BB est limitée à 25 articles pour chaque type.

  1. Formuler ce problème sous forme d'un programme linéaire (P)(P).
  2. Résoudre (P)(P). Donner la solution optimale XX^*, le bénéfice optimal et retrouver la base optimale.
  3. Que devient cette solution si la quantité de M1M_1 disponible était de 70 unités au lieu de 60 ? justifier. Quelle serait, dans ce cas, la quantité de matière M1M_1 restante si la production était optimale ?

B) L'entreprise se propose de maximiser la fabrication du produit BB, tout en maximisant le bénéfice.

  1. Formuler ce problème sous forme d'un programme linéaire bi-objectif (Q)(Q).
  2. Indiquer l'ensemble des sommets efficaces de (Q)(Q).
  3. Donner l'ensemble des solutions efficaces supportées de (Q)(Q).
الحل

A.1.

Variables : xAx_A = nombre d'articles A, xBx_B = nombre d'articles B.

maxZ=6xA+4xB\max Z = 6x_A + 4x_B s.c. {xA/3+xB/2152xA+xB603xA+2xB30xA25,  xB25xA,xB0\text{s.c. } \begin{cases} x_A/3 + x_B/2 \leq 15 \\ 2x_A + x_B \leq 60 \\ 3x_A + 2x_B \geq 30 \\ x_A \leq 25, \; x_B \leq 25 \\ x_A, x_B \geq 0 \end{cases}

A.2.

Résoudre graphiquement ou par le simplexe. La solution optimale est à un sommet du polyèdre.

A.3.

Analyse de sensibilité sur la contrainte M1M_1.

B.1.

Bi-objectif : max(xB,6xA+4xB)\max (x_B, 6x_A + 4x_B) sous les mêmes contraintes.

B.2-3.

Les sommets efficaces sont ceux non dominés. Les solutions efficaces supportées sont sur la frontière du polyèdre des objectifs réalisables.

Reˊsolution par meˊthode pondeˊreˊe ou ε-contrainte\boxed{\text{Résolution par méthode pondérée ou } \varepsilon\text{-contrainte}}