Archivo de la etiqueta: polinomio

Álgebra Lineal II: Polinomio característico de familias especiales

Por Julio Sampietro

Introducción

En la entrada anterior dimos la definición de polinomio característico. Vimos que siempre es un polinomio mónico y que su grado es exactamente del tamaño de la matriz. También, vimos cómo calcular el polinomio mínimo en algunos casos particulares. En esta entrada veremos varias propiedades que nos van a facilitar el calcular el polinomio característico (y por tanto los eigenvalores) en un amplio rango de matrices diferentes.

Comenzaremos estudiando el polinomio mínimo de las triangulares superiores. Luego, veremos cómo calcular el polinomio de matrices nilpotentes. No solo nos harán la vida más fácil los resultados a continuación, si no que los usaremos en la teoría más adelante.

Matrices triangulares superiores y transpuestas

El caso de las matrices triangulares superiores es muy sencillo, como veremos a través del siguiente problema.

Problema. Sea A=[aij] una matriz triangular superior. Demuestra que

χA(X)=i=1n(Xaii).

Solución. La matriz XInA sigue siendo triangular superior, y sus entradas diagonales son precisamente Xaii. Usando que el determinante de una matriz triangular superior es el producto de sus entradas diagonales y usando la definición se sigue que

χA(X)=det(XInA)=i=1n(Xaii).

◻

Ejemplo. Si queremos calcular el polinomio característico de la matriz

A=(1π2021010003).

entonces podemos aplicar el problema anterior y deducir inmediatamente que

χA(X)=(X1)(X+2)(X3).

¡Qué complicado hubiera sido calcular el determinante a pie!

Por otro lado, recordando la demostración que dice que los eigenvalores de la transpuesta de una matriz son iguales a los de la matriz original era de esperarse que el polinomio característico también «se portara bien» bajo transposición.

Problema. Demuestra que las matrices A y tA tienen el mismo polinomio característico para cualquier AMn(F).

Solución. Notamos que t(XInA)=XIn tA. Como una matriz y su transpuesta tienen el mismo determinante se tiene que

χA(X)=det(XInA)=det( t(XInA))=det(XIn tA)=χtA(X).

◻

Estrictamente hablando, estamos haciendo un poquito de trampa en la demostración anterior (y de hecho en varias que involucran a la variable X). Las propiedades de determinantes que hemos visto (como que una matriz y su transpuesta tienen el mismo determinante) las obtuvimos partiendo de la hipótesis de que las entradas vienen de un campo F. Pero cuando agregamos a la variable X, ahora las entradas vienen más bien de un anillo: el anillo de polinomios en F[X]. Aunque esto parezca un problema, en realidad no lo es. Las propiedades que usamos pueden mostrarse también en ese contexto.

Veamos ahora cómo podemos aplicar el resultado anterior en un ejemplo concreto.

Ejemplo. Queremos calcular el polinomio característico de la matriz

A=(000490112.)

Para esto notamos que

tA=(041091002)

que es triangular superior. Usando el primer problema

χtA(X)=X(X9)(X2).

Finalmente por el último problema χA(X)=χtA(X)=X(X9)(X2).

El término de la traza

Como vimos en la entrada anterior, en el polinomio det(XA+B) aparecen los términos det(A) y det(B). El siguiente problema aplica esto al polinomio característico e incluso deducimos otro término: la traza.

Problema. Demuestra que el polinomio característico de AMn(F) es de la forma

χA(X)=XnTr(A)Xn1++(1)ndetA.

Solución. Regresemos a la definición

det(XInA)=σSnsign(σ)(Xδ1σ(1)a1σ(1))(Xδnσ(n)anσ(n)).

Haciendo la expansión salvajemente podemos recuperar al menos los primeros términos de (Xδ1σ(1)a1σ(1))(Xδnσ(n)anσ(n)), que son Xni=1nδiσ(i)Xn1j=1n(kjδkσ(k))ajσ(j)+.

Más aún, nota cómo el producto j=1nδjσ(j) es distinto de cero si y sólo si j=σ(j) para todo j: es decir si σ es la identidad. Esto muestra que χA(X) es mónico de grado n, como ya habíamos mencionado en la entrada anterior.

Además, el término constante está dado por χA(0)=det(0InA)=det(A)=(1)ndet(A). Alternativamente pudimos haber usado la primera proposición de esta entrada para concluir estos hechos.

Nos falta estudiar el término de grado n1. Si j{1,2,,n}, entonces kjδjσ(j) es distinto de cero solo si σ(k)=k para todo kj: pero σ es una permutación, en particular una biyección, lo que fuerza que σ(j)=j también y entonces σ sea la identidad. Entonces el término de Xn1 en (Xδ1σ(1)a1σ(1))(Xδnσ(n)anσ(n)) es distinto de cero sólo cuando σ es la identidad. En ese caso es precisamente j=1najj=Tr(A).

◻

Ejemplo. Si A es la matriz del primer problema de esta entrada, tenemos que

χA(X)=(X1)(X+2)(X3)=X32X2++6.

Nota cómo el término de X2 es en efecto Tr(A)=(12+3) y el último es det(A).

Matrices nilpotentes

El caso de las matrices nilpotentes es todavía más sencillo.

Problema. Sea AMn(F) una matriz nilpotente. Es decir, existe k1 tal que Ak=On.

  1. Demuestra que
    χA(X)=Xn.
  2. Demuestra que TrAm=0 para todo m1.

Solución.

  1. Sea k1 tal que Ak=On (existe pues A es nilpotente). Entonces
    XkIn=XkInAk=(XInA)(Xk1In+Xk2A++Ak1).
    Tomando el determinante de ambos lados y recordando que abre productos llegamos a
    Xnk=det(XkIn)=χA(X)det(Xk1In++Ak1).
    De aquí, concluimos que χA(X) tiene que dividir a Xnk, pero sabemos que χA(X) es mónico y de grado n. Concluimos entonces que χA(X)=Xn.
  2. Puesto que Am también es una matriz nilpotente, el inciso anterior nos dice que
    χAm(X)=Xn.
    Pero sabemos por la sección sobre la traza que el término de Xn1 es Tr(Am). Como este término no aparece, concluimos que la traza es cero.

◻

Ejemplo. Para calcular el polinomio característico de la matriz

A=(53215961064)

podríamos notar (aunque no sea obvio a simple vista) que A2=O3. Luego, por el problema anterior, χA(X)=X3.

Un último caso particular

Acabamos con una última familia de matrices con polinomio característico simple. Esta familia está descrita por su forma, y será de particular importancia para el teorema de Cayley-Hamilton.

Problema. Para escalares a0,,an1F consideramos la matriz

A=(0000a01000a10100a20001an1).

en Mn(F).

Demuestra que

χA(X)=Xnan1Xn1a0.

Solución. Sea P(X)=Xnan1Xn1a0. Considera la matriz

B=XInA=(X000a01X00a101X0a20001Xan1).

Sumando el segundo renglón multiplicado por X al primer renglón, luego sumándole también al primer renglón el tercero multiplicado por X2, el cuarto por X3, y así sucesivamente hasta sumar el último renglón multiplicado por Xn1 llegamos a la matriz

C=(0000P(X)1X00a101X0a20001Xan1).

Recordamos que el determinante es invariante bajo sumas de renglones, por lo que

χA=detB=detC.

Expandiendo el determinante de C en el primer renglón obtenemos sencillamente

detC=(1)n+1P(X)|1X0010001|=(1)n+1P(X)(1)n1=P(X).

Para la segundaigualdad usamos que el determinante es el de una matriz triangular superior con puros 1 como entradas. Para la última, usamos que n+1+n1=2n siempre es un número par, así que queda 1 elevado a un número par. Esto concluye la prueba.

◻

Una de las consecuencias de la proposición anterior es que para cualquier polinomio mónico P de grado n en F[X], existe una matriz en Mn(F) tal que su polinomio característico es P.

Más adelante…

En la próxima entrada veremos unos últimos aspectos teóricos del polinomio característico antes de lanzarnos de lleno al teorema de Cayley-Hamilton y su demostración.

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. Encuentra una matriz A tal que χA(X)=X55X3+X22X+2. Sugerencia: Usa el último problema.
  2. Demuestra que el polinomio característico de una matriz A=[aij] triangular inferior está dado por i=1n(Xaii).
  3. Demuestra que 0 es eigenvalor de una matriz si y sólo si su determinante es cero.
  4. Calcula el polinomio característico de la siguiente matriz con entradas reales:
    A=(555666111111). Sugerencia: ¿Quién es A2?
  5. ¿Es cierto que si F es cualquier campo y A es una matriz con entradas en F, entonces el hecho de que Tr(A)=0 implica que A sea nilpotente? Sugerencia: Piensa en F2.
  6. Da una demostración alternativa al último problema de esta entrada usando inducción matemática sobre el tamaño de la matriz.

Entradas relacionadas

Agradecimientos

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

Álgebra Lineal II: Eigenvectores y eigenvalores

Por Julio Sampietro

Introducción

En esta entrada revisitamos los conceptos de eigenvalores y eigenvectores de una transformación lineal. Estos son esenciales para entender a las transformaciones lineales, y tienen un rango de aplicabilidad impresionante: aparecen en la física, las ecuaciones diferenciales parciales, la ciencia de datos, la topología algebraica y la probabilidad.

Primero enunciaremos la definición, después veremos un primer ejemplo para convencernos de que no son objetos imposibles de calcular. Luego daremos un método para vislumbrar una manera más sencilla de hacer dicho cálculo y concluiremos con unos ejercicios.

Eigen-definiciones

Comenzamos con V un espacio vectorial sobre F y T:VV una transformación lineal.

Definición. Un eigenvalor (también conocido como valor propio) de T es un escalar λF tal que λIdT no es invertible. Un eigenvector (también conocido como vector propio o λ-eigenvector) correspondiente a λ es un vector no-cero de ker(λIdT). A este kernel se le conoce como el eigenespacio correspondiente a λ (o λ-eigenespacio).

Entonces un λ-eigenvector es por definición distinto de cero y satisface

T(v)=λv.

Hay que tener cuidado. se permite que λ=0 sea eigenvalor, pero no se permite que v=0 sea eigenvector.

La colección de todos los eigenvectores, junto con el vector cero, es el eigenespacio asociado a λ. Podemos enunciar definiciones análogas con matrices.

Definición. Sea AMn(F) una matriz cuadrada. Un escalar λF es un eigenvalor de A si existe un vector XFn distinto de cero (un eigenvector) tal que AX=λX. En este caso el subespacio

ker(λInA):={XFnAX=λX}

es el λ-eigenespacio de A.

Puedes verificar que ambas definiciones se corresponden en el siguiente sentido:

Si V es un espacio de dimensión finita y T:VV es una transformación lineal, podemos escoger cualquier base de V y asociarle a T su forma matricial, digamos A, en esta base. Los eigenvalores de T son precisamente los eigenvalores de A. ¡Pero cuidado! Los eigenvectores de A dependerán de la base elegida.

Un primer ejemplo

Seguimos con un sencillo pero importante ejemplo.

Ejemplo 1. Considera la matriz

A=(0110).

Busquemos los eigenvectores y eigenvalores de A, pensando a A como una matriz con entradas complejas. Sea λC un eigenvalor y X un eigenvector asociado. Entonces se cumple la relación AX=λX. Si X=(x1,x2) entonces la condición mencionada es equivalente al par de ecuaciones

x2=λx1,x1=λx2.

Sustituyendo una en la otra obtenemos

x2=λ2x2.

Si x2=0 entonces x1=0 y así X es un vector nulo, lo que es imposible por definición (recuerda que pedimos que los eigenvectores sean distintos de cero). Entonces x20 y podemos dividir por x2 a la ecuación previa, de manera que λ2=1, o sea λ=±i. Conversamente, i y i son eigenvalores. En efecto, podemos tomar x2=1 y x1=λ como soluciones del problema anterior y obtener un vector propio asociado. De hecho, el eigenespacio está dado por

ker(λI2A)={(λx2,x2)x2C}

y esto no es más que la recta generada por el vector v=(λ,1)C2. Por lo tanto, vista como una matriz compleja, A tiene dos eigenvalores distintos ±i y dos eigenespacios, los generados por (i,1) y (i,1).

Por otro lado, veamos qué pasa si pensamos a A como una matriz con entradas reales. Haciendo las mismas cuentas llegamos a la misma ecuación, x2=λ2x2. Podemos reescribirla factorizando el término x2:

(λ2+1)x2=0.

Como λ esta vez es un número real, λ2+1 siempre es distinto de cero. Entonces para que el producto sea cero, tiene que ocurrir que x2=0, ¡pero entonces x1=0 y así X=0! En conclusión: vista como una matriz con entradas reales, A no tiene eigenvalores, y por tanto no tiene eigenespacios. La moraleja es que los eigenvalores y eigenvectores dependen mucho del campo en el que trabajemos.

¿Cómo calcularlos?

Si bien el ejemplo anterior resultó simple, no es difícil imaginar que matrices más complicadas y más grandes pueden resultar en procedimientos menos claros. En general:

  • ¿Cómo podemos calcular los eigenvalores?
  • ¿Cómo podemos calcular los eigenespacios de manera eficiente?
  • ¿Cómo podemos calcular los eigenvectores?

Una vez calculados los eigenvalores, calcular los eigenespacios se reduce a resolver el sistema de ecuaciones homogéneo (AλIn)X=0, lo cual ya hemos hecho muchas veces mediante reducción gaussiana. Luego, calcular los eigenvectores simplemente es tomar los elementos no cero del eigenespacio. Sin embargo, el cálculo de eigenvalores involucra encontrar raíces de polinomios lo cual de entrada no es obvio. Un primer paso es la siguiente observación que enunciamos como proposición.

Proposición. Un escalar λF es un eigenvalor de AMn(F) si y sólo si

det(λInA)=0.

Demostración. El sistema (λInA)X=0 tiene soluciones no triviales si y sólo si la matriz λInA no es invertible. A su vez, la matriz λInA no es invertible si y sólo si su determinante es nulo. El resultado se sigue.

◻

Regresemos a nuestra pregunta. Si

A=(a11a12a1na21a22a2nan1an2ann)

entonces la proposición nos dice que podemos calcular los valores propios de A resolviendo la ecuación polinomial

|λa11a12a1na21λa22a2nan1an2λann|=0

en F. Esta es una ecuación polinomial de grado n, y si el grado es mayor a 4 en general no existe una fórmula para resolverla en términos de radicales (aunque claro que hay casos particulares que si podemos resolver sin mucho problema).

Problema 2. Queremos calcular los eigenvalores de A, donde A está dada por

A=(100001010).

Solución. Como vimos en la proposición, esto se reduce a calcular las raíces del polinomio

|λ1000λ101λ|=0.

Calculando el determinante vemos que esto es de hecho

(λ1)(λ2+1)=0.

Sin embargo tenemos que recordar que las raíces dependen de nuestro campo de elección. Como no comentamos nada sobre el campo en el cual trabajamos, consideraremos dos casos. Si el campo es C entonces los eigenvalores son 1 y ±i. Si trabajamos sobre R entonces tenemos un único eigenvalor: 1.

Ejercicios

Acabamos esta entrada con unos ejercicios para reforzar lo que vimos.

Problema 1. Encuentra todos los números reales x tales que la matriz

A=(1x21)

tiene exactamente dos eigenvalores distintos. La misma pregunta para ningún eigenvalor.

Solución. El número de eigenvalores va a estar dado por el número de raíces del polinomio det(λI2A). Es decir, tenemos que trabajar la ecuación

det(λI2A)=|λ1x2λ1|=0.

Que a su vez se reduce a

(λ1)22x=0.

Y para que tenga dos soluciones basta con que 2x sea un número positivo. En efecto, en ese caso podemos despejar y resolver

λ=1±2x.

Como 2x es positivo solo si x lo es, podemos concluir que la condición necesaria y suficiente es que x sea un real positivo. Similarmente, si x es un número negativo no tendremos ningún eigenvalor.

Problema 2. Sea V el conjunto de todas las matrices AM2(C) tales que v=(12) es un eigenvector de A. Demuestra que V es un subespacio de M2(C) y da una base.

Solución. Supongamos que v es un eigenvector de A, con eigenvalor λ, y que es eigenvector de B, con eigenvalor μ. Entonces

(A+cB)(v)=Av+cBv=λv+cμv=(λ+cμ)v

por lo que v es eigenvector de A+cB con eigenvalor λ+cμ. Esto demuestra que V es un subespacio. Para darnos una idea de cómo podría ser una base para V, comencemos con una matriz genérica A=(abcd) tal que AV. Entonces A tiene que satisfacer Av=λv para algún λ. Escribamos esto más explícitamente

(abcd)(12)=(a+2bc+2d)=(λ2λ).

Así se desprenden dos ecuaciones

{a+2b=λc+2d=2λ.

Sabemos que λ es un parámetro libre, pues puede ser cualquier eigenvalor. Si conocemos a λ entonces necesitamos alguna de las variables, a o b para determinar a la otra y lo mismo con c y d. Entonces escojamos b y d como variables libres. Enseguida nuestra matriz es de la forma (reemplazando a a y c por sus valores en b y d):

A=(λ2bb2λ2dd)=b(2100)+d(0021)+λ(1020).

Entonces proponemos como base

β={(2100),(0021),(1020)}.

Ya vimos que β genera a V, y dejamos la independencia lineal como ejercicio.

◻

Más adelante…

En las próximas entradas desarrollaremos las propiedades relevantes de los eigenvalores y eigenvectores para eventualmente llegar al polinomio característico y establecer el puente con el polinomio mínimo.

Tarea moral

Aquí unos ejercicios para que repases el material de esta entrada.

  1. Encuentra todos los eigenvalores de la matriz A=(110021001)M3(C).
  2. Completa la demostración del último ejercicio de la sección de ejercicios, verificando que las soluciones encontradas son matrices linealmente independientes. ¿Puedes generalizar este ejercicio de alguna manera?
  3. Encuentra los eigenvalores de la matriz AMn(R) cuyas entradas son puros 2.
  4. Da contraejemplos para cada una de las siguientes afirmaciones:
    1. Si u y v son eigenvectores de A, entonces u+v es eigenvector de A.
    2. Si λ es eigenvalor de A y μ es eigenvalor de B, entonces λμ es eigenvalor de AB.
    3. Si A y B son formas matriciales de una misma transformación T y v es eigenvector de A, entonces v es eigenvector de B.
  5. Considera la transformación derivada en R[x]. ¿Quienes son sus eigenvectores y eigenvalores? Como sugerencia, estudia el coeficiente de mayor grado.

Entradas relacionadas

Agradecimientos

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

Álgebra Lineal II: Polinomio mínimo de transformaciones lineales y matrices

Por Julio Sampietro

Introducción

Anteriormente definimos qué quiere decir evaluar un polinomio en una matriz o en una transformación lineal. En esta entrada definiremos uno de los objetos más importantes del álgebra lineal: el polinomio mínimo. Si bien al principio nos va a costar un poco calcularlo, esto se compensa por la cantidad de propiedades teóricas que cumple. Comenzaremos dando su definición, y mostrando su existencia y unicidad. Luego exploraremos algunas propiedades y veremos ejemplos, seguido de un pequeño teorema de cambio de campos. Finalmente introduciremos un objeto similar (el polinomio mínimo puntual) y haremos unos ejercicios para cerrar.

El concepto de polinomio mínimo podría resultarle familiar a los más algebraicos de mente: ¡todo se debe a que trabajamos con dominios de ideales principales, o incluso euclidianos! Si has trabajado anteriormente con conceptos como el mínimo común múltiplo en enteros, puede que varios de los argumentos de esta entrada te suenen conocidos.

Existencia y unicidad

Comenzamos con un espacio vectorial V de dimensión n sobre un campo F. Fijando una transformación lineal T:VV, queremos entender para qué polinomios se cumple que P(T)=0. Nota como podríamos haber cambiado la pregunta: si fijamos un polinomio P, podríamos buscar todas las transformaciones T tales que P(T)=0. Ésta pregunta la estudiaremos más adelante.

Definimos el conjunto

I(T)={PF[X]P(T)=0}.

El polinomio cero pertenece a I(T) de manera trivial. Una cosa importante es que este conjunto I(T) que vamos a estudiar en verdad es «interesante», en el sentido de que debemos ver que hay más polinomios adentro y no es únicamente el conjunto {0}. Una manera de ver esto es sabiendo que el espacio de transformaciones lineales de V en V tiene dimensión n2 (lo puedes pensar como el espacio de matrices). Entonces, las n2+1 transformaciones Id,T,T2,,Tn2 no pueden ser todas linealmente independientes: uno de los corolarios del lema de Steinitz es que en un espacio de dimensión n a lo más se pueden tener n vectores linealmente independientes. Entonces existe una combinación lineal no trivial y nula

a0Id+a1T++an2Tn2=0.

Luego a0+a1X++an2Xn2 es un polinomio no cero tal que P(T)=0, es decir PI(T).

Con el argumento de arriba vimos que I(T) es «interesante» en el sentido de que tiene polinomios no cero. El siguiente teorema se puede entender como que I(T) se puede describir muy fácilmente.

Teorema. Existe un único polinomio mónico, distinto de cero μT tal que I(T) es precisamente el conjunto de múltiplos de μT. Es decir

I(T)=μTF[X]={μTP(X)P(X)F[X]}.

La demostración hará uso del algoritmo de la división para polinomios. Te lo compartimos aquí, sin demostración, por si no lo conoces o no lo recuerdas.

Teorema (algoritmo de la división en F[x]). Sean M(x) y N(x) polinomios en F[x], donde N(x) no es el polinomio cero. Entonces, existen únicos polinomios Q(x) y R(x) en F[x] tales que M(x)=Q(x)N(x)+R(x), en donde R(x) es el polinomio cero, o deg(R(x))<deg(G(x)).

Si te interesa saber cómo se demuestra, puedes seguir la teoría de polinomios disponible en la Unidad 4 del curso de Álgebra Superior II.

Demostración. Veamos primero que I(T) es un subespacio de F[X]. Para ello, tomemos polinomios P(x), Q(x) en I(T), y un escalar αF. Una de las proposiciones de la entrada pasada nos permite abrir la expresión (P+αQ)(T) como P(T)+αQ(T)=0+α0=0, de modo que P+αQ está en I(T) y por lo tanto I(T) es un subespacio de F[X].

Por otro lado si PI(T) y QF[X] entonces

(PQ)(T)=P(T)Q(T)=0Q(T)=0.

Lo que discutimos antes de enunciar el teorema nos dice que I(T){0}. Tomemos entonces PI(T) un polinomio no cero de grado mínimo. Podemos suponer sin perdida de generalidad que P es mónico, de no serlo, podemos dividir a P por su coeficiente principal sin cambiar el grado.

La ecuación previa nos indica que todos los múltiplos polinomiales de P también están en I(T). Veamos que todo elemento de I(T) es de hecho un múltiplo de P. Si SI(T), usamos el algoritmo de la división polinomial para escribir S=QP+R con Q,RF[X]. Aquí hay dos casos: que R sea el polinomio cero, o bien que no lo sea y entonces degR<degP. Nota que R=SQPI(T) dado que I(T) es un subespacio de F[X] y S,QPI(T). Si R0, entonces como degR<degP llegamos a una contradicción de la minimalidad del grado de P. Luego R=0 y por tanto S=QP. Entonces I(T) es precisamente el conjunto de todos los múltiplos de P y así podemos tomar μT=P.

Para verificar la unicidad de μT, si otro polinomio S tuviera las mismas propiedades, entonces S dividiría a μT y μT dividiría a S. Sin embargo, como ambos son mónicos se sigue que deben ser iguales: en efecto, si μT=SQ y S=μTR entonces degQ=degR=0, porlo tanto son constantes, y como el coeficiente principal de ambos es 1, se sigue que ambos son la constante 1 y así μT=S. Esto completa la demostración.

◻

Definición. Al polinomio μT se le conoce como el polinomio mínimo de T.

Primeras propiedades y ejemplos

Debido a su importancia, recalcamos las propiedades esenciales del polinomio mínimo μT:

  • Es mónico.
  • Cumple μT(T)=0.
  • Para cualquier otro polinomio PF[X], sucede que P(T)=0 si y sólo si μT divide a P.

Toda la teoría que hemos trabajado hasta ahora se traduce directamente a matrices usando exactamente los mismos argumentos. Lo enunciamos de todas maneras: si AMn(F) es una matriz cuadrada, entonces existe un único polinomio μAF[X] con las siguientes propiedades:

  • Es mónico.
  • Cumple μA(A)=On.
  • Si PF[X], entonces P(A)=On si y sólo si μA divide a P.

Como jerga, a veces diremos que un polinomio «anula T» si P(T)=0. En este sentido los polinomios que anulan a T son precisamente los múltiplos de μT.

Vimos antes de enunciar el teorema que podemos encontrar un polinomio P no cero de grado menor o igual a n2 tal que P(T)=0. Como μT divide a P se sigue que degμTn2. Esta cota resulta ser débil, y de hecho un objeto que hemos estudiado previamente nos ayudará a mejorarla: el polinomio característico. Este también va a anular a T y con ello obtendremos una mejor cota: degμTn.

Ejemplo 1. Si A=On, entonces μA=X. En efecto, μA(A)=0 y además es el polinomio de menor grado que cumple esto, pues ningún polinomio constante y no cero anula a On (¿por qué?). Nota como además I(A) es precisamente el conjunto de polinomios sin término constante.

Ejemplo 2. Considera la matriz AM2(R) dada por

A=(0110).

Nos proponemos calcular μA. Nota que A satisface A2=I2. Por tanto el polinomio P(X)=X2+1 cumple P(A)=0. Así, μA tiene que dividir a este polinomio ¡pero este es irreducible sobre los números reales! En efecto, si existiese un factor propio de P sobre R, tendríamos que la ecuación X2=1 tiene solución, y sabemos que este no es el caso. Entonces μA tiene que ser X2+1.

Ejemplo 3. Sean d1,,dnF escalares y A una matriz diagonal tal que [aii]=di. Los elementos pueden no ser distintos entre sí, así que escogemos una colección máxima di1,,dik de elementos distintos. Para cualquier polinomio P, tenemos que P(A) es simplemente la matriz diagonal con entradas P(di) (esto porque el producto An tiene como entradas a din). Entonces para que P(A)=0 se tiene que cumplir que P(di)=0, y para que esto pase es suficiente que P(dik)=0. Eso quiere decir que P tiene al menos a los dik como raíces, y entonces (Xdi1)(Xdi2)(Xdik) divide a P.

Nota como esto es suficiente: encontramos un polinomio mónico, (Xdi1)(Xdi2)(Xdik) que divide a cualquier P tal que P(A)=0. Así

μA(X)=(Xdi1)(Xdik).

Cambio de campos

En uno de los ejemplos argumentamos que el polinomio mínimo era X2+1 porque este es irreducible sobre R. Pero, ¿qué pasaría si cambiáramos nuestro campo a C? La situación puede ser incluso más delicada: a una matriz con entradas racionales la podemos considerar como una instancia particular de una matriz con entradas reales, que a su vez podemos considerar como una matriz compleja. ¿Hay tres polinomios mínimos distintos? El siguiente teorema nos da una respuesta tranquilizante.

Teorema. Sean F1F2 dos campos y AMn(F1) una matriz, entonces el polinomio mínimo de A vista como elemento de Mn(F1) y el polinomio mínimo de A vista como elemento de Mn(F2) son iguales.

Demostración. Sea μ1 el polinomio de AMn(F1) y μ2 el polinomio mínimo de AMn(F2). Puesto que F1[X]F2[X], se tiene que μ1F2[X] y además μ1(A)=0 por definición. Luego μ2 necesariamente divide a μ1. Sean d1=degμ1 y d2=degμ2, basta verificar que d2d1 y para que esto se cumpla basta con encontrar PF1[X] de grado a lo más d2 tal que P(A)=0 (entonces μ1 dividiría a este polinomio y se sigue la desigualdad).

Desarrollando que μ2(A)=0 en todas sus letras (o mejor dicho, en todos sus coeficientes) se tiene

a0In+a1A++ad2Ad2=On.

Esto es equivalente a tener n2 ecuaciones homogéneas en las variables a0,,ad2. Como A tiene entradas en F1 los coeficientes de estas ecuaciones todos pertenecen a F1. Tenemos un sistema de ecuaciones con coeficientes en F1 que tiene una solución no trivial en F2: tiene automáticamente una solución no trivial en F1 por un ejercicio de la entrada de Álgebra Lineal I de resolver sistemas de ecuaciones usando determinantes. Esto nos da el polinomio buscado.

◻

Mínimos puntuales

Ahora hablaremos (principalmente a través de problemas resueltos) de otro objeto muy parecido al polinomio mínimo: el polinomio mínimo puntual. Este es, esencialmente un «polinomio mínimo en un punto». Más específicamente si T:VV es lineal con polinomio mínimo μT y xV definimos

Ix={PF[X]P(T)(x)=0}.

Nota que la suma y diferencia de dos elementos en Ix también está en Ix.

Problema 1. Demuestra que existe un único polinomio mónico μxF[X] tal que Ix es el conjunto de múltiplos de μx en F[X]. Más aún, demuestra que μx divide a μT.

Solución. El caso x=0 se queda como ejercicio. Asumamos entonces que x0. Nota que μTIx puesto que μT(T)=0. Sea μx el polinomio mónico de menor grado en Ix. Demostraremos que Ix=μxF[X].

Primero si PμxF[X] entonces por definición P=μxQ para algún QF[X] y entonces

P(T)(x)=Q(T)(μx(T)(x))=Q(T)(0)=0.

Así PIx, y queda demostrado que μxF[X]Ix.

Conversamente, si PIx podemos usar el algoritmo de la división para llegar a una expresión de la forma P=Qμx+R para algunos polinomios Q,R con degR<degμx. Supongamos que R0. Similarmente a como procedimos antes, se cumple que R=PQμxIx dado que Ix es cerrado bajo sumas y diferencias. Dividiendo por el coeficiente principal de R, podemos asumir que R es mónico. Entonces R es un polinomio mónico de grado estrictamente menor que el grado de μx, una contradicción a nuestra suposición: μx es el polinomio de grado menor con esta propiedad. Luego R=0 y μx divide a P.

Así queda probado que si PIx entonces PμxF[X], lo que concluye la primera parte del problema. Para la segunda, vimos que μTIx y por tanto μx divide a μT.

◻

Problema 2. Sea Vx el subespacio generado por x,T(x),T2(x),. Demuestra que Vx es un subespacio de V de dimensión degμx, estable bajo T.

Solución. Es claro que Vx es un subespacio de V. Además, dado que T manda a generadores en generadores, también es estable bajo T. Sea d=degμx. Demostraremos que x,T(x),,Td1(x) forman una base de Vx, lo que concluiría el ejercicio.

Veamos que son linealmente independientes. Si a0x+a1T(x)+a2T2(x)++ad1Td1(x)=0 para algunos escalares ai no todos cero, entonces el polinomio

P=a0+a1X++ad1Xd1

es un elemento de Ix, pues P(T)(x)=0. Luego μx necesariamente divide a P, pero esto es imposible puesto que el grado de P es d1, estrictamente menor que el grado de μx. Luego los ai deben ser todos nulos, lo que muestra que x,T(x),T2(x),,Td1(x) es una colección linealmente independiente.

Sea W el espacio generado por x,T(x),,Td1(x). Afirmamos que W es invariante bajo T. Es claro que T(x)W, similarmente T(T(x))=T2(x)W y así sucesivamente. El único elemento «sospechoso» es Td1(x), para el cual basta verificar que T(Td1(x))=Td(x)W. Dado que μx(T)(x)=0 y μx es mónico de grado d, existen escalares bi (más precisamente, los coeficientes de μx) no todos cero tales que

Td(x)+bd1Td1(x)++b0x=0.

Esto nos muestra que podemos expresar a Td(x) en términos de x,T(x),,Td1(x) y por tanto Td(x) pertenece a W.

Ahora, dado que W es estable bajo T y contiene a x, se cumple que Tk(x)W para todo k0. En particular VxW. Luego Vx=W (la otra contención es clara) y x,T(x),,Td1(x) genera a W, o sea a Vx.

Mostramos entonces que x,T(x),,Td1(x) es una base para Vx y así dimVx=d.

◻

Unos ejercicios para terminar

Presentamos unos últimos ejercicios para calcular polinomios mínimos.

Problema 1. Calcula el polinomio mínimo de A donde

A=(010100001).

Solución. A estas alturas no tenemos muchas herramientas que usar. Comenzamos con calcular A2:

A2=(010100001)(010100001)=(100010001).

Entonces en particular A2=I3. Así, el polinomio mínimo μA tiene que dividir a X21. Este último se factoriza como (X1)(X+1), pero es claro que A no satisface ni AI3=0 ni A+I3=0. Entonces μA no puede dividir propiamente a X21, y por tanto tienen que ser iguales.

Problema 2. Calcula el polinomio mínimo de la matriz A con

A=(1201).

Solución. Nota como

AI2=(0200)

y es fácil verificar que el cuadrado de la matriz de la derecha es cero. Así (AI2)2=0, o sea, el polinomio P(X)=(X1)2 anula a A. Similarmente al problema anterior, μA tiene que dividir a P, pero P sólo tiene un factor: X1. Dado que A no satisface AI2=0 se tiene que μA no puede dividir propiamente a P, y entonces tienen que ser iguales. Luego μA=(X1)2=X22X+1.

Más adelante…

En las entradas subsecuentes repasaremos los eigenvalores y eigenvectores de una matriz, y (como mencionamos) ligaremos el polinomio característico de una matriz con su polinomio mínimo para entender mejor a ambos.

Tarea moral

Aquí unos ejercicios para practicar lo que vimos.

  1. Encuentra una matriz A cuyo polinomio mínimo sea X2. Para cada n, ¿puedes encontrar una matriz cuyo polinomio mínimo sea Xn?
  2. Encuentra una matriz A cuyo polinomio mínimo sea X21. Para cada n, ¿puedes encontrar una matriz cuyo polinomio mínimo sea Xn1?
  3. Encuentra el polinomio de la matriz A en Mn(F) cuyas entradas son todas 1.
  4. Si T:Mn(R)Mn(R) es la transformación que manda a cada matriz en su transpuesta, encuentra el polinomio mínimo de T.
  5. Sea V un espacio vectorial y x,y vectores linealmente independientes. Sea T:VV una transformación lineal. ¿Cómo son los polinomios P tales que P(T) se anula en todo el subespacio generado por x y y? ¿Cómo se relacionan con los polinomios mínimos puntuales de T para x y y?

Entradas relacionadas

Agradecimientos

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

Seminario de Resolución de Problemas: El teorema fundamental del cálculo

Por Fabian Ferrari

Introducción

Ya platicamos de continuidad, diferenciabilidad e integrales, así como de otros temas de cálculo. En esta sección reuniremos varias de estas ideas a través de uno de los resultados más importantes: el teorema fundamental del cálculo. Este teorema nos exhibe la relación que hay entre la derivada y la integral, distinguiéndolas como procedimientos inversos el uno del otro.

El teorema nos dice que si tenemos una función F(x) derivable sobre un intervalo [a,b], entonces

abF(t)dt=F(b)F(a).

Ahora bien, si nuestra función F(t) es derivable en [0,x], tenemos que

0xF(t)dt=F(x)F(0),

a lo que le sigue que

F(x)=0xF(t)dt+F(0).

Esto nos recuerda a la constante de integración

F(x)=0xF(t)dt+C.

Es decir, tenemos que C=F(0).

Aquí en el blog, en la entrada «Teoremas fundamentales de los cuadraditos» damos la intuición acerca de este teorema, comenzando con el caso discreto. Puedes leerlo antes de continuar.

Usar el teorema fundamental del cálculo para obtener una identidad trigonométrica

Veamos un ejemplo. Tenemos que la derivada de la función F(t)=sin2t es F(t)=2costsint. Por el teorema fundamental del cálculo, la integral de F(t) en el intervalo [0,x] está dada por

0x2sintcostdt=sin2x,

en donde usamos que F(0)=sin2(0)=0.

Por otro lado, resolviendo la integral utilizando el cambio de variable u=cost, tenemos que 0x2sintcostdt=cos2t|0x=cos2x+1.

Igualando ambos valores de la integral, tenemos que sin2x=cos2x+1. De aquí obtenemos la identidad trigonométrica pitagórica sin2x+cos2x=1 para toda x.

Veamos ahora un problema en el que, mediante el problema fundamental del cálculo,

Problema. Aplicando el teorema fundamental del calculo halla absecxdx.

Sugerencia pre-solución. Formula un problema equivalente multiplicando y dividiendo la expresión por secx+tanx. Intenta identificar la expresión resultante como la derivada de otra función.

Solución. Para resolver este problema tenemos que hallar una función F(x) de tal forma que F(x)=secx.

Para ello, tenemos que notar que

secx=secx(secx+tanxsecx+tanx)=sec2x+secxtanxsecx+tanx.

Y entonces la derivada de ln(secx+tanx) es igual a

(1secx+tanx)(sec2x+secxtanx)=sec2x+secxtanxsecx+tanx=secx.

Proponemos a la función

F(x)=ln(secx+tanx)

dado que

F(x)=secx.

Ahora, aplicando el teorema fundamental del cálculo tenemos que

absecxdx=F(b)F(a)=ln(secb+tanb)ln(seca+tana)

◻

Segundo teorema fundamental del cálculo

Veamos una implicación del teorema fundamental del cálculo, que también se le conoce como el «segundo teorema fundamental del cálculo».

Para una función f:[a,b]R continua en el intervalo [a,b] se tiene que:

ddx(axf(t)dt)=f(x)

Problema. Determina ddx(3x101t+4dt).

Sugerencia pre-solución. Usa el segundo teorema fundamental del cálculo y la regla de la cadena.

Solución. Como 3x101t+4dt=03x11t+4dt, tenemos entonces que

ddx(3x101t+4dt)=ddx(03x11t+4dt).

Por otro lado, consideremos las funciones

f(x)=0x1t+4dtyg(x)=3x1.

Aplicando el teorema fundamental del cálculo y derivando tenemos que

f(x)=1x+4yg(x)=3.

Notemos que

(fg)(x)=f(g(x))=f(3x1)=03x11t+4dt.

Así, aplicando la regla de la cadena, tenemos que

ddx(03x11t+4dt)=ddx(f(g(x))=f(g(x))g(x)=1(3x1)+43=1x+1.

◻

Veamos un último problema en el que se usa la segunda forma del teorema fundamental del cálculo.

Problema: Supongamos que f es una función continua para toda x, la cual satisface la ecuación

(1)0xf(t)dt=x1t2f(t)dt+x168+x189+C,

donde C es una constante. Encuentra la forma explícita de la función f(x) y determina el valor de la constante C.

Sugerencia pre-solución.

Solución. De la ecuación, tenemos lo siguiente

ddx(0xf(t)dt)=ddx(x1t2f(t)dt+x168+x189+C)

Como f es continua para toda x, por el teorema fundamental del cálculo en su segunda forma tenemos que

ddx(0xf(t)dt)=f(x)

y

ddx(x1t2f(t)dt)=ddx(1xt2f(t)dt)=x2f(x).

Entonces, derivando ambos lados de la expresión original nos resulta la ecuación

f(x)=x2f(x)+2x15+2x17,

de la cual se obtiene

f(x)(x2+1)=2x15+2x17=2x15(x2+1)

Así, tenemos que

f(x)=2x15.

Sustituyendo f(t)=2t15 en la ecuación (1), tenemos que

0x2t15dt=x1t2(2t15)dt+x168+x189+C

Así,

0x2t15dt=x1t2(2t15)dt+x168+x189+C0x2t15dt=1x2t17dt+x168+x189+C2t1616|0x=(2t1818)|1x+x168+x189+Cx168=(x18919)+x168+x1819+C

Con ello, tenemos que

C+19=0.

Por lo tanto la función que satisface la ecuación es f(x)=2x15 y el valor de la constante es C=19.

◻

Más problemas

Hay más ejemplos de problemas relacionados con la aplicación del teorema fundamental del cálculo en la Sección 6.9 del libro Problem Solving through Problems de Loren Larson.

Álgebra Lineal I: Formas bilineales, propiedades, ejemplos y aclaraciones

Por Leonardo Ignacio Martínez Sandoval

Introducción

En entradas anteriores hemos platicado de dualidad, ortogonalidad y transformaciones transpuestas. Es importante que repases esas entradas y nos escribas si tienes dudas, pues ahora pasaremos a un tema un poco diferente: formas bilineales y cuadráticas. Estas nociones nos permitirán seguir hablando acerca de la geometría de espacios vectoriales en general.

Para esta parte del curso, nos vamos a enfocar únicamente en espacios vectoriales sobre R. Se pueden definir los conceptos que veremos para espacios vectoriales en otros campos. Sobre todo, es posible definir conceptos análogos en C y obtener una teoría muy rica. Pero por ahora consideraremos sólo el caso de espacios vectoriales reales.

Aunque hablaremos de formas bilineales en general, una subfamilia muy importante de ellas son los productos interiores, que nos permiten hablar de espacios euclideanos. El producto interior es el paso inicial en una cadena muy profunda de ideas matemáticas:

  • Un producto interior nos permite definir la norma de un vector.
  • Con la noción de norma, podemos definir la distancia entre dos vectores.
  • A partir de un producto interior y su norma podemos mostrar la desigualdad de Cauchy-Schwarz, con la cual podemos definir ángulos entre vectores (por ejemplo, ¡podremos definir el ángulo entre dos polinomios!).
  • De la desigualdad de Cauchy-Schwarz, podemos probar que la noción de norma satisface la desigualdad del triángulo, y que por lo tanto la noción de distancia define una métrica.
  • Aunque no lo veremos en este curso, más adelante verás que una métrica induce una topología, y que con una topología se puede hablar de continuidad.

En resumen, a partir de un producto interior podemos hacer cálculo en espacios vectoriales en general.

Una forma bilineal con la cual probablemente estés familiarizado es el producto punto en Rn, que a dos vectores (x1,x2,,xn) y (y1,y2,,yn) los manda al real x1y1+x2y2++xnyn. Este es un ejemplo de una forma bilineal que es un producto interior. También puede que estés familiarizado con la norma en Rn, que a un vector (x1,,xn) lo manda al real x12+x22++xn2. Lo que está dentro de la raíz es un ejemplo de una forma cuadrática positiva definida. Incluyendo la raíz, este es un ejemplo de norma en espacios vectoriales.

Hay muchas otras formas bilineales y formas cuadráticas, pero los ejemplos mencionados arriba te pueden ayudar a entender la intuición detrás de algunos de los conceptos que mencionaremos. Para marcar algunas cosas en las que la intuición puede fallar, pondremos algunas «Aclaraciones» a lo largo de esta entrada.

En el futuro, tener una buena noción de la geometría de espacios vectoriales te ayudará a entender mucho mejor los argumentos de cursos de análisis matemático, de variable compleja y de optativas como geometría diferencial. Dentro de este curso, entender bien el concepto de forma bilineal te será de gran utilidad para cuando más adelante hablemos de formas multilineales y determinantes.

Formas bilineales

La definición fundamental para los temas que veremos en estas entradas es la siguiente, así que enunciaremos la definición, veremos varios ejemplos y haremos algunas aclaraciones.

Definición. Sea V un espacio vectorial sobre R. Una forma bilineal es una función b:V×VR tal que:

  • Para todo x en V, la función b(x,):VR que manda vV a b(x,v) es una forma lineal.
  • Para todo y en V, la función b(,y):VR que manda vV a b(v,y) es una forma lineal.

Ejemplo 1. Considera el espacio vectorial de polinomios R3[x] y considera la función b(p,q)=p(0)q(10)+p(1)q(11). Afirmamos que b es una forma bilineal. En efecto, fijemos un polinomio p y tomemos dos polinomios q1, q2 y un real r. Tenemos que
b(p,q1+rq2)=p(0)(q1+rq2)(10)+p(1)(q1+rq2)(11)=p(0)q1(10)+p(1)q1(11)+r(p(0)q2(10)+p(1)q2(11))=b(p,q1)+rb(p,q2),

De manera similar se puede probar que para q fijo y p1, p2 polinomios y r real tenemos que b(p1+rp2,q)=b(p1,q)+rb(p2,q). Esto muestra que b es una forma bilineal.

Si v=0, entonces por el primer inciso de la definición, b(x,v)=0 para toda x y por el segundo b(v,y)=0 para toda y, en otras palabras:

Proposición. Si b es una forma bilineal en b, y alguno de x o y es 0, entonces b(x,y)=0.

De la linealidad de ambas entradas de b, se tiene la siguiente proposición.

Proposición. Tomemos b:V×VR una forma bilineal, vectores x1,,xn, y1,,ym y escalares a1,,an,c1,,cm. Tenemos que b(i=1naixi,j=1mcjyj)=i=1nj=1maicjb(xi,yj).

La proposición anterior muestra, en particular, que para definir una forma bilineal en un espacio vectorial V de dimensión finita n, basta tomar una base {e1,,en} de V y definir b(ei,ej) para toda 1i,jn.

Hagamos algunas aclaraciones acerca de las formas bilineales.

Aclaración 1. No es lo mismo una forma bilineal en V, que una transformación lineal de V×V a R.

Ejemplo 2. La transformación b((w,x),(y,z))=w+x+y+z sí es una transformación lineal de R2×R2R, lo cual se puede verificar fácilmente a partir de la definición. Sin embargo, no es una forma bilineal. Una forma de verlo es notando que b((0,0),(1,1))=0+0+1+1=2. Aquí una de las entradas es el vector cero, pero el resultado no fue igual a cero.

Aclaración 2. Puede pasar que ninguna de las entradas de la forma bilineal sea 0, pero que evaluando en ella sí de 0.

Ejemplo 3. Consideremos la transformación b:R2×R2R tal que b((w,x),(y,z))=wyxz. Verificar que esta es una forma bilineal es sencillo y se deja como tarea moral. Además, se tiene que b((1,0),(0,1))=0.

Más adelante, cuando definamos producto interior, nos van a importar mucho las parejas de vectores v, w para las cuales b(v,w)=0.

Aclaración 3. Si b es una forma bilineal, no necesariamente es cierto que b(x,y)=b(y,x).

Ejemplo 4. Consideremos la transformación b:R2×R2R tal que b((w,x),(y,z))=wzxy. Verificar que esta es una forma bilineal es sencillo y se deja como tarea moral. Notemos que b((2,1),(2,3))=62=4, mientras que b((2,3),(2,1))=26=4.

Aquellas formas para las que sí sucede que b(x,y)=b(y,x) son importantes y merecen un nombre especial.

Definición. Una forma bilineal b:V×VR es simétrica si b(x,y)=b(y,x) para todo par de vectores x,y en V.

Para definir una forma bilineal b simétrica en un espacio V de dimensión finita n, basta tomar una base {e1,,en} y definir b en aquellas parejas b(ei,ej) con 1ijn.

Más ejemplos de formas bilineales

A continuación enunciamos más ejemplos de formas bilineales, sin demostración. Es un buen ejercicio verificar la definición para todas ellas.

Ejemplo 1. Si a1,a2,,an son números reales y V=Rn, entonces podemos definir b:V×VR que manda a x=(x1,,xn) y y=(y1,,yn) a b(x,y)=a1x1y1++anxnyn.

Este es un ejemplo de una forma bilineal simétrica. Si todos los ai son iguales a 1, obtenemos el producto punto o producto interior canónico de Rn.

Ejemplo 2. Tomemos V como el espacio vectorial de matrices Mn(R). La transformación b:V×VR tal que b(A,B)=tr(AB) es una forma bilineal. Además, es simétrica, pues la traza cumple la importante propiedad tr(AB)=tr(BA), cuya verificación queda como tarea moral.

Ejemplo 3. Tomemos V el conjunto de funciones continuas y de periodo 2π que van de R a sí mismo. Es decir, f:RR está en V si es continua y f(x)=f(x+2π) para todo real x. Se puede mostrar que V es un subespacio del espacio de funciones continuas, lo cual es sencillo y se queda como tarea moral. La transformación b:V×VR tal que b(f,g)=ππf(x)g(x)dx es una forma bilineal.

Ejemplo 4. Consideremos V=R[x], el espacio vectorial de polinomios con coeficientes reales. Para P y Q polinomios definimos b(P,Q)=n=1P(n)Q(2n)2n.

La serie de la derecha converge absolutamente, de modo que esta expresión está bien definida. Se tiene que b es una forma bilineal, pero no es simétrica.

Formas cuadráticas

Otra definición fundamental es la siguiente

Definición. Una forma cuadrática es una transformación q:VR que se obtiene tomando una forma bilineal b:V×VR y definiendo q(x)=b(x,x).

Aclaración 4. Es posible que la forma bilineal b que define a una forma cuadrática no sea única.

Ejemplo. Consideremos a la forma bilineal de R2 tal que b((x,y),(w,z))=xzyw. La forma cuadrática dada por b es q(x,y)=b((x,y),(x,y))=xyyx=0. Esta es la misma forma cuadrática que la dada por la forma bilineal b((x,y),(w,z))=ywxz. Pero b y b son formas bilineales distintas, pues b((1,0),(0,1))=1, mientras que b((1,0),(0,1))=1.

La aclaración anterior dice que puede que haya más de una forma bilineal que de una misma forma cuadrática. Sin embargo, resulta que la asignación es única si además pedimos a la forma bilineal ser simétrica. Este es el contenido del siguiente resultado importante.

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.

En la siguiente entrada mostraremos el teorema de la identidad de polarización. Por el momento, para tomar más intuición, observa como la identidad se parece mucho a la igualdad xy=(x+y)2x2y22 en números reales.

Más adelante…

En esta entrada estudiamos una extensión de la noción de transformaciones lineales que ya habíamos discutido en la unidad anterior. Enunciamos algunos teoremas muy importantes sobre las transformaciones bilineales e hicimos algunos ejemplos de cómo podemos verificar si una transformación es bilineal. La noción de transformación bilineal, nos permitirá abordar un concepto muy importante: el producto interior.

En las siguientes entradas hablaremos del producto interior y cómo éste nos ayuda a definir ángulos y distancias entre vectores de un espacio vectorial.

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.

  • Completa los detalles de la segunda parte del primer ejemplo.
  • Verifica que en efecto las transformaciones de los ejemplos de las aclaración 2 y 3 son formas bilineales.
  • Muestra que el subconjunto de funciones continuas R a R y de cualquier periodo p es un subespacio del espacio vectorial C(R) de funciones continuas reales.
  • Demuestra que para A y B matrices en Mn(F) se tiene que tr(AB)=tr(BA).
  • Encuentra una forma cuadrática en el espacio vectorial R3[x] que venga de más de una forma bilineal.
  • Muestra que el conjunto de formas bilineales de V es un subespacio del espacio de funciones V×VR. Muestra que el conjunto de formas bilineales simétricas de V es un subespacio del espacio de formas bilineales de V.
  • Piensa en cómo la igualdad xy=(x+y)2x2y22 de números reales está relacionada con la identidad de polarización para el producto punto en Rn.

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»