ИСПОЛЬЗОВАНИЕ  КВАТЕРНИОНОВ  ДЛЯ  ОПИСАНИЯ  ВРАЩАТЕЛЬНОГО  ДВИЖЕНИЯ

Камачкин  Александр  Михайлович

профессор,  д-р  физ.-мат.  наук,  профессор  Санкт-Петербургского  государственного  университета,  РФ,  г.  Санкт-Петербург

E-mail:  akamachkin@mail.ru

Свиркина  Лариса  Анатольевна

канд.  физ.-мат.  наук,  доцент  Санкт-Петербургского  государственного  университета,  РФ,  г.  Санкт-Петербург

E-mail:  lara_a@mail.ru

Хитров  Геннадий  Михайлович

доцент,  канд.  физ.-мат.  наук,  доцент  Санкт-Петербургского  государственного  университета,  РФ,  г.  Санкт-Петербург

E-mail: 

 

USAGE  OF  QUATERNIONS  FOR  THE  DESCRIPTION  OF  THE  ROTATIONAL  MOTION

Alexander  Kamachkin

professor,  doctor  of  physical  and  mathematical  sciences,  professor  of  St.  Petersburg  State  University,  Russia,  St.  Petersburg

Larisa  Svirkina

candidate  of  physical  and  mathematical  sciences,  associate  professor  of  St.  Petersburg  State  University,  Russia,  St.  Petersburg

Gennady  Khitrov

associate  professor,  candidate  of  physical  and  mathematical  sciences,  associate  professor  of  St.  Petersburg  State  University,  Russia,  St.  Petersburg

 

АННОТАЦИЯ

Статья  посвящена  описанию  свойств  кватернионов,  позволяющих  использовать  их  для  описания  вращательного  движения,  как  управляемого  вращательного  движения  космических  аппаратов,  так  и  в  компьютерных  играх.

ABSTRACT

The  paper  is  devoted  to  the  description  of  properties  of  the  quaternions,  allowing  us  to  use  them  for  the  description  of  a  rotational  motion  as  a  controlled  rotational  motion  of  spacecrafts,  and  in  computer  games.

 

Ключевые  слова:  кватернионы;  сложение  поворотов;  вращательное  движение;  кинематические  уравнения.

Keywords:  quaternions;  addition  of  turns;  rotational  motion;  kinematic  equations.

 

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

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

Пусть    и    —  кватернионы  или  «числа»  с  тремя  мнимыми  единицами  .  В  записи  кватернионов  наборы    и    —  вещественные  числа.  Числа    и    называются  скалярными  частями  кватернионов,  а    и  —  векторными  частями.  Умножение  кватернионов  как  чисел  с  тремя  мнимыми  единицами  производится  как  обычное  умножение  «скобки»  на  «скобку»,  т.  е.  как    с  учетом  перемножения  мнимых  единиц  между  собой:  .  Произведение  кватернионов  в  векторной  форме  имеет  вид:  ,  где  и    означают  скалярное  и  векторное  произведения    на    соответственно.  При  этом    и    могут  рассматриваться  как  векторы  трехмерного  пространства  с  декартовой  системой  координат  с  ортами  .  Отметим,  что    будет  скалярной,  а    —  векторной  частями  произведения  . 

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

 

.

 

Рассмотрим  функцию  .

 

 

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

 

,

 

где  ,  фиксированный  вектор  единичной  длины  (),  и  угол    —  параметры,  связанные  друг  с  другом  соотношением  ,  с  помощью  которых  задается  функция  .

Дадим  геометрическую  интерпретацию  функции  .  Для  этого  найдём  проекции  векторов    и    на  орты  ортогональной  системы  координат:  ,  ,  .  Начнем  с  вектора  .  Чтобы  найти  указанные  проекции,  достаточно  найти  скалярные  произведения  ,  ,  .  Поскольку    и    заданные  векторы,  то  будем  считать,  что    —  известно.  Найдем  :  ,  т.  е.    лежит  в  плоскости  проходящей  через  орты  p  и  r.  Найдем  теперь  :  .

Обозначим  проекцию  вектора    на  ось  с  ортом  r  через  ,  т.  е.  положим  ,  Очевидно,  поскольку  ,  что  проекция  вектора    на  плоскость  с  ортами  q  и  r  будет  совпадать  с  проекцией    на  ось  с  ортом  r,  т.  е.  будет  равна  .  Но  эта  же  проекция  будет  равна  ,  т.е.  .  (Мы  воспользовались  здесь  определениями  скалярного  и  векторного  произведения  векторов  из  аналитической  геометрии,  откуда  следует,  что    и  ,  где    -  угол  между  векторами  r  и  ).

 

Рисунок  1.  Взаимное  расположение  векторов 

 

Найдем  теперь  скалярные  произведения  ,    и  :

 

,

,

 

Сравним  теперь  координаты  векторов    и  h  в  декартовой  системе  координат  с  ортами  p,  q,  r:

 

,

 

Сравнивая  координаты  двух  указанных  векторов,  мы  видим,  что  вектор  h  получен  из  вектора    поворотом  последнего  на  угол    вокруг  оси  с  ортом  p.

 

Рисунок  2.  Поворот  вектора    на  угол 

 

Отметим,  что  обычно  единичный  кватернион    (кватернион  с  ),  используемый  в  преобразованиях  ,  записывают  в  виде  ,  где    и  p  определяются  из  соотношений:    ,  а  .  Действительно,  поскольку  ,  то  при  этих  предположениях  имеем:  ,  где  ,  а  .

Итак,  преобразование  ,  с  ,  поворачивает  вектор    вокруг  оси  определяемой  вектором    на  угол  ,  определяемый  из  условия:  .  Очевидно,  что  эта  интерпретация  сохраняет  смысл  и  для  предельных  случаев    и  .  В  первом  случае  не  происходит  ничего:  ,  поскольку  .  Во  втором  случае  -    происходит  поворот  на  угол    ()  безразлично  вокруг  какой  оси  ().

Предположим,  что  .  Совершим  теперь  поворот  вектора  h  вокруг  некоторой  оси  на  некоторый  угол,  определяемые  единичным  кватернионом    ().  Преобразование  можно  записать  в  виде:  .  То  есть,  два  последовательных  поворота  вектора  ,  определяемых  параметрами  единичных  кватернионов    и  ,  равносильны  одному  повороту,  определяемому  параметрами  (ось  и  угол)  единичного  кватерниона    ().  Операция  замены  двух  последовательных  поворотов  одним  поворотом  называется  сложением  поворотов.  Понятно,  что  эту  операцию  сложения  с  двух  поворотов  можно  обобщить  на  произвольное  число  поворотов. 

Кинематические  уравнения  движения  твердого  тела  в  декартовых  координатах.  При  изучении  вращательного  движения  твердого  тела,  как  правило,  требуется  вычислить  его  ориентацию,  то  есть  положение  системы  координат  связанной  с  телом  относительно  какой-либо  другой,  чаще  всего  инерциальной  системы  координат.  Это,  так  сказать,  «полная»  ориентация.  Иногда  требуется  знать  лишь  «частичную»  ориентацию,  например,  направление  какой-либо  оси  в  теле  относительно  выбранного  направления  в  пространстве.

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

Для  составления  кинематических  уравнений  заметим  следующее:  пусть    орт  i-той  оси  координат  в  инерциальном  пространстве,  а    —  вектор  угловой  скорости  тела.  Тогда  в  системе  координат    для  ортов    будем  иметь  уравнения  [3,  с.  410]

 

,                                              (1)

 

при  этом  на    и    в  равенстве  (1)  можно  смотреть  как  на  координатные  векторы  в  системе  координат  связанной  с  телом.

Итак    —  есть  не  что  иное,  как  столбцы  матрицы  ,  тем  самым  (1)  есть  система  дифференциальных  уравнений  для  определения  элементов  матрицы    как  функций  времени.  Система  уравнений  (1)  может  быть  записана  в  матричном  виде 

 

                                              (2)

 

(под  векторным  произведением  вектора    на  матрицу    понимается  векторное  произведение  указанного  вектора    на  каждый  столбец  матрицы  ).

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

 

,

 

то  матричные  соотношения  (2)  можно  записать  так: 

 

                                                (3)

 

Полученная  система  кинематических  уравнений  в  форме  (1)  (или,  что  то  же  самое,  в  форме  (2)  или  (3))  не  всегда  удобна,  так  как  она  имеет  размерность  равную  девяти.  С  другой  стороны  эти  девять  величин  связаны  шестью  независимыми  соотношениями  .  Следовательно,  положение  тела  определяется  всего  тремя  независимыми  параметрами.  Поэтому  представляют  интерес  системы  кинематических,  независимых  уравнений,  которые  имеют  меньшую  размерность.  Наименьшая  размерность  равна  трем  и  связана  она,  например,  с  углами  Эйлера.  Кинематические  уравнения  базирующиеся  на  использовании  углов  Эйлера  сугубо  нелинейные.  Поэтому  представляют  интерес  уравнения  обладающие  свойствами  уравнения  (3),  но  имеющие  меньшую  размерность.

Кинематические  уравнения  в  параметрах  Родрига–Гамильтона.  Известно,  что  для  твердого  тела,  имеющего  неподвижную  точку,  переход  от  одного  положения  к  другому  можно  совершить  с  помощью  поворота  вокруг  некоторой  оси  на  вполне  определенный  угол,  скажем  ,  в  одном  направлении,  и  на  угол    в  противоположном  направлении  вращения.  Достигают  ясности  в  этом  вопросе,  как  известно,  выбором  направляющего  вектора  оси  вращения  p  и  договоренности  о  выборе  положительного  направления  вращения,  то  есть  вращения  против  часовой  стрелки,  если  смотреть  с  конца  направляющего  вектора.  При  указанных  условиях  задание  орта  р  и  угла  поворота    однозначно  определяет  положение  тела.

В  теории  конечных  поворотов  и  в  механике  принято  вместо  р  и    рассматривать  следующие  величины:

 

  и  .

 

  и  координаты    вектора    в  какой  либо  правой  декартовой  системе  координат  называют  параметрами  Родрига-Гамильтона  [8,  с.  104].  Параметры  связаны  очевидным  соотношением:

 

,                                            (4)

 

и,  по  сути,  являются  компонентами  «единичного  кватерниона»  .

Орт  оси  поворота  р  и  угол  поворота    выражаются  через  эти  параметры  следующим  образом:

 

 

и  если  ,  то

 

.

 

Если  ,  то  это  соответствует  случаю  поворота  вокруг  произвольной  оси  на  углы  равные  нулю,  либо  кратные    радиан.  (В  этом  случае  из  соотношения  (4)  следует,  что    и  что  ось  поворота  не  определена;  или  начальное  и  конечное  положения  тела  совпадают.)

Таким  образом,  параметры  Родрига-Гамильтона  однозначно  определяют  конечное  положение  тела.

Замечание  1:  Обратим  внимание,  однако,  на  следующее:  если  конечное  положение  характеризуется  значениями  параметров  Родрига-Гамильтона

 

,                               (5)

 

то  тоже  положение  тела  характеризуется  и  значениями  параметров

 

.                       (6)

 

Действительно,  пусть  положение  тела  характеризуется  значениями  параметров  (5),  которым  соответствует  орт  р  и  угол  поворота  .  Но  то  же  положение  характеризуется  также  ортом    и  углом    и,  следовательно,  учитывая  что

 

,  ,

 

значениями  параметров  (6).

Рассмотрим  вновь  два  положения  твердого  тела,  точнее  два  положения  правой  декартовой  системы  координат,  связанной  с  ним.  Два  положения  системы  координат  можно  рассматривать  как  две  системы  координат    и  .  Пусть  S  матрица  перехода  от  системы    к  .  Матрица  S  полностью  определяет  как  ось  поворота,  так  и  угол  поворота  и,  следовательно,  параметры  Родрига-Гамильтона.

Чтобы  яснее  представлять  как  ось  и  угол  определяются  через  элементы  матрицы  S,  напомним  некоторые  ее  свойства:

Матрица  S  —  ортогональная  матрица,  то  есть  .  Далее,  так  как  рассматриваемая  декартова  система  координат  правая,  то    и  ,  откуда  следует,  что  .

Из  ортогональности  и  вещественности  матрицы  S  следует,  что  она  имеет  характеристические  корни  (числа)  по  модулю  равные  единице,  и  если  имеет  комплексное  характеристическое  число,  то  имеет  и  комплексное  сопряженное  к  нему.  Матрица  S  —  матрица  нечетного  (третьего)  порядка  и,  следовательно,  имеет,  по  крайней  мере,  одно  вещественное  характеристическое  число.  Из  того,  что  определитель  матрицы  S  равен  произведению  его  характеристических  чисел,  и  равен  единице,  следует,  что,  по  крайней  мере,  одно  характеристическое  число  матрицы  S  равно  единице.  Отсюда  следует,  что  характеристические  числа    матрицы  S  можно  представить  в  виде:

 

.                  (7)

 

Из  теории  канонической  структуры  нормальных  операторов  и  их  матриц  следует,  что  существуют  такие  ненулевые  вещественные  вектора  a,  b  и  c,  что

 

,                         (8)

 

Причем

 

.                                     (9)

 

Из  (8)  и  (9)  видно,  что,  не  нарушая  общности,  можно  считать  вектора  a,  b  и  c  ортами  правой  декартовой  системы  координат  .  (Если  это  не  так,  то  система  a,  с,  b  будет  правой,  либо,  заменив  вектор  а  на  ,  что  не  отразится  на  равенствах  (8),  получим  правую  тройку  .)

Построим  матрицу  Т  из  векторов  a,  b,  c  так,  что  .  Матрица  Т  будет  матрицей  перехода  от  системы  координат    к  .

Преобразование,  определяемое  матрицей  S,  переводит  систему    в  новое  положение,  скажем,  .  Очевидно,  что  Т  будет  также  матрицей  перехода  от  системы    к  .

Матрица  перехода    от  системы    к    будет  задавать  то  же  преобразование,  что  и  матрица  S.  Рассмотрим  цепочку  переходов

 

,

 

видим,  что  .  Найдем  явный  вид  матрицы  ,  для  чего  выпишем  вначале  следствие  из  формулы  (8):

 

                                             (8’)

 

Затем  вычислим  произведение  ST,  используя  формулы  (8)  и  (8’),  получим 

 

.

 

Теперь  подставляя  произведение  ST  в  выражение 

 

.                            (10)

 

Из  вида  матрицы    следует,  что  она  задает  преобразование  поворота  на  угол    вокруг  первой  координатной  оси,  то  есть  вокруг  оси  с  ортом  а.  (Если  a,  c,  b  образуют  правую  тройку  векторов,  то  берем    и  получаем  ,  которая  будет  транспонированной  к  выписанной  выше,  и  рассматриваемое  преобразование  будет  поворотом  вокруг  оси  с  ортом  а  на  угол  ,  или  на  угол    вокруг  оси  с  ортом  ).  Следовательно,  ось  поворота  задается  собственным  вектором  матрицы    отвечающим  собственному  числу  равному  единице.

Укажем  простые  способы  нахождения  собственного  вектора  d  матрицы  ,  соответствующего  собственному  числу  равному  единице  (d  —  не  обязательно  нормированный).

Из  равенств    и    следует,  что  .  Из  равенств  же    и    следует,  что  d  является  решением  системы

 

                                                  (11)

 

с  кососимметричной  матрицей  коэффициентов  равной  .

Если    (то  есть  S  несимметричная  матрица),  то  ненулевое  решение  системы  (11)  будет  собственным  вектором  матрицы  S,  соответствующим  собственному  числу,  равному  единице.

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

Домножив  (11)  слева  на  S,  убедимся,  что  ненулевое  решение  этой  системы  является  собственным  вектором  матрицы    и,  следовательно,  матрицы  S,  отвечающим  собственному  числу  равному  единице. 

Для  отыскания  же  решения  системы  (11)  заметим,  что  если  дана  ненулевая  кососимметрическая  матрица  вида  ,  то  вектор    будет  собственным  вектором  этой  матрицы,  отвечающим  нулевому  собственному  числу.

Таким  образом,  если  S  —  несимметричная  матрица,  то  в  качестве  вектора  d  можно  взять  вектор

 

                                                 (12)

 

Если  же  S  симметричная,  то  она  либо  имеет  собственные  числа  равные  ,  и  тогда  их  два  (определитель  матрицы  S  равен  единице),  либо  не  имеет,  и  тогда  ,  в  любом  случае  .

В  случае  симметрической  матрицы  S  не  совпадающей  с  E,  у  матрицы  S  только  три  элемента  отличны  от  нуля,  которые  равны  либо  +1,  либо  –1.  При  этом  есть  как  +1,  так  и  минус  1,  поскольку    и  .  В  этом  случае,  для  любой  конкретной  матрицы  S,  вектор  d  находится  весьма  просто.

Заметим,  что  преобразование  с  симметричной  матрицей  S,  имеющей  собственные  числа  равные  ,  будет  вращением  на  угол    или    вокруг  собственного  вектора  соответствующего  собственному  числу  равному  единице.

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

 

,

 

где  —  координаты  орта  оси  поворота  в  указанных  системах.

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

 

.

 

Здесь  орт  оси  поворота    в  указанных  системах  будет  такой:  .

Векторa    и    связаны  очевидным  соотношением:  ,  поскольку    -  координатный  вектор  в  системе  координат  ,  а    —  координатный  вектор  в  системе  координат  .  Аналогичным  соотношением  и  по  тем  же  соображениям  связаны  векторы  d  и  ,  т.  е.  ,  где    построен  по  матрице    таким  же  образом,  как  d  по  матрице  S. 

Используя  (10)  и  (12)  получаем,  что

 

.

 

Отсюда  уже  легко  получается

 

                                         (13)

 

Из  известного  факта  теории  матриц,  что  сумма  собственных  чисел  матрицы  S  равна    следует,  что

 

,  или  .                         (14)

 

Формулы  (13)  и  (14),  расписанные  покоординатно,  и  дают  выражение  параметров  Родрига-Гамильтона  через  элементы  матрицы  S  в  случае,  если  она  несимметричная.

 

                                             (15)

 

Замечание  2.  Формулы  (15)  справедливы  и  в  том  случае,  когда  S  симметричная  матрица.  Однако,  в  этом  случае  из  них  нельзя  получить  выражений  для  ,  если  S  не  совпадает  с  E.

Замечание  3.  Как  и  в  замечании  1,  при  выражении    из  (14)  знак  перед  квадратным  корнем  можно  брать  любой.

Если  матрица  S  симметричная  и  не  совпадает  с  единичной,  то  выражение  для  параметров  Родрига-Гамильтона  имеет  следующий  вид

 

,                                   (16)

 

где  d  нормированный  собственный  вектор,  отвечающий  собственному  числу  равному  единице. 

Если  ,  то  эти  параметры  следующие:

 

.                                                        (17)

 

Найдем  теперь  выражения  элементов  матрицы  S  через  параметры  Родрига-Гамильтона.  Для  этого  кроме  соотношений  (15),  которых  явно  недостаточно,  рассмотрим  так  же  при    следующие  соотношения:    или  вытекающие  из  него:

 

,                                                (18)

                     (19)

 

(соотношения  (18)  и  (19)  обращаются  при    в  тождество).

Соотношение  (19)  получено  из  наглядного  геометрического  соображения,  что  вектор  ,  ортогональный  к  вектору    и  i-тому  орту  системы  координат    при  повороте  вокруг  оси  с  направляющим  вектором    на  угол    перейдет  в  вектор  ортогональный  к    и  i-тому  орту  системы  .  Слева  в  (19)  стоит  скалярное  произведение  двух  векторов  одинаковой  длины  (один  получен  из  другого  поворотом  вокруг  заданной  оси).  Как  известно  скалярное  произведение  векторов  равно  произведению  длин  этих  векторов  (в  нашем  случае  квадрату  длины)  на  косинус  угла  между  этими  векторами  (выражение  справа).  Осталось  убедиться,  что  квадрат  длины  вектора    т.  е.  выражение    будет  равно  .  Действительно,

 

 

Теперь  из  соотношений  (18)  и  (19)  легко  находятся  .  Действительно,  (19)  в  развернутом  виде  дает:

 

                                           (19’)

 

Так  как    и  ,  то  (19)  в  развернутом  виде  дает:  левая  часть  – 

 

 

правая  часть  —  .  Подставляя  их  и  приравнивая,  получаем:

 

                                  (19’)

 

Из  (19’)  получаем:

 

 

Полагая    распишем  последние  соотношения  подробно:

 

 

Для  определения  остальных  элементов  матрицы  S  поступим  следующим  образом.  Разрешим  соотношение    относительно  ,  .  Учитывая  выражения  для  ,  получим:

 

                                            (20

 

Теперь  из  (15)  и  (20)  совсем  просто  получить  выражения  для,  .

Итак

 

    (21)

 

Соотношения  (3),  (4),  (15),  (21)  позволяют,  наконец,  получить  искомые  кинематические  уравнения  в  параметрах  Родрига-Гамильтона.  Для  этого  продифференцируем  (15)  в  силу  системы  (3),  затем  исключим  с  помощью  (15)  и  (21)  элементы  в  полученных  соотношениях.  В  итоге  получим  дифференциальные  уравнения,  разрешив  которые  относительно  ,  приведем  их  к  искомым  кинематическим  уравнениям.

Итак,  берем  первое  уравнение  системы  (15)  ,  получим:

 

 

Откуда

 

                             (22)

 

Далее  берем  второе  уравнение  системы  (15)  ,  имеем:

 

 

Используя  (21)  и  (22)  получим  далее 

 

 

Откуда

 

                                (23)

 

Совершенно  аналогично  получаем:

 

                                (24)

                                (25)

 

Соотношения  (22),  (23),  (24),  (25)  и  есть  искомые  кинематические  уравнения.  Используя  векторную  форму  записи,  эти  соотношения  записываются  более  компактно:

 

                                             (26)

 

Уравнения  (26)  задают  кинематику  движения  системы    относительно    в  координатах  системы  .

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

 

 

Если  учесть,  что  ,  (—  те  же,  что  в  уравнении  (26))  и  ввести  обозначение  ,  то  есть  ,  то  в  этих  обозначениях  получим:

 

                                     (27)

 

-координатный  в  системе    вектор,  соответствующий  вектору  угловой  скорости  системы    относительно  .

Системы  (26)  и  (27)  в  расписанном  покоординатном  виде  обычно  приводятся  во  многих  руководствах  по  аналитической  механике  [4,  с.  618—619]  (см.  также  [1]).  Некоторое  преимущество  кинематических  уравнений  в  форме  (26)  по  сравнению  с  (3)  очевидно,  то  есть  употребление  параметров  Родрига-Гамильтона  вполне  оправдано.  Покажем  также,  что  знание  параметров  Родрига-Гамильтона,  характеризующих  переход  системы    в    позволяет  также  просто  определять  конечное  положение  любого  вектора,  как  и  значение  матрицы  перехода  от    в  .  Кроме  того,  если  мы  совершаем  переход  системы    в    через  промежуточную  систему  ,  то  есть 

 

                                      (28)

 

где    —  соответствующие  матрицы  перехода,  то  знание  параметров    перехода    в    и    перехода    в  ,  позволяет  так  же  просто,  если  еще  не  проще,  определить  параметры    перехода    в    как  и  матрицу  перехода  S  от    к    через  .

Начнем  с  последнего  утверждения.

Из  (28)  видно,  что:

 

                                                (29)

 

Выразим  элементы  матриц    через  соответствующие  параметры  по  формулам  (21).  Тогда  из  (29)  видно,  что  элементы  матрицы  S  будут  алгебраическими  выражениями  от  .  Эти  выражения,  после  упрощений  на  основе  соотношения  (4),  подставим  в  (15).  Тогда  для  параметров    получим  следующие  решения:

 

                                                (30)

 

 

Другое  решение  (см.  замечание  3)  получится,  если  все  правые  части  системы  (30)  взять  с  обратным  знаком.  Соотношения  (30)  в  векторной  форме  можно  записать  следующим  образом:

 

                                    (31)

 

Решение  (30)  (или  (31),  что,  то  же  самое),  а  так  же  уравнения  (26),  (27)  можно  записать  в  еще  более  простом  виде,  если  использовать  кватернионы. 

Вернемся  к  формулам  (30),  предварительно  сопоставив  набору  параметров    кватернион  ,  набору  параметров    кватернион    и  параметрам  —  кватернион  .  Из  формул  (30)  или  (31)  непосредственно  следует,  что 

 

.

 

Вернемся  к  кинематическим  уравнениям  (26)  и  дополним  их  динамическими  уравнениями  вращательного  движения  твердого  тела  в  системе  координат,  связанной  с  телом,  и  центр  которой  находится  в  центре  инерции  тела  [3,  с.  410]

 

,                                 (32)

 

где:  —  тензор  инерции  твердого  тела, 

  —  угловая  скорость, 

  —  момент  внешних  сил  приложенных  к  телу.  Объединяя  (32)  и  (26)  получим  систему  уравнений

 

                                     (33)

 

 

 

полностью  описывающую  вращательное  движение  твердого  тела.  Система  уравнений  (33)  удобна  в  том  смысле,  что  позволяет  строить  для  неё  достаточно  простые  функции  Ляпунова.  Момент    в  уравнениях  (33),  рассматриваемый  как  управляющий,  можно  выбирать  в  силу  функций  Ляпунова  таким  образом,  чтобы  он  решал,  например,  задачу  полной  ориентации  твердого  тела  (см.,  например  [10]).

Заключение.  Использование  кватернионов  в  задачах  ориентации  и  управления  вращательным  движением  твердого  тела,  позволяет  усовершенствовать  результаты  работ  [2,  5,  6]  на  базе  уравнений  (26).

В  последнее  время  кватернионы  находят  все  большее  применение  в  компьютерных  технологиях  (см.  например  [7,  9,  11]).

       

Список  литературы:

1.Бранец  В.Н.,  Шмыглевский  И.П.  Применение  кватернионов  в  задачах  ориентации  твердого  тела.  М.:  Наука,  1973.  —  320  с.

2.Ермолин  В.С.,  Жабко  А.П.,  Камачкин  А.М.,  Овсянников  Д.А.  Школа  В.И.  Зубова  по  теории  устойчивости  и  теории  динамических  и  механических  систем  //  Труды  СВМО.  —  2005.  —  Т.  7.  —  №  1.  —  С.  426—432.

3.Зубов  В.И.  Лекции  по  теории  управления.  М.:  Наука,  1975.  —  496  с.

4.Ишлинский  А.Ю.  Ориентация,  гироскопы  и  инерциальная  навигация.  М.:  Наука,  1976.  —  672  с.

5.Камачкин  А.М.  Методы  оценки  точности  позиционирования  грузовой  платформы  при  ограниченных  возмущениях  //  Автоматизация  морских  судов  и  технологических  процессов  в  судостроении.  ВНТО  им.  акад.  А.Н.  Крылова.  Л.:  Судостроение.  —  1989.  —  Вып.  478.  —  С.  15—22. 

6.Камачкин  А.М.,  Шамберов  В.Н.  Метод  декомпозиции  в  многомерных  нелинейных  динамических  системах  //  Вестник  ВГУ.  Сер.:  Системный  анализ  и  информационные  технологии.  —  2012.  —  Вып.  1.  —  С.  47—55.

7.Кватернионы  в  геометрии,  механике,  релятивистской  физике,  теории  поля  //  fizteh.ru:  информационный  портал  [Электронный  ресурс]  —  Режим  доступа.  —  URL:  http://www.fizteh.ru/02-07-90327/index/qwat/  (дата  обращения  25.01.2014).

8.Лурье  А.И.  Аналитическая  механика.  М.:  Физматгиз,  1961.  —  826  с.

9.Побегайло  А.П.  Применение  кватернионов  в  компьютерной  графике.  Минск:  БГУ,  2010.  —  216  с.

10.Хитров  Г.М.  Применение  второго  метода  Ляпунова  к  решению  задачи  ориентации  твердого  тела  //  Дифференциальные  и  интегральные  уравнения:  межвуз.  сб.,  Горький:  ГГУ,  1981.  —  С.  152—155.

11.Цисарж  В.В.,  Марусик  Р.И.  Математические  методы  компьютерной  графики:  учебное  пособие.  К.:  Факт,  2004.  —  466  с.