Block matrix

From Wikipedia, the free encyclopedia

In the mathematical discipline of matrix theory, a block matrix or a partitioned matrix is a partition of a matrix into rectangular smaller matrices called blocks. Looking at it another way, the matrix is written in terms of smaller matrices written side-by-side. A block matrix must conform to a consistent way of splitting up the rows, and the columns: we group the rows into some adjacent 'bunches', and the columns likewise. The partition is into the rectangles described by one bunch of adjacent rows crossing one bunch of adjacent columns. In other words, the matrix is split up by some horizontal and vertical lines that go all the way across.

Contents

[edit] Example

The matrix

P = \begin{bmatrix} 1 & 1 & 2 & 2\\ 1 & 1 & 2 & 2\\ 3 & 3 & 4 & 4\\ 3 & 3 & 4 & 4\end{bmatrix}

can be partitioned into 4 2×2 blocks

P_{11} = \begin{bmatrix} 1 & 1 \\ 1 & 1 \end{bmatrix},   P_{12} = \begin{bmatrix} 2 & 2\\ 2 & 2\end{bmatrix},  P_{21} = \begin{bmatrix} 3 & 3 \\ 3 & 3 \end{bmatrix},   P_{22} = \begin{bmatrix} 4 & 4\\ 4 & 4\end{bmatrix}.

The partitioned matrix can then be written as

P_{\mathrm{partitioned}} = \begin{bmatrix} P_{11} & P_{12}\\ P_{21} & P_{22}\end{bmatrix}.

[edit] Block diagonal matrices

A block diagonal matrix is a block matrix which is a square matrix, and having main diagonal blocks square matrices, such that the off-diagonal blocks are zero matrices. A block diagonal matrix A has the form

\mathbf{A} = \begin{bmatrix}  A_{1} & 0 & \cdots & 0 \\ 0 & A_{2} & \cdots &  0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & A_{n}  \end{bmatrix}

where Ak is a square matrix; in other words, it is the direct sum of A1, …, An. It can also be indicated as A_1\oplus A_2\oplus\ldots\oplus A_n or \mbox{diag}\left(A_1, A_2,\ldots, A_n\right) (the latter being the same formalism used for a diagonal matrix). Any square matrix can trivially be considered a block diagonal matrix with only one block.

For the determinant and trace, the following properties hold

\operatorname{det} \mathbf{A} = \operatorname{det} A_1 \cdots \operatorname{det} A_n,
\operatorname{trace} \mathbf{A} = \operatorname{trace} A_1 +\cdots +\operatorname{trace} A_n.

[edit] Block tridiagonal matrices

A block tridiagonal matrix is a another special block matrix, which is just like the block diagonal matrix a square matrix, having square matrices (blocks) in the lower diagonal, main diagonal and upper diagonal, with all other blocks being zero matrices. It is essentially a tridiagonal matrix but has submatrices in places of scalars. A block tridiagonal matrix A has the form

\mathbf{A} = \begin{bmatrix} B_{1}  & C_{1}  &         &         & \cdots  &         & 0 \\ A_{2}  & B_{2}  & C_{2}   &         &         &         & \\        & \ddots & \ddots  & \ddots  &         &         & \vdots \\        &        & A_{k}   & B_{k}   & C_{k}   &         & \\ \vdots &        &         & \ddots  & \ddots  & \ddots  & \\        &        &         &         & A_{n-1} & B_{n-1} & C_{n-1}   \\ 0      &        & \cdots  &         &         & A_{n}   & B_{n} \end{bmatrix}

where Ak, Bk and Ck are square sub-matrices of the lower, main and upper diagonal respectively.

Block tridiagonal matricies are often encountered in numerical solutions of engineering problems (e.g. computational fluid dynamics). Optimized numerical methods for LU factorization are available and hence efficient solution algorithms for equation systems with a block tridiagonal matrix as coefficient matrix. The Thomas algorithm, used for efficient solution of equation systems involving a tridiagonal matrix can also be applied using matrix operations to block tridiagonal matrices (see also Block LU decomposition).

[edit] Direct sum

For any arbitrary matrices A (of size m × n) and B (of size p × q), we have the direct sum of A and B, denoted by A \oplus B and defined as

A \oplus B =   \begin{bmatrix}      a_{11} & \cdots & a_{1n} &      0 & \cdots &      0 \\      \vdots & \cdots & \vdots & \vdots & \cdots & \vdots \\     a_{m 1} & \cdots & a_{mn} &      0 & \cdots &      0 \\           0 & \cdots &      0 & b_{11} & \cdots &  b_{1q} \\      \vdots & \cdots & \vdots & \vdots & \cdots & \vdots \\           0 & \cdots &      0 & b_{p1} & \cdots &  b_{pq}    \end{bmatrix}

For instance,

\begin{bmatrix}     1 & 3 & 2 \\     2 & 3 & 1   \end{bmatrix} \oplus   \begin{bmatrix}     1 & 6 \\     0 & 1   \end{bmatrix} =   \begin{bmatrix}     1 & 3 & 2 & 0 & 0 \\     2 & 3 & 1 & 0 & 0 \\     0 & 0 & 0 & 1 & 6 \\     0 & 0 & 0 & 0 & 1   \end{bmatrix}

Note that any element in the direct sum of two vector spaces of matrices could be represented as a direct sum of two matrices.

[edit] Application

In linear algebra terms, the use of a block matrix corresponds to having a linear mapping thought of in terms of corresponding 'bunches' of basis vectors. That again matches the idea of having distinguished direct sum decompositions of the domain and range. It is always particularly significant if a block is the zero matrix; that carries the information that a summand maps into a sub-sum.

Given the interpretation via linear mappings and direct sums, there is a special type of block matrix that occurs for square matrices (the case m = n). For those we can assume an interpretation as an endomorphism of an n-dimensional space V; the block structure in which the bunching of rows and columns is the same is of importance because it corresponds to having a single direct sum decomposition on V (rather than two). In that case, for example, the diagonal blocks in the obvious sense are all square. This type of structure is required to describe the Jordan normal form.

This technique is used to cut down calculations of matrices, column-row expansions, and many computer science applications, including VLSI chip design. An example is the Strassen algorithm for fast matrix multiplication.