Archivo de la etiqueta: fibonacci

Álgebra Lineal I: Técnicas básicas de cálculo de determinantes

Por Leonardo Ignacio Martínez Sandoval

Introducción

Ya definimos a los determinantes para vectores, para transformaciones y para matrices. Además, mostramos algunas propiedades básicas de determinantes y las usamos para resolver varios problemas. Como hemos discutido, los determinantes guardan información importante sobre una transformación lineal o sobre una matriz. También ayudan a implementar la técnica de diagonalización la cual introdujimos hace algunas entradas y en la cual profundizaremos después. Es por esta razón que es importante tener varias técnicas para el cálculo de determinantes.

Fuera de este curso, los determinantes sirven en muchas otras áreas de las matemáticas. Cuando se hace cálculo de varias variables ayudan a enunciar el teorema del cambio de variable. En combinatoria ayudan a calcular el número de árboles generadores de una gráfica. Más adelante en tu formación matemática es probable que te encuentres con otros ejemplos.

Calculo de determinantes de 2×2

Como ya discutimos anteriormente, una matriz en M2(F), digamos A=(abcd) tiene determinante adbc.

Problema. Calcula el determinante de la matriz (0111)8.

Solución. Por la fórmula para el determinante de las matrices de 2×2, se tiene que |0111|=0111=1.

Como el determinante es multiplicativo, det(A2)=det(A)det(A)=(det(A))2, e inductivamente se puede mostrar que para todo entero positivo n se tiene que det(An)=(det(A))n. De esta forma, el determinante que buscamos es (1)8=1.

Observa que hubiera tomado más trabajo elevar la matriz a la octava potencia. Aunque esto usualmente no es recomendable, en este problema hay algo interesante que sucede con esta matriz. Llamémosla A=(0111). Haciendo las cuentas para las primeras potencias, se tiene que
A=(0111)A2=(1112)A3=(1223)A4=(2335)A5=(3558)

Aquí aparece la sucesión de Fibonacci, dada por F0=0, F1=1 y Fn+2=Fn+1+Fn para n0, cuyos primeros términos son 0,1,1,2,3,5,8,13,21,. De hecho se puede probar por inducción que An=(Fn1FnFnFn+1).

Así, por un lado el determinante de la matriz An es Fn1Fn+1Fn2, usando la fórmula de determinante de 2×2. Por otro lado, es (1)n, por el argumento del problema. Con esto hemos demostrado que para cualquier entero n tenemos la siguiente identidad para los números de Fibonacci: Fn1Fn+1Fn2=(1)n.

Cálculo de determinantes de 3×3

Para calcular el determinante de una matriz en M3(F) por definición, digamos de A=(abcdefghi), tenemos que hacer una suma de 3!=6 términos. Si se hacen las cuentas de manera explícita, el valor que se obtiene es aei+bfg+cdhcegafhbdi.

Esto se puede recordar mediante el siguiente diagrama, en el cual se ponen la primera y la segunda columna de nuevo, a la derecha. Las diagonales hacia abajo son términos positivos y las diagonales hacia arriba son términos negativos.

Cálculo de determinantes de matrices de 3x3
Cálculo de determinantes de 3×3

Veamos un ejemplo de un problema en el que se puede aprovechar esta técnica.

Problema. Determina para qué reales a,b,c se tiene que los vectores (a,b,0), (a,0,b) y (0,a,b) son una base de R3.

Solución. Para que estos vectores sean una base de R3, basta con que sean linealmente independientes, pues son 3. Como hemos visto en entradas anteriores, para que sean linealmente independientes, es necesario y suficiente que el determinante de la matriz (ab0a0b0ab) sea distinto de cero.

Usando la técnica de arriba, hacemos siguiente diagrama:

De aquí, vemos que el determinante es 0+0+00a2bab2=ab(a+b). Esta expresión es igual a cero si a=0, si b=0 o si a+b=0. En cualquier otro caso, el determinante no es cero, y por lo tanto los vectores forman una base.

Ten mucho cuidado. Esta técnica no funciona para matrices de 4×4 o más. Hay una forma sencilla de convencerse de ello. Por ejemplo, el determinante de una matriz de 4×4 debe tener 4!=24 sumandos. Si intentamos copiar la técnica de arriba, tendremos solamente 8 sumandos (4 en una diagonal y 4 en otra). Para cuando tenemos matrices de 4×4 o más, tenemos que recurrir a otras técnicas.

Reducción gaussiana para determinantes

Cuando vimos el tema de sistemas de ecuaciones hablamos del algoritmo de reducción gaussiana, y vimos que este siempre lleva una matriz en Mm,n(F) a su forma escalonada reducida mediante operaciones elementales. Cuando aplicamos el algoritmo a matrices en Mn(F), siempre llegamos a una matriz triangular, en donde sabemos fácilmente calcular el determinante: es simplemente el producto de las entradas en la diagonal. Nota cómo lo anterior también se cumple para las matrices diagonales, pues son un caso particular de matrices triangulares.

Por esta razón, es fundamental para el cálculo de determinantes saber qué le hacen las operaciones elementales al determinante de una matriz.

Teorema. Las operaciones elementales tienen el siguiente efecto en el determinante de una matriz A:

  1. Si todos los elementos de un renglón o columna de A se multiplican por λ, entonces el determinante se multiplica por λ.
  2. Cuando se intercambian dos renglones o columnas de A, el determinante se multiplica por 1.
  3. Si a un renglón de A se le suma un múltiplo escalar de otro renglón, entonces el determinante no cambia. Sucede algo análogo para columnas.

Demostración. El punto 1 ya lo demostramos en la entrada anterior, en donde vimos que el determinante es homogéneo.

Para los puntos 2 y 3, usemos que si e1,en es la base canónica de Fn, el determinante de una matriz con renglones R1,,Rn es det(e1,,en)(R1,,Rn).

Intercambiar los renglones i y j es hacer det(e1,,en)(Rσ(1),,Rσ(n)) para la transposición σ que intercambia i y j. Como el determinante es antisimétrico y σ tiene signo 1, obtenemos la conclusión.

Hagamos ahora el tercer punto. Tomemos ij y un escalar λ. Si al i-ésimo renglón de A le sumamos λ veces el j-ésimo renglón de A, esto es lo mismo que multiplicar a A por la izquierda por la matriz B que tiene unos en la diagonal y λ en la entrada (i,j). La matriz B es triangular, de modo que su determinante es el producto de las entradas, que es 1. De esta forma, det(BA)=det(B)det(A)=det(A).

◻

Así, una estrategia para calcular el determinante de una matriz es hacer reducción gaussiana hasta llegar a una matriz diagonal (incluso es suficiente que sea triangular superior) de determinante Δ. Si en el camino se hicieron r intercambios de renglones y se multiplicaron los renglones por escalares λ1,,λs, entonces el determinante de A será (1)rΔλ1λs.

Otras propiedades para calcular determinantes

Aquí recolectamos otras propiedades de determinantes que pueden ayudar a calcularlos. Ya mostramos todas ellas, salvo la número 2. Esta la mostramos después de la lista.

  1. Si se descompone una columna de una matriz como suma de dos columnas, entonces el determinantes es la suma de los determinantes en los que ponemos cada columna en vez de la original.
  2. Si A es una matriz en Mn(C), entonces el determinante de la matriz conjugada A es el conjugado del determinante de A.
  3. El determinante es multiplicativo.
  4. Si A es una matriz en Mn(F), el determinante de λA es λn veces el determinante de A.
  5. El determinante de una matriz triangular es el producto de sus entradas en la diagonal.
  6. El determinante de una matriz invertible es el inverso multiplicativo del determinante de la matriz.
  7. Una matriz tiene el mismo determinante que su transpuesta.

Proposición. Si A es una matriz en Mn(C), entonces el determinante de la matriz conjugada A es el conjugado del determinante de A.

Demostración. La conjugación compleja abre sumas y productos. Aplicando esto repetidas veces obtenemos la siguiente cadena de igualdades:

det(A)=σSnsign(σ)a1σ(1)anσ(n)=σSnsign(σ)a1σ(1)anσ(n)=σSnsign(σ)a1σ(1)anσ(n)=det(A).

◻

Hay una última técnica que es fundamental para el cálculo de determinantes: la expansión de Laplace. En algunos textos incluso se usa para definir el determinante. Probablemente la conoces: es la que consiste en hacer el determinante «con respecto a una fila o columna» y proceder de manera recursiva. Hablaremos de ella más adelante y veremos por qué funciona.

Dos problemas de cálculo de determinantes

Problema 1. Considera la matriz A=(5120011200521301). Calcula los siguientes determinantes:

  • detA
  • det(tA)
  • det(A1)
  • det(tAA)
  • det(2A)

Solución. Hagamos primero el determinante de la matriz A. Para ello, haremos operaciones elementales como sigue
(5120011200521301)(5120011200520145251)(51200112005200125335)(51200112005200018925).

En el primer paso sumamos 1/5 veces el primer renglón al último. Luego, sumamos 14/5 veces el segundo renglón al último. Finalmente, sumamos 12/25 veces el tercer renglón al último. De esta forma, nunca cambiamos el determinante de la matriz. Así, del determinante de A es el mismo que el de la matriz final, que por ser triangular superior es el producto de las entradas en su diagonal. De este modo, det(A)=5151895=189.

El determinante de una matriz es igual al de su transpuesta, así que det(tA)=det(A). El determinante det(A1) es el inverso multiplicativo de det(A), así que es 1189.

Como el determinante es multiplicativo, det(tAA)=det(tA)det(A)=189189=35721.

Finalmente, usando que el determinante es homogéneo y que estamos en M4(R), tenemos que
det(2A)=(2)4det(A)=16189=3024.

Problema 2. Sean a,b,c números complejos. Calculando el determinante de la matriz A=(abccabbca) en M3(C) de dos formas distintas, muestra que a3+b3+c33abc=(a+b+c)(a2+b2+c2abbcca).

Solución. Usando la técnica para determinantes de 33 tenemos que por un lado,
det(A)=a3+b3+c3abcbcacab=a3+b3+c33abc.

Por otro lado, el determinante no cambia si al primer renglón le sumamos los otros dos, así que el determinante de A también es |a+b+ca+b+ca+b+ccabbca|. Como el determinante es homogéneo, podemos factorizar a+b+c de la primera entrada para obtener que det(A)=(a+b+c)|111cabbca|.

Aplicando de nuevo la fórmula de determinantes de 3×3, tenemos que |111cabbca|=a2+b2+c2abbcca.

Concluimos entonces que det(A)=(a+b+c)(a2+b2+c2abbcca). Igualando ambas expresiones para det(A) obtenemos la identidad deseada.

Más adelante…

En esta entrada vimos varias formas para calcular el determinante de una matriz. Cuando nos enfrentemos con un problema que requiere el cálculo de un determinante, tenemos que elegir la que más nos convenga (o la que requiera menos pasos). La mejor forma de desarrollar un poco de «intuición» al momento de elegir el mejor método para calcular determinantes es haciendo ejercicios.

A continuación pondremos en práctica lo que aprendimos en esta entrada haciendo varios ejercicios de cálculo de determinantes.

Tarea moral

A continuación hay algunos ejercicios para que practiques los conceptos vistos en esta entrada. Te será de mucha utilidad intentarlos para entender más la teoría vista.

  • Sea α un número real. Encuentra el determinante de la matriz (sinαcosαcosαsinα).
  • Determina para qué valores de a la matriz (a0a0a0a0a000a000a0a0a0a0a) es invertible.
  • Encuentra el determinante de la matriz (2100002100002100002110002).
  • Sea x un número complejo. Muestra que el determinante de la matriz (3x26x+52x24x+2x22x2x24x+22x2+2x+1x2xx22xx2xx2) es x6. Sugerencia. Hay una solución simple, factorizando a la matriz como el producto de dos matrices triangulares, una superior y una inferior, una transpuesta de la otra.
  • Muestra que si A=(0111), entonces An=(Fn1FnFnFn+1), donde {Fn} es la sucesión de Fibonacci. Muestra que para los números de Fibonacci se satisface que F2n=Fn(Fn+1+Fn1).

Entradas relacionadas

Agradecimientos

Trabajo realizado con el apoyo del Programa UNAM-DGAPA-PAPIME PE104721 «Hacia una modalidad a distancia de la Licenciatura en Matemáticas de la FC-UNAM»

Elegir notación efectiva

Por Leonardo Ignacio Martínez Sandoval

HeuristicasPara que las matemáticas realmente simplifiquen las cosas y no las compliquen más, una de las cosas que se necesita es tener una notación adecuada.

En esta serie de videos veremos algunos ejemplos en los cuales elegir variables adecuadas o una representación adecuada del problema puede ayudar en la solución del problema o bien simplifica algunas cuentas.

Ir a los videos…