Matriz de Hilbert

En álgebra lineal, una Matriz de Hilbert es una matriz cuadrada cuyos campos constituyen una fracción de la unidad.

${\displaystyle H_{ij}={\frac {1}{i+j-1}},\ \ \ i,j=1...n}$

Por ejemplo, esta sería una matriz de Hilbert de 30 × 30:

${\displaystyle M={\begin{bmatrix}b_{1}&c_{1}&0&0&\cdots &\cdots &0\\[4pt]a_{2}&b_{2}&c_{2}&0&\cdots &\cdots &0\\[4pt]0&a_{3}&b_{3}&c_{3}&\cdots &\cdots &\vdots \\[4pt]0&0&a_{4}&b_{4}&c_{4}&\cdots &\vdots \\[4pt]\vdots &\ddots &\ddots &\ddots &\ddots &\ddots &\vdots \\[4pt]\vdots &\vdots &\vdots &\vdots &a_{n-1}&b_{n-1}&c_{n-1}\\[4pt]0&0&\cdots &0&\cdots &a_{n}&b_{n}\end{bmatrix}}.}$

La matriz de Hilbert se puede obtener a partir de la integral

${\displaystyle H_{ij}=\int _{0}^{1}x^{i+j-2}\,dx,}$

es decir, una Matriz de Gram para potencias de x. Representa la aproximación de mínimos cuadrados a funciones arbitrarias por polinomios.

Las matrices de Hilbert son ejemplos de matrices mal condicionadas, lo que las hace notoriamente difícil su uso en cálculo numérico. Por ejemplo, el número de condición de la matriz mostrada anteriormente es aproximadamente 4.8 · 10⁵.

Referencias[editar]