Диагонализация больших матриц Инвариантные подпространства. Матрица Рэлея.

Сентябрь 2, 2022

Главная
Алгебра
Диагонализация больших матриц Инвариантные подпространства. Матрица Рэлея.

Содержание

2. Постановки задач Размерность узельного базиса для бозе-системы из 12 узлов и 12 частиц: при решении квантовых
3. Пространства и инвариантные подпространства Скалярное произведение двух вещественных векторов: Линейная оболочка, натянутая на векторы s1, s2,
4. Матрица Рэлея Пусть Q=(q1,q2,…,qm) – векторы ортонормированного базиса инвариантного подпространства Σ, упорядоченные в виде матрицы Q
5. Процедура Рэлея – Ритца Процедура Рэлея – Ритца обеспечивает наилучшее приближение к собственным парам матрицы A
6. Подпространство Крылова Подпространство Крылова для произвольного ненулевого вектора b: При достаточно большом m подпространство Крылова будет
7. Алгоритм Ланцоша Стартовый вектор b выбирается либо в произвольном виде, либо с использованием априорной информации о
8. Алгоритм Ланцоша Пары Ритца, аппроксимирующие собственные пары матрицы A: Контроль сходимости – норма вектора невязки: Существуют
9. Пример Для случая m=4:
11. Скачать презентацию

Слайд 2

Постановки задач
Размерность узельного базиса для бозе-системы из 12 узлов и 12

частиц:
при решении квантовых многочастичных задач редко приходится вычислять все собственные пары (собственный вектор и отвечающее ему собственное значение) гамильтоновой матрицы
Обычные постановки задач:
Вычисление некоторых собственных значений, принадлежащих некоторому интервалу, или крайних собственных значений (например, низ спектра)
Нахождение некоторых собственных пар
Определение всех собственных значений или большого их количества без вычисления собственных векторов
Одним из самых мощных методов диагонализации симметричных матриц является алгоритм Ланцоша

Слайд 3

Пространства и инвариантные подпространства
Скалярное произведение двух вещественных векторов:
Линейная оболочка, натянутая на векторы

s1, s2, … :
Множество S называется образующим множеством подпространства Σ
Образующее множество, содержащее наименьшее количество векторов, называется базисом подпространства, а количество векторов базиса называется размерностью подпространства
Если базисные векторы ортонормированны, то базис называется ортонормированным
Подпространство Σ инвариантно относительно A, если для любого вектора x из Σ следует, что вектор Ax также принадлежит Σ

Слайд 4

Матрица Рэлея
Пусть Q=(q1,q2,…,qm) – векторы ортонормированного базиса инвариантного подпространства Σ, упорядоченные

в виде матрицы Q размера n×m
Матрица Рэлея C – сужение A на Σ:
Матрица Рэлея является симметричной
Определение собственных пар матрицы A может быть осуществлено путем диагонализации матрицы C меньшего размера:

Слайд 5

Процедура Рэлея – Ритца
Процедура Рэлея – Ритца обеспечивает наилучшее приближение к

собственным парам матрицы A
1. Ортонормализация матрицы Q, вычисление матрицы Рэлея
Матрица Рэлея H почти инвариантна относительно A
2. Определение требуемого количества собственны пар матрицы H:
3. Полученные значения Ритца ?? будут наилучшими приближениями к собственным значениям матрицы A, а векторы Ритца
– соответствующими приближениями к собственным векторам матрицы A

Слайд 6

Подпространство Крылова
Подпространство Крылова для произвольного ненулевого вектора b:
При достаточно большом m

подпространство Крылова будет посчти инвариантно относительно A
В этом случае матрица Рэлея является трехдиагональной, и вычисление пар Ритца существенно упрощается
Столбцы матрицы Q могут быть связаны друг с другом с помощью трехчленных рекуррентных соотношений, позволяющих вычислять новые векторы базиса Q
Алгоритм Ланцоша обычно использует последовательность подпространств Крылова
и вычисляет пары Ритца матрицы A, отвечающие каждому из этих подпространств. Сходимость значений собственных пар обычно оказывается достаточно быстрой:

Слайд 7

Алгоритм Ланцоша
Стартовый вектор b выбирается либо в произвольном виде, либо с

использованием априорной информации о собственных векторах матрицы A:
1. Пополнение ортонормированного базиса:
2. Вычисление промежуточной невязки:
3. Вычисление очередного диагонального элемента:
4. Завершение вычисления невязки:
5. Вычисление нормы невязки:

Слайд 8

Алгоритм Ланцоша
Пары Ритца, аппроксимирующие собственные пары матрицы A:
Контроль сходимости – норма

вектора невязки:
Существуют версии алгоритма Ланцоша, предназначенные для отыскания крайних собственных значений, некоторых внутренних собственных значений и даже всех собственных значений симметричных или эрмитовых матриц очень больших размеров, некоторые версии алгоритма позволяют вычислять и соответствующие собственные векторы
При практической реализации алгоритма Ланцоша возникают погрешности, связанные с машинной точностью обработки чисел двойной точности. Для решения этой проблемы разработаны специальные модификации алгоритма, производящие переортогонализацию векторов в зависимости от величины погрешности

Слайд 9

Диагонализация больших матриц Инвариантные подпространства. Матрица Рэлея.

Содержание

Постановки задачРазмерность узельного базиса для бозе-системы из 12 узлов и 12

Пространства и инвариантные подпространстваСкалярное произведение двух вещественных векторов:Линейная оболочка, натянутая на векторы

Матрица РэлеяПусть Q=(q1,q2,…,qm) – векторы ортонормированного базиса инвариантного подпространства Σ, упорядоченные

Процедура Рэлея – РитцаПроцедура Рэлея – Ритца обеспечивает наилучшее приближение к

Подпространство КрыловаПодпространство Крылова для произвольного ненулевого вектора b:При достаточно большом m

Алгоритм ЛанцошаСтартовый вектор b выбирается либо в произвольном виде, либо с

Алгоритм ЛанцошаПары Ритца, аппроксимирующие собственные пары матрицы A:Контроль сходимости – норма

ПримерДля случая m=4:

Похожие презентации