Théorème de Lax-Milgram

Un article de Wikipédia, l'encyclopédie libre.

Sommaire

1 Énoncé
2 Démonstration
3 Applications

[modifier] Énoncé

Soient :

$\mathcal{H}$ un espace de Hilbert réel muni de son produit scalaire noté $\langle.,.\rangle$
une forme bilinéaire qui est
- continue sur $\mathcal{H}\times\mathcal{H}$ : $\exists\,c>0, \forall (u,v)\in \mathcal{H}^2\,,\ |a(u,v)|\leq c\|u\|\|v\|\quad$
- coercitive sur $\mathcal{H}$ : $\exists\,\alpha>0, \forall u\in\mathcal{H}\,,\ a(u,u) \geq \alpha\|u\|^2 \quad$
L une forme linéaire continue sur $\mathcal{H}$ .

Sous ces hypothèses il existe un unique $u$ de $\mathcal{H}$ tel que l'équation $a (u, v) = L v$ soit vérifiée pour tout $v$ de $\mathcal{H}$ :

$(1) \quad \exists!\ u \in \mathcal{H},\ \forall v\in\mathcal{H},\quad a(u,v)=Lv$

Si de plus la forme bilinéaire $a$ est symétrique, alors $u$ est l'unique élément de $\mathcal{H}$ qui minimise la fonctionnelle $J:\mathcal{H}\rightarrow\R$ définie par $J(v)\ =\ \frac{1}{2}a(v,v)-Lv$ pour tout $v$ de $\mathcal{H}$ , c'est-à-dire :

$(2) \quad \exists!\ u \in \mathcal{H},\quad J(u) = \min_{v\in\mathcal{H}}\ J(v)$

[modifier] Démonstration

[modifier] Cas général

Par application du théorème de Riesz sur les formes linéaires continues, il existe un unique $f\in\mathcal{H}$ tel que $Lv=\langle f,v\rangle$ pour tout $v\in\mathcal{H}$ .

Pour tout $u\in\mathcal{H}$ , l'application $v\mapsto a(u,v)$ est une forme linéaire continue sur $\mathcal{H}$ et donc de la même manière, il existe un unique élément $Au\in\mathcal{H}$ tel que $a(u,v)=\langle Au,v\rangle$ pour tout $v\in\mathcal{H}$ . On montre facilement que l'opérateur $A$ ainsi défini est un endomorphisme linéaire continu sur $\mathcal{H}$ . La relation $(1)$ s'écrit donc de manière équivalente :

$\exists!\ u \in \mathcal{H},\ Au=f$

Pour prouver cette proposition, il suffit donc de montrer que $A$ est une bijection de $\mathcal{H}$ sur $\mathcal{H}$ . On montre dans un premier temps que l'opérateur est injectif, puis qu'il est surjectif.

Par la coercivité de $a$ et en appliquant l'inégalité de Cauchy-Schwarz, on a pour tout $v\in\mathcal{H}$

$\alpha\|v\|^2 \leq a(v,v) = \langle Av,v\rangle \leq \|Av\|\|v\|$

d'où $\|Av\| \geq \alpha\|v\|$ pour tout $v$ de $\mathcal{H}$ (*), ce qui montre que $A$ est injectif.

Soit $\mathcal{Z}$ l'image de l'opérateur $A$ dans $\mathcal{H}$ . Comme $A$ est endomorphisme linéaire continu satisfaisant à l'inégalité précédente (*) , on en déduit que $\mathcal{Z}$ un sous-espace fermé de $\mathcal{H}$ et que $\mathcal{H}$ est la somme directe de $\mathcal{Z}$ et $\mathcal{Z}^{\perp}$ où $\mathcal{Z}^{\perp}$ est le sous-espace vectoriel orthogonal de $\mathcal{Z}$ .

Soit ensuite un élément $w$ de $\mathcal{Z}^{\perp}$ , on a par définition $\langle Aw,w\rangle = 0$ et donc :

$\alpha\|w\|^2 \leq a(w,w) = \langle Aw,w\rangle = 0$

d'où $w = 0$ . Ainsi, $\mathcal{Z}^{\perp}$ est réduit à ${0}$ , ce qui montre que $A$ est surjectif.

L'endomorphisme $A$ est bijectif, il existe donc un unique $u$ de $\mathcal{H}$ tel que $Au=f\,$ et il est donné par $\,u=A^{-1}f$ .

[modifier] Remarque

Sans calculer $u$ on a l'inégalité

$\|u\| \leq \frac{\|L\|'}{\alpha}$

où $\|\cdot\|'$ désigne la norme de l'espace dual $\mathcal{H}'$ .

[modifier] Cas Symétrique

Si la forme bilinéaire $a$ est symétrique, on a pour tout $w$ de $\mathcal{H}$ :

$J(u+w) = J(u)+\Big(a(u,w)-Lw\Big)+\frac{1}{2}a(w,w)$

Comme $u$ est l'unique solution de la proposition (1), cela donne

$J(u+w) = J(u)+\frac{1}{2}a(w,w)$

Et comme $a$ est coercive, on a :

$J(u+w) \geq J(u) + \frac{\alpha}{2}\|w\|^2$

On a donc $J(u) \leq J(v)$ pour tout $v\in\mathcal{H}$ , d'où le résultat $(2)$ .

[modifier] Applications

Ce théorème est à la base des méthodes aux éléments finis, on peut en effet montrer que si au lieu de chercher u dans l'on cherche u_n dans , un sous espace de de dimension finie n, alors d'une part :
- Dans le cas où $a$ est symétrique $u n$ est le projeté de $u$ au sens du produit scalaire définit par $a$
- Si l'on se donne $(\varphi_i)$ une base de $\mathcal{H}_n$ , le problème se ramène alors à la résolution d'un système linéaire :