En matemáticas, la raíz cuadrada de una matriz extiende la noción de raíz cuadrada de los números a las matrices. Una matriz B se dice que es una raíz cuadrada de A si el producto matricial BB es igual a A.[1]
Introducción
La existencia de un producto de matrices permite definir la raíz de una matrices como aquella matriz que multiplicada por sí misma da la original.
Si
es una matriz definida positiva u operador, entonces existe exactamente una matriz definida positiva u operador
tal que
; entonces definimos
. Dada una matriz real su raíz cuadrada está definida si su espectro puntual está formado por números positivos, si el espectro no fuera estrictamente positivo la raíz cuadrada de una matriz involucrará matrices con coeficientes complejos.
Más generalmente, para cada matriz u operador normal
existen operadores normales
tales que
. En general, hay muchos de esos operadores
para cada
y entonces la función raíz cuadrada no puede ser definida satisfactoriamente para operadores normales. En cierta manera se puede decir que los operadores definidos positivos son similares a los números reales positivos, y los operadores normales son similares a los números complejos.
Algoritmos de cálculo
Método simplificado de Newton
Si A es una matriz n × n con valores complejos, el siguiente algoritmo — Método simplificado de Newton — aproxima la matriz
a la raíz cuadrada de A tras k iteraciones:[1]
Sea X0 = I, donde I es la matriz identidad. La iteración está definida por
![{\displaystyle X_{k+1}={\tfrac {1}{2}}(X_{k}+AX_{k}^{-1})\,.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/50a79cb4877bb73fea093c0b84432aa81c68629d)
La convergencia no está asegurada, pero si el proceso converge, la matriz
converge cuadráticamente a la raíz cuadrada A1/2. Este método es una extensión del algoritmo babilónico para el cálculo de raíces cuadradas de números positivos ordinarios.
Método IASMN
Para una matriz A real definida positiva existe un algoritmo denominado Iteración alternativa simplificada del Método de Newton, muy sencillo y computacionalmente eficiente para calcular la raíz cuadrada de este importante tipo de matrices. El algoritmo es el siguiente:[2]
El símbolo "\" representa el proceso de eliminación Gaussiana.
Véase también
Referencias
Enlaces externos