Neha Patil (Editor)

Vandermonde matrix

Updated on
Edit
Like
Comment
Share on FacebookTweet on TwitterShare on LinkedInShare on Reddit

In linear algebra, a Vandermonde matrix, named after Alexandre-Théophile Vandermonde, is a matrix with the terms of a geometric progression in each row, i.e., an m × n matrix

Contents

V = [ 1 α 1 α 1 2 … α 1 n − 1 1 α 2 α 2 2 … α 2 n − 1 1 α 3 α 3 2 … α 3 n − 1 ⋮ ⋮ ⋮ ⋱ ⋮ 1 α m α m 2 … α m n − 1 ] ,

or

V i , j = α i j − 1

for all indices i and j. (Some authors use the transpose of the above matrix.)

The determinant of a square Vandermonde matrix (where m = n) can be expressed as

det ( V ) = ∏ 1 ≤ i < j ≤ n ( α j − α i ) .

This is called the Vandermonde determinant or Vandermonde polynomial. If all the numbers α i are distinct, then it is non-zero.

The Vandermonde determinant is sometimes called the discriminant, although many sources, including this article, refer to the discriminant as the square of this determinant. Note that the Vandermonde determinant is alternating in the entries, meaning that permuting the α i by an odd permutation changes the sign, while permuting them by an even permutation does not change the value of the determinant. It thus depends on the order, while its square (the discriminant) does not depend on the order.

When two or more αi are equal, the corresponding polynomial interpolation problem (see below) is underdetermined. In that case one may use a generalization called confluent Vandermonde matrices, which makes the matrix non-singular while retaining most properties. If αi = αi + 1 = ... = αi+k and αi ≠ αi − 1, then the (i + k)th row is given by

V i + k , j = { 0 , if  j ≤ k ; ( j − 1 ) ! ( j − k − 1 ) ! α i j − k − 1 , if  j > k .

The above formula for confluent Vandermonde matrices can be readily derived by letting two parameters α i and α j go arbitrarily close to each other. The difference vector between the rows corresponding to α i and α j scaled to a constant yields the above equation (for k = 1). Similarly, the cases k > 1 are obtained by higher order differences. Consequently, the confluent rows are derivatives of the original Vandermonde row.

Properties

In the case of a square Vandermonde matrix, the Leibniz formula for the determinant gives

det ( V ) = ∑ σ ∈ S n sgn ⁡ ( σ ) ∏ i = 1 n α i σ ( i ) − 1 ,

where Sn denotes the set of permutations of { 1 , … , n } , and sgn ⁡ ( σ ) denotes the signature of the permutation σ. This determinant factors as

∑ σ ∈ S n sgn ⁡ ( σ ) ∏ i = 1 n α i σ ( i ) − 1 = ∏ 1 ≤ i < j ≤ n ( α j − α i ) .

Each of these factors must divide the determinant, because the latter is an alternating polynomial in the n variables. It also follows that the Vandermonde determinant divides any other alternating polynomial; the quotient will be a symmetric polynomial.

If m ≤ n, then the matrix V has maximum rank (m) if and only if all αi are distinct. A square Vandermonde matrix is thus invertible if and only if the αi are distinct; an explicit formula for the inverse is known.

Applications

The Vandermonde matrix evaluates a polynomial at a set of points; formally, it transforms coefficients of a polynomial a 0 + a 1 x + a 2 x 2 + ⋯ + a n − 1 x n − 1 to the values the polynomial takes at the points α i . The non-vanishing of the Vandermonde determinant for distinct points α i shows that, for distinct points, the map from coefficients to values at those points is a one-to-one correspondence, and thus that the polynomial interpolation problem is solvable with unique solution; this result is called the unisolvence theorem.

They are thus useful in polynomial interpolation, since solving the system of linear equations Vu = y for u with V an m × n Vandermonde matrix is equivalent to finding the coefficients uj of the polynomial(s)

P ( x ) = ∑ j = 0 n − 1 u j x j

of degree ≤ n − 1 which has (have) the property

P ( α i ) = y i for  i = 1 , … , m .

The Vandermonde matrix can be inverted in terms of Lagrange basis polynomials: each column is the coefficients of the Lagrange basis polynomial, with terms in increasing order going down. The resulting solution to the interpolation problem is called the Lagrange polynomial.

The Vandermonde determinant plays a central role in the Frobenius formula, which gives the character of conjugacy classes of representations of the symmetric group.

When the values α k range over powers of a finite field, then the determinant has a number of interesting properties: for example, in proving the properties of a BCH code.

Confluent Vandermonde matrices are used in Hermite interpolation.

A commonly known special Vandermonde matrix is the discrete Fourier transform matrix (DFT matrix), where the numbers αi are chosen to be the m different mth roots of unity.

The Vandermonde matrix diagonalizes the companion matrix.

The Vandermonde matrix is used in some forms of Reed–Solomon error correction codes.

References

Vandermonde matrix Wikipedia
(Text) CC BY-SA