ATLAS Ouvrir dans l'atlas →

ExerciceMaths Sup – MPSI, chapitre 19 — Espaces préhilbertiens réels · Bases orthonormées

Exercice 27

ProblèmeDémonstration

La matrice de Hilbert. Pour n∈ℕn \in \mathbb{N}, on note Hn=(1i+j+1)0⩽i,j⩽n∈ℳn+1(ℝ).H_{n} = \Big(\frac{1}{i + j + 1}\Big)_{0 \leqslant i, j \leqslant n} \in \mathcal{M}_{n + 1}(\mathbb{R}).

1. Montrer que HnH_{n} est la matrice de Gram de la famille (1,X,…,Xn)(1, X, \dots, X^{n}) pour le produit scalaire défini sur ℝ[X]\mathbb{R}[X] par ⟨P,Q⟩=∫01PQ\langle P, Q \rangle = \int_{0}^{1} PQ. En calculant XTHnX{X}^{T}H_{n}X pour une colonne X∈ℳn+1,1(ℝ)X \in \mathcal{M}_{n + 1,1}(\mathbb{R}), montrer que HnH_{n} est inversible. On comparera avec l'exercice 48 de la feuille du chapitre 17.

2. Montrer que, pour tout polynôme Q∈ℝ[X]Q \in \mathbb{R}[X], ∫−11Q(t) dt=−i∫0πQ(eiθ) eiθ dθ.\int_{-1}^{1} Q(t)\,\mathrm{d} t = -\mathrm{i}\int_{0}^{\pi} Q(\mathrm{e}^{\mathrm{i}\theta})\,\mathrm{e}^{\mathrm{i}\theta}\,\mathrm{d}\theta.

3. En déduire l'inégalité de Hilbert : pour tous réels a0a_{0}, …, ana_{n}, ∑0⩽k,ℓ⩽nakaℓk+ℓ+1⩽π∑k=0nak2.\sum_{0 \leqslant k, \ell \leqslant n} \frac{a_{k}a_{\ell}}{k + \ell + 1} \leqslant \pi\sum_{k = 0}^{n} a_{k}^{2}. On pourra utiliser l'égalité suivante, établie à l'exercice 52 de la feuille du chapitre 6, où l'on a noté P=∑k=0nakXkP = \sum_{k = 0}^{n} a_{k}X^{k} : ∫02π∣P(eiθ)∣2 dθ=2π∑k=0nak2.\int_{0}^{2\pi} \big|P(\mathrm{e}^{\mathrm{i}\theta})\big|^{2}\,\mathrm{d}\theta = 2\pi\sum_{k = 0}^{n} a_{k}^{2}.

Corrigé

1. Pour 0⩽i,j⩽n0 \leqslant i, j \leqslant n, la primitive usuelle de t↦ti+jt \mapsto t^{i + j} donne ⟨Xi,Xj⟩=∫01ti+j dt=1i+j+1,\langle X^{i}, X^{j} \rangle = \int_{0}^{1} t^{i + j}\,\mathrm{d} t = \frac{1}{i + j + 1}, et HnH_{n} est bien la matrice de Gram (⟨Xi,Xj⟩)i,j\big(\langle X^{i}, X^{j} \rangle\big)_{i, j} de la famille (1,X,…,Xn)(1, X, \dots, X^{n}). Soit X=(x0,…,xn)X = (x_{0}, \dots, x_{n}) une colonne, et P=∑ixiXiP = \sum_{i} x_{i}X^{i}. Par bilinéarité, XTHnX=∑i=0n∑j=0nxixj⟨Xi,Xj⟩=⟨∑ixiXi,∑jxjXj⟩=∫01P(t)2 dt.{X}^{T}H_{n}X = \sum_{i = 0}^{n}\sum_{j = 0}^{n} x_{i}x_{j}\langle X^{i}, X^{j} \rangle = \Big\langle \sum_{i} x_{i}X^{i}, \sum_{j} x_{j}X^{j}\Big\rangle = \int_{0}^{1} P(t)^{2}\,\mathrm{d} t. Si HnX=0H_{n}X = 0, alors XTHnX=0{X}^{T}H_{n}X = 0, donc ∥P∥2=0\lVert P \rVert^{2} = 0, et P=0P = 0, puisque la formule définit un produit scalaire sur ℝ[X]\mathbb{R}[X] (proposition 2 (ii)) : tous les xix_{i} sont nuls. Le noyau de HnH_{n} est réduit à {0}\{0\}, et HnH_{n} est inversible (chapitre 16, théorème 3). Hn=(⟨Xi,Xj⟩)0⩽i,j⩽n,XTHnX=∫01P2>0 si X≠0,Hn∈GLn+1(ℝ)\boxed{H_{n} = \big(\langle X^{i}, X^{j} \rangle\big)_{0 \leqslant i, j \leqslant n}, \qquad {X}^{T}H_{n}X = \int_{0}^{1} P^{2} > 0 \text{ si } X \neq 0, \qquad H_{n} \in \mathrm{GL}_{n + 1}(\mathbb{R})} L'exercice 48 de la feuille du chapitre 17 obtenait l'inversibilité en calculant le déterminant de Cauchy, après de lourdes manipulations de lignes et de colonnes. Ici, aucun calcul : l'inversibilité vient de ce que (1,X,…,Xn)(1, X, \dots, X^{n}) est libre et que l'intégrale de P2P^{2} est un produit scalaire. La proposition 9 (ii), appliquée dans une base orthonormée de l'espace euclidien ℝn[X]\mathbb{R}_{n}[X], donne même det⁡Hn>0\det H_{n} > 0 ; le calcul donne det⁡H1=112\det H_{1} = \frac{1}{12}, det⁡H2=12160\det H_{2} = \frac{1}{2160} et det⁡H3=16 048 000\det H_{3} = \frac{1}{6\,048\,000}, valeurs très petites qui expliquent la réputation numérique de cette matrice. ✓

2. Les deux membres dépendent linéairement de QQ, par linéarité de l'intégrale (chapitre 18, théorème 3, et proposition 8 pour les fonctions complexes) : il suffit de vérifier l'égalité pour Q=XkQ = X^{k}, k∈ℕk \in \mathbb{N}. On intègre au moyen des primitives usuelles : celle de tkt^{k} à gauche, celle de θ↦ei(k+1)θ\theta \mapsto \mathrm{e}^{\mathrm{i}(k + 1)\theta}, qui est θ↦ei(k+1)θi(k+1)\theta \mapsto \frac{\mathrm{e}^{\mathrm{i}(k + 1)\theta}}{\mathrm{i}(k + 1)}, à droite. ∫−11tk dt=1−(−1)k+1k+1=1+(−1)kk+1,−i∫0πei(k+1)θ dθ=−i⋅(−1)k+1−1i(k+1)=1+(−1)kk+1.\int_{-1}^{1} t^{k}\,\mathrm{d} t = \frac{1 - (-1)^{k + 1}}{k + 1} = \frac{1 + (-1)^{k}}{k + 1}, \qquad -\mathrm{i}\int_{0}^{\pi} \mathrm{e}^{\mathrm{i}(k + 1)\theta}\,\mathrm{d}\theta = -\mathrm{i} \cdot \frac{(-1)^{k + 1} - 1}{\mathrm{i}(k + 1)} = \frac{1 + (-1)^{k}}{k + 1}.

∀ Q∈ℝ[X],∫−11Q(t) dt=−i∫0πQ(eiθ) eiθ dθ\boxed{\forall\,Q \in \mathbb{R}[X], \quad \int_{-1}^{1} Q(t)\,\mathrm{d} t = -\mathrm{i}\int_{0}^{\pi} Q(\mathrm{e}^{\mathrm{i}\theta})\,\mathrm{e}^{\mathrm{i}\theta}\,\mathrm{d}\theta} Pour Q=(1+X)3+2X4Q = (1 + X)^{3} + 2X^{4}, les deux membres valent 245\frac{24}{5}, d'après un calcul symbolique. ✓ On peut lire cette égalité comme le fait qu'on va de −1-1 à 11 « par le segment » ou « par le demi-cercle supérieur » avec la même intégrale : c'est une première rencontre avec l'intégration le long d'un chemin du plan complexe.

3. Avec X=(a0,…,an)X = (a_{0}, \dots, a_{n}), la question 1 donne ∑k,ℓakaℓk+ℓ+1=∫01P2\sum_{k, \ell} \frac{a_{k}a_{\ell}}{k + \ell + 1} = \int_{0}^{1} P^{2}. Comme P2⩾0P^{2} \geqslant 0 sur [−1 ,0][-1\,,0], la relation de Chasles et la positivité de l'intégrale donnent ∫01P2⩽∫−11P2\int_{0}^{1} P^{2} \leqslant \int_{-1}^{1} P^{2}. On applique la question 2 au polynôme Q=P2Q = P^{2} ; le premier membre étant un réel positif, il est égal à son module, et l'inégalité triangulaire complexe (chapitre 18, théorème 5), avec ∣eiθ∣=1|\mathrm{e}^{\mathrm{i}\theta}| = 1, donne ∫−11P2=∣∫0πP(eiθ)2 eiθ dθ∣⩽∫0π∣P(eiθ)∣2 dθ.\int_{-1}^{1} P^{2} = \Big|\int_{0}^{\pi} P(\mathrm{e}^{\mathrm{i}\theta})^{2}\,\mathrm{e}^{\mathrm{i}\theta}\,\mathrm{d}\theta\Big| \leqslant \int_{0}^{\pi} \big|P(\mathrm{e}^{\mathrm{i}\theta})\big|^{2}\,\mathrm{d}\theta. Les coefficients aka_{k} étant réels, P(e−iθ)=P(eiθ)‾P(\mathrm{e}^{-\mathrm{i}\theta}) = \overline{P(\mathrm{e}^{\mathrm{i}\theta})} : la fonction h:θ↦∣P(eiθ)∣2h : \theta \mapsto |P(\mathrm{e}^{\mathrm{i}\theta})|^{2} est paire, et 2π2\pi-périodique. La proposition 6 du chapitre 6 donne alors ∫0πh=12∫−ππh=12∫02πh=12×2π∑k=0nak2,\int_{0}^{\pi} h = \frac{1}{2}\int_{-\pi}^{\pi} h = \frac{1}{2}\int_{0}^{2\pi} h = \frac{1}{2} \times 2\pi\sum_{k = 0}^{n} a_{k}^{2}, la dernière égalité étant celle de l'indication. En enchaînant les inégalités : ∑0⩽k,ℓ⩽nakaℓk+ℓ+1=∫01P2⩽∫−11P2⩽∫0π∣P(eiθ)∣2 dθ=π∑k=0nak2\boxed{\sum_{0 \leqslant k, \ell \leqslant n} \frac{a_{k}a_{\ell}}{k + \ell + 1} = \int_{0}^{1} P^{2} \leqslant \int_{-1}^{1} P^{2} \leqslant \int_{0}^{\pi} \big|P(\mathrm{e}^{\mathrm{i}\theta})\big|^{2}\,\mathrm{d}\theta = \pi\sum_{k = 0}^{n} a_{k}^{2}} Pour (a0,a1,a2,a3)=(1,−2,3,1)(a_{0}, a_{1}, a_{2}, a_{3}) = (1, -2, 3, 1), la chaîne s'écrit 8342⩽23021⩽15π\frac{83}{42} \leqslant \frac{230}{21} \leqslant 15\pi, ce que confirme un calcul symbolique. ✓

À retenir. La constante π\pi ne dépend pas de nn : c'est ce qui fait la force de l'inégalité. On peut montrer, ce que l'on admet, qu'elle ne peut pas être remplacée par une constante plus petite valable pour tout nn ; la plus grande valeur du quotient XTHnXXTX\frac{{X}^{T}H_{n}X}{{X}^{T}X}, calculée numériquement, vaut environ 1,271{,}27 pour n=1n = 1, 1,921{,}92 pour n=20n = 20 et 2,272{,}27 pour n=200n = 200 : elle s'approche de π\pi très lentement.

Dans l'édifice

Thèmes

bases orthonormées, produits scalaires, produits scalaires intégraux, polynômes orthogonaux, matrices et produit scalaire, matrices de gram

Notions

base orthogonale, produit scalaire

Ouvrir dans l'atlas interactif →