Обратная матрица. Примеры вычисления
Нахождение обратной матрицы является важной составляющей в разделе линейной алгебры. С помощью таких матриц, если они существуют, можно быстро найти решение системы линейных уравнений.
Матрицаназывается обратной к матрице,если выполняются следующие равенства.
.
Если определитель матрицыотличен от нуля, то матрицу называют не особо или невырожденной.
Для того, чтобы матрица имела обратную необходимо и достаточно, чтобы она была невырожденной
Алгоритм нахождения обратной матрицы
Пусть имеем квадратную матрицу
и нужно найти обратную к ней. Для этого нужно выполнить следующие действия:
1. Найти определитель матрицы. Если он не равен нулю то выполняем следующие действия. В противном случае данная матрица вырождена и для нее не существует обратной
2. Найти алгебраические дополнения элементов матрицы . Они равны минорам, умноженным на в степени суммы строки и столбца, для которого ищем.
3. Составить матрицу из алгебраических дополнений элементов матрицы матрицы и протранспонировать ее. Эта матрица называется присоединенной или союзной и обозначается .
4. Разделить присоединенную матрицу на детерминант . Полученная матрица будет обратной и иметь свойства, которые изложены в начале статьи.
———————————————
Пример 1.
Найти матрицу, обратную к матрице (Дубовик В.П., Юрик И.И. «Высшая математика. Сборник задач»)
1) (1.127)
2) (1.130)
3) (1.133)
Решение.
1)Находим определитель матрицы
Так как детерминант не равен нулю (), то обратная матрица существует. Находим матрицу, составленную из алгебраических дополнений
Матрица дополнений примет вид
Транспонируем ее и получаем присоединенную
Разделим ее на определитель и получим обратную
Видим, что в случае, когда определитель равен единице присоединена и обратная матрицы совпадают.
2) Вычисляем определитель матрицы
Находим матрицу алгебраических дополнений
Конечный вид матрицы дополнений
Транспонируем ее и находим союзную матрицу
Находим обратную матрицу
3) Вычислим детерминант матрицы. Для этого разложим его на первую строчку. В результате получим два отличны от нуля слагаемые
Находим матрицу алгебраических дополнений. Расписание определителя проводим по строкам и столбцам, в которых больше нулевых элементов (обозначены черным цветом).
Конечный вид матрицы дополнений следующий
Тран
Нахождение обратной матрицы: три алгоритма и примеры
Нахождение обратной матрицы — процесс, который состоит из достаточно простых действий. Но эти действия повторяются так часто, что процесс получается довольно продолжительным. Главное — не потерять внимание при решении.
При решении наиболее распространённым методом — алгебраических дополнений — потребуется:
При решении примеров мы разберём эти действия подробнее. А пока узнаем, что гласит теория об обратной матрице.
Для обратной матрицы существует уместная аналогия с обратным числом. Для каждого числа a, не равного нулю, существует такое число b, что произведение a и b равно единице: ab = 1. Число b называется обратным для числа b. Например, для числа 7 обратным является число 1/7, так как 7*1/7=1.
Обратной матрицей, которую требуется отыскать для данной квадратной матрицы А, называется такая матрица
,
произведение на которую матрицы А справа является единичной матрицей, т.е,
. (1)
Единичной матрицей называется диагональная матрица, у которой все диагональные элементы равны единице.
Нахождение обратной матрицы — задача, которая чаще решается двумя методами:
- методом алгебраических дополнений, при котором, как было замечено в начале урока, требуется находить определители, миноры и алгебраические дополнения и транспонировать матрицы;
- методом исключения неизвестных Гаусса, при котором требуется производить элементарные преобразования матриц (складывать строки, умножать строки на одно и то же число и т. д.).
Для особо любознательных существуют и другие методы, например, метод линейных преобразований. На этом уроке разберём три упомянутых метода и алгоритмы нахождения обратной матрицы этими методами.
Теорема. Для каждой неособенной (невырожденной, несингулярной) квадратной матрицы можно найти обратную матрицу, и притом только одну. Для особенной (вырожденной, сингулярной) квадратной матрицы обратная матрица не существует.
Квадратная матрица называется неособенной (или
Обратная матрица может быть найдена только для квадратной матрицы. Естественно, обратная матрица также будет квадратной и того же порядка, что и данная матрица. Матрица, для которой может быть найдена обратная матрица, называется обратимой матрицей.
На сайте есть онлайн калькулятор для нахождения обратной матрицы. Вы можете открыть его в новом окне уже сейчас, если держите перед собой ваши собственные задания. А мы разберём несколько разминочных.
Для неособенной квадратной матрицы А обратной является матрица
, (2)где — определитель матрицы А, а — матрица, союзная с матрицей А.
Разберём ключевые понятия, которые потребуются для решения задач — союзная матрица, алгебраические дополнения и транспонированная матрица.
Пусть существует квадратная матрица A:
Транспонированная относительно матрицы A матрица
Остановимся на минорах и алгебраических дополнениях.
Пусть есть квадратная матрица третьего порядка:
.
Её определитель:
Вычислим алгебраическое дополнение элемента
, то есть элемента 2, стоящего на пересечении первой строки и второго столбца.Для этого нужно сначала найти минор этого элемента. Он получается вычёркиванием из определителя строки и столбца, на пересечении которых стоит указанный элемент. В результате останется следующий определитель, который и является минором элемента :
.
Алгебраическое дополнение элемента
.
По этой инструкции нужно вычислить алгебраические дополнения всех элементов матрицы A’, транспонированной относительно матрицы матрица A.
И последнее из значимых для нахождение обратной матрицы понятий. Союзной с квадратной матрицей A называется матрица того же порядка, элементами которой являются алгебраические дополнения соответствующих элементов определителя матрицы , транспонированной относительно матрицы A. Таким образом, союзная матрица состоит из следующих элементов:Алгоритм нахождения обратной матрицы методом алгебраических дополнений
2. Найти матрицу, транспонированную относительно A.
3. Вычислить элементы союзной матрицы как алгебраические дополнения марицы, найденной на шаге 2.
4. Применить формулу (2): умножить число, обратное определителю матрицы A, на союзную матрицу, найденную на шаге 4.
5. Проверить полученный на шаге 4 результат, умножив данную матрицу A на обратную матрицу. Если произведение этих матриц равно единичной матрицы, значит обратная матрица была найдена верно. В противном случае начать процесс решения снова.
Пример 1. Для матрицы
найти обратную матрицу.
Решение. Для нахождения обратной матрицы необходимо найти определитель матрицы А . Находим по правилу треугольников:
Следовательно, матрица А – неособенная (невырожденная, несингулярная) и для неё существует обратная.
Найдём матрицу, союзную с данной матрицей А.
Найдём матрицу , транспонированную относительно матрицы A:
Вычисляем элементы союзной матрицы как алгебраические дополнения матрицы, транспонированной относительно матрицы A:
Следовательно, матрица , союзная с матрицей A, имеет вид
Замечание. Порядок вычисления элементов и транспонирования матрицы может быть иным. Можно сначала вычислить алгебраические дополнения матрицы A, а затем транспонировать матрицу алгебраических дополнений. В результате должны получиться те же элементы союзной матрицы.
Применяя формулу (2), находим матрицу, обратную матрице А:
Проверить решение можно с помощью онлайн калькулятора для нахождения обратной матрицы.
Алгоритм нахождения обратной матрицы методом исключения неизвестных Гаусса
1. К матрице A приписать единичную матрицу того же порядка.
2. Полученную сдвоенную матрицу преобразовать так, чтобы в левой её части получилась единичная матрица, тогда в правой части на месте единичной матрицы автоматически получится обратная матрица. Матрица A в левой части преобразуется в единичную матрицу путём элементарных преобразований матрицы.
2. Если в процессе преобразования матрицы A в единичную матрицу в какой-либо строке или в каком-либо столбце окажутся только нули, то определитель матрицы равен нулю, и, следовательно, матрица A будет вырожденной, и она не имеет обратной матрицы. В этом случае дальнейшее нахождение обратной матрицы прекращается.
Пример 2. Для матрицы
найти обратную матрицу.
Решение. Составляем сдвоенную матрицу
и будем её преобразовывать, так чтобы в левой части получилась единичная матрица. Начинаем преобразования.
Умножим первую строку левой и правой матрицы на (-3) и сложим её со второй строкой, а затем умножим первую строку на (-4) и сложим её с третьей строкой, тогда получим
.
Чтобы по возможности не было дробных чисел при последующих преобразованиях, создадим предварительно единицу во второй строке в левой части сдвоенной матрицы. Для этого умножим вторую строку на 2 и вычтем из неё третью строку, тогда получим
.
Сложим первую строку со второй, а затем умножим вторую строку на (-9) и сложим её с третьей строкой. Тогда получим
.
Разделим третью строку на 8, тогда
.
Умножим третью строку на 2 и сложим её со второй строкой. Получается:
.
Переставим местами вторую и третью строку, тогда окончательно получим:
.
Видим, что в левой части получилась единичная матрица, следовательно, в правой части получилась обратная матрица . Таким образом:
.
Можно проверить правильность вычислений, умножим исходную матрицу на найденную обратную матрицу:
.
В результате должна получиться обратная матрица.
Проверить решение можно с помощью онлайн калькулятора для нахождения обратной матрицы.
Пример 3. Для матрицы
найти обратную матрицу.
Решение. Составляем сдвоенную матрицу
и будем её преобразовывать.
Первую строку умножаем на 3, а вторую на 2, и вычитаем из второй, а затем первую строку умножаем на 5, а третью на 2 и вычитаем из третьей строки, тогда получим
.
Первую строку умножаем на 2 и складываем её со второй, а затем из третьей строки вычитаем вторую, тогда получим
.
Видим, что в третьей строке в левой части все элементы получились равными нулю. Следовательно, матрица вырожденная и обратной матрицы не имеет. Дальнейшее нахождение обратной марицы прекращаем.
Проверить решение можно с помощью онлайн калькулятора для нахождения обратной матрицы.
Матрицы теснейшим образом связаны с системами линейных уравнений. Каждой матрице соответствует система линейных уравнений, коэффициенты в которой есть элементы матрицы. И наоборот, системе линейных уравнений соответствует некоторая матрица.
Поэтому существует метод линейных преобразований для нахождения обратной матрицы. Для решения задач нам будет достаточно знать, что линейное преобразование — это система линейных уравнений, вид которой будет приведён ниже в алгоритме.
Алгоритм нахождения обратной матрицы методом линейных преобразований
1. Для данной невырожденной матрицы A составить линейное преобразование — систему линейных уравнений вида
,
где aij — элементы матрицы A.
2. Решить полученную систему относительно y — найти для предыдущего линейного преобразование обратное линейное преобразование
,
в котором Aij — алгебраические дополнения элементов матрицы A, Δ — определитель матрицы A. Внимание! Алгебраические дополнения располагаются как в транспонированной матрице, то есть для элементов строки — в столбце, а для элементов столбца — в строке.
3. Находим коэффициенты при y: , которые и будут элементами матрицы, обратной для матрицы A.
4. Пользуясь элементами, найденными на шаге 3, записать найденную обратную матрицу.
Наиболее наблюдательные могли заметить, что по сути метод линейных преобразований — это тот же метод алгебраических преобразований (союзной матрицы), но с другой формой записи. Для кого-то метод линейных преобразований может оказаться более удобным как более компактный.
Пример 4. Найти обратную матрицу для матрицы
.
Сначала проверим, не равен ли нулю определитель данной матрицы. Он не равен нулю, следовательно, обратная матрица существует.
Для данной матрицы записываем линейное преобразование:
.
Находим линейное преобразование, обратное предыдущему, для этого потребуется находить алгебраические дополнения (урок откроется в новом окне). Запишем обратное линейное преобразование:
Коэффициенты при иксах в обратном линейном преобразовании — это элементы обратной матрицы для матрицы A. Таким образом нашли обратную матрицу:
Проверить решение можно с помощью онлайн калькулятора для нахождения обратной матрицы.
Начало темы «Матрицы»
Другие темы линейной алгебры
Нахождение обратной матрицы с примеры решения. Описание основных методов вычисления обратной матрицы
Задание. Найти обратную матрицу к матрице $ A=\left( \begin{array}{rrr}{1} & {0} & {2} \\ {2} & {-1} & {1} \\ {1} & {3} & {-1}\end{array}\right) $
Решение. Вычисляем определитель матрицы:
$$ \Delta=\left| \begin{array}{rrr}{1} & {0} & {2} \\ {2} & {-1} & {1} \\ {1} & {3} & {-1}\end{array}\right|=1 \cdot(-1) \cdot(-1)+2 \cdot 3 \cdot 2+0 \cdot 1 \cdot 1- $$
$$ -1 \cdot(-1) \cdot 2-3 \cdot 1 \cdot 1-2 \cdot 0 \cdot(-1)=1+12+0+2-3+0=12 \neq 0 $$
Так как определитель не равен нулю, то матрица имеет обратную. Обратная матрица $A^{-1}$ к матрице $A$ находится по формуле:
$$ A^{-1}=\frac{1}{\Delta} \cdot \widetilde{A}^{T} $$
Найдем союзную матрицу $ \tilde{A} $ , для этого вычислим алгебраические дополнения к элементам матрицы $A$ :
$$ A_{11}=(-1)^{1+1} \left| \begin{array}{rr}{-1} & {1} \\ {3} & {-1}\end{array}\right|=(-1) \cdot(-1)-3 \cdot 1=1-3=-2 $$
$$ A_{12}=(-1)^{1+2} \left| \begin{array}{rr}{2} & {1} \\ {1} & {-1}\end{array}\right|=-[2 \cdot(-1)-1 \cdot 1]=-(-2-1)=3 $$
$$ A_{13}=(-1)^{1+3} \left| \begin{array}{rr}{2} & {-1} \\ {1} & {3}\end{array}\right|=2 \cdot 3-1 \cdot(-1)=6+1=7 $$
$$ A_{21}=(-1)^{2+1} \left| \begin{array}{rr}{0} & {2} \\ {3} & {-1}\end{array}\right|=-[0 \cdot(-1)-3 \cdot 2]=-(0-6)=6 $$
$$ A_{22}=(-1)^{2+2} \left| \begin{array}{rr}{1} & {2} \\ {1} & {-1}\end{array}\right|=1 \cdot(-1)-1 \cdot 2=-1-2=-3 $$
$$ A_{23}=(-1)^{2+3} \left| \begin{array}{cc}{1} & {0} \\ {1} & {3}\end{array}\right|=-[1 \cdot 3-1 \cdot 0]=-(3-0)=-3 $$
$$ A_{31}=(-1)^{3+1} \left| \begin{array}{rr}{0} & {2} \\ {-1} & {1}\end{array}\right|=0 \cdot 1-(-1) \cdot 2=0+2=2 $$
$$ A_{32}=(-1)^{3+2} \left| \begin{array}{cc}{1} & {2} \\ {2} & {1}\end{array}\right|=-[1 \cdot 1-2 \cdot 2]=-(1-4)=3 $$
$$ A_{33}=(-1)^{3+3} \left| \begin{array}{rr}{1} & {0} \\ {2} & {-1}\end{array}\right|=1 \cdot(-1)-2 \cdot 0=-1-0=-1 $$
Таким образом, $ \tilde{A}=\left( \begin{array}{rrr}{-2} & {3} & {7} \\ {6} & {-3} & {-3} \\ {2} & {3} & {-1}\end{array}\right) $
Транспонируем эту матрицу (т.е. строки матрицы делаем столбцами с тем же номером):
$$ \widetilde{A}^{T}=\left( \begin{array}{rrr}{-2} & {6} & {2} \\ {3} & {-3} & {3} \\ {7} & {-3} & {-1}\end{array}\right) $$
Итак, $ A^{-1}=\frac{1}{12} \left( \begin{array}{rrr}{-2} & {6} & {2} \\ {3} & {-3} & {3} \\ {7} & {-3} & {-1}\end{array}\right) $
Ответ. $ A^{-1}=\frac{1}{12} \left( \begin{array}{rrr}{-2} & {6} & {2} \\ {3} & {-3} & {3} \\ {7} & {-3} & {-1}\end{array}\right) $
Обратная матрица — Википедия
Обра́тная ма́трица — такая матрица A−1, при умножении на которую исходная матрица A даёт в результате единичную матрицу E:
- AA−1=A−1A=E{\displaystyle AA^{-1}=A^{-1}A=E}
Квадратная матрица обратима тогда и только тогда, когда она невырождена, то есть её определитель не равен нулю. Для неквадратных матриц и вырожденных матриц обратных матриц не существует. Однако возможно обобщить это понятие и ввести псевдообратные матрицы, похожие на обратные по многим свойствам.
Способы нахождения обратной матрицы[править | править код]
Если матрица обратима, то для нахождения обратной матрицы можно воспользоваться одним из следующих способов:
Точные (прямые) методы[править | править код]
Метод Жордана—Гаусса[править | править код]
Возьмём две матрицы: саму 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−a1m/amm0…0…0…1−am−1m/amm0…00…01/amm0…00…0−am+1m/amm1…0…0…0−anm/amm0…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(n3){\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)} — присоединенная матрица (матрица, составленная из алгебраических дополнений для соответствующих элементов транспонированной матрицы).
Сложность алгоритма зависит от сложности алгоритма расчета определителя Odet и равна O(n²)·Odet.
Использование LU/LUP-разложения[править | править код]
Матричное уравнение AX=In{\displaystyle AX=I_{n}} для обратной матрицы X{\displaystyle X} можно рассматривать как совокупность n{\displaystyle n} систем вида Ax=b{\displaystyle Ax=b}. Обозначим i{\displaystyle i}-й столбец матрицы X{\displaystyle X} через Xi{\displaystyle X_{i}}; тогда AXi=ei{\displaystyle AX_{i}=e_{i}}, i=1,…,n{\displaystyle i=1,\ldots ,n}, поскольку i{\displaystyle i}-м столбцом матрицы In{\displaystyle I_{n}} является единичный вектор ei{\displaystyle e_{i}}. другими словами, нахождение обратной матрицы сводится к решению n уравнений с одной матрицей и разными правыми частями. После выполнения LUP-разложения (время O(n³)) на решение каждого из n уравнений нужно время O(n²), так что и эта часть работы требует времени O(n³)[1].
Если матрица A невырождена, то для неё можно рассчитать LUP-разложение PA=LU{\displaystyle PA=LU}. Пусть PA=B{\displaystyle PA=B}, B−1=D{\displaystyle B^{-1}=D}. Тогда из свойств обратной матрицы можно записать: D=U−1L−1{\displaystyle D=U^{-1}L^{-1}}. Если умножить это равенство на U и L то можно получить два равенства вида UD=L−1{\displaystyle UD=L^{-1}} и DL=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−1P−1 = B−1 = D получаем равенство A−1=DP{\displaystyle A^{-1}=DP}.
В случае использования LU-разложения не требуется перестановки столбцов матрицы D, но решение может разойтись даже если матрица A невырождена.
Сложность алгоритма — O(n³).
Итерационные методы[править | править код]
Методы Шульца[править | править код]
{Ψk=E−AUk,Uk+1=Uk∑i=0nΨki{\displaystyle {\begin{cases}\Psi _{k}=E-AU_{k},\\U_{k+1}=U_{k}\sum _{i=0}^{n}\Psi _{k}^{i}\end{cases}}}
Оценка погрешности[править | править код]
Выбор начального приближения[править | править код]
Проблема выбора начального приближения U0{\displaystyle U_{0}} в рассматриваемых здесь процессах итерационного обращения матриц не позволяет относиться к ним как к самостоятельным универсальным методам, конкурирующими с прямыми методами обращения, основанными, например, на LU-разложении матриц. Имеются некоторые рекомендации по выбору U0{\displaystyle U_{0}}, обеспечивающие выполнение условия ρ(Ψ0)<1{\displaystyle \rho (\Psi _{0})<1} (спектральный радиус матрицы меньше единицы), являющегося необходимым и достаточным для сходимости процесса. Однако при этом, во-первых, требуется знать сверху оценку спектра обращаемой матрицы A либо матрицы AAT{\displaystyle AA^{T}} (а именно, если A — симметричная положительно определённая матрица и ρ(A)≤β{\displaystyle \rho (A)\leq \beta }, то можно взять U0=αE{\displaystyle U_{0}={\alpha }E}, где α∈(0,2β){\displaystyle \alpha \in \left(0,{\frac {2}{\beta }}\right)}; если же A — произвольная невырожденная матрица и ρ(AAT)≤β{\displaystyle \rho (AA^{T})\leq \beta }, то полагают U0=αAT{\displaystyle U_{0}={\alpha }A^{T}}, где также α∈(0,2β){\displaystyle \alpha \in \left(0,{\frac {2}{\beta }}\right)}; можно конечно упростить ситуацию и, воспользовавшись тем, что ρ(AAT)≤kAATk{\displaystyle \rho (AA^{T})\leq {\mathcal {k}}AA^{T}{\mathcal {k}}}, положить U0=AT‖AAT‖{\displaystyle U_{0}={\frac {A^{T}}{\|AA^{T}\|}}}). Во-вторых, при таком задании начальной матрицы нет гарантии, что ‖Ψ0‖{\displaystyle \|\Psi _{0}\|} будет малой (возможно, даже окажется ‖Ψ0‖>1{\displaystyle \|\Psi _{0}\|>1}), и высокий порядок скорости сходимости обнаружится далеко не сразу.
Матрица 2 × 2[править | править код]
- A−1=[abcd]−1=1detA[d−b−ca]=1ad−bc[d−b−ca]{\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 × 2 возможно только при условии, что ad−bc=detA≠0{\displaystyle ad-bc=\det A\neq 0}.
Обратная матрица и методы ее вычисления
Запишем вспомогательную матрицу
и приведем её, с помощью элементарных преобразований, к матрице, в которой единичная матрица будет слева. Переставим местами первую и вторую строки
Прибавим ко второй строке первую строку, умноженную на а к третьей строке первую, умноженную на
Прибавим ко второй строке третью, умноженную на
Умножим вторую строку на
Прибавим к первой строке вторую, умноженную на а к третьей вторую, умноженную на
Разделим третью строку на 3
К первой строке прибавим третью, умноженную на
Тогда обратная матрица равна
1.3.4. Примеры решения задач по теме «Обратная матрица»
Задача 1.
Найти обратную матрицу для матрицы
И проверить выполнение условий А А-1 = А-1А = Е.
Указание
Убедитесь, что матрица А – невырожденная, и примените способ вычисления обратной матрицы.
Решение
Убедимся, что матрица А – невырожденная. ΔА = 1·4 — 2·(-1) ≠ 0, следовательно, А-1 существует.
Вычислим алгебраические дополнения к элементам А:
Применим способ вычисления обратной матрицы:
.
Не забудьте, что обратная матрица образована из алгебраических дополнений к элементам Транспонированной матрицы! |
Найдем произведения А А-1 и А-1А:
Таким образом, найденная матрица А-1 отвечает определению обратной матрицы.
Ответ: .
Задача 2.
Найти обратную матрицу для матрицы
.
Указание
Убедитесь, что матрица А – невырожденная, и примените способ вычисления обратной матрицы.
Решение
Следовательно, матрица А невырожденная, и обратная матрица существует.
Вычислим алгебраические дополнения к элементам матрицы А:
Обратная матрица имеет вид:
Ответ: .
Задача 3.
Найти обратную матрицу для матрицы
.
Указание
Убедитесь, что матрица А – невырожденная, и примените способ вычисления обратной матрицы.
Решение
Вычислим определитель матрицы А разложением по первому столбцу:
.
Следовательно, обратная матрица для матрицы А существует.
Найдем алгебраические дополнения к элементам матрицы А:
Значит,
.
Ответ: .
Задача 4.
Найти обратную матрицу для матрицы
.
Указание
Убедитесь, что матрица А – невырожденная, и примените способ вычисления обратной матрицы.
Решение
.
Ответ:
Задача 5.
При каких X, Y, Z матрица
Является обратной к матрице
Указание
Необходимым условием того, что В = А-1, является требование АВ = Е.
Решение
Проверим невырожденность матрицы А:
Необходимым условием того, что В = А-1, является требование АВ = Е.
Найдем АВ:
Для того, чтобы выполнялось условие АВ = Е, X, Y, Z должны быть решением системы уравнений
Проверим, будет ли равно единичной матрице произведение ВА:
Значит, при найденных значениях X, Y, Z В = А-1.
Ответ: X = -3, Y = -3, Z = 4.
< Предыдущая | Следующая > |
---|
Обратная матрица
Эта тема является одной из самых ненавистных среди студентов. Хуже, наверное, только определители.
Фишка в том, что само понятие обратного элемента (и я сейчас не только о матрицах) отсылает нас к операции умножения. Даже в школьной программе умножение считается сложной операцией, а уж умножение матриц — вообще отдельная тема, которой у меня посвящён целый параграф и видеоурок.
Сегодня мы не будем вдаваться в подробности матричных вычислений. Просто вспомним: как обозначаются матрицы, как они умножаются и что из этого следует.
Повторение: умножение матриц
Прежде всего договоримся об обозначениях. Матрицей $A$ размера $\left[ m\times n \right]$ называется просто таблица из чисел, в которой ровно $m$ строк и $n$ столбцов:
\[A=\left[ m\times n \right]=\underbrace{\left[ \begin{matrix} {{a}_{11}} & {{a}_{12}} & … & {{a}_{1n}} \\ {{a}_{21}} & {{a}_{22}} & … & {{a}_{2n}} \\ … & … & … & … \\ {{a}_{m1}} & {{a}_{m2}} & … & {{a}_{mn}} \\\end{matrix} \right]}_{n}\]
Чтобы случайно не перепутать строки и столбцы местами (поверьте, на экзамене можно и единицу с двойкой перепутать — что уж говорить про какие-то там строки), просто взгляните на картинку:
Определение индексов для клеток матрицыЧто происходит? Если разместить стандартную систему координат $OXY$ в левом верхнем углу и направить оси так, чтобы они охватывали всю матрицу, то каждой клетке этой матрицы можно однозначно сопоставить координаты $\left( x;y \right)$ — это и будет номер строки и номер столбца.
Почему система координат размещена именно в левом верхнем углу? Да потому что именно оттуда мы начинаем читать любые тексты. Это очень просто запомнить.
А почему ось $x$ направлена именно вниз, а не вправо? Опять всё просто: возьмите стандартную систему координат (ось $x$ идёт вправо, ось $y$ — вверх) и поверните её так, чтобы она охватывала матрицу. Это поворот на 90 градусов по часовой стрелке — его результат мы и видим на картинке.
В общем, как определять индексы у элементов матрицы, мы разобрались. Теперь давайте разберёмся с умножением.
Определение. Матрицы $A=\left[ m\times n \right]$ и $B=\left[ n\times k \right]$, когда количество столбцов в первой совпадает с количеством строк во второй, называются согласованными.
Именно в таком порядке. Можно сумничать и сказать, мол, матрицы $A$ и $B$ образуют упорядоченную пару $\left( A;B \right)$: если они согласованы в таком порядке, то совершенно необязательно, что $B$ и $A$, т.е. пара $\left( B;A \right)$ — тоже согласована.
Умножать можно только согласованные матрицы.
Определение. Произведение согласованных матриц $A=\left[ m\times n \right]$ и $B=\left[ n\times k \right]$ — это новая матрица $C=\left[ m\times k \right]$, элементы которой ${{c}_{ij}}$ считаются по формуле:
\[{{c}_{ij}}=\sum\limits_{k=1}^{n}{{{a}_{ik}}}\cdot {{b}_{kj}}\]
Другими словами: чтобы получить элемент ${{c}_{ij}}$ матрицы $C=A\cdot B$, нужно взять $i$-строку первой матрицы, $j$-й столбец второй матрицы, а затем попарно перемножить элементы из этой строки и столбца. Результаты сложить.
Да, вот такое суровое определение. Из него сразу следует несколько фактов:
- Умножение матриц, вообще говоря, некоммутативно: $A\cdot B\ne B\cdot A$;
- Однако умножение ассоциативно: $\left( A\cdot B \right)\cdot C=A\cdot \left( B\cdot C \right)$;
- И даже дистрибутивно: $\left( A+B \right)\cdot C=A\cdot C+B\cdot C$;
- И ещё раз дистрибутивно: $A\cdot \left( B+C \right)=A\cdot B+A\cdot C$.
Дистрибутивность умножения пришлось отдельно описывать для левого и правого множителя-суммы как раз из-за некоммутативности операции умножения.
Если всё же получается так, что $A\cdot B=B\cdot A$, такие матрицы называются перестановочными.
Среди всех матриц, которые там на что-то умножаются, есть особые — те, которые при умножении на любую матрицу $A$ снова дают $A$:
Определение. Матрица $E$ называется единичной, если $A\cdot E=A$ или $E\cdot A=A$. В случае с квадратной матрицей $A$ можем записать:
\[A\cdot E=E\cdot A=A\]
Единичная матрица — частый гость при решении матричных уравнений. И вообще частый гость в мире матриц.:)
А ещё из-за этой $E$ кое-кто придумал всю ту дичь, которая будет написана дальше.
Что такое обратная матрица
Поскольку умножение матриц — весьма трудоёмкая операция (приходится перемножать кучу строчек и столбцов), то понятие обратной матрицы тоже оказывается не самым тривиальным. И требующим некоторых пояснений.
Ключевое определение
Что ж, пора познать истину.
Определение. Матрица $B$ называется обратной к матрице $A$, если
\[A\cdot B=B\cdot A=E\]
Обратная матрица обозначается через ${{A}^{-1}}$ (не путать со степенью!), поэтому определение можно переписать так:
\[A\cdot {{A}^{-1}}={{A}^{-1}}\cdot A=E\]
Казалось бы, всё предельно просто и ясно. Но при анализе такого определения сразу возникает несколько вопросов:
- Всегда ли существует обратная матрица? И если не всегда, то как определить: когда она существует, а когда — нет?
- А кто сказал, что такая матрица ровно одна? Вдруг для некоторой исходной матрицы $A$ найдётся целая толпа обратных?
- Как выглядят все эти «обратные»? И как, собственно, их считать?
Насчёт алгоритмов вычисления — об этом мы поговорим чуть позже. Но на остальные вопросы ответим прямо сейчас. Оформим их в виде отдельных утверждений-лемм.
Основные свойства
Начнём с того, как в принципе должна выглядеть матрица $A$, чтобы для неё существовала ${{A}^{-1}}$. Сейчас мы убедимся в том, что обе эти матрицы должны быть квадратными, причём одного размера: $\left[ n\times n \right]$.
Лемма 1. Дана матрица $A$ и обратная ей ${{A}^{-1}}$. Тогда обе эти матрицы — квадратные, причём одинакового порядка $n$.
Доказательство. Всё просто. Пусть матрица $A=\left[ m\times n \right]$, ${{A}^{-1}}=\left[ a\times b \right]$. Поскольку произведение $A\cdot {{A}^{-1}}=E$ по определению существует, матрицы $A$ и ${{A}^{-1}}$ согласованы в указанном порядке:
\[\begin{align} & \left[ m\times n \right]\cdot \left[ a\times b \right]=\left[ m\times b \right] \\ & n=a \end{align}\]
Это прямое следствие из алгоритма перемножения матриц: коэффициенты $n$ и $a$ являются «транзитными» и должны быть равны.
Вместе с тем определено и обратное умножение: ${{A}^{-1}}\cdot A=E$, поэтому матрицы ${{A}^{-1}}$ и $A$ тоже согласованы в указанном порядке:
\[\begin{align} & \left[ a\times b \right]\cdot \left[ m\times n \right]=\left[ a\times n \right] \\ & b=m \end{align}\]
Таким образом, без ограничения общности можем считать, что $A=\left[ m\times n \right]$, ${{A}^{-1}}=\left[ n\times m \right]$. Однако согласно определению $A\cdot {{A}^{-1}}={{A}^{-1}}\cdot A$, поэтому размеры матриц строго совпадают:
\[\begin{align} & \left[ m\times n \right]=\left[ n\times m \right] \\ & m=n \end{align}\]
Вот и получается, что все три матрицы — $A$, ${{A}^{-1}}$ и $E$ — являются квадратными размером $\left[ n\times n \right]$. Лемма доказана.
Что ж, уже неплохо. Мы видим, что обратимыми бывают лишь квадратные матрицы. Теперь давайте убедимся, что обратная матрица всегда одна.
Лемма 2. Дана матрица $A$ и обратная ей ${{A}^{-1}}$. Тогда эта обратная матрица — единственная.
Доказательство. Пойдём от противного: пусть у матрицы $A$ есть хотя бы два экземпляра обратных —$B$ и $C$. Тогда, согласно определению, верны следующие равенства:
\[\begin{align} & A\cdot B=B\cdot A=E; \\ & A\cdot C=C\cdot A=E. \\ \end{align}\]
Из леммы 1 мы заключаем, что все четыре матрицы — $A$, $B$, $C$ и $E$ — являются квадратными одинакового порядка: $\left[ n\times n \right]$. Следовательно, определено произведение:
\[B\cdot A\cdot C\]
Поскольку умножение матриц ассоциативно (но не коммутативно!), мы можем записать:
\[\begin{align} & B\cdot A\cdot C=\left( B\cdot A \right)\cdot C=E\cdot C=C; \\ & B\cdot A\cdot C=B\cdot \left( A\cdot C \right)=B\cdot E=B; \\ & B\cdot A\cdot C=C=B\Rightarrow B=C. \\ \end{align}\]
Получили единственно возможный вариант: два экземпляра обратной матрицы равны. Лемма доказана.
Приведённые рассуждения почти дословно повторяют доказательство единственность обратного элемента для всех действительных чисел $b\ne 0$. Единственное существенное дополнение — учёт размерности матриц.
Впрочем, мы до сих пор ничего не знаем о том, всякая ли квадратная матрица является обратимой. Тут нам на помощь приходит определитель — это ключевая характеристика для всех квадратных матриц.
Лемма 3. Дана матрица $A$. Если обратная к ней матрица ${{A}^{-1}}$ существует, то определитель исходной матрицы отличен от нуля:
\[\left| A \right|\ne 0\]
Доказательство. Мы уже знаем, что $A$ и ${{A}^{-1}}$ — квадратные матрицы размера $\left[ n\times n \right]$. Следовательно, для каждой из них можно вычислить определитель: $\left| A \right|$ и $\left| {{A}^{-1}} \right|$. Однако определитель произведения равен произведению определителей:
\[\left| A\cdot B \right|=\left| A \right|\cdot \left| B \right|\Rightarrow \left| A\cdot {{A}^{-1}} \right|=\left| A \right|\cdot \left| {{A}^{-1}} \right|\]
Но согласно определению $A\cdot {{A}^{-1}}=E$, а определитель $E$ всегда равен 1, поэтому
\[\begin{align} & A\cdot {{A}^{-1}}=E; \\ & \left| A\cdot {{A}^{-1}} \right|=\left| E \right|; \\ & \left| A \right|\cdot \left| {{A}^{-1}} \right|=1. \\ \end{align}\]
Произведение двух чисел равно единице только в том случае, когда каждое из этих чисел отлично от нуля:
\[\left| A \right|\ne 0;\quad \left| {{A}^{-1}} \right|\ne 0.\]
Вот и получается, что $\left| A \right|\ne 0$. Лемма доказана.
На самом деле это требование вполне логично. Сейчас мы разберём алгоритм нахождения обратной матрицы — и станет совершенно ясно, почему при нулевом определителе никакой обратной матрицы в принципе не может существовать.
Но для начала сформулируем «вспомогательное» определение:
Определение. Вырожденная матрица — это квадратная матрица размера $\left[ n\times n \right]$, чей определитель равен нулю.
Таким образом, мы можем утверждать, что всякая обратимая матрица является невырожденной.
Как найти обратную матрицу
Сейчас мы рассмотрим универсальный алгоритм нахождения обратных матриц. Вообще, существует два общепринятых алгоритма, и второй мы тоже сегодня рассмотрим.
Тот, который будет рассмотрен сейчас, очень эффективен для матриц размера $\left[ 2\times 2 \right]$ и — частично — размера $\left[ 3\times 3 \right]$. А вот начиная с размера $\left[ 4\times 4 \right]$ его лучше не применять. Почему — сейчас сами всё поймёте.
Алгебраические дополнения
Готовьтесь. Сейчас будет боль. Нет, не переживайте: к вам не идёт красивая медсестра в юбке, чулках с кружевами и не сделает укол в ягодицу. Всё куда прозаичнее: к вам идут алгебраические дополнения и Её Величество «Союзная Матрица».
Начнём с главного. Пусть имеется квадратная матрица размера $A=\left[ n\times n \right]$, элементы которой именуются ${{a}_{ij}}$. Тогда для каждого такого элемента можно определить алгебраическое дополнение:
Определение. Алгебраическое дополнение ${{A}_{ij}}$ к элементу ${{a}_{ij}}$, стоящего в $i$-й строке и $j$-м столбце матрицы $A=\left[ n\times n \right]$ — это конструкция вида
\[{{A}_{ij}}={{\left( -1 \right)}^{i+j}}\cdot M_{ij}^{*}\]
Где $M_{ij}^{*}$ — определитель матрицы, полученной из исходной $A$ вычёркиванием той самой $i$-й строки и $j$-го столбца.
Ещё раз. Алгебраическое дополнение к элементу матрицы с координатами $\left( i;j \right)$ обозначается как ${{A}_{ij}}$ и считается по схеме:
- Сначала вычёркиваем из исходной матрицы $i$-строчку и $j$-й столбец. Получим новую квадратную матрицу, и её определитель мы обозначаем как $M_{ij}^{*}$.
- Затем умножаем этот определитель на ${{\left( -1 \right)}^{i+j}}$ — поначалу это выражение может показаться мозговыносящим, но по сути мы просто выясняем знак перед $M_{ij}^{*}$.
- Считаем — получаем конкретное число. Т.е. алгебраическое дополнение — это именно число, а не какая-то новая матрица и т.д.
Саму матрицу $M_{ij}^{*}$ называют дополнительным минором к элементу ${{a}_{ij}}$. И в этом смысле приведённое выше определение алгебраического дополнения является частным случаем более сложного определения — того, что мы рассматривали в уроке про определитель.
Важное замечание. Вообще-то во «взрослой» математике алгебраические дополнения определяются так:
- Берём в квадратной матрице $k$ строчек и $k$ столбцов. На их пересечении получится матрица размера $\left[ k\times k \right]$ — её определитель называется минором порядка $k$ и обозначается ${{M}_{k}}$.
- Затем вычёркиваем эти «избранные» $k$ строчек и $k$ столбцов. Снова получится квадратная матрица — её определитель называется дополнительным минором и обозначается $M_{k}^{*}$.
- Умножаем $M_{k}^{*}$ на ${{\left( -1 \right)}^{t}}$, где $t$ — это (вот сейчас внимание!) сумма номеров всех выбранных строчек и столбцов. Это и будет алгебраическое дополнение.
Взгляните на третий шаг: там вообще-то сумма $2k$ слагаемых! Другое дело, что для $k=1$ мы получим лишь 2 слагаемых — это и будут те самые $i+j$ — «координаты» элемента ${{a}_{ij}}$, для которого мы ищем алгебраическое дополнение.
Таким образом сегодня мы используем слегка упрощённое определение. Но как мы увидим в дальнейшем, его окажется более чем достаточно. Куда важнее следующая штука:
Определение. Союзная матрица $S$ к квадратной матрице $A=\left[ n\times n \right]$ — это новая матрица размера $\left[ n\times n \right]$, которая получается из $A$ заменой ${{a}_{ij}}$ алгебраическими дополнениями ${{A}_{ij}}$:
\[A=\left[ \begin{matrix} {{a}_{11}} & {{a}_{12}} & … & {{a}_{1n}} \\ {{a}_{21}} & {{a}_{22}} & … & {{a}_{2n}} \\ … & … & … & … \\ {{a}_{n1}} & {{a}_{n2}} & … & {{a}_{nn}} \\\end{matrix} \right]\Rightarrow S=\left[ \begin{matrix} {{A}_{11}} & {{A}_{12}} & … & {{A}_{1n}} \\ {{A}_{21}} & {{A}_{22}} & … & {{A}_{2n}} \\ … & … & … & … \\ {{A}_{n1}} & {{A}_{n2}} & … & {{A}_{nn}} \\\end{matrix} \right]\]
Первая мысль, возникающая в момент осознания этого определения — «это сколько же придётся всего считать!» Расслабьтесь: считать придётся, но не так уж и много.:)
Что ж, всё это очень мило, но зачем это нужно? А вот зачем.
Основная теорема
Вернёмся немного назад. Помните, в Лемме 3 утверждалось, что обратимая матрица $A$ всегда не вырождена (т.е. её определитель отличен от нуля: $\left| A \right|\ne 0$).
Так вот, верно и обратное: если матрица $A$ не вырождена, то она всегда обратима. И даже существует схема поиска ${{A}^{-1}}$. Зацените:
Теорема об обратной матрице. Пусть дана квадратная матрица $A=\left[ n\times n \right]$, причём её определитель отличен от нуля: $\left| A \right|\ne 0$. Тогда обратная матрица ${{A}^{-1}}$ существует и считается по формуле:
\[{{A}^{-1}}=\frac{1}{\left| A \right|}\cdot {{S}^{T}}\]
А теперь — всё то же самое, но разборчивым почерком. Чтобы найти обратную матрицу, нужно:
- Посчитать определитель $\left| A \right|$ и убедиться, что он отличен от нуля.
- Составить союзную матрицу $S$, т.е. посчитать 100500 алгебраических дополнений ${{A}_{ij}}$ и расставить их на месте ${{a}_{ij}}$.
- Транспонировать эту матрицу $S$, а затем умножить её на некое число $q={1}/{\left| A \right|}\;$.
И всё! Обратная матрица ${{A}^{-1}}$ найдена. Давайте посмотрим на примеры:
Задача. Найдите обратную матрицу:
\[\left[ \begin{matrix} 3 & 1 \\ 5 & 2 \\\end{matrix} \right]\]
Решение. Проверим обратимость. Посчитаем определитель:
\[\left| A \right|=\left| \begin{matrix} 3 & 1 \\ 5 & 2 \\\end{matrix} \right|=3\cdot 2-1\cdot 5=6-5=1\]
Определитель отличен от нуля. Значит, матрица обратима. Составим союзную матрицу:
\[S=\left[ \begin{matrix} {{A}_{11}} & {{A}_{12}} \\ {{A}_{21}} & {{A}_{22}} \\\end{matrix} \right]\]
Посчитаем алгебраические дополнения:
\[\begin{align} & {{A}_{11}}={{\left( -1 \right)}^{1+1}}\cdot \left| 2 \right|=2; \\ & {{A}_{12}}={{\left( -1 \right)}^{1+2}}\cdot \left| 5 \right|=-5; \\ & {{A}_{21}}={{\left( -1 \right)}^{2+1}}\cdot \left| 1 \right|=-1; \\ & {{A}_{22}}={{\left( -1 \right)}^{2+2}}\cdot \left| 3 \right|=3. \\ \end{align}\]
Обратите внимание: определители |2|, |5|, |1| и |3| — это именно определители матриц размера $\left[ 1\times 1 \right]$, а не модули. Т.е. если в определителях стояли отрицательные числа, убирать «минус» не надо.
Итого наша союзная матрица выглядит так:
\[S=\left[ \begin{array}{*{35}{r}}2 & -5 \\ -1 & 3 \\\end{array} \right]\]
Осталось посчитать обратную:
\[{{A}^{-1}}=\frac{1}{\left| A \right|}\cdot {{S}^{T}}=\frac{1}{1}\cdot {{\left[ \begin{array}{*{35}{r}} 2 & -5 \\ -1 & 3 \\\end{array} \right]}^{T}}=\left[ \begin{array}{*{35}{r}} 2 & -1 \\ -5 & 3 \\\end{array} \right]\]
Ну вот и всё. Задача решена.
Ответ. $\left[ \begin{array}{*{35}{r}} 2 & -1 \\ -5 & 3 \\\end{array} \right]$
Задача. Найдите обратную матрицу:
\[\left[ \begin{array}{*{35}{r}} 1 & -1 & 2 \\ 0 & 2 & -1 \\ 1 & 0 & 1 \\\end{array} \right]\]
Решение. Опять считаем определитель:
\[\begin{align} & \left| \begin{array}{*{35}{r}} 1 & -1 & 2 \\ 0 & 2 & -1 \\ 1 & 0 & 1 \\\end{array} \right|=\begin{matrix} \left( 1\cdot 2\cdot 1+\left( -1 \right)\cdot \left( -1 \right)\cdot 1+2\cdot 0\cdot 0 \right)- \\ -\left( 2\cdot 2\cdot 1+\left( -1 \right)\cdot 0\cdot 1+1\cdot \left( -1 \right)\cdot 0 \right) \\\end{matrix}= \\ & =\left( 2+1+0 \right)-\left( 4+0+0 \right)=-1\ne 0. \\ \end{align}\]
Определитель отличен от нуля — матрица обратима. А вот сейчас будет самая жесть: надо посчитать аж 9 (девять, мать их!) алгебраических дополнений. И каждое из них будет содержать определитель $\left[ 2\times 2 \right]$. Полетели:
\[\begin{matrix} {{A}_{11}}={{\left( -1 \right)}^{1+1}}\cdot \left| \begin{matrix} 2 & -1 \\ 0 & 1 \\\end{matrix} \right|=2; \\ {{A}_{12}}={{\left( -1 \right)}^{1+2}}\cdot \left| \begin{matrix} 0 & -1 \\ 1 & 1 \\\end{matrix} \right|=-1; \\ {{A}_{13}}={{\left( -1 \right)}^{1+3}}\cdot \left| \begin{matrix} 0 & 2 \\ 1 & 0 \\\end{matrix} \right|=-2; \\ … \\ {{A}_{33}}={{\left( -1 \right)}^{3+3}}\cdot \left| \begin{matrix} 1 & -1 \\ 0 & 2 \\\end{matrix} \right|=2; \\ \end{matrix}\]
Короче, союзная матрица будет выглядеть так:
\[S=\left[ \begin{matrix} 2 & -1 & -2 \\ 1 & -1 & -1 \\ -3 & 1 & 2 \\\end{matrix} \right]\]
Следовательно, обратная матрица будет такой:
\[{{A}^{-1}}=\frac{1}{-1}\cdot \left[ \begin{matrix} 2 & -1 & -2 \\ 1 & -1 & -1 \\ -3 & 1 & 2 \\\end{matrix} \right]=\left[ \begin{array}{*{35}{r}}-2 & -1 & 3 \\ 1 & 1 & -1 \\ 2 & 1 & -2 \\\end{array} \right]\]
Ну и всё. Вот и ответ.
Ответ. $\left[ \begin{array}{*{35}{r}} -2 & -1 & 3 \\ 1 & 1 & -1 \\ 2 & 1 & -2 \\\end{array} \right]$
Как видите, в конце каждого примера мы выполняли проверку. В связи с этим важное замечание:
Не ленитесь выполнять проверку. Умножьте исходную матрицу на найденную обратную — должна получиться $E$.
Выполнить эту проверку намного проще и быстрее, чем искать ошибку в дальнейших вычислениях, когда, например, вы решаете матричное уравнение.
Альтернативный способ
Как я и говорил, теорема об обратной матрице прекрасно работает для размеров $\left[ 2\times 2 \right]$ и $\left[ 3\times 3 \right]$ (в последнем случае — уже не так уж и «прекрасно»), а вот для матриц больших размеров начинается прям печаль.
Но не переживайте: есть альтернативный алгоритм, с помощью которого можно невозмутимо найти обратную хоть для матрицы $\left[ 10\times 10 \right]$. Но, как это часто бывает, для рассмотрения этого алгоритма нам потребуется небольшая теоретическая вводная.
Элементарные преобразования
Среди всевозможных преобразований матрицы есть несколько особых — их называют элементарными. Таких преобразований ровно три:
- Умножение. Можно взять $i$-ю строку (столбец) и умножить её на любое число $k\ne 0$;
- Сложение. Прибавить к $i$-й строке (столбцу) любую другую $j$-ю строку (столбец), умноженную на любое число $k\ne 0$ (можно, конечно, и $k=0$, но какой в этом смысл? Ничего не изменится же).
- Перестановка. Взять $i$-ю и $j$-ю строки (столбцы) и поменять местами.
Почему эти преобразования называются элементарными (для больших матриц они выглядят не такими уж элементарными) и почему их только три — эти вопросы выходят за рамки сегодняшнего урока. Поэтому не будем вдаваться в подробности.
Важно другое: все эти извращения нам предстоит выполнять над присоединённой матрицей. Да, да: вы не ослышались. Сейчас будет ещё одно определение — последнее в сегодняшнем уроке.
Присоединённая матрица
Наверняка в школе вы решали системы уравнений методом сложения. Ну, там, вычесть из одной строки другую, умножить какую-то строку на число — вот это вот всё.
Так вот: сейчас будет всё то же, но уже «по-взрослому». Готовы?
Определение. Пусть дана матрица $A=\left[ n\times n \right]$ и единичная матрица $E$ такого же размера $n$. Тогда присоединённая матрица $\left[ A\left| E \right. \right]$ — это новая матрица размера $\left[ n\times 2n \right]$, которая выглядит так:
\[\left[ A\left| E \right. \right]=\left[ \begin{array}{rrrr|rrrr}{{a}_{11}} & {{a}_{12}} & … & {{a}_{1n}} & 1 & 0 & … & 0 \\{{a}_{21}} & {{a}_{22}} & … & {{a}_{2n}} & 0 & 1 & … & 0 \\… & … & … & … & … & … & … & … \\{{a}_{n1}} & {{a}_{n2}} & … & {{a}_{nn}} & 0 & 0 & … & 1 \\\end{array} \right]\]
Короче говоря, берём матрицу $A$, справа приписываем к ней единичную матрицу $E$ нужного размера, разделяем их вертикальной чертой для красоты — вот вам и присоединённая.:)
В чём прикол? А вот в чём:
Теорема. Пусть матрица $A$ обратима. Рассмотрим присоединённую матрицу $\left[ A\left| E \right. \right]$. Если с помощью элементарных преобразований строк привести её к виду $\left[ E\left| B \right. \right]$, т.е. путём умножения, вычитания и перестановки строк получить из $A$ матрицу $E$ справа, то полученная слева матрица $B$ — это обратная к $A$:
\[\left[ A\left| E \right. \right]\to \left[ E\left| B \right. \right]\Rightarrow B={{A}^{-1}}\]
Вот так всё просто! Короче говоря, алгоритм нахождения обратной матрицы выглядит так:
- Записать присоединённую матрицу $\left[ A\left| E \right. \right]$;
- Выполнять элементарные преобразования строк до тех пор, пока права вместо $A$ не появится $E$;
- Разумеется, слева тоже что-то появится — некая матрица $B$. Это и будет обратная;
- PROFIT!:)
Конечно, сказать намного проще, чем сделать. Поэтому давайте рассмотрим парочку примеров: для размеров $\left[ 3\times 3 \right]$ и $\left[ 4\times 4 \right]$.
Задача. Найдите обратную матрицу:
\[\left[ \begin{array}{*{35}{r}} 1 & 5 & 1 \\ 3 & 2 & 1 \\ 6 & -2 & 1 \\\end{array} \right]\]
Решение. Составляем присоединённую матрицу:
\[\left[ \begin{array}{rrr|rrr} 1 & 5 & 1 & 1 & 0 & 0 \\ 3 & 2 & 1 & 0 & 1 & 0 \\ 6 & -2 & 1 & 0 & 0 & 1 \\\end{array} \right]\]
Поскольку последний столбец исходной матрицы заполнен единицами, вычтем первую строку из остальных:
\[\begin{align} & \left[ \begin{array}{rrr|rrr} 1 & 5 & 1 & 1 & 0 & 0 \\ 3 & 2 & 1 & 0 & 1 & 0 \\ 6 & -2 & 1 & 0 & 0 & 1 \\\end{array} \right]\begin{matrix} \downarrow \\ -1 \\ -1 \\\end{matrix}\to \\ & \to \left[ \begin{array}{rrr|rrr} 1 & 5 & 1 & 1 & 0 & 0 \\ 2 & -3 & 0 & -1 & 1 & 0 \\ 5 & -7 & 0 & -1 & 0 & 1 \\\end{array} \right] \\ \end{align}\]
Больше единиц нет, кроме первой строки. Но её мы не трогаем, иначе в третьем столбце начнут «размножаться» только что убранные единицы.
Зато можем вычесть вторую строку дважды из последней — получим единицу в левом нижнем углу:
\[\begin{align} & \left[ \begin{array}{rrr|rrr} 1 & 5 & 1 & 1 & 0 & 0 \\ 2 & -3 & 0 & -1 & 1 & 0 \\ 5 & -7 & 0 & -1 & 0 & 1 \\\end{array} \right]\begin{matrix} \ \\ \downarrow \\ -2 \\\end{matrix}\to \\ & \left[ \begin{array}{rrr|rrr} 1 & 5 & 1 & 1 & 0 & 0 \\ 2 & -3 & 0 & -1 & 1 & 0 \\ 1 & -1 & 0 & 1 & -2 & 1 \\\end{array} \right] \\ \end{align}\]
Теперь можно вычесть последнюю строку из первой и дважды из второй — таким образом мы «занулим» первый столбец:
\[\begin{align} & \left[ \begin{array}{rrr|rrr} 1 & 5 & 1 & 1 & 0 & 0 \\ 2 & -3 & 0 & -1 & 1 & 0 \\ 1 & -1 & 0 & 1 & -2 & 1 \\\end{array} \right]\begin{matrix} -1 \\ -2 \\ \uparrow \\\end{matrix}\to \\ & \to \left[ \begin{array}{rrr|rrr} 0 & 6 & 1 & 0 & 2 & -1 \\ 0 & -1 & 0 & -3 & 5 & -2 \\ 1 & -1 & 0 & 1 & -2 & 1 \\\end{array} \right] \\ \end{align}\]
Умножим вторую строку на −1, а затем вычтем её 6 раз из первой и прибавим 1 раз к последней:
\[\begin{align} & \left[ \begin{array}{rrr|rrr} 0 & 6 & 1 & 0 & 2 & -1 \\ 0 & -1 & 0 & -3 & 5 & -2 \\ 1 & -1 & 0 & 1 & -2 & 1 \\\end{array} \right]\begin{matrix} \ \\ \left| \cdot \left( -1 \right) \right. \\ \ \\\end{matrix}\to \\ & \to \left[ \begin{array}{rrr|rrr} 0 & 6 & 1 & 0 & 2 & -1 \\ 0 & 1 & 0 & 3 & -5 & 2 \\ 1 & -1 & 0 & 1 & -2 & 1 \\\end{array} \right]\begin{matrix} -6 \\ \updownarrow \\ +1 \\\end{matrix}\to \\ & \to \left[ \begin{array}{rrr|rrr} 0 & 0 & 1 & -18 & 32 & -13 \\ 0 & 1 & 0 & 3 & -5 & 2 \\ 1 & 0 & 0 & 4 & -7 & 3 \\\end{array} \right] \\ \end{align}\]
Осталось лишь поменять местами строки 1 и 3:
\[\left[ \begin{array}{rrr|rrr} 1 & 0 & 0 & 4 & -7 & 3 \\ 0 & 1 & 0 & 3 & -5 & 2 \\ 0 & 0 & 1 & -18 & 32 & -13 \\\end{array} \right]\]
Готово! Справа — искомая обратная матрица.
Ответ. $\left[ \begin{array}{*{35}{r}}4 & -7 & 3 \\ 3 & -5 & 2 \\ -18 & 32 & -13 \\\end{array} \right]$
Задача. Найдите обратную матрицу:
\[\left[ \begin{matrix} 1 & 4 & 2 & 3 \\ 1 & -2 & 1 & -2 \\ 1 & -1 & 1 & 1 \\ 0 & -10 & -2 & -5 \\\end{matrix} \right]\]
Решение. Снова составляем присоединённую:
\[\left[ \begin{array}{rrrr|rrrr} 1 & 4 & 2 & 3 & 1 & 0 & 0 & 0 \\ 1 & -2 & 1 & -2 & 0 & 1 & 0 & 0 \\ 1 & -1 & 1 & 1 & 0 & 0 & 1 & 0 \\ 0 & -10 & -2 & -5 & 0 & 0 & 0 & 1 \\\end{array} \right]\]
Немного позалимаем, попечалимся от того, сколько сейчас придётся считать… и начнём считать. Для начала «обнулим» первый столбец, вычитая строку 1 из строк 2 и 3:
\[\begin{align} & \left[ \begin{array}{rrrr|rrrr} 1 & 4 & 2 & 3 & 1 & 0 & 0 & 0 \\ 1 & -2 & 1 & -2 & 0 & 1 & 0 & 0 \\ 1 & -1 & 1 & 1 & 0 & 0 & 1 & 0 \\ 0 & -10 & -2 & -5 & 0 & 0 & 0 & 1 \\\end{array} \right]\begin{matrix} \downarrow \\ -1 \\ -1 \\ \ \\\end{matrix}\to \\ & \to \left[ \begin{array}{rrrr|rrrr} 1 & 4 & 2 & 3 & 1 & 0 & 0 & 0 \\ 0 & -6 & -1 & -5 & -1 & 1 & 0 & 0 \\ 0 & -5 & -1 & -2 & -1 & 0 & 1 & 0 \\ 0 & -10 & -2 & -5 & 0 & 0 & 0 & 1 \\\end{array} \right] \\ \end{align}\]
Наблюдаем слишком много «минусов» в строках 2—4. Умножим все три строки на −1, а затем «выжжем» третий столбец, вычитая строку 3 из остальных:
\[\begin{align} & \left[ \begin{array}{rrrr|rrrr} 1 & 4 & 2 & 3 & 1 & 0 & 0 & 0 \\ 0 & -6 & -1 & -5 & -1 & 1 & 0 & 0 \\ 0 & -5 & -1 & -2 & -1 & 0 & 1 & 0 \\ 0 & -10 & -2 & -5 & 0 & 0 & 0 & 1 \\\end{array} \right]\begin{matrix} \ \\ \left| \cdot \left( -1 \right) \right. \\ \left| \cdot \left( -1 \right) \right. \\ \left| \cdot \left( -1 \right) \right. \\\end{matrix}\to \\ & \to \left[ \begin{array}{rrrr|rrrr} 1 & 4 & 2 & 3 & 1 & 0 & 0 & 0 \\ 0 & 6 & 1 & 5 & 1 & -1 & 0 & 0 \\ 0 & 5 & 1 & 2 & 1 & 0 & -1 & 0 \\ 0 & 10 & 2 & 5 & 0 & 0 & 0 & -1 \\\end{array} \right]\begin{matrix} -2 \\ -1 \\ \updownarrow \\ -2 \\\end{matrix}\to \\ & \to \left[ \begin{array}{rrrr|rrrr} 1 & -6 & 0 & -1 & -1 & 0 & 2 & 0 \\ 0 & 1 & 0 & 3 & 0 & -1 & 1 & 0 \\ 0 & 5 & 1 & 2 & 1 & 0 & -1 & 0 \\ 0 & 0 & 0 & 1 & -2 & 0 & 2 & -1 \\\end{array} \right] \\ \end{align}\]
Теперь самое время «поджарить» последний столбец исходной матрицы: вычитаем строку 4 из остальных:
\[\begin{align} & \left[ \begin{array}{rrrr|rrrr} 1 & -6 & 0 & -1 & -1 & 0 & 2 & 0 \\ 0 & 1 & 0 & 3 & 0 & -1 & 1 & 0 \\ 0 & 5 & 1 & 2 & 1 & 0 & -1 & 0 \\ 0 & 0 & 0 & 1 & -2 & 0 & 2 & -1 \\\end{array} \right]\begin{matrix} +1 \\ -3 \\ -2 \\ \uparrow \\\end{matrix}\to \\ & \to \left[ \begin{array}{rrrr|rrrr} 1 & -6 & 0 & 0 & -3 & 0 & 4 & -1 \\ 0 & 1 & 0 & 0 & 6 & -1 & -5 & 3 \\ 0 & 5 & 1 & 0 & 5 & 0 & -5 & 2 \\ 0 & 0 & 0 & 1 & -2 & 0 & 2 & -1 \\\end{array} \right] \\ \end{align}\]
Финальный бросок: «выжигаем» второй столбец, вычитая строку 2 из строки 1 и 3:
\[\begin{align} & \left[ \begin{array}{rrrr|rrrr} 1 & -6 & 0 & 0 & -3 & 0 & 4 & -1 \\ 0 & 1 & 0 & 0 & 6 & -1 & -5 & 3 \\ 0 & 5 & 1 & 0 & 5 & 0 & -5 & 2 \\ 0 & 0 & 0 & 1 & -2 & 0 & 2 & -1 \\\end{array} \right]\begin{matrix} 6 \\ \updownarrow \\ -5 \\ \ \\\end{matrix}\to \\ & \to \left[ \begin{array}{rrrr|rrrr} 1 & 0 & 0 & 0 & 33 & -6 & -26 & -17 \\ 0 & 1 & 0 & 0 & 6 & -1 & -5 & 3 \\ 0 & 0 & 1 & 0 & -25 & 5 & 20 & -13 \\ 0 & 0 & 0 & 1 & -2 & 0 & 2 & -1 \\\end{array} \right] \\ \end{align}\]
И снова слева единичная матрица, значит справа — обратная.:)
Ответ. $\left[ \begin{matrix} 33 & -6 & -26 & 17 \\ 6 & -1 & -5 & 3 \\ -25 & 5 & 20 & -13 \\ -2 & 0 & 2 & -1 \\\end{matrix} \right]$
Ну вот и всё. Проверку сделайте сами — мне в лом.:)
Смотрите также:
- Определитель
- Угол между двумя прямыми
- Тест к уроку «Сложные выражения с дробями» (легкий)
- Правила комбинаторики в задаче B6
- Пример решения задачи 15
- Тригонометрия в задаче B15: решаем без производных