Mathématiques 2: I Décompositions de Matrices

Mathématiques 2
2020
TSI
4 heures Calculatrice autorisée
Ce sujet se compose de deux problèmes indépendants l’un de l’autre.
I Décompositions de matrices
Notations et rappels
Dans tout ce problème, 𝑛 désigne un entier naturel supérieur ou égal à 2.
L’espace vectoriel ℝ𝑛 est muni de sa structure euclidienne canonique et on note ⟨⋅|⋅⟩ son produit scalaire.
ℳ𝑛 (ℝ) désigne l’espace vectoriel des matrices carrées de taille 𝑛 à coefficients réels.
𝐼𝑛 désigne la matrice identité de ℳ𝑛 (ℝ).
La transposée d’une matrice 𝑀 ∈ ℳ𝑛 (ℝ) est notée 𝑀 ⊤ .
On rappelle qu’une matrice 𝑀 ∈ ℳ𝑛 (ℝ) est orthogonale lorsque 𝑀 ⊤ 𝑀 = 𝑀 𝑀 ⊤ = 𝐼𝑛 .
Une matrice 𝐴 = (𝑎𝑖,𝑗 ) ∈ ℳ𝑛 (ℝ) est triangulaire supérieure lorsque 𝑎𝑖,𝑗 = 0 dès que 1 ⩽ 𝑗 < 𝑖 ⩽ 𝑛.
On pourra utiliser sans preuve les deux résultats suivants :
— le produit de deux matrices triangulaires supérieures est une matrice triangulaire supérieure ;
— l’inverse d’une matrice triangulaire supérieure inversible est triangulaire supérieure.
Objectifs
Ce problème étudie deux types de décompositions matricielles, d’abord pour une matrice inversible, puis pour
une matrice dont le polynôme caractéristique est scindé sur ℝ.
I.A – Un exemple introductif

1 −3 4
Soit 𝑃 = ⎛
⎜ 2 0 −1 ⎞ ⎟.
⎝ −2 3 1 ⎠
On note 𝑐1 , 𝑐2 et 𝑐3 les colonnes de 𝑃 considérées comme des vecteurs de ℝ3 .
Q 1. Justifier que la matrice 𝑃 est inversible. En déduire que la famille ℬ1 = (𝑐1 , 𝑐2 , 𝑐3 ) est une base de ℝ3 .
Q 2. Appliquer le procédé d’orthonormalisation de Gram-Schmidt à la base ℬ1 = (𝑐1 , 𝑐2 , 𝑐3 ) pour construire
une base orthonormée ℬ2 = (𝑢1 , 𝑢2 , 𝑢3 ) de ℝ3 .
Q 3. Soit 𝑄 la matrice de passage de la base canonique ℬ de ℝ3 à la base ℬ2 . Justifier que 𝑄−1 = 𝑄⊤ .
Q 4. Déterminer la matrice de passage 𝑅 de la base ℬ2 à la base ℬ1 .
On constate que 𝑅 est triangulaire supérieure à éléments diagonaux strictement positifs.
Q 5. Justifier que 𝑃 = 𝑄𝑅.
I.B – Cas général : décomposition 𝑄𝑅

Q 6. Soit 𝑃 ∈ ℳ𝑛 (ℝ) une matrice inversible. En s’inspirant de la démarche mise en place sur l’exemple,
montrer qu’il existe une matrice 𝑄 orthogonale et une matrice 𝑅 triangulaire supérieure à coefficients diagonaux
strictement positifs telles que 𝑃 = 𝑄𝑅.
Q 7. Soit 𝑏 un vecteur de ℝ𝑛 et 𝑃 une matrice inversible de ℳ𝑛 (ℝ). Expliquer l’intérêt de la décomposition
𝑃 = 𝑄𝑅, avec 𝑄 orthogonale et 𝑅 triangulaire supérieure à coefficients diagonaux strictement positifs, pour
résoudre le système linéaire 𝑃 𝑥 = 𝑏, d’inconnue 𝑥 ∈ ℝ𝑛 .
Les deux questions qui suivent permettent de démontrer l’unicité de la décomposition précédente.
Q 8. Soit 𝑀 ∈ ℳ𝑛 (ℝ) une matrice à la fois orthogonale et triangulaire supérieure à coefficients diagonaux
strictement positifs. En raisonnant de proche en proche de la première à la dernière colonne de 𝑀, montrer que
𝑀 = 𝐼𝑛 .
Q 9. On considère quatre matrices 𝑄1 , 𝑄2 , 𝑅1 , 𝑅2 de ℳ𝑛 (ℝ) telles que 𝑄1 et 𝑄2 sont orthogonales, 𝑅1 et
𝑅2 sont triangulaires supérieures à coefficients diagonaux strictement positifs et 𝑄1 𝑅1 = 𝑄2 𝑅2 . Montrer que
𝑄1 = 𝑄2 et 𝑅1 = 𝑅2 .
I.C – Décomposition d’une matrice dont le polynôme caractéristique est scindé sur ℝ
Q 10. Soit 𝐴 une matrice de ℳ𝑛 (ℝ) dont le polynôme caractéristique 𝜒𝐴 est scindé sur ℝ. En utilisant
la décomposition 𝑄𝑅 d’une matrice inversible bien choisie, démontrer qu’il existe une matrice 𝑄 ∈ ℳ𝑛 (ℝ)
orthogonale et une matrice 𝑇 ∈ ℳ𝑛 (ℝ) triangulaire supérieure telles que 𝐴 = 𝑄𝑇 𝑄⊤ .
2020-02-11 19:15:08 Page 1/4

Q 11. Donner un exemple de matrice 𝐴 ∈ ℳ3 (ℝ) diagonalisable et une décomposition 𝐴 = 𝑄𝑇 𝑄⊤ avec
𝑄 ∈ ℳ3 (ℝ) orthogonale et 𝑇 ∈ ℳ3 (ℝ) triangulaire supérieure à éléments diagonaux strictement positifs, mais
non diagonale.
Q 12. Donner un exemple de matrice 𝐴 ∈ ℳ3 (ℝ) diagonalisable et une décomposition 𝐴 = 𝑄𝑇 𝑄⊤ avec
𝑄 ∈ ℳ3 (ℝ) orthogonale et 𝑇 ∈ ℳ3 (ℝ) diagonale à éléments diagonaux strictement positifs.
II Calcul approché d’intégrales par quadrature

Notations et rappels
Dans tout ce problème, 𝑛 désigne un entier naturel.
1 si 𝑖 = 𝑗
Pour 0 ⩽ 𝑖, 𝑗 ⩽ 𝑛, le symbole de Kronecker 𝛿𝑖,𝑗 est défini par 𝛿𝑖,𝑗 = {
0 si 𝑖 ≠ 𝑗
On note ℝ[𝑋] l’ensemble des polynômes à coefficients réels et ℝ𝑛 [𝑋] l’ensemble des polynômes à coefficients
réels de degré inférieur ou égal à 𝑛. On identifie un polynôme et la fonction polynomiale associée, définie sur le
segment [−1, 1] à valeurs dans ℝ.
Pour toute fonction 𝑔 continue sur [−1, 1] et à valeurs dans ℝ, on note 𝑁 (𝑔) = max |𝑔(𝑥)|.
𝑥∈[−1,1]
Pour toute fonction 𝑔 de classe 𝒞 , on note 𝑔 sa dérivée 𝑘-ième.
𝑘 (𝑘)
Objectifs
On s’intéresse dans ce problème à une méthode numérique de calcul approché de
1
𝐼(𝑓) = ∫ 𝑓(𝑡) d𝑡
−1
pour une fonction 𝑓 continue sur [−1, 1] à valeurs dans ℝ. Le principe de cette méthode, dite par quadrature,
consiste à approcher 𝐼(𝑓) par une somme, pondérée par des poids (𝜔𝑖 )0⩽𝑖⩽𝑛 , de valeurs prises par la fonction 𝑓
en 𝑛 + 1 points distincts (𝑟0 , 𝑟1 , …, 𝑟𝑛 ) de l’intervalle [−1, 1].
𝑛
On note 𝑟 le (𝑛 + 1)-uplet (𝑟0 , 𝑟1 , …, 𝑟𝑛 ) et Σ𝑟 (𝑓) la somme ∑ 𝜔𝑘 𝑓(𝑟𝑘 ).
𝑘=0
Les réels (𝑟0 , …, 𝑟𝑛 ) ∈ [−1, 1]𝑛+1 sont appelés nœuds de la quadrature. La performance de cette méthode d’ap-
proximation dépend du nombre de nœuds, du choix de ces nœuds et de la régularité de la fonction 𝑓. Ces trois
aspects sont abordés dans ce problème.
II.A – Méthode de quadrature

Dans cette sous-partie, (𝑟0 , 𝑟1 , …, 𝑟𝑛 ) sont 𝑛 + 1 points distincts de l’intervalle [−1, 1].
II.A.1) Détermination des poids (𝜔𝑖 )
𝑋 − 𝑟𝑘
Pour tout entier 𝑖 ∈ ⟦0, 𝑛⟧, on note 𝐴𝑖 (𝑋) = ∏ .
𝑟 − 𝑟𝑘
𝑘≠𝑖 𝑖
Q 13. Vérifier que, pour tout couple (𝑖, 𝑗) ∈ ⟦0, 𝑛⟧2 , 𝐴𝑖 (𝑟𝑗 ) = 𝛿𝑖,𝑗 .
Q 14. Démontrer que (𝐴0 , 𝐴1 , …, 𝐴𝑛 ) est une base de ℝ𝑛 [𝑋] et donner la décomposition, dans cette base,
d’un polynôme 𝑃 quelconque de ℝ𝑛 [𝑋].
On se propose de démontrer qu’il existe un unique (𝑛 + 1)-uplet (𝜔0 , …, 𝜔𝑛 ) ∈ ℝ𝑛+1 pour lequel l’égalité
1
𝑛
∫ 𝑃 (𝑥) d𝑥 = ∑ 𝜔𝑘 𝑃 (𝑟𝑘 )
𝑘=0
−1
est valable pour toute fonction 𝑃 polynomiale de degré inférieur ou égal à 𝑛.

On raisonne par analyse-synthèse.
a) Analyse
On suppose l’existence d’un tel (𝑛 + 1)-uplet (𝜔0 , …, 𝜔𝑛 ) ∈ ℝ𝑛+1 .
Q 15. Pour tout entier 𝑖 ∈ ⟦0, 𝑛⟧, exprimer 𝜔𝑖 en fonction de 𝐴𝑖 .
Q 16. En déduire l’unicité de (𝜔0 , …, 𝜔𝑛 ).
b) Synthèse
Q 17. Démontrer que, pour les valeurs 𝜔0 , …, 𝜔𝑛 , déterminées à la question 15, l’égalité
1
𝑛
∫ 𝑃 (𝑥) d𝑥 = ∑ 𝜔𝑘 𝑃 (𝑟𝑘 )
𝑘=0
−1
est valable pour toute fonction 𝑃 polynomiale de degré inférieur ou égal à 𝑛.
2020-02-11 19:15:08 Page 2/4

Dans toute la suite, pour un jeu de nœuds 𝑟 = (𝑟0 , …, 𝑟𝑛 ) donné, on note 𝑄𝑟 le polynôme
𝑛
𝑄𝑟 = ∏ (𝑋 − 𝑟𝑘 ).
𝑘=0
Si (𝜔0 , …, 𝜔𝑛 ) sont les valeurs déterminées à la question 15, on rappelle que

𝑛
Σ𝑟 (𝑓) = ∑ 𝜔𝑘 𝑓(𝑟𝑘 ).
𝑘=0
Q 18. On choisit 𝑛 = 0 et 𝑟0 = 0. Donner l’interprétation géométrique, en termes d’aire, de Σ𝑟 (𝑓) lorsque 𝑓

est positive sur [−1, 1].
Q 19. On choisit 𝑛 = 1, 𝑟0 = −1 et 𝑟1 = 1. Donner l’interprétation géométrique, en termes d’aire, de Σ𝑟 (𝑓)
lorsque 𝑓 est positive sur [−1, 1].
II.A.2) Majoration de l’erreur
On suppose 𝑓 de classe 𝒞𝑛+1 sur [−1, 1] et on se propose de majorer l’erreur d’approximation |𝐼(𝑓) − Σ𝑟 (𝑓)|.
Q 20. Justifier l’existence de 𝑁 (𝑓 (𝑛+1) ) et de 𝑁 (𝑄𝑟 ).
𝑛
Q 21. Démontrer que 𝑇𝑟,𝑓 = ∑ 𝑓(𝑟𝑘 )𝐴𝑘 est l’unique polynôme de ℝ𝑛 [𝑋] vérifiant 𝑇𝑟,𝑓 (𝑟𝑖 ) = 𝑓(𝑟𝑖 ), pour
𝑘=0
tout entier 𝑖 ∈ ⟦0, 𝑛⟧.
On se propose de démontrer que, pour tout réel 𝑥 ∈ [−1, 1],
𝑁 (𝑓 (𝑛+1) )𝑁 (𝑄𝑟 )
|𝑓(𝑥) − 𝑇𝑟,𝑓 (𝑥)| ⩽ . (II.1)
(𝑛 + 1)!
On suppose que 𝑔 est une fonction de classe 𝒞𝑛+1 sur [−1, 1] s’annulant en au moins 𝑛 + 2 points distincts de
[−1, 1].
Q 22. Démontrer que 𝑔′ s’annule en au moins 𝑛 + 1 points distincts de ]−1, 1[.
Q 23. Démontrer que 𝑔(𝑛+1) s’annule en au moins un point de l’intervalle ]−1, 1[.
Q 24. Démontrer que l’inégalité (II.1) est vérifiée si 𝑥 ∈ {𝑟0 , …, 𝑟𝑛 }.
Q 25. On suppose que 𝑥 ∉ {𝑟0 , …, 𝑟𝑛 }. Montrer qu’il existe un réel 𝜆𝑥 pour lequel la fonction 𝑔𝑥 définie sur
[−1, 1] par
∀𝑡 ∈ [−1, 1], 𝑔𝑥 (𝑡) = 𝑓(𝑡) − 𝑇𝑟,𝑓 (𝑡) − 𝜆𝑥 𝑄𝑟 (𝑡)
vérifie 𝑔𝑥 (𝑥) = 0.
Q 26. En utilisant le résultat de la question 23, démontrer qu’il existe un réel 𝑐𝑥 ∈ ]−1, 1[ tel que
1
𝜆𝑥 = 𝑓 (𝑛+1) (𝑐𝑥 )
(𝑛 + 1)!
et conclure.
𝑁 (𝑓 (𝑛+1) )𝑁 (𝑄𝑟 )
Q 27. En déduire que |𝐼(𝑓) − Σ𝑟 (𝑓)| ⩽ 2 .
(𝑛 + 1)!
II.B – Choix d’un jeu de nœuds

II.B.1) Nœuds équidistants
On choisit un jeu 𝑠 = (𝑠0 , …, 𝑠𝑛 ) de 𝑛 + 1 nœuds équidistants dans l’intervalle [−1, 1] vérifiant 𝑠0 = −1 et 𝑠𝑛 = 1.
On cherche à minorer 𝑁 (𝑄𝑠 ).
𝑛
Q 28. On pose ℎ = 2/𝑛 et 𝑥 = −1 + 𝑡ℎ avec 𝑡 ∈ [0, 𝑛]. Exprimer |𝑄𝑠 (𝑥)| en fonction de 𝜑𝑛+1 (𝑡) = ∏ |𝑡 − 𝑘|.
𝑘=0
Q 29. Justifier que la fonction 𝜑𝑛+1 admet un maximum sur [0, 𝑛] et que ce maximum est atteint sur [0, 𝑛/2].
On admet dans la suite que ce maximum est atteint sur [0, 1].
2
Q 30. Démontrer que max 𝜑2 = √ .
[0,1] 3 3
1
Q 31. Démontrer par récurrence que, pour tout entier 𝑛 ⩾ 1, max 𝜑𝑛+1 ⩾ (𝑛 − 1)!.
[0,𝑛] 4
Q 32. En déduire une minoration de 𝑁 (𝑄𝑠 ).
2020-02-11 19:15:08 Page 3/4

II.B.2) Zéros des polynômes de Tchebychev
On cherche un nouveau jeu de nœuds 𝑐 = (𝑐0 , …, 𝑐𝑛 ) utilisant une famille de polynômes.
[0, 𝜋] → [−1, 1]
On rappelle que arccos est la fonction réciproque de la fonction ∣
𝜃 ↦ cos(𝜃)
Pour tout entier 𝑛 ∈ ℕ, on note 𝑇𝑛 : 𝑥 ↦ cos(𝑛 arccos(𝑥)).
Q 33. Préciser le domaine de définition de 𝑇𝑛 .
Q 34. Calculer, pour tout 𝑥 dans le domaine de définition, 𝑇0 (𝑥), 𝑇1 (𝑥), 𝑇2 (𝑥) et 𝑇3 (𝑥).
Q 35. Pour tout entier naturel 𝑛, calculer 𝑇𝑛 (−1), 𝑇𝑛 (0) et 𝑇𝑛 (1).
Q 36. Étudier la parité de 𝑇𝑛 en fonction de 𝑛.
Q 37. Pour tout 𝑥 dans le domaine de définition, démontrer que 𝑇𝑛+1 (𝑥) + 𝑇𝑛−1 (𝑥) = 2𝑥𝑇𝑛 (𝑥).
Q 38. En déduire que 𝑇𝑛 est une fonction polynomiale dont on précisera le degré et le coefficient dominant.
Pour tout entier naturel 𝑛, 𝑇𝑛 est le 𝑛-ième polynôme de Tchebychev.
Q 39. Montrer que 𝑇𝑛+1 admet 𝑛 + 1 racines distinctes dans l’intervalle [−1, 1].
On rappelle qu’après l’exécution de l’instruction Python import numpy as np, np.pi désigne la constante 𝜋 et
np.cos correspond à la fonction cosinus. La fonction np.cos peut s’appliquer à un tableau, elle produit alors un
nouveau tableau de même dimension dont les composantes sont les cosinus des composantes du tableau passé en
paramètre. Par ailleurs, l’expression np.linspace(x, y, n) construit un vecteur de n valeurs, régulièrement
espacées, la première valant x et la dernière y.
Q 40. Compléter la fonction Python Tchebychev(n) ci-dessous qui prend en argument un entier 𝑛 et renvoie
un couple de deux vecteurs (𝑢0 , …, 𝑢999 ) et (𝑦0 , …, 𝑦999 ), avec, pour tout 𝑘 ∈ ⟦0, 999⟧, 𝑦𝑘 = 𝑇𝑛 (𝑢𝑘 ).
import numpy as np
def Tchebychev(n):
T = np.linspace(np.pi/2, np.pi, 1000)
U = np.cos(T)
Y = . . .
return U, Y
En reliant les points de coordonnées (𝑢𝑘 , 𝑦𝑘 ) pour deux valeurs du paramètre 𝑛, on a obtenu les deux courbes
suivantes.
Courbe 1 Courbe 2
1 1
0,5 0,5
0 0
−0,5 −0,5
−1 −1
−1 −0,8 −0,6 −0,4 −0,2 0 −1 −0,8 −0,6 −0,4 −0,2 0
Q 41. Pour chacune de ces deux courbes, préciser, en la justifiant, la valeur utilisée pour le paramètre 𝑛.
Q 42. On choisit les racines obtenues en question 39 comme nœuds du jeu 𝑐 = (𝑐0 , …, 𝑐𝑛 ). Donner la valeur
de 𝑁 (𝑄𝑐 ).
II.B.3) Comparaison des ces deux jeux de nœuds
√ 𝑁 (𝑄𝑐 )
Q 43. En admettant que 𝑛! ∼ 2𝜋𝑛 e−𝑛 𝑛𝑛 , calculer la limite, lorsque 𝑛 tend vers +∞, du rapport .
𝑁 (𝑄𝑠 )
Q 44. À l’aide de ce résultat et de celui de la question 27, comparer la qualité des estimations de l’intégrale
sur [−1, 1] par quadrature, selon le choix de 𝑠 ou 𝑐 comme jeu de nœuds.
• • • FIN • • •
2020-02-11 19:15:08 Page 4/4

Mathématiques 2: I Décompositions de Matrices

Transféré par

Droits d'auteur :

Formats disponibles

Mathématiques 2: I Décompositions de Matrices

Transféré par

Informations du document

Description originale:

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Mathématiques 2: I Décompositions de Matrices

Transféré par

Droits d'auteur :

Formats disponibles

Mathématiques 2

I.A – Un exemple introductif

I.B – Cas général : décomposition 𝑄𝑅

2020-02-11 19:15:08 Page 1/4

II Calcul approché d’intégrales par quadrature

II.A – Méthode de quadrature

est valable pour toute fonction 𝑃 polynomiale de degré inférieur ou égal à 𝑛.

est valable pour toute fonction 𝑃 polynomiale de degré inférieur ou égal à 𝑛.

2020-02-11 19:15:08 Page 2/4

Si (𝜔0 , …, 𝜔𝑛 ) sont les valeurs déterminées à la question 15, on rappelle que

Q 18. On choisit 𝑛 = 0 et 𝑟0 = 0. Donner l’interprétation géométrique, en termes d’aire, de Σ𝑟 (𝑓) lorsque 𝑓

II.B – Choix d’un jeu de nœuds

2020-02-11 19:15:08 Page 3/4

2020-02-11 19:15:08 Page 4/4

Vous aimerez peut-être aussi