mtimes、*

四元数の乗算

R2019b 以降

構文

quatC = A*B

説明

quatC = A*B は、A または B のいずれかが四元数である場合の四元数の乗算を実装します。A または B のいずれかはスカラーでなければなりません。

四元数の乗算を使用して回転演算子を構成できます。

• 座標系の回転のシーケンスを構成するには、目的の回転シーケンスの順序で四元数を乗算します。たとえば、p 四元数の後に q 四元数を適用するには、pq の順序で乗算します。回転演算子は $${\left(pq\right)}^{\ast }v\left(pq\right)$$ になります。ここで、v は回転するオブジェクトであり、四元数形式で指定します。* は共役を表します。

• 点の回転のシーケンスを構成するには、目的の回転シーケンスと逆の順序で四元数を乗算します。たとえば、p 四元数の後に q 四元数を適用するには、逆の順序 qp で乗算します。回転演算子は $$\left(qp\right)v{\left(qp\right)}^{\ast }$$ になります。

例

Quaternion スカラーと Quaternion ベクトルの乗算

4 行 1 列の列ベクトル A とスカラー b を作成します。A と b の乗算を行います。

A = quaternion(randn(4,4))

A = 4x1 quaternion array
      0.53767 +  0.31877i +   3.5784j +   0.7254k
       1.8339 -   1.3077i +   2.7694j - 0.063055k
      -2.2588 -  0.43359i -   1.3499j +  0.71474k
      0.86217 +  0.34262i +   3.0349j -  0.20497k

b = quaternion(randn(1,4))

b = quaternion
    -0.12414 +  1.4897i +   1.409j +  1.4172k

C = A*b

C = 4x1 quaternion array
      -6.6117 +   4.8105i +  0.94224j -   4.2097k
      -2.0925 +   6.9079i +   3.9995j -   3.3614k
       1.8155 -   6.2313i -    1.336j -     1.89k
      -4.6033 +   5.8317i + 0.047161j -    2.791k

入力引数

A — 乗算する入力

quaternion オブジェクト | quaternion オブジェクトの配列 | 実数のスカラー | 実数の配列

乗算する入力。quaternion オブジェクト、任意の次元の quaternion オブジェクトの配列、実数スカラー、または任意の次元の実数の配列として指定します。数値のデータ型は single または double でなければなりません。

B が非スカラーの場合、A はスカラーでなければなりません。

B — 乗算する入力

quaternion オブジェクト | quaternion オブジェクトの配列 | 実数のスカラー | 実数の配列

乗算する入力。quaternion オブジェクト、任意の次元の quaternion オブジェクトの配列、実数スカラー、または任意の次元の実数の配列として指定します。数値のデータ型は single または double でなければなりません。

A が非スカラーの場合、B はスカラーでなければなりません。

出力引数

quatC — 四元数の積

quaternion オブジェクト | quaternion オブジェクトの配列

四元数の積。quaternion オブジェクト、または quaternion オブジェクトの配列として返されます。

アルゴリズム

四元数の実数スカラーによる乗算

次の四元数について考えます。

$$q=a_{\text{q}}+b_{\text{q}}\text{i}+c_{\text{q}}\text{j}+d_{\text{q}}\text{k,}$$

q と実数スカラー β の積は次のようになります。

$$\beta q=\beta a_{\text{q}}+\beta b_{\text{q}}\text{i}+\beta c_{\text{q}}\text{j}+\beta d_{\text{q}}\text{k}$$

四元数の Quaternion スカラーによる乗算

四元数の基本要素の定義は次のとおりです。

$${\text{i}}^2={\text{j}}^2={\text{k}}^2=\text{ijk}=-1,$$

これを展開すると、四元数の基本要素の乗算をまとめた表を次のように埋めることができます。

	1	i	j	k
1	1	i	j	k
i	i	-1	k	−j
j	j	−k	-1	i
k	k	j	−i	-1

この表を読むときは、ij = k や ji = −k のように行を先に読みます。

2 つの四元数 $$q=a_{\text{q}}+b_{\text{q}}\text{i}+c_{\text{q}}\text{j}+d_{\text{q}}\text{k,}$$ と $$p=a_{\text{p}}+b_{\text{p}}\text{i}+c_{\text{p}}\text{j}+d_{\text{p}}\text{k}$$ が与えられた場合、乗算を次のように展開できます。

$$\begin{array}{l}z=pq=\left(a_{\text{p}}+b_{\text{p}}\text{i}+c_{\text{p}}\text{j}+d_{\text{p}}\text{k}\right)\left(a_{\text{q}}+b_{\text{q}}\text{i}+c_{\text{q}}\text{j}+d_{\text{q}}\text{k}\right)\\ =a_{\text{p}}{a}_{\text{q}}+a_{\text{p}}{b}_{\text{q}}\text{i}+a_{\text{p}}{c}_{\text{q}}\text{j}+a_{\text{p}}{d}_{\text{q}}\text{k}\\ +b_{\text{p}}{a}_{\text{q}}\text{i}+b_{\text{p}}{b}_{\text{q}}{\text{i}}^2+b_{\text{p}}{c}_{\text{q}}\text{ij}+b_{\text{p}}{d}_{\text{q}}\text{ik}\\ +c_{\text{p}}{a}_{\text{q}}\text{j}+c_{\text{p}}{b}_{\text{q}}\text{ji}+c_{\text{p}}{c}_{\text{q}}{\text{j}}^2+c_{\text{p}}{d}_{\text{q}}\text{jk}\\ +d_{\text{p}}{a}_{\text{q}}k+d_{\text{p}}{b}_{\text{q}}\text{ki}+d_{\text{p}}{c}_{\text{q}}\text{kj}+d_{\text{p}}{d}_{\text{q}}{\text{k}}^2\end{array}$$

この方程式は、四元数の乗算の表を使用して次のように単純化できます。

$$\begin{array}{l}z=pq=a_{\text{p}}{a}_{\text{q}}+a_{\text{p}}{b}_{\text{q}}\text{i}+a_{\text{p}}{c}_{\text{q}}\text{j}+a_{\text{p}}{d}_{\text{q}}\text{k}\\ +b_{\text{p}}{a}_{\text{q}}\text{i}-b_{\text{p}}{b}_{\text{q}}+b_{\text{p}}{c}_{\text{q}}\text{k}-b_{\text{p}}{d}_{\text{q}}\text{j}\\ +c_{\text{p}}{a}_{\text{q}}\text{j}-c_{\text{p}}{b}_{\text{q}}\text{k}-c_{\text{p}}{c}_{\text{q}}+c_{\text{p}}{d}_{\text{q}}\text{i}\\ +d_{\text{p}}{a}_{\text{q}}k+d_{\text{p}}{b}_{\text{q}}\text{j}-d_{\text{p}}{c}_{\text{q}}\text{i}-d_{\text{p}}{d}_{\text{q}}\end{array}$$

参照

[1] Kuipers, Jack B. Quaternions and Rotation Sequences: A Primer with Applications to Orbits, Aerospace, and Virtual Reality. Princeton, NJ: Princeton University Press, 2007.

拡張機能

C/C++ コード生成
MATLAB® Coder™ を使用して C および C++ コードを生成します。

バージョン履歴

R2019b で導入

参考

関数

• .*,times

オブジェクト

• quaternion

トピック

• Rotations, Orientation, and Quaternions