Band matrix

Band matrix

In mathematics, particularly matrix theory, a band matrix is a sparse matrix whose non-zero entries are confined to a diagonal band, comprising the main diagonal and zero or more diagonals on either side.

Contents

Matrix bandwidth

Formally, consider an n×n matrix A=(ai,j ). If all matrix elements are zero outside a diagonally bordered band whose range is determined by constants k1 and k2:

a_{i,j}=0 \quad\mbox{if}\quad j<i-k_1 \quad\mbox{ or }\quad j>i+k_2; \quad k_1, k_2 \ge 0.\,

then the quantities k1 and k2 are called the left and right half-bandwidth, respectively (Golub & Van Loan 1996, §1.2.1). The bandwidth of the matrix is k1 + k2 + 1 (in other words, it is the smallest number of adjacent diagonals to which the non-zero elements are confined).

A matrix is called a band matrix or banded matrix if its bandwidth is reasonably small.

A band matrix with k1 = k2 = 0 is a diagonal matrix; a band matrix with k1 = k2 = 1 is a tridiagonal matrix; when k1 = k2 = 2 one has a pentadiagonal matrix and so on. If one puts k1 = 0, k2 = n−1, one obtains the definition of an upper triangular matrix; similarly, for k1 = n−1, k2 = 0 one obtains a lower triangular matrix.

Applications

In numerical analysis, matrices from finite element or finite difference problems are often banded. Such matrices can be viewed as descriptions of the coupling between the problem variables; the bandedness corresponds to the fact that variables are not coupled over arbitrarily large distances. Such matrices can be further divided - for instance, banded matrices exist where every element in the band is nonzero. These often arise when discretising one-dimensional problems.[citation needed]

Problems in higher dimensions also lead to banded matrices, in which case the band itself also tends to be sparse. For instance, a partial differential equation on a square domain (using central differences) will yield a matrix with a half-bandwidth equal to the square root of the matrix dimension, but inside the band only 5 diagonals are nonzero. Unfortunately, applying Gaussian elimination (or equivalently an LU decomposition) to such a matrix results in the band being filled in by many non-zero elements.

Band storage

Band matrices are usually stored by storing the diagonals in the band; the rest is implicitly zero.

For example, a tridiagonal matrix has bandwidth 3. The 6-by-6 matrix


\begin{bmatrix}
 B_{11} & B_{12} & 0      & \cdots & \cdots & 0 \\
 B_{21} & B_{22} & B_{23} & \ddots & \ddots & \vdots \\
  0     & B_{32} & B_{33} & B_{34} & \ddots & \vdots \\
 \vdots & \ddots & B_{43} & B_{44} & B_{45} & 0 \\
 \vdots & \ddots & \ddots & B_{54} & B_{55} & B_{56} \\
 0      & \cdots & \cdots & 0      & B_{65} & B_{66}
\end{bmatrix}

is stored as the 6-by-3 matrix


\begin{bmatrix}
 0 & B_{11} & B_{12}\\
 B_{21} & B_{22} & B_{23} \\
 B_{32} & B_{33} & B_{34} \\
 B_{43} & B_{44} & B_{45} \\
 B_{54} & B_{55} & B_{56} \\
 B_{65} & B_{66} & 0
\end{bmatrix}.

A further saving is possible when the matrix is symmetric. For example, consider a symmetric 6-by-6 matrix with a right bandwidth of 2:


\begin{bmatrix}
 A_{11} & A_{12} & A_{13} &   0  & \cdots & 0 \\
      & A_{22} & A_{23} & A_{24} & \ddots & \vdots \\
      &        & A_{33} & A_{34} & A_{35} & 0 \\
      &        &        & A_{44} & A_{45} & A_{46} \\
      & sym    &        &        & A_{55} & A_{56} \\
      &        &        &        &        & A_{66}
\end{bmatrix}.

This matrix is stored as the 6-by-3 matrix:


\begin{bmatrix}
 A_{11} & A_{12} & A_{13} \\
 A_{22} & A_{23} & A_{24} \\
 A_{33} & A_{34} & A_{35} \\
 A_{44} & A_{45} & A_{46} \\
 A_{55} & A_{56} & 0 \\
 A_{66} & 0 & 0
\end{bmatrix}.

Band form of sparse matrices

From a computational point of view, working with band matrices is always preferential to working with similarly dimensioned square matrices. A band matrix can be likened in complexity to a rectangular matrix whose row dimension is equal to the bandwidth of the band matrix. Thus the work involved in performing operations such as multiplication falls significantly, often leading to huge savings in terms of calculation time and complexity.

As sparse matrices lend themselves to more efficient computation than dense matrices, as well as in more efficient utilization of computer storage, there has been much research focused on finding ways to minimise the bandwidth (or directly minimise the fill in) by applying permutations to the matrix, or other such equivalence or similarity transformations.[citation needed]

The Cuthill–McKee algorithm can be used to reduce the bandwidth of a sparse symmetric matrix. There are, however, matrices for which the reverse Cuthill–McKee algorithm performs better. There are many other methods in use.

The problem of finding a representation of a matrix with minimal bandwidth by means of permutations of rows and columns is NP-hard.[1]

Examples and special cases

The following are special cases of band matrices:

The inverses of Lehmer matrices are constant tridiagonal matrices, and are thus band matrices.

See also

Graph bandwidth

References

  1. ^ "Coping with the NP-Hardness of the Graph Bandwidth Problem", Uriel Feige, Lecture Notes in Computer Science, Volume 1851, 2000, pp. 129-145, doi:10.1007/3-540-44985-X_2

External links


Wikimedia Foundation. 2010.

Игры ⚽ Поможем написать реферат

Look at other dictionaries:

  • Matrix-Metalloproteasen — (MMP) sind eine Gruppe verwandter Enzyme, die die Spaltung von Peptidbindungen in Proteinen katalysieren (so genannte Peptidasen). Sie kommen in allen Lebewesen außer Protozoen vor. In Eukaryoten sind sie außerhalb von Zellen in der… …   Deutsch Wikipedia

  • Matrix Reloaded — Données clés Titre original The Matrix Reloaded Réalisation Andy et Larry Wachowski Scénario Andy et Larry Wachowski Acteurs principaux Keanu Reeves Laurence Fishburne Carrie Anne Moss …   Wikipédia en Français

  • Matrix reloaded — Titre original The Matrix Reloaded Réalisation Andy et Larry Wachowski Acteurs principaux Keanu Reeves Laurence Fishburne Carrie Anne Moss Hugo Weaving Daniel Bernhardt Monica Bellucci Jada Pinkett Smith Lambert Wilson Anthony Wong …   Wikipédia en Français

  • Matrix Games — Тип Частная компания Год основания 2001 Расположение Статен …   Википедия

  • Matrix (jazz group) — Matrix is a jazz fusion group from Appleton, Wisconsin, that started in 1974, noted for tight brass ensemble lines and complex musical themes. It made its biggest impact on the music scene in the 1970s, including an appearance at the Monterey… …   Wikipedia

  • Matrix Quest — The Matrix Quest is the official collective name for issues #62 through #66 in Marvel s Transformers comics, written by UK writer Simon Furman. Each chapter of the story pays homage to a classic movie or book. Each issue of Matrix Quest was… …   Wikipedia

  • Matrix Games — For the Japanese company, see Matrix Software. Matrix Games Industry Strategy and Wargames for Personal Computers. Founded 1999 Headquarters Staten Island, New York Key …   Wikipedia

  • Matrix — Contents 1 Science and mathematics 2 Technology 3 Arts and entertainment …   Wikipedia

  • matrix — 1. [NA] The formative portion of a tooth or a nail. 2. The intercellular substance of a tissue. 3. A surrounding substance within which something is contained or embedded, e.g., the fatty tissue in which blood vessel s or …   Medical dictionary

  • band — 1. Any appliance or part of an apparatus that encircles or binds a part of the body. SEE ALSO: zone. 2. Any ribbon shaped or cordlike anatomic structure that encircles or binds another structure or that connects two or more parts. See fascia,… …   Medical dictionary

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”