Кватернионы и вращения в трехмерном пространстве: алгебра Гамильтона

Алгебраическая структура кватерниона

Комплексное число состоит из одной действительной и одной мнимой части. Гамильтон расширил эту концепцию, создав кватернион — число, состоящее из одной действительной (скалярной) части и трех мнимых (векторных) частей. Кватернион записывается как q = w + x*i + y*j + z*k, где i, j, k — фундаментальные мнимые единицы. Гениальность Гамильтона заключалась в правилах перемножения этих единиц: i^2 = j^2 = k^2 = i*j*k = -1. Из этих правил следует, что умножение кватернионов некоммутативно (порядок множителей меняет результат: i*j = k, но j*i = -k). Эта некоммутативность алгебры идеально отражает некоммутативность реальных физических вращений в пространстве (повернуть книгу по оси X, а затем по Y — это не то же самое, что повернуть сначала по Y, а затем по X).

Единичные кватернионы как операторы вращения

Для описания вращений в аналитической геометрии используются исключительно единичные (нормированные) кватернионы, сумма квадратов компонент которых строго равна единице: w^2 + x^2 + y^2 + z^2 = 1. Теорема Эйлера о вращении гласит, что любое сложное перемещение твердого тела с неподвижной точкой можно описать как одно единственное вращение вокруг некоторой оси на некоторый угол. Единичный кватернион кодирует эту информацию потрясающе компактно. Если мы хотим повернуть пространство на угол theta вокруг единичного вектора оси (nx, ny, nz), мы строим кватернион: q = cos(theta/2) + sin(theta/2) * (nx*i + ny*j + nz*k). Скалярная часть хранит информацию об угле, а мнимая векторная часть — направление оси вращения.

Аналитическое выполнение вращения (Сэндвич-умножение)

Как применить это вращение к конкретной точке пространства с координатами P(x, y, z)? В матричной алгебре мы бы умножили вектор на матрицу 3x3. В алгебре кватернионов точка P превращается в «чисто мнимый» кватернион (у которого скалярная часть равна нулю): p = 0 + x*i + y*j + z*k. Само вращение выполняется через операцию, известную как сэндвич-умножение: p_new = q * p * q^(-1). То есть мы умножаем кватернион вращения q на кватернион точки p, а затем умножаем результат на кватернион, обратный к q (сопряженный кватернион). Итоговый кватернион p_new также окажется чисто мнимым, и его коэффициенты при i, j, k дадут нам точные декартовы координаты точки после вращения. Эта операция лишена любых тригонометрических особенностей, оберегая систему от математических сингулярностей.

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

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