Существует ли обобщение закона инерции Сильвестра для симметричной обобщенной задачи на собственные значения?


9

Я знаю, что для решения симметричной задачи на собственные значения Ax=λxмы можем использовать закон инерции Сильвестра, то есть число собственных значений A меньше, чем aравно количеству отрицательных элементов в где диагональная матрица получается из факторизации ЛПНП . Затем с помощью метода деления пополам мы можем найти все или некоторые собственные значения по желанию. Я хотел бы знать, существует ли обобщение закона инерции Сильвестра для симметричных обобщенных задач на собственные значения, то есть решение , где и - симметричные матрицы. Спасибо.DDAaI=LDLTAx=λBxAB

Ответы:


5

Да, если карандаш определен, т. Е. Если и эрмитовы, а положительно определен. Затем подпись имеет ту же интерпретацию проблемы собственных значений , как и в случае . Более общий результат такого рода справедлив для любой определенной нелинейной задачи на собственные значения . Смотрите раздел 5.3 моей книгиABBAσB(AλB)x=0B=IA(λ)x=0

Арнольд Ноймайер, Введение в численный анализ, Кембридж Univ. Press, Cambridge 2001.

Для доказательство моего утверждения может быть выведено из аргумента, приведенного Джеком Полсоном, отметив, что и являются конгруэнтными, следовательно, имеют одинаковую инерцию.(AλB)x=0CσIAσB

В частности, можно непосредственно вычислить инерцию , и не нуждается в Холецкоге разложения с образованием . Действительно, если плохо обусловлен, то численное формирование ухудшает качество теста на инерцию.AσBBCBC


Хороший вопрос о плохой обусловленности B; Я думаю, что ваш подход лучше, если вы действительно заинтересованы только в вычислении инерции. Предложенный мною подход типичен для фактического решения проблемы собственных значений (в случае, когдаBв хорошем состоянии).
Джек Поулсон

@JackPoulson: тест инерции обычно применяется для получения собственных значений в определенном интервале, когда A а также B являются редкими, и их совместный образец разреженности генерирует не слишком много заполнения. Но ваш C будет плотным, когда Bявляется трехдиагональной, поэтому его использование никогда не подходит для нахождения собственных значений большой разреженной обобщенной задачи на собственные значения. (Принимая во внимание, что если проблема невелика, использование инерции не имеет смысла, так как поиск всех собственных значений обычно достаточно быстр.)
Арнольд Ноймайер,

Безусловно; кажется, что я ошибочно пропустил слово «плотный» из моего комментария.
Джек Полсон

3

В случае, когда B эрмитова и положительно-определенная, холесская факторизация B, сказать B=LLHдаёт

Ax=LLHxλ,

и этим уравнением можно манипулировать, чтобы показать, что

(L1ALH)(LHx)=(LHx)λ,

где должно быть ясно, что CL1ALH сохраняет симметрию A, а также имеет тот же спектр, что и карандаш (A,B), Таким образом, после формированияCс факторизацией Холецкого, за которой следует двустороннее треугольное решение , вы можете напрямую применить закон инерции Сильвестра кC собрать информацию о собственных значениях карандаша (A,B),

Отметим, что поскольку закон инерции Сильвестра инвариантен относительно конгруэнтных преобразований, например,SSHтогда матрица C соответствует A через преобразование L1LH, и так C имеет ту же инерцию, что и A, Однако, если инерцияCσI желательно, для некоторого ненулевого сдвига σто мы уже не можем просто считать A,


Понижение без какой-либо конструктивной критики?
Джек Поулсон,

2
Я не выходил из системы на компьютере в моем офисе, и мой товарищ по офису случайно наткнулся на эту вкладку в моем браузере и отклонил ответ. Я извиняюсь за недоразумение и спрошу его, почему он отказался от этого.
Шухао Цао

Вы были абсолютно правы, когда B спд матрица для пары (A,B)мы могли бы просто посмотреть на Aчтобы получить то, что мы хотим. Однако мой коллега сказал, что вы не ответили на вопрос, еслиBимеет симметричность Извините за путаницу.
Шухао Цао

@ Джон: вздох. Это не то, для чего понижательный голос.
Джек Полсон

Я знаю! Я уже сказал ему «пожалуйста, прочитайте правило» после того, как обнаружил, что он использовал мою учетную запись, чтобы понизить соответствующий ответ!
Шухао Цао
Используя наш сайт, вы подтверждаете, что прочитали и поняли нашу Политику в отношении файлов cookie и Политику конфиденциальности.
Licensed under cc by-sa 3.0 with attribution required.