Фактор Рэлея

редактировать

конструкция для эрмитовых матриц

В математике фактор Рэлея () для данной комплексной эрмитовой матрицы M и ненулевого вектора x определяется как:

R (M, x) = х * М хх * х. {\ displaystyle R (M, x) = {x ^ {*} Mx \ over x ^ {*} x}.}{\ displaystyle R (M, x) = {x ^ {*} Mx \ over x ^ {*} x}.}

Для вещественных матриц и векторов условие эрмитовости сводится к тому, чтобы быть симметричный, а сопряжение транспонирует x ∗ {\ displaystyle x ^ {*}}x ^ {*} в обычное transpose x ′ {\ displaystyle x '}x'. Обратите внимание, что R (M, c x) = R (M, x) {\ displaystyle R (M, cx) = R (M, x)}R (M, cx) = R (M, x) для любого ненулевого скаляра c. Напомним, что эрмитова (или вещественная симметричная) матрица диагонализуема только с действительными собственными значениями. Можно показать, что для данной матрицы коэффициент Рэлея достигает минимального значения λ min {\ displaystyle \ lambda _ {\ min}}\ lambda_ \ min (наименьшее собственное значение of M), когда x равен v min {\ displaystyle v _ {\ min}}v_ \ min (соответствующий собственный вектор ). Аналогично, R (M, x) ≤ λ max {\ displaystyle R (M, x) \ leq \ lambda _ {\ max}}R (M, x) \ leq \ lambda_ \ max и R (M, v max) = λ max {\ displaystyle R (M, v _ {\ max}) = \ lambda _ {\ max}}R (M, v_ \ max) = \ lambda_ \ max .

Коэффициент Рэлея используется в теореме min-max для получения точных значений все собственные значения. Он также используется в алгоритмах собственных значений (таких как итерация коэффициента Рэлея ) для получения аппроксимации собственного значения из аппроксимации собственного вектора.

Диапазон отношения Рэлея (для любой матрицы, не обязательно эрмитовой) называется числовым диапазоном и содержит его спектр. Когда матрица эрмитова, числовой диапазон равен спектральной норме. Еще в функциональном анализе λ max {\ displaystyle \ lambda _ {\ max}}\ lambda_ \ max известен как спектральный радиус. В контексте C * -алгебр или алгебраической квантовой механики функция, которая связывает с M фактор Рэлея – Ритца R (M, x) для фиксированного x и M, изменяющегося через алгебру, будет называться «векторным состоянием» алгебра.

В квантовой механике коэффициент Рэлея дает математическое ожидание наблюдаемой, соответствующей оператору M для системы, состояние которой задается x.

Если мы зафиксируем комплексную матрицу M, то результирующая карта отношений Рэлея (рассматриваемая как функция от x) полностью определяет M через поляризационное тождество ; действительно, это остается верным, даже если мы позволим M быть неэрмитовым. (Однако, если мы ограничим поле скаляров действительными числами, то фактор Рэлея определяет только симметричную часть M.)

Содержание
  • 1 Границы для эрмитова M { \ displaystyle M}M
  • 2 Частный случай ковариационных матриц
    • 2.1 Формулировка с использованием множителей Лагранжа
  • 3 Использование в теории Штурма – Лиувилля
  • 4 Обобщения
  • 5 См. также
  • 6 Ссылки
  • 7 Дополнительная литература
Границы эрмитова M {\ displaystyle M}M

Как указано во введении, для любого вектора x R (M, x) ∈ [λ min, λ max ] {\ displaystyle R (M, x) \ in \ left [\ lambda _ {\ min}, \ lambda _ {\ max} \ right]}R (M, x) \ in \ left [\ lambda _ {\ min}, \ lambda _ {\ max} \ right] , где λ min, λ max {\ displaystyle \ lambda _ {\ min}, \ lambda _ {\ max}}{\ displaystyle \ lambda _ {\ min}, \ lambda _ {\ max}} - соответственно наименьшее и наибольшее собственные значения M {\ displaystyle M}M . Это происходит сразу после наблюдения, что фактор Рэлея представляет собой средневзвешенное значение собственных значений M:

R (M, x) = x ∗ M xx ∗ x = ∑ i = 1 n λ iyi 2 ∑ i = 1 nyi 2 { \ Displaystyle R (M, x) = {x ^ {*} Mx \ over x ^ {*} x} = {\ frac {\ sum _ {i = 1} ^ {n} \ lambda _ {i} y_ { i} ^ {2}} {\ sum _ {i = 1} ^ {n} y_ {i} ^ {2}}}}R (M, x) = {x ^ {{*}} Mx \ over x ^ {{*}} x} = {\ frac {\ sum _ {{i = 1}} ^ {n} \ lambda _ {i} y_ {i} ^ {2}} {\ sum _ {{i = 1}} ^ {n} y_ {i} ^ {2}}}

где (λ i, vi) {\ displaystyle (\ lambda _ {i}, v_ {i})}(\ lambda _ {i}, v_ {i}) - это i {\ displaystyle i}i -я собственная пара после ортонормировки и yi = vi ∗ x {\ displaystyle y_ {i} = v_ {i} ^ {*} x}y_ {i} = v_ {i} ^ {*} x - это i {\ displaystyle i}i -я координата x в собственном базисе. Затем легко проверить, что границы достигаются на соответствующих собственных векторах v min, v max {\ displaystyle v _ {\ min}, v _ {\ max}}v _ {\ min}, v _ {\ max} .

Тот факт, что частное является средневзвешенным собственных значений можно использовать для идентификации второго, третьего,... наибольших собственных значений. Пусть λ max = λ 1 ≥ λ 2 ≥ ⋯ ≥ λ N = λ min {\ displaystyle \ lambda _ {\ max} = \ lambda _ {1} \ geq \ lambda _ {2} \ geq \ cdots \ geq \ lambda _ {n} = \ lambda _ {\ min}}{\ displaystyle \ lambda _ {\ max} = \ lambda _ {1} \ geq \ lambda _ {2} \ geq \ cdots \ geq \ lambda _ {n} = \ lambda _ {\ min}} - собственные значения в порядке убывания. Если n = 2 {\ displaystyle n = 2}п = 2 и x {\ displaystyle x}x ограничено, чтобы быть ортогональным v 1 {\ displaystyle v_ {1}}v_ {1} , в этом случае y 1 = v 1 ∗ x = 0 {\ displaystyle y_ {1} = v_ {1} ^ {*} x = 0}y_ {1} = v_ {1} ^ {*} x = 0 , тогда R (M, x) {\ displaystyle R (M, x)}R (M, x) имеет максимальное значение λ 2 {\ displaystyle \ lambda _ {2}}\ lambda _ {2} , который достигается, когда x = v 2 {\ displaystyle x = v_ {2}}x = v_ {2 } .

Особый случай ковариационных матриц

Эмпирическая ковариационная матрица M {\ displaystyle M}M может быть представлен как продукт A ′ A {\ displaystyle A'A}{\displaystyle A'A}матрицы данных A {\ displaystyle A}A , предварительно умноженное на его транспонирование A ′ {\ displaystyle A '}A'. Будучи положительной полуопределенной матрицей, M {\ displaystyle M}M имеет неотрицательные собственные значения и ортогональные (или ортогонализуемые) собственные векторы, что можно продемонстрировать следующим образом.

Во-первых, собственные значения λ i {\ displaystyle \ lambda _ {i}}\ lambda _ {i} неотрицательны:

M vi = A ′ A vi = λ ivi {\ displaystyle Mv_ {i} = A'Av_ {i} = \ lambda _ {i} v_ {i}}M v_i = A' A v_i = \lambda_i v_i
⇒ vi ′ A ′ A vi = vi ′ λ ivi {\ displaystyle \ Rightarrow v_ {i} 'A'Av_ {i} = v_ {i}' \ lambda _ {i} v_ {i}}\Rightarrow v_i' A' A v_i = v_i' \lambda_i v_i
⇒ ‖ A vi ‖ 2 = λ я ‖ vi ‖ 2 {\ displaystyle \ Rightarrow \ left \ | Av_ {i} \ right \ | ^ {2} = \ lambda _ {i} \ left \ | v_ {i} \ right \ | ^ {2}}\ Rightarrow \ left \ | A v_i \ right \ | ^ 2 = \ lambda_i \ left \ | v_i \ right \ | ^ 2
⇒ λ i = ‖ A vi ‖ 2 ‖ vi ‖ 2 ≥ 0. {\ displaystyle \ Rightarrow \ lambda _ {i} = {\ frac {\ left \ | Av_ {i} \ right \ | ^ {2}} {\ left \ | v_ {i} \ right \ | ^ {2}}} \ geq 0.}\ Rightarrow \ lambda_i = \ frac {\ left \ | A v_i \ right \ | ^ 2} {\ left \ | v_i \ right \ | ^ 2} \ geq 0.

Во-вторых, собственные векторы vi {\ displaystyle v_ {i}}v_ {i} ортогональны друг другу:

M vi = λ ivi ⇒ vj ′ M vi = vj ′ λ ivi ⇒ (M vj) ′ vi = λ jvj ′ vi ⇒ λ jvj ′ vi = λ ivj ′ vi ⇒ (λ j - λ i) vj ′ vi = 0 ⇒ vj ′ vi = 0 {\ displaystyle {\ begin {выровнено} Mv_ {i} = \ lambda _ {i} v_ {i} \\ \ Rightarrow v_ {j} 'Mv_ {i} = v_ {j}' \ lambda _ { i} v_ {i} \\ \ Rightarrow \ left (Mv_ {j} \ right) 'v_ {i} = \ lambda _ {j } v_ {j} 'v_ {i} \\ \ Rightarrow \ lambda _ {j} v_ {j}' v_ {i} = \ lambda _ {i} v_ {j} 'v_ {i} \\ \ Rightarrow \ left (\ lambda _ {j} - \ lambda _ {i} \ right) v_ {j} 'v_ {i} = 0 \\ \ Rightarrow v_ {j}' v_ {i} = 0 \ end { выровнен}}}{\displaystyle {\begin{aligned}Mv_{i}=\lambda _{i}v_{i}\\\Rightarrow v_{j}'Mv_{i}=v_{j}'\lambda _{i}v_{i}\\\Rightarrow \left(Mv_{j}\right)'v_{i}=\lambda _{j}v_{j}'v_{i}\\\Rightarrow \lambda _{j}v_{j}'v_{i}=\lambda _{i}v_{j}'v_{i}\\\Rightarrow \left(\lambda _{j}-\lambda _{i}\right)v_{j}'v_{i}=0\\\Rightarrow v_{j}'v_{i}=0\end{aligned}}}

если собственные значения разные - в случае множественности базис может быть ортогонализирован.

Теперь, чтобы установить, что коэффициент Рэлея максимизируется собственным вектором с наибольшим собственным значением, рассмотрим разложение произвольного вектора x {\ displaystyle x}x на основе собственных векторов vi {\ displaystyle v_ {i}}v_ {i} :

x = ∑ i = 1 n α ivi, {\ displaystyle x = \ sum _ {i = 1} ^ {n} \ alpha _ {i} v_ {i },}x = \ sum _ {i = 1} ^ n \ alpha _i v_i,

где

α i = x ′ vivi ′ vi = ⟨x, vi⟩ ‖ vi ‖ 2 {\ displaystyle \ alpha _ {i} = {\ frac {x'v_ {i}} { {v_ {i}} '{v_ {i}}}} = {\ frac {\ langle x, v_ {i} \ rangle} {\ left \ | v_ {i} \ right \ | ^ {2}}} }{\displaystyle \alpha _{i}={\frac {x'v_{i}}{{v_{i}}'{v_{i}}}}={\frac {\langle x,v_{i}\rangle }{\left\|v_{i}\right\|^{2}}}}

- координата x {\ displaystyle x}x , ортогонально проецируемая на vi {\ displaystyle v_ {i}}v_ {i} . Следовательно, имеем:

R (M, x) = x ′ A ′ A xx ′ x = (∑ j = 1 n α jvj) ′ (A ′ A) (∑ i = 1 n α ivi) (∑ j = 1 n α jvj) ′ (∑ i = 1 n α ivi) = (∑ j = 1 n α jvj) ′ (∑ i = 1 n α i (A ′ A) vi) (∑ i = 1 n α я 2 vi ′ vi) знак равно (∑ j = 1 n α jvj) ′ (∑ i = 1 n α i λ ivi) (∑ i = 1 n α i 2 ‖ vi ‖ 2) {\ displaystyle {\ begin {выровнено } R (M, x) = {\ frac {x'A'Ax} {x'x}} \\ = {\ frac {{\ Bigl (} \ sum _ {j = 1} ^ {n} \ alpha _ {j} v_ {j} {\ Bigr)} '\ left (A'A \ right) {\ Bigl (} \ sum _ {i = 1} ^ {n} \ alpha _ {i} v_ { i} {\ Bigr)}} {{\ bigl (} \ sum _ {j = 1} ^ {n} \ alpha _ {j} v_ {j} {\ Bigr)} '{\ Bigl (} \ sum _ {i = 1} ^ {n} \ alpha _ {i} v_ {i} {\ Bigr)}}} \\ = {\ frac {{\ Bigl (} \ sum _ {j = 1} ^ {n } \ alpha _ {j} v_ {j} {\ Bigr)} '{\ Bigl (} \ sum _ {i = 1} ^ {n} \ alpha _ {i} (A'A) v_ {i} { \ Bigr)}} {{\ bigl (} \ sum _ {i = 1} ^ {n} \ alpha _ {i} ^ {2} {v_ {i}} '{v_ {i}} {\ Bigr) }}} \\ = {\ frac {{\ Bigl (} \ sum _ {j = 1} ^ {n} \ alpha _ {j} v_ {j} {\ Bigr)} '{\ Bigl (} \ сумма _ {i = 1} ^ {n} \ alpha _ {i} \ lambda _ {i} v_ {i} {\ Bigr)}} {{\ bigl (} \ sum _ {i = 1} ^ {n } \ alpha _ {i} ^ {2} \ | {v_ {i}} \ | ^ {2} {\ Bigr)}}} \ end {a ligned}}}{\displaystyle {\begin{aligned}R(M,x)={\frac {x'A'Ax}{x'x}}\\={\frac {{\Bigl (}\sum _{j=1}^{n}\alpha _{j}v_{j}{\Bigr)}'\left(A'A\right){\Bigl (}\sum _{i=1}^{n}\alpha _{i}v_{i}{\Bigr)}}{{\bigl (}\sum _{j=1}^{n}\alpha _{j}v_{j}{\Bigr)}'{\Bigl (}\sum _{i=1}^{n}\alpha _{i}v_{i}{\Bigr)}}}\\={\frac {{\Bigl (}\sum _{j=1}^{n}\alpha _{j}v_{j}{\Bigr)}'{\Bigl (}\sum _{i=1}^{n}\alpha _{i}(A'A)v_{i}{\Bigr)}}{{\bigl (}\sum _{i=1}^{n}\alpha _{i}^{2}{v_{i}}'{v_{i}}{\Bigr)}}}\\={\frac {{\Bigl (}\sum _{j=1}^{n}\alpha _{j}v_{j}{\Bigr)}'{\Bigl (}\sum _{i=1}^{n}\alpha _{i}\lambda _{i}v_{i}{\Bigr)}}{{\bigl (}\sum _{i=1}^{n}\alpha _{i}^{2}\|{v_{i}}\|^{2}{\Bigr)}}}\end{aligned}}}

который в силу ортонормированности собственных векторов становится:

R (M, x) = ∑ i = 1 n α i 2 λ i ∑ i = 1 n α i 2 знак равно ∑ я знак равно 1 N λ я (x ′ vi) 2 (x ′ x) (vi ′ vi) {\ displaystyle R (M, x) = {\ frac {\ sum _ {i = 1} ^ {n } \ alpha _ {i} ^ {2} \ lambda _ {i}} {\ sum _ {i = 1} ^ {n} \ alpha _ {i} ^ {2}}} = \ sum _ {i = 1} ^ {n} \ lambda _ {i} {\ frac {(x'v_ {i}) ^ {2}} {(x'x) (v_ {i} 'v_ {i})}}}R(M,x) = \frac{\sum _{i=1} ^n \alpha_i^2 \lambda _i}{\sum _{i=1} ^n \alpha_i^2} = \sum_{i=1}^n \lambda_i \frac{(x'v_i)^2}{ (x'x)( v_i' v_i)}

Последнее представление устанавливает, что фактор Рэлея представляет собой сумму квадратов косинусов углов, образованных вектором x {\ displaystyle x}x и каждым собственным вектором vi {\ displaystyle v_ {i}}v_ {i} , взвешенный по соответствующим собственным значениям.

Если вектор x {\ displaystyle x}x максимизирует R (M, x) {\ displaystyle R (M, x)}R (M, x) , то любое ненулевое скалярное кратное kx {\ displaystyle kx}kx также максимизирует R {\ displaystyle R}R , поэтому проблема может быть сведена к Задача Лагранжа максимизации ∑ i = 1 n α i 2 λ i {\ displaystyle \ sum _ {i = 1} ^ {n} \ alpha _ {i} ^ {2} \ lambda _ {i}}\ sum _ {i = 1} ^ n \ alpha_i ^ 2 \ lambda _i при ограничении ∑ i = 1 n α i 2 = 1 {\ displaystyle \ sum _ {i = 1} ^ {n} \ alpha _ {i} ^ { 2} = 1}\ sum _ {i = 1} ^ n \ alpha _i ^ 2 = 1 .

Определите: β i = α i 2 {\ displaystyle \ beta _ {i} = \ alpha _ {i} ^ {2}}{\ displaystyle \ beta _ {i} = \ alpha _ {i} ^ {2}} . Тогда это становится линейной программой, которая всегда достигает своего максимума в одном из углов области. Максимальная точка будет иметь α 1 = ± 1 {\ displaystyle \ alpha _ {1} = \ pm 1}\ alpha_1 = \ pm 1 и α i = 0 {\ displaystyle \ alpha _ {i} = 0}\ альфа _i = 0 для всех i>1 {\ displaystyle i>1}i>1 (когда собственные значения упорядочены по убывающей величине).

Таким образом, собственный вектор Рэлея максимизируется с помощью наибольшее собственное значение.

Формулировка с использованием множителей Лагранжа

В качестве альтернативы этот результат может быть получен методом множителей Лагранжа. Первая часть - показать, что частное постоянна при масштабировании x → cx {\ displaystyle x \ to cx}{\ displaystyle x \ to cx} , где c {\ displaystyle c}c - скаляр

R (M, сх) знак равно (сх) * М сх (сх) * сх = с * сс * сх * М хх * х = р (М, х). {\ Displaystyle R (M, сх) = {\ гидроразрыва {(сх) ^ {*} Mcx} {(cx) ^ {*} cx}} = {\ frac {c ^ {*} c} {c ^ {*} c}} {\ frac {x ^ {*} Mx} {x ^ {*} x}} = R (M, x).}{\ displaystyle R (M, cx) = {\ frac {(cx) ^ { *} Mcx} {(cx) ^ {*} cx}} = {\ frac {c ^ {*} c} {c ^ {*} c}} {\ frac {x ^ {*} Mx} {x ^ {*} x}} = R (M, x).}

Из-за этой инвариантности достаточно изучить частный случай ‖ x ‖ 2 = Икс T Икс знак равно 1 {\ Displaystyle \ | х \ | ^ {2} = х ^ {T} х = 1}{\ displaystyle \ | x \ | ^ {2} = x ^ {T} x = 1} . Тогда задача состоит в том, чтобы найти критические точки функции

R (M, x) = x TM x {\ displaystyle R (M, x) = x ^ {T} Mx}R (M, x) = x ^ TM x ,

с учетом ограничения ‖ x ‖ 2 = x T x = 1. {\ displaystyle \ | x \ | ^ {2} = x ^ {T} x = 1.}\ | x \ | ^ 2 = x ^ Tx = 1. В другом словами, это найти критические точки

L (x) = x TM x - λ (x T x - 1), {\ displaystyle {\ mathcal {L}} (x) = x ^ {T} Mx- \ lambda \ left (x ^ {T} x-1 \ right),}\ mathcal {L} (x) = x ^ TM x - \ лямбда \ left (x ^ Tx - 1 \ right),

, где λ {\ displaystyle \ lambda}\ lambda - множитель Лагранжа. Стационарные точки L (x) {\ displaystyle {\ mathcal {L}} (x)}\ mathcal {L} (x) находятся в

d L (x) dx = 0 ⇒ 2 x TM - 2 λ x T знак равно 0 ⇒ 2 M x - 2 λ x = 0 ⇒ M x = λ x {\ displaystyle {\ begin {align} {\ frac {d {\ mathcal {L}} (x)} {dx} } = 0 \\ \ Rightarrow 2x ^ {T} M-2 \ lambda x ^ {T} = 0 \\ \ Rightarrow 2Mx-2 \ lambda x = 0 \\ \ Rightarrow Mx ​​= \ lambda x \ end {выровнено}}}{\ displaystyle {\ begin {align} {\ frac {d {\ mathcal {L}} (x)} {dx}} = 0 \\ \ Стрелка вправо 2x ^ {T} M-2 \ lambda x ^ {T} = 0 \\ \ Rightarrow 2Mx-2 \ lambda x = 0 \\ \ Rightarrow Mx ​​= \ lambda x \ end {align}}}

и

∴ R (M, x) = x TM xx T x = λ x T xx T x = λ. {\ displaystyle \, следовательно, R (M, x) = {\ frac {x ^ {T} Mx} {x ^ {T} x}} = \ lambda {\ frac {x ^ {T} x} {x ^ { T} x}} = \ lambda.}{\ displaystyle \, следовательно, R (M, x) = {\ frac { x ^ {T} Mx} {x ^ {T} x}} = \ lambda {\ frac {x ^ {T} x} {x ^ {T} x}} = \ lambda.}

Следовательно, собственные векторы x 1,…, xn {\ displaystyle x_ {1}, \ ldots, x_ {n}}x_ {1}, \ ldots, x_ {n} of M {\ displaystyle M}M - критические точки отношения Рэлея и соответствующие им собственные значения λ 1,…, λ n {\ displaystyle \ lambda _ {1}, \ ldots, \ лямбда _ {n}}\ lambda _ {1}, \ ldots, \ lambda _ {n} - это стационарные значения L {\ displaystyle {\ mathcal {L}}}{\ mathcal {L}} . Это свойство является основой для анализа главных компонент и канонической корреляции.

Использование в теории Штурма – Лиувилля

теория Штурма – Лиувилля касается действия линейного оператора

L (y) = 1 вес (Икс) (- ddx [p (x) dydx] + q (x) y) {\ displaystyle L (y) = {\ frac {1} {w (x)} } \ left (- {\ frac {d} {dx}} \ left [p (x) {\ frac {dy} {dx}} \ right] + q (x) y \ right)}L (y) = \ frac {1} {w (x)} \ left (- \ frac {d} {dx} \ left [p (x) \ frac {dy} {dx} \ right] + q (x) y \ right)

на внутреннее пространство продукта определяется как

⟨y 1, y 2⟩ = ∫ abw (x) y 1 (x) y 2 (x) dx {\ displaystyle \ langle {y_ {1}, y_ {2}} \ rangle = \ int _ {a} ^ {b} w (x) y_ {1} (x) y_ {2} (x) \, dx}\ langle {y_1, y_2} \ rangle = \ int_a ^ bw (x) y_1 (x) y_2 (x) \, dx

функций, удовлетворяющих некоторому заданному граничные условия в точках a и b. В этом случае фактор Рэлея равен

⟨y, L y⟩ ⟨y, y⟩ = ∫ aby (x) (- ddx [p (x) dydx] + q (x) y (x)) dx ∫ abw (х) у (х) 2 дх. {\ displaystyle {\ frac {\ langle {y, Ly} \ rangle} {\ langle {y, y} \ rangle}} = {\ frac {\ int _ {a} ^ {b} y (x) \ left (- {\ frac {d} {dx}} \ left [p (x) {\ frac {dy} {dx}} \ right] + q (x) y (x) \ right) dx} {\ int _ {a} ^ {b} {w (x) y (x) ^ {2}} dx}}.}\ frac {\ langle {y, Ly} \ rangle} { \ langle {y, y} \ rangle} = \ frac {\ int_a ^ by (x) \ left (- \ frac {d} {dx} \ left [p (x) \ frac {dy} {dx} \ right ] + q (x) y (x) \ right) dx} {\ int_a ^ b {w (x) y (x) ^ 2} dx}.

Иногда это представляется в эквивалентной форме, полученной путем разделения интеграла в числителе и использования интегрирование по частям :

⟨y, L y⟩ ⟨y, y⟩ = {∫ aby (x) (- ddx [p (x) y ′ (x)]) dx} + {∫ abq (x) y (x) 2 dx} ∫ abw (x) y (x) 2 dx = {- y (x) [p (x) y ′ (x)] | ab} + {∫ aby ′ (x) [p (x) y ′ (x)] dx} + {∫ abq (x) y (x) 2 dx} ∫ abw (x) y (x) 2 dx = { - p (x) y (x) y ′ (x) | a b} + {∫ a b [p (x) y ′ (x) 2 + q (x) y (x) 2] d x} ∫ a b w (x) y (x) 2 d x. {\ displaystyle {\ begin {align} {\ frac {\ langle {y, Ly} \ rangle} {\ langle {y, y} \ rangle}} = {\ frac {\ left \ {\ int _ {a } ^ {b} y (x) \ left (- {\ frac {d} {dx}} \ left [p (x) y '(x) \ right] \ right) dx \ right \} + \ left \ {\ int _ {a} ^ {b} {q (x) y (x) ^ {2}} \, dx \ right \}} {\ int _ {a} ^ {b} {w (x) y (x) ^ {2}} \, dx}} \\ = {\ frac {\ left \ {\ left.-y (x) \ left [p (x) y '(x) \ right] \ right | _ {a} ^ {b} \ right \} + \ left \ {\ int _ {a} ^ {b} y '(x) \ left [p (x) y' (x) \ right] \, dx \ right \} + \ left \ {\ int _ {a} ^ {b} {q (x) y (x) ^ {2}} \, dx \ right \}} {\ int _ {a} ^ {b} w (x) y (x) ^ {2} \, dx}} \\ = {\ frac {\ left \ {\ left.-p (x) y (x) y '(x) \ right | _ {a} ^ {b} \ right \} + \ left \ {\ int _ {a} ^ {b} \ left [p (x) y '(x) ^ {2} + q (x) y (x) ^ {2} \ right] \, dx \ right \}} {\ int _ {a} ^ {b} {w (x) y (x) ^ {2}} \, dx}}. \ end {align}}}{\displaystyle {\begin{aligned}{\frac {\langle {y,Ly}\rangle }{\langle {y,y}\rangle }}={\frac {\left\{\int _{a}^{b}y(x)\left(-{\frac {d}{dx}}\left[p(x)y'(x)\right]\right)dx\right\}+\left\{\int _{a}^{b}{q(x)y(x)^{2}}\,dx\right\}}{\int _{a}^{b}{w(x)y(x)^{2}}\,dx}}\\={\frac {\left\{\left.-y(x)\left[p(x)y'(x)\right]\right|_{a}^{b}\right\}+\left\{\int _{a}^{b}y'(x)\left[p(x)y'(x)\right]\,dx\right\}+\left\{\int _{a}^{b}{q(x)y(x)^{2}}\,dx\right\}}{\int _{a}^{b}w(x)y(x)^{2}\,dx}}\\={\frac {\left\{\left.-p(x)y(x)y'(x)\right|_{a}^{b}\right\}+\left\{\int _{a}^{b}\left[p(x)y'(x)^{2}+q(x)y(x)^{2}\right]\,dx\right\}}{\int _{a}^{b}{w(x)y(x)^{2}}\,dx}}.\end{aligned}}}
Обобщения
  1. Для данной пары (A, B) матриц и данного ненулевого вектора x обобщенное частное Рэлея определяется как:
    R (A, B; x): = x ∗ A xx ∗ B x. {\ displaystyle R (A, B; x): = {\ frac {x ^ {*} Ax} {x ^ {*} Bx}}.}R (A, B; x): = \ frac {x ^ * A x} {x ^ * B x}.
    Обобщенный коэффициент Рэлея можно уменьшить до коэффициента Рэлея R (D, C ∗ x) {\ displaystyle R (D, C ^ {*} x)}R (D, C ^ * x) через преобразование D = C - 1 AC ∗ - 1 {\ displaystyle D = C ^ {- 1} A {C ^ {*}} ^ {- 1}}D = C ^ {- 1} A {C ^ *} ^ {- 1} где CC ∗ {\ displaystyle CC ^ {*}}CC ^ * - это Разложение Холецкого эрмитовой положительно определенной матрицы B.
  2. Для данной пары (x, y) ненулевых векторов и данной эрмитовой матрицы H обобщенное частное Рэлея можно определить как:
    R (H; x, y): = y ∗ H xy ∗ y ⋅ x ∗ x {\ displaystyle R (H; x, y): = {\ frac {y ^ { *} Hx} {\ sqrt {y ^ {*} y \ cdot x ^ {*} x}}}}R (H; x, y): = {\ frac {y ^ {*} Hx} {\ sqrt {y ^ {*} y \ cdot x ^ {*} x} }}
    , который совпадает с R (H, x), когда x = y. В квантовой механике эта величина называется «матричным элементом» или иногда «амплитудой перехода».
См. Также
Ссылки
  1. ^Также известно как отношение Рэлея – Ритца ; назван в честь Вальтера Ритца и лорда Рэлея.
  2. ^Хорна, Р. А.; Джонсон, К. А. (1985). Матричный анализ. Издательство Кембриджского университета. С. 176–180. ISBN 0-521-30586-1.
  3. ^Парлетт Б. Н. (1998). Симметричная проблема собственных значений. Классика прикладной математики. СИАМ. ISBN 0-89871-402-8.
  4. ^Костин, Родика Д. (2013). «Среднесрочные заметки» (PDF). Математика 5102 Линейная математика в бесконечных измерениях, конспекты лекций. Университет штата Огайо.
Дополнительная литература
Последняя правка сделана 2021-06-03 09:25:27
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте