📚 الرئيسية

مسابقة دكتوراه 2021Université Mohamed El Bachir El Ibrahimi - Bordj Bou Arréridj — الموضوع 02

مسابقة تخصص · Analyse Numérique & Optimisation · المدة: 2سا

Concours d'accès à la formation doctorale LMD — Optimisation et analyse convexe — Épreuve de spécialité, Université Mohammed El Bachir El Ibrahimi, Faculté MI, Département de Mathématiques — Samedi 6 mars 2021 — Durée 2 heures.

التمرين 1

Exercice 1 — Coercivité et points critiques d'une fonction sur R²

#optimization#coercivity#critical-points#convexity

On considère la fonction ff définie sur R2\mathbb{R}^2 par :

f(x,y)=x2+y2+4(x2+y2116)2+2xy.f(x,y) = x^2 + y^2 + 4\left(x^2 + y^2 - 1 - \frac{1}{6}\right)^2 + 2xy.
  1. (2 pts) Montrer qu'il existe (α,β)R+2(\alpha, \beta) \in \mathbb{R}_+^2 tels que f(x,y)α(x,y)2+βf(x,y) \geq \alpha \|(x,y)\|^2 + \beta pour tous (x,y)R2(x,y) \in \mathbb{R}^2 et déduire que ff est coercive.
  2. (2 pts) Déterminer les points critiques et la nature de chacun.
  3. (1 pt) La fonction ff est-elle convexe ?
الحل

1.

On développe : f(x,y)=x2+y2+2xy+4(x2+y27/6)2f(x,y) = x^2+y^2+2xy + 4(x^2+y^2-7/6)^2. Le terme 4(x2+y27/6)204(x^2+y^2-7/6)^2 \geq 0. Par ailleurs x2+y2+2xy=(x+y)20x^2+y^2+2xy = (x+y)^2 \geq 0, mais on a aussi x2+y2(x+y)2/2x^2+y^2 \geq (x+y)^2/2. Pour la coercivité : quand x,y\|x,y\| \to \infty, le terme quartique 4(x2+y2)24(x^2+y^2)^2 domine. Donc f(x,y)+f(x,y) \to +\infty.

2.

f=0\nabla f = 0 donne un système. Les points critiques se trouvent en résolvant les équations f/x=0\partial f/\partial x = 0 et f/y=0\partial f/\partial y = 0.

3.

La matrice hessienne contient des termes quartiques, donc ff n'est pas quadratique. On vérifie si HfH_f est semi-définie positive partout. Le terme 4(x2+y27/6)24(x^2+y^2-7/6)^2 rend la hessienne non constante. En général, ff n'est pas convexe (à cause de la structure du polynôme de degré 4).

التمرين 2

Exercice 2 — Opérateur proximal et minimisation convexe

#proximal-operator#convex-optimization#projection#l1-norm

On considère Rn\mathbb{R}^n muni de la norme euclidienne. Soit CRnC \subset \mathbb{R}^n un ensemble convexe fermé non vide et f:RnRf : \mathbb{R}^n \to \mathbb{R} une fonction convexe. Pour xRnx \in \mathbb{R}^n on définit

yRn,Φf,x(y)=f(y)+12xy2.\forall y \in \mathbb{R}^n, \quad \Phi_{f,x}(y) = f(y) + \frac{1}{2}\|x - y\|^2.
  1. Montrer qu'il existe xCx^* \in C unique tel que Φf,x(x)=inf{Φf,x(y),yC}\Phi_{f,x}(x^*) = \inf\{\Phi_{f,x}(y), y \in C\}. On note Proxf(x)=x\text{Prox}_f(x) = x^*.
  2. À quoi correspond Proxf\text{Prox}_f lorsque f(x)=0f(x) = 0 pour tout xCx \in C ?
  3. Calculer Proxf\text{Prox}_f lorsque C=RnC = \mathbb{R}^n et f(x)=x1f(x) = \|x\|_1.
  4. Montrer que si xˉC\bar{x} \in C est un minimiseur de ff sur CC, alors Proxf(xˉ)=xˉ\text{Prox}_f(\bar{x}) = \bar{x}.
  5. On suppose f=g+hf = g + h avec gg convexe différentiable et hh convexe continue. Si xˉ\bar{x} minimise ff sur CC, alors xˉ=Proxτh(xˉτg(xˉ))\bar{x} = \text{Prox}_{\tau h}(\bar{x} - \tau \nabla g(\bar{x})) pour tout τ>0\tau \gt 0.
الحل

1.

Φf,x\Phi_{f,x} est strictement convexe (car 12xy2\frac{1}{2}\|x-y\|^2 est strictement convexe) et coercive. Donc elle admet un unique minimum sur CC (fermé convexe).

2.

Si f0f \equiv 0 : Proxf(x)=argminyC12xy2=PC(x)\text{Prox}_f(x) = \arg\min_{y \in C} \frac{1}{2}\|x-y\|^2 = P_C(x), la projection orthogonale sur CC.

Prox0=PC (projection sur C)\boxed{\text{Prox}_0 = P_C \text{ (projection sur } C)}

3.

Avec C=RnC = \mathbb{R}^n et f=1f = \|\cdot\|_1 : on minimise yi+12(xiyi)2\sum|y_i| + \frac{1}{2}\sum(x_i-y_i)^2 composante par composante. La solution est le seuillage doux :

(Prox1(x))i=sign(xi)max(xi1,0)\boxed{(\text{Prox}_{\|\cdot\|_1}(x))_i = \text{sign}(x_i)\max(|x_i|-1, 0)}

4.

xˉ\bar{x} minimise ff sur CC, donc f(xˉ)f(y)f(\bar{x}) \leq f(y) pour tout yCy \in C. Alors Φf,xˉ(xˉ)=f(xˉ)f(y)f(y)+12xˉy2=Φf,xˉ(y)\Phi_{f,\bar{x}}(\bar{x}) = f(\bar{x}) \leq f(y) \leq f(y) + \frac{1}{2}\|\bar{x}-y\|^2 = \Phi_{f,\bar{x}}(y).

5.

La condition d'optimalité de xˉ\bar{x} pour f=g+hf = g+h s'écrit 0g(xˉ)+h(xˉ)+NC(xˉ)0 \in \nabla g(\bar{x}) + \partial h(\bar{x}) + N_C(\bar{x}). Cela équivaut à xˉ=Proxτh(xˉτg(xˉ))\bar{x} = \text{Prox}_{\tau h}(\bar{x} - \tau\nabla g(\bar{x})) (point fixe de l'itération forward-backward).