📚 الرئيسية

مسابقة دكتوراه 2022Université Abderrahmane Mira - Béjaïa — الموضوع 02

مسابقة تخصص · Recherche Opérationnelle · المدة: 1سا 30د

Concours d'entrée en Doctorat en Mathématiques Appliquées, Épreuve de Recherche Opérationnelle Fondamentale, Université de Béjaïa, Faculté des Sciences Exactes, Département de Recherche Opérationnelle, 24 février 2022, durée 01h30.

التمرين 1

Exercice 1 — Programmation quadratique (KKT) et inégalité valide de sac-à-dos

#quadratic-programming#kkt-conditions#convex-optimization#integer-programming

Partie 1. Une ressource disponible en quantité dd doit être affectée à trois activités en quantités x1,x2,x3x_1, x_2, x_3 respectivement. On suppose que la totalité de la quantité disponible dd doit être utilisée. L'allocation de xkx_k unités de ressource à l'activité kk procure un gain net évalué par zk(xk)=8xkkxk2z_k(x_k) = 8x_k - k x_k^2, k=1,2,3k=1,2,3. On souhaite déterminer une allocation qui procure un gain net total maximal.

  1. Donner la formulation mathématique du problème décrit, qu'on notera (P)(P).
  2. Le problème (P)(P) est-il quadratique ? Si c'est le cas donner ses paramètres.
  3. Est-ce que (P)(P) est un problème convexe ? Justifier.
  4. Écrire les conditions de KKT pour (P)(P).
  5. Résoudre le système KKT.
  6. La solution obtenue x0x^0 est-elle un point optimal local ou global ? Justifier.

Partie 2. Soit un problème du sac-à-dos (noté (Q)(Q)) formulé comme suit :

Z(x)=8x1+x2+5x3+2x4maxZ(x) = 8x_1 + x_2 + 5x_3 + 2x_4 \to \max

s.c.3x1+2x2+3x3+4x45,xj{0,1}, j=1,4.\text{s.c.}\quad 3x_1 + 2x_2 + 3x_3 + 4x_4 \le 5,\qquad x_j\in\{0,1\},\ j=\overline{1,4}.

  • Montrer que l'inégalité x1+x3+x41x_1 + x_3 + x_4 \le 1 est valide pour ce problème.
الحل

1.

Le gain total est k=13zk(xk)=8(x1+x2+x3)(x12+2x22+3x32)\sum_{k=1}^3 z_k(x_k) = 8(x_1+x_2+x_3) - (x_1^2 + 2x_2^2 + 3x_3^2), sous la contrainte d'utilisation totale et de positivité :

(P):max f(x)=8(x1+x2+x3)(x12+2x22+3x32)s.c.x1+x2+x3=d, xk0.(P):\quad \max\ f(x) = 8(x_1+x_2+x_3) - (x_1^2+2x_2^2+3x_3^2)\quad \text{s.c.}\quad x_1+x_2+x_3 = d,\ x_k\ge 0.

2.

Oui, (P)(P) est quadratique : f(x)=cxxDxf(x) = c^\top x - x^\top D x avec

c=(8,8,8),D=diag(1,2,3),2f=2D=diag(2,4,6).c = (8,8,8)^\top,\qquad D = \operatorname{diag}(1,2,3),\qquad \nabla^2 f = -2D = \operatorname{diag}(-2,-4,-6).

La matrice hessienne est constante, ce qui caractérise un programme quadratique.

3.

La hessienne 2f=diag(2,4,6)\nabla^2 f = \operatorname{diag}(-2,-4,-6) est définie négative, donc ff est strictement concave. On maximise une fonction concave sur un ensemble convexe (hyperplan xk=d\sum x_k = d intersecté avec l'orthant positif). C'est donc un problème convexe (maximisation concave).

 (P) est convexe : f strictement concave, domaine convexe \boxed{\ (P)\ \text{est convexe : } f\ \text{strictement concave, domaine convexe}\ }

4.

Lagrangien avec multiplicateur λ\lambda (égalité) et μk0\mu_k\ge 0 (contraintes xk0x_k\ge 0) :

fxkλ+μk=0,82kxkλ+μk=0,\frac{\partial f}{\partial x_k} - \lambda + \mu_k = 0,\qquad 8 - 2k x_k - \lambda + \mu_k = 0,

μk0,xk0,μkxk=0 (k=1,2,3),x1+x2+x3=d.\mu_k \ge 0,\quad x_k\ge 0,\quad \mu_k x_k = 0\ (k=1,2,3),\qquad x_1+x_2+x_3 = d.

5.

Cherchons la solution intérieure xk>0x_k>0, donc μk=0\mu_k = 0 et 82kxkλ=08 - 2k x_k - \lambda = 0, soit xk=8λ2kx_k = \dfrac{8-\lambda}{2k}. La contrainte donne

8λ2(1+12+13)=d  8λ2116=d  8λ=12d11.\frac{8-\lambda}{2}\left(1 + \tfrac12 + \tfrac13\right) = d\ \Longrightarrow\ \frac{8-\lambda}{2}\cdot\frac{11}{6} = d\ \Longrightarrow\ 8-\lambda = \frac{12d}{11}.

D'où λ=812d11\lambda = 8 - \dfrac{12d}{11} et

 x0=(6d11, 3d11, 2d11),f(x0)=8d6d211 \boxed{\ x^0 = \left(\tfrac{6d}{11},\ \tfrac{3d}{11},\ \tfrac{2d}{11}\right),\qquad f(x^0) = 8d - \tfrac{6d^2}{11}\ }

Ces valeurs sont positives pour d>0d>0 et μk=00\mu_k=0\ge0 : les conditions KKT sont satisfaites.

6.

ff étant strictement concave et le domaine convexe, tout point KKT est l'unique maximum global. Donc x0x^0 est le maximum global (et unique).

 x0 est l’optimum global unique de (P) \boxed{\ x^0\ \text{est l'optimum global unique de } (P)\ }

Partie 2.

Les poids des variables x1,x3,x4x_1, x_3, x_4 dans la contrainte sont respectivement 3,3,43, 3, 4. La somme des deux plus petits vaut 3+3=6>53+3 = 6 > 5 : deux quelconques de ces trois variables mises à 11 dépassent déjà la capacité 55. Donc dans toute solution admissible x{0,1}4x\in\{0,1\}^4, au plus une des variables x1,x3,x4x_1,x_3,x_4 peut valoir 11. Formellement, {x1,x3,x4}\{x_1,x_3,x_4\} forme une couverture minimale (tout sous-ensemble de taille 22 viole la contrainte), d'où l'inégalité de couverture

 x1+x3+x41 (valide pour (Q)) \boxed{\ x_1 + x_3 + x_4 \le 1\ \text{(valide pour } (Q))\ }

التمرين 2

Exercice 2 — Concurrence monopolistique : équilibre, coopération et jeu répété

#game-theory#nash-equilibrium#repeated-games#bertrand-competition

Deux entreprises sont en concurrence monopolistique sur le marché d'un bien différencié. Chaque entreprise choisit son prix de vente et la demande des consommateurs à l'entreprise i{1,2}i\in\{1,2\} est qi=1002pi+pjq_i = 100 - 2p_i + p_j, jij\ne i, où pkp_k est le prix choisi par l'entreprise k{1,2}k\in\{1,2\}. Les coûts de production sont supposés nuls.

Partie 1. On suppose que chaque entreprise choisit son prix de vente, sans communication avec sa concurrente. L'objectif de chacune est de maximiser son profit.

  1. Construire le jeu, qu'on notera Γ\Gamma, modélisant cette situation de concurrence entre les 2 entreprises. Préciser ses caractéristiques.
  2. Donner la définition du concept de solution du jeu.
  3. Déterminer l'équilibre de ce jeu, ainsi que les profits correspondants.

Partie 2.

  1. On suppose que les entreprises décident de coopérer pour former la grande coalition. Elles décident alors de fixer en commun un prix qui maximiserait la somme de leurs profits et de partager par la suite de manière égalitaire le profit de la coalition. Déterminer les stratégies et les profits de la solution « coopérative ».
  2. Supposons qu'une entreprise décide de ne pas respecter le prix issu de la coopération, pendant que l'autre le respecte. Dans ce cas, quel serait le prix qui serait le plus profitable pour la firme ii ?
  3. En comparant les profits de concurrence et de coopération pour chaque entreprise, quel est le comportement qui leur est le plus favorable ?

Partie 3. Nous considérons le cas où le marché peut continuer sans fin, et les deux entreprises sont appelées au début de chaque année à choisir simultanément et sans concertation leurs prix de vente pour l'année en cours. Supposons que chacune des deux entreprises adopte la stratégie (dite déclic) suivante : je commence la première année par coopérer en fixant le prix de l'issue de coopération ; je continue à coopérer tant que l'autre entreprise coopère en fixant son prix de l'issue de coopération. Dès qu'elle dévie en fixant un prix qui lui est plus favorable, alors je pratiquerai le prix de l'équilibre pour toute la suite du marché. On suppose que les deux entreprises ont un même taux d'escompte δ[0,1]\delta\in[0,1]. Sous quelles conditions sur δ\delta, l'issue de coopération peut-elle émerger comme un équilibre parfait en sous-jeu du jeu répété indéfiniment ?

الحل

Partie 1 — 1.

Le jeu Γ\Gamma est un jeu statique à information complète à la Bertrand (prix) différencié :

Γ=({1,2}, (Si)i, (πi)i),Si=[0,+) (prix pi0),\Gamma = \big(\{1,2\},\ (S_i)_{i},\ (\pi_i)_i\big),\quad S_i = [0,+\infty)\ (\text{prix } p_i\ge 0),

πi(pi,pj)=piqi=pi(1002pi+pj).\pi_i(p_i,p_j) = p_i\,q_i = p_i\,(100 - 2p_i + p_j).

Caractéristiques : 2 joueurs, stratégies continues (prix), jeu non coopératif à coups simultanés, information complète, fonctions de gain (profits) πi\pi_i.

Partie 1 — 2.

Le concept de solution est l'équilibre de Nash : un profil (p1,p2)(p_1^*,p_2^*) tel que pour chaque ii, pip_i^* maximise πi(pi,pj)\pi_i(p_i,p_j^*) étant donné pjp_j^* :

πi(pi,pj)πi(pi,pj)pi0, i.\pi_i(p_i^*,p_j^*)\ge \pi_i(p_i,p_j^*)\quad \forall p_i\ge 0,\ \forall i.

Partie 1 — 3.

Condition du premier ordre (fonction de meilleure réponse) :

πipi=1004pi+pj=0  pi=100+pj4.\frac{\partial \pi_i}{\partial p_i} = 100 - 4p_i + p_j = 0\ \Longrightarrow\ p_i = \frac{100 + p_j}{4}.

Par symétrie p1=p2=pp_1^* = p_2^* = p^*, d'où 4p=100+p4p^* = 100 + p^*, soit p=1003p^* = \tfrac{100}{3}. Alors qi=10021003+1003=2003q_i = 100 - 2\cdot\tfrac{100}{3} + \tfrac{100}{3} = \tfrac{200}{3} et

 p=100333,33,πN=10032003=2000092222,2 \boxed{\ p^* = \tfrac{100}{3}\approx 33{,}33,\qquad \pi^N = \tfrac{100}{3}\cdot\tfrac{200}{3} = \tfrac{20000}{9}\approx 2222{,}2\ }

Partie 2 — 4.

On maximise π1+π2\pi_1+\pi_2. Par symétrie p1=p2=pp_1=p_2=p donne qi=100pq_i = 100 - p et π1+π2=2p(100p)\pi_1+\pi_2 = 2p(100-p), maximal en

ddp2p(100p)=2004p=0  pc=50.\frac{d}{dp}\,2p(100-p) = 200 - 4p = 0\ \Longrightarrow\ p^c = 50.

(On vérifie l'optimum sans symétrie : (π1+π2)/p1=1004p1+2p2=0\partial(\pi_1+\pi_2)/\partial p_1 = 100 - 4p_1 + 2p_2 = 0 donne aussi p=50p=50.) Alors qi=50q_i = 50 et

 pc=50,πc=5050=2500 par entreprise \boxed{\ p^c = 50,\qquad \pi^c = 50\cdot 50 = 2500\ \text{par entreprise}\ }

Partie 2 — 5.

Si la firme jj reste à pj=50p_j = 50, la meilleure réponse de ii est pi=100+504=37,5p_i = \dfrac{100 + 50}{4} = 37{,}5. Alors qi=10075+50=75q_i = 100 - 75 + 50 = 75 et

 pid=37,5,πid=37,575=2812,5 \boxed{\ p_i^{d} = 37{,}5,\qquad \pi_i^{d} = 37{,}5\cdot 75 = 2812{,}5\ }

(la firme loyale obtient alors πj=5037,5=1875\pi_j = 50\cdot 37{,}5 = 1875).

Partie 2 — 6.

Comparaison des profits par entreprise :

πN=2222,2 < πc=2500 < πd=2812,5.\pi^N = 2222{,}2\ \lt\ \pi^c = 2500\ \lt\ \pi^{d} = 2812{,}5.

La coopération domine la concurrence (2500>22222500 > 2222), donc collectivement le comportement le plus favorable est de coopérer. Mais chaque firme est individuellement tentée de dévier (2812,5>25002812{,}5 > 2500) : c'est une structure de dilemme du prisonnier, et l'unique équilibre du jeu à un coup reste la concurrence.

Partie 3.

Avec la stratégie déclic (« grim trigger »), coopérer à jamais rapporte

VC=πc1δ=25001δ.V_C = \frac{\pi^c}{1-\delta} = \frac{2500}{1-\delta}.

Dévier une fois (gain πd=2812,5\pi^d = 2812{,}5) puis subir l'équilibre de Nash à jamais rapporte

VD=πd+δπN1δ.V_D = \pi^d + \frac{\delta\,\pi^N}{1-\delta}.

La coopération est un équilibre parfait en sous-jeu ssi VCVDV_C\ge V_D, c'est-à-dire πcπdδ(πNπd)\pi^c - \pi^d \ge \delta(\pi^N - \pi^d), soit

δπdπcπdπN=2812,525002812,52222,2=625/210625/18=917.\delta \ge \frac{\pi^d - \pi^c}{\pi^d - \pi^N} = \frac{2812{,}5 - 2500}{2812{,}5 - 2222{,}2} = \frac{625/2}{10625/18} = \frac{9}{17}.

 Coopeˊration soutenable en EPSJ    δ9170,53 \boxed{\ \text{Coopération soutenable en EPSJ}\iff \delta \ge \tfrac{9}{17}\approx 0{,}53\ }