Обратная матрица 4 порядка пример. Нахождение обратной матрицы

Нахождение обратной матрицы - задача, которая чаще решается двумя методами:

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

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

Обратной матрицей А , называется такая матрица

А
. (1)

Обратной матрицей , которую требуется отыскать для данной квадратной матрицы А , называется такая матрица

произведение на которую матрицы А справа является единичной матрицей, т.е,
. (1)

Единичной матрицей называется диагональная матрица, у которой все диагональные элементы равны единице.

Теорема. Для каждой неособенной (невырожденной, несингулярной) квадратной матрицы можно найти обратную матрицу, и притом только одну. Для особенной (вырожденной, сингулярной) квадратной матрицы обратная матрица не существует.

Квадратная матрица называется неособенной (или невырожденной , несингулярной ), если её определитель не равен нулю, и особенной (или вырожденной , сингулярной ), если её определитель равен нулю.

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

Для обратной матрицы существует уместная аналогия с обратным числом. Для каждого числа a , не равного нулю, существует такое число b , что произведение a и b равно единице: ab = 1 . Число b называется обратным для числа b . Например, для числа 7 обратным является число 1/7, так как 7*1/7=1.

Нахождение обратной матрицы методом алгебраических дополнений (союзной матрицы)

Для неособенной квадратной матрицы А обратной является матрица

где - определитель матрицы А , а - матрица, союзная с матрицей А .

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

то

Алгоритм нахождения обратной матрицы методом алгебраических дополнений

1. Найти определитель данной матрицы A . Если определитель равен нулю, нахождение обратной матрицы прекращается, так как матрица вырожденная и обратная для неё не существует.

2. Найти матрицу, транспонированную относительно A .

3. Вычислить элементы союзной матрицы как алгебраические дополнения марицы, найденной на шаге 2.

4. Применить формулу (2): умножить число, обратное определителю матрицы A , на союзную матрицу, найденную на шаге 4.

5. Проверить полученный на шаге 4 результат, умножив данную матрицу A на обратную матрицу. Если произведение этих матриц равно единичной матрицы, значит обратная матрица была найдена верно. В противном случае начать процесс решения снова.

Пример 1. Для матрицы

найти обратную матрицу.

Решение. Для нахождения обратной матрицы необходимо найти определитель матрицы А . Находим по правилу треугольников:

Следовательно, матрица А – неособенная (невырожденная, несингулярная) и для неё существует обратная.

Найдём матрицу, союзную с данной матрицей А .

Найдём матрицу , транспонированную относительно матрицы A :

Вычисляем элементы союзной матрицы как алгебраические дополнения матрицы, транспонированной относительно матрицы A :

Следовательно, матрица , союзная с матрицей A , имеет вид

Замечание. Порядок вычисления элементов и транспонирования матрицы может быть иным. Можно сначала вычислить алгебраические дополнения матрицы A , а затем транспонировать матрицу алгебраических дополнений. В результате должны получиться те же элементы союзной матрицы.

Применяя формулу (2), находим матрицу, обратную матрице А :

Нахождение обратной матрицы методом исключения неизвестных Гаусса

Первый шаг для нахождения обратной матрицы методом исключения неизвестных Гаусса - приписать к матрице A единичную матрицу того же порядка, отделив их вертикальной чертой. Мы получим сдвоенную матрицу . Умножим обе части этой матрицы на , тогда получим

Алгоритм нахождения обратной матрицы методом исключения неизвестных Гаусса

1. К матрице A приписать единичную матрицу того же порядка.

2. Полученную сдвоенную матрицу преобразовать так, чтобы в левой её части получилась единичная матрица, тогда в правой части на месте единичной матрицы автоматически получится обратная матрица. Матрица A в левой части преобразуется в единичную матрицу путём элементарных преобразований матрицы.

2. Если в процессе преобразования матрицы A в единичную матрицу в какой-либо строке или в каком-либо столбце окажутся только нули, то определитель матрицы равен нулю, и, следовательно, матрица A будет вырожденной, и она не имеет обратной матрицы. В этом случае дальнейшее нахождение обратной матрицы прекращается.

Пример 2. Для матрицы

найти обратную матрицу.

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

Умножим первую строку левой и правой матрицы на (-3) и сложим её со второй строкой, а затем умножим первую строку на (-4) и сложим её с третьей строкой, тогда получим

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

Сложим первую строку со второй, а затем умножим вторую строку на (-9) и сложим её с третьей строкой. Тогда получим

Разделим третью строку на 8, тогда

Умножим третью строку на 2 и сложим её со второй строкой. Получается:

Переставим местами вторую и третью строку, тогда окончательно получим:

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

Можно проверить правильность вычислений, умножим исходную матрицу на найденную обратную матрицу:

В результате должна получиться обратная матрица.

Пример 3. Для матрицы

найти обратную матрицу.

Решение. Составляем сдвоенную матрицу

и будем её преобразовывать.

Первую строку умножаем на 3, а вторую на 2, и вычитаем из второй, а затем первую строку умножаем на 5, а третью на 2 и вычитаем из третьей строки, тогда получим

Первую строку умножаем на 2 и складываем её со второй, а затем из третьей строки вычитаем вторую, тогда получим

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

Способы нахождения обратной матрицы, . Рассмотрим квадратную матрицу

Обозначим Δ =det A.

Квадратная матрица А называется невырожденной, или неособенной , если ее определитель отличен от нуля, и вырожденной, или особенной , если Δ = 0.

Квадратная матрица В есть для квадратной матрицы А того же порядка, если их произведение А В = В А = Е, где Е - единичная матрица того же порядка, что и матрицы А и В.

Теорема . Для того, чтобы матрица А имела обратную матрицу, необходимо и достаточно, чтобы ее определитель был отличен от нуля.

Обратная матрица матрице А, обозначается через А - 1 , так что В = А - 1 и вычисляется по формуле

, (1)

где А i j - алгебраические дополнения элементов a i j матрицы A..

Вычисление A -1 по формуле (1) для матриц высокого порядка очень трудоемко, поэтому на практике бывает удобно находить A -1 с помощью метода элементарных преобразований (ЭП). Любую неособенную матрицу А путем ЭП только столбцов (или только строк) можно привести к единичной матрице Е. Если совершенные над матрицей А ЭП в том же порядке применить к единичной матрице Е, то в результате получится обратная матрица. Удобно совершать ЭП над матрицами А и Е одновременно, записывая обе матрицы рядом через черту. Отметим еще раз, что при отыскании канонического вида матрицы с целью нахождения можно пользоваться преобразованиями строк и столбцов. Если нужно найти обратную матрицу, в процессе преобразований следует использовать только строки или только столбцы.

Пример 2.10 . Для матрицы найти A -1 .

Решение. Находим сначала детерминант матрицы А
значит, обратная матрица существует и мы ее можем найти по формуле: , где А i j (i,j=1,2,3) - алгебраические дополнения элементов а i j исходной матрицы.

Откуда .

Пример 2.11 . Методом элементарных преобразований найти A -1 для матрицы: А= .

Решение. Приписываем к исходной матрице справа единичную матрицу того же порядка: . С помощью элементарных преобразований столбцов приведем левую “половину” к единичной, совершая одновременно точно такие преобразования над правой матрицей.
Для этого поменяем местами первый и второй столбцы: ~ . К третьему столбцу прибавим первый, а ко второму - первый, умноженный на -2: . Из первого столбца вычтем удвоенный второй, а из третьего - умноженный на 6 второй; . Прибавим третий столбец к первому и второму: . Умножим последний столбец на -1: . Полученная справа от вертикальной черты квадратная матрица является обратной матрицей к данной матрице А. Итак,
.

Определение 1: матрица называется вырожденной, если её определитель равен нулю.

Определение 2: матрица называется невырожденной, если её определитель не равен нулю.

Матрица "A" называется обратной матрицей , если выполняется условие A*A-1 = A-1 *A = E (единичной матрице).

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

Схема вычисления обратной матрицы:

1) Вычислить определитель матрицы "A", если ∆ A = 0, то обратной матрицы не существует.

2) Найти все алгебраические дополнения матрицы "A".

3) Составить матрицу из алгебраических дополнений (Aij )

4) Транспонировать матрицу из алгебраических дополнений (Aij )T

5) Умножить транспонированную матрицу на число, обратное определителю данной матрицы.

6) Выполнить проверку:

На первый взгляд может показаться, что это сложно, но на самом деле всё очень просто. Все решения основаны на простых арифметических действиях, главное при решении не путаться со знаками "-" и "+", и не терять их.

А теперь давайте вместе с Вами решим практическое задание, вычислив обратную матрицу.

Задание: найти обратную матрицу "A", представленную на картинке ниже:

Решаем всё в точности так, как это указано в план-схеме вычисления обратной матрицы.

1. Первое, что нужно сделать, это найти определитель матрицы "A":

Пояснение:

Мы упростили наш определитель, воспользовавшись его основными функциями. Во первых, мы прибавили ко 2 и 3 строке элементы первой строки, умноженные на одно число.

Во-вторых, мы поменяли 2 и 3 столбец определителя, и по его свойствам поменяли знак перед ним.

В-третьих, мы вынесли общий множитель (-1) второй строки, тем самым, снова поменяв знак, и он стал положительным. Также мы упростили 3 строку также, как в самом начале примера.

У нас получилась треугольный определитель, у которого элементы ниже диагонали равны нулю, и по 7 свойству он равен произведению элементов диагонали. В итоге мы получили ∆ A = 26, следовательно обратная матрица существует.

А11 = 1*(3+1) = 4

А12 = -1*(9+2) = -11

А13 = 1*1 = 1

А21 = -1*(-6) = 6

А22 = 1*(3-0) = 3

А23 = -1*(1+4) = -5

А31 = 1*2 = 2

А32 = -1*(-1) = -1

А33 = 1+(1+6) = 7

3. Следующий шаг - составление матрицы из получившихся дополнений:

5. Умножаем эту матрицу на число, обратное определителю, то есть на 1/26:

6. Ну а теперь нам просто нужно выполнить проверку:

В ходе проверки мы получили единичную матрицу, следовательно, решение было выполнено абсолютно верно.

2 способ вычисления обратной матрицы.

1. Элементарное преобразование матриц

2. Обратная матрица через элементарный преобразователь.

Элементарное преобразование матриц включает:

1. Умножение строки на число, не равное нулю.

2. Прибавление к любой строке другой строки, умноженной на число.

3. Перемена местами строк матрицы.

4. Применяя цепочку элементарных преобразований, получаем другую матрицу.

А-1 = ?

1. (A|E) ~ (E|A-1 )

2. A-1 * A = E

Рассмотрим это на практическом примере с действительными числами.

Задание: Найти обратную матрицу.

Решение:

Выполним проверку:

Небольшое разъяснение по решению:

Сперва мы переставили 1 и 2 строку матрицы, затем умножили первую строку на (-1).

После этого умножили первую строку на (-2) и сложили со второй строкой матрицы. После чего умножили 2 строку на 1/4.

Заключительным этапом преобразований стало умножение второй строки на 2 и прибавлением с первой. В результате слева у нас получилась единичная матрица, следовательно, обратная матрица - это матрица справа.

После проверки мы убедились в правильности решения.

Как вы видите, вычисление обратной матрицы - это очень просто.

В заключении данной лекции хотелось бы также уделить немного времени свойствам такой матрицы.

Нахождение обратной матрицы.

В этой статье разберемся с понятием обратной матрицы, ее свойствами и способами нахождения. Подробно остановимся на решении примеров, в которых требуется построить обратную матрицу для заданной.

Навигация по странице.

Обратная матрица - определение.

Нахождение обратной матрицы с помощью матрицы из алгебраических дополнений.

Свойства обратной матрицы.

Нахождение обратной матрицы методом Гаусса-Жордана.

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

Обратная матрица - определение.

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

Определение.

Матрица называется обратной для матрицы , определитель которой отличен от нуля , если справедливы равенства , где E – единичная матрица порядка n на n .

Нахождение обратной матрицы с помощью матрицы из алгебраических дополнений.

Как же находить обратную матрицу для данной?

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

Определение.

Минор k-ого порядка матрицы A порядка m на n – это определитель матрицы порядка k на k , которая получается из элементов матрицы А , находящихся в выбранныхk строках и k столбцах. (k не превосходит наименьшего из чисел m или n ).

Минор (n-1)-ого порядка, который составляется из элементов всех строк, кроме i-ой , и всех столбцов, кроме j-ого , квадратной матрицы А порядка n на n обозначим как .

Иными словами, минор получается из квадратной матрицы А порядка n на n вычеркиванием элементов i-ой строки и j-ого столбца.

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

Определение.

Алгебраическим дополнением элемента квадратной матрицы называют минор (n-1)-ого порядка, который получается из матрицы А , вычеркиванием элементов ее i-ой строки и j-ого столбца, умноженный на .

Алгебраическое дополнение элемента обозначается как . Таким обрзом, .

Например, для матрицы алгебраическое дополнение элемента есть .

Во-вторых, нам пригодятся два свойства определителя, которые мы разобрали в разделевычисление определителя матрицы :

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

Матрица действительно является обратной для матрицы А , так как выполняются равенства . Покажем это

Составим алгоритм нахождения обратной матрицы с использованием равенства .

Разберем алгоритм нахождения обратной матрицы на примере.

Пример.

Дана матрица . Найдите обратную матрицу.

Решение.

Вычислим определитель матрицы А , разложив его по элементам третьего столбца:

Определитель отличен от нуля, так что матрица А обратима.

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

Поэтому

Выполним транспонирование матрицы из алгебраических дополнений:

Теперь находим обратную матрицу как :

Проверяем полученный результат:

Равенства выполняются, следовательно, обратная матрица найдена верно.

Свойства обратной матрицы.

Понятие обратной матрицы, равенство , определения операций над матрицами и свойства определителя матрицы позволяют обосновать следующие свойства обратной матрицы :

Рассмотрим еще один способ нахождения обратной матрицы для квадратной матрицы А порядка n на n .

Этот метод основан на решении n систем линейных неоднородных алгебраических уравнений с n неизвестными. Неизвестными переменными в этих системах уравнений являются элементы обратной матрицы.

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

Приравнивая соответствующие элементы по столбцам, получим n систем линейных уравнений

Решаем их любым способом и из найденных значений составляем обратную матрицу.

Разберем этот метод на примере.

Пример.

Дана матрица . Найдите обратную матрицу.

Решение.

Примем . Равенство дает нам три системы линейных неоднородных алгебраических уравнений:

Не будем расписывать решение этих систем, при необходимости обращайтесь к разделурешение систем линейных алгебраических уравнений .

Из первой системы уравнений имеем , из второй - , из третьей - . Следовательно, искомая обратная матрица имеет вид . Рекомендуем сделать проверку, чтобы убедиться в правильности результата.

Подведем итог.

Мы рассмотрели понятие обратной матрицы, ее свойства и три метода ее нахождения.

Пример решений методом обратной матрицы

Задание 1. Решить СЛАУ методом обратной матрицы. 2 x 1 + 3x 2 + 3x 3 + x 4 = 1 3 x 1 + 5x 2 + 3x 3 + 2x 4 = 2 5 x 1 + 7x 2 + 6x 3 + 2x 4 = 3 4 x 1 + 4x 2 + 3x 3 + x 4 = 4

Начало формы

Конец формы

Решение . Запишем матрицу в виде: Вектор B: B T = (1,2,3,4) Главный определитель Минор для (1,1): = 5 (6 1-3 2)-7 (3 1-3 2)+4 (3 2-6 2) = -3 Минор для (2,1): = 3 (6 1-3 2)-7 (3 1-3 1)+4 (3 2-6 1) = 0 Минор для (3,1): = 3 (3 1-3 2)-5 (3 1-3 1)+4 (3 2-3 1) = 3 Минор для (4,1): = 3 (3 2-6 2)-5 (3 2-6 1)+7 (3 2-3 1) = 3 Определитель минора ∆ = 2 (-3)-3 0+5 3-4 3 = -3

Транспонированная матрица Алгебраические дополнения ∆ 1,1 = 5 (6 1-2 3)-3 (7 1-2 4)+2 (7 3-6 4) = -3 ∆ 1,2 = -3 (6 1-2 3)-3 (7 1-2 4)+1 (7 3-6 4) = 0 ∆ 1,3 = 3 (3 1-2 3)-3 (5 1-2 4)+1 (5 3-3 4) = 3 ∆ 1,4 = -3 (3 2-2 6)-3 (5 2-2 7)+1 (5 6-3 7) = -3 ∆ 2,1 = -3 (6 1-2 3)-3 (5 1-2 4)+2 (5 3-6 4) = 9 ∆ 2,2 = 2 (6 1-2 3)-3 (5 1-2 4)+1 (5 3-6 4) = 0 ∆ 2,3 = -2 (3 1-2 3)-3 (3 1-2 4)+1 (3 3-3 4) = -6 ∆ 2,4 = 2 (3 2-2 6)-3 (3 2-2 5)+1 (3 6-3 5) = 3 ∆ 3,1 = 3 (7 1-2 4)-5 (5 1-2 4)+2 (5 4-7 4) = -4 ∆ 3,2 = -2 (7 1-2 4)-3 (5 1-2 4)+1 (5 4-7 4) = 1 ∆ 3,3 = 2 (5 1-2 4)-3 (3 1-2 4)+1 (3 4-5 4) = 1 ∆ 3,4 = -2 (5 2-2 7)-3 (3 2-2 5)+1 (3 7-5 5) = 0 ∆ 4,1 = -3 (7 3-6 4)-5 (5 3-6 4)+3 (5 4-7 4) = -12 ∆ 4,2 = 2 (7 3-6 4)-3 (5 3-6 4)+3 (5 4-7 4) = -3 ∆ 4,3 = -2 (5 3-3 4)-3 (3 3-3 4)+3 (3 4-5 4) = 9 ∆ 4,4 = 2 (5 6-3 7)-3 (3 6-3 5)+3 (3 7-5 5) = -3 Обратная матрица Вектор результатов X X = A -1 ∙ B X T = (2,-1,-0.33,1) x 1 = 2 x 2 = -1 x 3 = -0.33 x 4 = 1

см. также решений СЛАУ методом обратной матрицы online. Для этого введите свои данные и получите решение с подробными комментариями.

Задание 2 . Систему уравнений записать в матричной форме и решить ее с помощью обратной матрицы. Сделать проверку полученного решения. Решение :xml :xls

Пример 2 . Записать систему уравнений в матричной форме и решить с помощью обратной матрицы. Решение :xml :xls

Пример . Дана система трех линейных уравнений с тремя неизвестными. Требуется: 1) найти ее решение с помощью формул Крамера ; 2) записать систему в матричной форме и решить ее средствами матричного исчисления. Методические рекомендации . После решения методом Крамера, найдите кнопку "Решение методом обратной матрицы для исходных данных". Вы получите соответствующее решение. Таким образом, данные вновь заполнять не придется. Решение . Обозначим через А - матрицу коэффициентов при неизвестных; X - матрицу-столбец неизвестных; B - матрицу-столбец свободных членов:

Вектор B: B T =(4,-3,-3) С учетом этих обозначений данная система уравнений принимает следующую матричную форму: А*Х = B. Если матрица А - невырожденная (ее определитель отличен от нуля, то она имеет обратную матрицу А -1 . Умножив обе части уравнения на А -1 , получим: А -1 *А*Х = А -1 *B, А -1 *А=Е. Это равенство называется матричной записью решения системы линейных уравнений . Для нахождения решения системы уравнений необходимо вычислить обратную матрицу А -1 . Система будет иметь решение, если определитель матрицы A отличен от нуля. Найдем главный определитель. ∆=-1 (-2 (-1)-1 1)-3 (3 (-1)-1 0)+2 (3 1-(-2 0))=14 Итак, определитель 14 ≠ 0, поэтому продолжаем решение. Для этого найдем обратную матрицу через алгебраические дополнения. Пусть имеем невырожденную матрицу А:

Вычисляем алгебраические дополнения.

∆ 1,1 =(-2 (-1)-1 1)=1

∆ 1,2 =-(3 (-1)-0 1)=3

∆ 1,3 =(3 1-0 (-2))=3

∆ 2,1 =-(3 (-1)-1 2)=5

∆ 2,2 =(-1 (-1)-0 2)=1

∆ 2,3 =-(-1 1-0 3)=1

∆ 3,1 =(3 1-(-2 2))=7

∆ 3,2 =-(-1 1-3 2)=7

X T =(-1,1,2) x 1 = -14 / 14 =-1 x 2 = 14 / 14 =1 x 3 = 28 / 14 =2 Проверка . -1 -1+3 1+0 2=4 3 -1+-2 1+1 2=-3 2 -1+1 1+-1 2=-3 doc :xml :xls Ответ: -1,1,2.

Похожие на обратные по многим свойствам.

Энциклопедичный YouTube

1 / 5

✪ Как находить обратную матрицу - bezbotvy

✪ Обратная матрица (2 способа нахождения)

✪ Обратная матрица #1

✪ 2015-01-28. Обратная матрица 3x3

✪ 2015-01-27. Обратная матрица 2х2

Субтитры

Свойства обратной матрицы

det A − 1 = 1 det A {\displaystyle \det A^{-1}={\frac {1}{\det A}}} , где det {\displaystyle \ \det } обозначает определитель .
(A B) − 1 = B − 1 A − 1 {\displaystyle \ (AB)^{-1}=B^{-1}A^{-1}} для двух квадратных обратимых матриц A {\displaystyle A} и B {\displaystyle B} .
(A T) − 1 = (A − 1) T {\displaystyle \ (A^{T})^{-1}=(A^{-1})^{T}} , где (. . .) T {\displaystyle (...)^{T}} обозначает транспонированную матрицу.
(k A) − 1 = k − 1 A − 1 {\displaystyle \ (kA)^{-1}=k^{-1}A^{-1}} для любого коэффициента k ≠ 0 {\displaystyle k\not =0} .
E − 1 = E {\displaystyle \ E^{-1}=E} .
Если необходимо решить систему линейных уравнений , (b - ненулевой вектор) где x {\displaystyle x} - искомый вектор, и если A − 1 {\displaystyle A^{-1}} существует, то x = A − 1 b {\displaystyle x=A^{-1}b} . В противном случае либо размерность пространства решений больше нуля, либо их нет вовсе.

Способы нахождения обратной матрицы

Если матрица обратима, то для нахождения обратной матрицы можно воспользоваться одним из следующих способов:

Точные (прямые) методы

Метод Гаусса-Жордана

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

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

Λ 1 ⋅ ⋯ ⋅ Λ n ⋅ A = Λ A = E ⇒ Λ = A − 1 {\displaystyle \Lambda _{1}\cdot \dots \cdot \Lambda _{n}\cdot A=\Lambda A=E\Rightarrow \Lambda =A^{-1}} . Λ m = [ 1 … 0 − a 1 m / a m m 0 … 0 … 0 … 1 − a m − 1 m / a m m 0 … 0 0 … 0 1 / a m m 0 … 0 0 … 0 − a m + 1 m / a m m 1 … 0 … 0 … 0 − a n m / a m m 0 … 1 ] {\displaystyle \Lambda _{m}={\begin{bmatrix}1&\dots &0&-a_{1m}/a_{mm}&0&\dots &0\\&&&\dots &&&\\0&\dots &1&-a_{m-1m}/a_{mm}&0&\dots &0\\0&\dots &0&1/a_{mm}&0&\dots &0\\0&\dots &0&-a_{m+1m}/a_{mm}&1&\dots &0\\&&&\dots &&&\\0&\dots &0&-a_{nm}/a_{mm}&0&\dots &1\end{bmatrix}}} .

Вторая матрица после применения всех операций станет равна Λ {\displaystyle \Lambda } , то есть будет искомой. Сложность алгоритма - O (n 3) {\displaystyle O(n^{3})} .

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

Матрица, обратная матрице A {\displaystyle A} , представима в виде

A − 1 = adj (A) det (A) {\displaystyle {A}^{-1}={{{\mbox{adj}}(A)} \over {\det(A)}}}

где adj (A) {\displaystyle {\mbox{adj}}(A)} - присоединенная матрица ;

Сложность алгоритма зависит от сложности алгоритма расчета определителя O det и равна O(n²)·O det .

Использование LU/LUP-разложения

Матричное уравнение A X = I n {\displaystyle AX=I_{n}} для обратной матрицы X {\displaystyle X} можно рассматривать как совокупность n {\displaystyle n} систем вида A x = b {\displaystyle Ax=b} . Обозначим i {\displaystyle i} -ый столбец матрицы X {\displaystyle X} через X i {\displaystyle X_{i}} ; тогда A X i = e i {\displaystyle AX_{i}=e_{i}} , i = 1 , … , n {\displaystyle i=1,\ldots ,n} ,поскольку i {\displaystyle i} -м столбцом матрицы I n {\displaystyle I_{n}} является единичный вектор e i {\displaystyle e_{i}} . другими словами, нахождение обратной матрицы сводится к решению n уравнений с одной матрицей и разными правыми частями. После выполнения LUP-разложения (время O(n³)) на решение каждого из n уравнений нужно время O(n²), так что и эта часть работы требует времени O(n³) .

Если матрица A невырождена, то для неё можно рассчитать LUP-разложение P A = L U {\displaystyle PA=LU} . Пусть P A = B {\displaystyle PA=B} , B − 1 = D {\displaystyle B^{-1}=D} . Тогда из свойств обратной матрицы можно записать: D = U − 1 L − 1 {\displaystyle D=U^{-1}L^{-1}} . Если умножить это равенство на U и L то можно получить два равенства вида U D = L − 1 {\displaystyle UD=L^{-1}} и D L = U − 1 {\displaystyle DL=U^{-1}} . Первое из этих равенств представляет собой систему из n² линейных уравнений для n (n + 1) 2 {\displaystyle {\frac {n(n+1)}{2}}} из которых известны правые части (из свойств треугольных матриц). Второе представляет также систему из n² линейных уравнений для n (n − 1) 2 {\displaystyle {\frac {n(n-1)}{2}}} из которых известны правые части (также из свойств треугольных матриц). Вместе они представляют собой систему из n² равенств. С помощью этих равенств можно реккурентно определить все n² элементов матрицы D. Тогда из равенства (PA) −1 = A −1 P −1 = B −1 = D. получаем равенство A − 1 = D P {\displaystyle A^{-1}=DP} .

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

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

Итерационные методы

Методы Шульца

{ Ψ k = E − A U k , U k + 1 = U k ∑ i = 0 n Ψ k i {\displaystyle {\begin{cases}\Psi _{k}=E-AU_{k},\\U_{k+1}=U_{k}\sum _{i=0}^{n}\Psi _{k}^{i}\end{cases}}}

Оценка погрешности

Выбор начального приближения

Проблема выбора начального приближения в рассматриваемых здесь процессах итерационного обращения матриц не позволяет относиться к ним как к самостоятельным универсальным методам, конкурирующими с прямыми методами обращения, основанными, например, на LU-разложении матриц. Имеются некоторые рекомендации по выбору U 0 {\displaystyle U_{0}} , обеспечивающие выполнение условия ρ (Ψ 0) < 1 {\displaystyle \rho (\Psi _{0})<1} (спектральный радиус матрицы меньше единицы), являющегося необходимым и достаточным для сходимости процесса. Однако при этом, во-первых, требуется знать сверху оценку спектра обращаемой матрицы A либо матрицы A A T {\displaystyle AA^{T}} (а именно, если A - симметричная положительно определённая матрица и ρ (A) ≤ β {\displaystyle \rho (A)\leq \beta } , то можно взять U 0 = α E {\displaystyle U_{0}={\alpha }E} , где ; если же A - произвольная невырожденная матрица и ρ (A A T) ≤ β {\displaystyle \rho (AA^{T})\leq \beta } , то полагают U 0 = α A T {\displaystyle U_{0}={\alpha }A^{T}} , где также α ∈ (0 , 2 β) {\displaystyle \alpha \in \left(0,{\frac {2}{\beta }}\right)} ; можно конечно упростить ситуацию и, воспользовавшись тем, что ρ (A A T) ≤ k A A T k {\displaystyle \rho (AA^{T})\leq {\mathcal {k}}AA^{T}{\mathcal {k}}} , положить U 0 = A T ‖ A A T ‖ {\displaystyle U_{0}={\frac {A^{T}}{\|AA^{T}\|}}} ). Во-вторых, при таком задании начальной матрицы нет гарантии, что ‖ Ψ 0 ‖ {\displaystyle \|\Psi _{0}\|} будет малой (возможно, даже окажется ‖ Ψ 0 ‖ > 1 {\displaystyle \|\Psi _{0}\|>1} ), и высокий порядок скорости сходимости обнаружится далеко не сразу.

Примеры

Матрица 2х2

A − 1 = [ a b c d ] − 1 = 1 det (A) [ d − b − c a ] = 1 a d − b c [ d − b − c a ] . {\displaystyle \mathbf {A} ^{-1}={\begin{bmatrix}a&b\\c&d\\\end{bmatrix}}^{-1}={\frac {1}{\det(\mathbf {A})}}{\begin{bmatrix}\,\,\,d&\!\!-b\\-c&\,a\\\end{bmatrix}}={\frac {1}{ad-bc}}{\begin{bmatrix}\,\,\,d&\!\!-b\\-c&\,a\\\end{bmatrix}}.}

Обращение матрицы 2х2 возможно только при условии, что a d − b c = det A ≠ 0 {\displaystyle ad-bc=\det A\neq 0} .