Алгебра кватернионов: математика надежных трехмерных вращений

От комплексных чисел к 4D пространству

Если комплексное число имеет одну мнимую единицу (i), то кватернион обладает сразу тремя мнимыми единицами (i, j, k) и одной вещественной частью. Кватернион записывается как q = a + bi + cj + dk. Фундаментальный закон умножения кватернионов, вырезанный Гамильтоном на мосту Брум-Бридж в Дублине, гласит: i^2 = j^2 = k^2 = ijk = -1. Из этого закона немедленно вытекает, что умножение кватернионов антикоммутативно (ij = k, но ji = -k). Именно эта некоммутативность (зависимость от порядка действий) идеально описывает физическую реальность трехмерных поворотов, ведь если мы повернем объект сначала по оси X, а затем по оси Y, мы получим совершенно иной результат, чем если сделаем это в обратном порядке.

Версоры: вращение через сопряжение

Как четырехмерное число описывает трехмерное вращение? Любой вектор трехмерного пространства (x, y, z) записывается как чисто мнимый кватернион v = xi + yj + zk (с нулевой вещественной частью). Для самого вращения используется кватернион единичной длины (версор) q = cos(theta/2) + u*sin(theta/2), где u — единичный направляющий вектор оси вращения, а theta — угол поворота вокруг этой оси. Операция вращения вектора v осуществляется с помощью так называемого двойного (сопряженного) умножения: v' = q * v * q^(-1). Эта элегантная алгебраическая формула берет трехмерный вектор, переносит его в 4D пространство, осуществляет идеальное вращение и гарантированно возвращает обратно в трехмерный мир, сохраняя его длину до микрона.

Преимущества над матрицами вращения

С точки зрения вычислительной линейной алгебры кватернионы многократно превосходят матрицы. Матрица поворота содержит 9 чисел, которые связаны 6 условиями ортогональности. При серии из тысяч поворотов (в симуляции физики) ошибки плавающей запятой накапливаются, и матрицу приходится принудительно ортогонализировать алгоритмом Грама-Шмидта, что требует огромных процессорных мощностей. Кватернион же хранит всего 4 числа (одно уравнение связи — единичная длина). Чтобы исправить накопившуюся ошибку кватерниона, его достаточно просто нормализовать — поделить на его собственную длину. Это выполняется за доли миллисекунды, экономя память и такты GPU в современных системах рендеринга.

Сферическая линейная интерполяция (SLERP)

Главным триумфом кватернионов в анимации является алгоритм SLERP (Spherical Linear Interpolation), разработанный Кеном Шумейком. Если персонажу в игре нужно плавно повернуть голову от объекта А к объекту Б, использование матриц или углов Эйлера приведет к неестественным рывкам или «проскальзыванию» скорости вращения. Поскольку единичные кватернионы геометрически лежат на поверхности четырехмерной гиперсферы, алгоритм SLERP вычисляет кратчайший путь между двумя кватернионами по дуге большого круга этой гиперсферы. Это обеспечивает абсолютно равномерную, математически идеальную угловую скорость вращения. Анимация получается кинематографически плавной, что сделало кватернионы святым Граалем компьютерной мультипликации и виртуальной реальности.