Mines-Ponts 2025 — Mathématiques II — MP
Critère de Schur-Cohn et généralisation au cas non inversible
Notations et objectifs
On fixe {n\geq1} et l’on note {\mathcal M_n(\mathbb R)}, {\mathcal S_n(\mathbb R)}, {\mathcal D_n(\mathbb R)} et {\operatorname{GL}_n(\mathbb R)} les ensembles des matrices réelles carrées, symétriques, diagonales et inversibles. Une matrice d’ordre {1} est identifiée à son coefficient.
On note {M^{\mathsf T}} la transposée, {M_{ij}} le coefficient de ligne {i}, colonne {j}, et {\pi(M)} le nombre de valeurs propres réelles strictement positives, comptées avec multiplicité. On utilise la notation {\operatorname{Diag}(u_1,\ldots,u_n)}.
Le PGCD unitaire de deux polynômes non simultanément nuls est noté {f\wedge g}. On identifie un polynôme à sa fonction associée. On note {\sigma(f)} le nombre de ses racines réelles dans {]-1,1[}, comptées avec multiplicité ; par exemple {\sigma(X^2(X-1)(X+1))=2}.
Une racine réelle {\alpha} de {f} est dite stable si {\alpha\ne0} et {f(\alpha)=f(\alpha^{-1})=0}.
Pour {f=\sum\limits_{k=0}^ma_kX^k} de degré {m}, son réciproque est {f_0=\sum\limits_{k=0}^ma_{m-k}X^k}.
On pose {U=(1,0,\ldots,0)^{\mathsf T}} et {S_{ij}=\delta_{i+1,j}} : {S} a des {1} juste au-dessus de la diagonale et des {0} ailleurs. On définit
{J(f)=f_0(S)^{\mathsf T}f_0(S)-f(S)^{\mathsf T}f(S).}Le polynôme {p=\sum\limits_{k=0}^na_kX^k}, avec {a_n\ne0}, est scindé sur {\mathbb R}. Ses racines, comptées avec multiplicité, sont {\alpha_1\leq\cdots\leq\alpha_n}.
On veut prouver {\sigma(p)=\pi(J(p))} lorsque {J(p)} est inversible, puis compter les racines dans {]-1,1[} dans le cas général. Ces résultats interviennent dans l’étude de la stabilité des systèmes dynamiques.
A — Polynôme réciproque et stabilité des racines
1. Montrer que {p_0(x)=x^np(1/x)} pour {x\ne0}, puis que {p_0=a_n\prod\limits_{j=1}^n(1-\alpha_jX)}.
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. Montrer que {p\wedge p_0=1} si et seulement si {p} n’a aucune racine stable.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez :
Pour poursuivre votre exploration, vous pouvez :
Jusqu’à la fin de A, toutes les racines de {p} sont supposées stables et simples.
3. Justifier qu’il existe {\lambda\in\{-1,1\}} tel que {p=\lambda p_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 :
On pose {h=Xp'}, de degré {n}, et l’on note {h_0} et {(p')_0} les réciproques respectifs de {h} et de {p'}.
4. Montrer que {h=np-\lambda(p')_0}, puis que {h_0=\lambda(np-Xp')}.
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. Vérifier que {p'} est scindé sur {\mathbb R}, montrer que {h\wedge h_0=1}, puis que {p'} n’a aucune racine stable.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez :
Pour poursuivre votre exploration, vous pouvez :
B — Liberté d’une famille de polynômes
Pour {1\leq j\leq n}, on pose
{f_j=a_n\prod_{k=j+1}^n(1-\alpha_kX)\prod_{k=1}^{j-1}(X-\alpha_k),}avec la convention du produit vide égal à {1}.
6. Si {1\leq i\lt k\leq n} et {\alpha_i\alpha_k=1}, montrer que {\alpha_i} est racine de chaque {f_j} et que la famille {(f_1,\ldots,f_n)} est lié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 :
Jusqu’à la fin de B, aucune racine de {p} n’est supposée stable.
On note {E} l’espace des fractions rationnelles réelles dont les pôles éventuels sont des inverses de racines de {p}. Leurs dénominateurs peuvent être choisis comme produits finis de facteurs {1-\alpha_iX}, le produit vide étant admis.
On pose
{\begin{array}{rl}\displaystyle g_j&=\displaystyle\frac{f_j}{\prod_{i=1}^n(1-\alpha_iX)},\\[12pt]\displaystyle P_j(f)&=\displaystyle\frac{(1-\alpha_jX)f-(1-\alpha_j^2)f(\alpha_j)}{X-\alpha_j}\end{array}}
7. Montrer que {P_j} est un endomorphisme de {E} et déterminer son noyau.
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. Pour {g\in E}, calculer {P_j((X-\alpha_j)g/(1-\alpha_jX))}.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez :
Pour poursuivre votre exploration, vous pouvez :
9. En déduire que {(f_1,\ldots,f_n)} est libre.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez :
Pour poursuivre votre exploration, vous pouvez :
C — Expression de {J(p)}
10. Montrer que {((S^{\mathsf T})^iU)_{0\leq i\leq n-1}} est une base de l’espace des colonnes réelles de taille {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 :
Pour {1\leq j\leq n}, on pose {B_j=S-\alpha_jI_n} et {C_j=I_n-\alpha_jS}.
11. Montrer que
{J(p)=\sum_{j=1}^nf_j(S)^{\mathsf T}(C_j^{\mathsf T}C_j-B_j^{\mathsf T}B_j)f_j(S).}
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. Montrer que {C_j^{\mathsf T}C_j-B_j^{\mathsf T}B_j=(1-\alpha_j^2)UU^{\mathsf T}}.
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. On pose {D=\operatorname{Diag}(1-\alpha_1^2,\ldots,1-\alpha_n^2)} et l’on définit {V} par ses colonnes {V_j=f_j(S^{\mathsf T})U}. Montrer que {J(p)=VDV^{\mathsf T}}.
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. En déduire, à l’aide de Q6, que si {p} possède une racine stable, alors {J(p)} n’est pas inversible.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez :
Pour poursuivre votre exploration, vous pouvez :
D — Cas inversible : critère de Schur-Cohn
On munit l’espace des colonnes de son produit scalaire euclidien. Un sous-espace {F} vérifie {(C_M)} si {X^{\mathsf T}MX\gt0} pour tout {X\in F\setminus\{0\}}.
On note {d(M)} la dimension maximale d’un tel sous-espace. Ce maximum existe : le sous-espace nul convient et les dimensions possibles appartiennent à {\{0,\ldots,n\}}.
15. Si {A=P^{\mathsf T}BP} avec {P} inversible, montrer que {d(B)\geq d(A)}, puis que {d(B)=d(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. Pour {M} symétrique, construire un sous-espace {F_M} de dimension {\pi(M)} vérifiant {(C_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 :
17. Si un sous-espace {G} de dimension strictement supérieure à {\pi(M)} vérifiait {(C_M)}, montrer que {\dim(F_M^\perp\cap G)\geq1} et en déduire que {d(M)=\pi(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 :
18. Démontrer le critère de Schur-Cohn : si {J(p)} est inversible, alors {p} n’a aucune racine stable et {\sigma(p)=\pi(J(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 :
E — Condition nécessaire et suffisante d’inversibilité
19. Si {p} n’a aucune racine stable mais {J(p)} n’est pas inversible, montrer à l’aide de Q9 et Q13 qu’il existe un polynôme non nul {q} de degré au plus {n-1} tel que {q(S^{\mathsf T})U=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 :
20. En déduire que {J(p)} est inversible si et seulement si {p} n’a aucune racine stable.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez :
Pour poursuivre votre exploration, vous pouvez :
F — Un cas particulier
Toutes les racines de {p} sont maintenant supposées stables et simples. On rappelle {h=Xp'} et {p=\lambda p_0} avec {\lambda\in\{-1,1\}}.
21. Montrer que {J(h)} est inversible.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez :
Pour poursuivre votre exploration, vous pouvez :
22. Montrer qu’il existe {\eta\gt0} tel que, pour {1-\eta\lt r\lt1}, le polynôme {p(rX)} est scindé, possède exactement {\sigma(p)} racines dans {]-1,1[} et n’a aucune racine stable.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez :
Pour poursuivre votre exploration, vous pouvez :
Pour {r\gt0}, on pose {F(r)=J(p(rX))}.
23. Montrer que
{\lim_{r\to1^-}\pi\left(\frac n{2(r-1)}F(r)\right)=n-\sigma(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 :
24. Justifier que
{F} est dérivable sur
{\mathbb R_+^*} et que
{\begin{array}{l}F'(1)=2n\,p(S)^{\mathsf T}p(S)\\[8pt]\quad-2S^{\mathsf T}p'(S)^{\mathsf T}p(S)-2p(S)^{\mathsf T}p'(S)S.\end{array}}
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez :
Pour poursuivre votre exploration, vous pouvez :
25. En déduire, à l’aide de Q4, que
{\frac n{2(r-1)}F(r)=J(h)+o(1)\quad(r\to1).}
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 admet que les valeurs propres d’une matrice symétrique, rangées dans l’ordre décroissant et comptées avec multiplicité, dépendent continûment de la matrice.
26. En déduire que {\sigma(p)=n-1-\pi(J(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 :
G — Méthode générale
On revient au cas général. On construit {f=p\wedge p_0} et {g} défini par {p=fg}.
27. Montrer que {\sigma(g)=\pi(J(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 :
28. Proposer une méthode construisant un nombre fini, éventuellement nul, de polynômes
{g_1,\ldots,g_\ell} dont les racines sont stables et simples, tels que
{f=g_1\cdots g_\ell}.
Exprimer {\sigma(p)} à l’aide de {n}, {\deg g}, {\pi(J(g))}, {\ell} et des {\pi(J(g_i'))}.
Cliquer ici pour voir (ou cacher) la réponse
Pour voir la suite de ce contenu, vous devez :
Pour poursuivre votre exploration, vous pouvez :