Логарифм матрицы

редактировать

Математические операции с обратимыми матрицами

В математике - логарифм матрицы представляет собой другую матрицу, такую, что экспоненциальная матрица последней матрицы равна исходной матрице. Таким образом, это обобщение скалярного логарифма и в некотором смысле обратная функция для экспоненциальной матрицы. Не все матрицы имеют логарифм, и те матрицы, у которых есть логарифм, могут иметь более одного логарифма. Изучение логарифмов матриц приводит к теории Ли, поскольку, когда матрица имеет логарифм, то она находится в группе Ли, а логарифм является соответствующим элементом векторного пространства Алгебра Ли.

Содержание
  • 1 Определение
  • 2 Выражение степенного ряда
  • 3 Пример: логарифм вращений на плоскости
  • 4 Существование
  • 5 Свойства
  • 6 Дальнейший пример: логарифм вращений в трехмерном пространстве
  • 7 Вычисление логарифма диагонализуемой матрицы
  • 8 Логарифм недиагонализуемой матрицы
  • 9 Перспектива функционального анализа
  • 10 Перспектива теории групп Ли
  • 11 Ограничения в случае 2 × 2
  • 12 См. также
  • 13 Примечания
  • 14 Ссылки
Определение

Экспонента матрицы A определяется как

e A ≡ ∑ n = 0 ∞ A nn! {\ displaystyle e ^ {A} \ Equiv \ sum _ {n = 0} ^ {\ infty} {\ frac {A ^ {n}} {n!}}}{\ displaystyle e ^ {A} \ Equiv \ sum _ {n = 0} ^ {\ infty} {\ frac {A ^ {n} } {n!}}} .

Дана матрица B, другая матрица A называется матричным логарифмом числа B, если e = B. Поскольку экспоненциальная функция не является взаимно однозначной для комплексных чисел (например, e π i = e 3 π i = - 1 {\ displaystyle e ^ {\ pi i} = e ^ {3 \ pi i} = - 1}{\ displaystyle e ^ {\ pi i} = e ^ {3 \ pi i} = - 1} ), числа могут иметь несколько комплексных логарифмов, и, как следствие этого, некоторые матрицы могут иметь больше чем один логарифм, как описано ниже.

Выражение степенного ряда

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

log (B) = ∑ к знак равно 1 ∞ (- 1) k + 1 (B - I) kk = (B - I) - (B - I) 2 2 + (B - I) 3 3 ⋯ {\ displaystyle \ mathrm {log} (B) = \ sum _ {k = 1} ^ {\ infty} {(- 1) ^ {k + 1} {\ frac {(BI) ^ {k}} {k}}} = (BI) - {\ frac {(BI) ^ {2}} {2}} + {\ frac {(BI) ^ {3}} {3}} \ cdots}{\ displaystyle \ mathrm {log} (B) = \ sum _ {k = 1} ^ {\ infty} {(- 1) ^ {k + 1} { \ frac {(BI) ^ {k}} {k}}} = (BI) - {\ frac {(BI) ^ {2}} {2}} + {\ frac {(BI) ^ {3}} {3}} \ cdots} .

В частности, если ‖ B - I ‖ < 1 {\displaystyle \left\|{B-I}\right\|<1}{\ displaystyle \ left \ | {BI} \ right \ | <1} , тогда предыдущий ряд сходится и elog (B) = B {\ displaystyle e ^ {\ mathrm {log} (B)} = B}{\ displaystyle e ^ {\ mathrm {log} (B)} = B} .

Пример: логарифм поворотов в плоскости

Вращения в плоскости дают простой пример. Поворот на угол α вокруг начала координат представлен матрицей 2 × 2

A = (cos ⁡ (α) - sin ⁡ (α) sin ⁡ (α) cos ⁡ (α)). {\ displaystyle A = {\ begin {pmatrix} \ cos (\ alpha) - \ sin (\ alpha) \\\ sin (\ alpha) \ cos (\ alpha) \\\ end {pmatrix}}.}A = \ begin {pmatrix} \ cos (\ alpha) - \ sin (\ alpha) \\ \ sin (\ alpha) \ cos (\ alpha) \\ \ конец {pmatrix}.

Для любого целого n матрица

B n = (α + 2 π n) (0–1 1 0), {\ displaystyle B_ {n} = (\ alpha +2 \ pi n) {\ begin {pmatrix} 0 -1 \\ 1 0 \\\ end {pmatrix}},}B_n = (\ alpha + 2 \ pi n) \ begin {pmatrix} 0 -1 \\ 1 0 \\ \ конец {pmatrix},

является логарифмом A. Таким образом, матрица A имеет бесконечно много логарифмов. Это соответствует тому факту, что угол поворота определяется только с точностью до 2π.

На языке теории Ли матрицы вращения A являются элементами группы Ли SO (2). Соответствующие логарифмы B являются элементами алгебры Ли so (2), которая состоит из всех кососимметричных матриц. Матрица

(0 1 - 1 0) {\ displaystyle {\ begin {pmatrix} 0 1 \\ - 1 0 \\\ end {pmatrix}}}\ begin {pmatrix} 0 1 \\ -1 0 \\ \ end {pmatrix}

является генератором алгебры Ли так (2).

Существование

На вопрос, имеет ли матрица логарифм, есть самый простой ответ при рассмотрении в сложных условиях. Сложная матрица имеет логарифм тогда и только тогда, когда является обратимым. Логарифм не уникален, но если матрица не имеет отрицательных вещественных собственных значений, то существует уникальный логарифм, все собственные значения которого лежат в полосе {z ∈ C | −π < Im z < π}. This logarithm is known as the principal logarithm.

Ответ больше связан с реальной обстановкой. Вещественная матрица имеет действительный логарифм тогда и только тогда, когда она обратима и каждый жорданов блок, принадлежащий отрицательному собственному значению, встречается четное число раз. Если обратимая вещественная матрица не удовлетворяет условию с жордановыми блоками, то она имеет только ненастоящие логарифмы. Это уже можно увидеть в скалярном случае: никакая ветвь логарифма не может быть действительной в -1. Существование вещественных матричных логарифмов вещественных матриц 2 × 2 рассматривается в следующем разделе.

Свойства

Если A и B обе являются положительно определенными матрицами, то

tr ⁡ log ⁡ (AB) = tr ⁡ log ⁡ (A) + тр ⁡ журнал ⁡ (В), {\ displaystyle \ operatorname {tr} {\ log {(AB)}} = \ operatorname {tr} {\ log {(A)}} + \ operatorname {tr} {\ log { (B)}},}{\ displaystyle \ operatorname {tr} {\ log {(AB)}} = \ operatorname {tr} {\ log {(A)}} + \ operatorname {tr } {\ log {(B)}},}

и если A и B коммутируют, т. Е. AB = BA, то

log ⁡ (AB) = log ⁡ (A) + log ⁡ (B). {\ displaystyle \ log {(AB)} = \ log {(A)} + \ log {(B)}. \,}{\ displaystyle \ log {(AB)} = \ log {(A)} + \ log {(B)}. \,}

Подставляя в это уравнение B = A, получаем

log ⁡ (A - 1) = - журнал ⁡ (A). {\ displaystyle \ log {(A ^ {- 1})} = - \ log {(A)}.}{\ displaystyle \ log {(A ^ {- 1})} = - \ log {(A)}.}

Аналогично, теперь для некоммутирующих A и B

log ⁡ (A + t B) = журнал ⁡ (A) + t ∫ 0 ∞ dz 1 A + z IB 1 A + z I + O (t 2). {\ displaystyle \ log {(A + tB)} = \ log {(A)} + t \ int _ {0} ^ {\ infty} \! \! \! dz ~~ {\ frac {1} {A + zI}} B {\ frac {1} {A + zI}} + O (t ^ {2}).}{\ displaystyle \ log {(A + tB)} = \ log {(A)} + t \ int _ {0} ^ {\ infty} \! \! \! dz ~~ {\ frac {1} {A + zI}} B {\ frac {1} {A + zI}} + O (t ^ {2}).}
Дополнительный пример: логарифм вращений в трехмерном пространстве

Поворот R ∈ SO (3) в ³ задается ортогональной матрицей 3 × 3 .

. Логарифм такой матрицы вращения R может быть легко вычислен из антисимметричной части формулы вращения Родригеса (см. также Угол оси ). Он дает логарифм минимальной нормы Фробениуса, но терпит неудачу, когда R имеет собственные значения, равные -1, где это не уникально.

Также обратите внимание, что для заданных матриц вращения A и B

dg (A, B): = ‖ log ⁡ (A ⊤ B) ‖ F {\ displaystyle d_ {g} (A, B) : = \ | \ log (A ^ {\ top} B) \ | _ {F}}d_g (A, B): = \ | \ log (A ^ \ top B) \ | _F

- геодезическое расстояние на трехмерном многообразии матриц вращения.

Вычисление логарифма диагонализуемой матрицы

Метод нахождения ln A для диагонализуемой матрицы A следующий:

Найти матрицу V из собственные векторы матрицы A (каждый столбец матрицы V является собственным вектором матрицы A).
Найдите обратный V числа V.
Пусть
A ′ = V - 1 AV. {\ displaystyle A '= V ^ {- 1} AV. \,} A' = V^{-1} A V.\,
Тогда A ′ будет диагональной матрицей, диагональные элементы которой являются собственными значениями A.
Заменить каждый диагональный элемент A ′ на его (натуральный) логарифм, чтобы получить log ⁡ A ′ {\ displaystyle \ log A '}{\displaystyle \log A'}.
Затем
log ⁡ A = V (log ⁡ A ′) V - 1. {\ displaystyle \ log A = V (\ log A ') V ^ {- 1}. \,}{\displaystyle \log A=V(\log A')V^{-1}.\,}

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

Логарифм недиагонализуемой матрицы

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

[1 1 0 1]. {\ displaystyle {\ begin {bmatrix} 1 1 \\ 0 1 \ end {bmatrix}}.}\ begin {bmatrix} 1 1 \\ 0 1 \ end {bmatrix}.

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

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

B = (λ 1 0 0 ⋯ 0 0 λ 1 0 ⋯ 0 0 0 λ 1 ⋯ 0 ⋮ ⋮ ⋮ ⋱ ⋱ ⋮ 0 0 0 0 λ 1 0 0 0 0 0 λ) = λ (1 λ - 1 0 0 ⋯ 0 0 1 λ - 1 0 ⋯ 0 0 0 1 λ - 1 ⋯ 0 ⋮ ⋮ ⋮ ⋱ ⋱ ⋮ 0 0 0 0 1 λ - 1 0 0 0 0 0 1) = λ (I + K) {\ displaystyle B = {\ begin {pmatrix} \ lambda 1 0 0 \ cdots 0 \\ 0 \ lambda 1 0 \ cdots 0 \\ 0 0 \ lambda 1 \ cdots 0 \\\ vdots \ vdots \ vdots \ ddots \ ddots \ vdots \\ 0 0 0 0 \ lambda 1 \\ 0 \\ 0 0 0 end { pmatrix}} = \ lambda {\ begin {pmatrix} 1 \ lambda ^ {- 1} 0 0 \ cdots 0 \\ 0 1 \ lambda ^ {- 1} 0 \ cdots 0 \\ 0 0 1 \ lambda ^ {- 1} \ cdots 0 \\\ vdots \ vdots \ vdots \ ddots \ ddots \ vdots \\ 0 0 0 0 1 \ lambda ^ {- 1} \\ 0 0 0 0 0 0 1 \\\ end {pmatrix}} = \ lambda (I + K) }B = \ begin {pmatrix} \ lambda 1 0 0 \ cdots 0 \\ 0 \ lambda 1 0 \ cdots 0 \\ 0 0 \ lambda 1 \ cdots 0 \\ \ vdots \ vdots \ vdots \ ddots \ ddots \ vdots \\ 0 0 0 0 \ lambda 1 \\ 0 0 0 0 0 0 \ lambda \\\ end {pmatrix} = \ lambda \ begin {pmatrix} 1 \ lambda ^ {- 1} 0 0 \ cdots 0 \\ 0 1 \ lambda ^ {- 1} 0 \ cdots 0 \\ 0 0 1 \ lambda ^ {- 1} \ cdots 0 \\ \ vdots \ vdots \ vdots \ ddots \ ddots \ vdots \\ 0 0 0 0 1 \ lambda ^ {- 1 } \\ 0 0 0 0 0 1 \\\ end {pmatrix} = \ lambda (I + K)

где K - матрица с нулями на главной диагонали и под ней. (Число λ отлично от нуля в предположении, что матрица, логарифм которой пытаются найти, обратима.)

Тогда по ряду Меркатора

log ⁡ (1 + x) = x - x 2 2 + x 3 3 - x 4 4 + ⋯ {\ displaystyle \ log (1 + x) = x - {\ frac {x ^ {2}} {2}} + {\ frac {x ^ {3} } {3}} - {\ frac {x ^ {4}} {4}} + \ cdots}{ \ displaystyle \ log (1 + x) = x - {\ frac {x ^ {2}} {2}} + {\ frac {x ^ {3}} {3}} - {\ frac {x ^ {4 }} {4}} + \ cdots}

получается

log ⁡ B = log ⁡ (λ (I + K)) = log ⁡ ( λ I) + журнал ⁡ (I + K) знак равно (журнал ⁡ λ) I + K - K 2 2 + K 3 3 - K 4 4 + ⋯ {\ displaystyle \ log B = \ log {\ big (} \ lambda (I + K) {\ big)} = \ log (\ lambda I) + \ log (I + K) = (\ log \ lambda) I + K - {\ frac {K ^ {2}} {2} } + {\ frac {K ^ {3}} {3}} - {\ frac {K ^ {4}} {4}} + \ cdots}{\ displaystyle \ log B = \ log {\ big (} \ lambda (I + K) {\ big)} = \ log (\ lambda I) + \ log (I + K) = ( \ log \ lambda) I + K - {\ frac {K ^ {2}} {2}} + {\ frac {K ^ {3}} {3}} - {\ frac {K ^ {4}} { 4}} + \ cdots}

Эта серия имеет конечное число членов (K равно нулю, если m - размерность K), и поэтому его сумма хорошо определена.

Используя этот подход, можно найти

log ⁡ [1 1 0 1] = [0 1 0 0]. {\ displaystyle \ log {\ begin {bmatrix} 1 1 \\ 0 1 \ end {bmatrix}} = {\ begin {bmatrix} 0 1 \\ 0 0 \ end {bmatrix}}.}{\ displaystyle \ log {\ begin {bmatrix} 1 1 \\ 0 1 \ end {bmatrix}} = {\ begin {bmatrix} 0 1 \\ 0 0 \ end {bmatrix}}.}
Перспектива функционального анализа

Квадратная матрица представляет собой линейный оператор в евклидовом пространстве R, где n - размерность матрицы. Поскольку такое пространство конечномерно, этот оператор фактически ограничен.

Используя инструменты голоморфного функционального исчисления, задана голоморфная функция f (z), определенная на открытое множество в комплексной плоскости и ограниченный линейный оператор T, можно вычислить f (T), пока f (z) определена на спектре of T.

Функция f (z) = log z может быть определена на любом односвязном открытом множестве в комплексной плоскости, не содержащем начала координат, и она голоморфна на таком домен. Это означает, что можно определить ln T до тех пор, пока спектр T не содержит начало координат и существует путь, идущий от начала координат к бесконечности, не пересекающий спектр T (например, если спектр T представляет собой круг с начало внутри него, определить ln T невозможно).

Спектр линейного оператора на R - это набор собственных значений его матрицы, а значит, и конечный набор. Пока начало координат не находится в спектре (матрица обратима), условие пути из предыдущего абзаца выполняется и ln T четко определено. Неединственность логарифма матрицы следует из того факта, что можно выбрать более одной ветви логарифма, которая определена на множестве собственных значений матрицы.

Теория групп Ли

В теории групп Ли существует экспоненциальное отображение из алгебры Ли g в соответствующую группу Ли G

exp: g → G. {\ displaystyle \ exp: g \ rightarrow G.}\ exp: g \ rightarrow G.

Для матричных групп Ли элементы g и G являются квадратными матрицами, а экспоненциальное отображение задается экспоненциальной матрицей . Обратное отображение log = exp - 1 {\ displaystyle \ log = \ exp ^ {- 1}}\ log = \ exp ^ {- 1} является многозначным и совпадает с обсуждаемым здесь матричным логарифмом. Логарифм отображает группу Ли G в алгебру Ли g. Обратите внимание, что экспоненциальное отображение является локальным диффеоморфизмом между окрестностью U нулевой матрицы 0 _ ∈ g {\ displaystyle {\ underline {0}} \ in g}\ underline {0} \ in g и окрестностью V матрицы единичная матрица 1 _ ∈ G {\ displaystyle {\ underline {1}} \ in G}\ underline {1} \ in G . Таким образом, (матричный) логарифм корректно определяется как отображение

log: V ⊂ G → U ⊂ g. {\ displaystyle \ log: V \ subset G \ rightarrow U \ subset g. \,}\ log: V \ subset G \ rightarrow U \ subset g. \,

Важное следствие формулы Якоби, тогда

log ⁡ (det (A)) = tr (журнал ⁡ A). {\ displaystyle \ log (\ det (A)) = \ mathrm {tr} (\ log A) ~.}\ log (\ det (A)) = \ mathrm {tr} (\ log A) ~.
Ограничения в случае 2 × 2

Если вещественная матрица 2 × 2 имеет отрицательный определитель, он не имеет действительного логарифма. Прежде всего отметим, что любую вещественную матрицу 2 × 2 можно рассматривать как один из трех типов комплексного числа z = x + y ε, где ε² ∈ {−1, 0, +1}. Эта z - точка на сложной подплоскости кольца матриц.

Случай, когда определитель отрицательный, возникает только в плоскости с ε² = + 1, то есть в плоскости комплексных чисел с разбиением. Только одна четверть этой плоскости является изображением экспоненциальной карты, поэтому логарифм определяется только в этой четверти (квадранте). Остальные три квадранта являются изображениями этого квадранта под четырьмя группами Клейна , порожденными ε и −1.

Например, пусть a = log 2; тогда ch a = 5/4 и sh a = 3/4. Для матриц это означает, что

A = exp ⁡ (0 aa 0) = (cosh ⁡ a sinh ⁡ a sinh ⁡ a cosh ⁡ a) = (1,25 0,75 0,75 1,25) {\ displaystyle A = \ exp { \ begin {pmatrix} 0 a \\ a 0 \ end {pmatrix}} = {\ begin {pmatrix} \ cosh a \ sinh a \\\ sinh a \ cosh a \ end {pmatrix}} = {\ begin {pmatrix} 1.25.75 \\. 75 1.25 \ end {pmatrix}}}A = \ exp \ begin {pmatrix} 0 a \\ a 0 \ end {pmatrix} = \ begin {pmatrix} \ cosh a \ sinh a \\ \ sinh a \ cosh a \ end {pmatrix} = \ begin {pmatrix} 1.25.75 \\.75 1.25 \ end {pmatrix} .

Итак, эта последняя матрица имеет логарифм

log ⁡ A = (0 log ⁡ 2 log ⁡ 2 0) {\ displaystyle \ log A = { \ begin {pmatrix} 0 \ log 2 \\\ log 2 0 \ end {pmatrix}}}{\ displaystyle \ log A = {\ begin {pmatrix} 0 \ log 2 \\\ log 2 0 \ end {pmatrix}}} .

Эти матрицы, однако, не имеют логарифма:

(3/4 5/4 5/4 3 / 4), (- 3/4 - 5/4 - 5/4 - 3/4), (- 5/4 - 3/4 - 3/4 - 5/4) {\ displaystyle {\ begin {pmatrix} 3 / 4 и 5/4 \\ 5/4 и 3/4 \ end {pmatrix}}, \ {\ begin {pmatrix} -3 / 4 -5 / 4 \\ - 5/4 -3 / 4 \ end {pmatrix}}, \ {\ begin {pmatrix} -5 / 4 -3 / 4 \\ - 3/4 -5 / 4 \ end {pmatrix}}}\ begin {pmatrix} 3/4 и 5/4 \\ 5/4 и 3/4 \ end { pmatrix}, \ \ begin {pmatrix} -3/4 -5/4 \\ -5/4 -3/4 \ end {pmatrix}, \ \ begin {pmatrix} -5/4 -3/4 \\ -3/4 -5/4 \ end {pmatrix} .

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

Неособая матрица 2 x 2 не обязательно имеет логарифм, но она сопряжена четырехгруппой с матрицей, которая имеет логарифм.

Отсюда также следует, что, например, квадратный корень из этой матрицы A можно получить непосредственно возведением в степень (logA) / 2,

A = (ch ⁡ ((log ⁡ 2) / 2) sinh ⁡ ((log ⁡ 2) / 2) sinh ⁡ ((log ⁡ 2) / 2) ch ⁡ ((log 2) / 2)) = (1,06 0,35 0,35 1,06). {\ Displaystyle {\ sqrt {A}} = {\ begin {pmatrix} \ cosh ((\ log 2) / 2) \ sinh ((\ log 2) / 2) \\\ sinh ((\ log 2) / 2) \ cosh ((\ log 2) / 2) \ end {pmatrix}} = {\ begin {pmatrix} 1.06.35 \\. 35 1.06 \ end {pmatrix}} ~.}{\ displaystyle {\ sqrt {A}} = {\ begin {pmatrix} \ cosh ((\ log 2) / 2) \ sinh ( (\ log 2) / 2) \\\ sinh ((\ log 2) / 2) \ ch ((\ log 2) / 2) \ end {pmatrix}} = {\ begin {pmatrix} 1.06 и.35 \\. 35 1.06 \ end {pmatrix}} ~.}

Для более богатого примера начните с тройки Пифагора (p, q, r) и пусть a = log (p + r) - log q. Тогда

ea = p + rq = cosh ⁡ a + sinh ⁡ a {\ displaystyle e ^ {a} = {\ frac {p + r} {q}} = \ cosh a + \ sinh a}e ^ a = \ frac {p + r} {q} = \ ch a + \ sinh a .

Теперь

ехр ⁡ (0 аа 0) = (г / qp / qp / qr / q) {\ displaystyle \ exp {\ begin {pmatrix} 0 a \\ a 0 \ end {pmatrix}} = {\ begin {pmatrix} r / q p / q \\ p / q r / q \ end {pmatrix}}}\ exp \ begin {pmatrix} 0 a \\ a 0 \ end {pmatrix} = \ begin {pmatrix} r / q p / q \\ p / q r / q \ end {pmatrix} .

Таким образом,

1 q (rppr) {\ displaystyle {\ tfrac {1} {q}} {\ begin {pmatrix} r p \\ p r \ end {pmatrix}}}\ tfrac {1} {q} \ begin {pmatrix} r p \\ p r \ end {pmatrix}

имеет матрицу логарифма

(0 aa 0) {\ displaystyle {\ begin {pmatrix} 0 a \\ a 0 \ end {pmatrix}}}\ begin {pmatrix} 0 a \\ a 0 \ end {pmatrix} ,

где а = журнал (р + г) - журнал q.

См. Также
Примечания
  1. ^Холл 2015 Теорема 2.8
  2. ^Хайэм (2008), Теорема 1.27
  3. ^Хайэм (2008), Теорема 1.31
  4. ^Калвер (1966)
  5. ^Энгё (2001)
  6. ^Холл 2015 Теорема 3.42
Литература
Последняя правка сделана 2021-05-28 05:30:58
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте