正方行列の三連続積の展開

小ネタです.まだとくに計算が楽になる応用もありません. 「こんな表現ができるんだ.展開の時,対称性が崩れた過程を経ないね」程度に思っていただければ成功です. もしよかったら何かに使えたら報告していただけると嬉しいです. 欲しい記事に投稿するか,私クロメルの紹介ページのアドレスからメールでお願いします.

展開法

行列 AVB の三連続積を考えます. 簡単のため三次正方行列の場合を考えます.

まず,行列 AB の見方を変えます.

A=\begin{pmatrix}a_{11} & a_{12} & a_{13} \\a_{21} & a_{22} & a_{23} \\a_{31} & a_{32} & a_{33} \end{pmatrix}=\begin{pmatrix} \bm{a}_1^T &\bm{a}_2^T &\bm{a}_3^T \end{pmatrix}
B=\begin{pmatrix}b_{11} & b_{12} & b_{13} \\b_{21} & b_{22} & b_{23} \\b_{31} & b_{32} & b_{33} \end{pmatrix}=\begin{pmatrix}\bm{b}_1 \\ \bm{b}_2 \\ \bm{b}_3\end{pmatrix}

ここで行列を,ベクトルの集まりとして見方を変えました. \bm{P} ^T は列ベクトル, \bm{Q} は行ベクトルを表します. これで準備は完了です.三連続積を展開しましょう!

AVB =\begin{pmatrix} \bm{a}_1^T &\bm{a}_2^T &\bm{a}_3^T \end{pmatrix}\begin{pmatrix}v_{11} & v_{12} & v_{13} \\v_{21} & v_{22} & v_{23} \\v_{31} & v_{32} & v_{33} \end{pmatrix}\begin{pmatrix}\bm{b}_1 \\ \bm{b}_2 \\ \bm{b}_3\end{pmatrix}=\sum_{i , j}v_{ij} \bm{a}_{i}  \bm{b}_{j}

ここで上の式変形の最後の行に使ったダイアドを知らない方の為に説明しておくと, \bm{a}_{i} \bm{b}_{j} は, ダイアドと呼ばれるもので,

\bm{a}_{i} \bm{b}_{j} = \begin{pmatrix} a_{1i}  \\  a_{2i}  \\  a_{3i} \end{pmatrix}\begin{pmatrix} b_{j1} &b_{j2} &b_{j3} \end{pmatrix}=\begin{pmatrix}a_{1i} b_{j1} & a_{1i} b_{j2} & a_{1i} b_{j3} \\a_{2i} b_{j1} & a_{2i} b_{j2} & a_{2i} b_{j3} \\a_{3i} b_{j1} & a_{3i} b_{j2} & a_{3i} b_{j3}  \end{pmatrix}

という積です.

この記法は,もっと高次の場合も同様に,四次なら16個のダイアドの和,五次なら25個のダイアドの和となっていくのです.

ひとつ注意しておくと,これは展開する際に組になるベクトルと行列要素の位置を 実際に一度は展開してみて確認しておくことをお勧めします. 真ん中の V の非対角項を含む組がちとやっかいです.

追記:真ん中の V を単位行列 I に変えてもおもしろいですね.

AB=AIB =\begin{pmatrix} \bm{a}_1^T &\bm{a}_2^T &\bm{a}_3^T \end{pmatrix}\begin{pmatrix}1 & 0 & 0 \\0 & 1 & 0 \\0 & 0 & 1  \end{pmatrix}\begin{pmatrix}\bm{b}_1 \\ \bm{b}_2 \\ \bm{b}_3\end{pmatrix}=\sum_{i} \bm{a}_{i} \bm{b}_{i}