ATLAS Ouvrir dans l'atlas →

Proposition 2Extrema liés

Trois applications des extrema liés

On applique le théorème des extrema liés (Multiplicateurs de Lagrange) à trois problèmes classiques. Dans chacun, l'existence d'un extremum vient d'un argument de compacité, ou se suppose, et les multiplicateurs désignent le point où il est atteint.

Énoncé

(i) Inégalité arithmético-géométrique. Pour tout entier n⩾1n \geqslant 1 et tous réels x1,…,xn⩾0x_1, \dots, x_n \geqslant 0, x1x2⋯xn⩽(x1+⋯+xnn)n,x_1 x_2 \cdots x_n \leqslant \Big(\frac{x_1 + \cdots + x_n}{n}\Big)^n, avec égalité si et seulement si x1=⋯=xnx_1 = \cdots = x_n.

(ii) Forme quadratique sur la sphère. Soit A∈ℳn(ℝ)A \in {\mathcal M_n}(\mathbb R) une matrice symétrique, n⩾1n \geqslant 1, et q(x)=txAxq(x) = {}^txAx. Sur la sphère unité S={x∈ℝn:x12+⋯+xn2=1}S = \{x \in \mathbb R^n : x_1^2 + \cdots + x_n^2 = 1\}, la fonction qq atteint un maximum MM et un minimum mm, qui sont des valeurs propres de AA, et toute valeur propre réelle de AA est comprise entre mm et MM.

Plus précisément, un point a∈Sa \in S vérifie ∇q(a)=λ∇g(a)\nabla q(a) = \lambda \nabla g(a), pour g(x)=∑ixi2−1g(x) = \sum_i x_i^2 - 1, si et seulement si Aa=λaAa = \lambda a, et alors λ=q(a)\lambda = q(a).

(iii) Distance à un ensemble de niveau. Avec les notations du nœud Ensembles de niveau, soit b∈ℝnb \in \mathbb R^n et a∈Γa \in \Gamma tel que dg(a)dg(a) soit surjective. Si la restriction à Γ\Gamma de x↦∥x−b∥2x \mapsto \|x - b\|_2 admet en aa un minimum local, alors a−ba - b est combinaison linéaire des gradients ∇g1(a),…,∇gp(a)\nabla g_1(a), \dots, \nabla g_p(a), et ⟨a−b,v⟩=0pour tout v∈TaΓ.\langle a - b, v \rangle = 0 \quad \text{pour tout } v \in T_a\Gamma.

Démonstration

(i) Le cas n=1n = 1 est une égalité. Soit n⩾2n \geqslant 2, ω=(1,…,1)\omega = (1, \dots, 1), f(x)=x1⋯xnf(x) = x_1 \cdots x_n et g(x)=x1+⋯+xn−ng(x) = x_1 + \cdots + x_n - n, fonctions polynomiales, donc de classe C∞C^\infty (Classe Cᵏ, (iii)).

L'ensemble U= ]0,+∞[nU = \,]0, +\infty[^n est ouvert, comme intersection des images réciproques de ]0,+∞[]0, +\infty[ par les projections. Sur UU, dg(x)h=h1+⋯+hndg(x)h = h_1 + \cdots + h_n est surjective, puisqu'elle vaut 11 en (1,0,…,0)(1, 0, \dots, 0), et l'on note Γ={x∈U:g(x)=0}\Gamma = \{x \in U : g(x) = 0\}.

Existence d'un maximum. La partie K={x∈ℝn:xi⩾0 pour tout i, g(x)=0}K = \{x \in \mathbb R^n : x_i \geqslant 0 \text{ pour tout } i,\ g(x) = 0\} est fermée, comme intersection d'images réciproques de fermés par des applications continues, et bornée, puisque 0⩽xi⩽n0 \leqslant x_i \leqslant n ; elle est donc compacte (L'espace ℝⁿ, (iii)), et elle contient ω\omega. La fonction continue ff y atteint un maximum MM en un point cc (Compacts de ℝ, (iii)).

On a M⩾f(ω)=1>0M \geqslant f(\omega) = 1 > 0, donc tous les cic_i sont strictement positifs, et c∈Γc \in \Gamma. Comme Γ⊂K\Gamma \subset K, la fonction ff admet en cc un maximum lié.

Le point de maximum. Par le théorème des extrema liés (Multiplicateurs de Lagrange, (ii)), il existe λ\lambda tel que ∇f(c)=λ ω\nabla f(c) = \lambda\, \omega, c'est-à-dire ∏j≠icj=λpour tout i.\prod_{j \neq i} c_j = \lambda \quad \text{pour tout } i. En multipliant par cic_i, on obtient M=λciM = \lambda c_i pour tout ii. Comme M>0M > 0, λ\lambda n'est pas nul, et tous les cic_i valent M/λM/\lambda ; leur somme étant nn, c=ωc = \omega et M=1M = 1.

Ainsi f⩽1f \leqslant 1 sur KK. Si x∈Kx \in K vérifie f(x)=1f(x) = 1, c'est un point de maximum de ff sur KK, et le raisonnement précédent, appliqué à xx au lieu de cc, donne x=ωx = \omega.

Cas général. Soient x1,…,xn⩾0x_1, \dots, x_n \geqslant 0 et s=x1+⋯+xns = x_1 + \cdots + x_n. Si s=0s = 0, tous les xix_i sont nuls et l'égalité a lieu. Sinon, y=(n/s) xy = (n/s)\,x appartient à KK, et f(y)=(n/s)nf(x)⩽1f(y) = (n/s)^n f(x) \leqslant 1 est l'inégalité annoncée. L'égalité a lieu si et seulement si y=ωy = \omega, c'est-à-dire si tous les xix_i valent s/ns/n.

(ii) La fonction gg est polynomiale, et dg(x)h=2⟨x,h⟩dg(x)h = 2\langle x, h \rangle. Pour x∈Sx \in S, dg(x)x=2≠0dg(x)x = 2 \neq 0 : la forme linéaire dg(x)dg(x) est surjective sur ℝ\mathbb R.

La sphère SS est fermée, comme image réciproque de {1}\{1\}, bornée, puisque ∣xi∣⩽1|x_i| \leqslant 1, et non vide : elle est compacte (L'espace ℝⁿ, (iii)), et la fonction polynomiale qq y atteint un maximum MM en un point aa et un minimum mm en un point a′a' (Compacts de ℝ, (iii)).

Comme q(x)=∑j,kajkxjxkq(x) = \sum_{j,k} a_{jk} x_j x_k, on a ∂jq(x)=∑k(ajk+akj)xk=2(Ax)j\partial_j q(x) = \sum_k (a_{jk} + a_{kj}) x_k = 2(Ax)_j par symétrie de AA, donc ∇q(x)=2Ax\nabla q(x) = 2Ax ; et ∇g(x)=2x\nabla g(x) = 2x. L'égalité ∇q(a)=λ∇g(a)\nabla q(a) = \lambda \nabla g(a) équivaut donc à Aa=λaAa = \lambda a, et alors, aa étant dans SS, λ=λ⟨a,a⟩=⟨Aa,a⟩=q(a).\lambda = \lambda\langle a, a \rangle = \langle Aa, a \rangle = q(a).

Le point aa est un maximum lié de qq ; le théorème des extrema liés (Multiplicateurs de Lagrange, (ii)) fournit λ\lambda tel que ∇q(a)=λ∇g(a)\nabla q(a) = \lambda \nabla g(a). Ainsi Aa=λaAa = \lambda a avec a≠0a \neq 0 : λ\lambda est une valeur propre de AA, et λ=q(a)=M\lambda = q(a) = M. De même pour mm, au point a′a'.

Enfin, si Av=μvAv = \mu v avec v≠0v \neq 0, le vecteur u=v/⟨v,v⟩u = v/\sqrt{\langle v, v \rangle} appartient à SS, et q(u)=μ⟨u,u⟩=μq(u) = \mu\langle u, u \rangle = \mu est compris entre mm et MM.

(iii) La fonction f(x)=∥x−b∥22=∑i(xi−bi)2f(x) = \|x - b\|_2^2 = \sum_i (x_i - b_i)^2 est polynomiale, et ∇f(x)=2(x−b)\nabla f(x) = 2(x - b). Comme t↦t2t \mapsto t^2 est croissante sur [0,+∞[[0, +\infty[, ff admet en aa un minimum local lié.

Le théorème des extrema liés donne des réels λi\lambda_i tels que 2(a−b)=∑iλi∇gi(a)2(a - b) = \sum_i \lambda_i \nabla g_i(a). Pour v∈TaΓ=Ker⁡dg(a)v \in T_a\Gamma = \operatorname{Ker} dg(a) (Ensembles de niveau, (ii)), on a ⟨∇gi(a),v⟩=dgi(a)v=0\langle \nabla g_i(a), v \rangle = dg_i(a)v = 0 pour tout ii, d'où ⟨a−b,v⟩=0\langle a - b, v \rangle = 0.

□
Remarque

Pour n=3n = 3 et la contrainte x+y+z=3x + y + z = 3, le point (1,1,1)(1, 1, 1), de multiplicateur λ=1\lambda = 1, est le maximum de xyzxyz sur la partie KK de (i) ; le second ordre confirme qu'il est un maximum local lié strict sur le plan, et montre pourquoi il faut restreindre la hessienne à l'espace tangent. Le lagrangien L=xyz−(x+y+z−3)L = xyz - (x + y + z - 3) a pour hessienne en (1,1,1)(1, 1, 1) HL=(011101110),H_L = \begin{pmatrix} 0 & 1 & 1 \\ 1 & 0 & 1 \\ 1 & 1 & 0 \end{pmatrix}, indéfinie, avec Q(1,1,1)=6>0Q(1, 1, 1) = 6 > 0 et Q(1,−1,0)=−2<0Q(1, -1, 0) = -2 < 0.

Mais sur l'espace tangent T={h:h1+h2+h3=0}T = \{h : h_1 + h_2 + h_3 = 0\}, on a Q(h)=2(h1h2+h1h3+h2h3)=(h1+h2+h3)2−(h12+h22+h32)=−(h12+h22+h32),\begin{aligned} Q(h) &= 2(h_1h_2 + h_1h_3 + h_2h_3) \\ &= (h_1 + h_2 + h_3)^2 - (h_1^2 + h_2^2 + h_3^2) \\ &= -(h_1^2 + h_2^2 + h_3^2), \end{aligned} forme définie négative, et le nœud Second ordre, (iii), donne un maximum local lié strict.

Sur ℝ3\mathbb R^3 tout entier, et non sur l'ouvert U= ]0,+∞[3U = \,]0, +\infty[^3, le système de Lagrange a d'autres solutions, comme (3,0,0)(3, 0, 0) avec λ=0\lambda = 0 ; elles ne sont pas des extrema liés, puisque xyzxyz change de signe au voisinage de ces points sur le plan x+y+z=3x + y + z = 3. C'est en se restreignant à UU que la preuve de (i) les écarte.

Remarque

Le second ordre départage aussi un minimum et un col. Pour A=(100020003)A = \begin{pmatrix} 1 & 0 & 0 \\ 0 & 2 & 0 \\ 0 & 0 & 3 \end{pmatrix} et q(x)=x12+2x22+3x32q(x) = x_1^2 + 2x_2^2 + 3x_3^2 sur la sphère de ℝ3\mathbb R^3, les points de la sphère où le système de Lagrange est satisfait sont, par (ii), les vecteurs propres unitaires ±δ1\pm\delta_1, ±δ2\pm\delta_2, ±δ3\pm\delta_3, de multiplicateurs 11, 22 et 33.

Au point aa de multiplicateur λ\lambda, le lagrangien q−λgq - \lambda g a pour hessienne 2(A−λI3)2(A - \lambda I_3), et l'espace tangent est le plan {h:⟨a,h⟩=0}\{h : \langle a, h \rangle = 0\}.

En ±δ1\pm\delta_1, la restriction de QQ au plan engendré par δ2\delta_2 et δ3\delta_3 est 2h22+4h322h_2^2 + 4h_3^2, définie positive : c'est un minimum local lié strict, qui vaut 11. En ±δ3\pm\delta_3, la restriction au plan de δ1\delta_1 et δ2\delta_2 est −4h12−2h22-4h_1^2 - 2h_2^2, et c'est un maximum, qui vaut 33.

En ±δ2\pm\delta_2, la restriction au plan de δ1\delta_1 et δ3\delta_3 est −2h12+2h32-2h_1^2 + 2h_3^2, indéfinie : par (iv) du nœud Second ordre, il n'y a pas d'extremum lié, ce qui se voit directement, puisque qq vaut 2−x12<22 - x_1^2 < 2 sur le grand cercle x3=0x_3 = 0 et 2+x32>22 + x_3^2 > 2 sur le grand cercle x1=0x_1 = 0. Par analogie avec le nœud Extrema, on dit, de façon imagée, que δ2\delta_2 est un col de qq sur la sphère.

Remarque

Le point (ii) contient l'étape décisive du théorème spectral : toute matrice symétrique réelle a une valeur propre réelle, et la plus grande est le maximum de txAx{}^txAx sur la sphère. Le nœud Théorème spectral de l'analyse fonctionnelle l'obtient sans calcul différentiel, par le quotient de Rayleigh, puis en déduit par récurrence une base orthonormée de vecteurs propres. L'exemple de xyxy sur le cercle x2+y2=2x^2 + y^2 = 2 (Multiplicateurs de Lagrange) en est un cas : la matrice (01/21/20)\begin{pmatrix} 0 & 1/2 \\ 1/2 & 0 \end{pmatrix} a pour valeurs propres ±12\pm\frac12, que le carré du rayon, 22, multiplie par 22.

Remarque

Le point (iii) dit que le segment qui joint bb à un point le plus proche de Γ\Gamma est orthogonal à l'espace tangent. Pour un hyperplan affine Γ={x:⟨u,x⟩=c}\Gamma = \{x : \langle u, x \rangle = c\}, avec u≠0u \neq 0, il donne a=b+tua = b + tu, et la contrainte impose t=(c−⟨u,b⟩)/⟨u,u⟩t = (c - \langle u, b \rangle)/\langle u, u \rangle.

Ce candidat est bien le point le plus proche : pour x∈Γx \in \Gamma, le vecteur x−ax - a vérifie ⟨u,x−a⟩=0\langle u, x - a \rangle = 0, donc il est orthogonal à a−ba - b, et ∥x−b∥22=∥x−a∥22+∥a−b∥22.\|x - b\|_2^2 = \|x - a\|_2^2 + \|a - b\|_2^2. On retrouve la distance ∣c−⟨u,b⟩∣/∥u∥2|c - \langle u, b \rangle| / \|u\|_2 d'un point à un hyperplan.

Repose sur

Sert à

Aucune fiche de ce graphe ne l'utilise encore.
Ouvrir dans l'atlas interactif →