Archivo de la etiqueta: sylvester

Álgebra Lineal I: Formas cuadráticas, propiedades, polarización y Gauss

Por Leonardo Ignacio Martínez Sandoval

Introducción

En la entrada anterior hablamos acerca de formas bilineales y comenzamos a hablar de formas cuadráticas. Discutimos cómo a partir de estas nociones a la larga podremos hablar de geometría y cálculo en espacios vectoriales. El objetivo de esta entrada es entender mejor a las formas cuadráticas y su relación con formas bilineales.

Lo primero que haremos es demostrar la identidad de polarización, que a grandes rasgos dice que hay una biyección entre las formas bilineales simétricas y las formas cuadráticas. Veremos algunos ejemplos concretos de esta biyección. A partir de ella demostraremos algunas propiedades de formas cuadráticas. Finalmente, hablaremos brevemente de un bello resultado de Gauss que caracteriza las formas cuadráticas en Rn en términos de formas lineales, de las cuales discutimos mucho cuando hablamos de espacio dual.

Como pequeño recordatorio de la entrada anterior, una forma bilineal de un espacio vectorial V es una transformación b:V×VR tal que cada que fijamos una coordenada, es lineal en la otra. Esta forma es simétrica si b(x,y)=b(y,x) para cada par de vectores x,y en V. Una forma cuadrática de V es una transformación q:VR tal que q(x)=b(x,x) para alguna forma bilineal b.

Formas cuadráticas y polarización

En la entrada anterior enunciamos el siguiente teorema, que mostraremos ahora.

Teorema (identidad de polarización). Sea q:VR una forma cuadrática. Existe una única forma bilineal simétrica b:V×VR tal que q(x)=b(x,x) para todo vector x. Esta forma bilineal está determinada mediante la identidad de polarización b(x,y)=q(x+y)q(x)q(y)2.

Demostración. Tomemos una forma cuadrática q de V. Por definición, está inducida por una forma bilineal B de V, es decir, q(x)=B(x,x). Definamos la transformación b mediante b(x,y)=q(x+y)q(x)q(y)2. Comencemos probando que b es una transformación bilineal simétrica. Notemos que:
b(x,y)=q(x+y)q(x)q(y)2=B(x+y,x+y)B(x,x)B(y,y)2=B(x,x)+B(x,y)+B(y,x)+B(y,y)B(x,x)B(y,y)2=B(x,y)+B(y,x)2.

De aquí es muy claro que b es forma bilineal, pues fijando x, set tiene que b(x,y) es combinación lineal de dos formas lineales en y; y fijando y, se tiene que b(x,y) es combinación lineal de dos formas lineales en x. Además, de esta igualdad (o directo de la definición de b) es claro que b(x,y)=b(y,x).

También de esta igualdad obtenemos que b(x,x)=B(x,x)=q(x).

Para mostrar la unicidad, notemos que cualquier forma bilineal simétrica b tal que b(x,x)=q(x) debe satisfacer, como en las cuentas que hicimos arriba, que
q(x+y)=b(x+y,x+y)=q(x)+q(y)+b(x,y)+b(y,x)=q(x)+q(y)+2b(x,y).

De aquí, despejando b, se obtiene que debe tener la forma de b.

◻

El teorema anterior justifica la siguiente definición.

Definición. Dada una forma cuadrática q de V, a la única forma bilineal simétrica b de V tal que q(x)=b(x,x) le llamamos la forma polar de q.

Ejemplo 1. En el espacio vectorial Rn, la transformación q:RnR dada por q(x1,,xn)=x12++xn2. es una forma cuadrática. Su forma polar es la forma bilineal producto punto que manda a x=(x1,,xn) y y=(y1,,yn) a b(x,y)=x1y1++xnyn.

Esto coincide con la construcción dada por la identidad de polarización, ya que q(x+y)q(x)q(y)=i=1n(xi+yi)2xi2yi2=i=1nxiyi

Ejemplo 2. En el espacio vectorial R[x] de polinomios con coeficientes reales, la transformación Q dada por Q(p)=p(0)p(1)+p(2)2 es una forma cuadrática. Para encontrar a su forma bilineal polar, usamos la identidad de polarización
B(p,q)=Q(p+q)Q(p)Q(q)2=(p+q)(0)(p+q)(1)+(p+q)(2)2p(0)p(1)p(2)2q(0)q(1)q(2)22=p(0)q(1)+q(0)p(1)+2p(2)q(2)2=p(0)q(1)2+p(1)q(0)2+p(2)q(2).

Propiedades de formas cuadráticas

Si q es una forma cuadrática, x es un vector y c es un real, tenemos que q(cx)=c2q(x), pues sale una c por cada una de las coordenadas de la forma bilineal asociada. En particular, q(x)=q(x).

La identidad de polarización nos permite probar otras propiedades de formas bilineales y formas cuadráticas.

Proposición. Sea q una forma cuadrática en V con forma polar b. Entonces:

  • Para todo par de vectores x y y en V, se tiene que b(x,y)=q(x+y)q(xy)4.
  • (Ley del paralelogramo) Para todo par de vectores x y y en V, se tiene que q(x+y)+q(xy)=2(q(x)+q(y)).
  • (Teorema de Pitágoras) Para vectores x y y tales que b(x,y)=0, se tiene que q(x+y)=q(x)+q(y).
  • (Diferencia de cuadrados) Para todo par de vectores x y y en V, se tiene que b(x+y,xy)=q(x)q(y).

Demostración. Por la identidad de polarización tenemos que b(x,y)=q(x+y)q(x)q(y)2, y como q(y)=q(y), tenemos también por la identidad de polarización que b(x,y)=b(x,y)=q(xy)q(x)q(y)2.

Restando la segunda ecuación de la primera, obtenemos la primer propiedad. Sumando ambas obtenemos la ley del paralelogramo.

El teorema de Pitágoras es una consecuencia directa de la identidad de polarización.

La identidad de diferencia de cuadrados es una consecuencia de la primer propiedad aplicada a los vectores x+y y xy, y de usar que q(2x)=4q(x) y que q(2y)=4q(y).

◻

Forma de las formas cuadráticas

Otra consecuencia de la identidad de polarización es que establece una biyección entre las formas cuadráticas y las formas simétricas bilineales. Esta asociación nos permite decir cómo se ven exactamente las formas cuadráticas en espacios vectoriales de dimensión finita.

Toda forma cuadrática viene de una forma bilineal simétrica. En la entrada anterior, mencionamos que para definir una forma bilineal simétrica en un espacio vectorial V de dimensión n, basta tomar una base {e1,,en} de V y decidir los valores bij de b(ei,ej) para 1ijn. Como b es simétrica, para j<i se tendría que b(ei,ej)=b(ej,ei), es decir, que bji=bij.

De esta forma, para todo vector v en V podemos encontrar el valor de q(v) expresando v en la base {e1,,en}, digamos, v=a1e1++anen, de donde q(v)=i=1nbiiai2+21i<jnbijaiaj.

Ejemplo. Toda forma cuadrática en R3 se obtiene de elegir reales a,b,c,d,e,f y definir q(x,y,z)=ax2+by2+cz2+2dxy+2eyz+2fzx. La forma polar de q es la forma bilineal B tal que para la base canónica e1,e2,e3 de R3 hace lo siguiente

B(e1,e1)=aB(e2,e2)=bB(e3,e3)=cB(e1,e2)=B(e2,e1)=dB(e2,e3)=B(e3,e2)=eB(e3,e1)=B(e1,e3)=f.

Teorema de Gauss de formas cuadráticas (opcional)

Para esta sección, fijemos al espacio vectorial como Rn. Hay una forma muy natural de construir formas cuadráticas a partir de formas lineales. Tomemos números reales α1,,αr y formas lineales l1,,lr. Consideremos q(x)=α1l1(x)2++αrlr(x)2. Se tiene que q es una forma cuadrática. La demostración de ello es sencillo y se queda como tarea moral.

Lo que descubrió Gauss es que todas las formas cuadráticas se pueden expresar de esta forma, y de hecho, es posible hacerlo usando únicamente formas lineales que sean linealmente independientes y coeficientes 1 y 1.

Teorema (clasificación de Gauss de formas cuadráticas). Sea q una forma cuadrática en Rn. Entonces, existen enteros no negativos r y s, y formas lineares l1,,lr,m1,,ms en (Rn), todas ellas linealmente independientes, tales que q=l12++lr2m12ms2.

Hay un pequeño refinamiento de este teorema, demostrado por Sylvester.

Teorema (teorema de la inercia de Sylverster). Los números r y s en el teorema de clasificación de Gauss de formas cuadráticas son únicos.

Ejemplo. Tomemos la forma cuadrática en R3 dada por q(x,y,z)=xy+yz+zx. Por el teorema de Gauss, esta forma se debe de poder poner como combinación lineal de cuadrados de formas lineales independientes. En efecto, tenemos que: xy+yz+zx=(2x+y+z2)2(yz2)2x2, en donde
(x,y,z)2x+y+z2,(x,y,z)yz2 y (x,y,z)x
son formas lineales linealmente independientes.

Más adelante…

En esta entrada estudiamos a fondo la identidad de polarización; esto nos permitió concluir que existe una biyección entre las funciones bilineales simétricas y las formas cuadráticas. También, pusimos mucho énfasis en ejemplos concretos de esta biyección.

Con esto estamos listos para empezar a pensar en cómo haríamos geometría o cálculo en espacios vectoriales. Abordaremos estos temas al final de esta unidad. En la siguiente entrada hablaremos del producto interior.

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.

  • Verifica que las formas cuadráticas de los ejemplos del teorema de polarización en efecto son formas cuadráticas.
  • Muestra que q(x,y)=3x2y2+7y no es una forma cuadrática.
  • Muestra que si α1,,αr son reales y tomamos formas lineales l1,,lr en Rn, entonces q(x)=a1l1(x)2++αrlr(x)2 es una forma cuadrática.
  • ¿Quién es la forma polar de la forma cuadrática Q(f)=01f2(x)dx en el espacio vectorial de funciones continuas en el intervalo [0,1]?

Una demostración algorítmica del teorema de Gauss se puede encontrar en la Sección 10.1 del libro de Álgebra Lineal de Titu Andreescu.

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: Problemas de rango de transformaciones y matrices.

Por Ayax Calderón

Introducción

Con anterioridad vimos el concepto de rango de una matriz y rango de una transformación lineal, además del muy importante teorema de rango-nulidad y la desigualdad de Sylvester. Vimos también, como contenido optativo, el versátil teorema de la factorización PJQ. En esta ocasión nos enfocaremos en resolver problemas de rango que nos servirán para repasar dichos conceptos.

Problemas resueltos

Problema 1. Encuentra el kernel y el rango de la transformación lineal T:R2[x]R3[x] definida por T(f(x))=2f(x)+0x3f(t)dt.

Antes de comenzar a leer la solución, es conveniente que te convenzas de que T es una transformación lineal y que está bien definida, es decir, que en efecto toma un polinomio de grado a lo más dos con coeficientes reales y lo lleva a un polinomio de grado a lo más tres con coeficientes reales.

Solución. Consideremos B={1,x,x2} la base canónica de R2[x].
Entonces
Im(T)=span({T(1),T(x),T(x2)})=span({3x,2+32x2,4x+x3}).

Para determinar el rango de ImT, colocamos a las coordenadas de estas imágenes en la siguiente matriz A,

A=(0300203200401)

y con el algoritmo de reducción gaussiana llegamos a que

Ared=(1034001000001)

Como Ared tiene 3 pivotes se sigue que rank(T)=3.

Luego, por el teorema de rango nulidad se tiene que

3=dim(R2[x])=dim(ker(T))+rank(T)=dim(ker(T))+3.

Así, dim(ker(T))=0, por lo tanto ker(T)={0}.

La desigualdad de Sylvester nos ayuda a acotar el rango de una suma de matrices por abajo. La desigualdad rank(A+B)rank(A)+rank(B) nos ayuda a acotarlo por arriba. Combinar ambas ideas puede ser útil en problemas de rango de matrices.

Problema 2. Sea AMn(C) una matriz idempotente. Prueba que rank(A)+rank(InA)=n.

Recuerda que una matriz es idempotente si A2=A.

Solución. Como A2=A, entonces A(InA)=On.
Luego, por la desigualdad de Sylvester se tiene que
0=rank(On)=rank(A(InA))rank(A)+rank(InA)n,

entonces rank(A)+rank(InA)n.

Por otro lado, como para cualesquiera matrices X,Y se tiene
rank(X+Y)rank(X)+rank(Y), entonces
n=rank(In)rank(A)+rank(InA),
de modo que nrank(A)+rank(InA).

Combinando ambas desigualdades, rank(A)+rank(InA)=n.

◻

Problema 3. Encuentra el rango de la transformación lineal T:R2[x]M2(R) definida por
T(f(x))=(f(1)f(2)00f(0)).

Solución. Para determinar el rango, basta tomar una base, encontrar la imagen de sus elementos bajo T y determinar cuántos de estos elementos son linealmente independientes. Considera B={1,x,x2} la base canónica de R2[x]. Tenemos que

Im(T)=span(T(B))=span({T(1),T(x),T(x2)})=span({(0001),(1000),(3000)})=span({(0001),(1000)}).

Notemos también que Extra close brace or missing open brace es linealmente independiente.

Por lo tanto C es una base para Im(T) y así rank(T)=2.

Problema 4. Sean AM3,2(R) y BM2,3(R) matrices tales que
AB=(224134123)

Muestra que BA es la identidad.

El enunciado no parece mostrar que este sea uno de los problemas de rango de matrices. Sin embargo, para poder resolverlo usaremos las herramientas que hemos desarrollado hasta ahora.

Partiremos el problema en los siguientes pasos.

  1. Verificar que (AB)2=AB y que rank(AB)=2.
  2. Probar que BA es invertible.
  3. Probar que (BA)3=(BA)2 y deducir que BA=I2.

Solución.

1. Realizamos la operación matricial:

(224134123)(224134123)=(224134123)

Ahora, aplicando reducción gaussiana en AB obtenemos que (AB)red=(101011000).

Como (AB)red tiene sólo dos pivotes, entonces rank(AB)=2.

2. Usando la desigualdad de rango para producto de matrices, obtenemos que
rank(BA)rank(A(BA)B)=rank((AB)2)=rank(AB)=2.

Entonces, rank(BA)2. Por otro lado, como BAM2(R), entonces rank(BA)2. Así, rank(BA)=2 y BA es una matriz en M2(R), así que es invertible.

3. Como (AB)2=AB, entonces B(AB)2A=B(AB)A=(BA)2. Por consiguiente BABABA=(BA)2 y así (BA)3=(BA)2 y como BA es invertible, podemos multiplicar en ambos lados de esta última igualdad por ((BA)1)2 para obtener BA=I2.

◻

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: Rango de transformaciones lineales y matrices

Por Leonardo Ignacio Martínez Sandoval

Introducción

En entradas anteriores hablamos de transformaciones lineales, cómo actúan en conjuntos especiales de vectores y de cómo se pueden representar con matrices. Hablamos también de cómo cambiar de una base a otra y cómo usar esto para entender transformaciones en varias bases. Estamos listos para introducir un concepto fundamental de álgebra lineal, el de rango de una transformación lineal y de una matriz.

Antes de entrar en las definiciones formales, vale la pena hablar un poco de rango de manera intuitiva. Supongamos que V es un espacio vectorial de dimensión n y que W es un espacio vectorial sobre el mismo campo que V. Una transformación lineal T:VW puede «guardar mucha independencia lineal» o «muy poquita». Si T es inyectiva, ya vimos antes que T manda linealmente independientes a linealmente independientes. Si T es la transformación 0, entonces se «pierde toda la independencia».

El rango mide algo intermedio entre estos dos extremos. Mientras mayor sea el rango, más independencia lineal se preserva y viceversa. Si mantienes esta intuición en mente, varias de las proposiciones te resultarán más naturales.

Otro buen ejemplo para tener en mente es tomar una transformación lineal T:R3R3. Si es la transformación identidad, la base canónica se preserva. Si es la proyección al plano xy, entonces «perdemos» al vector (0,0,1), pues se va al (0,0,0). Si es la proyección al eje x, «perdemos» al (0,1,0) y al (0,0,1) pues ambos se van a (0,0,0). Y si es la transformación 0, perdemos a todos. El rango precisamente va a medir esto, y para estos ejemplos tendremos rango 3, 2, 1 y 0 respectivamente.

Rango para transformaciones lineales

Como en otras ocasiones, cuando hablemos de transformaciones lineales entre espacios vectoriales, serán sobre un mismo campo F.

Definición. Sean V y W espacios de dimensión finita. El rango de una transformación lineal T:VW es la dimensión de la imagen de T, es decir, rank(T)=dimImT.

Si B es una base de V, entonces genera a V. La transformación T es suprayectiva de V a ImT, de modo que T(B) es generador de ImT. De esta forma, para encontrar el rango de una transformación lineal T:VW basta:

  • Tomar una base B de V.
  • Aplicar T a cada elemento de B.
  • Determinar un conjunto linealmente independiente máximo en T(B).

Para hacer este último paso, podemos poner a los vectores coordenada de T(B) con respecto a una base de W como los vectores fila de una matriz A y usar reducción gaussiana. Las operaciones elementales no cambian el espacio generado por las filas, así que el rango de T es el número de vectores fila no cero en la forma escalonada reducida Ared de A.

Ejemplo. Encuentra el rango de la transformación lineal T:R3M2(R) que manda (x,y,z) a (x+yz2x2y2zx+zy).

Solución. Tomemos e1,e2,e3 la base canónica de R3. Tenemos que T(e1)=(1201), T(e2)=(1021) y T(e3)=(1021).

Tomando la base canónica E11,E12,E21,E22 de M2(R), podemos entonces poner a las coordenadas de T(e1),T(e2),T(e2) como vectores fila de una matriz (120110211021). Sumando la segunda fila a la tercera, y después restando la primera a la segunda,obtenemos la matriz (120102220000). De aquí, sin necesidad de terminar la reducción gaussiana, podemos ver que habrá exactamente dos filas no cero. De este modo, el rango de la transformación es 2.

Propiedades del rango

Demostremos ahora algunas propiedades teóricas importantes acerca del rango de una transfromación lineal.

Proposición. Sean U, V y W espacios de dimensión finita. Sean S:UV, T:VW, T:VW transformaciones lineales. Entonces:

  1. rank(T)dimV
  2. rank(T)dimW
  3. rank(TS)rank(T)
  4. rank(TS)rank(S)
  5. rank(T+T)rank(T)+rank(T)

Demostración. (1) Pensemos a T como una transformación T:VIm(T). Haciendo esto, T resulta ser suprayectiva, y por un resultado anterior tenemos que dimVdimImT=rank(T).

(2) Sabemos que Im(T) es un subespacio de W, así que rank(T)=dimImTdimW.

(3) La imagen de T contiene a la imagen de TS, pues cada vector de la forma T(S(v)) es de la forma T(w) (para w=S(v)). Así, rank(T)=dimImTdimImTS=rank(TS).

(4) La función TS coincide con la restricción TImS de T a ImS. Por el inciso (1), rank(TImS)dimImS=rank(S), así que rank(TS)rank(S).

(5) Tenemos que Im(T+T)ImT+ImT. Además, por un corolario de la fórmula de Grassman, sabemos que
dim(ImT+ImT)dimImT+dimImT=rank(T)+rank(T).

Así,
rank(T+T)rank(ImT+ImT)rank(T)+rank(T).

◻

Proposición. Sean R:UV, T:VW y S:WZ transformaciones lineales con R suprayectiva y S inyectiva. Entonces rank(STR)=rank(T).

Dicho de otra forma «composición por la izquierda con transformaciones inyectivas no cambia el rango» y «composición por la derecha con transformaciones suprayectivas no cambia el rango». Un corolario es «composición con transformaciones invertibles no cambia el rango».

Demostración. De la proposición anterior, tenemos que rank(ST)rank(T). La restricción SImT de S a la imagen de T es una transformación lineal de ImT a Im(ST) que es inyectiva, de modo que dimImTdimIm(ST), que es justo rank(T)rank(ST), de modo que tenemos la igualdad rank(ST)=rank(T).

Como R es suprayectiva, ImR=V, de modo que Im(STR)=Im(ST). Así, rank(STR)=rank(ST)=rank(T).

◻

Teorema de rango-nulidad

Una transformación lineal T:VW determina automáticamente dos subespacios de manera natural: el kernel kerT y la imagen ImT. Resulta que las dimensiones de kerT, de ImT y de V están fuertemente relacionadas entre sí.

Teorema. Sean V y W espacios de dimensión finita. Sea T:VW una transformación lineal. Entonces dimkerT+rank(T)=dimV.

Demostración. Supongamos que dimV=n y dimkerT=k. Queremos mostrar que rank(T)=nk. Para ello, tomemos una base B de kerT y tomemos B={v1,,vnk} tal que BB sea base de V. Basta mostrar que T(B)={T(v1),,T(vnk)}ImT es base de ImT. Sea U el generado por B, de modo que V=UkerT.

Veamos que T(B) es generador de ImT. Tomemos T(v) en ImT. Podemos escribir v=z+u con zkerT y uU. Así, T(v)=T(z)+T(u)=T(u), y este último está en el generado por T(B).

Ahora veamos que T(B) es linealmente independiente. Si α1T(v1)++αnkT(vnk)=0, entonces T(α1v1++αnkvnk)=0, de modo que α1v1++αnkvnk está en U y en kerT, pero la intersección de estos espacios es {0}. Como esta combinación lineal es 0 y B es linealmente independiente, α1==αn=0.

De esta forma, T(B) es linealmente independiente y genera a ImT, de modo que rank(T)=|B|=nk.

◻

Ejemplo. Consideremos de nuevo la transformación lineal T:R3M2(R) que manda (x,y,z) a (x+yz2x2y2zx+zy). Muestra que T no es inyectiva.

Solución. Ya determinamos previamente que esta transformación tiene rango 2. Por el teorema de rango-nulidad, su kernel tiene dimensión 1. Así, hay un vector v(0,0,0) en el kernel, para el cual T(v)=0=T(0), de modo que T no es inyectiva.

◻

Problema. Demuestra que para cualquier entero n existe una terna (a,b,c)(0,0,0) con a+b+c=0 y tal que 01at2n+btn+cdt=0.

Solución. Podríamos hacer la integral y plantear dos ecuaciones lineales. Sin embargo, daremos argumentos dimensionales para evitar la integral. Consideremos las transformaciones lineales T:R3R y S:R3R dadas por
T(x,y,z)=01xt2n+ytn+zdtS(x,y,z)=x+y+z.
Notemos que T(0,0,1)=011dt=1=S(0,0,1), de modo que ni T ni S son la transformación 0. Como su rango puede ser a lo más dimR=1, entonces su rango es 1. Por el teorema de rango-nulidad, dimkerS=dimkerT=2. Como ambos son subespacios de R3, es imposible que kerSkerT={0}, de modo que existe (a,b,c) no cero tal que T(a,b,c)=S(a,b,c)=0. Esto es justo lo que buscábamos.

◻

Rango para matrices

Definición. El rango de una matriz A en Mm,n(F) es el rango de la transformación lineal asociada de Fn a Fm dada por XAX. Lo denotamos por rank(A).

A partir de esta definición y de las propiedades de rango para transformaciones lineales obtenemos directamente las siguientes propiedades para rango de matrices.

Proposición. Sean m, n y p enteros. Sea B una matriz en Mn,p(F) y A, A matrices en Mm,n(F). Sea P una matriz en Mn,p(F) cuya transformación lineal asociada es suprayectiva y Q una matriz en Mr,m(F) cuya transformación lineal asociada es inyectiva. Entonces:

  1. rank(A)min(m,n)
  2. rank(AB)min(rank(A),rank(B))
  3. rank(A+A)rank(A)+rank(A)
  4. rank(QAP)=rank(A)

Como discutimos anteriormente, el rango de una transformación se puede obtener aplicando la transformación a una base y viendo cuál es el máximo subconjunto de imágenes de elementos de la base que sea linealmente independiente. Si tomamos una matriz A en Mm,n(F), podemos aplicar esta idea con los vectores e1,,en de la base canónica de Fn. Como hemos visto con anterioridad, para cada i=1,,n tenemos que el vector Aei es exactamente la i-ésima columna de A. Esto nos permite determinar el rango de una matriz en términos de sus vectores columna.

Proposición. El rango de una matriz en Mm,n(F) es igual a la dimensión del subespacio de Fm generado por sus vectores columna.

Problema. Determina el rango de la matriz (310500829001042).

Solución. Como es una matriz con 3 filas, el rango es a lo más 3. Notemos que entre las columnas están los vectores (3,0,0), (0,2,0) y (0,0,2), que son linealmente independientes. De esta forma, el rango de la matriz es 3.

A veces queremos ver que el rango de un producto de matrices es grande. Una herramienta que puede servir en estos casos es la desigualdad de Sylvester.

Problema (Desigualdad de Sylvester). Muestra que para todas las matrices A, B en Mn(F) se tiene que rank(AB)rank(A)+rank(B)n.

Solución. Tomemos T1:FnFn y T2:FnFn tales que T1(X)=AX y T2(X)=BX. Lo que tenemos que probar es que rank(T1T2)rank(T1)+rank(T2)n.

Consideremos S1 como la restricción de T1 a ImT2. Tenemos que kerS1kerT1, así que dimkerS1dimkerT1. Por el teorema de rango-nulidad en S1, tenemos que
rank(T2)=dimImT2=dimkerS1+rank(S1)=dimkerS1+rank(T1T2)dimkerT1+rank(T1T2), así que rank(T2)dimkerT1+rank(T1T2).

Por el teorema de rango-nulidad en T1 tenemos que dimkerT1+rank(T1)=n.

Sumando la desigualdad anterior con esta igualdad obtenemos el resultado.

◻

El teorema PJQ (opcional)

El siguiente resultado no se encuentra en el temario usual de Álgebra Lineal I. Si bien no formará parte de la evaluación del curso, recomendamos fuertemente conocerlo y acostumbrarse a usarlo pues tiene amplias aplicaciones a través del álgebra lineal.

Teorema (Teorema PJQ). Sea A una matriz en Mm,n(F) y r un entero en {0,,min(m,n)}. El rango de A es igual a r si y sólo si existen matrices invertibles PMm(F) y QMn(F) tales que A=PJrQ, en donde Jr es la matriz en Mm,n cuyas primeras r entradas de su diagonal principal son 1 y todas las demás entradas son cero, es decir, en términos de matrices de bloque, Jr=(Ir000).

No damos la demostración aquí. Se puede encontrar en el libro de Titu Andreescu, Teorema 5.68. Veamos algunas aplicaciones de este teorema.

Problema 1. Muestra que una matriz tiene el mismo rango que su transpuesta.

Solución. Llamemos r al rango de A. Escribimos A=PJrQ usando el teorema PJQ, con P y Q matrices invertibles. Tenemos que tA=tQtJrtP, con tQ y tP matrices invertibles. Además, tJr es de nuevo de la forma de Jr. Así, por el teorema PJQ, tenemos que tA es de rango r.

Combinando el problema anterior con el resultado del rango de una matriz en términos de sus vectores columna obtenemos lo siguiente.

Proposición. El rango de una matriz en Mm,n(F) es igual a la dimensión del subespacio de Fn generado por sus vectores renglón.

Terminamos esta entrada con una aplicación más del teorema PJQ.

Problema 2. Muestra que una matriz A de rango r se puede escribir como suma de r matrices de rango 1. Muestra que es imposible hacerlo con menos matrices.

Solución. Expresamos A=PJrQ usando el teorema PJQ. Si definimos Ai=PEiiQ para i=1,,r, donde Eii es la matriz cuya entrada (i,i) es uno y las demás cero, claramente tenemos que Jr=E11+E22++Err, por lo que A=PJrQ=A1+A2++Ar. Además, como Eii es de rango 1, por el teorema PJQ cada matriz Ai es de rango 1.

Veamos que es imposible con menos. Si B1,,Bs son matrices de rango 1, como el rango es subaditivo tenemos que rank(B1++Bs)s. Así, si sumamos menos de r matrices, no podemos obtener a A.

◻

Más adelante…

Esta entrada es solamente una breve introducción al concepto de rango y a algunas propiedades que pueden ser de utilidad al momento de calcular el rango de una matriz o una transformación lineal. Más adelante, veremos que el rango de una matriz está también relacionado con las soluciones de su sistema lineal homogéneo asociado.

El teorema de rango-nulidad es fundamental para el álgebra lineal. Muchas veces necesitamos calcular el rango de la imagen de una transformación lineal, pero es mucho más fácil calcular la dimensión de su kernel. O viceversa. En estas situaciones es muy importante recordar la forma en la que dicho teorema las relaciona.

Con este tema termina la segunda unidad del curso. Ahora estudiaremos aspectos un poco más geométricos de espacios vectoriales. En la siguiente unidad, hablaremos de dualidad, ortogonalidad, formas bilineales y productos interiores.

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.

  • Termina de hacer la reducción gaussiana del primer ejemplo.
  • Sea T una transformación de un espacio vectorial V de dimensión finita a si mismo. Usa el teorema de rango-nulidad para mostrar que si T es inyectiva o suprayectiva, entonces es biyectiva.
  • Determina el rango de la matriz (000837812031449).
  • Demuestra que aplicar operaciones elementales a una matriz no cambia su rango.
  • Demuestra que matrices similares tienen el mismo rango.
  • Demuestra por inducción que para matrices A1,,An del mismo tamaño tenemos que rank(A1++An)i=1nrank(Ai).
  • Escribe la demostración de la última proposición de la sección del teorema PJQ
  • Revisa la demostración del teorema de descomposición PJQ en el libro de Titu Andreescu.

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»