next up previous contents
Next: Symmetry Up: Linear Algebra Previous: Second neighbor influences   Contents

Block Matrices

Matrices sometimes have a structure in which some submatrix is repeated over and over, usually reappearing in the form of different multiples. For example, the product of a column by a row is a rectangular matrix, in which all the rows are proportional, just as the columns are likewise proportional. When the repetition of one matrix is guided by another, the result is sometimes called a KroneckerKronecker product product of the two matrices, and at other times a tensortensor product product. For example, given matrices $P$ and $Q$,

\begin{displaymath}P = \left[ \begin{array}{cc} a & b \\ c & d \end{array} \righ...
... \left[ \begin{array}{cc} A & B \\ C & D \end{array} \right],
\end{displaymath}

their Kronecker product $P \otimes Q$ would be defined by

\begin{eqnarray*}
P \otimes Q & = & \left[ \begin{array}{cc}
a Q & b Q \\ c Q...
...cA & cB & dA & dB \\
cC & cD & dC & dD
\end{array} \right].
\end{eqnarray*}



In complete generality, give $P \otimes Q$ four indices in the form of two pairs, and set

\begin{eqnarray*}[P \otimes Q]_{ij,kl} & = & P_{ik} Q_{jl}
\end{eqnarray*}



These indices are to be run out in lexicographic order, which means that $P \otimes Q$ is not the same thing as $Q \otimes P$, although it is equivalent by a change of basis which permutes the indices.

\begin{eqnarray*}
\left[ \begin{array}{cccc}
aA & aB & bA & bB \\
aC & aD &...
...
aC & bC & aD & bD \\
cC & dC & cD & dD
\end{array} \right]
\end{eqnarray*}



Calculation with Kronecker products is facilitated by verifying all the algebraic rules, such as distributivity or associativity, which they obey, thereafter using them automatically. Let $I$ stand for a unit matrix, $O$ zero matrices of the appropriate size and shape, $P$, $Q$ and $R$ generic matrices. As needed, $\alpha$ and $\beta$ could be scalar factors. Then

\begin{eqnarray*}
I \otimes I & = & I \\
O \otimes P & = & O \\
P \otimes O...
..._2) \\
(P \otimes Q) \otimes R & = & P \otimes (Q \otimes R).
\end{eqnarray*}



Apparently a Kronecker product resembles a cartesian product, in the sense that the factors combine fairly independently when it comes to mingling the Kronecker product with other operations such as linear combination or matrix multiplication. Some further results are useful; if the factors are invertible,

\begin{eqnarray*}
P^{-1} \otimes Q^{-1} & = & (P \otimes Q)^{-1},
\end{eqnarray*}



wherein it should be noted that the order of the factors is preserved.

If there are eigenvectors $PX=\lambda X$ and $QY=\mu Y$,

\begin{eqnarray*}
(P \otimes Q) (X \otimes Y) & = & (\lambda \mu) (X \otimes Y),
\end{eqnarray*}



This means that the Kronecker product has all possible products of eigenvalues associated with all possible Kronecker products of eigenvectors, due regard being given to the possibility of Jordan normal forms, and with similar results whenever they occur. Beyond that, under the same assumptions,

\begin{eqnarray*}
(P \otimes I + I \otimes Q) (X \otimes Y) & = &
(\lambda + ...
...\otimes Q) (X \otimes Y) & = &
(\lambda - \mu) (X \otimes Y),
\end{eqnarray*}



which provides a matrix with all the possible sums, or all the possible differences, of the eigenvalues of two given matrices as its own eigenvalues. Checking whether the second combination were singular would tell immediately whether a pair of matrices had common eigenvalues.


next up previous contents
Next: Symmetry Up: Linear Algebra Previous: Second neighbor influences   Contents
Pedro Hernandez 2004-02-28