Обратная матрица,ее свойства.обращение матриц методом жордана-гаусса,условие обратимости матриц

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

Понятие обратной матрицы существует только для квадратных матриц, матриц «два на два», «три на три» и т.д.

Обозначения: Как вы уже, наверное, заметили, обратная матрица обозначается надстрочным индексом

Начнем с простейшего случая – матрицы «два на два». Чаще всего, конечно, требуется найти обратную матрицу для матрицы «три на три», но, тем не менее, настоятельно рекомендую изучить более простое задание, для того чтобы усвоить общий принцип решения.

Пример:

Найти обратную матрицу для матрицы

Решаем. Последовательность действий удобно разложить по пунктам.

1) Сначала находим определитель матрицы.

Если с пониманием сего действа плоховато, ознакомьтесь с материалом Как вычислить определитель?

Важно! В том случае, если определитель матрицы равен НУЛЮ – обратной матрицы НЕ СУЩЕСТВУЕТ.

В рассматриваемом примере, как выяснилось, , а значит, всё в порядке.

2) Находим матрицу миноров .

Для решения нашей задачи не обязательно знать, что такое минор, однако, желательно ознакомиться со статьей Как вычислить определитель.

Матрица миноров имеет такие же размеры, как и матрица , то есть в данном случае .
Дело за малым, осталось найти четыре числа и поставить их вместо звездочек.

Возвращаемся к нашей матрице
Сначала рассмотрим левый верхний элемент:

Как найти его минор?
А делается это так: МЫСЛЕННО вычеркиваем строку и столбец, в котором находится данный элемент:

Оставшееся число и является минором данного элемента, которое записываем в нашу матрицу миноров:

Рассматриваем следующий элемент матрицы :

Мысленно вычеркиваем строку и столбец, в котором стоит данный элемент:

То, что осталось, и есть минор данного элемента, который записываем в нашу матрицу:

Аналогично рассматриваем элементы второй строки и находим их миноры:


Готово.

– матрица миноров соответствующих элементов матрицы .

3) Находим матрицу алгебраических дополнений .

Это просто. В матрице миноров нужно ПОМЕНЯТЬ ЗНАКИ у двух чисел:

Именно у этих чисел, которые я обвел в кружок!

– матрица алгебраических дополнений соответствующих элементов матрицы .

И всего-то лишь…

4) Находим транспонированную матрицу алгебраических дополнений.

Что такое транспонирование матрицы, и с чем это едят, смотрите в лекции Действия с матрицами.

– транспонированная матрица алгебраических дополнений соответствующих элементов матрицы .

5) Ответ.

Вспоминаем нашу формулу
Всё найдено!

Таким образом, обратная матрица:

Ответ лучше оставить в таком виде. НЕ НУЖНО делить каждый элемент матрицы на 2, так как получатся дробные числа. Более подробно данный нюанс рассмотрен в той же статье Действия с матрицами.

Как проверить решение?

Необходимо выполнить матричное умножение либо

Проверка:

Получена так называемая единичная матрица (с единицами по главной диагонали и нулями в остальных местах).

Таким образом, обратная матрица найдена правильно.

Если провести действие , то в результате тоже получится единичная матрица. Это один из немногих случаев, когда умножение матриц перестановочно, более подробную информацию можно найти в статье Свойства операций над матрицами. Матричные выражения. Также заметьте, что в ходе проверки константа (дробь) выносится вперёд и обрабатывается в самом конце – после матричного умножения. Это стандартный приём.

Переходим к более распространенному на практике случаю – матрице «три на три».

Пример:

Найти обратную матрицу для матрицы

Алгоритм точно такой же, как и для случая «два на два».

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

1) Находим определитель матрицы.

Здесь определитель раскрыт по первой строке.

Также не забываем, что , а значит, всё нормально – обратная матрица существует.

2) Находим матрицу миноров .

Матрица миноров имеет размерность «три на три» , и нам нужно найти девять чисел.

Я подробно рассмотрю парочку миноров:

Рассмотрим следующий элемент матрицы:

МЫСЛЕННО вычеркиваем строку и столбец, в котором находится данный элемент:

Оставшиеся четыре числа записываем в определитель «два на два»

Этот определитель «два на два» и является минором данного элемента. Его нужно вычислить:

Всё, минор найден, записываем его в нашу матрицу миноров:

Как вы, наверное, догадались, необходимо вычислить девять определителей «два на два». Процесс, конечно, муторный, но случай не самый тяжелый, бывает хуже.

Ну и для закрепления – нахождение еще одного минора в картинках:

Остальные миноры попробуйте вычислить самостоятельно.

Окончательный результат:
– матрица миноров соответствующих элементов матрицы .

То, что все миноры получились отрицательными – чистая случайность.

3) Находим матрицу алгебраических дополнений .

В матрице миноров необходимо СМЕНИТЬ ЗНАКИ строго у следующих элементов:

В данном случае:
– матрица алгебраических дополнений соответствующих элементов матрицы .

4) Находим транспонированную матрицу алгебраических дополнений.

– транспонированная матрица алгебраических дополнений соответствующих элементов матрицы .

5) Ответ:

Проверка:

Свойства обратной матрицы[править | править исходный текст]

· , где обозначает определитель.

· для любых двух обратимых матриц и .

· где обозначает транспонированную матрицу.

· для любого коэффициента .

· Если необходимо решить систему линейных уравнений , (b — ненулевой вектор) где — искомый вектор, и если существует, то . В противном случае либо размерность пространстварешений больше нуля, либо их нет вовсе.

Метод Гаусса—Жордана[править | править исходный текст]

Возьмём две матрицы: саму A и единичную E. Приведём матрицу A к единичной матрице методом Гаусса—Жордана. После применения каждой операции к первой матрице применим ту же операцию ко второй. Когда приведение первой матрицы к единичному виду будет завершено, вторая матрица окажется равной A⁻¹.

При использовании метода Гаусса первая матрица будет умножаться слева на одну из элементарных матриц (трансвекцию или диагональную матрицу с единицами на главной диагонали, кроме одной позиции):

.

.

Вторая матрица после применения всех операций станет равна , то есть будет искомой. Сложность алгоритма — .

С помощью матрицы алгебраических дополнений[править | править исходный текст]

где — присоединенная матрица;

Полученная матрица A ⁻¹ и будет обратной. Сложность алгоритма зависит от сложности алгоритма расчета определителя O_det и равна O(n²)·O_det.

Использование LU/LUP-разложения[править | править исходный текст]

Матричное уравнение для обратной матрицы можно рассматривать как совокупность систем вида . Обозначим -ый столбец матрицы через ; тогда , ,поскольку -м столбцом матрицы является единичный вектор . другими словами, нахождение обратной матрицы сводится к решению n уравнений с одной матрицей и разными правыми частями. После выполнения LUP-разложения (время O(n³)) на решение каждого из n уравнений нужно время O(n²), так что и эта часть работы требует времени O(n³)^[1].

Если матрица A невырождена, то для неё можно рассчитать LUP-разложение . Пусть , . Тогда из свойств обратной матрицы можно записать: . Если умножить это равенство на U и L то можно получить два равенства вида и . Первое из этих равенств представляет собой систему из n² линейных уравнений для из которых известны правые части (из свойств треугольных матриц). Второе представляет также систему из n² линейных уравнений для из которых известны правые части (также из свойств треугольных матриц). Вместе они представляют собой систему из n² равенств. С помощью этих равенств можно реккурентно определить все n² элементов матрицы D. Тогда из равенства (PA)⁻¹ = A⁻¹P⁻¹ = B⁻¹ = D. получаем равенство .

В случае использования LU-разложения не требуется перестановки столбцов матрицы D но решение может разойтись даже если матрица A невырождена.

Сложность алгоритма — O(n³).