Ортогональная матрица

Для матриц с ортогональностью по полю комплексных чисел см. унитарную матрицу .

В линейной алгебре или ортогональная матрица ортонормированная матрица — это действительная квадратная матрица , столбцы и строки которой являются ортонормированными векторами .

Один из способов выразить это $$Q^{\mathrm {T} }Q=QQ^{\mathrm {T} }=I,$$где Q ^Т — транспонирование Q , а I — единичная матрица .

Это приводит к эквивалентной характеристике: матрица Q ортогональна, если ее транспонирование равно обратному : $$Q^{\mathrm {T} }=Q^{-1},$$где Q ⁻¹ является обратным Q .

Ортогональная матрица Q обязательно обратима (с обратным Q ⁻¹ = К ^Т ), унитарный ( Q ⁻¹ = К ^∗ ), где Q ^∗ является эрмитовым сопряженным ( сопряженным транспонированием ) к Q и, следовательно, нормальным ( Q ^∗ Q = QQ ^∗ ) над действительными числами . Определитель любой ортогональной матрицы равен +1 или -1. В качестве линейного преобразования ортогональная матрица сохраняет скалярное произведение векторов и, следовательно, действует как изометрия евклидова пространства , например вращение , отражение или роторное отражение . Другими словами, это унитарное преобразование .

Набор ортогональных матриц размера n × n при умножении образует группу O( n ) , известную как ортогональная группа . Подгруппа состоящая SO( n ), из ортогональных матриц с определителем +1, называется специальной ортогональной группой , а каждый ее элемент является специальной ортогональной матрицей. В качестве линейного преобразования каждая специальная ортогональная матрица действует как вращение.

Обзор

Визуальное понимание умножения путем транспонирования матрицы. Если A — ортогональная матрица, а B — ее транспонированная, ij-й элемент произведения AA ^Т исчезнет, ​​если i≠j, потому что i-я строка A ортогональна j-й строке A.

Ортогональная матрица — это реальная специализация унитарной матрицы и, следовательно, всегда нормальная матрица . Хотя мы рассматриваем здесь только вещественные матрицы, определение можно использовать для матриц с элементами из любого поля . Однако ортогональные матрицы естественным образом возникают из скалярных произведений , а для матриц комплексных чисел это вместо этого приводит к унитарному требованию. Ортогональные матрицы сохраняют скалярное произведение, Итак, для векторов u и v в n -мерном вещественном евклидовом пространстве $${\mathbf {u} }\cdot {\mathbf {v} }=\left(Q{\mathbf {u} }\right)\cdot \left(Q{\mathbf {v} }\right)$$где Q — ортогональная матрица. Чтобы увидеть связь внутреннего произведения, рассмотрим вектор v в n -мерном реальном евклидовом пространстве . Если записать относительно ортонормированного базиса, квадрат длины v равен v ^Т в . Если линейное преобразование в матричной форме Q v сохраняет длины векторов, то $${\mathbf {v} }^{\mathrm {T} }{\mathbf {v} }=(Q{\mathbf {v} })^{\mathrm {T} }(Q{\mathbf {v} })={\mathbf {v} }^{\mathrm {T} }Q^{\mathrm {T} }Q{\mathbf {v} }.$$

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

Ортогональные матрицы важны по ряду причин, как теоретических, так и практических. × Ортогональные матрицы размера n n образуют группу при умножении матриц, ортогональную группу , обозначаемую O( n ) , которая вместе со своими подгруппами широко используется в математике и физических науках. Например, точечная группа молекулы является подгруппой O(3). Поскольку версии ортогональных матриц с плавающей запятой обладают выгодными свойствами, они являются ключевыми для многих алгоритмов числовой линейной алгебры, таких как QR- разложение . Другой пример: при соответствующей нормализации дискретное косинусное преобразование (используемое при сжатии MP3 ) представляется ортогональной матрицей.

Примеры

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

• ${\begin{bmatrix}1&0\\0&1\\\end{bmatrix}}$ (трансформация личности)
• ${\begin{bmatrix}\cos \theta &-\sin \theta \\\sin \theta &\cos \theta \\\end{bmatrix}}$ (вращение вокруг начала координат)
• ${\begin{bmatrix}1&0\\0&-1\\\end{bmatrix}}$ (отражение по оси X )
• ${\begin{bmatrix}0&0&0&1\\0&0&1&0\\1&0&0&0\\0&1&0&0\end{bmatrix}}$ (перестановка осей координат)

Элементарные конструкции

Нижние размеры

Простейшими ортогональными матрицами являются матрицы 1 × 1 [1] и [−1], которые мы можем интерпретировать как тождество и отражение реальной линии через начало координат.

Матрицы 2 × 2 имеют вид $${\begin{bmatrix}p&t\\q&u\end{bmatrix}},$$требования ортогональности удовлетворяют трем уравнениям $${\begin{aligned}1&=p^{2}+t^{2},\\1&=q^{2}+u^{2},\\0&=pq+tu.\end{aligned}}$$

При рассмотрении первого уравнения без ограничения общности положим p = cos θ , q = sin θ ; тогда либо т знак равно - q , ты знак равно п , либо т знак равно q , ты знак равно - п . Мы можем интерпретировать первый случай как поворот на θ (где θ = 0 — единица), а второй — как отражение через линию под углом ⁠ θ / 2 ⁠ .

$${\begin{bmatrix}\cos \theta &-\sin \theta \\\sin \theta &\cos \theta \\\end{bmatrix}}{\text{ (rotation), }}\qquad {\begin{bmatrix}\cos \theta &\sin \theta \\\sin \theta &-\cos \theta \\\end{bmatrix}}{\text{ (reflection)}}$$

Особый случай матрицы отражения с θ = 90° генерирует отражение относительно линии под углом 45°, заданной y = x, и, следовательно, меняет местами x и y ; это матрица перестановок с одной единицей в каждом столбце и строке (в противном случае 0): $${\begin{bmatrix}0&1\\1&0\end{bmatrix}}.$$

Идентичность также является матрицей перестановок.

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

Высшие измерения

Независимо от размерности всегда можно классифицировать ортогональные матрицы как чисто вращательные или нет, но для матриц 3 × 3 и больше невращательные матрицы могут быть более сложными, чем отражения. Например, $${\begin{bmatrix}-1&0&0\\0&-1&0\\0&0&-1\end{bmatrix}}{\text{ and }}{\begin{bmatrix}0&-1&0\\1&0&0\\0&0&-1\end{bmatrix}}$$

представляют собой инверсию через начало координат и ротоинверсию соответственно вокруг оси z .

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

Однако у нас есть элементарные строительные блоки для перестановок, отражений и вращений, которые применимы в целом.

Примитивы

Самая элементарная перестановка — это транспозиция, полученная из единичной матрицы путем замены двух строк. Любая матрица перестановок размера n × n может быть построена как произведение не более чем n − 1 транспозиций.

Отражение Хаусхолдера строится из ненулевого вектора v как $$Q=I-2{\frac {{\mathbf {v} }{\mathbf {v} }^{\mathrm {T} }}{{\mathbf {v} }^{\mathrm {T} }{\mathbf {v} }}}.$$

Здесь числитель — симметричная матрица, а знаменатель — число, квадрат величины v . Это отражение в гиперплоскости, перпендикулярной v (отрицающее любую компоненту вектора, параллельную v ). Если v — единичный вектор, то Q = I − 2 vv ^Т достаточно. Отражение Хаусхолдера обычно используется для одновременного обнуления нижней части колонны. Любая ортогональная матрица размера n × n может быть построена как произведение не более n таких отражений.

действует Вращение Гивенса на двумерное (плоское) подпространство, охватываемое двумя координатными осями, вращающееся на выбранный угол. Обычно он используется для обнуления одной субдиагональной записи. Любая матрица вращения размера n × n может быть построена как произведение не более ⁠ n ( n − 1) / 2 ⁠ таких вращений. В случае 3×3 матриц достаточно трех таких поворотов; и, зафиксировав последовательность, мы можем, таким образом, описать все матрицы вращения 3 × 3 (хотя и не однозначно) в терминах трех используемых углов, часто называемых углами Эйлера .

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

Характеристики

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

Действительная квадратная матрица ортогональна тогда и только тогда, когда ее столбцы образуют ортонормированный базис евклидова пространства R. ^н с обычным евклидовым скалярным произведением , что имеет место тогда и только тогда, когда его строки образуют ортонормированный базис R ^н . Может возникнуть соблазн предположить, что матрица с ортогональными (не ортонормированными) столбцами будет называться ортогональной матрицей, но такие матрицы не представляют особого интереса и не имеют специального названия; они удовлетворяют только М ^Т M = D , где D — матрица диагональная .

Определитель . любой ортогональной матрицы равен +1 или -1 Это следует из следующих основных фактов об детерминантах: $$1=\det(I)=\det \left(Q^{\mathrm {T} }Q\right)=\det \left(Q^{\mathrm {T} }\right)\det(Q)={\bigl (}\det(Q){\bigr )}^{2}.$$

Обратное неверно; наличие определителя ±1 не является гарантией ортогональности даже для ортогональных столбцов, как показывает следующий контрпример. $${\begin{bmatrix}2&0\\0&{\frac {1}{2}}\end{bmatrix}}$$

В матрицах перестановок определитель соответствует сигнатуре и равен +1 или -1, поскольку четность перестановки четная или нечетная, поскольку определитель является альтернативной функцией строк.

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

Свойства группы

Обратная каждая ортогональная матрица снова ортогональна, как и матричное произведение двух ортогональных матриц. Фактически, набор всех ортогональных матриц размера n × n удовлетворяет всем аксиомам группы . Это компактная группа Ли размерности ⁠ n ( n - 1) / 2 ⁠ , называемая ортогональной группой и обозначаемая O( n ) .

Ортогональные матрицы, определитель которых равен +1, образуют -связную нормальную подгруппу O ( n ) индекса линейно 2, специальную ортогональную группу SO( n ) вращений. Факторгруппа с проекционным отображением, выбирающим O( n )/SO( n ) изоморфна O(1) [+1] или [−1] в соответствии с определителем. Ортогональные матрицы с определителем −1 не содержат единицы и поэтому не образуют подгруппу, а только смежный класс ; он также (отдельно) подключен. Таким образом, каждая ортогональная группа распадается на две части; и поскольку карта проекции разделяется , O( n ) является полупрямым произведением SO ( n ) на O(1) . С практической точки зрения аналогичное утверждение состоит в том, что любую ортогональную матрицу можно получить, взяв матрицу вращения и, возможно, отрицая один из ее столбцов, как мы видели на примере матриц 2 × 2 . Если n нечетно, то полупрямое произведение на самом деле является прямым произведением , и любую ортогональную матрицу можно получить, взяв матрицу вращения и, возможно, отрицая все ее столбцы. Это следует из свойства определителей, согласно которому отрицание столбца отрицает определитель, и, таким образом, отрицание нечетного (но не четного) числа столбцов отрицает определитель.

Теперь рассмотрим ( n + 1) × ( n + 1) ортогональные матрицы с нижней правой записью, равной 1. Остаток последнего столбца (и последней строки) должен быть нулевым, а произведение любых двух таких матриц имеет тот же вид . Остальная часть матрицы представляет собой ортогональную матрицу размера n × n ; таким образом, O( n ) является подгруппой O( n + 1) (и всех более высоких групп).

$${\begin{bmatrix}&&&0\\&\mathrm {O} (n)&&\vdots \\&&&0\\0&\cdots &0&1\end{bmatrix}}$$

Поскольку элементарное отражение в форме матрицы Хаусхолдера может привести любую ортогональную матрицу к этой ограниченной форме, серия таких отражений может привести любую ортогональную матрицу к единице; таким образом, ортогональная группа является группой отражения . Последний столбец может быть привязан к любому единичному вектору, и каждый выбор дает другую копию O( n ) в O( n + 1) ; таким образом, O( n + 1) является расслоением над единичной сферой S ^н с волокном O( n ) .

Аналогично, SO( n ) является подгруппой SO( n + 1) ; и любая специальная ортогональная матрица может быть создана путем вращения плоскости Гивенса с использованием аналогичной процедуры. Структура расслоения сохраняется: SO( n ) ↪ SO( n + 1) → S ^н . Одно вращение может дать ноль в первой строке последнего столбца, а серия из n - 1 вращений обнулит все, кроме последней строки последнего столбца матрицы вращения n × n . Поскольку плоскости фиксированы, каждое вращение имеет только одну степень свободы — угол. Следовательно , по индукции SO( n ) имеет $$(n-1)+(n-2)+\cdots +1={\frac {n(n-1)}{2}}$$степеней свободы, а также O( n ) .

Матрицы перестановок еще проще; они образуют не группу Ли, а только конечную группу порядка n ! симметричная группа S _n . тем же соображениям Sn ₊ является подгруппой Sn _{По 1} . Четные перестановки образуют подгруппу матриц перестановок определителя +1 порядка ⁠ н ! / 2 ⁠ чередующаяся группа .

Каноническая форма

В более широком смысле, эффект любой ортогональной матрицы разделяется на независимые действия на ортогональных двумерных подпространствах. То есть, если Q является специальной ортогональной, то всегда можно найти ортогональную матрицу P , (вращательную) замену базиса , которая приводит Q в блочно-диагональную форму:

$$P^{\mathrm {T} }QP={\begin{bmatrix}R_{1}&&\\&\ddots &\\&&R_{k}\end{bmatrix}}\ (n{\text{ even}}),\ P^{\mathrm {T} }QP={\begin{bmatrix}R_{1}&&&\\&\ddots &&\\&&R_{k}&\\&&&1\end{bmatrix}}\ (n{\text{ odd}}).$$

где матрицы R ₁ , ..., R _k представляют собой матрицы вращения 2 × 2 , а остальные элементы равны нулю. В исключительных случаях блок вращения может быть диагональным, ± I . Таким образом, отрицая при необходимости один столбец и отмечая, что отражение 2 × 2 диагонализируется до +1 и -1, любую ортогональную матрицу можно привести к виду $$P^{\mathrm {T} }QP={\begin{bmatrix}{\begin{matrix}R_{1}&&\\&\ddots &\\&&R_{k}\end{matrix}}&0\\0&{\begin{matrix}\pm 1&&\\&\ddots &\\&&\pm 1\end{matrix}}\\\end{bmatrix}},$$

Матрицы R1 _; ,..., Rk _, лежащих на единичной окружности комплексной плоскости дают сопряженные пары собственных значений таким образом, это разложение подтверждает, что все собственные значения имеют абсолютное значение 1. Если n нечетно, существует хотя бы одно действительное собственное значение, +1 или -1; для вращения 3 × 3 собственный вектор, связанный с +1, является осью вращения.

Алгебра Ли

что элементы Q являются дифференцируемыми функциями t , и что t = 0 дает Q = I. Предположим , Дифференцирование условия ортогональности $$Q^{\mathrm {T} }Q=I$$урожайность $${\dot {Q}}^{\mathrm {T} }Q+Q^{\mathrm {T} }{\dot {Q}}=0$$

Тогда оценка при t = 0 ( Q = I ) подразумевает $${\dot {Q}}^{\mathrm {T} }=-{\dot {Q}}.$$

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

Например, физика трехмерных объектов называет угловую скорость дифференциальным вращением, то есть вектором в алгебре Ли. ${\mathfrak {so}}(3)$ касательная к SO(3) . Учитывая ω = ( xθ , yθ , zθ ) , где v = ( x , y , z ) является единичным вектором, правильная кососимметричная матричная форма ω равна $$\Omega ={\begin{bmatrix}0&-z\theta &y\theta \\z\theta &0&-x\theta \\-y\theta &x\theta &0\end{bmatrix}}.$$

Экспонента этого является ортогональной матрицей для вращения вокруг оси v на угол θ ; установка c = cos ⁠ θ / 2 ⁠ , s = грех ⁠ θ / 2 ⁠ , $$\exp(\Omega )={\begin{bmatrix}1-2s^{2}+2x^{2}s^{2}&2xys^{2}-2zsc&2xzs^{2}+2ysc\\2xys^{2}+2zsc&1-2s^{2}+2y^{2}s^{2}&2yzs^{2}-2xsc\\2xzs^{2}-2ysc&2yzs^{2}+2xsc&1-2s^{2}+2z^{2}s^{2}\end{bmatrix}}.$$

Численная линейная алгебра

Преимущества

Численный анализ использует многие свойства ортогональных матриц числовой линейной алгебры, и они возникают естественным образом. Например, часто желательно вычислить ортонормированный базис пространства или ортогональную замену базисов; оба принимают форму ортогональных матриц. Наличие определителя ±1 и всех собственных значений величины 1 очень полезно для числовой стабильности . Одним из следствий является то, что число обусловленности равно 1 (что является минимумом), поэтому ошибки не увеличиваются при умножении на ортогональную матрицу. многие алгоритмы используют ортогональные матрицы, такие как отражения Хаусхолдера и вращения Гивенса По этой причине . Полезно также то, что ортогональная матрица не только обратима, но и то, что ее обратная матрица доступна практически бесплатно путем обмена индексами.

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

Аналогичным образом, алгоритмы, использующие матрицы Хаусхолдера и Гивенса, обычно используют специализированные методы умножения и хранения. Например, поворот Гивенса влияет только на две строки матрицы, которую он умножает, изменяя полное умножение порядка n. ³ до гораздо более эффективного порядка n . Когда использование этих отражений и вращений приводит к появлению нулей в матрице, освободившегося пространства достаточно для хранения достаточного количества данных для воспроизведения преобразования и для надежного выполнения этого преобразования. (Следуя Стюарту (1976) , мы не сохраняем угол поворота, что дорого и плохо ведет себя.)

Разложения

Ряд важных матричных разложений ( Golub & Van Loan 1996 ) используют ортогональные матрицы, в том числе, в частности:

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

M = QR , Q ортогональный, R верхний треугольный

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

М = U Σ V ^Т , U и V ортогональны, Σ диагональная матрица

Собственное разложение симметричной матрицы (разложение по спектральной теореме )

S = Q Λ Q ^Т , S симметричный, Q ортогональный, Λ диагональный

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

M = QS , Q ортогональный, S симметричный положительно-полуопределенный

Примеры

Рассмотрим переопределенную систему линейных уравнений , что может произойти при повторных измерениях физического явления для компенсации экспериментальных ошибок. Запишите A x знак равно b , где A равно m × n , m > n .QR - разложение сводит A к верхнему R. треугольному Например, если A равно 5 × 3 , то R имеет вид $$R={\begin{bmatrix}\cdot &\cdot &\cdot \\0&\cdot &\cdot \\0&0&\cdot \\0&0&0\\0&0&0\end{bmatrix}}.$$

Линейная задача наименьших квадратов состоит в том, чтобы найти x , который минимизирует ‖ A x − b ‖ , что эквивалентно проецированию b в подпространство, охватываемое столбцами A . Предполагая, что столбцы A (и, следовательно, R ) независимы, проекционное решение находится из A ^Т А х = А ^Т б . Теперь А ^Т A квадратный ( n × n ) и обратимый, а также равен R ^Т Р. ​Но нижние строки нулей в R являются лишними в произведении, которое, таким образом, уже находится в нижнетреугольной верхнетреугольной факторизованной форме, как при методе исключения Гаусса ( разложение Холецкого ). Здесь ортогональность важна не только для уменьшения A ^Т А = ( Р ^Т вопрос ^Т ) QR - R ^Т R , но и за возможность решения без увеличения числовых задач.

В случае линейной системы, которая недоопределена, или в противном случае необратимой матрицы , разложение по сингулярным значениям (SVD) одинаково полезно. С учетом A как U Σ V ^Т , удовлетворительное решение использует псевдообратную функцию Мура-Пенроуза , V Σ ⁺ В ^Т , где Σ ⁺ просто заменяет каждую ненулевую диагональную запись на обратную. Установите x на V Σ ⁺ В ^Т б .

Интересен также случай квадратной обратимой матрицы. Предположим, например, что A — это матрица вращения размером 3 × 3 , которая была вычислена как композиция многочисленных поворотов и поворотов. Число с плавающей запятой не соответствует математическому идеалу действительных чисел, поэтому A постепенно утратило свою истинную ортогональность. Процесс Грама -Шмидта мог бы ортогонализировать столбцы, но это не самый надежный, не самый эффективный и не самый инвариантный метод. Полярное разложение разлагает матрицу на пару, одна из которых является уникальной, ближайшей ортогональной матрицей к данной матрице, или одной из ближайших, если данная матрица сингулярна. (Близость может быть измерена с помощью любой матричной нормы, инвариантной относительно ортогональной замены базиса, такой как спектральная норма или норма Фробениуса.) Для почти ортогональной матрицы быстрая сходимость к ортогональному множителю может быть достигнута с помощью « метода Ньютона ». подход Хайэма (1986) ( 1990 ), многократно усредняющий матрицу с ее обратным транспонированием. Дубрюль (1999) опубликовал ускоренный метод с удобным тестом сходимости.

Например, рассмотрим неортогональную матрицу, для которой простой алгоритм усреднения занимает семь шагов. $${\begin{bmatrix}3&1\\7&5\end{bmatrix}}\rightarrow {\begin{bmatrix}1.8125&0.0625\\3.4375&2.6875\end{bmatrix}}\rightarrow \cdots \rightarrow {\begin{bmatrix}0.8&-0.6\\0.6&0.8\end{bmatrix}}$$и какое ускорение сокращается до двух ступеней (при γ = 0,353553, 0,565685).

$${\begin{bmatrix}3&1\\7&5\end{bmatrix}}\rightarrow {\begin{bmatrix}1.41421&-1.06066\\1.06066&1.41421\end{bmatrix}}\rightarrow {\begin{bmatrix}0.8&-0.6\\0.6&0.8\end{bmatrix}}$$

Грамм-Шмидт дает худшее решение, демонстрируемое расстоянием Фробениуса 8,28659 вместо минимального 8,12404.

$${\begin{bmatrix}3&1\\7&5\end{bmatrix}}\rightarrow {\begin{bmatrix}0.393919&-0.919145\\0.919145&0.393919\end{bmatrix}}$$

Рандомизация

Некоторые численные приложения, такие как методы Монте-Карло и исследование многомерных пространств данных, требуют создания равномерно распределенных случайных ортогональных матриц. В этом контексте «равномерность» определяется с точки зрения меры Хаара , которая по существу требует, чтобы распределение не менялось при умножении на любую свободно выбранную ортогональную матрицу. Ортогонализация матриц с независимыми равномерно распределенными случайными элементами не приводит к получению равномерно распределенных ортогональных матриц. , но QR- разложение независимых нормально распределенных случайных элементов делает это, пока диагональ R содержит только положительные элементы ( Mezzadri 2006 ). Стюарт (1980) заменил это более эффективной идеей, которую Диаконис и Шахшахани (1987) позже обобщили как «алгоритм подгруппы» (в этой форме он работает так же хорошо для перестановок и вращений). Чтобы сгенерировать ортогональную матрицу ( n + 1) × ( n + 1) , возьмите матрицу n × n и равномерно распределенный единичный вектор размерности n + 1 . Постройте отражение Хаусхолдера из вектора, затем примените его к меньшей матрице (встроенной в матрицу большего размера с цифрой 1 в правом нижнем углу).

Ближайшая ортогональная матрица

Проблема нахождения ортогональной матрицы Q, ближайшей к данной матрице M, связана с ортогональной проблемой Прокруста . Существует несколько различных способов получить уникальное решение, самый простой из которых — разложить сингулярным значениям M по и заменить сингулярные значения единицами. Другой метод явно выражает R , но требует использования квадратного корня матрицы : $$Q=M\left(M^{\mathrm {T} }M\right)^{-{\frac {1}{2}}}$$

Это можно объединить с вавилонским методом извлечения квадратного корня из матрицы, чтобы получить рекурсию, которая квадратично сходится к ортогональной матрице: $$Q_{n+1}=2M\left(Q_{n}^{-1}M+M^{\mathrm {T} }Q_{n}\right)^{-1}$$где Q ₀ = М .

Эти итерации стабильны, если обусловленности M число меньше трех.

Использование обратного приближения первого порядка и той же инициализации приводит к модифицированной итерации:

$$N_{n}=Q_{n}^{\mathrm {T} }Q_{n}$$$$P_{n}={\frac {1}{2}}Q_{n}N_{n}$$$$Q_{n+1}=2Q_{n}+P_{n}N_{n}-3P_{n}$$

Вращение и булавка

В некоторых случаях использования ортогональных матриц возникает тонкая техническая проблема. Компоненты группы с определителем +1 и −1 не только не связаны друг с другом, но даже компонент +1, SO( n ) , не является односвязным (за исключением SO(1), что тривиально). Таким образом, иногда полезно или даже необходимо работать с покрывающей группой SO( n ), группой спиновой Spin( n ) . Аналогично, O( n ) имеет покрывающие группы, группы выводов Pin( n ). Для n > 2 Spin ( n ) односвязен и, следовательно, является универсальной накрывающей группой для SO( n ) . Безусловно, самым известным примером спиновой группы является Spin(3) , которая представляет собой не что иное, как SU(2) или группу единичных кватернионов .

Группы Pin и Spin находятся в алгебрах Клиффорда , которые сами могут быть построены из ортогональных матриц.

Прямоугольные матрицы

Основная статья: Полуортогональная матрица

Если Q не является квадратной матрицей, то условия Q ^Т Q = I и QQ ^Т = Я не эквивалентен. Условие Q ^Т Q = I говорит, что столбцы Q ортонормированы. Это может произойти только в том случае, если Q — матрица размера m × n с n ≤ m (из-за линейной зависимости). Аналогично, QQ ^Т = Я говорю, что строки Q ортонормированы, что требует n ≥ m .

Для этих матриц не существует стандартной терминологии. Их по-разному называют «полуортогональными матрицами», «ортонормированными матрицами», «ортогональными матрицами», а иногда просто «матрицами с ортонормированными строками/столбцами».

В случае n ≤ m матрицы с ортонормированными столбцами можно назвать ортогональными k-шкалами и они являются элементами многообразия Штифеля .

См. также

• Биортогональная система

Ссылки

• Диаконис, Персия ; Шахшахани, Мехрдад (1987), «Алгоритм подгруппы для генерации однородных случайных величин», Вероятность в инженерных и информационных науках , 1 : 15–32, doi : 10.1017/S0269964800000255 , ISSN   0269-9648 , S2CID   122752374
• Дюбрюлль, Огюстен А. (1999), «Оптимальная итерация для полярного разложения матрицы» , Электронные транзакции по численному анализу , 8 : 21–25
• Голуб, Джин Х .; Ван Лоан, Чарльз Ф. (1996), Матричные вычисления (3-е изд.), Балтимор: Издательство Университета Джонса Хопкинса, ISBN  978-0-8018-5414-9
• Хайэм, Николас (1986), «Вычисление полярного разложения — с приложениями» (PDF) , SIAM Journal on Scientific and Statistical Computing , 7 (4): 1160–1174, doi : 10.1137/0907079 , ISSN   0196-5204
• Хайэм, Николас ; Шрайбер, Роберт (июль 1990 г.), «Быстрое полярное разложение произвольной матрицы», SIAM Journal on Scientific and Statistical Computing , 11 (4): 648–655, CiteSeerX   10.1.1.230.4322 , doi : 10.1137/0911038 , ISSN   0196 -5204 , S2CID   14268409 [1]
• Стюарт, Г.В. (1976), «Экономическое хранение вращений плоскостей», Numerische Mathematik , 25 (2): 137–138, doi : 10.1007/BF01462266 , ISSN   0029-599X , S2CID   120372682
• Стюарт, Г.В. (1980), «Эффективное создание случайных ортогональных матриц с применением к средствам оценки условий», SIAM Journal on Numerical Analysis , 17 (3): 403–409, Bibcode : 1980SJNA...17..403S , doi : 10.1137/0717034 , ISSN   0036-1429
• Меццадри, Франческо (2006), «Как генерировать случайные матрицы из классических компактных групп», Уведомления Американского математического общества , 54 , arXiv : math-ph/0609050 , Bibcode : 2006math.ph...9050M

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

• «Ортогональная матрица» , Математическая энциклопедия , EMS Press , 2001 [1994]
• Учебно-интерактивная программа по ортогональной матрице