Quotient de Rayleigh

Quotient de Rayleigh

Le quotient de Rayleigh est un nombre réel caractérisant l'effet d'une matrice symétrique (respectivement hermitienne) sur un vecteur, et offrant les deux propriétés fondamentales suivantes :

  • le quotient de Rayleigh atteint un extremum relatif pour les vecteurs propres de la matrice
  • l’application du quotient de Rayleigh à un vecteur propre donne la valeur propre correspondante.

Ces deux propriétés peuvent servir à déterminer numériquement les Valeur propre, vecteur propre et espace propre d'un opérateur hermitien ou symétrique.

Le quotient de Rayleigh, dont la propriété d'extremum peut être reliée au principe du minimum de l'énergie potentielle en mécanique, a été étudié pour la première fois par Rayleigh (1877). Walter Ritz reprit l'idée en 1909 pour en faire la base d’une méthode d’approximation variationnelle.

Sommaire

Définitions et propriétés caractéristiques

En mathématiques, pour une matrice hermitienne A à coefficients complexes et un vecteur x non nul, on appelle quotient de Rayleigh R(A,x) le scalaire :

R(A,x) = \frac{x^{*} A x}{x^{*} x}.

x * désigne le vecteur adjoint de x, c'est-à-dire le conjugué du vecteur transposé.

Pour des matrices et des vecteurs à coefficients réels, le caractère hermitien se traduit par la caractère symétrique, et il faut substituer à la matrice adjointe x * la transposée familière x'. Notons que

pour tout nombre réel c, R(A,cx) = R(A,x).

Rappelons également que les valeurs propres d'une matrice symétrique (respectivement hermitienne) sont toutes réelles. On peut montrer[1] que, sous ces conditions,

  1. le quotient de Rayleigh atteint un minimum λmin (qui n'est autre que la plus petite valeur propre de A) lorsque x est un vecteur propre vmin associé à cette valeur.
  2. En outre, R(A, x) \leq \lambda_{\text{max}} et R(A,vmax) = λmax. Le quotient de Rayleigh, combiné au théorème du minimax, permet de déterminer une à une toutes les valeurs propres d'une matrice. On peut également l'employer pour calculer une valeur approchée d'une valeur propre à partir d'une approximation d'un vecteur propre. Ces idées forment d'ailleurs la base de l’algorithme d’itération de Rayleigh.

Cas particulier des matrices de covariance

On peut factoriser une matrice de covariance Σ sous la forme A'A. Les valeurs propres d'une telle matrice sont positives, car :

Σvi = λivi
A'Avi = λivi
vi'A'Avi = viivi
 \left\| A v_i \right\|^2 = \lambda _i \left\| v_i \right\|^2
 \lambda _i = \frac{\left\| A v_i \right\|^2}{\left\| v_i \right\|^2} \geq 0.

Les vecteurs propres sont deux à deux orthogonaux :

Σvi = λivi
vjvi = λivj'vi
vj)'vi = λivj'vi
λjvj'vi = λivj'vi
j − λi)vj'vi = 0
vj'vi = 0 (pour des valeurs propres différentes ; en cas de multiplicité, la base de Im(Σ) peut être orthogonalisée).

Les vecteurs propres formant une base, on peut exprimer le quotient de Rayleigh en fonction des valeurs propres en décomposant un vecteur quelconque x sur cette base des vecteurs propres :

x = \sum _{i=1} ^n \alpha _i v_i
\rho = \frac{x' A' A x}{x' x}
\rho = \frac{(\sum _{j=1} ^n \alpha _j v_j)' A' A (\sum _{i=1} ^n \alpha _i v_i)}{(\sum _{j=1} ^n \alpha _j v_j)' (\sum _{i=1} ^n \alpha _i v_i)}

relation qui, de par l’orthogonalité des vecteurs propres, devient :

\rho = \frac{\sum _{i=1} ^n \alpha _i ^2 \lambda _i}{\sum _{i=1} ^n \alpha _i ^2}.

Si un vecteur x rend maximal le scalaire ρ, alors tout vecteur k.x (k \ne 0) produira également ce maximum : il suffit donc de rechercher le vecteur propre de norme égale à 1 ; on peut donc imposer \textstyle\sum _{i=1} ^n \alpha _i ^2 = 1.

Par la relation ci-dessus, il apparaît que la recherche d'un vecteur propre se ramène à la recherche d'un extremum de Lagrange, à savoir \textstyle\sum _{i=1} ^n \alpha _i ^2 \lambda _i \rightarrow \max sous la contrainte \textstyle<\sum _{i=1} ^n \alpha _i ^2 = 1.

Puisque toutes les valeurs propres sont positives, le problème est convexe et le maximum est atteint à la frontière du domaine, à savoir quand α1 = 1 et \forall i > 1, \alpha _i = 0 (en supposant que les valeurs propres soit numérotées par ordre de valeur décroissante).

On aboutit au même résultat par la méthode des multiplicateur de Lagrange. Le problème consiste alors à déterminer les points critiques de la fonction

ρ(x) = xTΣx, sous la contrainte \|x\|^2 = x^Tx = 1.

c’est-à-dire trouver les points critiques de

\mathcal{L}(x) = x^T\Sigma x  -\lambda (x^Tx - 1) ,

λ est un multiplicateur de Lagrange. Les points stationnaires de \mathcal{L}(x) surviennent pour

\frac{\mathrm d\mathcal{L}(x)}{\mathrm dx} = 0 \Longrightarrow 2x^T\Sigma - 2\lambda x^T = 0 \Longrightarrow \Sigma x = \lambda x

et  \rho(x) = \frac{x^T \Sigma x}{x^T x} = \lambda \frac{x^Tx}{x^T x} = \lambda.

Par conséquent, les vecteurs propres x1...xn de Σ ne sont autre que les points critiques du quotient de Rayleigh et les valeurs propres correspondantes λ1...λn sont les valeurs rendant ρ(x) stationnaire.

Cette propriété remarquable est à la base de l’analyse en composantes principales et des corrélations canoniques.

Méthode de Rayleigh-Ritz

La Théorie de Sturm-Liouville a trait à l’action de l’application linéaire

L(y) = \frac{1}{w(x)}\left(-\frac{\mathrm d}{\mathrm dx}\left[p(x)\frac{\mathrm dy}{\mathrm dx}\right] + q(x)y\right)

sur l’espace préhilbertien des fonctions y(x) vérifiant des conditions aux limites particulières en x = a et b, muni du produit scalaire : \langle{y_1,y_2}\rangle = \int_a^b{w(x)y_1(x)y_2(x)}\mathrm dx.

Dans ce cas, le quotient de Rayleigh est

\rho(x) =\frac{\langle{y,Ly}\rangle}{\langle{y,y}\rangle} = \frac{\int_a^b{y(x)\left(-\frac{\mathrm d}{\mathrm dx}\left[p(x)\frac{\mathrm dy}{\mathrm dx}\right] + q(x)y(x)\right)}\mathrm dx}{\int_a^b{w(x)y(x)^2}\mathrm dx}.

On le présente parfois sous une forme équivalente, obtenue en découpant l'intégrale du numérateur et en intégrant par parties :

\rho(x) =\frac{\langle{y,Ly}\rangle}{\langle{y,y}\rangle} = \frac{\int_a^b{y(x)\left(-\frac{\mathrm d}{\mathrm dx}\left[p(x)y'(x)\right]\right)}\mathrm dx + \int_a^b{q(x)y(x)^2}\mathrm dx}{\int_a^b{w(x)y(x)^2}\mathrm dx}
= \frac{-y(x)\left[p(x)y'(x)\right]|_a^b + \int_a^b{y'(x)\left[p(x)y'(x)\right]}\mathrm dx + \int_a^b{q(x)y(x)^2}\mathrm dx}{\int_a^b{w(x)y(x)^2}\mathrm dx}
= \frac{-p(x)y(x)y'(x)|_a^b + \int_a^b\left[p(x)y'(x)^2 + q(x)y(x)^2\right]\mathrm dx}{\int_a^b{w(x)y(x)^2}\mathrm dx}.

Pour déterminer une solution approchée \bar y (x) de l’équation

-\frac{\mathrm d}{\mathrm dx}\left[p(x)\frac{\mathrm dy}{\mathrm dx}\right] + q(x)y = 0

vérifiant les conditions aux limites, on choisit un certain nombre de fonctions u1,u2,...,up vérifiant elles-mêmes les conditions aux limites, et on cherche la solution approchée comme une combinaison linéaire des p modes retenus : \bar y (x) = \textstyle\sum_{i=1}^p \alpha_i u_i(x). Les coefficients inconnus αi s’obtiennent en écrivant la stationnarité du quotient de Rayleigh : \tfrac{\partial \rho}{\partial \alpha_i}=0, qui détermine p équations linéaires d'inconnues i)i = 1...p

Généralisation

On peut étendre la notion de quotient de Rayleigh à deux matrices symétriques définies positives réelles (A,B), et à un vecteur non-nul x, selon :

R(A,B; x) := \frac{x^T A x}{x^T B x}.

Ce « quotient de Rayleigh généralisé » se réduit au quotient de Rayleigh R(D,Cx) par la transformation D = C TAC − 1C est la factorisation de Cholesky de la matrice B.

Notes et références

  1. Cf. par ex. Ciarlet pp.12-13.

Voir aussi

Articles connexes

Bibliographie

  • Philippe Ciarlet, Introduction à l’analyse numérique matricielle et à l’optimisation, Masson, coll. « Math. Appl. pour la Maîtrise », 1985 (réimpr. 2001) (ISBN 2-225-68893-1) 
  • P. Lascaux, R. Théodor, Analyse numérique matricielle appliquée à l’art de l’ingénieur, vol. I, Masson, 1985 (réimpr. 2001) (ISBN 2-225-84546-8), « 1.4 - Forme hermitienne associée... » 
  • John William Strutt Rayleigh, The theory of Sound, McMillan Co., 1877 (réimpr. 1945), 2 vol. (ISBN 0-486-60292-3), « IV-Vibrating systems in general », p. 106-129 

Wikimedia Foundation. 2010.

Contenu soumis à la licence CC-BY-SA. Source : Article Quotient de Rayleigh de Wikipédia en français (auteurs)

Игры ⚽ Поможем решить контрольную работу

Regardez d'autres dictionnaires:

  • Rayleigh — Cette page d’homonymie répertorie les différents sujets et articles partageant un même nom. John William Strutt Rayleigh est un physicien anglais qui a donné son nom : au Rayleigh : une unité de mesure d intensité lumineuse, à la… …   Wikipédia en Français

  • Rayleigh quotient iteration — is an eigenvalue algorithm which extends the idea of the inverse iteration by using the Rayleigh quotient to obtain increasingly accurate eigenvalue estimates.Rayleigh quotient iteration is an iterative method, that is, it must be repeated until… …   Wikipedia

  • Rayleigh — may refer to:*Rayleigh (unit), named after the son of Lord Rayleigh *Rayleigh criterion in Angular resolution *Rayleigh distribution *Rayleigh fading *Rayleigh law on low field magnetization *Rayleigh length *Rayleigh number *Rayleigh quotient… …   Wikipedia

  • Rayleigh — bezeichnet einen Ort in der englischen Grafschaft Essex, siehe Rayleigh (Essex) den Adelstitel Baron Rayleigh, darunter John William Strutt, 3. Baron Rayleigh, britischer Physiker und Nobelpreisträger, Entdecker der Duplex Theorie Siehe auch:… …   Deutsch Wikipedia

  • Rayleigh quotient — In mathematics, for a given complex Hermitian matrix A and nonzero vector x, the Rayleigh quotient R(A, x) is defined as::{x^{*} A x over x^{*} x}For real matrices and vectors, the condition of being Hermitian reduces to that of being symmetric,… …   Wikipedia

  • Rayleigh-Koeffizient — Der Rayleigh Quotient des n Vektors x zur quadratischen Matrix A ist die Zahl mit . Der Rayleigh Quotient hat eine enge Beziehung zu den Eigenwerten von A. Ist x ein Eigenvektor der Matrix A und λ der zugehörige Eigenwert, dann gilt …   Deutsch Wikipedia

  • Rayleigh-Quotient — Der Rayleigh Quotient eines Vektors x zur quadratischen Matrix A ist die Zahl mit . Hierbei bezeichnet das hermitesch Transponierte von x. Der Rayleigh Quotient hat eine enge Beziehung zu den Eigenwerten von A. Ist x …   Deutsch Wikipedia

  • John William Strutt Rayleigh — John W. Rayleigh Naissance 12 novembre 1842 Landford Grove (Essex) (Royaume Uni) Décès 30 juin 1919 Witham (Essex) (Royaume Uni) …   Wikipédia en Français

  • John Rayleigh — John William Strutt Rayleigh John W. Rayleigh Naissance 12 novembre 1842 Landford Grove (Essex) (Royaume Uni) Décès …   Wikipédia en Français

  • Lord Rayleigh — John William Strutt Rayleigh John W. Rayleigh Naissance 12 novembre 1842 Landford Grove (Essex) (Royaume Uni) Décès …   Wikipédia en Français

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”