Теорема Бауэра – Фике

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

В математике теорема Бауэра – Фике является стандартным результатом в теория возмущений для собственного значения комплексной диагонализуемой матрицы. По сути, он устанавливает абсолютную верхнюю границу отклонения одного возмущенного собственного значения матрицы от правильно выбранного собственного значения точной матрицы. Неформально говоря, здесь говорится, что чувствительность собственных значений оценивается числом обусловленности матрицы собственных векторов.

Теорема была доказана Фридрихом Л. Бауэром и CT Fike в 1960 году.

Содержание
  • 1 Установка
  • 2 Теорема Бауэра – Фике
  • 3 Альтернативная формулировка
  • 4 Относительная граница
  • 5 Случай нормальных матриц
  • 6 Ссылки
Схема

Далее мы предполагаем, что:

κ p (X) = ‖ X ‖ p ‖ X - 1 ‖ p. {\ displaystyle \ kappa _ {p} (X) = \ | X \ | _ {p} \ left \ | X ^ {- 1} \ right \ | _ {p}.}\ kappa _ {p} (X) = \ | X \ | _ {p} \ left \ | X ^ {{- 1}} \ right \ | _ {p}.
Теорема Бауэра – Фике
Теорема Бауэра – Фике. Пусть μ - собственное значение оператора A + δA. Тогда существует λ ∈ Λ (A) такое, что:
| λ - μ | ≤ κ п (V) ‖ δ A ‖ p {\ displaystyle | \ lambda - \ mu | \ leq \ kappa _ {p} (V) \ | \ delta A \ | _ {p}}| \ lambda - \ mu | \ leq \ kappa _ {p} (V) \ | \ delta A \ | _ {p}

Доказательство. Можно предположить, что μ ∉ Λ (A), в противном случае возьмем λ = μ, и результат тривиально верен, поскольку κ p (V) ≥ 1. Поскольку μ является собственным значением A + δA, имеем det (A + δA - μI) = 0, поэтому

0 = det (A + δ A - μ I) = det (V - 1) det (A + δ A - μ I) det (V) = det (V - 1 (A + δ A - μ I) V) = det (V - 1 AV + V - 1 δ AV - V - 1 μ IV) = det (Λ + V - 1 δ AV - μ I) = det (Λ - μ I) det ((Λ - μ I) - 1 V - 1 δ AV + I) {\ displaystyle {\ begin {align} 0 = \ det (A + \ delta A- \ mu I) \\ = \ det (V ^ {- 1}) \ det (A + \ delta A- \ mu I) \ det (V) \\ = \ det \ left (V ^ {- 1} (A + \ delta A- \ mu I) V \ right) \\ = \ det \ left (V ^ {- 1} AV + V ^ {- 1} \ delta AV-V ^ {- 1} \ mu IV \ right) \\ = \ det \ left (\ Lambda + V ^ {- 1} \ delta AV- \ mu I \ right) \\ = \ det (\ Lambda - \ mu I) \ det \ left ((\ Lambda - \ mu I) ^ {- 1} V ^ {- 1} \ delta AV + I \ right) \\\ end {align}}}{\ begin {align} 0 = \ det (A + \ delta A- \ mu I) \\ = \ det (V ^ {{- 1}}) \ det (A + \ delta A- \ mu I) \ det (V) \\ = \ det \ left (V ^ {{- 1}} (A + \ delta A- \ mu I) V \ right) \\ = \ det \ left (V ^ {{- 1}} AV + V ^ {{- 1}} \ delta AV-V ^ {{- 1}} \ mu IV \ right) \\ = \ det \ left (\ Lambda + V ^ {{- 1}} \ delta AV- \ mu I \ right) \\ = \ det (\ Lambda - \ mu I) \ det \ left ( (\ Lambda - \ mu I) ^ {{- 1}} V ^ {{- 1}} \ delta AV + I \ right) \\\ конец {выровнено}}

Однако наше предположение, μ ∉ Λ (A), означает, что: det (Λ - μI) ≠ 0, поэтому мы можем записать:

det (( Λ - μ I) - 1 В - 1 δ AV + I) знак равно 0. {\ Displaystyle \ det \ left ((\ Lambda - \ mu I) ^ {- 1} V ^ {- 1} \ delta AV + I \ right) = 0.}\ det \ left ((\ Lambda - \ mu I) ^ {{- 1}} V ^ {{- 1}} \ delta AV + I \ right) = 0.

Это показывает, что −1 является собственным значением

(Λ - μ I) - 1 V - 1 δ AV. {\ displaystyle (\ Lambda - \ mu I) ^ {- 1} V ^ {- 1} \ delta AV.}(\ Лямбда - \ му I) ^ {{- 1}} V ^ {{- 1}} \ дельта AV.

Поскольку все p-нормы являются согласованными матричными нормами, мы имеем | λ | ≤ || A || p, где λ - собственное значение A. В данном случае это дает нам:

1 = | - 1 | ≤ ‖ (Λ - μ I) - 1 V - 1 δ AV ‖ p ≤ ‖ (Λ - μ I) - 1 ‖ p ‖ V - 1 ‖ p ‖ V ‖ p ‖ δ A ‖ p = ‖ (Λ - μ I) - 1 ‖ п κ п (V) ‖ δ A ‖ p {\ displaystyle 1 = | -1 | \ leq \ left \ | (\ Lambda - \ mu I) ^ {- 1} V ^ {- 1} \ delta AV \ right \ | _ {p} \ leq \ left \ | (\ Lambda - \ mu I) ^ {- 1} \ right \ | _ {p} \ left \ | V ^ {- 1} \ right \ | _ {p} \ | V \ | _ {p} \ | \ delta A \ | _ {p} = \ left \ | (\ Lambda - \ mu I) ^ {- 1} \ right \ | _ { p} \ \ kappa _ {p} (V) \ | \ delta A \ | _ {p}}1 = | -1 | \ leq \ left \ | (\ Lambda - \ mu I) ^ {{- 1}} V ^ { {-1}} \ delta AV \ right \ | _ {p} \ leq \ left \ | (\ Lambda - \ mu I) ^ {{- 1}} \ right \ | _ {p} \ left \ | V ^ {{- 1}} \ right \ | _ {p} \ | V \ | _ {p} \ | \ delta A \ | _ {p} = \ left \ | (\ Lambda - \ mu I) ^ { {-1}} \ right \ | _ {p} \ \ kappa _ {p} (V) \ | \ delta A \ | _ {p}

Но (Λ - μI) - диагональная матрица, p-норма которой легко вычисляется:

‖ (Λ - μ I) - 1 ‖ p = max ‖ x ‖ p ≠ 0 ‖ (Λ - μ I) - 1 x ‖ p ‖ x ‖ p = max λ ∈ Λ (A) 1 | λ - μ | = 1 min λ ∈ Λ (A) | λ - μ | {\ displaystyle \ left \ | \ left (\ Lambda - \ mu I \ right) ^ {- 1} \ right \ | _ {p} \ = \ max _ {\ | {\ boldsymbol {x}} \ | _ {p} \ neq 0} {\ frac {\ left \ | \ left (\ Lambda - \ mu I \ right) ^ {- 1} {\ boldsymbol {x}} \ right \ | _ {p}} {\ | {\ boldsymbol {x}} \ | _ {p}}} = \ max _ {\ lambda \ in \ Lambda (A)} {\ frac {1} {| \ lambda - \ mu |}} \ = { \ frac {1} {\ min _ {\ lambda \ in \ Lambda (A)} | \ lambda - \ mu |}}}{\ displaystyle \ left \ | \ left (\ Lambda - \ mu I \ right) ^ {- 1} \ right \ | _ {p} \ = \ max _ {\ | {\ boldsymbol {x}} \ | _ {p} \ neq 0} {\ frac {\ left \ | \ left (\ Lambda - \ mu I \ right) ^ {- 1} {\ boldsymbol {x}} \ right \ | _ {p}} {\ | {\ boldsymbol {x}} \ | _ { p}}} = \ max _ {\ lambda \ in \ Lambda (A)} {\ frac {1} {| \ lambda - \ mu |}} \ = {\ frac {1} {\ min _ {\ lambda \ in \ Lambda (A)} | \ lambda - \ mu |}}}

откуда:

min λ ∈ Λ (A) | λ - μ | ≤ κ p (V) ‖ δ A ‖ p. {\ displaystyle \ min _ {\ lambda \ in \ Lambda (A)} | \ lambda - \ mu | \ leq \ \ kappa _ {p} (V) \ | \ delta A \ | _ {p}.}{\ displaystyle \ min _ {\ lambda \ in \ Lambda (A)} | \ lambda - \ mu | \ leq \ \ kappa _ {p} (V) \ | \ delta A \ | _ {p}.}
Альтернативная формулировка

Теорема также может быть переформулирована, чтобы лучше подходить для численных методов. Фактически, имея дело с реальными проблемами собственной системы, можно часто иметь точную матрицу A, но знать только приблизительную пару собственное значение-собственный вектор (λ, v ) и необходимо ограничить ошибку. Следующая версия поможет вам.

Теорема Бауэра – Фике (альтернативная формулировка). Пусть (λ, v ) будет приближенной парой собственное значение-собственный вектор, и r = A v - λ v . Тогда существует λ ∈ Λ (A) такое, что:
| λ - λ a | ≤ κ п (В) ‖ р ‖ п ‖ ва ‖ п {\ displaystyle \ left | \ lambda - \ lambda ^ {a} \ right | \ leq \ kappa _ {p} (V) {\ frac {\ | { \ boldsymbol {r}} \ | _ {p}} {\ left \ | {\ boldsymbol {v}} ^ {a} \ right \ | _ {p}}}}\ left | \ lambda - \ lambda ^ {a} \ right | \ leq \ kappa _ {p} (V) {\ frac {\ | {\ boldsymbol {r}} \ | _ {p}} {\ left \ | {\ boldsymbol {v}} ^ {a} \ right \ | _ {p}}}

Доказательство. Мы можем предположим, что λ ∉ Λ (A), в противном случае возьмем λ = λ, и результат тривиально верен, поскольку κ p (V) ≥ 1. Значит (A - λI) существует, поэтому мы можем написать:

va знак равно (A - λ a I) - 1 r = V (D - λ a I) - 1 V - 1 r {\ displaystyle {\ boldsymbol {v}} ^ {a} = \ left (A- \ lambda ^ {a} I \ right) ^ {- 1} {\ boldsymbol {r}} = V \ left (D- \ lambda ^ {a} I \ right) ^ {- 1} V ^ {- 1} {\ boldsymbol {r}}}{\ boldsymbol {v}} ^ {a} = \ left (A- \ lambda ^ {a} I \ right) ^ {{- 1}} {\ boldsymbol {r}} = V \ left (D- \ lambda ^ {a} I \ right) ^ {{- 1}} V ^ {{- 1}} {\ boldsymbol {r}}

поскольку A диагонализуема; взяв p-норму обеих частей, получаем:

‖ va ‖ p = ‖ V (D - λ a I) - 1 V - 1 r ‖ p ≤ ‖ V ‖ p ‖ (D - λ a I) - 1 ‖ p ‖ V - 1 ‖ p ‖ r ‖ p = κ p (V) ‖ (D - λ a I) - 1 ‖ p ‖ r ‖ p. {\ displaystyle \ left \ | {\ boldsymbol {v}} ^ {a} \ right \ | _ {p} = \ left \ | V \ left (D- \ lambda ^ {a} I \ right) ^ {- 1} V ^ {- 1} {\ boldsymbol {r}} \ right \ | _ {p} \ leq \ | V \ | _ {p} \ left \ | \ left (D- \ lambda ^ {a} I \ right) ^ {- 1} \ right \ | _ {p} \ left \ | V ^ {- 1} \ right \ | _ {p} \ | {\ boldsymbol {r}} \ | _ {p} = \ kappa _ {p} (V) \ left \ | \ left (D- \ lambda ^ {a} I \ right) ^ {- 1} \ right \ | _ {p} \ | {\ boldsymbol {r}} \ | _ {p}.}\ left \ | {\ boldsymbol {v} } ^ {a} \ right \ | _ {p} = \ left \ | V \ left (D- \ lambda ^ {a} I \ right) ^ {{- 1}} V ^ {{- 1}} { \ boldsymbol {r}} \ right \ | _ {p} \ leq \ | V \ | _ {p} \ left \ | \ left (D- \ lambda ^ {a} I \ right) ^ {{- 1} } \ right \ | _ {p} \ left \ | V ^ {{- 1}} \ right \ | _ {p} \ | {\ boldsymbol {r}} \ | _ {p} = \ kappa _ {p } (V) \ left \ | \ left (D- \ lambda ^ {a} I \ right) ^ {{- 1}} \ right \ | _ {p} \ | {\ boldsymbol {r}} \ | _ {p}.

Однако

(D - λ a I) - 1 {\ displaystyle \ left (D- \ lambda ^ {a} I \ right) ^ {- 1}}\ left (D- \ lambda ^ {a} I \ right) ^ {{- 1}}

- диагональная матрица, и ее p-норма легко вычисляется:

‖ (D - λ a I) - 1 ‖ p = max ‖ x ‖ p ≠ 0 ‖ (D - λ a I) - 1 x ‖ p ‖ X ‖ p = max λ ∈ σ (A) 1 | λ - λ a | = 1 min λ ∈ σ (A) | λ - λ a | {\ displaystyle \ left \ | \ left (D- \ lambda ^ {a} I \ right) ^ {- 1} \ right \ | _ {p} = \ max _ {\ | {\ boldsymbol {x}} \ | _ {p} \ neq 0} {\ frac {\ left \ | \ left (D- \ lambda ^ {a} I \ right) ^ {- 1} {\ boldsymbol {x}} \ right \ | _ { p}} {\ | {\ boldsymbol {x}} \ | _ {p}}} = \ max _ {\ lambda \ in \ sigma (A)} {\ frac {1} {\ left | \ lambda - \ lambda ^ {a} \ right |}} = {\ frac {1} {\ min _ {\ lambda \ in \ sigma (A)} \ left | \ lambda - \ lambda ^ {a} \ right |}}}\ left \ | \ left (D- \ lambda ^ {a} I \ right) ^ {{ -1}} \ right \ | _ {p} = \ max _ {{\ | {\ boldsymbol {x}} \ | _ {p} \ neq 0}} {\ frac {\ left \ | \ left (D - \ lambda ^ {a} I \ right) ^ {{- 1}} {\ boldsymbol {x}} \ right \ | _ {p}} {\ | {\ boldsymbol {x}} \ | _ {p} }} = \ max _ {{\ lambda \ in \ sigma (A)}} {\ frac {1} {\ left | \ lambda - \ lambda ^ {a} \ right |}} = {\ frac {1} {\ min _ {{\ lambda \ in \ sigma (A)}} \ left | \ lambda - \ lambda ^ {a} \ right |}}

откуда:

min λ ∈ λ (A) | λ - λ a | ≤ κ p (V) ‖ r ‖ p ‖ v a ‖ p. {\ displaystyle \ min _ {\ lambda \ in \ lambda (A)} \ left | \ lambda - \ lambda ^ {a} \ right | \ leq \ kappa _ {p} (V) {\ frac {\ | { \ boldsymbol {r}} \ | _ {p}} {\ left \ | {\ boldsymbol {v}} ^ {a} \ right \ | _ {p}}}.}{\ displaystyle \ min _ {\ lambda \ in \ lambda (A)} \ left | \ lambda - \ lambda ^ {a} \ right | \ leq \ kappa _ {p} (V) {\ frac {\ | {\ boldsymbol {r}} \ | _ {p}} {\ left \ | { \ boldsymbol {v}} ^ {a} \ right \ | _ {p}}}.}
Относительная граница

Обе формулировки теоремы Бауэра – Фике дают абсолютную оценку. Следующее следствие полезно всякий раз, когда требуется относительная оценка:

Следствие. Предположим, что A обратима и что μ является собственным значением A + δA. Тогда существует λ ∈ Λ (A) такое, что:
| λ - μ | | λ | ≤ κ п (V) ‖ A - 1 δ A ‖ p {\ displaystyle {\ frac {| \ lambda - \ mu |} {| \ lambda |}} \ leq \ kappa _ {p} (V) \ left \ | A ^ {- 1} \ delta A \ right \ | _ {p}}{\ frac {| \ lambda - \ mu |} {| \ lambda |}} \ leq \ kappa _ {p } (V) \ left \ | A ^ {{- 1}} \ delta A \ right \ | _ {p}

Примечание. || AδA || формально можно рассматривать как относительную вариацию A, так же как | λ - μ | / | λ | - относительное изменение λ.

Доказательство. Так как μ является собственным значением оператора A + δA и det (A) ≠ 0, умножая его на −A слева, получаем:

- A - 1 (A + δ A) v = - μ A - 1 v. {\ displaystyle -A ^ {- 1} (A + \ delta A) {\ boldsymbol {v}} = - \ mu A ^ {- 1} {\ boldsymbol {v}}.}-A ^ {{ -1}} (A + \ delta A) {\ boldsymbol {v}} = - \ mu A ^ {{- 1}} {\ boldsymbol {v}}.

Если мы установим:

A a = μ A - 1, (δ A) a = - A - 1 δ A {\ displaystyle A ^ {a} = \ mu A ^ {- 1}, \ qquad (\ delta A) ^ {a } = - A ^ {- 1} \ delta A}A ^ {a} = \ mu A ^ {{-1}}, \ qquad (\ delta A) ^ {a} = - A ^ {{- 1}} \ delta A

тогда имеем:

(A a + (δ A) a - I) v = 0 {\ displaystyle \ left (A ^ {a} + (\ delta A) ^ {a} -I \ right) {\ boldsymbol {v}} = {\ boldsymbol {0}}}\ left (A ^ {a} + (\ delta A) ^ {a} -I \ right) {\ boldsymbol {v}} = {\ boldsymbol {0}}

, что означает, что 1 является собственным значением A + (δA), с v как собственный вектор. Теперь собственные значения A равны μ / λ i, тогда как у него такое же, как у A. Применение теоремы Бауэра – Фике к A + (δA) с собственным значением 1 дает нам:

min λ ∈ Λ (A) | μ λ - 1 | = min λ ∈ Λ (A) | λ - μ | | λ | ≤ κ п (В) ‖ A - 1 δ A ‖ p {\ displaystyle \ min _ {\ lambda \ in \ Lambda (A)} \ left | {\ frac {\ mu} {\ lambda}} - 1 \ right | = \ min _ {\ lambda \ in \ Lambda (A)} {\ frac {| \ lambda - \ mu |} {| \ lambda |}} \ leq \ kappa _ {p} (V) \ left \ | A ^ {- 1} \ delta A \ right \ | _ {p}}{\ displaystyle \ min _ {\ lambda \ in \ Lambda (A)} \ left | {\ frac {\ mu} {\ lambda}} - 1 \ right | = \ min _ {\ lambda \ in \ Лямбда (A)} {\ frac {| \ lambda - \ mu |} {| \ lambda |}} \ leq \ kappa _ {p} (V) \ left \ | A ^ {- 1} \ delta A \ right \ | _ {p}}
Случай нормальных матриц

Если A нормальный, V является унитарным матрица, поэтому:

‖ V ‖ 2 = ‖ V - 1 ‖ 2 = 1, {\ displaystyle \ | V \ | _ {2} = \ left \ | V ^ {- 1} \ right \ | _ {2} = 1,}\ | V \ | _ {2} = \ left \ | V ^ {{- 1}} \ right \ | _ {2} = 1,

, так что κ 2 (V) = 1. Тогда теорема Бауэра – Фике принимает вид:

∃ λ ∈ Λ (A): | λ - μ | ≤ ‖ δ A ‖ 2 {\ displaystyle \ exists \ lambda \ in \ Lambda (A): \ quad | \ lambda - \ mu | \ leq \ | \ delta A \ | _ {2}}{\ displaystyle \ exists \ lambda \ in \ Lambda (A): \ quad | \ lambda - \ mu | \ leq \ | \ delta A \ | _ {2} }

Или альтернативно формулировка:

∃ λ ∈ Λ (A): | λ - λ a | ≤ ‖ р ‖ 2 ‖ ва ‖ 2 {\ displaystyle \ exists \ lambda \ in \ Lambda (A): \ quad \ left | \ lambda - \ lambda ^ {a} \ right | \ leq {\ frac {\ | { \ boldsymbol {r}} \ | _ {2}} {\ left \ | {\ boldsymbol {v}} ^ {a} \ right \ | _ {2}}}}{\ di splaystyle \ exists \ lambda \ in \ Lambda (A): \ quad \ left | \ lambda - \ lambda ^ {a} \ right | \ leq {\ frac {\ | {\ boldsymbol {r}} \ | _ {2 }} {\ left \ | {\ boldsymbol {v}} ^ {a} \ right \ | _ {2}}}}

что, очевидно, остается верным, если A a Эрмитова матрица. В этом случае, однако, имеет место гораздо более сильный результат, известный как теорема Вейля о собственных значениях. В эрмитовом случае можно также переформулировать теорему Бауэра – Фике в виде, что отображение A ↦ Λ (A), которое отображает матрицу в ее спектр, является нерасширяющей функцией относительно расстояния Хаусдорфа на множестве компактных подмножеств C.

Ссылки
  • Bauer, FL; Фике, К. Т. (1960). «Нормы и теоремы исключения». Нумер. Математика. 2 (1): 137–141. doi : 10.1007 / BF01386217.
  • Eisenstat, S.C.; Ипсен, И. К. Ф. (1998). «Три абсолютных границы возмущения для собственных значений матрицы подразумевают относительные границы». Журнал SIAM по матричному анализу и приложениям. 20 (1): 149–158. CiteSeerX 10.1.1.45.3999. doi :10.1137/S0895479897323282.
Последняя правка сделана 2021-05-12 07:38:10
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте