Gramian matrix

In linear algebra, the Gram matrix (or Gramian matrix, Gramian) of a set of vectors ${\displaystyle v_{1},\dots ,v_{n}}$ in an inner product space is the Hermitian matrix of inner products, whose entries are given by ${\displaystyle G_{ij}=\langle v_{i},v_{j}\rangle }$.[1]

An important application is to compute linear independence: a set of vectors are linearly independent if and only if the Gram determinant (the determinant of the Gram matrix) is non-zero.

It is named after Jørgen Pedersen Gram.

Examples

For finite-dimensional real vectors in ${\displaystyle \mathbb {R} ^{n}}$ with the usual Euclidean dot product, the Gram matrix is simply ${\displaystyle G=V^{\mathrm {T} }V}$, where ${\displaystyle V}$ is a matrix whose columns are the vectors ${\displaystyle v_{k}}$. For complex vectors in ${\displaystyle \mathbb {C} ^{n}}$ , ${\displaystyle G=V^{H}V}$, where ${\displaystyle V^{H}}$ is the conjugate transpose of ${\displaystyle V}$.

Given square-integrable functions ${\displaystyle \{\ell _{i}(\cdot ),\,i=1,\dots ,n\}}$ on the interval ${\displaystyle [t_{0},t_{f}]}$, the Gram matrix ${\displaystyle G=[G_{ij}]}$ is:

${\displaystyle G_{ij}=\int _{t_{0}}^{t_{f}}\ell _{i}(\tau ){\bar {\ell _{j}}}(\tau )\,d\tau .}$

For any bilinear form ${\displaystyle B}$ on a finite-dimensional vector space over any field we can define a Gram matrix ${\displaystyle G}$ attached to a set of vectors ${\displaystyle v_{1},\dots ,v_{n}}$ by ${\displaystyle G_{ij}=B(v_{i},v_{j})}$. The matrix will be symmetric if the bilinear form ${\displaystyle B}$ is symmetric.

Applications

• In Riemannian geometry, given an embedded ${\displaystyle k}$-dimensional Riemannian manifold ${\displaystyle M\subset \mathbb {R} ^{n}}$ and a coordinate chart ${\displaystyle \phi :U\to M}$ for ${\displaystyle (x_{1},\ldots ,x_{k})\in U\subset \mathbb {R} ^{k}}$, the volume form ${\displaystyle \omega }$ on ${\displaystyle M}$ induced by the embedding may be computed using the Gramian of the coordinate tangent vectors:

${\displaystyle \omega ={\sqrt {\det G}}\ dx_{1}\cdots dx_{k},\quad G=\left[\left\langle {\tfrac {\partial \phi }{\partial x_{i}}},{\tfrac {\partial \phi }{\partial x_{j}}}\right\rangle \right].}$

This generalizes the classical surface integral of a parametrized surface ${\displaystyle \phi :U\to S\subset \mathbb {R} ^{3}}$ for ${\displaystyle (x,y)\in U\subset \mathbb {R} ^{2}}$:

${\displaystyle \int _{S}f\ dA\ =\ \iint _{U}f(\phi (x,y))\,\left|{\tfrac {\partial \phi }{\partial x}}\,{\times }\,{\tfrac {\partial \phi }{\partial y}}\right|\,dx\,dy.}$

• If the vectors are centered random variables, the Gramian is approximately proportional to the covariance matrix, with the scaling determined by the number of elements in the vector.
• In quantum chemistry, the Gram matrix of a set of basis vectors is the overlap matrix.
• In control theory (or more generally systems theory), the controllability Gramian and observability Gramian determine properties of a linear system.
• Gramian matrices arise in covariance structure model fitting (see e.g., Jamshidian and Bentler, 1993, Applied Psychological Measurement, Volume 18, pp. 79–94).
• In the finite element method, the Gram matrix arises from approximating a function from a finite dimensional space; the Gram matrix entries are then the inner products of the basis functions of the finite dimensional subspace.
• In machine learning, kernel functions are often represented as Gram matrices.[2]
• Since the Gram matrix over the reals is a symmetric matrix, it is diagonalizable and its eigenvalues are non-negative. The diagonalization of the Gram matrix is the singular value decomposition.

Properties

Positive-semidefiniteness

The Gram matrix is symmetric in the case the real product is real-valued; it is Hermitian in the general, complex case by definition of an inner product.

The Gram matrix is positive semidefinite, and every positive semidefinite matrix is the Gramian matrix for some set of vectors. The fact that the Gramian matrix is positive-semidefinite can be seen from the following simple derivation:

${\displaystyle x^{\textsf {T}}\mathbf {G} x=\sum _{i,j}x_{i}x_{j}\left\langle v_{i},v_{j}\right\rangle =\sum _{i,j}\left\langle x_{i}v_{i},x_{j}v_{j}\right\rangle =\left\langle \sum _{i}x_{i}v_{i},\sum _{j}x_{j}v_{j}\right\rangle =\left\|\sum _{i}x_{i}v_{i}\right\|^{2}\geq 0.}$

The first equality follows from the definition of matrix multiplication, the second and third from the bi-linearity of the inner-product, and the last from the positive definiteness of the inner product. Note that this also shows that the Gramian matrix is positive definite if and only if the vectors ${\displaystyle v_{i}}$ are linearly independent (that is, ${\displaystyle \textstyle \sum _{i}x_{i}v_{i}\neq 0}$ for all ${\displaystyle x}$).[1]

Finding a vector realization

Given any positive semidefinite matrix ${\displaystyle M}$, one can decompose it as:

${\displaystyle M=B^{*}B}$,

where ${\displaystyle B^{*}}$ is the conjugate transpose of ${\displaystyle B}$ (or ${\displaystyle M=B^{\textsf {T}}B}$ in the real case). Here ${\displaystyle B}$ is a ${\displaystyle k\times n}$ matrix, where ${\displaystyle k}$ is the rank of ${\displaystyle M}$. Various ways to obtain such a decomposition include computing the Cholesky decomposition or taking the non-negative square root of ${\displaystyle M}$.

The columns ${\displaystyle b^{(1)},\dots ,b^{(n)}}$ of ${\displaystyle B}$ can be seen as n vectors in ${\displaystyle \mathbb {C} ^{k}}$ (or k-dimensional Euclidean space ${\displaystyle \mathbb {R} ^{k}}$, in the real case). Then

${\displaystyle M_{ij}=b^{(i)}\cdot b^{(j)}}$

where the dot product ${\displaystyle a\cdot b=\sum _{\ell =1}^{k}{\overline {a_{\ell }}}b_{\ell }}$ is the usual inner product on ${\displaystyle \mathbb {C} ^{k}}$.

Thus a Hermitian matrix ${\displaystyle M}$ is positive semidefinite if and only if it is the Gram matrix of some vectors ${\displaystyle b^{(1)},\dots ,b^{(n)}}$. Such vectors are called a vector realization of ${\displaystyle M}$. The infinite-dimensional analog of this statement is Mercer's theorem.

Uniqueness of vector realizations

If ${\displaystyle M}$ is the Gram matrix of vectors ${\displaystyle v_{1},\dots ,v_{n}}$ in ${\displaystyle \mathbb {R} ^{k}}$, then applying any rotation or reflection of ${\displaystyle \mathbb {R} ^{k}}$ (any orthogonal transformation, that is, any Euclidean isometry preserving 0) to the sequence of vectors results in the same Gram matrix. That is, for any ${\displaystyle k\times k}$ orthogonal matrix ${\displaystyle Q}$, the Gram matrix of ${\displaystyle Qv_{1},\dots ,Qv_{n}}$ is also ${\displaystyle M}$.

This is the only way in which two real vector realizations of ${\displaystyle M}$ can differ: the vectors ${\displaystyle v_{1},\dots ,v_{n}}$ are unique up to orthogonal transformations. In other words, the dot products ${\displaystyle v_{i}\cdot v_{j}}$ and ${\displaystyle w_{i}\cdot w_{j}}$ are equal if and only if some rigid transformation of ${\displaystyle \mathbb {R} ^{k}}$ transforms the vectors ${\displaystyle v_{1},\dots ,v_{n}}$ to ${\displaystyle w_{1},\dots ,w_{n}}$ and 0 to 0.

The same holds in the complex case, with unitary transformations in place of orthogonal ones. That is, if the Gram matrix of vectors ${\displaystyle v_{1},\dots ,v_{n}}$ is equal to the Gram matrix of vectors ${\displaystyle w_{1},\dots ,w_{n}}$ in ${\displaystyle \mathbb {C} ^{k}}$, then there is a unitary ${\displaystyle k\times k}$ matrix ${\displaystyle U}$ (meaning ${\displaystyle U^{*}U=I}$) such that ${\displaystyle v_{i}=Uw_{i}}$ for ${\displaystyle i=1,\dots ,n}$.[3]

Other properties

• The Gram matrix of any orthonormal basis is the identity matrix.
• The rank of the Gram matrix of vectors in ${\displaystyle \mathbb {R} ^{k}}$ or ${\displaystyle \mathbb {C} ^{k}}$ equals the dimension of the space spanned by these vectors.[1]

Gram determinant

The Gram determinant or Gramian is the determinant of the Gram matrix:

${\displaystyle G(x_{1},\dots ,x_{n})={\begin{vmatrix}\langle x_{1},x_{1}\rangle &\langle x_{1},x_{2}\rangle &\dots &\langle x_{1},x_{n}\rangle \\\langle x_{2},x_{1}\rangle &\langle x_{2},x_{2}\rangle &\dots &\langle x_{2},x_{n}\rangle \\\vdots &\vdots &\ddots &\vdots \\\langle x_{n},x_{1}\rangle &\langle x_{n},x_{2}\rangle &\dots &\langle x_{n},x_{n}\rangle \end{vmatrix}}.}$

If ${\displaystyle x_{1},\cdots ,x_{n}}$ are vectors in ${\displaystyle \mathbb {R} ^{n}}$, then it is the square of the n-dimensional volume of the parallelotope formed by the vectors. In particular, the vectors are linearly independent if and only if the parallelotope has nonzero n-dimensional volume, if and only if Gram determinant is nonzero, if and only if the Gram matrix is nonsingular.

The Gram determinant can also be expressed in terms of the exterior product of vectors by

${\displaystyle G(x_{1},\dots ,x_{n})=\|x_{1}\wedge \cdots \wedge x_{n}\|^{2}.}$

See also

• Controllability Gramian
• Observability Gramian

References

1. Horn & Johnson 2013, p. 441, p.441, Theorem 7.2.10
2. Lanckriet, G. R. G.; Cristianini, N.; Bartlett, P.; Ghaoui, L. E.; Jordan, M. I. (2004). "Learning the kernel matrix with semidefinite programming". Journal of Machine Learning Research. 5: 27–72 [p. 29].
3. Horn & Johnson (2013), p. 452, Theorem 7.3.11

• Horn, Roger A.; Johnson, Charles R. (2013). Matrix Analysis (2nd ed.). Cambridge University Press. ISBN 978-0-521-54823-6.CS1 maint: ref=harv (link)

External links

• "Gram matrix", Encyclopedia of Mathematics, EMS Press, 2001 [1994]
• Volumes of parallelograms by Frank Jones