Разложение матрицы

редактировать
Не путать с матричной факторизацией многочлена.

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

СОДЕРЖАНИЕ

  • 1 Пример
  • 2 Разложения, связанные с решением систем линейных уравнений
    • 2.1 LU разложение
    • 2.2 сокращение LU
    • 2.3 Блочная декомпозиция LU
    • 2.4 Факторизация рангов
    • 2.5 Разложение Холецкого
    • 2.6 QR-разложение
    • 2.7 Факторизация RRQR
    • 2.8 Интерполяционная декомпозиция
  • 3 Разложения на основе собственных значений и связанных понятий
    • 3.1 Собственное разложение
    • 3.2 Жорданов разложение
    • 3.3 Разложение Шура
    • 3.4 Вещественное разложение Шура
    • 3.5 QZ разложение
    • 3.6 Факторизация Такаги
    • 3.7 Разложение по сингулярным числам
    • 3.8 Масштабно-инвариантные разложения
  • 4 Другие разложения
    • 4.1 Полярное разложение
    • 4.2 Алгебраическое полярное разложение
    • 4.3 Разложение Мостова
    • 4.4 Нормальная форма синкхорна
    • 4.5 Секторальная декомпозиция
    • 4.6 нормальная форма Вильямсона
    • 4.7 Матричный квадратный корень
  • 5 Обобщения
  • 6 См. Также
  • 7 ссылки
    • 7.1 Примечания
    • 7.2 Цитаты
    • 7.3 Библиография
  • 8 Внешние ссылки

Пример

В численном анализе используются различные декомпозиции для реализации эффективных матричных алгоритмов.

Например, при решении системы линейных уравнений матрица A может быть разложена посредством LU-разложения. Разложение LU факторизация матрицы в нижний треугольную матрице L и верхнюю треугольную матрицу U. Системы и требуют меньшего количества сложений и умножений для решения по сравнению с исходной системой, хотя может потребоваться значительно больше цифр в неточной арифметике, такой как числа с плавающей запятой. А Икс знак равно б {\ Displaystyle А \ mathbf {x} = \ mathbf {b}} L ( U Икс ) знак равно б {\ Displaystyle L (U \ mathbf {x}) = \ mathbf {b}} U Икс знак равно L - 1 б {\ Displaystyle U \ mathbf {x} = L ^ {- 1} \ mathbf {b}} А Икс знак равно б {\ Displaystyle А \ mathbf {x} = \ mathbf {b}}

Точно так же разложение QR выражает A как QR, где Q - ортогональная матрица, а R - верхнетреугольная матрица. Система Q ( R x ) = b решается с помощью R x = Q T b = c, а система R x = c решается с помощью « обратной подстановки ». Количество требуемых сложений и умножений примерно в два раза больше, чем при использовании решателя LU, но в неточной арифметике больше не требуется цифр, поскольку QR-разложение численно стабильно.

Разложения, связанные с решением систем линейных уравнений

LU разложение

Основная статья: LU разложение

Сокращение LU

Основная статья: сокращение LU

Блочная декомпозиция LU

Основная статья: Блочная декомпозиция LU

Факторизация рангов

Основная статья: Факторизация рангов

Разложение Холецкого

Основная статья: разложение Холецкого
  • Применимо к: квадратной, эрмитовой, положительно определенной матрице A
  • Разложение:, где - верхний треугольник с вещественными положительными диагональными элементами А знак равно U * U {\ displaystyle A = U ^ {*} U} U {\ displaystyle U}
  • Комментарий: если матрица эрмитова и положительно полуопределенная, то она имеет разложение вида, если диагональные элементы матрицы могут быть равны нулю. А {\ displaystyle A} А знак равно U * U {\ displaystyle A = U ^ {*} U} U {\ displaystyle U}
  • Единственность: для положительно определенных матриц разложение Холецкого единственно. Однако в положительном полуопределенном случае он не уникален.
  • Комментарий: если A действительный и симметричный, имеет все действительные элементы U {\ displaystyle U}
  • Комментарий: Альтернативой является разложение LDL, которое позволяет избежать извлечения квадратных корней.

QR-разложение

Основная статья: QR-разложение
  • Применимо к: m -by- n матрице A с линейно независимыми столбцами
  • Разложение: где является унитарной матрицей размера М матрицы с размерностью м, и является верхней треугольной матрицей размера м матрицы с размерностью п А знак равно Q р {\ displaystyle A = QR} Q {\ displaystyle Q} р {\ displaystyle R}
  • Уникальность: Как правило, он не уникален, но если он полного ранга, то существует единственный, который имеет все положительные диагональные элементы. Если квадрат, то и уникален. А {\ displaystyle A} р {\ displaystyle R} А {\ displaystyle A} Q {\ displaystyle Q}
  • Комментарий: QR-разложение обеспечивает эффективный способ решения системы уравнений. Тот факт, что это ортогональный означает, что, таким образом, что эквивалентно, что очень легко решить, так как это треугольный. А Икс знак равно б {\ Displaystyle А \ mathbf {x} = \ mathbf {b}} Q {\ displaystyle Q} Q Т Q знак равно я {\ Displaystyle Q ^ {\ mathrm {T}} Q = I} А Икс знак равно б {\ Displaystyle А \ mathbf {x} = \ mathbf {b}} р Икс знак равно Q Т б {\ Displaystyle R \ mathbf {x} = Q ^ {\ mathsf {T}} \ mathbf {b}} р {\ displaystyle R}

RRQR факторизация

Основная статья: факторизация RRQR

Интерполяционная декомпозиция

Основная статья: Интерполяционная декомпозиция

Разложения на основе собственных значений и связанных понятий

Собственное разложение

Основная статья: Собственное разложение (матрица)
  • Также называется спектральным разложением.
  • Применимо к: квадратной матрице A с линейно независимыми собственными векторами (не обязательно различными собственными значениями).
  • Разложение:, где D представляет собой диагональная матрица формируется из собственных значений из А, а столбцы V являются соответствующим собственными векторами из А. А знак равно V D V - 1 {\ displaystyle A = VDV ^ {- 1}}
  • Существование: п матрицы с размерностью п матрица всегда имеет п (комплексные) собственные значения, которые могут быть заказаны (в более чем одним способе) с образованием п матрица с размерностью п диагональной матрицей D и соответствующей матрицей ненулевых столбцов V, который удовлетворяет уравнение на собственные значения. обратима тогда и только тогда, когда n собственных векторов линейно независимы (т. е. каждое собственное значение имеет геометрическую кратность, равную его алгебраической кратности ). Достаточным (но не необходимым) условием для этого является то, что все собственные значения различны (в этом случае геометрическая и алгебраическая кратность равны 1) А V знак равно V D {\ displaystyle AV = VD} V {\ displaystyle V}
  • Комментарий: всегда можно нормализовать собственные векторы, чтобы они имели длину, равную единице (см. Определение уравнения для собственных значений)
  • Комментарий: Каждая нормальная матрица A (т. Е. Матрица, для которой, где - сопряженное транспонирование ) может быть разложена на собственные. Для нормальной матрицы A (и только для нормальной матрицы) собственные векторы также можно сделать ортонормированными (), а собственное разложение читается как. В частности, все унитарные, эрмитовы или кососимметричные (в вещественном случае все ортогональные, симметричные или кососимметричные соответственно) матрицы нормальны и, следовательно, обладают этим свойством. А А * знак равно А * А {\ displaystyle AA ^ {*} = A ^ {*} A} А * {\ displaystyle A ^ {*}} V V * знак равно я {\ displaystyle VV ^ {*} = I} А знак равно V D V * {\ displaystyle A = VDV ^ {*}}
  • Комментарий: Для любой вещественной симметричной матрицы A всегда существует собственное разложение, и его можно записать как, где и D, и V являются действительными значениями. А знак равно V D V Т {\ Displaystyle A = VDV ^ {\ mathsf {T}}}
  • Комментарий: Собственное разложение полезно для понимания решения системы линейных обыкновенных дифференциальных уравнений или линейных разностных уравнений. Например, разностное уравнение, начиная с начальным условием решается, что эквивалентно, где V и D являются матрицами, образованные из собственных векторов и собственных значений А. Поскольку D диагональный, возведение его в степень просто включает в себя возведение каждого элемента по диагонали в степень t. Это намного проще сделать и понять, чем возвести A в степень t, поскольку A обычно не диагональна. Икс т + 1 знак равно А Икс т {\ displaystyle x_ {t + 1} = Ax_ {t}} Икс 0 знак равно c {\ displaystyle x_ {0} = c} Икс т знак равно А т c {\ displaystyle x_ {t} = A ^ {t} c} Икс т знак равно V D т V - 1 c {\ displaystyle x_ {t} = VD ^ {t} V ^ {- 1} c} D т {\ displaystyle D ^ {t}}

Разложение Жордана

Жордановой нормальной формы и разложение Жордана-Шевалье

  • Применимо к: квадратной матрице A
  • Комментарий: нормальная форма Жордана обобщает собственное разложение на случаи, когда есть повторяющиеся собственные значения и не могут быть диагонализованы, разложение Жордана – Шевалле делает это без выбора базиса.

Разложение Шура

Основная статья: разложение Шура

Реальное разложение Шура

QZ разложение

Основная статья: QZ разложение
  • Также называется: обобщенное разложение Шура
  • Применимо к: квадратным матрицам A и B
  • Комментарий: есть две версии этой декомпозиции: сложная и реальная.
  • Разложение (комплексная версия): и где Q и Z - унитарные матрицы, верхний индекс * представляет сопряженное транспонирование, а S и T - верхние треугольные матрицы. А знак равно Q S Z * {\ displaystyle A = QSZ ^ {*}} B знак равно Q Т Z * {\ displaystyle B = QTZ ^ {*}}
  • Комментарий: в комплексном разложении QZ отношения диагональных элементов S к соответствующим диагональным элементам T, являются обобщенными собственными значениями, которые решают обобщенную проблему собственных значений (где - неизвестный скаляр, а v - неизвестный ненулевой вектор). λ я знак равно S я я / Т я я {\ displaystyle \ lambda _ {i} = S_ {ii} / T_ {ii}} А v знак равно λ B v {\ Displaystyle A \ mathbf {v} = \ lambda B \ mathbf {v}} λ {\ displaystyle \ lambda}
  • Разложение (реальная версия): и где A, B, Q, Z, S и T - матрицы, содержащие только действительные числа. В этом случае Q и Z - ортогональные матрицы, верхний индекс T представляет транспонирование, а S и T - блочные верхнетреугольные матрицы. Блоки на диагонали S и T имеют размер 1 × 1 или 2 × 2. А знак равно Q S Z Т {\ displaystyle A = QSZ ^ {\ mathsf {T}}} B знак равно Q Т Z Т {\ Displaystyle B = QTZ ^ {\ mathsf {T}}}

Факторизация Такаги

  • Применимо к: квадрат, сложной, симметричной матрицы A.
  • Разложение:, где D - вещественная неотрицательная диагональная матрица, а V - унитарная. обозначает матрицу транспонирование из V. А знак равно V D V Т {\ Displaystyle A = VDV ^ {\ mathsf {T}}} V Т {\ Displaystyle V ^ {\ mathsf {T}}}
  • Комментарий: диагональные элементы D являются неотрицательными квадратными корнями из собственных значений. А А * {\ displaystyle AA ^ {*}}
  • Комментарий: V может быть сложным, даже если A реально.
  • Комментарий: это не частный случай собственного разложения (см. Выше), в котором вместо. Более того, если A не является реальным, это не эрмитово, и форма using также не применяется. V - 1 {\ displaystyle V ^ {- 1}} V Т {\ Displaystyle V ^ {\ mathsf {T}}} V * {\ Displaystyle V ^ {*}}

Разложение по сингулярным числам

Основная статья: Разложение по сингулярным числам
  • Применимо к: м матрицу с размерностью п матрицы А.
  • Разложение:, где D - неотрицательная диагональная матрица, а U и V удовлетворяют. Вот это сопряженное транспонирование из V (или просто транспонированная, если V содержит только действительные числа), и я обозначает единичную матрицу (некоторой размерности). А знак равно U D V * {\ displaystyle A = UDV ^ {*}} U * U знак равно я , V * V знак равно я {\ Displaystyle U ^ {*} U = I, V ^ {*} V = I} V * {\ Displaystyle V ^ {*}}
  • Комментарий: Диагональные элементы D называются сингулярные значения из A.
  • Комментарий: Как и собственное разложение выше, разложение по сингулярным числам включает в себя поиск базисных направлений, по которым умножение матриц эквивалентно скалярному умножению, но оно имеет большую общность, поскольку рассматриваемая матрица не обязательно должна быть квадратной.
  • Уникальность: сингулярные значения всегда определяются однозначно. и не обязательно должен быть уникальным в целом. А {\ displaystyle A} U {\ displaystyle U} V {\ displaystyle V}

Масштабно-инвариантные разложения

Относится к вариантам существующих матричных разложений, таких как SVD, которые инвариантны относительно диагонального масштабирования.

  • Применимо к: м матрицу с размерностью п матрицы А.
  • Блок-масштабно-инвариантная болванка-Значение разложение:, где S является уникальной неотрицательной диагональной матрицей из масштабно-инвариантных сингулярных значений, U и V являются унитарными матрицами, является сопряженной транспозицией из V, и положительные диагональные матрицы D и E. А знак равно D U S V * E {\ displaystyle A = DUSV ^ {*} E} V * {\ Displaystyle V ^ {*}}
  • Комментарий: Аналогичен SVD, за исключением того, что диагональные элементы S инвариантны относительно левого и / или правого умножения A на произвольные невырожденные диагональные матрицы, в отличие от стандартного SVD, для которого сингулярные значения инвариантны относительно левого и / или правое умножение A на произвольные унитарные матрицы.
  • Комментарий: Является альтернативой стандартному СВД, когда требуется инвариантность по отношению к диагонали, а не унитарных преобразований A.
  • Уникальность: масштабно-инвариантные особые значения (заданные диагональными элементами S) всегда определяются однозначно. Диагональные матрицы D и E и унитарные U и V, вообще говоря, не обязательно уникальны. А {\ displaystyle A}
  • Комментарий: матрицы U и V отличаются от матриц SVD.

Аналогичные масштабно-инвариантные разложения могут быть получены из других матричных разложений, например, для получения масштабно-инвариантных собственных значений.

Другие разложения

Полярное разложение

Основная статья: Полярное разложение
  • Применимо к: любой квадратной комплексной матрицы A.
  • Разложение: (правое полярное разложение) или (левое полярное разложение), где U - унитарная матрица, а P и P ' - положительные полуопределенные эрмитовы матрицы. А знак равно U п {\ displaystyle A = UP} А знак равно п U {\ displaystyle A = P'U}
  • Уникальность: всегда уникальна и равна (что всегда эрмитово и положительно полуопределено). Если обратимо, то единственно. п {\ displaystyle P} А * А {\ displaystyle {\ sqrt {A ^ {*} A}}} А {\ displaystyle A} U {\ displaystyle U}
  • Комментарий: Поскольку любая эрмитова матрица допускает спектральное разложение с унитарной матрицей, ее можно записать как. Поскольку положительно полуопределенный, все элементы в неотрицательны. Поскольку произведение двух унитарных матриц унитарно, взяв одну, можно написать, что является разложением по сингулярным значениям. Следовательно, существование полярного разложения эквивалентно существованию разложения по сингулярным значениям. п {\ displaystyle P} п знак равно V D V * {\ displaystyle P = VDV ^ {*}} п {\ displaystyle P} D {\ displaystyle D} W знак равно U V {\ displaystyle W = UV} А знак равно U ( V D V * ) знак равно W D V * {\ Displaystyle A = U (VDV ^ {*}) = WDV ^ {*}}

Алгебраическое полярное разложение

  • Применимо к: квадратной, сложной, невырожденной матрице А.
  • Разложение:, где Q - комплексная ортогональная матрица, а S - комплексная симметричная матрица. А знак равно Q S {\ displaystyle A = QS}
  • Уникальность: если нет отрицательных действительных собственных значений, то разложение уникально. А Т А {\ displaystyle A ^ {\ mathsf {T}} A}
  • Комментарий: Существование этого разложения эквивалентно тому, чтобы быть похожим на. А А Т {\ Displaystyle AA ^ {\ mathsf {T}}} А Т А {\ displaystyle A ^ {\ mathsf {T}} A}
  • Комментарий: вариант этого разложения: где R - вещественная матрица, а C - круговая матрица. А знак равно р C {\ Displaystyle A = RC}

Разложение Мостова

Основная статья: разложение Мостова
  • Применимо к: квадратной, сложной, невырожденной матрице А.
  • Разложение:, где U унитарно, M вещественно антисимметрично, а S вещественно симметрично. А знак равно U е я M е S {\ displaystyle A = Ue ^ {iM} e ^ {S}}
  • Комментарий: Матрица A также может быть разложена как, где U 2 является унитарным, M 2 вещественно антисимметричным и S 2 вещественно симметричным. А знак равно U 2 е S 2 е я M 2 {\ displaystyle A = U_ {2} e ^ {S_ {2}} e ^ {IM_ {2}}}

Рога обыкновенная

Основная статья: теорема Синкхорна
  • Применимо к: квадратной вещественной матрице A со строго положительными элементами.
  • Разложение:, где S - дважды стохастический, а D 1 и D 2 - вещественные диагональные матрицы со строго положительными элементами. А знак равно D 1 S D 2 {\ displaystyle A = D_ {1} SD_ {2}}

Секторальная декомпозиция

  • Применимо к: квадратной комплексной матрице A с числовым диапазоном, содержащимся в секторе. S α знак равно { р е я θ C р gt; 0 , | θ | α lt; π 2 } {\ displaystyle S _ {\ alpha} = \ left \ {re ^ {i \ theta} \ in \ mathbb {C} \ mid rgt; 0, | \ theta | \ leq \ alpha lt;{\ frac {\ pi} { 2}} \ right \}}
  • Разложение:, где C - обратимая комплексная матрица и со всеми. А знак равно C Z C * {\ displaystyle A = CZC ^ {*}} Z знак равно диагональ ( е я θ 1 , , е я θ п ) {\ displaystyle Z = \ operatorname {diag} \ left (e ^ {i \ theta _ {1}}, \ ldots, e ^ {i \ theta _ {n}} \ right)} | θ j | α {\ displaystyle \ left | \ theta _ {j} \ right | \ leq \ alpha}

Нормальная форма Вильямсона

  • Применимо к квадратной положительно определенной вещественной матрице A порядка 2 n × 2 n.
  • Распад:, где есть симплектическая матрица и D является неотрицательным п матрицы с размерностью п диагональной матрицей. А знак равно S Т диагональ ( D , D ) S {\ displaystyle A = S ^ {\ mathsf {T}} \ operatorname {diag} (D, D) S} S Sp ( 2 п ) {\ Displaystyle S \ in {\ text {Sp}} (2n)}

Матричный квадратный корень

Основная статья: квадратный корень из матрицы
  • Разложение:, в общем, не уникальное. А знак равно B B {\ displaystyle A = BB}
  • В случае положительного полуопределенного существует единственное положительное полуопределенное значение, такое что. А {\ displaystyle A} B {\ displaystyle B} А знак равно B * B знак равно B B {\ displaystyle A = B ^ {*} B = BB}

Обобщения

Существуют аналоги SVD, QR, LU и факторизации Холецкого для квазиматриц и c- матриц или непрерывных матриц. «Квазиматрица», как и матрица, представляет собой прямоугольную схему, элементы которой индексированы, но один дискретный индекс заменен непрерывным индексом. Точно так же «c-матрица» непрерывна по обоим индексам. В качестве примера c-матрицы можно представить ядро интегрального оператора.

Эти факторизации основаны на ранних работах Фредгольма (1903), Гильберта (1904) и Шмидта (1907). Отчет и перевод основополагающих статей на английский см. В Stewart (2011).

Смотрите также

Рекомендации

Заметки

Цитаты

Библиография

Внешние ссылки

Последняя правка сделана 2024-01-01 11:50:54
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте