Álgebra Lineal
Tema 2. Matriz de Vandermonde
Grado en Ingeniería Informática
Doble Grado en Ingeniería Informática y Administración
de Empresas
AUTORES: J. S ALAS , A. T ORRENTE Y E.J.S. V ILLASEÑOR
Índice general
2. Matriz de Vandermonde 1
2.1. Cálculo del determinante de una matriz de Vandermonde cuadrada . . . . . 2
I
Tema 2
Matriz de Vandermonde
La matriz de Vandermonde de dimensión m × n viene dada por la expresión
1 α1 α21 ··· αn−1
1
1 α2 α22 ··· αn−1
2
An (α1 , . . . , αm ) =
.. .. .. ..
. . . .
1 αm α2m · · · αn−1
m
donde los αi con 1 6 i 6 m son números reales arbitrarios. El nombre de este tipo especial
de matrices hace referencia al matemático francés Alexandre-Théophile Vandermonde,
1735–1796. En algunos libros se define la matriz de Vandermonde como la transpuesta de
la anterior.
Este tipo de matrices aparece de manera natural cuando queremos ajustar m datos
experimentales (αi , yi ) a un polinomio en x de grado n − 1
X
n−1
y= ck xk = c0 + c1 x + · · · + cn−1 xn−1
k=0
usando el método de mínimos cuadrados que se estudiará en el Tema 13 del curso.
1
2.1. Cálculo del determinante de una matriz de Vandermon-
de cuadrada
Supongamos ahora que n = m de manera que la matriz de Vandermonde es cuadrada
1 α1 α21 ··· αn−1
1
1 α2 α22 · · · αn−1
2
An (α1 , . . . , αn ) = .
.. .. .. ..
. . . .
1 αn α2n · · · αn−1
n
Un problema clásico consiste en evaluar su determinante
Vn (α1 , . . . , αn ) = det(An (α1 , . . . , αn )) .
De las propiedades de los determinantes es obvio que si hay al menos dos variables αi
que son iguales, entonces dicho determinante es cero:
Vn (α1 , . . . , αn ) = 0 , si αj = αk para algún par j 6= k.
Veamos algunos ejemplos sencillos:
Caso n = 2
Este caso es trivial:
1 α1
V2 (α1 , α2 ) = = α2 − α1 .
1 α
2
Obviamente si α1 = α2 = α, V2 (α, α) = 0.
2
Caso n = 3
Este caso es algo más complicado
2
1 α1 α1
V3 (α1 , α2 , α3 ) = 1 α2 α22 .
1 α3 α23
Podemos calcular este determinante directamente, pero es más ilustrativo usar las propie-
dades de los determinantes. Si restamos la primera fila al resto de las filas, el determinante
no varía:
α21
1 α1
V3 (α1 , α2 , α3 ) = 0 α2 − α1 α22 − α21 .
0 α3 − α1 α23 − α21
Si desarrollamos ahora por la primera columna (método de Laplace) obtenemos la fór-
mula
α2 − α1 α2 − α2 α2 − α1 (α2 − α1 )(α2 + α1 )
2 1
V3 (α1 , α2 , α3 ) = =
α3 − α1 α23 − α21
α3 − α1 (α3 − α1 )(α3 + α1 )
1 α2 + α1
= (α3 − α1 )(α2 − α1 )
1 α3 + α1
= (α3 − α2 )(α3 − α1 )(α2 − α1 ) .
Para pasar de la primera a la segunda línea hemos sacado factor común α2 − α1 en la
primera fila y α3 − α1 en la segunda. Es obvio que, si alguna de las αi coinciden, este
determinante se anula. También se ve que el resultado tiene una estructura parecida a la
del caso n = 2.
3
Caso n = 4
Este caso requiere bastante más álgebra. Sea
1 α1 α21 α31
1 2 3
α2 α2 α2
V4 (α1 , α2 , α3 , α4 ) = .
2 3
1 α3 α3 α3
2 3
1 α4 α4 α4
Al igual que en el caso anterior, restamos la primera fila al resto de las filas y desarrolla-
mos por la primera columna. Obtenemos que
α2 − α1 α22 − α21 α32 − α31
V4 (α1 , α2 , α3 , α4 ) = α3 − α1 α23 − α21 α33 − α31 .
α4 − α1 α24 − α21 α34 − α31
Como a3 − b3 = (a − b)(a2 + ab + b2 ), podemos sacar factor común αk − α1 de la fila
k-ésima y obtener:
1 α2 + α1 α22 + α2 α1 + α21
V4 (α1 , α2 , α3 , α4 ) = (α4 − α1 )(α3 − α1 )(α2 − α1 ) 1 α3 + α1 α23 + α3 α1 + α21 .
1 α4 + α1 α24 + α4 α1 + α21
Ahora volvemos a restar la primera fila al resto de las filas y obtenemos
α22 + α2 α1 + α21
1 α2 + α1
V4 (α1 , α2 , α3 , α4 ) = (α4 − α1 )(α3 − α1 )(α2 − α1 ) 0 α3 − α2 (α3 − α2 )(α1 + α2 + α3 ) .
0 α4 − α2 (α4 − α2 )(α1 + α2 + α4 )
Volvemos a desarrollar por la primera columna, sacamos los factores comunes de cada
fila y obtenemos
1 α1 + α2 + α3
V4 (α1 , α2 , α3 , α4 ) = (α4 − α1 )(α3 − α1 )(α2 − α1 )(α3 − α2 )(α4 − α2 ) .
1 α1 + α2 + α4
4
De aquí se deduce la fórmula final
V4 (α1 , α2 , α3 , α4 ) = (α4 − α3 )(α4 − α2 )(α4 − α1 )(α3 − α2 )(α3 − α1 )(α2 − α1 ) .
Es importante darse cuenta de que la estructura de este resultado es muy parecida a la
del caso n = 3.
Caso general
Los anteriores ejemplos sugieren que la fórmula general para un determinante de Van-
dermonde es
Y
Vn (α1 , . . . , αn ) = (αj − αi )
16i<j6n
= (αn − αn−1 ) · · · (αn − α1 ) (αn−1 − αn−2 ) · · · (αn−1 − α1 ) · · ·
×(α3 − α2 )(α3 − α1 ) (α2 − α1 ) .
Una prueba sencilla de este resultado se obtiene por inducción (que es un método general
que se estudiará en las asignaturas de Cálculo y Matemática Discreta).
