Archivo de la etiqueta: determinante

Á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»

Álgebra Lineal I: Determinantes de matrices y transformaciones lineales

Por Leonardo Ignacio Martínez Sandoval

Introducción

En la entrada anterior dimos la definición de determinante para ciertos vectores con respecto a una base. En esta entrada continuamos con la construcción de determinantes. Primero, basados en la teoría que desarrollamos anteriormente, definiremos determinantes de transformaciones lineales. Luego, mediante la cercanía entre transformaciones lineales y matrices, definimos determinantes de matrices.

Determinantes de transformaciones lineales

Ahora definiremos el determinante para transformaciones lineales. Antes de esto, necesitamos hacer algunas observaciones iniciales y demostrar un resultado.

Si tomamos un espacio vectorial V de dimensión finita n1 sobre un campo F, una transformación lineal T:VV y una forma n-lineal f:VnF, se puede mostrar que la transformación Tf:VnF dada por Tf(x1,,xn)=f(T(x1),,T(xn)) también es una forma n-lineal. Además, se puede mostrar que si f es alternante, entonces Tf también lo es. Mostrar ambas cosas es relativamente sencillo y queda como tarea moral.

Teorema. Sea V un espacio vectorial de dimensión finita n1 sobre el campo F. Para cualquier transformación lineal T:VV existe un único escalar detT en F tal que f(T(x1),,T(xn))=detTf(x1,,xn) para cualquier forma n-lineal alternante f:VnF y cualquier elección x1,,xn de vectores en V.

Demostración. Fijemos una base B=(b1,,bn) cualquiera de V. Llamemos g a la forma n-lineal alternante det(b1,,bn). Por la discusión de arriba, la asignación Tg:VnF dada por (x1,,xn)g(T(x1),,T(xn)) es una forma n-lineal y alternante.

Por el teorema que mostramos en la entrada de determinantes de vectores, se debe cumplir que Tg=Tg(b1,,bn)g. Afirmamos que detT:=Tg(b1,,bn) es el escalar que estamos buscando.

En efecto, para cualquier otra forma n-lineal alternante f, tenemos por el mismo teorema que f=f(b1,,bn)g. Usando la linealidad de T y la igualdad anterior, se tiene que

Tf=f(b1,,bn)Tg=f(b1,,bn)detTg=detTf.

Con esto se prueba que detT funciona para cualquier forma lineal f. La unicidad sale eligiendo (x1,,xn)=(b1,,bn) y f=g en el enunciado del teorema, pues esto forza a que detT=g(T(b1),,T(bn)).

◻

Ahora sí, estamos listos para definir el determinante de una transformación lineal.

Definición. El escalar detT del teorema anterior es el determinante de la transformación lineal T.

Para obtener el valor de detT, podemos entonces simplemente fijar una base B=(b1,,bn) y el determinante estará dado por detT=det(b1,,bn)(T(b1),,T(bn)). Como el teorema también prueba unicidad, sin importar que base B elijamos este número siempre será el mismo.

Ejemplo 1. Vamos a encontrar el determinante de la transformación lineal T:R3R3 dada por T(x,y,z)=(2z,2y,2x). Para ello, usaremos la base canónica de R3. Tenemos que
T(1,0,0)=(0,0,2)=2e3T(0,1,0)=(0,2,0)=2e2T(0,0,1)=(2,0,0)=2e1.

De acuerdo al teorema anterior, podemos encontrar al determinante de T como detT=det(e1,e2,e3)(2e3,2e2,2e1).

Como el determinante (para vectores) es antisimétrico, al intercambiar las entradas 1 y 3 su signo cambia en 1. Usando la 3-linealidad en cada entrada, podemos sacar un factor 2 de cada una. Así, tenemos:
detT=det(e1,e2,e3)(2e3,2e2,2e1)=det(e1,e2,e3)(2e1,2e2,2e3)=8det(e1,e2,e3)(e1,e2,e3)=8.

Concluimos entonces que el determinante de T es 8.

Ejemplo 2. Vamos ahora a encontrar el determinante de la transformación T:Rn[x]Rn[x] que deriva polinomios, es decir, tal que T(p)=p. Tomemos q0=1,q1=x,,qn=xn la base canónica de Rn[x].

Notemos que, T(1)=0, de modo que los vectores T(1),,T(xn) son linealmente dependientes. Así, sin tener que hacer el resto de los cálculos, podemos deducir ya que det(q0,,qn)(T(q0),,T(qn))=0. Concluimos entonces que detT=0.

Determinantes de matrices

La expresión detT=det(b1,,bn)(T(b1),,T(bn)) para una transformación lineal T también nos permite poner al determinante en términos de las entradas de la matriz de T con respecto a la base B. Recordemos que dicha matriz AT=[aij] tiene en la columna i las coordenadas de bi en la base B. En otras palabras, para cada i se cumple que T(bi)=j=1najibi.

Usando esta notación, obtenemos que detT=σSnsign(σ)a1σ(1)anσ(n), de manera que podemos expresar a detT en términos únicamente de su matriz en la base B.

Esto nos motiva a definir el determinante de una matriz en general.

Definición. Para una matriz A en Mn(F) de entradas A=[aij], el determinante de A es detA=σSnsign(σ)a1σ(1)anσ(n). A detA también lo escribimos a veces en notación de «matriz con barras verticales» como sigue:

detA=|a11a12a1na21a22a2nan1an2ann.|

Ejemplo. Si queremos calcular el determinante de una matriz en M2(F), digamos A=(abcd), debemos considerar dos permutaciones: la identidad y la transposición (1,2).

La identidad tiene signo 1 y le corresponde el sumando ad. La transposición tiene signo 1 y le corresponde el sumando bc. Así, |abcd|=adbc.

Retomando la discusión antes de la definición, tenemos entonces que detT=detAT, en donde a la izquierda hablamos de un determinante de transformaciones lineales y a la derecha de uno de matrices. La matriz de T depende de la base elegida, pero como vimos, el determinante de T no. Esta es una conclusión muy importante, y la enunciamos como teorema en términos de matrices.

Teorema. Sean A y P matrices en Mn(F) con P invertible. El determinante de A y el de P1AP son iguales.

Determinantes de matrices triangulares

Terminamos esta entrada con un problema que nos ayudará a repasar la definición y que más adelante servirá para calcular determinantes.

Problema. Muestra que el determinante de una matriz triangular superior o triangular inferior es igual al producto de las entradas de su diagonal.

Solución. En una matriz triangular superior tenemos que aij=0 si i>j. Vamos a estudiar la expresión σSnsign(σ)a1σ(1)anσ(n).

Si una permutación σ no es la identidad, entonces hay un entero i que no deja fijo, digamos σ(i)i. Tomemos a i como el mayor entero que σ no deja fijo. Notemos que σ(i) tampoco queda fijo por σ pues σ(σ(i))=σ(i) implica σ(i)=i, ya que σ es biyectiva, y estamos suponiendo σ(i)i. Por la maximalidad de i, concluimos que σ(i)<i.Entonces el sumando correspondiente a σ es 0 pues tiene como factor a la entrada aiσ(i)=0.

En otras palabras, la única permutación a la que le puede corresponder un sumando no cero es la identidad, cuyo signo es 1. De esta forma,
det(A)=σSnsign(σ)a1σ(1)anσ(n)=a11ann.

◻

Más adelante…

En esta entrada planteamos cómo se define el concepto de matriz para transformaciones lineales y cómo esta definición se extiende naturalmente a la definición del determinante de una matriz, recordando que a cada transformación lineal se le puede asociar una matriz y viceversa.

En las siguientes entradas vamos a ver qué propiedades que cumplen los determinantes y aprenderemos diferentes técnicas para calcularlos. A lo largo de la unidad, desarrollaremos bastante práctica en el cálculo y la manipulación de los determinantes, ya sea el determinante de un conjunto de vectores, de una transformación lineal o de una matriz.

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.

  • Muestra que la transformación Tf definida en la entrada es n-lineal y alternante.
  • Usando la definición de determinante para transformaciones lineales, encuentra el determinante de la transformación lineal T:RnRn dada por T(x1,x2,,xn)=(x2,x3,,x1).
  • Calcula por definición el determinante de las matrices (3241) y (1241391416).
  • Calcula por definición el determinante de la matriz (1112344916) y compáralo con el de la matriz de 3×3 del inciso anterior. ¿Qué notas?
  • Completa el argumento para mostrar que el determinante de una matriz triangular inferior es el producto de las entradas en su diagonal.

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»

Álgebra Lineal I: Determinantes de vectores e independencia lineal

Por Leonardo Ignacio Martínez Sandoval

Introducción

En este cuarto y último bloque del curso comenzamos hablando de transformaciones multilineales y de permutaciones. Luego, nos enfocamos en las transformaciones multilineales antisimétricas y alternantes. Con la teoría que hemos desarrollado hasta ahora, estamos listos para definir determinantes de vectores, de transformaciones lineales y de matrices.

En esta entrada comenzaremos con la definición de determinantes de vectores. En la siguiente entrada hablaremos acerca de determinantes de matrices y de transformaciones lineales. Después de definir determinantes, probaremos varias de las propiedades que satisfacen. Posteriormente, hablaremos de varias técnicas que nos permitirán calcular una amplia variedad de determinantes para tipos especiales de matrices.

Determinantes de vectores

Para empezar, definiremos qué es el determinante de un conjunto de vectores en un espacio de dimensión finita con respecto a una base.

Definición. Sea B=(b1,,bn) una base de un espacio vectorial V de dimensión finita n y x1,,xn vectores de V. Cada uno de los xi se puede escribir como xi=j=1najibj.

El determinante de x1,,xn con respecto a (b1,,bn) es σSnsign(σ)a1σ(1)anσ(n), y lo denotamos por det(b1,,bn)(x1,,xn).

Observa que estamos sumando tantos términos como elementos en Sn. Como existen n! permutaciones de un conjunto de n elementos, entonces la suma de la derecha tiene n! sumandos.

Ejemplo. Consideremos la base b1=1, b2=1+x y b3=1+x+x2 del espacio vectorial R2[x] de polinomios con coeficientes reales y grado a lo más 2. Tomemos los polinomios v1=1, v2=2x y v3=3x2. Vamos a calcular el determinante de v1,v2,v3 con respecto a la base (b1,b2,b3).

Para hacer eso, lo primero que tenemos que hacer es expresar a v1,v2,v3 en términos de la base. Hacemos esto a continuación:
v1=1b1+0b2+0b3v2=2b1+2b2+0b3v3=0b13b2+3b3.

De aquí, obtenemos
a11=1,a21=0,a31=0,a12=2,a22=2,a32=0,a13=0,a23=3,a33=3.

Si queremos calcular el determinante, tenemos que considerar las 3!=321=6 permutaciones en S3. Estas permutaciones son

σ1=(123123)σ2=(123132)σ3=(123213)σ4=(123231)σ5=(123321)σ6=(123312).

Los signos de σ1,,σ6 son, como puedes verificar, 1, 1, 1, 1, 1 y 1, respectivamente.

El sumando correspondiente a σ1 es
(1)sign(σ1)a1σ1(1)a2σ1(2)a3σ1(3)(2)=1a11a22a33(3)=1123=6.

El sumando correspondiente a σ2 es
(4)sign(σ2)a1σ2(1)a2σ2(2)a3σ2(3)(5)=(1)a11a23a32(6)=(1)1(3)0=0.

Continuando de esta manera, se puede ver que los sumandos correspondientes a σ1,,σ6 son +6,0,0,+0,0,+0, respectivamente de modo que el determinante es 6.

La expresión de determinante puede parecer algo complicada, pero a través de ella podemos demostrar fácilmente algunos resultados. Consideremos como ejemplo el siguiente resultado.

Proposición. Sea B=(b1,,bn) una base de un espacio vectorial V de dimensión finita n. El determinante de B con respecto a sí mismo es 1.

Demostración. Cuando escribimos a bi en términos de la base b, tenemos que bi=j=1najibj. Como la expresión en una base es única, debemos tener aii=1 y aji=0 si ji. Ahora, veamos qué le sucede al determinante σSnsign(σ)a1σ(1)anσ(n).

Si σ es una permutación tal que σ(i)i para alguna i, entonces en el producto del sumando correspondiente a σ aparece aiσ(i)=0, de modo que ese sumando es cero. En otras palabras, el único sumando no cero es cuando σ es la permutación identidad.

Como el signo de la identidad es 1 y cada aii es 1, tenemos que el determinante es
σSnsign(σ)a1σ(1)anσ(n)=a11ann=11=1.

◻

El determinante es una forma n-lineal alternante

La razón por la cual hablamos de transformaciones n-lineales antisimétricas y alternantes antes de hablar de determinantes es que, en cierto sentido, los determinantes de vectores son las únicas transformaciones de este tipo. Los siguientes resultados formalizan esta intuición.

Teorema. Sea B=(b1,,bn) una base de un espacio vectorial V sobre F. Entonces la transformación det(b1,,bn):VnF es una forma n-lineal y alternante.

Demostración. La observación clave para demostrar este resultado es que det(b1,,bn) se puede reescribir en términos de la base dual b1,,bn. En efecto, recuerda que bi es la forma lineal que «lee» la coordenada de un vector v escrito en la base B. De esta forma,

det(b1,,bn)(v1,,vn)=σSn(sign(σ)j=1nbj(vσ(j)))

Para cada permutación σ, el sumando correspondiente es una forma n-lineal, pues es producto de n formas lineales evaluadas en los distintos vectores. Así que det(b1,,bn) es suma de formas n-lineales y por lo tanto es forma n-lineal.

Para mostrar que el determinante es alternante, tenemos que mostrar que es igual a 0 cuando algún par de sus entradas son iguales. Supongamos que ij y que vi=vj. Tomemos τ a la transposición que intercambia a i y a j. Cuando se compone una permutación con una transposición, su signo cambia. Así, para cualquier permutación σ, tenemos que στ tiene signo diferente.

Además, para cualquier σ tenemos que a1σ(1)anσ(n) y a1στ(1)anστ(n) son iguales, pues vi=vj. Combinando ambas ideas, podemos emparejar a cada sumando del determinante con otro con el cual sume cero. Esto muestra que el determinante es 0.

◻

Usando la teoría que desarrollamos en la entrada anterior, tenemos el siguiente corolario.

Corolario. La forma n-lineal det(b1,,bn) es antisimétrica.

Los determinantes de vectores son las «únicas» formas n-lineales alternantes

Ya vimos que el determinante es una forma n-lineal alternante. Veamos ahora por qué decimos que es «la única». El siguiente resultado dice que cualquier otra forma n-lineal alternante varía de det(b1,,bn) únicamente por un factor multiplicativo.

Teorema. Sea B=(b1,,bn) una base de un espacio vectorial V. Si f:VnF es cualquier forma n-lineal y alternante, entonces f=f(b1,,bn)det(b1,,bn).

Demostración. Para mostrar la igualdad del teorema, que es una igualdad de transformaciones, tenemos que ver que es cierta al evaluar en cualesquiera vectores x1,,xn. Escribamos a cada xi en términos de la base B: xi=j=1naijbj.

Usando la n-linealidad de f en cada una de las entradas, tenemos que
f(x1,,xn)=i=1na1if(bi,x2,,xn)=i,j=1na1ia2if(bi,bj,x3,,xn)==i1,,in=1na1i1aninf(bi1,,bin).

Aquí hay muchos términos, pero la mayoría de ellos son 0. En efecto, si bik=bil, como f es alternante tendríamos que ese sumando es 0. Así, los únicos sumandos que pueden ser no cero son cuando la elección de subíndices es una permutación, es decir cuando existe σ en Sn tal que para ik=σ(k).

Por lo tanto, podemos simplificar la expresión anterior a
f(x1,,xn)=σSna1σ(1)anσ(n)f(bσ(1),,bσ(n)).

Como f es alternante, entonces es antisimétrica. De este modo, podemos continuar la igualdad anterior como
=σSnsign(σ)a1σ(1)anσ(n)f(b1,,bn)=f(b1,,bn)det(b1,,bn)(x1,,xn).

Esto es justo lo que queríamos probar.

◻

Los determinantes de vectores caracterizan bases

Como consecuencia del último teorema de la sección anterior, los determinantes de vectores caracterizan totalmente a los conjuntos de vectores que son bases. A continuación enunciamos esto formalmente.

Corolario. En un espacio vectorial V de dimensión n son equivalentes las siguientes tres afirmaciones para vectores x1,,xn de V:

  1. El determinante de x1,,xn con respecto a toda base es distinto de 0.
  2. El determinante de x1,,xn con respecto a alguna base es distinto de 0.
  3. x1,,xn es una base de V.

Demostración. La afirmación (1) es más fuerte que la (2) y por lo tanto la implica.

Ahora, probemos que la afirmación (2) implica la afirmación (3). Como x1,,xn son n vectores y n es la dimensión de V, para mostrar que forman una base basta mostrar que son linealmente independientes. Anteriormente, vimos que cualquier forma alternante manda vectores linealmente dependientes a 0. Como la hipótesis de (2) es que existe alguna forma alternante que no se anula en x1,,xn, entonces deben ser linealmente independientes y por lo tanto formar una base.

Finalmente, probemos que (3) implica (1). Tomemos B=(b1,,bn) otra base de V. Como det(x1,,xn) es una forma n-lineal, podemos aplicar el teorema anterior y evaluar en x1,,xn para concluir que
det(x1,,xn)(x1,,xn)=det(x1,,xn)(b1,,bn)det(b1,,bn)(x1,,xn).

El término de la izquierda es igual a 1, de modo que ambos factores a la derecha deben ser distintos de 0.

◻

Ejemplo. En el ejemplo que dimos de polinomios vimos que el determinante de 1, 2x y 3x2 con respecto a la base 1, 1+x y 1+x+x2 es igual a 6. De acuerdo al teorema anterior, esto implica que 1, 2x y 3x2 es un conjunto linealmente independiente de polinomios, y de hecho una base.

Además, el teorema anterior también implica que sin importar que otra base B de R2[x] tomemos, el determinante de 1, 2x y 3x2 con respecto a B también será distinto de 0.

Más adelante…

A lo largo de esta entrada estudiamos la definición de determinantes para un conjunto de vectores y enunciamos sus principales propiedades. En las siguientes entradas vamos a hablar cómo se define el determinante para matrices y para transformaciones lineales. Después de las definiciones, pasaremos a estudiar cómo se calculan los determinantes y veremos cómo se aplican a diferentes problemas de álgebra lineal.

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.

  • ¿Cuántos sumandos tendrá el determinante de 5 vectores en un espacio vectorial de dimensión 5 con respecto a cualquier base? Da el número de manera explícita.
  • Verifica que en el primer ejemplo de determinantes de esta entrada, en efecto los sumandos correspondientes a σ1,,σ6 son los que se enuncian.
  • Encuentra el determinante de los vectores (3,1) y (2,4) con respecto a la base ((5,1),(2,3)) de R2.
  • Muestra que los vectores (1,4,5,2), (0,3,2,1), (0,0,1,4) y (0,0,0,1) son linealmente independientes calculando por definición su determinante con respecto a la base canónica de R4.
  • Usa un argumento de determinantes para mostrar que los vectores (1,4,3), (2,2,9), (7,8,27) de R3 no son linealmente independientes. Sugerencia. Calcula su determinante con respecto a la base canónica.

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»

Seminario de Resolución de Problemas: Identidad de Gauss e identidad de suma de cubos

Por Leonardo Ignacio Martínez Sandoval

[latexpage]

Introducción

En la entrada anterior comenzamos a platicar acerca de identidades algebraicas útiles en la resolución de problemas matemáticos. Vimos algunas identidades básicas y platicamos acerca del teorema del binomio de Newton. En esta entrada veremos dos identidades más: la identidad de Gauss para suma de cuadrados y la identidad para factorizar a3+b3+c33abc. Damos más de una demostración de cada una de ellas para seguir explorando ideas algebraicas.

Identidad de cuadrados de Gauss

Proposición. Para a,b,c,d números reales se cumple que (a2+b2)(c2+d2)=(acbd)2+(ad+bc)2.

Demostración 1. Simplemente desarrollamos. Por un lado,
(a2+b2)(c2+d2)=a2c2+a2d2+b2c2+b2d2.

Por otro lado, (acbd)2+(ad+bc)2 es
a2c22abcd+b2d2+a2d2+2abcd+b2c2=a2c2+a2d2+b2c2+b2d2.

◻

La siguiente demostración nos ayuda a entender un poco mejor la identidad y tiene una idea que se puede aplicar en varios contextos.

Demostración 2. Vamos a dar un pequeño brinco a los números complejos, pues ahí podemos hacer la factorización x2+y2=(x+yi)(xyi).

Usando esa identidad:
(a2+b2)(c2+d2)=(a+bi)(abi)(c+di)(cdi)=(a+bi)(c+di)(abi)(cdi)=((acbd)+(ad+bc)i)((acbd)(ad+bc)i)=(acbd)2+(ad+bc)2.

◻

La idea que se puede recuperar de la demostración anterior es la siguiente: a veces una identidad no se puede factorizar en los números reales (racionales, enteros, etc), pero sí en los números complejos (otro sistema numérico más grande). Aunque el problema hable de números reales, es posible que podamos ir a los complejos y regresar a los reales con información.

Problema ejemplo para identidad de Gauss

Problema. Muestra que si tienes un número x de la forma r2+7s2, con r y s números enteros, entonces el número x2020 también es de esa forma.

Sugerencia pre-solución. Aquí, el exponente 2020 es sospechoso, y sugiere que en realidad el problema debe ser más general. Haz algunos casos pequeños para buscar un patrón de cómo se comporta el producto de dos números de esa forma. Después, para estudiar las potencias, usa el principio de inducción.

Solución. Notemos que x=r2+7s2=(r+7si)(r7si) Tomemos otro número de esa forma, digamos y=t2+7u2=(t+7ui)(t7ui). Al hacer el producto de x y y, aparecerá un factor (r+7si)(t+7ui)=((rt7su)+(ru+st)7i) y un factor (r7si)(t7ui)=((rt7su)(ru+st)7i), que multiplicados son iguales a (rt7su)2+7(ru+st)2. Con todo esto, concluimos que el producto de cualesquiera dos números de la forma buscada, también es de la forma buscada. De aquí, x2 es de la forma buscada, e inductivamente xn es de la forma buscada para todo entero n1. En particular, x2020 es de la forma que se quiere.

◻

Identidad para a3+b3+c33abc

Proposición. Para a,b,c números reales, se tiene que a3+b3+c33abc es igual a (a+b+c)(a2+b2+c2abbcca).

Esta identidad también tiene varias demostraciones, que en conjunto guardan varias ideas. Veamos dos de ellas.

Demostración 1. Simplemente hacemos el producto de la segunda expresión para verificar que nos de la primera. Claramente aparece un único a3 y por simétría aparecen b3 y c3 exactamente una vez. También, claramente aparece tres veces la expresión abc. Todas las expresiones que aparecen son cúbicas y ya contamos las «de la forma» x3 y xyz, así que por simetría basta ver qué pasa con cada expresión de la forma x2y. Estas se obtienen ya sea de elegir x en la primera y xy en la segunda, o bien y en la primera y x2 en la segunda, de modo que todas ellas se cancelan.

Sólo para asegurarnos que hicimos todo bien, deberíamos haber contado 36=18 monomios. Hay tres de la forma x3, tres de la forma xyz y cada uno de los seis la forma x2y ya lo encontramos 2 veces, una vez positivo y una vez negativo. Así, nuestra cuenta abarca 3+3+63=18 monomios, así que ya contamos todos los términos.

◻

Hay una segunda demostración, que usa ideas de álgebra lineal. Daremos la idea general, y más adelante, cuando hablemos de matrices y determinantes, platicaremos de estas ideas más a detalle.

Demostración. Calculemos el determinante D de la matriz (abccabbca) de dos formas distintas. Por un lado, podemos sumar los renglones 2 y 3 al primer renglón sin que cambie el determinante, así, D=|a+b+ca+b+ca+b+ccabbca|. De aquí, podemos factorizar a+b+c pues está en cada entrada del primer renglón D=(a+b+c)|111cabbca|.

Finalmente, desarrollando el determinante que queda usando el primer renglón, tenemos que
D=(a+b+c)((a2bc)(cab2)+(c2ab))=(a+b+c)(a2+b2+c2abbcca).

Por otro lado, usando el truco para desarrollar un determinante de 3×3 por diagonales,
D=a3+b3+c3abcabcabc=a3+b3+c33abc.

Igualando ambas expresiones para D, obtenemos la identidad deseada.

◻

Problema ejemplo de factorización de a3+b3+c33abc

Problema. Sean a,b,c números reales. Muestra que a3+b3+c3=3abc si y sólo si a+b+c=0 o a=b=c.

Sugerencia pre-solución. Necesitarás la identidad anterior y un análisis de casos. También, para uno de los casos necesitarás usar la factorización de x22xy+y2 algunas veces.

Solución. De acuerdo a la identidad de la sección anterior, a3+b3+c3=3abc si y sólo si (a+b+c)(a2+b2+c2abbcca)=0.

Notemos que a2+b2+c2abbcca=(ab)2+(bc)2+(ca)22, que siempre es mayor o igual que cero y es igual a 0 si y sólo si ab=bc=ca=0, si y sólo si a=b=c.

Así, a3+b3+c3=3abc si y sólo si alguno de los factores que lo conforman es cero, lo cual pasa si y sólo si a+b+c=0 o a=b=c.

◻

Más problemas

Puedes ver más problemas que usan identidades algebraicas en la entrada anterior de este tema. Además, puedes encontrar más problemas de identidades algebraicas en la Sección 4.1 del libro Problem Solving through Problems de Loren Larson.

Álgebra Superior II: Sistemas de ecuaciones lineales complejos

Por Leonardo Ignacio Martínez Sandoval

Introducción

En la entrada anterior comenzamos a hablar acerca de resolver, en los complejos, ecuaciones de distintos tipos. Además, profundizamos en cómo resolver las ecuaciones cuadráticas complejas. En esta entrada platicaremos acerca de los sistemas de ecuaciones lineales complejos.

Resolveremos a detalle el caso de dos variables y dos ecuaciones. Después, hablaremos un poco acerca de sistemas de ecuaciones con más variables. Un estudio cuidadoso de los sistemas de ecuaciones lineales con más variables se hace en los cursos de álgebra lineal. Un muy buen texto para aprender estos temas es el libro Essential Linear Algebra de Titu Andreescu.

Sistemas de ecuaciones lineales complejos con dos incógnitas

Si a,b son elementos de C y a0, la ecuación lineal ax=b tiene una única solución, dada por x=ba, la cual está bien definida pues todo complejo distinto de 0 tiene inverso multiplicativo.

Si tenemos los números complejos a,b,c,d,e y f, el sistema de ecuaciones lineales en los complejos

ax+by=cdx+ey=f

puede comportarse de tres formas distintas:

  • Su solución existe y es única.
  • Tiene una infinidad de soluciones.
  • No tiene solución.

Si tiene al menos soluciones distintas, tenemos entonces que tiene una infinidad. Cuando la solución del sistema es única, el sistema se puede resolver por los métodos básicos con los que se resuelve un sistema en R:

  • Por substitución: de la primera ecuación se despeja la variable x y su valor se pone en la segunda ecuación. De ahí, obtenemos una ecuación en y. Se despeja y para obtener su valor y con ello se obtiene el valor de x.
  • Igualando coeficientes: multiplicamos la primer ecuación por d y la segunda por a. Al sumar ambas ecuaciones resultantes, queda una ecuación lineal en y.

Ejemplos de sistemas de ecuaciones lineales complejos

Ejemplo 1. Determina todas las soluciones del sistema
2x+iy=3+4iix+5y=94i.

Solución. Para empezar, multiplicamos la segunda ecuación por 2i, de donde obtenemos el sistema
2x+iy=3+4i2x+10iy=8+18i.

Sumando ambas ecuaciones, obtenemos que 11iy=11+22i. Multiplicando por i11 de ambos lados, obtenemos y=2i.

Substituyendo en la segunda ecuación, notamos que 2x=3+4ii(2i)=2+2i, de donde x=1+i. De aquí, la única solución puede ser x=1+i y y=2i, que se puede verificar que en efecto satisfacen la ecuación.

Ejemplo 2. Determina todas las soluciones del sistema
(3+2i)x+iy=3+3i(4+6i)x2y=6+6i.

Solución. Multiplicando la primer ecuación por 2i obtenemos que es equivalente a la ecuación (4i+6i)x2y=6+6i, es decir, ambas ecuaciones difieren sólo por un factor 2i, así que son la misma. Si elegimos cualquier valor de y, podemos encontrar un valor de x que cumpla con la ecuación. Por ejemplo, tomando y=1, de la ecuación obtenemos que x=1. Así, esta ecuación tiene una infinidad de soluciones, dadas por elegir un y y definir x=3+3iiy3+2i.

Ejemplo 3. Determina todas las soluciones del sistema
(1+2i)x+(2+i)y=3+6i3x+3iy=8.

Solución. Supongamos que existe alguna solución para x y y. Multipliquemos la primer ecuación por 3 y la segunda por 1+2i. Obtenemos que
(3+6i)x+(6+3i)y=9+18i(3+6i)x+(6+3i)y=8+16i.

De aquí, 9+18i=8+16i, lo cual es una contradicción. Así, esta ecuación no tiene soluciones.

Método del determinante

Un método más general para resolver sistemas de ecuaciones lineales complejos con dos incógnitas, que nos dice todo lo que puede suceder, es el siguiente. De hecho, exactamente el mismo teorema funciona para R.

Teorema. Sean a,b,c,d,e y f en C. Para el sistema ax+by=cdx+ey=f

definimos a su determinante como el número complejo aebd. Entonces:

  • Si el determinante es distinto de 0, el sistema tiene una solución única para x y y dada por
    x=cebfaebdy=afcdaebd.
  • Si el determinante es 0, entonces el sistema no tiene solución, o tiene una infinidad.

Demostración. Cuando el determinante no es 0, resolvemos el sistema por igualación de coeficientes. Multiplicando la primer ecuación por d, la segunda por a y sumando, obtenemos que (aebd)y=afcd. Como el determinante no es cero, y=afcdaebd. Así mismo, multiplicando la primer ecuación por e, la segunda por b y sumando, obtenemos de manera análoga que x=cebfaebd. Así, si existe una solución, debe tener estos valores. Queda como tarea moral verificar que estos valores cumplen.

Cuando el determinante es 0, tenemos que ae=bd. Si a=b=e=d=0, para que exista una solución se necesita forzosamente que c=f=0, y de hecho en este caso cualquier pareja x,y funciona. Si en este caso alguno de c o f no es 0, el sistema no tiene solución.

Así, continuando el análisis podemos suponer sin pérdida de generalidad que a0. De este modo, e=bda, por lo que la segunda ecuación es equivalente a dx+bday=f, que es adx+bdy=af.

Si d=0, tenemos, de la ecuación anterior, que af=0 y del determinante que ae=bd=0. Como a0, se necesita que e=f=0, de modo que en realidad sólo tenemos una ecuación, la primera. Como a0, podemos elegir cualquier valor de y y de ahí despejar el valor de x, obteniendo una infinidad de soluciones.

Si d0, entonces la ecuación adx+bdy=af es equivalente a la ecuación ax+by=afd. La primer ecuación y esta implican que si hay solución, entonces afd=c. De ser así ,sólo tenemos una ecuación, pero repetida. Por el mismo argumento de arriba, hay una infinidad de soluciones.

◻

Sistemas de ecuaciones lineales complejos con más incógnitas

Los sistemas lineales complejos con más incógnitas se pueden resolver con las mismas técnicas que aquellos en los reales. En cursos como álgebra lineal verás cómo resolver un sistema lineal en general y cómo saber cómo se ven todas sus soluciones. Sin embargo, puedes aprovechar lo que ya sabes del álgebra de los complejos para resolver distintos sistemas lineales.

Problema 1. Resuelve en los complejos el sistema de ecuaciones

3a+(2+i)b+(1+2i)c=1+i3b+(2+i)c=2+2i3c=3+3i.

Solución. Resolvemos el sistema por substitución. Nos conviene empezar con la tercer ecuación, que tiene únicamente una variable. De ella obtenemos que c=1+i. Substituyendo en la segunda ecuación, obtenemos que 3b+(2+i)(1+i)=2+2i, de donde 3b+1+3i=2+2i, así que 3b=1i, entonces b=1313i.

Con los valores de b y c podemos substituir en la primer ecuación. Notando que
(2+i)(1313i)=113i(1+2i)(1+i)=1+3i(1+i)(113i)(1+3i)=153i,

obtenemos que a=1359i.

En resumen,
a=1359ib=1313ic=1+i

es la única posible solución, y se puede mostrar que en efecto satisface las tres ecuaciones.

Problema 2. Resuelve en los complejos el sistema de ecuaciones

(1+5i)a+b+c+d+e=2a+(1+5i)b+c+d+e=2a+b+(1+5i)c+d+e=2a+b+c+(1+5i)d+e=2a+b+c+d+(1+5i)e=2.

Solución. Sumando todas las ecuaciones, tenemos que (5+5i)(a+b+c+d+e)=10, de donde obtenemos que
a+b+c+d+e=21+i=1i.

De la primera ecuación, obtenemos que 2=(a+b+c+d+e)+5ia=1i+5ia, por lo que a=1+i5i=1515i. Por simetría, el resto de las variables también tiene este valor, de modo que a=b=c=d=e=1515i es la única solución.

Más adelante…

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.

  1. Verifica que las soluciones de los ejemplos de sistemas de ecuaciones lineales complejos de dos variables en efecto son soluciones.
  2. Resuelve en los complejos el sistema de ecuaciones 2x+(1+i)y=4(5i)x+(3+2i)y=0.
  3. En el teorema del método del determinante, cuando el determinante no es cero, encontramos una solución. Verifica que en efecto satisface el sistema original.
  4. Verifica que las soluciones de los ejemplos en varias variables en efecto satisfacen el sistema original.
  5. Resuelve en los complejos el sistema de ecuaciones x+(1+i)y=4y+(2+i)z=5z+(3+i)x=6.

Entradas relacionadas

Agradecimientos

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