Двойной кватернион

редактировать
Мемориальная доска на мосту Брум (Дублин) в память об изобретении Гамильтоном кватернионов

В математике двойные кватернионы представляют собой 8-мерную действительную алгебру, изоморфную тензорному произведению кватернионов и двойных чисел. Таким образом, они могут быть построены так же, как кватернионы, за исключением использования двойных чисел вместо действительных чисел в качестве коэффициентов. Двойственный кватернион может быть представлен в виде A + εB, где A и B - обычные кватернионы, а ε - двойственная единица, которая удовлетворяет ε = 0 и коммутирует с каждым элементом алгебры. В отличие от кватернионов, двойные кватернионы не образуют алгебру деления.

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

Подобно тому, как вращения в трехмерном пространстве могут быть представлены кватернионами единичной длины, жесткие движения в трехмерном пространстве могут быть представлены двойными кватернионами единичной длины. Этот факт используется в теоретической кинематике (см. Маккарти), а также в приложениях к 3D компьютерной графике, робототехнике и компьютерному зрению.

Содержание
  • 1 История
  • 2 Формулы
    • 2.1 Сложение
    • 2.2 Умножение
    • 2.3 Сопряжение
    • 2.4 Сопряжение с двойным числом
    • 2,5 Норма
    • 2,6 Обратное
  • 3 Двойные кватернионы и пространственные смещения
  • 4 Матричная форма умножения двойных кватернионов
  • 5 Подробнее о пространственных смещениях
  • 6 Двойные кватернионы и однородные преобразования 4 × 4
  • 7 Связь с алгебрами Клиффорда
  • 8 Эпонимы
  • 9 См. Также
  • 10 Список литературы
  • 11 Дополнительная литература
  • 12 Внешние ссылки
История

W. Р. Гамильтон ввел кватернионы в 1843 году, а к 1873 году W. К. Клиффорд получил широкое обобщение этих чисел, которое он назвал бикватернионами, что является примером того, что сейчас называется алгеброй Клиффорда.

. В 1898 году Александр МакОлей использовал Ω с Ω = 0, чтобы сгенерировать двойственную кватернионную алгебру. Однако его терминология «октонионов» не сохранилась, поскольку сегодняшние октонионы представляют собой другую алгебру.

В России Александр Котельников разработал двойственные векторы и двойные кватернионы для использования в изучении механики.

В 1891 Эдуард Штюч понял, что эта ассоциативная алгебра идеальна для описания группы движений трехмерного пространства. Он далее развил идею в Geometrie der Dynamen в 1901 году. Б. Л. ван дер Варден назвал структуру «Изучение бикватернионов», одну из трех восьмимерных алгебр, называемых бикватернионами.

Формулы

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

Кватернион - это линейная комбинация базисных элементов 1, i, j и k. Правило произведения Гамильтона для i, j и k часто записывается как

i 2 = j 2 = k 2 = ijk = - 1. {\ displaystyle i ^ {2} = j ^ {2} = k ^ {2 } = ijk = -1.}{ \ displaystyle i ^ {2} = j ^ {2} = k ^ {2} = ijk = -1.}

Вычислить i (ijk) = −jk = −i, чтобы получить jk = i и (ijk) k = −ij = −k или ij = k. Теперь, поскольку j (j k) = j i = −k, мы видим, что это произведение дает i j = −j i, что связывает кватернионы со свойствами определителей.

Удобный способ работы с кватернионным произведением - записать кватернион как сумму скаляра и вектора, то есть A = a 0+ A, где a 0 - это действительное число и A = A 1 i + A 2 j + A 3 k - трехмерный вектор. Операции с векторной точкой и перекрестием теперь можно использовать для определения кватернионного произведения A = a 0+ Aи C = c 0+ Cкак

G = AC = (a 0 + A) (c 0 + C) = (a 0 c 0 - A ⋅ C) + (c 0 A + a 0 C + A × C). {\ Displaystyle G = AC = (a_ {0} + \ mathbf {A}) (c_ {0} + \ mathbf {C}) = (a_ {0} c_ {0} - \ mathbf {A} \ cdot \ mathbf {C}) + (c_ {0} \ mathbf {A} + a_ {0} \ mathbf {C} + \ mathbf {A} \ times \ mathbf {C}).}G = AC = (a_0 + \ mathbf {A}) (c_0 + \ mathbf { C}) = (a_0 c_0 - \ mathbf {A} \ cdot \ mathbf {C}) + (c_0 \ mathbf {A} + a_0 \ mathbf {C} + \ mathbf {A} \ times \ mathbf {C}).

Двойной кватернион обычно описывается как кватернион с двойными числами в качестве коэффициентов. Двойное число - это упорядоченная пара â = (a, b). Два двойных числа покомпонентно складываются и умножаются по правилу â ĉ = (a, b) (c, d) = (a c, a d + b c). Двойные числа часто записываются в виде â = a + εb, где ε - двойственная единица, которая коммутирует с i, j, k и обладает свойством ε = 0.

В результате двойственный кватернион может записывается как упорядоченная пара кватернионов (A, B). Два двойных кватерниона добавляются покомпонентно и умножаются по правилу

A ^ C ^ = (A, B) (C, D) = (A C, A D + B C). {\ displaystyle {\ hat {A}} {\ hat {C}} = (A, B) (C, D) = (AC, AD + BC).}{\displaystyle {\hat {A}}{\hat {C}}=(A,B)(C,D)=(AC,AD+BC).}

Двойной кватернион удобно записать как сумма двойственного скаляра и двойственного вектора, Â = â 0 + A, где â 0 = (a, b) и A = (A, B) - двойственное вектор, определяющий винт . Эти обозначения позволяют нам записать произведение двух двойственных кватернионов как

G ^ = A ^ C ^ = (a ^ 0 + A) (c ^ 0 + C) = (a ^ 0 c ^ 0 - A ⋅ C) + (c ^ 0 A + a ^ 0 C + A × C). {\ displaystyle {\ hat {G}} = {\ hat {A}} {\ hat {C}} = ({\ hat {a}} _ {0} + {\ mathsf {A}}) ({\ шляпа {c}} _ {0} + {\ mathsf {C}}) = ({\ hat {a}} _ {0} {\ hat {c}} _ {0} - {\ mathsf {A}} \ cdot {\ mathsf {C}}) + ({\ hat {c}} _ {0} {\ mathsf {A}} + {\ hat {a}} _ {0} {\ mathsf {C}} + {\ mathsf {A}} \ times {\ mathsf {C}}).}\ hat {G} = \ hat {A} \ hat {C} = (\ hat {a} _0 + \ mathsf {A}) (\ шляпа {c} _0 + \ mathsf {C}) = (\ hat {a} _0 \ hat {c} _0 - \ mathsf {A} \ cdot \ mathsf {C}) + (\ hat {c} _0 \ mathsf {A} + \ hat {a} _0 \ mathsf {C} + \ mathsf {A} \ times \ mathsf {C}).

Сложение

Сложение двойных кватернионов определяется покомпонентно, так что задано,

A ^ = (A, В) знак равно a 0 + a 1 я + a 2 j + a 3 К + б 0 ε + b 1 ε я + b 2 ε j + b 3 ε К, {\ displaystyle {\ hat {A}} = (A, B) = a_ {0} + a_ {1} i + a_ {2} j + a_ {3} k + b_ {0} \ varepsilon + b_ {1} \ varepsilon i + b_ {2} \ varepsilon j + b_ {3} \ varepsilon k,}{\ displaystyle {\ шляпа {A}} = (A, B) = a_ {0} + a_ {1} i + a_ {2} j + a_ {3} k + b_ {0} \ varepsilon + b_ {1} \ varepsilon i + б_ {2} \ varepsi lon j + b_ {3} \ varepsilon k,}

и

C ^ = (C, D) = c 0 + c 1 i + c 2 j + c 3 k + d 0 ε + d 1 ε i + d 2 ε J + d 3 ε К, {\ Displaystyle {\ Hat {C}} = (C, D) = c_ {0} + c_ {1} я + c_ {2} j + c_ {3} k + d_ {0} \ varepsilon + d_ {1} \ varepsilon i + d_ {2} \ varepsilon j + d_ {3} \ varepsilon k,}{\ displaystyle {\ hat {C}} = (C, D) = c_ {0} + c_ {1} i + c_ {2} j + c_ {3} k + d_ {0} \ varepsilon + d_ {1 } \ varepsilon i + d_ {2} \ varepsilon j + d_ {3} \ varepsilon k,}

, затем

A ^ + C ^ = (A + C, B + D) = (a 0 + c 0) + (a 1 + c 1) i + (a 2 + c 2) j + (a 3 + c 3) k + (b 0 + d 0) ε + (b 1 + d 1) ε i + (b 2 + d 2) ε j + (б 3 + d 3) ε К, {\ displaystyle {\ hat {A}} + {\ hat {C}} = (A + C, B + D) = (a_ {0} + c_ {0}) + (a_ {1} + c_ {1}) i + (a_ {2} + c_ {2}) j + (a_ {3} + c_ {3}) k + (b_ {0} + d_ {0}) \ varepsilon + (b_ {1} + d_ {1}) \ varepsilon i + (b_ {2} + d_ {2}) \ varepsilon j + (b_ {3} + d_ {3}) \ varepsilon k,}{\ displaystyle {\ hat { A}} + {\ hat {C}} = (A + C, B + D) = (a_ {0} + c_ {0}) + (a_ {1} + c_ {1}) i + (a_ {2 } + c_ {2}) j + (a_ {3} + c_ {3}) k + (b_ {0} + d_ {0}) \ varepsilon + (b_ {1} + d_ {1}) \ varepsilon i + (b_ {2} + d_ {2}) \ varepsilon j + (b_ {3} + d_ {3}) \ varepsilon k,}

Умножение

Умножение двух двойственных кватернионов следует из правил умножения для кватернионных единиц i, j, k и коммутативного умножения на двойственную единицу ε. В частности, при

A ^ = (A, B) = A + ε B, {\ displaystyle {\ hat {A}} = (A, B) = A + \ varepsilon B,}{\ displaystyle {\ hat {A}} = (A, B) = A + \ varepsilon B,}

и

C ^ = (C, D) = C + ε D, {\ displaystyle {\ hat {C}} = (C, D) = C + \ varepsilon D,}{\ displaystyle {\ hat {C}} = (C, D) = C + \ varepsilon D, }

, затем

A ^ C ^ = (A + ε B) (C + ε D) = AC + ε (AD + BC). {\ displaystyle {\ hat {A}} {\ hat {C}} = (A + \ varepsilon B) (C + \ varepsilon D) = AC + \ varepsilon (AD + BC). \!}{\ displaystyle {\ hat {A}} {\ hat {C} } = (A + \ varepsilon B) (C + \ varepsilon D) = AC + \ varepsilon (AD + BC). \!}

Обратите внимание, что есть нет термина BD, потому что определение двойных чисел требует, чтобы ε = 0.

Это дает нам таблицу умножения (обратите внимание, что порядок умножения - это строка, умноженная на столбец):

Таблица умножения для двойных кватернионов
(Строка x столбец)1ijkεεiεjεk
11ijkεεiεjεk
ii−1k−jεi−εεk−εj
jj−k−1iεj-εk−εεi
kkj−i−1εkεj−εi−ε
εεεiεjεk0000
εiεi−εεk−εj0000
εjεj−εk−εεi0000
εkεkεj−εi−ε0000

Conjugate

Сопряжение двойного кватерниона является расширение сопряженного кватерниона, то есть

A ^ ∗ = (A ∗, B ∗) = A ∗ + ε B ∗. {\ displaystyle {\ hat {A}} ^ {*} = (A ^ {*}, B ^ {*}) = A ^ {*} + \ varepsilon B ^ {*}. \!}{\ displaystyle {\ hat {A}} ^ {*} = (A ^ {*}, B ^ {*}) = A ^ {*} + \ varepsilon B ^ {*}. \!}

Как с кватернионами, сопряжение произведения двойственных кватернионов, Ĝ = ÂĈ, является произведением их сопряженных в обратном порядке,

G ^ ∗ = (A ^ C ^) ∗ = C ^ ∗ A ^ ∗. {\ displaystyle {\ hat {G}} ^ {*} = ({\ hat {A}} {\ hat {C}}) ^ {*} = {\ hat {C}} ^ {*} {\ hat {A}} ^ {*}. \!}\ hat {G} ^ * = (\ hat {A} \ hat {C}) ^ * = \ hat {C} ^ * \ hat {A} ^ *. \!

Полезно ввести функции Sc (∗) и Vec (∗), которые выбирают скалярную и векторную части кватерниона или двойственные скалярные и двойственные векторные части двойного кватерниона. В частности, если Â = â 0 + A, то

Sc (A ^) = a ^ 0, Vec (A ^) = A. {\ displaystyle {\ t_dv {Sc}} ({\ hat {A}}) = {\ hat {a}} _ {0}, {\ t_dv {Vec}} ({\ hat {A}}) = { \ mathsf {A}}. \!}\ t_dv {Sc} (\ hat {A}) = \ hat {a} _0, \ t_dv {Vec} (\ hat {A}) = \ mathsf {A}. \!

Это позволяет определить сопряжение Â как

A ^ ∗ = Sc (A ^) - Vec (A ^). {\ displaystyle {\ hat {A}} ^ {*} = {\ t_dv {Sc}} ({\ hat {A}}) - {\ t_dv {Vec}} ({\ hat {A}}). \ !}\ hat {A} ^ * = \ t_dv {Sc} (\ hat {A}) - \ t_dv {Vec} (\ hat {A}). \!

или,

(a ^ 0 + A) ∗ = a ^ 0 - A. {\ displaystyle ({\ hat {a}} _ {0} + {\ mathsf {A}}) ^ {*} = {\ hat {a}} _ {0} - {\ mathsf {A}}. \ !}(\ hat {a} _0 + \ mathsf {A}) ^ * = \ hat {a} _0 - \ mathsf {A}. \!

Произведение двойственного кватерниона на сопряженный ему дает

A ^ A ^ ∗ = (a ^ 0 + A) (a ^ 0 - A) = a ^ 0 2 + A ⋅ A. {\ displaystyle {\ hat {A}} {\ hat {A}} ^ {*} = ({\ hat {a}} _ {0} + {\ mathsf {A}}) ({\ hat {a}) } _ {0} - {\ mathsf {A}}) = {\ hat {a}} _ {0} ^ {2} + {\ mathsf {A}} \ cdot {\ mathsf {A}}. \! }\ hat {A} \ hat {A} ^ * = (\ hat {a} _0 + \ mathsf {A}) (\ hat {a} _0 - \ mathsf {A}) = \ hat {a} _0 ^ 2 + \ mathsf {A} \ cdot \ mathsf {A}. \!

Это двойной скаляр, который представляет собой квадрат величины двойного кватерниона.

Сопряжение двойного числа

Второй тип сопряжения двойственного кватерниона дается взятием сопряженного двойственного числа, заданного как

A ^ ¯ = (A, - B) = A - ε B. {\ displaystyle {\ overline {\ hat {A}}} = (A, -B) = A- \ varepsilon B. \!}{\ displaystyle {\ overline {\ hat {A}}} = (A, -B) = A- \ varepsilon B. \!}

Кватернионное и двойное числовое сопряжение можно объединить в третью форму сопряженного данного по

A ^ ∗ ¯ = (A ∗, - B ∗) = A ∗ - ε B ∗. {\ displaystyle {\ overline {{\ hat {A}} ^ {*}}} = (A ^ {*}, - B ^ {*}) = A ^ {*} - \ varepsilon B ^ {*}. \!}{\ displaystyle {\ overline {{\ hat {A}} ^ {*}}} = (A ^ {*}, - B ^ {*}) = A ^ {*} - \ варепсилон B ^ {*}. \!}

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

Норма

Норма двойного кватерниона |  | вычисляется с использованием конъюгата для вычисления |  | = √ Â. Это двойное число, называемое величиной двойного кватерниона. Двойные кватернионы с |  | = 1 - единичные двойственные кватернионы.

Двойные кватернионы величины 1 используются для представления пространственных евклидовых смещений. Обратите внимание, что требование, что Â Â = 1, вводит два алгебраических ограничения на компоненты Â, то есть

A ^ A ^ ∗ = (A, B) (A ∗, B ∗) = (AA ∗, AB ∗ + BA ∗) = (1, 0). {\ displaystyle {\ hat {A}} {\ hat {A}} ^ {*} = (A, B) (A ^ {*}, B ^ {*}) = (AA ^ {*}, AB ^ {*} + BA ^ {*}) = (1,0). \!}\ hat {A} \ шляпа {A} ^ * = (A, B) (A ^ *, B ^ *) = (AA ^ *, AB ^ * + BA ^ *) = (1, 0). \!

Обратный

Если p + ε q - двойственный кватернион, а p не равен нулю, то обратный двойственный кватернион задается как

p (1 - ε qp).

Таким образом, элементы подпространства {ε q: q ∈ H} не имеют обратных. Это подпространство называется идеалом в теории колец. Это единственный максимальный идеал кольца двойственных чисел.

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

Двойные кватернионы и пространственные смещения

Преимущество формулировки двойных кватернионов композиции двух пространственных смещений D B = ([R B ], b ) и D A = ([R A],a) состоит в том, что полученный двойной кватернион дает непосредственно винт ось и двойной угол составного смещения D C = D BDA.

В общем, двойной кватернион, связанный с пространственным смещением D = ([A], d ) состоит из его оси винта S = (S, V) и двойного угла (φ, d), где φ - вращение вокруг, а d - скольжение вдоль этой оси, что определяет смещение D. Соответствующий двойственный кватернион определяется как,

S ^ = cos ⁡ ϕ ^ 2 + sin ⁡ ϕ ^ 2 S. {\ displaystyle {\ hat {S}} = \ cos {\ frac {\ hat {\ phi}} {2}} + \ sin {\ frac {\ hat {\ phi}} {2}} {\ mathsf { S}}.}\ hat {S} = \ cos \ frac {\ hat {\ phi}} {2} + \ sin \ frac {\ hat {\ phi}} {2} \ mathsf {S}.

Пусть композиция смещения D B с D A будет смещением D C = D BDA. Ось винта и двойной угол D C получают из произведения двойных кватернионов D A и D B, определяемых как

A ^ = cos ⁡ (α ^ / 2) + sin ⁡ (α ^ / 2) A и B ^ = cos ⁡ (β ^ / 2) + sin ⁡ (β ^ / 2) B. {\ displaystyle {\ hat {A}} = \ cos ({\ hat {\ alpha}} / 2) + \ sin ({\ hat {\ alpha}} / 2) {\ mathsf {A}} \ quad { \ text {and}} \ quad {\ hat {B}} = \ cos ({\ hat {\ beta}} / 2) + \ sin ({\ hat {\ beta}} / 2) {\ mathsf {B }}.}\ hat {A} = \ cos (\ hat {\ alpha} / 2) + \ sin (\ hat {\ alpha} / 2) \ mathsf {A} \ quad \ text {and} \ quad \ hat {B} = \ cos (\ hat {\ beta} / 2) + \ sin (\ hat {\ beta} / 2) \ mathsf {B}.

То есть составное смещение D C=DBDAимеет связанный двойственный кватернион, задаваемый формулой

C ^ = cos ⁡ γ ^ 2 + sin ⁡ γ ^ 2 C = (cos ⁡ β ^ 2 + sin ⁡ β ^ 2 B) (cos ⁡ α ^ 2 + sin ⁡ α ^ 2 A). {\ displaystyle {\ hat {C}} = \ cos {\ frac {\ hat {\ gamma}} {2}} + \ sin {\ frac {\ hat {\ gamma}} {2}} {\ mathsf { C}} = {\ Big (} \ cos {\ frac {\ hat {\ beta}} {2}} + \ sin {\ frac {\ hat {\ beta}} {2}} {\ mathsf {B} } {\ Big)} {\ Big (} \ cos {\ frac {\ hat {\ alpha}} {2}} + \ sin {\ frac {\ hat {\ alpha}} {2}} {\ mathsf { A}} {\ Big)}.}\ hat {C} = \ cos \ frac {\ hat {\ gamma}} {2} + \ sin \ frac {\ hat {\ gamma}} {2} \ mathsf {C} = \ Big (\ cos \ frac {\ hat {\ beta}} {2} + \ sin \ frac {\ hat {\ beta}} {2} \ mathsf {B} \ Big) \ Big (\ cos \ frac {\ hat {\ alpha}} {2} + \ sin \ frac {\ hat {\ alpha}} {2} \ mathsf {A} \ Big).

Разверните этот продукт, чтобы получить

cos ⁡ γ ^ 2 + sin ⁡ γ ^ 2 C = (cos ⁡ β ^ 2 cos ⁡ α ^ 2 - sin ⁡ β ^ 2 sin ⁡ α ^ 2 B ⋅ A) + (sin ⁡ β ^ 2 cos ⁡ α ^ 2 B + sin ⁡ α ^ 2 cos ⁡ β ^ 2 A + sin ⁡ β ^ 2 sin ⁡ α ^ 2 B × А). {\ displaystyle \ cos {\ frac {\ hat {\ gamma}} {2}} + \ sin {\ frac {\ hat {\ gamma}} {2}} {\ mathsf {C}} = {\ Big ( } \ cos {\ frac {\ hat {\ beta}} {2}} \ cos {\ frac {\ hat {\ alpha}} {2}} - \ sin {\ frac {\ hat {\ beta}} { 2}} \ sin {\ frac {\ hat {\ alpha}} {2}} {\ mathsf {B}} \ cdot {\ mathsf {A}} {\ Big)} + {\ Big (} \ sin { \ frac {\ hat {\ beta}} {2}} \ cos {\ frac {\ hat {\ alpha}} {2}} {\ mathsf {B}} + \ sin {\ frac {\ hat {\ alpha }} {2}} \ cos {\ frac {\ hat {\ beta}} {2}} {\ mathsf {A}} + \ sin {\ frac {\ hat {\ beta}} {2}} \ sin {\ frac {\ hat {\ alpha}} {2}} {\ mathsf {B}} \ times {\ mathsf {A}} {\ Big)}.}\ cos \ frac {\ hat {\ gamma} } {2} + \ sin \ frac {\ hat {\ gamma}} {2} \ mathsf {C} = \ Big (\ cos \ frac {\ hat {\ beta}} {2} \ cos \ frac {\ шляпа {\ alpha}} {2} - \ sin \ frac {\ hat {\ beta}} {2} \ sin \ frac {\ hat {\ alpha}} {2} \ mathsf {B} \ cdot \ mathsf { A} \ Big) + \ Big (\ sin \ frac {\ hat {\ beta}} {2} \ cos \ frac {\ hat {\ alpha}} {2} \ mathsf {B} + \ sin \ frac { \ hat {\ alpha}} {2} \ cos \ frac {\ hat {\ beta}} {2} \ mathsf {A} + \ sin \ frac {\ hat {\ beta}} {2} \ sin \ frac {\ hat {\ alpha}} {2} \ mathsf {B} \ times \ mathsf {A} \ Big).

Разделите обе части этого уравнения на тождество

соз ⁡ γ ^ 2 = соз ⁡ β ^ 2 соз ⁡ α ^ 2 - грех ⁡ β ^ 2 грех ⁡ α ^ 2 В ⋅ A {\ displaystyle \ cos {\ frac {\ hat {\ gamma}} { 2}} = \ cos {\ frac {\ hat {\ beta}} {2}} \ cos {\ frac {\ hat {\ alpha}} {2}} - \ sin {\ frac {\ hat {\ beta }} {2}} \ sin {\ frac {\ hat {\ alpha}} {2}} {\ mathsf {B}} \ cdot {\ mathsf {A}}}\ cos \ frac {\ hat {\ gamma}} {2} = \ cos \ frac {\ hat {\ beta}} {2} \ cos \ frac {\ hat {\ alpha}} {2} - \ sin \ frac {\ hat {\ beta}} {2} \ sin \ frac {\ hat {\ alpha}} {2} \ mathsf {B} \ cdot \ mathsf {A}

, чтобы получить

tan ⁡ γ ^ 2 C = tan ⁡ β ^ 2 B + tan ⁡ α ^ 2 A + tan ⁡ β ^ 2 tan ⁡ α ^ 2 B × A 1 - tan ⁡ β ^ 2 tan ⁡ α ^ 2 Б А. {\ displaystyle \ tan {\ frac {\ hat {\ gamma}} {2}} {\ mathsf {C}} = {\ frac {\ tan {\ frac {\ hat {\ beta}} {2}} { \ mathsf {B}} + \ tan {\ frac {\ hat {\ alpha}} {2}} {\ mathsf {A}} + \ tan {\ frac {\ hat {\ beta}} {2}} \ загар {\ frac {\ hat {\ alpha}} {2}} {\ mathsf {B}} \ times {\ mathsf {A}}} {1- \ tan {\ frac {\ hat {\ beta}} { 2}} \ tan {\ frac {\ hat {\ alpha}} {2}} {\ mathsf {B}} \ cdot {\ mathsf {A}}}}.}\ tan \ frac {\ hat {\ gamma}} {2} \ mathsf {C} = \ frac {\ tan \ frac {\ hat {\ beta}} {2} \ mathsf {B} + \ tan \ frac {\ hat {\ alpha}} {2} \ mathsf {A} + \ tan \ frac {\ hat {\ beta}} {2} \ tan \ frac {\ hat {\ alpha}} { 2} \ mathsf {B} \ times \ mathsf {A}} {1 - \ tan \ frac {\ hat {\ beta}} {2} \ tan \ frac {\ hat {\ alpha}} {2} \ mathsf {B} \ cdot \ mathsf {A}}.

Это Родригес 'формула для оси винта составного смещения, определенного в терминах осей винта двух смещений. Он вывел эту формулу в 1840 году.

Три оси винта A, B и C образуют пространственный треугольник и двойные углы в этих вершинах между общими нормалями, которые образуют стороны этого треугольник напрямую связаны с двойными углами трех пространственных перемещений.

Матричная форма двойного умножения кватернионов

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

Кватернионное произведение AC - это линейное преобразование оператором A компонентов кватерниона C, поэтому существует матричное представление A, работающее с вектором, сформированным из компонентов C.

Соберите компоненты кватерниона C = c 0+ Cв массив C = (C 1, C 2, C 3, c 0). Обратите внимание, что компоненты векторной части кватерниона перечислены первыми, а скаляр - последним. Это произвольный выбор, но как только это соглашение будет выбрано, мы должны его соблюдать.

Кватернионный продукт AC теперь может быть представлен как матричный продукт

AC = [A +] C = [a 0 A 3 - A 2 A 1 - A 3 a 0 A 1 A 2 A 2 - A 1 a 0 A 3 - A 1 - A 2 - A 3 a 0] {C 1 C 2 C 3 c 0}. {\ displaystyle AC = [A ^ {+}] C = {\ begin {bmatrix} a_ {0} A_ {3} - A_ {2} A_ {1} \\ - A_ {3} a_ {0} A_ {1} A_ {2} \\ A_ {2} - A_ {1} a_ {0} A_ {3} \\ - A_ {1} - A_ {2} - A_ {3} a_ {0} \ end {bmatrix}} {\ begin {Bmatrix} C_ {1} \\ C_ {2} \\ C_ {3} \\ c_ {0} \ end {Bmatrix}}.}{\ displaystyle AC = [A ^ {+}] C = {\ begin {bmatrix} a_ {0} A_ {3} - A_ {2} A_ {1} \\ - A_ {3} a_ { 0} A_ {1} A_ {2} \\ A_ {2} - A_ {1} a_ {0} A_ {3} \\ - A_ {1} - A_ {2} - A_ {3} a_ {0} \ end {bmatrix}} {\ begin {Bmatrix} C_ {1} \\ C_ {2} \ C_ {3} \\ c_ {0} \ end {Bmatrix}}.}

Продукт AC также может можно рассматривать как операцию C над компонентами A, и в этом случае мы имеем

AC = [C -] A = [c 0 C 3 - C 2 C 1 - C 3 c 0 C 1 C 2 C 2 - C 1 c 0 C 3 - C 1 - C 2 - C 3 c 0] {A 1 A 2 A 3 a 0}. {\ displaystyle AC = [C ^ {-}] A = {\ begin {bmatrix} c_ {0} C_ {3} - C_ {2} C_ {1} \\ - C_ {3} c_ {0} C_ {1} C_ {2} \\ C_ {2} - C_ {1} c_ {0} C_ {3} \\ - C_ {1} - C_ {2} - C_ {3} c_ {0} \ end {bmatrix}} {\ begin {Bmatrix} A_ {1} \\ A_ {2} \\ A_ {3} \\ a_ {0} \ end {Bmatrix}}.}AC = [C ^ -] A = \ begin {bmatrix} c_0 C_3 -C_2 C_1 \\ -C_3 c_0 C_1 C_2 \\ C_2 -C_1 c_0 C_3 \\ -C_1 -C_2 -C_3 c_0 \ end {bmatrix} \ begin {Bmatrix} A_1 \\ A_2 \\ A_3 \\ a_0 \ end {Bmatrix}.

Двойное кватернионное произведение ÂĈ = (A, B) (C, D) = (AC, AD + BC) можно сформулировать как матричную операцию следующим образом. Соберите компоненты Ĉ в восьмеричный массив Ĉ = (C 1, C 2, C 3, c 0, D 1, D 2, D 3, d 0), тогда ÂĈ задается матричным произведением 8x8

A ^ C ^ = [A ^ +] C ^ = [A + 0 B + A +] {CD}. {\ displaystyle {\ hat {A}} {\ hat {C}} = [{\ hat {A}} ^ {+}] {\ hat {C}} = {\ begin {bmatrix} A ^ {+} 0 \\ B ^ {+} A ^ {+} \ end {bmatrix}} {\ begin {Bmatrix} C \\ D \ end {Bmatrix}}.}\ hat {A} \ hat {C} = [\ hat {A} ^ +] \ hat {C} = \ begin {bmatrix} A ^ + 0 \\ B ^ + A ^ + \ end {bmatrix} \ begin {Bmatrix} C \\ D \ end {Bmatrix}.

Как мы видели для кватернионов, произведение ÂĈ может можно рассматривать как операцию Ĉ над координатным вектором Â, что означает, что ÂĈ также может быть сформулировано как,

A ^ C ^ = [C ^ -] A ^ = [C - 0 D - C -] {AB}. {\ displaystyle {\ hat {A}} {\ hat {C}} = [{\ hat {C}} ^ {-}] {\ hat {A}} = {\ begin {bmatrix} C ^ {-} 0 \\ D ^ {-} C ^ {-} \ end {bmatrix}} {\ begin {Bmatrix} A \\ B \ end {Bmatrix}}.}\ hat {A} \ hat {C} = [\ hat {C} ^ -] \ hat {A} = \ begin {bmatrix} C ^ - 0 \\ D ^ - C ^ - \ end {bmatrix} \ begin {Bmatrix} A \\ B \ end {Bmatrix}.
Подробнее о пространственных смещениях

Двойной кватернион смещения D = ([A], d ) может быть построен из кватерниона S = cos (φ / 2) + sin (φ / 2) S, который определяет вращение [A] и векторный кватернион, построенный из вектора перемещения d, заданного как D = d 1 i + d 2 j + d 3 к. Используя это обозначение, двойственный кватернион для смещения D = ([A], d ) задается как

S ^ = S + ε 1 2 D S. {\ displaystyle {\ hat {S}} = S + \ varepsilon {\ frac {1} {2}} DS.}\ hat {S} = S + \ varepsilon \ frac {1} {2} DS.

Пусть координаты Плюккера линии в направлении x через точку p в движущемся теле и его координаты в неподвижной системе отсчета, которая находится в направлении X через точку P, будут заданы как,

x ^ = x + ε p × x и X ^ = X + ε P × X. {\ displaystyle {\ hat {x}} = \ mathbf {x} + \ varepsilon \ mathbf {p} \ times \ mathbf {x} \ quad {\ text {and}} \ quad {\ hat {X}} = \ mathbf {X} + \ varepsilon \ mathbf {P} \ times \ mathbf {X}.}\ hat {x} = \ mathbf {x} + \ varepsilon \ mathbf {p} \ times \ mathbf {x} \ quad \ text {and} \ quad \ hat {X} = \ mathbf {X} + \ varepsilon \ mathbf {P} \ times \ mathbf {X}.

Затем двойной кватернион смещения этого тела преобразует координаты Плюккера в подвижной системе отсчета в координаты Плюккера в фиксированной системе отсчета формула

X ^ = S ^ x ^ S ^ ∗ ¯. {\ displaystyle {\ hat {X}} = {\ hat {S}} {\ hat {x}} {\ overline {{\ hat {S}} ^ {*}}}.}{\ displaystyle {\ hat {X}} = {\ hat {S}} {\ hat {x}} {\ overline {{\ hat {S}} ^ {*}}}.}

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

X ^ = [S ^ +] [S ^ -] ∗ x ^. {\ displaystyle {\ hat {X}} = [{\ hat {S}} ^ {+}] [{\ hat {S}} ^ {-}] ^ {*} {\ hat {x}}.}\ hat {X} = [\ hat {S} ^ +] [\ hat {S} ^ -] ^ * \ hat {x}.

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

Двойные кватернионы и однородные преобразования 4 × 4

Может быть полезно, особенно при движении твердого тела, представить двойные кватернионы как однородные матрицы. Как указано выше, двойной кватернион можно записать как: q ^ = r + d ε {\ displaystyle {\ hat {q}} = r + d \ varepsilon}\ hat q = r + d \ varepsilon где r и d оба кватернионы. Кватернион r известен как действительная или вращательная часть, а кватернион d {\ displaystyle d}d известен как двойная часть или часть смещения.

Часть вращения может быть задана как

r = rw + rxi + ryj + rzk = cos ⁡ (θ 2) + sin ⁡ (θ 2) (a → ⋅ (i, j, k)) {\ displaystyle r = r_ {w} + r_ {x} i + r_ {y} j + r_ {z} k = \ cos \ left ({\ frac {\ theta} {2}} \ right) + \ sin \ left ({\ frac {\ theta} {2}} \ right) \ left ({\ vec {a}} \ cdot (i, j, k) \ right)}{\ displaystyle r = r_ {w} + r_ {x} i + r_ {y} j + r_ {z} k = \ cos \ left ({\ frac {\ theta} {2}} \ right) + \ sin \ left ({\ frac {\ theta} {2}} \ right) \ left ({\ vec {a}} \ cdot (i, j, k) \ right)}

где θ { \ displaystyle \ theta}\ theta - угол поворота относительно направления, заданного единичным вектором a → {\ displaystyle {\ vec {a}}}{\ vec {a}} . Смещение может быть записано как

d = 0 + Δ x 2 i + Δ y 2 j + Δ z 2 k {\ displaystyle d = 0 + {\ frac {\ Delta x} {2}} i + {\ frac {\ Delta y} {2}} j + {\ frac {\ Delta z} {2}} k}{\ displaystyle d = 0 + {\ frac {\ Delta x} {2}} i + {\ frac {\ Delta y} {2}} j + {\ frac {\ Delta z} {2}} k} .

Двойной кватернионный эквивалент 3D-вектора:

v ^: = 1 + ε ( vxi + vyj + vzk) {\ displaystyle {\ hat {v}}: = 1+ \ varepsilon (v_ {x} i + v_ {y} j + v_ {z} k)}{\ displaystyle {\ hat {v}}: = 1+ \ varepsilon (v_ {x} i + v_ {y} j + v_ {z} k)}

и его преобразование с помощью q ^ {\ displaystyle {\ hat {q}}}\ hat q задается выражением

v ^ ′ = q ^ ⋅ v ^ ⋅ q ^ ∗ ¯ {\ displaystyle {\ hat {v} } '= {\ hat {q}} \ cdot {\ hat {v}} \ cdot {\ overline {{\ hat {q}} ^ {*}}}}{\displaystyle {\hat {v}}'={\hat {q}}\cdot {\hat {v}}\cdot {\overline {{\hat {q}}^{*}}}}.

Эти двойные кватернионы (или фактически их преобразования на 3D-векторах) может быть представлена ​​однородной матрицей преобразования

T = (1 0 0 0 Δ x Δ y R Δ z) {\ displaystyle T = {\ begin {pmatrix} 1 0 0 0 \\\ Delta x \\ \ Delta y R \\\ Delta z \\\ end {pmatrix}}}{\ displaystyle T = {\ begin {pmatrix} 1 0 0 0 \\\ Delta x \\\ Delta y R \\\ Delta z \\\ end {pmatrix}}}

где ортогональная матрица 3 × 3 задается как

R = (rw 2 + rx 2 - ry 2 - rz 2 2 rxry - 2 rwrz 2 rxrz + 2 rwry 2 rxry + 2 rwrzrw 2 - rx 2 + ry 2 - р з 2 2 р й р з - 2 р ш р х 2 р х р з - 2 р ш р у 2 р й р з + 2 р ш р х р ш 2 - г х 2 - р у 2 + г г 2). {\ displaystyle R = {\ begin {pmatrix} r_ {w} ^ {2} + r_ {x} ^ {2} -r_ {y} ^ {2} -r_ {z} ^ {2} 2r_ {x} r_ {y} -2r_ {w} r_ {z} 2r_ {x} r_ {z} + 2r_ {w} r_ {y} \\ 2r_ {x} r_ {y} + 2r_ {w} r_ {z} r_ {w} ^ {2} -r_ {x} ^ {2} + r_ {y} ^ {2} -r_ {z} ^ {2} 2r_ {y} r_ {z} -2r_ {w} r_ {x } \\ 2r_ {x} r_ {z} -2r_ {w} r_ {y} 2r_ {y} r_ {z} + 2r_ {w} r_ {x} r_ {w} ^ {2} -r_ {x} ^ {2} -r_ {y} ^ {2} + r_ {z} ^ {2} \\\ end {pmatrix}}.}R =\begin{pmatrix} r_w^2+r_x^2-r_y^2-r_z^2 2r_xr_y-2r_wr_z 2r_xr_z+2r_wr_y \\ 2r_xr_y+2r_wr_z r_w^2-r_x^2+r_y^2-r_z^2 2r_yr_z-2r_wr_x \\ 2r_xr_z-2r_wr_y 2r_yr_z+2r_wr_x r_w^2-r_x^2-r_y^2+r_z^2\\ \end{pmatrix}.

Для 3D-вектора

v = (1 vxvyvz) { \ displaystyle v = {\ begin {pmatrix} 1 \\ v_ {x} \\ v_ {y} \\ v_ {z} \\\ end {pmatrix}}}{\ displaystyle v = {\ begin {pmatrix} 1 \\ v_ {x} \\ v_ {y } \\ v_ {z} \\\ конец {pmatrix}}}

преобразование по T задается как

v → ′ = T ⋅ v → {\ displaystyle {\ vec {v}} '= T \ cdot {\ vec {v}}}{\displaystyle {\vec {v}}'=T\cdot {\vec {v}}}
Связь с алгебрами Клиффорда

Помимо тензорного произведения двух алгебр Клиффорда, кватернионов и двойственных чисел, двойственные кватернионы имеют две другие формулировки в терминах алгебр Клиффорда.

Во-первых, двойственные кватернионы изоморфны алгебре Клиффорда, порожденной тремя антикоммутирующими элементами i, j, e с i = j = -1 и e = 0. Если мы определим k = ij и ε = k, то из них подразумеваются соотношения, определяющие двойственные кватернионы, и наоборот. Во-вторых, двойственные кватернионы изоморфны четной части алгебры Клиффорда, порожденной 4 антикоммутирующими элементами e 1, e 2, e 3, e 4 {\ displaystyle e_ {1}, e_ {2}, e_ {3 }, e_ {4}}{\ displaystyle e_ {1}, e_ {2}, e_ {3}, e_ {4}} с

e 1 2 = e 2 2 = e 3 2 = - 1, e 4 2 = 0. {\ displaystyle e_ {1} ^ {2} = e_ {2} ^ {2} = e_ {3} ^ {2} = - 1, \, \, e_ {4} ^ {2} = 0.}{\ displaystyle e_ {1 } ^ {2} = e_ {2} ^ {2} = e_ {3} ^ {2} = - 1, \, \, e_ {4} ^ {2} = 0.}

Подробнее см. Алгебры Клиффорда: двойные кватернионы.

Эпонимы

Поскольку и Эдуард Этюд, и Уильям Кингдон Клиффорд использовали и писали о двойных кватернионах, иногда авторы называют двойные кватернионы «Исследовать бикватернионы. "или" Клиффорд бикватернионы ". Последний эпоним также использовался для обозначения сплит-бикватернионов. Прочтите статью Джо Руни, ссылка на которую приведена ниже, чтобы увидеть сторонника W.K. Утверждение Клиффорда. Поскольку утверждения Клиффорда и Стюда являются спорными, во избежание конфликта удобно использовать текущее обозначение двойного кватерниона.

См. Также
Ссылки

Примечания

Источники

Дополнительная литература
Внешние ссылки
Последняя правка сделана 2021-05-18 05:08:27
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте