Rayleigh quotient

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, and the conjugate transpose x^{*} to the usual transpose x'. Note that R(A, c . x) = R(A,x) for any real scalar c. Recall that a Hermitian (or real symmetric) matrix has real eigenvalues. It can be shown that the Rayleigh quotient reaches its minimum value lambda_{operatorname{min (the smallest eigenvalue of A) when x is v_{operatorname{min (the corresponding eigenvector). Similarly, R(A, x) leq lambda_{operatorname{max and R(A, v_{operatorname{max) = lambda_{operatorname{max. The Rayleigh quotient is used in Min-max theorem to get exact values of all eigenvalues. It is also used in eigenvalue algorithms to obtain an eigenvalue approximation from an eigenvector approximation. Specifically, this is the basis for Rayleigh quotient iteration.

pecial case of covariance matrices

A covariance matrix Sigma can be represented as the product A' A. Its eigenvalues are positive:

:Sigma v_i = lambda _i v_i

:A' A v_i = lambda _i v_i

:v_i' A' A v_i = v_i' lambda _i v_i

: left| A v_i ight|^2 = lambda _i left| v_i ight|^2

: lambda _i = frac{left| A v_i ight|^2}{left| v_i ight|^2} geq 0

The eigenvectors are orthogonal to one another:

:A' A v_i = lambda _i v_i

:v _j' A' A v_i = lambda _i v_j' v_i

:(A' A v_j )' v_i = lambda _i v_j' v_i

:lambda _j v_j ' v_i = lambda _i v_j' v_i

:(lambda _j - lambda _i) v_j ' v_i = 0

:v_j ' v_i = 0 (different eigenvalues, in case of multiplicity, the basis can be orthogonalized)

The Rayleigh quotient can be expressed as a function of the eigenvalues by decomposing any vector x on the basis of eigenvectors::x = sum _{i=1} ^n alpha _i v_i

: ho = frac{x' A' A x}{x' x}

: ho = 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)}

Which, by orthogonality of the eigenvectors, becomes:

: ho = frac{sum _{i=1} ^n alpha _i ^2 lambda _i}{sum _{i=1} ^n alpha _i ^2}

If a vector x maximizes ho, then any vector k . x (for k e 0) also maximizes it, one can reduce to the Lagrange problem of maximizing sum _{i=1} ^n alpha _i ^2 lambda _i under the constraint that sum _{i=1} ^n alpha _i ^2 = 1.

Since all the eigenvalues are non-negative, the problem is convex and the maximum occurs on the edges of the domain, namely when alpha _1 = 1 and forall i > 1, alpha _i = 0 (when the eigenvalues are ordered in decreasing magnitude).

Alternatively, this result can be arrived at by the method of Lagrange multipliers. The problem is to find the critical points of the function

: ho(x) = x^TSigma x , subject to the constraint |x|^2 = x^Tx = 1.I.e. to find the critical points of :mathcal{L}(x) = x^TSigma x -lambda (x^Tx - 1) , where lambda is a Lagrange multiplier. The stationary points of mathcal{L}(x) occur at

:frac{dmathcal{L}(x)}{dx} = 0 : herefore 2x^TSigma - 2lambda x^T = 0 : herefore Sigma x = lambda x and ho(x) = frac{x^T Sigma x}{x^T x} = lambda frac{x^Tx}{x^T x} = lambda

Therefore, the eigenvectors x_1 ldots x_n of Sigma are the critical points of the Raleigh Quotient and their corresponding eigenvalues lambda_1 ldots lambda_n are the stationary values of ho(x).

This property is the basis for principal components analysis and canonical correlation.

Use in Sturm-Liouville theory

Sturm-Liouville theory concerns the action of the linear operator:L(y) = frac{1}{w(x)}left(-frac{d}{dx}left [p(x)frac{dy}{dx} ight] + q(x)y ight)on the inner product space defined by:langle{y_1,y_2} angle = int_a^b{w(x)y_1(x)y_2(x)}dxof functions satisfying some specified boundary conditions at "a" and "b". In this case the Rayleigh quotient is:frac{langle{y,Ly} angle}{langle{y,y} angle} = frac{int_a^b{y(x)left(-frac{d}{dx}left [p(x)frac{dy}{dx} ight] + q(x)y(x) ight)}dx}{int_a^b{w(x)y(x)^2}dx}This is sometimes presented in an equivalent form, obtained by separating the integral in the numerator and using integration by parts::frac{langle{y,Ly} angle}{langle{y,y} angle} = frac{int_a^b{y(x)left(-frac{d}{dx}left [p(x)y'(x) ight] ight)}dx + int_a^b{q(x)y(x)^2}dx}{int_a^b{w(x)y(x)^2}dx}:= frac{-y(x)left [p(x)y'(x) ight] |_a^b + int_a^b{y'(x)left [p(x)y'(x) ight] }dx + int_a^b{q(x)y(x)^2}dx}{int_a^b{w(x)y(x)^2}dx}:= frac{-p(x)y(x)y'(x)|_a^b + int_a^bleft [p(x)y'(x)^2 + q(x)y(x)^2 ight] dx}{int_a^b{w(x)y(x)^2}dx}

ee also

* Field of values


Wikimedia Foundation. 2010.

Игры ⚽ Нужно сделать НИР?

Look at other dictionaries:

  • 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

  • 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

  • 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… …   Wikipédia en Français

  • 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 — 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

  • John William Strutt, 3. Lord Rayleigh — John William Strutt, 3. Baron Rayleigh, (* 12. November 1842 in Langford Grove, Meldon, England; † 30. Juni 1919 in Terlins Place bei Witham, England) war ein englischer Physiker. Er erhielt 1904 den Nobelpreis für Physik. John William Strutt, 3 …   Deutsch Wikipedia

  • John William Strutt Rayleigh — John William Strutt, 3. Baron Rayleigh, (* 12. November 1842 in Langford Grove, Meldon, England; † 30. Juni 1919 in Terlins Place bei Witham, England) war ein englischer Physiker. Er erhielt 1904 den Nobelpreis für Physik. John William Strutt, 3 …   Deutsch Wikipedia

  • Lord Rayleigh — John William Strutt, 3. Baron Rayleigh, (* 12. November 1842 in Langford Grove, Meldon, England; † 30. Juni 1919 in Terlins Place bei Witham, England) war ein englischer Physiker. Er erhielt 1904 den Nobelpreis für Physik. John William Strutt, 3 …   Deutsch Wikipedia

Share the article and excerpts

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