Algèbre linéaire/Devoir/Décomposition QR

De testwiki
Aller à la navigation Aller à la recherche

Modèle:Devoir Modèle:Wikipédia Soit A une matrice réelle inversible de taille n. n est muni de sa structure euclidienne canonique.

a) Méthode de Schmidt

Soient ' de matrice A dans la base canonique de n, ' la b.o.n. de n construite à partir de ' par l'Modèle:W[1], R la matrice de ' dans ' et Q la matrice de ' dans .

i) Montrer que Q est orthogonale.

ii) Montrer que R est triangulaire supérieure.

iii) Quelle équation relie A,Q,R ?

b) Méthode de Householder

Pour tout vecteur non nul v, on pose sv(x)=x2v,xv2v, et Hv= la matrice de sv dans la base canonique. Les matrices Hv ainsi obtenues sont appelées Modèle:W de taille n.

i) Démontrer que sv est une réflexion[2].

ii) Démontrer que si H est une matrice de Householder de taille nk alors (Ik00H) est une matrice de Householder de taille n.

iii) Soient e un vecteur unitaire et x un vecteur non nul. On pose v=x+αe avec α=±x. Montrer que sv(x)=αe. En choisissant convenablement e, x et α, montrer qu'il existe une matrice de Householder Q1 (de taille n) telle que Q1A soit de la forme (αL0A1), avec A1 matrice inversible de taille n1.

iv) En déduire (en précisant t) un algorithme permettant de construire des matrices de Householder Q1,,Qt (de taille n) telles que R:=QtQ1A soit triangulaire supérieure.

v) Quelle équation relie alors A,R, et Q:=Q1Qt ?

c)

i) Quelles matrices sont à la fois triangulaires et orthogonales ? Montrer qu'elles sont produits de matrices de réflexions[2]. En déduire (par b) que tout endomorphisme orthogonal est un produit de réflexions.

ii) Quelles matrices sont à la fois triangulaires et orthogonales et à termes diagonaux >0 ? En déduire (en améliorant b.iv) que tout endomorphisme orthogonal est produit d'un nombre n de réflexions. Modèle:Corrigé

Notes

Modèle:Références

  1. Rappel sur Gram-Schmidt pour la question a : pour toute base (u1,,un) d'un e.v. euclidien E, il existe une unique base orthonormée (v1,,vn) de E telle que pour tout k, vkVect(u1,,uk) et vk,uk>0.
  2. 2,0 et 2,1 Précision pour les questions b.i et c : une réflexion est une symétrie orthogonale par rapport à un hyperplan.