The square root of a matrix, or matrix root, is a linear algebra term that generalizes the concept of the square root of a real number . A square root of a square matrix is a matrix that, when multiplied by itself, results in the output matrix . A unique square root can be defined for symmetric positive semidefinite matrices. In general, however, a square root does not need to exist, nor does it need to be unique if it does exist.
Square root of a positive semi-definite matrix
definition
For a symmetric positive semidefinite matrix , a symmetric positive semidefinite matrix is called the square root or, for short, the root of if
applies. The square root of is uniquely determined and is denoted by.
presentation
The square root of is given as follows. According to the spectral theorem, there is an orthogonal matrix
with pairwise orthonormal eigenvectors of as columns and a diagonal matrix
with the eigenvalues associated with these eigenvectors on the diagonal , so that
applies. The square root of then results in
where the diagonal matrix
has the square roots of the eigenvalues of on the diagonal. Since the eigenvalues of a positive semidefinite matrix are always real and nonnegative, their square roots can also be chosen to be real and nonnegative.
example
The matrix
has the eigenvalues and and form the corresponding orthonormal basis from eigenvectors. So it applies
summarized to
and thus
properties
The square of the matrix is the matrix
The matrix is symmetrical:
The matrix is positive semidefinite (the displacement property of the standard scalar product is used ):
for all where . If positive is definite, then positive is definite.
Square roots of arbitrary matrices
definition
The root of a square matrix is any matrix that, when multiplied by itself, results in:
There are also sources in which a root of is mentioned when applies.
A root of is also written . In this notation, however, it is unclear which root is meant, since several can exist.
Number of existing roots
As with the square root of real or complex numbers, the square root of a matrix is generally ambiguous. Is about a root from then too
Unlike the root of a complex number, matrices can have more than two roots.
For example, matrices, whose characteristic polynomial breaks down into different linear factors in pairs, have up to different roots.
There are even matrices with an infinite number of roots. The identity matrix , for example, has a root for every complex number .
There are also matrices for which there is no root at all: An example is
Geometric interpretation of roots
Looking at the array as a linear transformation , that is, as a mapping between vector spaces through which a vector , a vector is assigned, then a root is a transformation that is twice in succession must perform in order to convert.
Example:
be the two-dimensional rotation matrix with the angle
Then every rotation matrix belonging to an angle with an integer is a square root of For with the first multiplication of a vector with a rotation by half the angle and with the second multiplication again.
Calculating a Root
One can easily determine the roots of a matrix of size if it is a diagonal matrix or can at least be converted into a diagonal form (see diagonalization ).
Case 1: diagonal matrix
If the diagonal entries of a diagonal matrix are different in pairs, then all the roots of the diagonal matrix can be determined simply by determining a root from each entry on the main diagonal . If you denote the diagonal entries of as usual with , you get the roots of the matrices
For each of the diagonal elements, the sign can be chosen arbitrarily, so that one obtains different roots in pairs, if all diagonal entries are different from zero. If a diagonal entry is zero, then one obtains roots that are different in pairs. Since the matrix can also have negative values on the diagonal, the roots can also contain complex numbers.
It should be noted that there can be further roots if the diagonal entries are not different in pairs. However, these are then not diagonal matrices. For example, the identity matrix has an infinite number of roots, as already explained above . Diagonal matrices with negative diagonal entries can in this case also have real roots. For example:
Case 2: Diagonalizable matrix
If the matrix can be diagonalized, the roots of can be found in the following way :
First one determines an invertible matrix and a diagonal matrix so that . The matrix then has eigenvectors of the matrix as columns and the matrix has the associated eigenvalues as diagonal entries .
If now is a root of then is a root of the matrix , because the following applies:
Since is a diagonal matrix, one obtains possible roots as in case 1. Here, too, it should be noted that some eigenvalues of the diagonal matrix can be negative, the roots of which are then complex. If the matrix pairwise distinct eigenvalues, are also obtained as in case 1 or different solutions.
Case 3: Non-diagonalizable matrix
If the matrix cannot be diagonalized, the method shown cannot be used to calculate a root. However, this does not mean that it does not have a root: For example, the shear matrix cannot be diagonalized, but it does have the root
If we allow complex numbers in arithmetic, every matrix can be transformed into Jordanian normal form , even if it cannot be diagonalized.
One determines matrices their inverse and with where has the following block diagonal form:
They are Jordan blocks of shape
A square root of is calculated according to
The root of is to be taken from each Jordan block individually.
If true, the power of a Jordan block is through
given with where the -th derivative of the power function is. This results explicitly and where the size of the Jordan block is denoted by (in the illustration ), the subdiagonals by ( is the diagonal) and the gamma function by . Set for the square root .
For example, for
If and at the same time , the root of the Jordan block does not exist.
Outside the Jordan blocks there are all zeros.
If so, the number has two roots, so you get two different roots for each Jordan block this way . Combination creates roots, where the number of Jordan blocks denotes.
This technique generally gives only some, and not all, of the square roots of a matrix.
See also
Web links
Individual evidence
-
↑ a b Christian Kanzow: Numerics of linear systems of equations. Direct and iterative procedures. Springer-Verlag, Berlin / Heidelberg / New York 2005, pp. 13–15.