Матрица поворота

В двумерном пространстве поворот можно описать одним углом ${\displaystyle \;\theta }$ со следующей матрицей линейного преобразования в декартовой системе координат:

${\displaystyle M(\theta )={\begin{pmatrix}\cos {\theta }&\mp \sin {\theta }\\\pm \sin {\theta }&\cos {\theta }\end{pmatrix}}}$

или

${\displaystyle M(\theta )=\exp \left({\theta }\cdot {\begin{bmatrix}0&\mp 1\\\pm 1&0\end{bmatrix}}\right)}$.

Поворот выполняется путём умножения матрицы поворота на вектор-столбец, описывающий вращаемую точку:

${\displaystyle {\begin{bmatrix}x'\\y'\\\end{bmatrix}}={\begin{bmatrix}\cos \theta &\mp \sin \theta \\\pm \sin \theta &\cos \theta \\\end{bmatrix}}{\begin{bmatrix}x\\y\\\end{bmatrix}}.}$

Координаты (x′,y′) в результате поворота точки (x, y) имеют вид:

${\displaystyle x'=x\cos \theta \mp y\sin \theta ,}$

${\displaystyle y'=\pm x\sin \theta +y\cos \theta .}$

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

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

Матрицами вращения вокруг оси декартовой системы координат на угол ${\displaystyle \alpha }$ в трёхмерном пространстве с неподвижной системой координат являются:

• Вращение вокруг оси ${\displaystyle x}$:

${\displaystyle M_{x}(\alpha )={\begin{pmatrix}1&0&0\\0&\cos \alpha &-\sin \alpha \\0&\sin \alpha &\cos \alpha \end{pmatrix}}.}$

• Вращение вокруг оси ${\displaystyle y}$:

${\displaystyle M_{y}(\alpha )={\begin{pmatrix}\cos \alpha &0&\sin \alpha \\0&1&0\\-\sin \alpha &0&\cos \alpha \end{pmatrix}}.}$

• Вращение вокруг оси ${\displaystyle z}$:

${\displaystyle M_{z}(\alpha )={\begin{pmatrix}\cos \alpha &-\sin \alpha &0\\\sin \alpha &\cos \alpha &0\\0&0&1\end{pmatrix}}.}$

Положительным углам при этом соответствует вращение вектора против часовой стрелки в правой системе координат, и по часовой стрелке в левой системе координат, если смотреть против направления соответствующей оси[2]. Например, при повороте на угол ${\displaystyle \alpha =90^{\circ }}$ вокруг оси ${\displaystyle z}$ ось ${\displaystyle x}$ переходит в ${\displaystyle y}$: ${\displaystyle M_{z}(90^{\circ })\cdot \mathbf {e} _{x}=\mathbf {e} _{y}}$. Аналогично, ${\displaystyle M_{y}(90^{\circ })\cdot \mathbf {e} _{z}=\mathbf {e} _{x}}$ и ${\displaystyle M_{x}(90^{\circ })\cdot \mathbf {e} _{y}=\mathbf {e} _{z}}$. Правая система координат связана с выбором правого базиса (см. правило буравчика).

Совершенно аналогично могут быть записаны матрицы поворота конечномерного пространства любой более высокой размерности.

Надо только иметь в виду, что для размерностей пространства, не равных трём, невозможно указать единственную прямую, ортогональную двум данным прямым, а поэтому нельзя говорить о вращении вокруг какой-то оси, можно же говорить о вращении в какой-то плоскости[3]. Все точки при повороте в пространстве любой размерности, начиная с 2, всегда движутся параллельно некоторой (двумерной) плоскости.

Итак, совершенно аналогично трёхмерному случаю (с приведенной оговоркой) можем написать матрицу поворота в любой координатной плоскости для любой размерности пространства.

Например:

${\displaystyle M_{1,2}(\alpha )={\begin{pmatrix}\cos \alpha &-\sin \alpha &0&0&0\\\sin \alpha &\cos \alpha &0&0&0\\0&0&1&0&0\\0&0&0&1&0\\0&0&0&0&1\end{pmatrix}}}$

— матрица поворота в 5-мерном пространстве в плоскости ${\displaystyle x_{1}x_{2}}$,

${\displaystyle M_{2,4}(\alpha )={\begin{pmatrix}1&0&0&0&0&0&0\\0&\cos \alpha &0&-\sin \alpha &0&0&0\\0&0&1&0&0&0&0\\0&\sin \alpha &0&\cos \alpha &0&0&0\\0&0&0&0&1&0&0\\0&0&0&0&0&1&0\\0&0&0&0&0&0&1\end{pmatrix}}}$

— матрица поворота в 7-мерном пространстве в плоскости ${\displaystyle x_{2}x_{4}}$.

• При таком подходе знаки перед синусами расставлять даже легче, поскольку они определяются порядком перечисления осей плоскости вращения: какая названа первой, в той строке перед синусом минус.
• Легко видеть, что матрица поворота в плоскости ${\displaystyle x_{1}x_{2}}$ совпадает (что естественно) с матрицей поворота в плоскости ${\displaystyle x_{2}x_{1}}$ и т. д. с точностью до замены угла поворота на противоположный.
• Поэтому такие матрицы с переставленными индексами очевидно не независимы, и для получения произвольного поворота достаточно включить в композицию каждую плоскость только один раз, то есть, скажем, только ${\displaystyle M_{1,2}(\alpha _{1,2})}$, а не ${\displaystyle M_{1,2}(\alpha _{1,2})}$ и ${\displaystyle M_{2,1}(\alpha _{2,1})}$.
• Исходя из этого, нетрудно сосчитать их общее количество: ${\displaystyle n(n-1)/2}$, где n — размерность пространства.

Пусть ${\displaystyle \;M}$ — матрица поворота вокруг оси с ортом ${\displaystyle \;n}$ на угол ${\displaystyle \;\alpha }$, ${\displaystyle \;M'}$ — матрица поворота вокруг оси с ортом ${\displaystyle \;n'}$ на тот же угол, причем

${\displaystyle \;n'=M''\cdot \;n,}$

где ${\displaystyle \;M''}$ — матрица поворота, изменяющая орт оси поворота ${\displaystyle \;n}$. Тогда

${\displaystyle M'=M''\cdot M\cdot M''^{\;T},}$

где ${\displaystyle \;M''^{\;T}}$ — транспонированная матрица ${\displaystyle \;M''}$.

Если ${\displaystyle \;M_{1}}$ — матрица поворота вокруг оси с ортом ${\displaystyle \;n}$ на угол ${\displaystyle \;\alpha }$, ${\displaystyle \;M_{2}}$ — матрица поворота вокруг оси с ортом ${\displaystyle \;m}$ на угол ${\displaystyle \;\beta }$, то ${\displaystyle M_{2}\cdot M_{1}}$ — матрица, описывающая поворот, являющийся результатом двух последовательно осуществленных поворотов (${\displaystyle \;M_{1}}$ и ${\displaystyle \;M_{2}}$), поскольку

${\displaystyle (M_{2}\cdot M_{1})\cdot r=M_{2}\cdot (M_{1}\cdot r).}$

При этом последовательность поворотов можно поменять, видоизменив поворот ${\displaystyle \;M_{1}}$:

${\displaystyle M_{2}\cdot M_{1}=M'_{1}\cdot M_{2},}$

где матрица ${\displaystyle \;M'_{1}}$ — матрица поворота на угол ${\displaystyle \;\alpha }$ вокруг оси c ортом ${\displaystyle \;n',}$ повернутым с помощью поворота ${\displaystyle \;M_{2}}$:

${\displaystyle n'=M_{2}\cdot n,\qquad M'_{1}=M_{2}\cdot M_{1}\cdot M_{2}^{T},}$

поскольку ${\displaystyle M_{2}^{T}\cdot M_{2}=E}$, так как матрица поворота является ортогональной матрицей (${\displaystyle \;E}$ — единичная матрица). Заметим, что коммутативности поворотов в обычном смысле нет, то есть

${\displaystyle M_{2}\cdot M_{1}\neq M_{1}\cdot M_{2}.}$

Последовательные повороты около осей ${\displaystyle \;Z,X',Z''}$ на угол прецессии (${\displaystyle \alpha }$), угол нутации (${\displaystyle \beta }$) и на угол собственного вращения (${\displaystyle \gamma }$) приводят к следующему выражению для матрицы поворота:

${\displaystyle M(\alpha ,\beta ,\gamma )={\begin{pmatrix}\cos \alpha \cos \gamma -\sin \alpha \cos \beta \sin \gamma &-\cos \alpha \sin \gamma -\sin \alpha \cos \beta \cos \gamma &\sin \alpha \sin \beta \\\sin \alpha \cos \gamma +\cos \alpha \cos \beta \sin \gamma &-\sin \alpha \sin \gamma +\cos \alpha \cos \beta \cos \gamma &-\cos \alpha \sin \beta \\\sin \beta \sin \gamma &\sin \beta \cos \gamma &\cos \beta \end{pmatrix}}.}$

Ось ${\displaystyle \;X'}$ — ось X, повёрнутая первым поворотом (на ${\displaystyle \;\alpha }$), ${\displaystyle \;Z''}$ — ось Z, повёрнутая первым и вторым поворотом (на ${\displaystyle \;\alpha }$ и ${\displaystyle \;\beta }$). Вследствие перестановочности поворотов приведённая матрица соответствует поворотам на углы ${\displaystyle \;\gamma }$, ${\displaystyle \;\beta }$, ${\displaystyle \;\alpha }$ вокруг осей Z, X, Z:

${\displaystyle M(\alpha ,\beta ,\gamma )=M_{Z}(\alpha )\cdot M_{X}(\beta )\cdot M_{Z}(\gamma )}$.

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

• 1) Поворот около осей: ${\displaystyle X,Y,X}$
• 2) Соответственно: ${\displaystyle X,Y,Z}$
• 3) ${\displaystyle X,Z,X}$
• 4) ${\displaystyle X,Z,Y}$
• 5) ${\displaystyle Y,X,Y}$
• 6) ${\displaystyle Y,X,Z}$
• 7) ${\displaystyle Y,Z,X}$
• 8) ${\displaystyle Y,Z,Y}$
• 9) ${\displaystyle Z,X,Y}$
• 10) ${\displaystyle Z,X,Z}$
• 11) ${\displaystyle Z,Y,X}$
• 12) ${\displaystyle Z,Y,Z}$

Пусть ось вращения задана единичным вектором ${\displaystyle {\hat {\mathbf {v} }}=(x,y,z)}$, а угол поворота ${\displaystyle \theta }$.

Тогда матрица поворота в декартовых координатах имеет вид:

${\displaystyle M({\hat {\mathbf {v} }},\theta )={\begin{pmatrix}\cos \theta +(1-\cos \theta )x^{2}&(1-\cos \theta )xy-(\sin \theta )z&(1-\cos \theta )xz+(\sin \theta )y\\(1-\cos \theta )yx+(\sin \theta )z&\cos \theta +(1-\cos \theta )y^{2}&(1-\cos \theta )yz-(\sin \theta )x\\(1-\cos \theta )zx-(\sin \theta )y&(1-\cos \theta )zy+(\sin \theta )x&\cos \theta +(1-\cos \theta )z^{2}\end{pmatrix}}}$

Если задан кватернион ${\displaystyle q=(w,x,y,z)}$, то соответствующая матрица поворота имеет вид:

${\displaystyle Q={\begin{bmatrix}1-2y^{2}-2z^{2}&2xy-2zw&2xz+2yw\\2xy+2zw&1-2x^{2}-2z^{2}&2yz-2xw\\2xz-2yw&2yz+2xw&1-2x^{2}-2y^{2}\end{bmatrix}}.}$

Если ${\displaystyle \mathbf {M} }$ — матрица, задающая поворот вокруг оси ${\displaystyle {\vec {n}}}$ на угол ${\displaystyle \varphi }$, то:

• ${\displaystyle |\mathbf {M} {\vec {v}}|=|{\vec {v}}|}$ ${\displaystyle \forall {\vec {v}}}$
• ${\displaystyle \mathbf {M} {\vec {n}}={\vec {n}}}$
• ${\displaystyle (\mathbf {M} {\vec {v}},{\vec {v}})=(1-\cos \varphi )({\vec {n}}\cdot {\vec {v}})^{2}+({\vec {v}})^{2}\cos \varphi }$
• ${\displaystyle \operatorname {Tr} (\mathbf {M} )=n-2+2\cos \varphi }$ (след матрицы вращения), где n — размерность пространства (размер матрицы).
• ${\displaystyle \det \mathbf {M} =1}$ (матрица имеет единичный определитель).
• Матрица обратного поворота получается обычным транспонированием матрицы прямого поворота, т. о. ${\displaystyle \mathbf {M} ^{\mathrm {-1} }=\mathbf {M} ^{\mathrm {T} }}$.
• Для трёхмерного пространства (матриц ${\displaystyle 3\times 3}$): если строки (или столбцы матрицы) рассматривать как координаты векторов ${\displaystyle {\vec {a}},{\vec {b}},{\vec {c}}}$, то верны следующие соотношения):
  • ${\displaystyle |{\vec {a}}|=|{\vec {b}}|=|{\vec {c}}|=1}$
  • ${\displaystyle {\vec {a}}\cdot {\vec {b}}=0,{\vec {b}}\cdot {\vec {c}}=0,{\vec {c}}\cdot {\vec {a}}=0}$
  • ${\displaystyle {\vec {a}}\times {\vec {b}}={\vec {c}},{\vec {b}}\times {\vec {c}}={\vec {a}},{\vec {c}}\times {\vec {a}}={\vec {b}}}$
• Первые два свойства[4], означающие условие ортогональности матрицы, верны и для произвольной размерности пространства (размера матрицы).

• Матрица перехода
• Поворот

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

• Поворот плоскости. Матрица поворота