2023-X-Ens-Psi-Corrigé

X-ENS — PSI — 2023

Notations et rappels

On note {\mathbb R_+} l’ensemble des réels positifs ou nuls, et {\mathbb R_+^*} celui des réels strictement positifs.
Une fonction continue sur un fermé borné non vide de {\mathbb R^n} y atteint son minimum.
Un sous-ensemble {C} d’un espace vectoriel réel est convexe si {(1-t)x+ty\in C} pour tous {x,y\in C} et {t\in[0,1]}.
Une fonction {f:C\to\mathbb R} est convexe si
{f((1-t)x+ty)\le(1-t)f(x)+tf(y).}Elle est strictement convexe si l’inégalité est stricte lorsque {x\ne y} et {0\lt t\lt1}.
Une fonction {f:A\times B\to\mathbb R} admet un point selle en {(a_*,b_*)} si, pour tous {a\in A} et {b\in B},
{f(a_*,b)\le f(a_*,b_*)\le f(a,b_*).}Toutes les variables aléatoires sont définies sur un même espace probabilisé {(\Omega,\mathcal F,P)}. Les parties I et II sont indépendantes.

I. Convexité et points selles

Soient {m,n\ge1} et {E} un espace vectoriel réel.
(1) Soit {C\subset E} convexe ; soient {f,g:C\to\mathbb R} convexes.

(1) (a) Montrer que {f+g} est convexe, et strictement convexe si l’une des deux fonctions est strictement convexe.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :
(1) (b) On suppose {f} strictement convexe. Vérifier que son minimum sur {C} est atteint en au plus un point.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :

(2) Soit {A\in\mathcal M_{m,n}(\mathbb R)}. On utilise les produits scalaires usuels sur {\mathbb R^n} et {\mathbb R^m} ; {A^\top} désigne la transposée.

(2) (a) Montrer que, pour {(x,\nu)\in\mathbb R^n\times\mathbb R^m},
{\langle Ax,\nu\rangle_{\mathbb R^m}=\langle x,A^\top\nu\rangle_{\mathbb R^n}.}
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :
(2) (b) En déduire {\ker A\subset(\operatorname{Im}A^\top)^\perp}, l’orthogonal étant pris dans {\mathbb R^n}.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :
(2) (c) Montrer que {\ker A=(\operatorname{Im}A^\top)^\perp}.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :

(3) Soit {U\subset\mathbb R^n} ouvert, {h:U\to\mathbb R} de classe {C^1}, et {b\in\mathbb R^m}. On suppose que {x_*\in U} minimise {h} sur
{V_b=\{x\in U:Ax+b=0\}.}

(3) (a) Montrer que, pour tout {u\in\ker A}, {\langle\nabla h(x_*),u\rangle=0}.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :
(3) (b) Montrer qu’il existe {\nu_*\in\mathbb R^m} tel que {\nabla h(x_*)-A^\top\nu_*=0}.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :
(3) (c) On définit {L(x,\nu)=h(x)-\langle\nu,Ax+b\rangle} sur {U\times\mathbb R^m}. En déduire que {\partial L/\partial x_k(x_*,\nu_*)=0} pour tout {1\le k\le n}.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :
(3) (d) Si {U} est convexe et {h} convexe sur {U}, conclure que {L} admet un point selle en {(x_*,\nu_*)}, c’est-à-dire
{L(x_*,\nu)\le L(x_*,\nu_*)\le L(x,\nu_*)}pour tout {(x,\nu)\in U\times\mathbb R^m}.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :

II. Entropie et codage

Soit {\mathcal X} un ensemble fini et {p=(p_x)_{x\in\mathcal X}} une loi de probabilité telle que {p_x\gt0} pour tout {x}. Son entropie est
{H(p)=-\sum_{x\in\mathcal X}p_x\ln p_x.}On pose {Q_{\mathcal X}=\{q\in\mathbb R^{\mathcal X}:q_x\ge0\}}.
Pour {q,q'\in Q_{\mathcal X}}, avec {q'_x\gt0} pour tout {x}, on définit
{\operatorname{KL}(q,q')=\sum_{x\in\mathcal X}q'_x\varphi(q_x/q'_x),}où {\varphi(t)=t\ln t-t+1} pour {t\gt0}, prolongée en {0} par continuité.

(4) (a) Préciser {\varphi(0)}.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :
(4) (b) Vérifier que {\varphi} est continue, strictement convexe, positive, et que {\varphi(t)=0} si et seulement si {t=1}.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :
(4) (c) Montrer que {Q_{\mathcal X}} est convexe et que {q\mapsto\operatorname{KL}(q,q')} est strictement convexe, positive, et s’annule si et seulement si {q=q'}.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :

Soit {A} un ensemble fini. Un mot {u=u_1\cdots u_n} sur {A} est une suite finie d’éléments de {A}, de longueur {|u|=n}. Le mot vide {\epsilon} a longueur {0}.
On note {A^*} l’ensemble des mots et {A^+=A^*\setminus\{\epsilon\}}. La concaténation {u\cdot v} juxtapose les lettres, le mot vide étant neutre. Un mot {u} est un préfixe de {v} s’il existe {w\in A^*} tel que {v=u\cdot w}.
Soit {\mathcal X} fini non vide et {c:\mathcal X\to\{0,1\}^+} injective : c’est un code binaire. On suppose ce code préfixe : si {x\ne y}, {c(x)} n’est pas un préfixe de {c(y)}.
(5) On définit {\bar c:\mathcal X\to\{0,1\}^*} en retirant la première lettre : {c(x)=c(x)_1\cdot\bar c(x)}.

(5) (a) Vérifier que, si {x\ne y} et {c(x)_1=c(y)_1}, alors {\bar c(x)\ne\bar c(y)} et {\bar c(x)} n’est pas un préfixe de {\bar c(y)}.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :
(5) (b) Pour {a\in\{0,1\}}, on pose {\mathcal X_a=\{x:c(x)_1=a\}}. Montrer que, si {\mathcal X_a} contient au moins deux éléments, la restriction de {\bar c} à {\mathcal X_a} est un code préfixe.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :
(5) (c) En déduire {\sum\limits_{x\in\mathcal X}2^{-|c(x)|}\le1}. On pourra séparer les sommes sur {\mathcal X_0} et {\mathcal X_1} et raisonner par récurrence sur {L(c)=\max_x|c(x)|}.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :

On pose {q_x=2^{-|c(x)|}}, et {X} est une variable aléatoire à valeurs dans {\mathcal X}, de loi {p}.

(6) (a) Vérifier que {\ln(2)E(|c(X)|)=-\sum\limits_{x\in\mathcal X}p_x\ln q_x}.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :
(6) (b) En déduire {E(|c(X)|)\ge H(p)/\ln2}. On pourra exprimer {\ln2\,E(|c(X)|)} à l’aide de {H(p)} et de {\operatorname{KL}(p,q)}.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :

III. Transport régularisé

Dans la suite, {I,J} sont deux ensembles finis. Les familles {\alpha=(\alpha_i)_{i\in I}} et {\beta=(\beta_j)_{j\in J}} ont des coordonnées strictement positives et vérifient
{\sum_{i\in I}\alpha_i=\sum_{j\in J}\beta_j=1.}Elles définissent des lois de probabilités. On pose {E=\mathbb R^{I\times J}} et {Q=\{q\in E:q_{ij}\ge0\}}.
L’ensemble {F(\alpha,\beta)} est formé des {q\in Q} tels que
{\sum_{j\in J}q_{ij}=\alpha_i,\quad\sum_{i\in I}q_{ij}=\beta_j.}La famille {p}, définie par {p_{ij}=\alpha_i\beta_j\gt0}, appartient à {F(\alpha,\beta)}.

(7) Vérifier que {F(\alpha,\beta)} est convexe dans {E}.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :

(8) Soient {X_1,X_2} deux variables aléatoires à valeurs respectives dans {I,J}.

(8) (a) Vérifier que, si {q\in F(\alpha,\beta)}, alors {\sum\limits_{i\in I}\sum\limits_{j\in J}q_{ij}=1}.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :
(8) (b) On suppose {P(X_1=i,X_2=j)=q_{ij}} avec {q\in F(\alpha,\beta)}. Calculer les lois de {X_1} et de {X_2}.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :
(8) (c) Que dire de {X_1} et {X_2} lorsque {q=p} ?
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :

Soient {C=(C_{ij})\in\mathbb R_+^{I\times J}} et {\varepsilon\gt0}. On définit sur {Q}{J_\varepsilon(q)=\sum_{i,j}q_{ij}C_{ij}+\varepsilon\operatorname{KL}(q,p),}où la divergence est celle de la partie II, avec {\mathcal X=I\times J}.

(9) Montrer que {J_\varepsilon} est strictement convexe sur {Q}.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :
(10) (a) Vérifier que {F(\alpha,\beta)} est un fermé borné de {\mathbb R^{I\times J}}.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :
(10) (b) Montrer qu’il existe un unique {q(\varepsilon)} minimisant {J_\varepsilon} sur {F(\alpha,\beta)}.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :
(10) (c) Donner un contre-exemple simple montrant que l’unicité n’est plus garantie si {\varepsilon=0}.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :
(11) (a) Vérifier que {q(\varepsilon)_{ij}\gt0} pour tout {(i,j)}. On pourra raisonner par l’absurde, considérer {q(\varepsilon,t)=(1-t)q(\varepsilon)+tp}, et observer {\varphi} au voisinage de {0}.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :
(11) (b) Montrer que cette propriété n’est plus garantie si {\varepsilon=0}.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :

IV. Dualité

On pose {Q_{\gt0}=(\mathbb R_+^*)^{I\times J}} et, sur {Q_{\gt0}\times(\mathbb R^I\times\mathbb R^J)},
{\begin{array}{rl}\mathcal L(q,(f,g))&=J_\varepsilon(q)\\[8pt]&\displaystyle\quad+\sum_i f_i\left(\alpha_i-\sum_jq_{ij}\right)\\[8pt]&\displaystyle\quad+\sum_j g_j\left(\beta_j-\sum_iq_{ij}\right).\end{array}}

(12) (a) Vérifier que {Q_{\gt0}} est un ouvert convexe de {\mathbb R^{I\times J}}.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :
(12) (b) Montrer qu’il existe {(f(\varepsilon),g(\varepsilon))} tel que {(q(\varepsilon),(f(\varepsilon),g(\varepsilon)))} soit un point selle de {\mathcal L}. On pourra identifier les espaces à {\mathbb R^n} et {\mathbb R^m}, avec {n=|I\times J|}, {m=|I|+|J|}, et utiliser Q3.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :
(13) (a) Montrer que, pour tout {(f,g)}, le minimum de {q\mapsto\mathcal L(q,(f,g))} sur {Q_{\gt0}} est atteint en
{q(f,g)_{ij}=p_{ij}e^{(f_i+g_j-C_{ij})/\varepsilon}.}
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :
(13) (b) Calculer {G(f,g)=\mathcal L(q(f,g),(f,g))}.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :
(13) (c) Vérifier que {G} est concave sur {\mathbb R^I\times\mathbb R^J}.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :
(14) On définit {f_*:\mathbb R^J\to\mathbb R^I} et {g_*:\mathbb R^I\to\mathbb R^J} par
{\begin{array}{l}\displaystyle f_*(g)_i=-\varepsilon\ln\left(\sum_j\beta_je^{(g_j-C_{ij})/\varepsilon}\right),\\[8pt]\displaystyle g_*(f)_j=-\varepsilon\ln\left(\sum_i\alpha_ie^{(f_i-C_{ij})/\varepsilon}\right).\end{array}}Vérifier, pour tout {(f,g)} et tout {(i,j)},
{\frac{\partial G}{\partial f_i}(f_*(g),g)=\frac{\partial G}{\partial g_j}(f,g_*(f))=0.}
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :

Soit {(f^0,g^0)\in\mathbb R^I\times\mathbb R^J}. Pour {k\ge0}, on définit
{g^{k+1}=g_*(f^k),\quad f^{k+1}=f_*(g^{k+1}).}La mention {\mathbb R^{I\times J}} pour le couple initial dans le sujet est une coquille de domaine : {f^0} est indexé par {I} et {g^0} par {J}.

(15) Montrer que la suite {(G(f^k,g^k))_{k\ge0}} est croissante.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :

(16) On suppose que {f^k_i\to f_i^\infty} et {g^k_j\to g_j^\infty} pour tous {i,j}. On note
{G_* =\sup_{(f,g)\in\mathbb R^I\times\mathbb R^J}G(f,g).}

(16) (a) Montrer que {G(f^\infty,g^\infty)=G_*}.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :
(16) (b) Montrer que {G(f(\varepsilon),g(\varepsilon))=G_*}.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :
(16) (c) Montrer qu’il existe {a\in\mathbb R} tel que, pour tous {i,j},
{f(\varepsilon)_i=f_i^\infty+a,\quad g(\varepsilon)_j=g_j^\infty-a.}
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :
(16) (d) En déduire que {q(f^k,g^k)\to q(\varepsilon)}.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez : Pour poursuivre votre exploration, vous pouvez :