Álgebra Superior I: Producto de matrices con matrices

Por Eduardo García Caballero

Introducción

Hasta ahora hemos conocido varias operaciones que involucran escalares, vectores y matrices. En esta entrada aprenderemos sobre una de las operaciones más importantes en el álgebra lineal: el producto de matrices con matrices.

Definición de producto de matrices

Para poder efectuar el producto de dos matrices, hay que asegurarnos de que el número de columnas de la primera matriz sea igual al número de filas de la segunda matriz.

El resultado de una matriz A de tamaño m×n por una matriz B de tamaño n× será la matriz C=AB de tamaño m×, donde la entrada cij de C está dada por la fórmula
cij=ai1b1j+ai2b2j++ainbnj.

A primera vista esta fórmula puede parecer complicada, sin embargo, practicando con algunos ejemplos verás que es muy fácil de implementar.

  • Producto de matrices de tamaño 2×2:

Sean
A=(1357)yB=(2468).

Como estamos multiplicando una matriz de tamaño 2×2 por una matriz de tamaño 2×2, sabemos que el resultado será otra matriz de tamaño 2×2. Ahora, iremos calculando una por una sus entradas.

Sea C=AB. Para calcular la entrada c11 observamos la primera fila de A y la primera columna de B, las cuales son
A=(1357)yB=(2468),
de modo que c11=(1)(2)+(3)(6)=20:
AB=(20285276).

Para la entrada c12, nos fijamos en la primera columna de A y en la segunda columna de B, que son
A=(1357)yB=(2468),
obteniendo c12=(1)(4)+(3)(8)=28:
AB=(20285276).

De manera similar, observemos la segunda fila de A y la primera columna de B,
A=(1357),B=(2468),
obteniendo c21=(5)(2)+(7)(6)=52, mientras que la segunda fila de A y la segunda columna de B son
A=(1357),B=(2468),
obteniendo c22=(5)(4)+(7)(8)=76.

Por lo tanto,
AB=(20285276).

En general, el resultado del producto de las matrices
A=(a11a12a21a22)yB=(b11b12b21b22)
es
AB=(a11a12a21a22)(b11b12b21b22)=(a11b11+a12b21a11b12+a12b22a21b11+a22b21a21b12+a22b22).

  • Producto de matriz de 3×2 por matriz de 2×2:

Supongamos que
A=(351043)yB=(7852).

En este caso, como estamos multiplicando una matriz de tamaño 3×2 por una matriz de tamaño 2×2, la matriz resultante tendrá tamaño 3×2.

Podemos obtener sus entradas de manera similar al caso anterior. Si C=AB, entonces la entrada c12 la podemos encontrar revisando la primera fila de A y la segunda columna de B,
A=(351043),B=(7852).
de modo que c12=(3)(8)+(5)(2)=34. Por su parte, para obtener la entrada c31 nos fijamos en la tercera fila de A y la primera columna de B,
A=(351043),B=(7852).
obteniendo c31=(4)(7)+(3)(5)=43.

¿Podrías comprobar que
AB=(4634784338)?

Así, para el caso general de matrices de 3×2 por 2×2, obtendremos
(a11a12a21a22a31a32)(b11b12b21b22)=(a11b11+a12b21a11b12+a12b22a21b11+a22b21a21b12+a22b22a31b11+a32b21a31b12+a32b22).

  • Producto de matriz de 4×2 por matriz de 2×3:

¿Podrías verificar que la siguiente fórmula es correcta?
(a11a12a21a22a31a32a41a42)(b11b12b13b21b22b23)=(a11b11+a12b21a11b12+a12b22a11b13+a12b23a21b11+a22b21a21b12+a22b22a21b13+a22b23a31b11+a32b21a31b12+a32b22a31b13+a32b23a41b11+a42b21a41b12+a42b22a41b13+a42b23).

Propiedades del producto de matrices

A continuación revisaremos algunas de las propiedades que cumple la multiplicación de matrices. Para demostrar las siguientes propiedades, consideraremos la matriz A de tamaño 3×2 y las matrices B y C de tamaño 2×2, aunque se pueden probar para matrices de cualesquier otro tamaño entre las cuales se puedan efectuar las operaciones.

Veamos que si efectuamos la multiplicación de una matriz de tamaño m×n por una matriz de tamaño n×1 siguiendo el algoritmo descrito anteriormente, el resultado coincide con el de multiplicar la matriz de tamaño m×n por un vector de tamaño n. Por ejemplo, si multiplicamos A por una matriz U de tamaño 2×1, obtendremos
(a11a12a21a22a31a32)(u11u12)=(a11u11+a12u21a21u11+a22u21a31u11+a32u21).

Esta es una observación importante pues todo lo que demostremos para el producto de matrices también lo tendremos para el producto de matriz por vector.

Veamos que la multiplicación de matrices es asociativa:

(AB)C=((a11a12a21a22a31a32)(b11b12b21b22))(c11c12c21c22)=(a11b11+a12b21a11b12+a12b22a21b11+a22b21a21b12+a22b22a31b11+a32b21a31b12+a32b22)(c11c12c21c22)=((a11b11+a12b21)c11+(a11b12+a12b22)c21(a11b11+a12b21)c12+(a11b12+a12b22)c22(a21b11+a22b21)c11+(a21b12+a22b22)c21(a21b11+a22b21)c12+(a21b12+a22b22)c22(a31b11+a32b21)c11+(a31b12+a32b22)c21(a31b11+a32b21)c12+(a31b12+a32b22)c22)=(a11(b11c11+b12c21)+a12(b21c11+b22c21)a11(b11c12+b12c22)+a12(b21c12+b22c22)a21(b11c11+b12c21)+a22(b21c11+b22c21)a21(b11c12+b12c22)+a22(b21c12+b22c22)a31(b11c11+b12c21)+a32(b21c11+b22c21)a31(b11c12+b12c22)+a32(b21c12+b22c22))=(a11a12a21a22a31a32)(b11c11+b12c21b11c12+b12c22b21c11+b22c21b21c12+b22c22)=(a11a12a21a22a31a32)((b11b12b21b22)(c11c12c21c22))=A(BC).

De manera muy similar, si u es un vector de tamaño 2, podemos ver que se cumple que A(Bu)=(AB)u. ¿Puedes demostrarlo? Hazlo por lo menos para matrices A y B ambas de 2×2.

Quizás tengas la impresión de que hay que hacer demasiadas cuentas y que sería sumamente difícil demostrar estas propiedades para matrices más grandes. Sin embargo, en cursos posteriores verás cómo trabajar apropiadamente con la notación para poder hacer estas demostraciones más fácilmente.

El producto de matrices es asociativo. Sin embargo, no es conmutativo. Por ejemplo, consideremos las matrices
E=(5730)yF=(1291).


Veamos que
EF=(68336)(174863)=FE.

En términos de combinar el producto de matrices con otras operaciones, tenemos que el producto de matrices por la izquierda se distribuye sobre la suma de matrices:
A(B+C)=(a11a12a21a22a31a32)((b11b12b21b22)+(c11c12c21c22))=(a11a12a21a22a31a32)(b11+c11b12+c12b21+c21b22+c22)=(a11(b11+c11)+a12(b21+c21)a11(b12+c21)+a12(b22+c22)a21(b11+c11)+a22(b21+c21)a21(b12+c21)+a22(b22+c22)a31(b11+c11)+a32(b21+c21)a31(b12+c21)+a32(b22+c22))=(a11b11+a11c11+a12b21+a12c21a11b12+a11c11+a12b22+a12c22a21b11+a21c11+a22b21+a22c21a21b12+a21c12+a22b22+a22c22a31b11+a31c11+a32b21+a32c21a31b12+a31c12+a32b22+a32c22)=(a11b11+a12b21a11b12+a12b22a21b11+a22b21a21b12+a22b22a31b11+a32b21a31b12+a32b22)+(a11c11+a12c21a11c12+a12c22a21c11+a22c21a21c12+a22c22a31c11+a32c21a31c12+a32c22)=(a11a12a21a22a31a32)(b11b12b21b22)+(a11a12a21a22a31a32)(c11c12c21c22)=AB+AC.

El producto también se distribuye sobre la suma cuando la suma aparece a la izquierda. ¿Podrías probar que si D es una matriz de tamaño 3×2, entonces se cumple (A+D)B=AB+DB?

En entradas anteriores vimos que In tiene la propiedad de ser neutro al multiplicarla por un vector de tamaño n. Resulta que In también tiene esta propiedad al multiplicarla por la izquierda por una matriz de tamaño n×m. Por ejemplo, veamos que al multiplicar I3 por la izquierda por A, obtenemos
I3A=(100010001)(a11a12a21a22a31a32)=(1a11+0a21+0a311a12+0a22+0a320a11+1a21+0a310a12+1a22+0a320a11+0a21+1a310a12+0a22+1a32)=(a11a12a21a22a31a32)=A.

¿Podrías probar que AI2=A (es decir, que I2 es neutro por la derecha para A)?

Habiendo visto que el producto de matrices es asociativo, conmutativo y tiene neutros, probablemente te estarás preguntando si existen inversos en la multiplicación de matrices. Este cuestionamiento lo dejaremos para la siguiente entrada.

Relación con la composición de transformaciones

Como vimos en la entrada anterior, una forma de visualzar el producto de una matriz A por un vector u es como una transformación que envía el vector u a un único vector Au.

Teniendo en mente esto, veamos que la propiedad de que A(Bu)=(AB)u resulta aún más interesante. Para esto, veamos que el siguiente ejemplo: sean
A=(0211),B=(1230),yu=(12).

Si multiplicamos B por u, vemos que corresponde a la transformación que envía u=(12) al vector Bu=(53).

Ahora, si multiplicamos A por el vector Bu, vemos que corresponde a la transformación que envía Bu al vector A(Bu)=(68) (Acabamos de obtener el resultado de aplicar a u la composición de las transformaciones B y A).

Por otra parte, si realizamos la multiplicación
AB=(0211)(1230)=(6042),
la transformación asociada a AB envía u al vector (AB)u=(68).

¡La composición de las transformaciones asociadas a B y A aplicada al vector u coincide con la transformación asociada a la matriz AB aplicada al mismo vector!

Si probamos esto para un vector arbitrario, nos daremos cuenta de que en todos los casos se cumple lo mismo. En realidad, esto no es una coincidencia: como aprenderás en tus cursos de álgebra lineal, la composición de transformaciones lineales está directamente asociada al producto de matrices.

Potencias de matrices

Podemos ver que si una matriz A es cuadrada, al tener el mismo número de filas que de columnas, entonces podemos realizar la multiplicaciones AA, AAA, AAAA, etc., que por asociatividad no importa en qué orden multipliquemos. Esto nos sugiere que podemos cacular potencias de matrices.

Para una matriz cuadrada A, definiremos de manera recursiva la potencia An:

  • Definimos A0=I.
  • Dada An, con n un número natural, definimos An+1=AnA.

Por ejemplo, si
A=(2134),
calculemos A3 empleando la definición recursiva. Para esto, iremos calculando una por una las potencias de A, hasta llegar a A3:
A0=I=(1001),A1=A0A=(1001)(2134)=(2134),A2=A1A=(2134)(2134)=((2)(2)+(1)(3)(2)(1)+(1)(4)(3)(2)+(4)(3)(3)(1)+(4)(4))=(761819),A3=A2A=(761819)(2134)=((7)(2)+(6)(3)(7)(1)+(6)(4)(18)(2)+(19)(3)(18)(1)+(19)(4))=(32319394).

Prueba calcular algunas potencias de la matriz (2003). ¿Notas algún patrón especial?

Más adelante…

En esta entrada aprendimos sobre el producto de matrices con matrices y conocimos algunas de sus propiedades. En la siguiente entrada abordaremos la pregunta sobre si existen los inversos en la multiplicación de matrices.

Tarea moral

  1. Realiza el producto de matrices (123012113)(111111111).
  2. Considera la matriz A=(3445). Realiza las siguientes operaciones por separado, sin usar la asociatividad del producto de matrices. ¿Cuál de las dos operaciones te resultó más fácil de hacer?
    • A(A(A(A(23)))).
    • (((AA)A)A)(23).
  3. Completa las pruebas faltantes de las propiedades de la multiplicación de matrices.
  4. Demuestra la siguiente ley de exponentes para matrices: AmAn=Am+n.
  5. Prueba que si
    A=(a1100a22),
    y k es un entero mayor o igual que 0, entonces
    Ak=(a11k00a22k)
    (Sugerencia: realizarlo por inducción sobre k, utilizando la definición recursiva).
  6. Encuentra matrices A y B de 2×2 para las cuales A2B2(A+B)(AB).

Entradas relacionadas

Deja una respuesta

Tu dirección de correo electrónico no será publicada. Los campos obligatorios están marcados con *

Este sitio usa Akismet para reducir el spam. Aprende cómo se procesan los datos de tus comentarios.