В линейной алгебре , п матрицы с размерностью п квадратной матрицы называется обратимым (также неособым или невырожденным ), если существует п матрица с размерностью п квадратной матрица B таким образом, что
где я п обозначает в н матрице с размерностью п единичной матрицы , а умножение используется обычное умножение матриц . Если это так, то матрица B однозначно определяется A , и называется (мультипликативный) обратный из А , обозначаемый A -1 . [1] [2] Матрица инверсии является процесс нахождения матрицы B , которая удовлетворяет уравнению перед для данной обратимой матрицы А .
Квадратная матрица, не обратимы называется сингулярным или вырожденным . Квадратная матрица сингулярна тогда и только тогда, когда ее определитель равен нулю. [3] Сингулярные матрицы встречаются редко в том смысле, что если элементы квадратной матрицы выбираются случайным образом из любой конечной области числовой прямой или комплексной плоскости, вероятность того, что матрица является сингулярной, равна 0, то есть «почти никогда» не будет. быть единичным. Non-квадратные матрицы ( т матрицу с размерностью п матриц , для которых м ≠ п ) не имеет обратного. Однако в некоторых случаях такая матрица может иметь левую обратную или правую обратную . Если является т матрицу с размерностью п и ранга из А равна п ( п ≤ м ), то имеет левый обратный, An в н матрицу с размерностью м матрицы B таким образом, что BA = I н . Если A имеет ранг m ( m ≤ n ), то у него есть правая обратная матрица B размером n на m, такая что AB = I m .
Хотя наиболее распространенным случаем является матрица над действительными или комплексными числами, все эти определения могут быть даны для матриц над любым кольцом . Однако в случае коммутативности кольца условием обратимости квадратной матрицы является то, что ее определитель обратим в кольце, что в общем случае является более строгим требованием, чем ненулевое значение. Для некоммутативного кольца обычный определитель не определен. Условия существования левообратных или правообратных более сложны, так как понятие ранга не существует над кольцами.
Набор обратимых матриц размера n × n вместе с операцией умножения матриц (и элементами из кольца R ) образуют группу , общую линейную группу степени n , обозначаемую GL n ( R ) . [1]
Характеристики
Теорема об обратимой матрице
Пусть A будет квадратной матрицей n на n над полем K (например, полем R действительных чисел). Следующие утверждения эквивалентны (т. Е. Либо все они истинны, либо полностью ложны для любой данной матрицы): [4]
- A обратимо, то есть A имеет обратный, невырожденный или невырожденный.
- Является строкой-эквивалентно к N матрицы с размерностью п единичной матрицы I н .
- Является колонным эквивалентны к N матрицы с размерностью п единичной матрицы I н .
- A имеет n опорных позиций .
- det A ≠ 0 . В общем случае квадратная матрица над коммутативным кольцом обратима тогда и только тогда, когда ее определитель является единицей в этом кольце.
- А имеет полный ранг; то есть ранг A = n .
- Уравнение Ax = 0 имеет только тривиальное решение x = 0 .
- Ядро из А тривиально, то есть, он содержит только вектор нуль в качестве элемента, кег ( ) = { 0 }.
- Уравнение Ax = b имеет ровно одно решение для каждого b в K n .
- Столбцы А являются линейно независимыми .
- Столбцы A образуют K n .
- Col A = K n .
- Столбцы А образуют базис в K н .
- Отображение линейного преобразования x в Ax является биекцией из K n в K n .
- Существует матрица B размером n x n такая, что AB = I n = BA .
- Транспонированная Т является обратимой матрицей (следовательно , ряды А являются линейно независимыми , пролет К п , и образуют основу из К п ).
- Число 0 не является собственным значением из A .
- Матрица A может быть выражена как конечное произведение элементарных матриц .
- Матрица A имеет левую обратную (то есть существует B такая, что BA = I ) или правую обратную (то есть существует C такая, что AC = I ), и в этом случае существуют как левая, так и правая инверсия, и В = С = А -1 .
Прочие свойства
Кроме того, для обратимой матрицы A выполняются следующие свойства :
- ( А -1 ) -1 = А ;
- ( k A ) −1 = k −1 A −1 для ненулевого скаляра k ;
- ( Ax ) + = x + A −1, если A имеет ортонормированные столбцы, где + обозначает обратное преобразование Мура – Пенроуза, а x - вектор;
- ( A T ) −1 = ( A −1 ) T ;
- Для любого обратимого N матрицы с размерностью п матриц A и B , ( AB ) -1 = B -1 A -1 . В более общем смысле , если А 1 , ..., к обратимым п матрица с размерностью N матрицы, то ( 1 2 ⋯ к -1 к ) -1 =−1
кА−1
к −1⋯ А−1
2А−1
1; - det A −1 = (det A ) −1 .
Строки обратной матрицы V матриц U являются ортонормированными к колонкам U (и наоборот перестановка строк для столбцов). Чтобы увидеть это, предположим, что UV = VU = I, где строки V обозначены кака столбцы U как для . Тогда ясно, что евклидов скалярное произведение любых двух. Это свойство также может быть полезно при построении обратной квадратной матрицы в некоторых случаях, когда известен набор ортогональных векторов (но не обязательно ортонормированных векторов) столбцам U. В этом случае, можно применить итерационный процесс Грама-Шмидт к этому исходному набору для определения строк обратной V .
Матрица, которая является своей собственной обратной (т. Е. Матрица A такая, что A = A −1 и A 2 = I ), называется инволютивной матрицей .
По отношению к его адъюгату
Adjugate матрицы можно использовать, чтобы найти обратное следующим образом:
Если обратимая матрица, то
По отношению к единичной матрице
Из ассоциативности умножения матриц следует, что если
для конечных квадратных матриц A и B , то также
- [5]
Плотность
Над поля действительных чисел, множества особых п матрицы с размерностью п матриц, рассматриваемым как подмножество R п × п , является множество нуля , то есть, имеет лебегову меру нуля . Это верно, потому что особые матрицы являются корнями детерминантной функции. Это непрерывная функция, поскольку она является полиномом от элементов матрицы. Таким образом , на языке теории меры , почти все п матрицы с размерностью п матрицы обратимы.
Кроме того, п матрицы с размерностью п обратимых матрицами являются плотным открытым множеством в топологическом пространстве из всех N матрицы с размерностью п матриц. Эквивалентно, множество особых матриц закрыто , и нигде не плотно в пространстве п матрицы с размерностью п матриц.
Однако на практике можно встретить необратимые матрицы. И при численных расчетах матрицы, которые являются обратимыми, но близкими к необратимым, все еще могут быть проблематичными; такие матрицы называются плохо обусловленными .
Примеры
Рассмотрим следующую матрицу 2 на 2:
Матрица обратимо. Чтобы проверить это, можно вычислить, что, которая не равна нулю.
В качестве примера необратимой или сингулярной матрицы рассмотрим матрицу
Определитель равен 0, что является необходимым и достаточным условием необратимости матрицы.
Методы обращения матриц
Гауссово исключение
Исключение Гаусса – Жордана - это алгоритм, который можно использовать для определения, является ли данная матрица обратимой, и для нахождения обратной. Альтернативой является разложение LU , которое генерирует верхние и нижние треугольные матрицы, которые легче инвертировать.
Метод Ньютона
Обобщение метода Ньютона, используемого для мультипликативного обратного алгоритма, может быть удобным, если удобно найти подходящее начальное начальное число:
Виктор Пэн и Джон Рейф проделали работу, которая включает способы создания начального семени. [6] [7] Журнал Byte резюмировал один из их подходов. [8]
Метод Ньютона особенно полезен при работе с семействами связанных матриц, которые ведут себя достаточно похоже на последовательность, созданную для гомотопии выше: иногда хорошей отправной точкой для уточнения приближения для новой обратной матрицы может быть уже полученная обратная матрица предыдущей, которая почти соответствует текущая матрица, например, пара последовательностей обратных матриц, используемая для получения квадратного корня матрицы с помощью итерации Денмана – Биверса ; для этого может потребоваться более одного прохода итерации в каждой новой матрице, если они не достаточно близко друг к другу, чтобы хватило только одного. Метод Ньютона также полезен для «доработки» поправок к алгоритму Гаусса – Джордана, который был загрязнен небольшими ошибками из-за несовершенной компьютерной арифметики .
Метод Кэли – Гамильтона
Теорема Кэли – Гамильтона допускает обратное быть выраженным в терминах , следы и силы : [9]
где это измерение , а также это след матрицыдается суммой главной диагонали. Сумма принимается и наборы всех удовлетворяющий линейному диофантову уравнению
Формулу можно переписать в терминах полных многочленов Белла от аргументов в виде
Собственное разложение
Если матрица может быть eigendecomposed, и если ни один из его собственных значений равны нулю, то не обратим и обратный дается
где квадратная ( N × N ) матрица, i -й столбец которой является собственным вектором из , а также - диагональная матрица , диагональные элементы которой являются соответствующими собственными значениями, т. е.. Если симметрично, гарантированно будет ортогональной матрицей , поэтому. Кроме того, поскольку - диагональная матрица, обратная к ней легко вычисляется:
Разложение Холецкого
Если матрица является положительно определена , то его обратная можно получить как
где L является нижним треугольным Холецким разложением на А , а L * обозначает сопряженные транспонированную L .
Аналитическое решение
Запись транспонированной матрицы кофакторов , известной как сопряженная матрица , также может быть эффективным способом вычисления обратного значения малых матриц, но этот рекурсивный метод неэффективен для больших матриц. Чтобы определить обратное, вычисляем матрицу сомножителей:
чтобы
где | А | является фактором , определяющим из А , С является матрица кофакторов , и С Т представляет собой матрицу транспонирование .
Инверсия матриц 2 × 2
Уравнение кофактора перечисленных выше , дает следующий результат для 2 × 2 матриц. Инверсия этих матриц может быть произведена следующим образом: [10]
Это возможно, потому что 1 / ( ad - bc ) является обратной величиной определителя рассматриваемой матрицы, и та же стратегия может быть использована для других размеров матриц.
Метод Кэли – Гамильтона дает
Инверсия матриц 3 × 3
Вычислительно эффективная инверсия матрицы 3 × 3 дается выражением
(где скаляр A не следует путать с матрицей A ).
Если определитель не равен нулю, матрица обратима, а элементы промежуточной матрицы в правой части выше заданы формулой
Определитель A можно вычислить, применив правило Сарруса следующим образом:
Разложение Кэли – Гамильтона дает
Общее обратное 3 × 3 может быть кратко выражено в терминах перекрестного произведения и тройного произведения . Если матрица (состоящий из трех векторов-столбцов, , , а также ) обратима, обратная ему дается формулой
Определитель A, , равна тройному произведению , , а также - объем параллелепипеда, образованный строками или столбцами:
Правильность формулы можно проверить, используя свойства перекрестного и тройного произведения, а также отметив, что для групп всегда совпадают левые и правые инверсии. Интуитивно из-за перекрестных произведений каждая строка ортогонален двум несоответствующим столбцам матрицы (в результате чего недиагональные члены быть нулевым). Деление на
вызывает диагональные элементы быть единством. Например, первая диагональ:
Инверсия матриц 4 × 4
С увеличением размерности выражения для обратной величины A усложняются. При n = 4 метод Кэли – Гамильтона приводит к выражению, которое все еще остается приемлемым:
Блочная инверсия
Матрицы также можно инвертировать поблочно , используя следующую формулу аналитического обращения:
( 1 )
где A , B , C и D - матричные субблоки произвольного размера. ( A должен быть квадратным, чтобы его можно было инвертировать. Кроме того, A и D - CA −1 B должны быть невырожденными. [11] ) Эта стратегия особенно выгодна, если A диагональна, а D - CA −1 B ( алгоритм Шура дополнение к A ) является маленькой матрицей, так как это единственные матрицы, требующие инверсии.
Этот метод был изобретен несколько раз заново, и он принадлежит Гансу Больцу (1923), [ необходима цитата ], который использовал его для обращения геодезических матриц, и Тадеушу Банахевичу (1937), который обобщил его и доказал его правильность.
Теорема о нулевом значении говорит, что нулевое значение A равно нулю субблока в правом нижнем углу обратной матрицы, а нулевое значение B равно нулю субблока в верхнем правом углу обратной матрицы.
Процедура инверсии, которая привела к уравнению ( 1 ), выполняла операции матричного блока, которые сначала работали с C и D. Вместо этого, если сначала оперировать A и B , и при условии, что D и A - BD −1 C неособые, [12] результат будет
( 2 )
Приравнивание уравнений ( 1 ) и ( 2 ) приводит к
( 3 )
где уравнение ( 3 ) представляет собой матричное тождество Вудбери , что эквивалентно биномиальной обратной теореме .
Если A и D оба обратимы, то две вышеупомянутые инверсии блочной матрицы могут быть объединены, чтобы обеспечить простую факторизацию
( 2 )
Согласно тождеству Вайнштейна – Ароншайна одна из двух матриц в блочно-диагональной матрице обратима в точности тогда, когда другая обратима.
Поскольку для поблочного обращения матрицы размера n × n требуется инверсия двух матриц половинного размера и 6 умножений между двумя матрицами половинного размера, можно показать, что алгоритм разделяй и властвуй, который использует поблочное обращение для инвертирования матрицы, работает с тем же временная сложность как алгоритм умножения матриц, который используется внутри компании. [13] Исследование сложности матричного умножения показывает, что существуют алгоритмы матричного умножения со сложностью O ( n 2.3727 ) операций, в то время как наилучшая доказанная нижняя граница - Ω ( n 2 log n ) . [14]
Эта формула значительно упрощается, когда верхняя правая блочная матрица - нулевая матрица. Эта формулировка полезна, когда матрицы а также имеют относительно простые обратные формулы (или псевдообратные в случае, когда блоки не все квадратные. В этом частном случае формула обращения блочной матрицы, изложенная выше в полной общности, становится
По серии Neumann
Если матрица A обладает свойством, что
тогда A невырожден, и его обратное может быть выражено рядом Неймана : [15]
Усечение суммы приводит к «приблизительному» обратному результату, который может быть полезен в качестве предобуславливателя . Обратите внимание, что усеченный ряд можно ускорить экспоненциально, если учесть, что ряд Неймана представляет собой геометрическую сумму . Таким образом, он удовлетворяет
- .
Поэтому только умножения матриц необходимы для вычисления условия суммы.
В более общем смысле, если A находится «рядом» с обратимой матрицей X в том смысле, что
то A невырожден и обратный ему равен
Если также верно, что A - X имеет ранг 1, то это упрощается до
p -адическое приближение
Если A - матрица с целыми или рациональными коэффициентами, и мы ищем решение в рациональных числах произвольной точности , то метод p -адической аппроксимации сходится к точному решению в, предполагая стандартный используется матричное умножение. [16] Метод основан на решении n линейных систем с помощью метода p -адической аппроксимации Диксона (каждая в) и доступен как таковой в программном обеспечении, специализирующемся на матричных операциях произвольной точности, например, в IML. [17]
Метод взаимных базисных векторов
Учитывая квадратная матрица , , с участием строки интерпретируются как векторов ( Предполагается суммирование Эйнштейна ), гдеявляются стандартной ортонормированный базис в евклидовом пространстве (), затем, используя алгебру Клиффорда (или геометрическую алгебру ), мы вычисляем обратные (иногда называемые двойственными ) векторами-столбцами как столбцы обратной матрицы . Обратите внимание, что место "" указывает, что ""удалено с этого места в приведенном выше выражении для . Тогда у нас есть, где - дельта Кронекера . У нас также есть, как требуется. Если векторы не являются линейно независимыми, то и матрица не обратима (не имеет обратного).
Производная обратной матрицы
Предположим, что обратимая матрица A зависит от параметра t . Тогда производная, обратная к A по t, равна [18]
Чтобы вывести указанное выше выражение для производной обратной матрицы A , можно дифференцировать определение матрицы, обратнойа затем решите для обратного к A :
Вычитание с обеих сторон вышеуказанного и умножая справа на дает правильное выражение для производной обратного:
Аналогично, если это небольшое число, тогда
В более общем смысле, если
тогда,
Учитывая положительное целое число ,
Следовательно,
Обобщенная обратная
Некоторые свойства обратных матриц являются общими для обобщенных обратных матриц (например, обратная матрица Мура – Пенроуза ), которая может быть определена для любой матрицы размера m на n .
Приложения
Для большинства практических приложений нет необходимости инвертировать матрицу для решения системы линейных уравнений ; однако, для однозначного решения, то есть необходимо, чтобы матрица участвует обратимым.
Методы декомпозиции, такие как LU-декомпозиция , намного быстрее, чем инверсия, также были разработаны различные быстрые алгоритмы для специальных классов линейных систем.
Регрессия / наименьшие квадраты
Хотя явное обратное не требуется для оценки вектора неизвестных, это самый простой способ оценить их точность, находящуюся на диагонали обратной матрицы (матрица апостериорной ковариации вектора неизвестных). Однако во многих случаях известны более быстрые алгоритмы вычисления только диагональных элементов обратной матрицы. [19]
Обращение матриц в симуляциях в реальном времени
Инверсия матриц играет важную роль в компьютерной графике , особенно в рендеринге трехмерной графики и трехмерном моделировании. Примеры включают преобразование лучей из экрана в мир , преобразования объектов из мира в подпространство в мир и физические симуляции.
Матрица инверсия в беспроводной связи MIMO
Инверсия матрицы также играет важную роль в технологии MIMO (Multiple-Input, Multiple-Output) в беспроводной связи. Система MIMO состоит из N передающих и M приемных антенн. Уникальные сигналы, занимающие одну и ту же полосу частот, отправляются через N передающих антенн и принимаются через M приемных антенн. Сигнала , поступающего на каждой приемной антенны будет представлять собой линейную комбинацию из N передаваемых сигналов , образующих N × M матрицу передачи H . Крайне важно, чтобы матрица H была обратимой, чтобы приемник мог вычислить передаваемую информацию.
Смотрите также
- Биномиальная обратная теорема
- LU разложение
- Разложение матрицы
- Матричный квадратный корень
- Минор (линейная алгебра)
- Частичная обратная матрица
- Псевдообратный
- Разложение по сингулярным числам
- Тождество матрицы Вудбери
Рекомендации
- ^ a b «Исчерпывающий список символов алгебры» . Математическое хранилище . 2020-03-25 . Проверено 8 сентября 2020 .
- ^ «Обратимые матрицы» . www.sosmath.com . Проверено 8 сентября 2020 .
- ^ Вайсштейн, Эрик В. «Матрица инверсная» . mathworld.wolfram.com . Проверено 8 сентября 2020 .
- ^ Вайсштейн, Эрик В. «Теорема об обратимой матрице» . mathworld.wolfram.com . Проверено 8 сентября 2020 .
- ^ Хорн, Роджер А .; Джонсон, Чарльз Р. (1985). Матричный анализ . Издательство Кембриджского университета . п. 14. ISBN 978-0-521-38632-6..
- ^ Пан, Виктор; Рейф, Джон (1985), Эффективное параллельное решение линейных систем , Труды 17-го ежегодного симпозиума ACM по теории вычислений, Providence: ACM
- ^ Пан, Виктор; Рейф, Джон (1985), Центр исследований в области вычислительной техники Гарвардского университета, отчет TR-02-85 , Кембридж, Массачусетс: вычислительная лаборатория Эйкена
- ^ «Обращение больших матриц». Байт Журнал . 11 (4): 181–190. Апрель 1986 г.
- ^ Доказательство можно найти в Приложении B к Кондратюк, Л.А.; Криворученко М.И. (1992). «Сверхпроводящее кварковое вещество в цветовой группе SU (2)» . Zeitschrift für Physik . 344 : 99–115. DOI : 10.1007 / BF01291027 . S2CID 120467300 .
- ^ Стрэнг, Гилберт (2003). Введение в линейную алгебру (3-е изд.). СИАМ. п. 71. ISBN 978-0-9614088-9-3., Глава 2, страница 71
- ^ Бернштейн, Деннис (2005). Матричная математика . Издательство Принстонского университета. п. 44. ISBN 978-0-691-11802-4.
- ^ Бернштейн, Деннис (2005). Матричная математика . Издательство Принстонского университета. п. 45. ISBN 978-0-691-11802-4.
- ^ TH Cormen, CE Leiserson, RL Rivest, C. Stein, Введение в алгоритмы , 3-е изд., MIT Press, Кембридж, Массачусетс, 2009, §28.2.
- ↑ Ran Raz . О сложности матричного произведения. В материалах тридцать четвертого ежегодного симпозиума ACM по теории вычислений. ACM Press, 2002. DOI : 10,1145 / 509907,509932 .
- ^ Стюарт, Гилберт (1998). Матричные алгоритмы: основные разложения . СИАМ. п. 55. ISBN 978-0-89871-414-2.
- ^ Haramoto, H .; Мацумото, М. (2009). «P-адический алгоритм вычисления инверсии целочисленных матриц» . Журнал вычислительной и прикладной математики . 225 : 320–322. DOI : 10.1016 / j.cam.2008.07.044 .
- ^ «IML - Библиотека целочисленных матриц» . cs.uwaterloo.ca . Проверено 14 апреля 2018 года .
- ^ Магнус, Ян Р .; Neudecker, Хайнц (1999). Матричное дифференциальное исчисление: с приложениями в статистике и эконометрике (пересмотренное издание). Нью-Йорк: Джон Вили и сыновья. С. 151–152. ISBN 0-471-98633-X.
- ^ Линь, Линь; Лу, Цзяньфэн; Инь, Лексинг; Автомобиль, Роберто; E, Weinan (2009). «Быстрый алгоритм выделения диагонали обратной матрицы с приложением к анализу электронной структуры металлических систем» . Сообщения в математических науках . 7 (3): 755–777. DOI : 10.4310 / CMS.2009.v7.n3.a12 .
дальнейшее чтение
- "Обращение матрицы" , Энциклопедия математики , EMS Press , 2001 [1994]
- Кормен, Томас Х .; Лейзерсон, Чарльз Э .; Ривест, Рональд Л .; Стейн, Клиффорд (2001) [1990]. «28.4: Инвертирование матриц». Введение в алгоритмы (2-е изд.). MIT Press и McGraw-Hill. С. 755–760. ISBN 0-262-03293-7.
- Бернштейн, Деннис С. (2009). Матричная математика: теория, факты и формулы (2-е изд.). Издательство Принстонского университета. ISBN 978-0691140391- через Google Книги .
- Петерсен, Кааре Брандт; Педерсен, Майкл Сискинд (15 ноября 2012 г.). "Поваренная книга Матрицы" (PDF) . С. 17–23.
Внешние ссылки
- Сандерсон, Грант (15 августа 2016 г.). «Обратные матрицы, пространство столбцов и пустое пространство» . Суть линейной алгебры - через YouTube .
- Стрэнг, Гилберт. «Лекция по линейной алгебре по обратным матрицам» . MIT OpenCourseWare .
- Символьный калькулятор обратной матрицы с указанными шагами
- Обратная матрица Мура-Пенроуза