Формула поворота Родрига

Формула поворота Родрига — формула, связывающая два вектора с общим началом, один из которых получен поворотом другого на известный угол вокруг оси, проходящей через их общее начало:

${\vec {R}}_{2}-\tan(\chi /2)[{\vec {e}}\times {\vec {R}}_{2}]={\vec {R}}_{1}+\tan(\chi /2)[{\vec {e}}\times {\vec {R}}_{1}]$

где ${\vec {R}}_{1}$ — исходный вектор, ${\vec {R}}_{2}$ — результирующий вектор, ${\vec {e}}$ — единичный вектор оси поворота, $\chi$ — угол поворота. Также формула может быть записана в виде:

${\vec {R}}_{2}=({\vec {e}}\cdot {\vec {R}}_{1})(1-\cos \chi ){\vec {e}}+({\vec {e}}\times {\vec {R}}_{1})\sin \chi +{\vec {R}}_{1}\cdot \cos \chi$

Лежит в основе векторной теории конечных поворотов и сложения вращений. Получена О. Родригом в 1840 г.

Вывод

Без потери общности, направим ось ${\vec {z}}$ вдоль единичного вектора ${\vec {e}}$, а вектор ${\vec {R}}_{1}$ — лежащим в плоскости OXZ, тогда:

${\vec {R}}_{1x}={\vec {R}}_{1}-{\vec {R}}_{1z}$

${\vec {R}}_{1y}=0$

${\vec {R}}_{1z}=({\vec {e}}\cdot {\vec {R}}_{1}){\vec {e}}$

Откуда:

${\vec {R}}_{1x}={\vec {R}}_{1}-({\vec {e}}\cdot {\vec {R}}_{1}){\vec {e}}$

Положим вектор ${\vec {w}}$, равный:

${\vec {w}}={\vec {e}}\times {\vec {R}}_{1}$

Заметим, что:

$|{\vec {w}}|=|{\vec {e}}\times {\vec {R}}_{1}|=|{\vec {e}}|\,|{\vec {R}}_{1}|\sin \phi \ =|{\vec {R}}_{1}|\sin \phi$

$|{\vec {R}}_{1x}|=|{\vec {R}}_{1}|\cos(\pi /2-\phi )=|{\vec {R}}_{1}|\sin \phi .$

Тогда вектор ${\vec {R}}_{2x}$ можно выразить через векторы ${\vec {w}}$ и ${\vec {R}}_{1x}$ и угол $\chi$:

${\vec {R}}_{2x}={\vec {R}}_{1x}\cos \chi +{\vec {w}}\sin \chi =({\vec {R}}_{1}-({\vec {e}}\cdot {\vec {R}}_{1}){\vec {e}})\cos \chi +({\vec {e}}\times {\vec {R}}_{1})\sin \chi$

Результирующий вектор ${\vec {R}}_{2}$ выражается через векторы ${\vec {R}}_{2x}$ и ${\vec {R}}_{1z}$:

${\vec {R}}_{2}={\vec {R}}_{2x}+{\vec {R}}_{1z}=({\vec {R}}_{1}-({\vec {e}}\cdot {\vec {R}}_{1}){\vec {e}})\cos \chi +({\vec {e}}\times {\vec {R}}_{1})\sin \chi +({\vec {e}}\cdot {\vec {R}}_{1}){\vec {e}}$

Приведя подобные, получим формулу поворота Родрига:

${\vec {R}}_{2}=({\vec {e}}\cdot {\vec {R}}_{1})(1-\cos \chi ){\vec {e}}+({\vec {e}}\times {\vec {R}}_{1})\sin \chi +{\vec {R}}_{1}\cdot \cos \chi$

В матричной форме

Векторное умножение на вектор k можно представить в виде умножения на матрицу K:

$\mathbf {k} \times \mathbf {v} ={\begin{bmatrix}(\mathbf {k} \times \mathbf {v} )_{x}\\(\mathbf {k} \times \mathbf {v} )_{y}\\(\mathbf {k} \times \mathbf {v} )_{z}\end{bmatrix}}={\begin{bmatrix}k_{y}v_{z}-k_{z}v_{y}\\k_{z}v_{x}-k_{x}v_{z}\\k_{x}v_{y}-k_{y}v_{x}\end{bmatrix}}={\begin{bmatrix}0&-k_{z}&k_{y}\\k_{z}&0&-k_{x}\\-k_{y}&k_{x}&0\end{bmatrix}}{\begin{bmatrix}v_{x}\\v_{y}\\v_{z}\end{bmatrix}}=\mathbf {K} \mathbf {v} \,.$

Вектор v при повороте вокруг единичного вектора k перейдет в вектор

$\mathbf {v} _{\mathrm {rot} }=\mathbf {v} +(\sin \theta )\mathbf {K} \mathbf {v} +(1-\cos \theta )\mathbf {K} ^{2}\mathbf {v} =\mathbf {R} \mathbf {v} \,,$

где $\mathbf {K} (\mathbf {K} \mathbf {v} )=\mathbf {K} ^{2}\mathbf {v} =\mathbf {k} \times (\mathbf {k} \times \mathbf {v} )\,.$

Таким образом получается, что матрица поворота вокруг единичного вектора k на угол $\theta$

$\mathbf {R} =\mathbf {I} +(\sin \theta )\mathbf {K} +(1-\cos \theta )\mathbf {K} ^{2}~.$

где

$\mathbf {K} =\left[{\begin{array}{ccc}0&-k_{z}&k_{y}\\k_{z}&0&-k_{x}\\-k_{y}&k_{x}&0\end{array}}\right]~.$

Литература

• Лурье А. И.  Аналитическая механика. — М.: Физматгиз, 1961. — 824 с. — С. 101—103.
• Rodrigues O.  Des lois géométriques qui régissent les déplacements d’une système solide dans l’espace et de la variation des coordonnées provenant de ces déplacements considérés indépendamment des causes qui peuvent les produire // Liouvillés Journ. Math.. — 1840. — Vol. 5. — P. 380—440.