Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
BMLA.DOC
Скачиваний:
19
Добавлен:
09.11.2019
Размер:
3.4 Mб
Скачать

Якобиевы матрицы

Вещественная матрица

,

называется якобиевой (у нас ).

Лемма 3.

Пусть – якобиева матрица, тогда

  1. если , то , если , то .

Док–во

оставляется читателю в качестве упражнения.

Лемма 4.

Собственные значения якобиевой матрицы попарно различные (простые).

Док–во.

Т.к. размерность ядра симметричной матрицы совпадает с кратностью , а из леммы 3 следует, что у вырожденной якобиевой матрицы минор , то и простое собственное значение матрицы .

Теорема.

Пусть – якобиева матрица, тогда

,

если приписать знак .

Док–во.

1. Если , то это лемма 1.

2. Пусть . Пусть .

Определим и рассмотрим якобиевы матрицы .

Т.к. , то

а) (т.к. определитель матрицы равен произведению ее собственных значений),

б) ,

в) ,

(т.к. из леммы 4 следует, что простое и отрицательных собственных значений у матрицы на одно больше, чем у матрицы ),

г) .

Из леммы 1, а) и г) следует, что

,

,

Подсчитаем эти числа:

Из б) следует, что если и , то перемена знака происходит (или нет) одновременно в этих последовательностях.

Случай .

Из леммы 3 имеем , отсюда и из б) следует и на участках

по одной перемене знака.

Случай . Отсюда, из в) и г) следует, что

, ,

.

Следовательно, (если приписать знак ) последовательности миноров матриц и имеют одинаковые знаки. Теорема доказана.

О вычислении чпз

Для вычисления якобиевой матрицы достаточно знать знак каждого . Если

(обычно выбирают ), то и . Нормировку можно применять неоднократно, что позволит избежать быстрого роста (переполнения) чисел .

О вычислении собственного вектора

Лемма 5.

Последняя компонента собственного вектора якобиевой матрицы не равна нулю.

Док–во.

Пусть . Предположим, что . Тогда

– противоречие, значит .

Собственный вектор якобиевой матрицы мы можем, положив , вычислить по формулам

или решив систему

с неособенной матрицей.

Лекция 13. Метод вращений (Якоби)

Для самосопряженной матрицы существует унитарная матрица (столбцы которой – собственные векторы матрицы ):

, где .

Идея:

построить : , тогда на диагональные элементы будут приближать собственные значения, а столбцы – собственные векторы матрицы .

Определим .

Лемма 1.

Для любых квадратной матрицы и унитарной матрицы имеем

.

Док–во.

Если , то

В качестве матриц будем выбирать элементарные матрицы вращения.

Лемма 2.

Пусть , ,

где – элементарная матрица вращения, тогда

.

Док–во.

Заметим, что изменились только строки и столбцы с номерами , .

Тогда, используя лемму 1, получим

откуда следует утверждение леммы.

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]