ATLAS Ouvrir dans l'atlas →

ExerciceMaths Sup – MPSI, chapitre 21 — Dénombrement et probabilités · Variance et covariance

Exercice 93

Calcul

Soit n⩾1n \geqslant 1 et (Xi,j)1⩽i,j⩽n(X_{i,j})_{1 \leqslant i,j \leqslant n} des variables aléatoires mutuellement indépendantes, de loi uniforme sur {−1,1}\{-1, 1\}. On considère la matrice aléatoire M=(Xi,j)1⩽i,j⩽nM = (X_{i,j})_{1 \leqslant i,j \leqslant n}, c'est-à-dire la variable aléatoire qui à ω\omega associe la matrice (Xi,j(ω))\big(X_{i,j}(\omega)\big).

1. Déterminer la loi de la trace tr⁡(M)\operatorname{tr}(M), puis calculer son espérance et sa variance.

2. Calculer la probabilité que MM soit de rang 11.

3. Déterminer l'espérance et la variance de det⁡(M)\det(M). On utilisera l'expression du déterminant comme somme sur 𝔖n\mathfrak{S}_{n} du chapitre 17.

Corrigé

Les n2n^{2} variables sont centrées : E(Xi,j)=0\mathrm{E}(X_{i,j}) = 0, et Xi,j2=1X_{i,j}^{2} = 1, donc V(Xi,j)=1\mathrm{V}(X_{i,j}) = 1.

1. On a tr⁡(M)=∑i=1nXi,i\operatorname{tr}(M) = \sum_{i=1}^{n} X_{i,i}. Les variables Bi=Xi,i+12B_{i} = \frac{X_{i,i} + 1}{2} sont mutuellement indépendantes (proposition 16), de loi ℬ(12)\mathcal{B}(\frac{1}{2}), et tr⁡(M)+n2=∑iBi\frac{\operatorname{tr}(M) + n}{2} = \sum_{i} B_{i} suit la loi ℬ(n,12)\mathcal{B}(n,\frac{1}{2}) (théorème 8). D'où P(tr⁡(M)=2k−n)=(nk)12n,k∈⟦0 ,n⟧\boxed{\mathrm{P}\big(\operatorname{tr}(M) = 2k - n\big) = \binom{n}{k}\frac{1}{2^{n}}, \qquad k \in \llbracket 0\,,n \rrbracket} Par linéarité et par le théorème 13, les Xi,iX_{i,i} étant indépendantes, E(tr⁡(M))=0,V(tr⁡(M))=n\boxed{\mathrm{E}\big(\operatorname{tr}(M)\big) = 0, \qquad \mathrm{V}\big(\operatorname{tr}(M)\big) = n} C'est la marche aléatoire de la section 21.5.7 : la trace est la position d'un marcheur après nn pas.

2. Comme tous les coefficients sont non nuls, M≠0M \neq 0 et rg(M)⩾1\mathrm{rg}(M) \geqslant 1. La matrice est de rang 11 si et seulement si toutes ses colonnes sont proportionnelles à la première, C1C_{1}, qui est non nulle (chapitre 16). Une colonne CjC_{j} à coefficients ±1\pm 1 s'écrit λC1\lambda C_{1} si et seulement si λ=±1\lambda = \pm 1, puisque les coefficients de λC1\lambda C_{1} sont ±λ\pm\lambda. Ainsi rg(M)=1\mathrm{rg}(M) = 1 si et seulement si Cj∈{C1,−C1}C_{j} \in \{C_{1}, -C_{1}\} pour tout j∈⟦2 ,n⟧j \in \llbracket 2\,,n \rrbracket. La matrice MM suit la loi uniforme sur les 2n22^{n^{2}} matrices à coefficients ±1\pm 1, par indépendance, et il y a 2n2^{n} choix de C1C_{1} et 22 choix pour chacune des n−1n - 1 autres colonnes : P(rg(M)=1)=2n×2n−12n2,P(rg(M)=1)=12(n−1)2\mathrm{P}\big(\mathrm{rg}(M) = 1\big) = \frac{2^{n} \times 2^{n-1}}{2^{n^{2}}}, \qquad \boxed{\mathrm{P}\big(\mathrm{rg}(M) = 1\big) = \frac{1}{2^{(n-1)^{2}}}} Pour n=1n = 1, on trouve 11, et pour n=2n = 2, 12\frac{1}{2} : une matrice 2×22 \times 2 à coefficients ±1\pm 1 est inversible une fois sur deux.

3. Par la définition du déterminant (chapitre 17), det⁡(M)=∑σ∈𝔖nε(σ) Tσ,Tσ=∏j=1nXσ(j),j.\det(M) = \sum_{\sigma \in \mathfrak{S}_{n}} \varepsilon(\sigma)\,T_{\sigma}, \qquad T_{\sigma} = \prod_{j=1}^{n} X_{\sigma(j),j}. Le produit TσT_{\sigma} porte sur nn coefficients de positions distinctes (σ(j),j)(\sigma(j), j), donc mutuellement indépendants : par la proposition 20, E(Tσ)=0\mathrm{E}(T_{\sigma}) = 0, et, par linéarité, E(det⁡(M))=0\boxed{\mathrm{E}\big(\det(M)\big) = 0}.

Soit σ≠τ\sigma \neq \tau. Dans le produit TσTτT_{\sigma}T_{\tau}, le facteur de la colonne jj vaut Xσ(j),j2=1X_{\sigma(j),j}^{2} = 1 si σ(j)=τ(j)\sigma(j) = \tau(j), et Xσ(j),jXτ(j),jX_{\sigma(j),j}X_{\tau(j),j} sinon. Ainsi TσTτT_{\sigma}T_{\tau} est le produit des coefficients situés aux positions (σ(j),j)(\sigma(j), j) et (τ(j),j)(\tau(j), j), pour les jj tels que σ(j)≠τ(j)\sigma(j) \neq \tau(j) ; ces positions sont deux à deux distinctes, et il y en a au moins deux, puisque σ≠τ\sigma \neq \tau. Par la proposition 20, E(TσTτ)=0\mathrm{E}(T_{\sigma}T_{\tau}) = 0, donc Cov⁡(Tσ,Tτ)=0\operatorname{Cov}(T_{\sigma}, T_{\tau}) = 0, alors que V(Tσ)=E(Tσ2)=1\mathrm{V}(T_{\sigma}) = \mathrm{E}(T_{\sigma}^{2}) = 1. Les n!n! variables ε(σ)Tσ\varepsilon(\sigma)T_{\sigma} sont deux à deux décorrélées, et le théorème 13 donne V(det⁡(M))=∑σ∈𝔖nε(σ)2V(Tσ)=n!\boxed{\mathrm{V}\big(\det(M)\big) = \sum_{\sigma \in \mathfrak{S}_{n}} \varepsilon(\sigma)^{2}\mathrm{V}(T_{\sigma}) = n!} Pour n=2n = 2, det⁡(M)\det(M) vaut 00 avec la probabilité 12\frac{1}{2}, et ±2\pm 2 avec la probabilité 14\frac{1}{4} chacun : la variance vaut bien 22. Pour n=3n = 3, l'énumération des 512512 matrices donne P(det⁡M=0)=58\mathrm{P}(\det M = 0) = \frac{5}{8} et P(det⁡M=±4)=316\mathrm{P}(\det M = \pm 4) = \frac{3}{16} chacun, de variance 66 ✓.

Piège. Les variables TσT_{\sigma} ne sont pas indépendantes : deux d'entre elles partagent en général des coefficients. Le théorème 13 n'a besoin que de leur décorrélation deux à deux, qui se vérifie par la proposition 20 appliquée aux coefficients, eux indépendants. L'erreur la plus répandue est d'invoquer une indépendance qui n'existe pas.

À retenir. L'écart type de det⁡(M)\det(M) vaut n!\sqrt{n!} : le déterminant d'une matrice aléatoire à coefficients ±1\pm 1 est en moyenne nul, mais il est typiquement de l'ordre de n!\sqrt{n!}, bien plus petit que n!n!, nombre de termes de la somme, en raison des compensations entre termes de signes opposés.

Ouvrir dans l'atlas interactif →