2022-10-01
Свёртка глобального вектора с помощью кватернионов
Ещё однин красивый вариант свёртки глобального вектора получается, если использовать для этого кватернионы [1].
В этом случае уже необязательно, чтобы точка двигалась прямолинейно и равномерно, как это было необходимо ранее.
Напомним, что кватернион состоит из действительной и мнимой части также, как и комплексное число, но мнимая часть состоит из трёх мнимых значений:
\[q = a + \mathbf{i} b + \mathbf{j} c + \mathbf{k} d \tag{2.1}\]
где: \(a, b, c, d\) — вещественные числа, \(\mathbf{i}, \mathbf{j}, \mathbf{k}\) — мнимые единицы, которые задаются по следующему правилу
\[\mathbf{i}^2 = \mathbf{j}^2 = \mathbf{k}^2 = \mathbf{i} \mathbf{j} \mathbf{k} = -1\]
\[\mathbf{i} \mathbf{j} = -\mathbf{j} \mathbf{i} = \mathbf{k} ,\quad \mathbf{j} \mathbf{k} = -\mathbf{k} \mathbf{j} = \mathbf{i} ,\quad \mathbf{k} \mathbf{i} = -\mathbf{i} \mathbf{k} = \mathbf{j} \tag{2.2}\]
Хотя здесь мнимые части и называются векторами, но с точки зрения векторной алгебры это не совсем векторы.
Поэтому мы будем их называть мнимыми векторами.
Тогда глобальный многомерный вектор скорости (GVV), в общем случае, может быть свёрнут до четырёхмерного пространства таким образом:
\[\mathbf{V} = \frac{c}{\gamma} \sum \limits_{n=0}^{\infty} \mathbf{j_n} \beta^n \to c\, e^{\mathbf{i}a} e^{\mathbf{j}b} e^{\mathbf{k}d} \tag{2.3}\]
Как и в случае c многомерным вектором, модуль мнимого вектора также равен скорости света:
\[ c\, | e^{\mathbf{i}a} e^{\mathbf{j}b} e^{\mathbf{k}d} | = c \tag{2.4}\]
Давайте разберём частные случаи.
Двумерное пространство
Этот случай возникает, если нам необходимо отразить (свернуть) GVV на двумерное пространство, на две его координаты.
Тогда GVV можно представить так:
\[\mathbf{V} = c\,e^{\mathbf{i}a} \tag{2.5}\]
где угол \(a\) находится через арксинус \(\beta\):
\[a = \arcsin(\beta), \quad \beta= {v \over c} \tag{2.6}\]
Причём в более общем случае все эти величины могут зависеть и меняться от времени \(t\)
\[a(t) = \arcsin(\beta(t)), \quad \beta(t)= {v(t) \over c} \tag{2.7}\]
но далее мы всё же будем применять более простой вариант формулы (2.6), подразумевая (2.7).
Тогда, раскрывая полученное выражение, и вспоминая, что \(\cos(\arcsin(\beta)) = \sqrt{1 - \beta^2}\), мы получим:
\[\mathbf{V} = c \left( \sqrt{1 - \beta^2} + \mathbf{i}\beta \right) \tag{2.8}\]
А теперь мы просим наших читателей обратить внимание на всю красоту момента, в котором мы сможем переписать эту же формулу, но не для мнимых векторов, а для настоящих:
\[\mathbf{V} = c \left( \mathbf{j_0} \sqrt{1 - \beta^2} + \mathbf{j_1} \beta \right) \tag{2.9}\]
Здесь \(\mathbf{j_0}, \mathbf{j_1}\) — единичные векторы.
Мы можем проверить правильность такого перехода так:
\[ \left| \mathbf{j}_0 \sqrt{1 - \beta^2} + \mathbf{j}_1 \beta \right| = \left| \sqrt{1 - \beta^2} + \mathbf{i}\beta \right| = 1 \tag{2.10}\]
Это же правило перехода от мнимых координат к векторным можно распространить и для бо́льших мерностей.
Четырёхмерное пространство
Здесь мы можем воспользоваться те же принципом и разложить формулу (2.3) сначала на кватернионы, а затем — перейти к векторам.
Четырёхмерный кватернион GVV выглядит так:
\[\mathbf{V} = c\, e^{\mathbf{i}a} e^{\mathbf{j}b} e^{\mathbf{k}d}, \quad a = \arcsin(\beta_1),\, b = \arcsin(\beta_2),\, d = \arcsin(\beta_3) \tag{2.11}\]
Здесь
\[ \beta_1= \beta_1(t) = {v_1(t) \over c},\, \beta_2= \beta_2(t) = {v_2(t) \over c},\, \beta_3= \beta_3(t) = {v_3(t) \over c} \tag{2.12}\]
Раскрывая формулу (2.11) мы получим такой кватернион:
\[\mathbf{V} = c\, (A + \mathbf{i} B + \mathbf{j} C + \mathbf{k} D) \tag{2.13}\]
в котором
\[A = \sqrt{1 - \beta_1^2} \sqrt{1 - \beta_2^2} \sqrt{1 - \beta_3^2} - \beta_1 \beta_2 \beta_3 \]
\[B = \beta_1 \sqrt{1 - \beta_2^2} \sqrt{1 - \beta_3^2} + \beta_2 \beta_3 \sqrt{1 - \beta_1^2} \]
\[C = \beta_2 \sqrt{1 - \beta_1^2} \sqrt{1 - \beta_3^2} - \beta_1 \beta_3 \sqrt{1 - \beta_2^2} \]
\[D = \beta_3 \sqrt{1 - \beta_1^2} \sqrt{1 - \beta_2^2} + \beta_1 \beta_2 \sqrt{1 - \beta_3^2} \tag{2.14}\]
Переходя к векторам мы получим:
\[\mathbf{V} = c\, (\mathbf{j_0} A + \mathbf{j_1} B + \mathbf{j_2} C + \mathbf{j_3} D) \tag{2.15}\]
Здесь \(\mathbf{j_0}, \mathbf{j_1}, \mathbf{j_2}, \mathbf{j_3}\) — единичные векторы.
Проверка на модуль от функции позволяет сделать вывод о правильности наших выкладок:
\[ |\mathbf{j_0} A + \mathbf{j_1} B + \mathbf{j_2} C + \mathbf{j_3} D| = |A + \mathbf{i} B + \mathbf{j} C + \mathbf{k} D| = 1 \tag{2.16}\]
Примеры с картинками
Для примера свёртки многомеронго единичного вектора в четырёхмерный кватернион, мы возьмём вращение точки по окружности вокруг оси,
плюс дополнительное движение с постоянной скоростью вдоль этой оси.
Математически отразить это можно так:
\[ \beta_1 = {v \over c} \sin(\omega t),\, \beta_2 = {v \over c} \cos(\omega t),\, \beta_3 = {v_3 \over c} \tag{2.17}\]
где \(v\) — скорость движения точки по окружности, \(v_3\) — скорость движения точки вдоль оси, \(\omega\) — частота вращения, радиан в секунду, \(t\) — время.
Свёрнутый вектор считается по формуле (2.14-2.15), куда мы и подставляем (2.17).
На картинках ниже будут приведены примеры реализации полученной формулы.
При этом нужно учитывать, что координата X — это направление \(\beta_1\), Y — это направление \(\beta_2\), а Z — это направление \(\beta_3\) (движение вдоль оси).
Время на рисунках проявляется параметрическим образом.
![]() Рис.9. Направление глобальной скорости из (2.17) при: \(v/c = 0.1,\, v_3/c = 0.1\)
|
![]() Рис.10. Другой ракурс графика (рис.9)
|
Если теперь проинтегрировать каждое направление из (2.14), которые из себя представляют проекции скорости на координаты, то мы получим график движение точки в пространстве.
Что и будут отражать следующие графики.
![]() Рис.11. Движение точки в пространстве координат XYZ при \(v/c = 0.1,\, v_3/c = 0.1\)
|
![]() Рис.12. Другой ракурс графика (рис.11). Два оборота точки
|
На рисунках 9 и 10 мы получили направление скорости точки, а на рисунках 11 и 12 — график движения точки в пространстве, в координатах XYZ.
Было правильно отобразить и координату времени, которая в формуле (2.15) представлена действительной частью \(A\), но это невозможно сделать на плоскости.
На следующих рисунках будут представлены направление скоростей и движение точки (их интегралы) при световых сокростях по одной из координат.
![]() Рис.13. Направление глобальной скорости из (2.17) при: \(v/c = 1,\, v_3/c = 0.1\)
|
![]() Рис.14. Другой ракурс графика (рис.13)
|
Интегрируем вектор глобальной скорости.
![]() Рис.15. Движение точки в пространстве координат XYZ при \(v/c = 1,\, v_3/c = 0.1\)
|
![]() Рис.16. Другой ракурс графика (рис.15). Два оборота точки
|
Еще один вариант
![]() Рис.13. Направление глобальной скорости из (2.17) при: \(v/c = 0.1,\, v_3/c = 1\)
|
![]() Рис.14. Другой ракурс графика (рис.13)
|
Интегрируем вектор глобальной скорости.
![]() Рис.15. Движение точки в пространстве координат XYZ при \(v/c = 0.1,\, v_3/c = 1\)
|
![]() Рис.16. Другой ракурс графика (рис.15). Два оборота точки
|
Еще один вариант, в котором все скорости равны скорости света.
![]() Рис.17. Направление глобальной скорости из (2.17) при: \(v/c = 1,\, v_3/c = 1\)
|
![]() Рис.18. Другой ракурс графика (рис.17)
|
Интегрируем вектор глобальной скорости.
![]() Рис.19. Движение точки в пространстве координат XYZ при \(v/c = 1,\, v_3/c = 1\)
|
![]() Рис.20. Другой ракурс графика (рис.19). Два оборота точки
|
Согласно классической теории ОТО, в последнем варианте, где линейная и круговая скорости равны скорости света, время в подвижной системе отсчёта должно полностью остановиться.
Но если расписать по формуле (2.14) элемент \(A\), отвечающий за координату времени в подвижной системе координат,
то получится, что эта координата ненулевая, и зависит от времени в неподвижной системе координат следующим образом:
\[ A = - \int \limits_0^t \sin(\omega t) \cos(\omega t)\, dt = - {\sin(\omega t)^2 \over 4\pi} \tag{2.18}\]
Движение точки в пространстве при этом можно увидеть на рисунках 19-20.
Используемые материалы
- Википедия. Кватернион.