Обратная матрица
Обра́тная ма́трица — такая матрица , при умножении которой на исходную матрицу получается единичная матрица :
Обратную матрицу можно определить как:
- где — соответствующая присоединённая матрица,
- — определитель матрицы .
Из этого определения следует критерий обратимости: матрица обратима тогда и только тогда, когда она невырождена, то есть её определитель не равен нулю. Для неквадратных матриц и вырожденных матриц обратных матриц не существует. Однако возможно обобщить это понятие и ввести псевдообратные матрицы, похожие на обратные по многим свойствам.
Свойства обратной матрицы
правитьПусть квадратные матрицы — невырожденные. Тогда:
- Если обратима, то единственна.
- , где обозначает определитель.
- , где обозначает операцию транспонирования.
- для любого коэффициента .
- .
- Пусть необходимо решить систему линейных уравнений , где — искомый вектор, — ненулевой вектор. Если существует, то . В противном случае либо размерность пространства решений больше нуля, либо их нет вовсе.
Способы нахождения обратной матрицы
правитьЕсли матрица обратима, то для нахождения обратной матрицы можно воспользоваться одним из следующих способов:
Точные (прямые) методы
правитьМетод Жордана—Гаусса
правитьВозьмём две матрицы: саму и единичную матрицу . Приведём матрицу к единичной методом Гаусса—Жордана, применяя преобразования по строкам (можно также применять преобразования и по столбцам). После применения каждой операции к первой матрице применим ту же операцию ко второй. Когда приведение первой матрицы к единичному виду будет завершено, вторая матрица окажется равной .
При использовании метода Гаусса первая матрица будет умножаться слева на одну из элементарных матриц (трансвекцию или диагональную матрицу с единицами на главной диагонали, кроме одной позиции):
Вторая матрица после применения всех операций станет равна , то есть будет искомой. Сложность алгоритма — .
С помощью матрицы алгебраических дополнений
правитьМатрица, обратная матрице , представима в виде:
- где — присоединенная матрица (матрица, составленная из алгебраических дополнений для соответствующих элементов транспонированной матрицы).
Сложность алгоритма зависит от сложности алгоритма расчета определителя и равна .
Использование LU- или LUP-разложения
правитьМатричное уравнение для обратной матрицы можно рассматривать как совокупность систем вида . Обозначим -й столбец матрицы через ; тогда , , поскольку -м столбцом матрицы является единичный вектор . Иными словами, нахождение обратной матрицы сводится к решению уравнений с одной матрицей и разными правыми частями. Решение этих уравнений может быть упрощено с помощью LU- или LUP-разложения матрицы . После выполнения LUP-разложения за время на решение каждого из уравнений нужно время , так что и этот алгоритм требует времени [1].
Матрицу, обратную к заданной невырожденной матрице , можно также вычислить непосредственно с помощью матриц, полученных в результате разложения.
Результатом LUP-разложения матрицы является равенство . Пусть , . Тогда из свойств обратной матрицы можно записать: . Если умножить это равенство на и то можно получить два равенства вида и . Первое из этих равенств представляет собой систему из линейных уравнений, для из которых известны правые части (из свойств треугольных матриц). Второе также представляет систему из линейных уравнений, для из которых известны правые части (также из свойств треугольных матриц). Вместе они представляют собой систему из равенств. С их помощью можно рекуррентно определить все элементов матрицы . Тогда из равенства получаем равенство .
В случае использования LU-разложения ( ) не требуется перестановки столбцов матрицы , но решение может разойтись даже если матрица невырождена.
Сложность обоих алгоритмов — .
Итерационные методы
правитьМатрицу можно вычислить с произвольной точностью в результате выполнения следующего итерационного процесса, называющегося методом Шульца и являющегося обобщением классического метода Ньютона:
Последовательность матриц сходится к при . Существует также так называемый обобщённый метод Шульца, который описывается следующими рекуррентными соотношениями[2]:
Выбор начального приближения
правитьПроблема выбора начального приближения в рассматриваемых здесь процессах итерационного обращения матриц не позволяет относиться к ним как к самостоятельным универсальным методам, конкурирующими с прямыми методами обращения, основанными, например, на -разложении матриц. Имеются некоторые рекомендации по выбору , обеспечивающие выполнение условия (спектральный радиус матрицы меньше единицы), являющегося необходимым и достаточным для сходимости итерационного процесса. Однако при этом, во-первых, требуется знать оценку сверху спектра обращаемой матрицы либо матрицы (а именно, если — симметричная положительно определённая матрица и , то можно взять , где ; если же — произвольная невырожденная матрица и , то полагают , где также ; можно, конечно, упростить ситуацию и, воспользовавшись тем, что , положить ). Во-вторых, при таком задании начальной матрицы нет гарантии, что будет малой (возможно, даже окажется ), и высокий порядок скорости сходимости обнаружится далеко не сразу.
Для метода Ньютона в качестве начального приближения можно выбрать , где верхний индекс обозначает эрмитово сопряжение, и — соответствующие матричные нормы. Такое вычисляется всего за операций, где — порядок матрицы, и обеспечивает сходимость алгоритма[3].
Примеры
правитьМатрица 2 × 2
править- [4].
Обращение матрицы 2 × 2 возможно только при условии, что .
Примечания
править- ↑ Кормен Т., Лейзерсон Ч., Ривест Р., Штайн К. Алгоритмы: построение и анализ, — М.: Вильямс, 2006 (с. 700).
- ↑ Petković, M. D.. Generalized Schultz iterative methods for the computation of outer inverses (англ.) // Computers & Mathematics with Applications. — 2014. — June (vol. 67, iss. 10). — P. 1837—1847. — doi:10.1016/j.camwa.2014.03.019.
- ↑ Pan, V., Reif, J.. Fast and efficient parallel solution of dense linear systems (англ.) // Computers & Mathematics with Applications. — 1989. — Vol. 17, iss. 11. — P. 1481—1491. — doi:10.1016/0898-1221(89)90081-3.
- ↑ Как найти обратную матрицу? mathprofi.ru. Дата обращения: 18 октября 2017. Архивировано 17 октября 2017 года.
Ссылки
правитьВ статье не хватает ссылок на источники (см. рекомендации по поиску). |