Transpose(Redirected from Matrix transpose)
In linear algebra, the transpose of a matrix is an operator which flips a matrix over its diagonal, that is it switches the row and column indices of the matrix by producing another matrix denoted as AT (also written A′, Atr, tA or At). It is achieved by any one of the following equivalent actions:
- reflect A over its main diagonal (which runs from top-left to bottom-right) to obtain AT,
- write the rows of A as the columns of AT,
- write the columns of A as the rows of AT.
Formally, the i-th row, j-th column element of AT is the j-th row, i-th column element of A:
If A is an m × n matrix, then AT is an n × m matrix.
For matrices A, B and scalar c we have the following properties of transpose:
- The transpose respects addition.
- Note that the order of the factors reverses. From this one can deduce that a square matrix A is invertible if and only if AT is invertible, and in this case we have (A−1)T = (AT)−1. By induction this result extends to the general case of multiple matrices, where we find that (A1A2...Ak−1Ak)T = AkTAk−1T…A2TA1T.
- The determinant of a square matrix is the same as the determinant of its transpose.
- The dot product of two column vectors a and b can be computed as the single entry of the matrix product:
- If A has only real entries, then ATA is a positive-semidefinite matrix.
- The transpose of an invertible matrix is also invertible, and its inverse is the transpose of the inverse of the original matrix. The notation A−T is sometimes used to represent either of these equivalent expressions.
- If A is a square matrix, then its eigenvalues are equal to the eigenvalues of its transpose, since they share the same characteristic polynomial.
Special transpose matricesEdit
A square matrix whose transpose is equal to itself is called a symmetric matrix; that is, A is symmetric if
A square matrix whose transpose is equal to its negative is called a skew-symmetric matrix; that is, A is skew-symmetric if
A square complex matrix whose transpose is equal to the matrix with every entry replaced by its complex conjugate (denoted here with an overline) is called a Hermitian matrix (equivalent to the matrix being equal to its conjugate transpose); that is, A is Hermitian if
A square matrix whose transpose is equal to its inverse is called an orthogonal matrix; that is, A is orthogonal if
A square complex matrix whose transpose is equal to its conjugate inverse is called a unitary matrix; that is, A is unitary if
Transpose of a linear mapEdit
The transpose may be defined more generally:
Equivalently, the transpose tf is defined by the relation
The definition of the transpose may be seen to be independent of any bilinear form on the vector spaces, unlike the adjoint (below).
If the matrix A describes a linear map with respect to bases of V and W, then the matrix AT describes the transpose of that linear map with respect to the dual bases.
Transpose of a bilinear formEdit
Every linear map to the dual space f : V → V∗ defines a bilinear form B : V × V → F, with the relation B(v, w) = f(v)(w). By defining the transpose of this bilinear form as the bilinear form tB defined by the transpose tf : V∗∗ → V∗ i.e. tB(w, v) = tf(Ψ(w))(v), we find that B(v, w) = tB(w, v). Here, Ψ is the natural homomorphism V → V∗∗ into the double dual.
These bilinear forms define an isomorphism between V and V∗, and between W and W∗, resulting in an isomorphism between the transpose and adjoint of f. The matrix of the adjoint of a map is the transposed matrix only if the bases are orthonormal with respect to their bilinear forms. In this context, many authors use the term transpose to refer to the adjoint as defined here.
The adjoint allows us to consider whether g : W → V is equal to f −1 : W → V. In particular, this allows the orthogonal group over a vector space V with a quadratic form to be defined without reference to matrices (nor the components thereof) as the set of all linear maps V → V for which the adjoint equals the inverse.
Over a complex vector space, one often works with sesquilinear forms (conjugate-linear in one argument) instead of bilinear forms. The Hermitian adjoint of a map between such spaces is defined similarly, and the matrix of the Hermitian adjoint is given by the conjugate transpose matrix if the bases are orthonormal.
Implementation of matrix transposition on computersEdit
On a computer, one can often avoid explicitly transposing a matrix in memory by simply accessing the same data in a different order. For example, software libraries for linear algebra, such as BLAS, typically provide options to specify that certain matrices are to be interpreted in transposed order to avoid the necessity of data movement.
However, there remain a number of circumstances in which it is necessary or desirable to physically reorder a matrix in memory to its transposed ordering. For example, with a matrix stored in row-major order, the rows of the matrix are contiguous in memory and the columns are discontiguous. If repeated operations need to be performed on the columns, for example in a fast Fourier transform algorithm, transposing the matrix in memory (to make the columns contiguous) may improve performance by increasing memory locality.
Ideally, one might hope to transpose a matrix with minimal additional storage. This leads to the problem of transposing an n × m matrix in-place, with O(1) additional storage or at most storage much less than mn. For n ≠ m, this involves a complicated permutation of the data elements that is non-trivial to implement in-place. Therefore, efficient in-place matrix transposition has been the subject of numerous research publications in computer science, starting in the late 1950s, and several algorithms have been developed.
- Arthur Cayley (1858) "A memoir on the theory of matrices", Philosophical Transactions of the Royal Society of London, 148 : 17–37. The transpose (or "transposition") is defined on page 31.
- Bourbaki, "II §2.5", Algebra I