Mathematics · Libro 4 · Bachelor Year 2

Matemáticas universitarias — Grado 2

Matemáticas universitarias — Grado 2 · Bachelor Year 2

2Álgebra lineal

El álgebra lineal del volumen del primer año trabajaba sobre R\R o C\C en dimensión finita y admitía el determinante general. Este capítulo levanta las tres restricciones: la teoría se enuncia sobre un cuerpo KK arbitrario, se desarrolla sistemáticamente el juego entre un espacio y su dual (bases duales, anuladores, traspuestas) y el determinante queda por fin construido a partir de las formas multilineales alternadas y de la signatura del Capítulo 1, saldando así todas las deudas del primer año.

En todo el capítulo, KK es un cuerpo (Q\Q, R\R, C\C o Z/pZ\Z/p\Z; a la teoría le da igual) y, salvo mención expresa, los espacios son de dimensión finita sobre KK. Los resultados del primer año (bases, dimensión, teorema del rango, matrices) se trasladan literalmente: sus demostraciones nunca usaron nada más que los axiomas de cuerpo.

2.1 Espacio dual

Definición 2.1 (Espacio dual, base dual)

El dual de EE es E=L(E,K)E^* = \mathcal{L}(E, K), el espacio de las formas lineales. Si B=(e1,,en)\mathcal{B} = (e_1, \dots, e_n) es una base de EE, las formas coordenadas e1,,ene_1^*, \dots, e_n^* definidas por ei(ej)=δije_i^*(e_j) = \delta_{ij} (delta de Kronecker: 11 si i=ji = j, y 00 en otro caso) forman la base dual B\mathcal{B}^* de EE^*; en particular dimE=dimE\dim E^* = \dim E, y

x=i=1nei(x)ei(xE),φ=i=1nφ(ei)ei(φE).x = \sum_{i=1}^{n} e_i^*(x)\, e_i \quad (x \in E), \qquad \varphi = \sum_{i=1}^{n} \varphi(e_i)\, e_i^* \quad (\varphi \in E^*).

Demostración de que B\mathcal{B}^* es una base. Libre: aplicando una combinación nula λiei=0\sum \lambda_i e_i^* = 0 a eje_j se obtiene λj=0\lambda_j = 0. Generadora: para φE\varphi \in E^*, la forma φiφ(ei)ei\varphi - \sum_i \varphi(e_i) e_i^* anula todos los eje_j y por tanto es nula (una aplicación lineal que se anula sobre una base es nula). Las dos fórmulas destacadas son esos mismos cálculos leídos en sentido directo.

Ejemplo 2.2

En Kn[X]K_n[X] con la base (1,X,,Xn)(1, X, \dots, X^n), la base dual es PP(k)(0)k!P \mapsto \frac{P^{(k)}(0)}{k!} (coeficientes de Taylor). Otra base del dual: las evaluaciones PP(xi)P \mapsto P(x_i) en n+1n + 1 puntos distintos; su base “predual” en Kn[X]K_n[X] es exactamente la familia de los polinomios de Lagrange LiL_i (volumen del primer año), ya que Li(xj)=δijL_i(x_j) = \delta_{ij}. La interpolación es dualidad.

Método 2.3 (Bases duales y anteduales en la práctica)

Para desarrollar una forma φ\varphi sobre una base (ei)(e_i) de EE: las coordenadas son los valores φ(ei)\varphi(e_i); no hay ningún sistema que resolver. Para hallar la base (uj)(u_j) de EE cuya dual es una base dada (φ1,,φn)(\varphi_1, \dots, \varphi_n) de EE^* (la antedual), se resuelven los nn sistemas lineales

φi(uj)=δij(1in),\varphi_i(u_j) = \delta_{ij} \qquad (1 \leq i \leq n),

una columna uju_j cada vez; en términos matriciales, si las filas de MM dan los coeficientes de los φi\varphi_i en una base conocida de EE^*, las columnas de M1M^{-1} son los uju_j. La existencia y la unicidad de la antedual se demuestran en el problema de fin de semana de este capítulo; el cálculo es siempre esta inversión.

Ejemplo 2.4 (Una base dual de R2\R^2, calculada por completo)

Para la base b1=(1,1)b_1 = (1, 1), b2=(1,1)b_2 = (1, -1) de R2\R^2, la base dual (b1,b2)(b_1^*, b_2^*) ha de cumplir bi(bj)=δijb_i^*(b_j) = \delta_{ij}. Escribiendo b1(x,y)=αx+βyb_1^*(x, y) = \alpha x + \beta y, las condiciones α+β=1\alpha + \beta = 1 y αβ=0\alpha - \beta = 0 dan

b1(x,y)=x+y2,y anaˊlogamenteb2(x,y)=xy2.b_1^*(x, y) = \frac{x + y}{2}, \qquad\text{y análogamente}\qquad b_2^*(x, y) = \frac{x - y}{2} .

Comprobaciones de sensatez: b1b_1^* no es e1+e2e_1^* + e_2^* evaluado ingenuamente; la base dual depende de la base entera y no de cada vector por separado (al sustituir b2b_2 por (0,1)(0, 1), b1b_1^* pasa a ser xxx \mapsto x). Y la fórmula de desarrollo funciona: (x,y)=x+y2b1+xy2b2(x, y) = \frac{x+y}2\,b_1 + \frac{x-y}2\,b_2, la descomposición en parte par e impar de un par; las bases duales son extractores de coordenadas, y esta extrae las partes simétrica y antisimétrica.

Definición 2.5 (Anulador)

Para un subespacio FEF \subseteq E, el anulador es

F={φE:φF=0},F^{\circ} = \{\varphi \in E^* : \varphi|_F = 0\},

un subespacio de EE^*.

Teorema 2.6 (Dimensión del anulador)

dimF=dimEdimF\dim F^{\circ} = \dim E - \dim F. Además, FFF \mapsto F^\circ invierte las inclusiones, y FF se recupera a partir de su anulador:

F={xE:φF, φ(x)=0}.F = \{x \in E : \forall\varphi \in F^\circ,\ \varphi(x) = 0\}.

En consecuencia, todo subespacio de dimensión pp en dimensión nn es el conjunto de soluciones de npn - p ecuaciones lineales independientes, y recíprocamente.

Demostración. Elijamos una base (e1,,ep)(e_1, \dots, e_p) de FF completada hasta una base de EE. Una forma φ=φ(ei)ei\varphi = \sum \varphi(e_i) e_i^* anula FF si y solo si se anulan sus pp primeros coeficientes: F=Vect(ep+1,,en)F^\circ = \operatorname{Vect}(e_{p+1}^*, \dots, e_n^*), de dimensión npn - p. La inversión de las inclusiones es inmediata. Para la recuperación: el miembro derecho contiene a FF; recíprocamente, si xFx \notin F, complétese una base de FF con xx y con otros vectores; la forma coordenada de xx en esa base anula FF pero no xx. La lectura en términos de “ecuaciones” se obtiene tomando una base (φ1,,φnp)(\varphi_1, \dots, \varphi_{n-p}) de FF^\circ: entonces F=kerφjF = \bigcap \ker\varphi_j, intersección de npn - p hiperplanos independientes.

Ejemplo 2.7 (Un anulador, en los dos sentidos)

Sea F=Vect((1,2,1), (1,0,1))R3F = \operatorname{Vect}\bigl((1, 2, 1),\ (1, 0, -1)\bigr) \subseteq \R^3. Una forma φ=ae1+be2+ce3\varphi = a\,e_1^* + b\,e_2^* + c\,e_3^* anula FF si y solo si

a+2b+c=0yac=0,a + 2b + c = 0 \qquad\text{y}\qquad a - c = 0 ,

es decir, c=ac = a y b=ab = -a: F=R(e1e2+e3)F^\circ = \R\,(e_1^* - e_2^* + e_3^*), de dimensión 32=13 - 2 = 1, como exige el Teorema 2.6. Leído al revés: F={(x,y,z):xy+z=0}F = \{(x, y, z) : x - y + z = 0\}, el plano recuperado como núcleo de la única forma que genera FF^\circ. Pasar de una familia generadora a unas ecuaciones es calcular un anulador; pasar de las ecuaciones a una parametrización es calcular un preanulador. (Comprobación: ambos vectores generadores cumplen xy+z=0x - y + z = 0.)

Definición 2.8 (Aplicación traspuesta)

Para uL(E,F)u \in \mathcal{L}(E, F), la traspuesta uTL(F,E)u^{\mathsf T} \in \mathcal{L}(F^*, E^*) es

uT(ψ)=ψu.u^{\mathsf T}(\psi) = \psi \circ u .

Cumple (vu)T=uTvT(v \circ u)^{\mathsf T} = u^{\mathsf T} \circ v^{\mathsf T} y, en bases duales, la matriz de uTu^{\mathsf T} es la matriz traspuesta de la de uu, lo que por fin explica la traspuesta del primer año.

Ejemplo 2.9 (La traspuesta, entrada a entrada)

Sea u ⁣:R2R3u \colon \R^2 \to \R^3 de matriz A=(120130)A = \left(\begin{smallmatrix} 1 & 2\\ 0 & 1\\ 3 & 0\end{smallmatrix}\right) en las bases canónicas. Para ψ=b1f1+b2f2+b3f3(R3)\psi = b_1f_1^* + b_2f_2^* + b_3f_3^* \in (\R^3)^*, calculemos uT(ψ)=ψuu^{\mathsf T}(\psi) = \psi \circ u sobre la base de R2\R^2:

(ψu)(e1)=ψ(1,0,3)=b1+3b3,(ψu)(e2)=ψ(2,1,0)=2b1+b2.(\psi \circ u)(e_1) = \psi(1, 0, 3) = b_1 + 3b_3, \qquad (\psi \circ u)(e_2) = \psi(2, 1, 0) = 2b_1 + b_2 .

Luego uT(ψ)=(b1+3b3)e1+(2b1+b2)e2u^{\mathsf T}(\psi) = (b_1 + 3b_3)\,e_1^* + (2b_1 + b_2)\,e_2^*, y en las bases duales la matriz de uTu^{\mathsf T} es

(103210)=AT:\begin{pmatrix} 1 & 0 & 3\\ 2 & 1 & 0\end{pmatrix} = A^{\mathsf T} :

la traspuesta abstracta es la matriz volteada, sin dejar ningún cálculo a la fe. Obsérvese el mecanismo: la columna jj-ésima de AA pasó a ser la fila jj-ésima de la nueva matriz porque ψu\psi \circ u lee las salidas de uu a través de los coeficientes de ψ\psi.

Proposición 2.10

keruT=(imu)\ker u^{\mathsf T} = (\operatorname{im} u)^{\circ} y imuT=(keru)\operatorname{im} u^{\mathsf T} = (\ker u)^{\circ}. En consecuencia, rk(uT)=rk(u)\operatorname{rk}(u^{\mathsf T}) = \operatorname{rk}(u): el rango por filas es igual al rango por columnas, demostrado estructuralmente.

Demostración. ψkeruT    ψu=0    ψ\psi \in \ker u^{\mathsf T} \iff \psi \circ u = 0 \iff \psi anula imu\operatorname{im} u: esa es la primera identidad. Para la segunda: uT(ψ)=ψuu^{\mathsf T}(\psi) = \psi \circ u anula siempre keru\ker u, luego imuT(keru)\operatorname{im} u^{\mathsf T} \subseteq (\ker u)^\circ; las dimensiones coinciden por el teorema del rango y el Teorema 2.6:

rkuT=dimFdimkeruT=dimF(dimFrku)=rku=dim(keru).\operatorname{rk} u^{\mathsf T} = \dim F^* - \dim\ker u^{\mathsf T} = \dim F - \bigl(\dim F - \operatorname{rk} u\bigr) = \operatorname{rk} u = \dim (\ker u)^{\circ} . \qedhere

Ejemplo 2.11 (El rango leído por los dos lados)

Sea

A=(120101111312).A = \begin{pmatrix} 1 & 2 & 0 & 1\\ 0 & 1 & 1 & 1\\ 1 & 3 & 1 & 2 \end{pmatrix} .

Rango por columnas: la tercera fila es la suma de las dos primeras, luego rkA2\operatorname{rk} A \leq 2; las columnas 11 y 22 son libres: rkA=2\operatorname{rk} A = 2. Núcleo de la traspuesta: resolviendo ATy=0A^{\mathsf T}y = 0 se obtiene yR(1,1,1)y \in \R\,(1, 1, -1), de modo que kerAT\ker A^{\mathsf T} tiene dimensión 1=321 = 3 - 2: exactamente (imA)(\operatorname{im} A)^\circ bajo la identificación de (R3)(\R^3)^* con los vectores fila, tal como afirma la Proposición 2.10; la única relación “fila3_3 = fila1_1 + fila2_2es el anulador del espacio de columnas. El rango por filas (22 filas libres) y el rango por columnas coinciden no por casualidad, sino porque ambos valen rkA=rkAT\operatorname{rk} A = \operatorname{rk} A^{\mathsf T}.

Ejemplo 2.12 (La dualidad lee una fórmula de cuadratura)

¿Por qué existe una regla como la de Simpson (Ejercicio 2.4) y por qué es única? La dualidad responde antes de todo cálculo. Sobre E=R2[X]E = \R_2[X], la integral P01PP \mapsto \int_0^1 P es un vector concreto del dual EE^*, de dimensión tres; las evaluaciones en 00, 12\frac12 y 11 forman una base de EE^*; por tanto la integral se desarrolla de manera única sobre ellas, y ese desarrollo es la regla de Simpson, con sus coeficientes incluidos. Un recuento de dimensiones también calibra las expectativas: sobre R3[X]\R_3[X], cuatro dimensiones de formas no pueden en general ser generadas por tres evaluaciones, de modo que la exactitud sobre las cúbicas no se la debemos a la dualidad; que Simpson integre las cúbicas exactamente pese a todo es una simetría de regalo (cancelación de los grados impares en torno a 12\frac12), y hay que comprobarla a mano. Las reglas con n+1n + 1 nodos son desarrollos de la forma integración en una base de evaluaciones de Rn[X]\R_n[X]^*: la existencia y la unicidad cuestan un teorema sobre bases duales; solo los grados de regalo cuestan trabajo.

2.2 Formas multilineales alternadas

Definición 2.13

Una aplicación f ⁣:EnKf \colon E^n \to K es nn-lineal cuando es lineal en cada variable, y alternada cuando se anula siempre que dos argumentos son iguales. Ser alternada implica ser antisimétrica: intercambiar dos argumentos cambia el signo (desarróllese f(,x+y,,x+y,)=0f(\dots, x + y, \dots, x + y, \dots) = 0); más en general, para σSn\sigma \in \mathfrak{S}_n,

f(xσ(1),,xσ(n))=ε(σ)f(x1,,xn),f(x_{\sigma(1)}, \dots, x_{\sigma(n)}) = \varepsilon(\sigma)\, f(x_1, \dots, x_n),

descomponiendo σ\sigma en transposiciones (Teorema 1.21).

Teorema 2.14 (Teorema fundamental de los determinantes)

Sea dimE=n\dim E = n y B=(e1,,en)\mathcal{B} = (e_1, \dots, e_n) una base. El espacio de las formas nn-lineales alternadas sobre EE tiene dimensión 11: toda forma así es múltiplo de

detB(x1,,xn)=σSnε(σ)i=1naσ(i),i,xj=iaijei,\det{}_{\mathcal{B}}(x_1, \dots, x_n) = \sum_{\sigma \in \mathfrak{S}_n} \varepsilon(\sigma) \prod_{i=1}^{n} a_{\sigma(i),\,i}, \qquad x_j = \sum_{i} a_{ij} e_i ,

y detB\det_{\mathcal{B}} es la única que toma el valor 11 sobre B\mathcal{B}.

Demostración. Sea ff una forma nn-lineal alternada. Desarrollando cada argumento sobre B\mathcal{B} por multilinealidad,

f(x1,,xn)=i1,,inai1,1ain,nf(ei1,,ein).f(x_1, \dots, x_n) = \sum_{i_1, \dots, i_n} a_{i_1,1}\cdots a_{i_n,n}\, f(e_{i_1}, \dots, e_{i_n}).

Los términos con un índice repetido se anulan (por ser alternada); las tuplas (i1,,in)(i_1, \dots, i_n) que sobreviven son las inyectivas, es decir, ik=σ(k)i_k = \sigma(k) para una permutación σ\sigma, y la antisimetría reordena f(eσ(1),,eσ(n))=ε(σ)f(e1,,en)f(e_{\sigma(1)}, \dots, e_{\sigma(n)}) = \varepsilon(\sigma) f(e_1, \dots, e_n). Por tanto

f=f(e1,,en)detB:f = f(e_1, \dots, e_n) \cdot \det{}_{\mathcal{B}} :

toda forma alternada es ese múltiplo, siempre que el propio detB\det_{\mathcal{B}} (la suma destacada) sea nn-lineal alternada y tome el valor 11 sobre B\mathcal B. La multilinealidad es clara (cada sumando es lineal en cada columna). Valor sobre B\mathcal B: el único término no nulo es el de σ=id\sigma = \mathrm{id}. Alternada: supongamos xj=xkx_j = x_k (jkj \neq k), de modo que las columnas de coordenadas cumplen aij=aika_{i j} = a_{i k} para todo ii. Emparejemos cada σ\sigma con σ=σ(jk)\sigma' = \sigma\circ(j\,k), una involución sin puntos fijos de Sn\mathfrak{S}_n. Los productos emparejados coinciden:

iaσ(i),i=aσ(k),j  aσ(j),kij,kaσ(i),i=aσ(k),k  aσ(j),jij,kaσ(i),i=iaσ(i),i,\prod_i a_{\sigma'(i),\,i} = a_{\sigma(k),\,j}\; a_{\sigma(j),\,k} \prod_{i \neq j,k} a_{\sigma(i),\,i} = a_{\sigma(k),\,k}\; a_{\sigma(j),\,j} \prod_{i \neq j,k} a_{\sigma(i),\,i} = \prod_i a_{\sigma(i),\,i},

usando la igualdad de las columnas jj y kk; mientras que ε(σ)=ε(σ)\varepsilon(\sigma') = -\varepsilon(\sigma). Cada pareja aporta cero: la suma se anula.

Ejemplo 2.15 (Sarrus, deducida y demolida)

Para n=3n = 3, la fórmula de las permutaciones tiene exactamente 3!=63! = 6 términos. Listando S3\mathfrak{S}_3 por signatura —id\mathrm{id}, (123)(1\,2\,3), (132)(1\,3\,2) pares; (12)(1\,2), (13)(1\,3), (23)(2\,3) impares— se obtiene

detA=a11a22a33+a21a32a13+a31a12a23a21a12a33a31a22a13a11a32a23:\det A = a_{11}a_{22}a_{33} + a_{21}a_{32}a_{13} + a_{31}a_{12}a_{23} - a_{21}a_{12}a_{33} - a_{31}a_{22}a_{13} - a_{11}a_{32}a_{23} :

precisamente la regla de las “diagonales” de Sarrus que se enseña en secundaria, ahora un teorema y con los signos misteriosos identificados como signaturas. La demolición: para n=4n = 4 hay 2424 permutaciones, de las cuales cualquier esquema de diagonales solo recoge 88; Sarrus no tiene versión en grado 44, y el desarrollo por cofactores (Teorema 2.17 (4)) toma el relevo. Contar términos es además una advertencia: la fórmula de las permutaciones tiene n!n! sumandos, así que es una definición, no un algoritmo; el método de eliminación calcula det\det en O(n3)O(n^3) operaciones.

Definición 2.16 (Determinantes)

El determinante de una familia en una base es detB(x1,,xn)\det_{\mathcal{B}}(x_1, \dots, x_n); el determinante de una matriz AA es el determinante de sus columnas en la base canónica — la fórmula de las permutaciones de más arriba—; el determinante de un endomorfismo uu es el escalar detu\det u tal que

detB(u(x1),,u(xn))=detudetB(x1,,xn)para todos los xi\det{}_{\mathcal{B}}\bigl(u(x_1), \dots, u(x_n)\bigr) = \det u \cdot \det{}_{\mathcal{B}}(x_1, \dots, x_n) \quad \text{para todos los } x_i

(el miembro izquierdo es nn-lineal alternado, luego múltiplo de detB\det_\mathcal{B} por el Teorema 2.14; el factor no depende de B\mathcal{B}).

Teorema 2.17 (El cálculo de determinantes, demostrado)

  1. det(uv)=detudetv\det(uv) = \det u\,\det v;   det(AB)=detAdetB\;\det(AB) = \det A \det B.
  2. uu es invertible     detu0\iff \det u \neq 0; una familia es base     \iff su determinante en alguna base es no nulo.
  3. det(AT)=detA\det(A^{\mathsf T}) = \det A.
  4. El desarrollo por cofactores a lo largo de cualquier fila o columna, tal como se enunció en el volumen del primer año, es válido; y dos matrices semejantes tienen el mismo determinante.

Demostración. (1) Aplíquese dos veces la relación de definición: detB(uv(xi))=detudetB(v(xi))=detudetvdetB(xi)\det_{\mathcal B}(uv(x_i)) = \det u \cdot \det_{\mathcal B}(v(x_i)) = \det u \det v \cdot \det_{\mathcal B}(x_i).

(2) Si uu es invertible, detudetu1=detid=10\det u \det u^{-1} = \det \mathrm{id} = 1 \neq 0. Si no lo es, las imágenes u(ei)u(e_i) están ligadas; expresando una a través de las demás y desarrollando, detB(u(ei))=0\det_{\mathcal B}(u(e_i)) = 0 (el carácter alternado mata las direcciones repetidas), luego detu=0\det u = 0. El criterio de base es el mismo enunciado para familias.

(3) En la fórmula de las permutaciones, reindexemos cada producto mediante j=σ(i)j = \sigma(i), es decir, i=τ(j)i = \tau(j) con τ=σ1\tau = \sigma^{-1}: los factores son los mismos números en otro orden, así que

i=1naσ(i),i=j=1naj,τ(j),\prod_{i=1}^{n} a_{\sigma(i),\,i} = \prod_{j=1}^{n} a_{j,\,\tau(j)} ,

y ε(τ)=ε(σ)1=ε(σ)\varepsilon(\tau) = \varepsilon(\sigma)^{-1} = \varepsilon(\sigma) (los valores son ±1\pm1 y ε\varepsilon es un morfismo). Sumar sobre σ\sigma es lo mismo que sumar sobre τ\tau (la inversión es una biyección de Sn\mathfrak{S}_n):

detA=τε(τ)jaj,τ(j)=det(AT),\det A = \sum_{\tau}\varepsilon(\tau)\prod_j a_{j,\tau(j)} = \det(A^{\mathsf T}),

siendo la última suma la fórmula de las permutaciones aplicada a las entradas traspuestas (AT)ij=aji(A^{\mathsf T})_{ij} = a_{ji}.

(4) Fijemos la columna jj y separemos xj=iaijeix_j = \sum_i a_{ij} e_i por linealidad: detA=iaijdet(,ei,)\det A = \sum_i a_{ij}\, \det(\dots, e_i, \dots), y llevar eie_i a la última posición (nin - i transposiciones de filas y njn - j de columnas, vía (3)) identifica det(,ei,)=(1)i+jΔij\det(\dots, e_i, \dots) = (-1)^{i+j}\Delta_{ij} con el menor: exactamente la regla de los cofactores del primer año. Semejanza: det(P1AP)=detP1detAdetP=detA\det(P^{-1}AP) = \det P^{-1}\det A \det P = \det A por (1).

Ejemplo 2.18 (Desarrollo por cofactores, ejecutado)

Calculemos

det(213041120)\det\begin{pmatrix} 2 & 1 & 3\\ 0 & 4 & 1\\ 1 & 2 & 0 \end{pmatrix}

a lo largo de la primera columna (con la pereza que permite un cero: uno). Los signos siguen el tablero (1)i+j(-1)^{i+j}:

2det(4120)0+1det(1341)=2(02)+(112)=15.2\,\det\begin{pmatrix}4 & 1\\ 2 & 0\end{pmatrix} - 0 + 1\cdot\det\begin{pmatrix}1 & 3\\ 4 & 1\end{pmatrix} = 2(0 - 2) + (1 - 12) = -15 .

Comprobación por Sarrus (Ejemplo 2.15): 0+1+01204=150 + 1 + 0 - 12 - 0 - 4 = -15. Estrategia, no dogma: desarróllese por la línea con más ceros y, cuando ninguna tenga, créense primero mediante operaciones elementales; una ronda de eliminación cuesta menos que dos capas de cofactores.

Ejemplo 2.19 (Un determinante mediante las reglas)

Sea JMn(K)J \in \mathcal{M}_n(K) la matriz de unos y aKa \in K; calculemos det(aIn+J)\det(aI_n + J) con las herramientas recién demostradas. Todas las columnas de aIn+JaI_n + J suman igual: sumemos todas las filas a la primera (el determinante no varía, pues sumar un múltiplo de una fila a otra añade un término con una dirección repetida, que la alternancia mata). La primera fila pasa a ser (a+n,a+n,,a+n)(a + n, a + n, \dots, a + n); saquemos a+na + n por linealidad en esa fila y restemos después la primera columna de todas las demás: lo que queda es triangular con diagonal (1,a,,a)(1, a, \dots, a). Por tanto

det(aIn+J)=(a+n)an1.\det(aI_n + J) = (a + n)\,a^{\,n-1}.

La moraleja: las raíces a=0a = 0 (de multiplicidad n1n - 1) y a=na = -n dicen que JJ tiene el valor propio 00 con multiplicidad n1n - 1 y el valor propio nn una vez, es decir, el espectro de la matriz JJ de rango uno, un capítulo antes de tiempo (el Capítulo 3 lo hará sistemático).

Ejemplo 2.20 (Un determinante por la fórmula de las permutaciones)

Para una matriz con muchos ceros la fórmula resulta práctica por sí sola: en

A=(0a0000b0000cd000),A = \begin{pmatrix} 0 & a & 0 & 0\\ 0 & 0 & b & 0\\ 0 & 0 & 0 & c\\ d & 0 & 0 & 0 \end{pmatrix},

la única permutación que selecciona entradas no nulas es el ciclo de longitud 44 σ=(1234)\sigma = (1\,2\,3\,4), que envía la columna 11 a la fila 44, etc.; como ε(σ)=(1)3=1\varepsilon(\sigma) = (-1)^3 = -1, resulta detA=abcd\det A = -abcd. (Compruébese mediante tres intercambios de columnas hasta llegar a una matriz diagonal.)

Ejemplo 2.21 (Un Vandermonde por la fórmula del producto)

Para los nodos 0,1,20, 1, 2 (usados por reglas de cuadratura como la del Ejercicio 2.4), el determinante de Vandermonde del Ejercicio 2.11 se evalúa de un vistazo:

det(111012014)=(10)(20)(21)=2,\det\begin{pmatrix} 1 & 1 & 1\\ 0 & 1 & 2\\ 0 & 1 & 4 \end{pmatrix} = (1 - 0)(2 - 0)(2 - 1) = 2 ,

y por desarrollo directo a lo largo de la primera columna: 1(42)=21\cdot(4 - 2) = 2; coinciden. Que no se anule cuando los nodos son distintos es toda la teoría de la interpolación en un solo determinante: las formas de evaluación PP(ai)P \mapsto P(a_i) son una base del dual exactamente cuando este determinante no es nulo, es decir, siempre que los aia_i sean distintos; el Ejemplo 2.2, cuantificado.

2.3 La traza, revisitada

Proposición 2.22

La traza tr ⁣:Mn(K)K\operatorname{tr} \colon \mathcal{M}_n(K) \to K es la única forma lineal con tr(AB)=tr(BA)\operatorname{tr}(AB) = \operatorname{tr}(BA) y tr(In)=n\operatorname{tr}(I_n) = n (si charK=0\operatorname{char} K = 0); la traza de un endomorfismo está bien definida a través de cualquier representación matricial, y

tr(u)=iei(u(ei))\operatorname{tr}(u) = \sum_{i} e_i^*\bigl(u(e_i)\bigr)

en cualquier base: la dualidad escribe la traza sin recurrir a una base.

Demostración. Las identidades tr(AB)=tr(BA)\operatorname{tr}(AB) = \operatorname{tr}(BA) y la invariancia por cambio de base se demostraron en el primer año. Unicidad: una forma lineal tt con t(AB)=t(BA)t(AB) = t(BA) anula todo conmutador ABBAAB - BA. Afirmamos que los conmutadores generan el hiperplano de traza nula, de dimensión n21n^2 - 1. Bastan dos familias de conmutadores. La regla de multiplicación de las matrices elementales es EabEcd=δbcEadE_{ab}E_{cd} = \delta_{bc}E_{ad}. Para iji \neq j da

EiiEijEijEii=Eij0=EijE_{ii}E_{ij} - E_{ij}E_{ii} = E_{ij} - 0 = E_{ij}

(el segundo producto es EijEii=δjiEii=0E_{ij}E_{ii} = \delta_{ji}E_{ii} = 0, pues jij \neq i): toda EijE_{ij} fuera de la diagonal es un conmutador. Y

EijEjiEjiEij=EiiEjj.E_{ij}E_{ji} - E_{ji}E_{ij} = E_{ii} - E_{jj} .

Las EijE_{ij} (iji \neq j; hay n2nn^2 - n) junto con las E11EjjE_{11} - E_{jj} (j2j \geq 2; hay n1n - 1) son n21n^2 - 1 matrices de traza nula linealmente independientes: generan el hiperplano kertr\ker\operatorname{tr}. Así pues, tt se anula donde lo hace tr\operatorname{tr} y factoriza a través de ella: t=ctrt = c\operatorname{tr}; después, t(I)=nt(I) = n obliga a c=1c = 1. En cuanto a la fórmula destacada: la ii-ésima entrada diagonal de la matriz de uu es precisamente ei(u(ei))e_i^*(u(e_i)).

Observación 2.23 (Errores frecuentes)

(i) El determinante es nn-lineal en las columnas, no lineal en la matriz: det(A+B)detA+detB\det(A + B) \neq \det A + \det B en general, y det(λA)=λndetA\det(\lambda A) = \lambda^n\det A, no λdetA\lambda\det A. (ii) La trasposición invierte los productos: (vu)T=uTvT(vu)^{\mathsf T} = u^{\mathsf T}v^{\mathsf T}; olvidar la inversión arruina todo cálculo con inversas. (iii) El anulador FF^\circ vive en EE^*, no en EE: solo se convierte en el familiar “complemento ortogonal” una vez que un producto escalar identifica EE con EE^* (Capítulo 12); ninguna identificación de ese tipo es canónica. (iv) “El rango por filas es igual al rango por columnas” no significa que el espacio de filas sea igual al de columnas: ambos viven en espacios distintos (KnK^n y KmK^m) y se relacionan mediante la Proposición 2.10, no son iguales. (v) La fórmula de las permutaciones es un instrumento de demostración: para calcular con números, úsense operaciones elementales y cofactores (Ejemplo 2.15).

Ejemplo 2.24 (El emparejamiento traza parte el espacio de matrices)

Sobre M2(R)\mathcal{M}_2(\R), con el emparejamiento A,B=tr(AB)\langle A, B\rangle = \operatorname{tr}(AB) del Ejercicio 2.9, descompongamos M=(1423)M = \left(\begin{smallmatrix}1 & 4\\ 2 & 3\end{smallmatrix}\right) en sus partes simétrica y antisimétrica:

M=S+A,S=12(M+MT)=(1333),A=12(MMT)=(0110).M = S + A, \qquad S = \tfrac12(M + M^{\mathsf T}) = \begin{pmatrix}1 & 3\\ 3 & 3\end{pmatrix}, \qquad A = \tfrac12(M - M^{\mathsf T}) = \begin{pmatrix}0 & 1\\ -1 & 0\end{pmatrix}.

Entonces tr(SA)=tr(3133)=0\operatorname{tr}(SA) = \operatorname{tr} \left(\begin{smallmatrix}-3 & 1\\ -3 & 3\end{smallmatrix}\right) = 0: las dos partes son “ortogonales” para el emparejamiento traza, caso particular del hecho general (demostrado en el problema de fin de semana de este capítulo) de que las matrices antisimétricas forman exactamente el anulador de las simétricas. La dualidad ve la descomposición Mn=SnAn\mathcal{M}_n = \mathcal{S}_n \oplus \mathcal{A}_n antes de que se elija ningún producto escalar.

Observación 2.25 (Perspectivas dentro de este volumen)

Conviene observar cómo las tres construcciones de este capítulo se cambian de traje más adelante. La traspuesta vuelve en el Capítulo 3: uu y uTu^{\mathsf T} comparten valores propios con multiplicidades geométricas iguales (problema de fin de semana de este capítulo, pregunta 15), y por eso los análisis por filas y por columnas de una matriz nunca discrepan. El determinante pasa a ser función de un parámetro en el Capítulo 3 (χu(X)=det(Xidu)\chi_u(X) = \det(X\,\mathrm{id} - u)) y un jacobiano en el Capítulo 20, donde su multilinealidad se convierte en el factor del cambio de variable. La traza siembra los invariantes de semejanza: es el segundo coeficiente de χu\chi_u, la suma de los valores propios y, con el tiempo, la integral de la diagonal en identidades del estilo del Capítulo 14. Un capítulo de álgebra lineal, tres largas sombras.

Observación 2.26 (Dónde se usa este capítulo)

El espacio dual no es una abstracción gratuita: los anuladores y las traspuestas gobiernan la teoría de resolubilidad de los sistemas lineales (el problema de fin de semana de este capítulo demuestra a partir de ellos la alternativa de Fredholm en dimensión finita), los emparejamientos no degenerados reaparecen como la forma polar en el Capítulo 12 y el adjunto en el Capítulo 13, y el determinante construido aquí mueve todo el Capítulo 3. En el volumen del tercer año, esa misma dualidad, transportada a dimensión infinita, se convierte en el teorema de representación de Riesz y en la teoría de Fredholm sobre espacios de Hilbert, con la compacidad ocupando el lugar de los recuentos de dimensiones usados aquí.

2.4 Ejercicios

Ejercicio 2.1

En R3\R^3, sean φ1(x,y,z)=x+y\varphi_1(x,y,z) = x + y, φ2=y+z\varphi_2 = y + z, φ3=x+z\varphi_3 = x + z. Demuestra que (φ1,φ2,φ3)(\varphi_1, \varphi_2, \varphi_3) es una base de (R3)(\R^3)^* y halla la base de R3\R^3 de la que es la dual.

Solución

Solución de Ejercicio 2.1.

Tres formas en un dual de dimensión 33: basta con la libertad. Una relación αφ1+βφ2+γφ3=0\alpha\varphi_1 + \beta\varphi_2 + \gamma\varphi_3 = 0 evaluada en (1,0,0),(0,1,0),(0,0,1)(1,0,0), (0,1,0), (0,0,1) da α+γ=0\alpha + \gamma = 0, α+β=0\alpha + \beta = 0, β+γ=0\beta + \gamma = 0, de donde α=β=γ=0\alpha = \beta = \gamma = 0.

Base predual (u1,u2,u3)(u_1, u_2, u_3): resolvemos φi(uj)=δij\varphi_i(u_j) = \delta_{ij}. Escribiendo uj=(x,y,z)u_j = (x, y, z): para u1u_1, las ecuaciones x+y=1x + y = 1, y+z=0y + z = 0, x+z=0x + z = 0 dan u1=(12,12,12)u_1 = \bigl(\tfrac12, \tfrac12, -\tfrac12\bigr); simétricamente, u2=(12,12,12)u_2 = \bigl(-\tfrac12, \tfrac12, \tfrac12\bigr) y u3=(12,12,12)u_3 = \bigl(\tfrac12, -\tfrac12, \tfrac12\bigr).

Ejercicio 2.2

Calcula mediante la fórmula de las permutaciones los determinantes de

(00a0b0c00),(ab00cd0000ef00gh),\begin{pmatrix} 0 & 0 & a\\ 0 & b & 0\\ c & 0 & 0 \end{pmatrix}, \qquad \begin{pmatrix} a & b & 0 & 0\\ c & d & 0 & 0\\ 0 & 0 & e & f\\ 0 & 0 & g & h \end{pmatrix},

y enuncia la regla para matrices diagonales por bloques que sugiere el segundo.

Solución

Solución de Ejercicio 2.2.

Primera matriz: la única permutación con producto no nulo envía 131 \mapsto 3, 222 \mapsto 2, 313 \mapsto 1; es la transposición (13)(1\,3), de signatura 1-1: el determinante vale abc-abc.

Segunda: una permutación con producto no nulo no puede mezclar los dos bloques (toda entrada que los une es 00), luego se descompone como una permutación de {1,2}\{1,2\} por otra de {3,4}\{3,4\}, y la signatura es el producto de las dos signaturas: la suma se factoriza como

(adbc)(ehfg).(ad - bc)(eh - fg) .

Regla general que esto sugiere (y que es cierta, con la misma demostración): el determinante de una matriz diagonal por bloques es el producto de los determinantes de los bloques.

Ejercicio 2.3

Sea F={(x,y,z,t)R4:x+y=z+t y x=2y}F = \{(x,y,z,t) \in \R^4 : x + y = z + t \text{ y } x = 2y\}. Da una base de FF^\circ y comprueba el Teorema 2.6 sobre las dimensiones.

Solución

Solución de Ejercicio 2.3.

FF viene definido por las dos ecuaciones independientes φ1(x,y,z,t)=x+yzt=0\varphi_1(x,y,z,t) = x + y - z - t = 0 y φ2=x2y=0\varphi_2 = x - 2y = 0: por el Teorema 2.6 leído al revés, F=Vect(φ1,φ2)F^\circ = \operatorname{Vect}(\varphi_1, \varphi_2) —pertenecen a FF^\circ por construcción, son libres (no son proporcionales) y dimF=4dimF=42=2\dim F^\circ = 4 - \dim F = 4 - 2 = 2, ya que dimF=2\dim F = 2 (dos ecuaciones independientes en R4\R^4)—. Base: (φ1,φ2)(\varphi_1, \varphi_2); dimensiones: 2+2=42 + 2 = 4, como exige el teorema.

Ejercicio 2.4 ★★

Sean a0,,ana_0, \dots, a_n puntos distintos de KK y φi ⁣:PP(ai)\varphi_i \colon P \mapsto P(a_i) sobre Kn[X]K_n[X]. Demuestra que (φ0,,φn)(\varphi_0, \dots, \varphi_n) es una base de Kn[X]K_n[X]^*, identifica su base predual y desarrolla la forma P01P(t) ⁣dtP \mapsto \int_0^1 P(t)\,\dd t (para K=RK = \R, n=2n = 2, ai=0,12,1a_i = 0, \frac12, 1) en esa base, reconociendo la regla de Simpson.

Solución

Solución de Ejercicio 2.4.

Los φi\varphi_i son n+1n + 1 formas sobre un espacio de dimensión n+1n + 1: basta con la libertad. Si iλiφi=0\sum_i \lambda_i \varphi_i = 0, evaluando en el polinomio de Lagrange LjL_j de los nodos se obtiene λj=0\lambda_j = 0. La base predual es (L0,,Ln)(L_0, \dots, L_n), puesto que φi(Lj)=Lj(ai)=δij\varphi_i(L_j) = L_j(a_i) = \delta_{ij}.

Para la forma integral con nodos 0,12,10, \frac12, 1 sobre R2[X]\R_2[X]: 01P=iciP(ai)\int_0^1 P = \sum_i c_i P(a_i) con ci=01Lic_i = \int_0^1 L_i. Calculemos: L0=2(X12)(X1)L_0 = 2(X - \tfrac12)(X - 1), 01L0=16\int_0^1 L_0 = \frac16; L1=4X(X1)L_1 = -4X(X-1), 01L1=46\int_0^1 L_1 = \frac46; L2=2X(X12)L_2 = 2X(X - \tfrac12), 01L2=16\int_0^1 L_2 = \frac16. Por tanto

01P=16(P(0)+4P(12)+P(1))(PR2[X]):\int_0^1 P = \frac{1}{6}\Bigl(P(0) + 4P\bigl(\tfrac12\bigr) + P(1)\Bigr) \quad (P \in \R_2[X]) :

la regla de Simpson, exacta sobre las cuadráticas; un enunciado sobre bases duales.

Ejercicio 2.5 ★★

Sea uL(E)u \in \mathcal{L}(E) con dimE=n\dim E = n y rku=1\operatorname{rk} u = 1. Demuestra que u=φ()au = \varphi(\cdot)\, a para cierto vector aa y cierta forma φ\varphi; que tru=φ(a)\operatorname{tr} u = \varphi(a); y que u2=(tru)uu^2 = (\operatorname{tr} u)\, u. Deduce det(I+u)=1+tru\det(I + u) = 1 + \operatorname{tr} u.

Solución

Solución de Ejercicio 2.5.

imu=Ka\operatorname{im} u = Ka para algún a0a \neq 0; entonces u(x)=φ(x)au(x) = \varphi(x)\,a, donde φ(x)\varphi(x) es la coordenada de u(x)u(x) sobre aa, lineal en xx. Traza: complétese a=e1a = e_1 hasta una base; la matriz de uu tiene columnas φ(ej)e1\varphi(e_j)\,e_1, de modo que su única entrada diagonal es φ(e1)=φ(a)\varphi(e_1) = \varphi(a): tru=φ(a)\operatorname{tr} u = \varphi(a). Entonces

u2(x)=φ(x)u(a)=φ(x)φ(a)a=(tru)u(x).u^2(x) = \varphi(x)\, u(a) = \varphi(x)\varphi(a)\, a = (\operatorname{tr} u)\, u(x).

Determinante, en dos casos. Si φ(a)0\varphi(a) \neq 0: tómese cualquier base del hiperplano kerφ\ker\varphi y añádase aa. Entonces uu anula kerφ\ker\varphi (allí u(x)=φ(x)a=0u(x) = \varphi(x)a = 0) y u(a)=φ(a)au(a) = \varphi(a)\,a: la matriz de I+uI + u es diagonal, (1,,1,1+φ(a))(1, \dots, 1,\, 1 + \varphi(a)), luego det(I+u)=1+φ(a)=1+tru\det(I + u) = 1 + \varphi(a) = 1 + \operatorname{tr} u. Si φ(a)=0\varphi(a) = 0: entonces akerφa \in \ker\varphi; tómese una base de kerφ\ker\varphi cuyo primer vector sea aa y añádase un vector bb con φ(b)=1\varphi(b) = 1. Entonces I+uI + u deja fija la base de kerφ\ker\varphi y envía bb+ab \mapsto b + a: triangular con diagonal de unos, det(I+u)=1=1+tru\det(I + u) = 1 = 1 + \operatorname{tr} u. Ambos casos concuerdan con la fórmula.

Ejercicio 2.6 ★★

Demuestra que todo hiperplano de Mn(K)\mathcal{M}_n(K) (n2n \geq 2) contiene una matriz invertible. Indicación: un hiperplano es {M:tr(AM)=0}\{M : \operatorname{tr}(AM) = 0\} para cierta A0A \neq 0 (Ejercicio 2.9). Si AA es escalar, exhibe una matriz invertible de traza nula; en caso contrario, halla una MM invertible que haga que AMAM tenga diagonal nula —una matriz de tipo permutación sirve—.

Solución

Solución de Ejercicio 2.6.

Por el Ejercicio 2.9, el hiperplano es HA={M:tr(AM)=0}H_A = \{M : \operatorname{tr}(AM) = 0\} con A0A \neq 0.

Si A=λIA = \lambda I: HAH_A es el hiperplano de traza nula; la matriz de la permutación cíclica de longitud nn (con unos en las posiciones (i,i+1)(i, i+1) y (n,1)(n, 1)) es invertible (su determinante vale ±1\pm 1 por el cálculo del Ejemplo 2.20) y tiene traza nula.

Si AA no es escalar: hallemos primero una PP invertible tal que B=P1APB = P^{-1}AP tenga una entrada no diagonal bji0b_{ji} \neq 0 (jij \neq i). En efecto, si AA ya tiene una, tómese P=IP = I; si AA es diagonal con dos entradas distintas d1d2d_1 \neq d_2, conjugar por la transvección P=I+E12P = I + E_{12} produce la entrada no diagonal d1d20d_1 - d_2 \neq 0 (calcúlese: P1AP=A+(d1d2)E12P^{-1}AP = A + (d_1 - d_2)E_{12}); y una matriz diagonal con todas las entradas iguales es escalar, caso excluido. Pongamos ahora M=I+tEijM' = I + tE_{ij} con t=tr(B)/bjit = -\operatorname{tr}(B)/b_{ji}: entonces

tr(BM)=trB+tbji=0,\operatorname{tr}(BM') = \operatorname{tr} B + t\,b_{ji} = 0,

y MM' es invertible (triangular con diagonal de unos). Deshaciendo la conjugación, M=PMP1M = PM'P^{-1} es invertible y tr(AM)=tr(BM)=0\operatorname{tr}(AM) = \operatorname{tr}(BM') = 0: MHAM \in H_A.

Ejercicio 2.7 ★★

(Derivada del determinante) Para AMn(R)A \in \mathcal{M}_n(\R), demuestra a partir de la multilinealidad que

 ⁣d ⁣dtt=0det(In+tA)=trA,\frac{\dd}{\dd t}\Big|_{t=0} \det(I_n + tA) = \operatorname{tr} A ,

y deduce det(etA)=ettrA\det(\eu^{tA}) = \eu^{t\operatorname{tr} A} suponiendo la derivabilidad de tdet(etA)t \mapsto \det(\eu^{tA}) y la propiedad de grupo e(s+t)A=esAetA\eu^{(s+t)A} = \eu^{sA}\eu^{tA} (establecida en el Capítulo 16).

Solución

Solución de Ejercicio 2.7.

Por la fórmula de las permutaciones, det(I+tA)\det(I + tA) es un polinomio en tt; su término constante vale 11 (para t=0t = 0). Su coeficiente en tt: desarróllese det\det como forma alternada de las columnas ej+tcj(A)e_j + t\,c_j(A); por multilinealidad, los términos lineales en tt sustituyen exactamente un eje_j por cj(A)c_j(A):

jdet(e1,,cj(A),,en)=jajj=trA,\sum_{j} \det(e_1, \dots, c_j(A), \dots, e_n) = \sum_j a_{jj} = \operatorname{tr} A ,

(el determinante con todas las columnas canónicas salvo cj(A)c_j(A) en el lugar jj recoge la jj-ésima entrada diagonal). Así pues, la derivada en 00 vale trA\operatorname{tr} A.

Sea g(t)=det(etA)g(t) = \det(\eu^{tA}). La propiedad de grupo da g(s+t)=g(s)g(t)g(s + t) = g(s)g(t) (multiplicatividad de det\det), gg es derivable y g(0)=trAg'(0) = \operatorname{tr} A por lo anterior (etA=I+tA+O(t2)\eu^{tA} = I + tA + O(t^2)). Un morfismo derivable (R,+)(R,×)(\R, +) \to (\R^*, \times) cumple g=g(0)gg' = g'(0)\,g (derívese g(s+t)g(s+t) respecto de ss en 00), luego g(t)=ettrAg(t) = \eu^{t\operatorname{tr} A} por la unicidad de las soluciones de y=cyy' = cy con y(0)=1y(0) = 1 (volumen del primer año).

Ejercicio 2.8 ★★

(Circulante 3×33 \times 3) Sea j=e2iπ/3j = \eu^{2\iu\pi/3} y

C=(abccabbca)M3(C).C = \begin{pmatrix} a & b & c\\ c & a & b\\ b & c & a \end{pmatrix} \in \mathcal{M}_3(\C).

Comprueba que las columnas de la matriz de Vandermonde de 1,j,j21, j, j^2 son vectores propios de CC y deduce

detC=(a+b+c)(a+bj+cj2)(a+bj2+cj).\det C = (a + b + c)(a + bj + cj^2)(a + bj^2 + cj).
Solución

Solución de Ejercicio 2.8.

Sea vk=(1,jk,j2k)Tv_k = (1, j^k, j^{2k})^{\mathsf T} para k=0,1,2k = 0, 1, 2. Usando 1+j+j2=01 + j + j^2 = 0 y j3=1j^3 = 1:

Cvk=(a+bjk+cj2kc+ajk+bj2kb+cjk+aj2k)=(a+bjk+cj2k)(1jkj2k),C v_k = \begin{pmatrix} a + b j^k + c j^{2k}\\ c + a j^k + b j^{2k}\\ b + c j^k + a j^{2k} \end{pmatrix} = (a + b j^k + c j^{2k}) \begin{pmatrix} 1\\ j^k\\ j^{2k}\end{pmatrix},

(compruébese la segunda fila: jk(a+bjk+cj2k)=ajk+bj2k+cj3k=c+ajk+bj2kj^k(a + bj^k + cj^{2k}) = aj^k + bj^{2k} + cj^{3k} = c + aj^k + bj^{2k}). Luego vkv_k es un vector propio de valor propio λk=a+bjk+cj2k\lambda_k = a + bj^k + cj^{2k}. Los vkv_k forman una base (Vandermonde de los valores distintos 1,j,j21, j, j^2), así que CC es diagonalizable con esos valores propios y

detC=λ0λ1λ2=(a+b+c)(a+bj+cj2)(a+bj2+cj).\det C = \lambda_0\lambda_1\lambda_2 = (a+b+c)(a + bj + cj^2)(a + bj^2 + cj).

Ejercicio 2.9 ★★★

Demuestra que toda forma lineal tt sobre Mn(K)\mathcal{M}_n(K) es Mtr(AM)M \mapsto \operatorname{tr}(AM) para una única AA: la aplicación Atr(A)A \mapsto \operatorname{tr}(A\,\cdot) es un isomorfismo de Mn(K)\mathcal{M}_n(K) sobre su dual. Deduce de nuevo el enunciado de unicidad de la Proposición 2.22.

Solución

Solución de Ejercicio 2.9.

La aplicación Θ ⁣:Atr(A)\Theta \colon A \mapsto \operatorname{tr}(A\,\cdot) es lineal de Mn(K)\mathcal{M}_n(K) en su dual, entre espacios de la misma dimensión n2n^2: basta con la inyectividad. Si tr(AM)=0\operatorname{tr}(AM) = 0 para toda MM, tómese M=EjiM = E_{ji}: tr(AEji)=aij=0\operatorname{tr}(A E_{ji}) = a_{ij} = 0 para todos i,ji, j, luego A=0A = 0. Así pues, Θ\Theta es un isomorfismo.

Unicidad de la traza (Proposición 2.22): una forma tt que anula todos los conmutadores es tr(A)\operatorname{tr}(A\,\cdot) para cierta AA con tr(A(MNNM))=0\operatorname{tr}(A(MN - NM)) = 0 para todas M,NM, N, es decir, tr((AMMA)N)=0\operatorname{tr}((AM - MA)N) = 0 para toda NN (por ciclicidad), es decir, AM=MAAM = MA para toda MM (por inyectividad de Θ\Theta): AA conmuta con todo y por tanto es escalar (conmutar con todas las EijE_{ij} obliga a que las entradas no diagonales sean 00 y las diagonales iguales), luego t=ctrt = c \operatorname{tr}.

Ejercicio 2.10 ★★★

Sean u,vL(E)u, v \in \mathcal{L}(E) con uvvu=uu \circ v - v \circ u = u. Demuestra que uu es nilpotente. Indicación: prueba que tr(uk)=0\operatorname{tr}(u^k) = 0 para todo k1k \geq 1 (calcula ukvvuku^k v - v u^k por inducción) y usa después el hecho siguiente, que puede demostrarse con las identidades de Newton o por inducción sobre la dimensión: un endomorfismo de un C\C-espacio vectorial cuyas potencias tienen todas traza nula es nilpotente. Trabaja sobre C\C.

Solución

Solución de Ejercicio 2.10.

Trabajemos sobre C\C (una matriz real es nilpotente si y solo si lo es como matriz compleja: la nilpotencia es un=0u^n = 0).

Paso 1: tr(uk)=0\operatorname{tr}(u^k) = 0 para k1k \geq 1. Por inducción, ukvvuk=kuku^k v - v u^k = k\, u^k: para k=1k = 1 es la hipótesis; para el paso inductivo,

uk+1vvuk+1=uk(uvvu)+(ukvvuk)u=uk+1+kuk+1.u^{k+1}v - vu^{k+1} = u^k(uv - vu) + (u^k v - v u^k)u = u^{k+1} + k\,u^{k+1} .

Tomando trazas: 0=tr(ukv)tr(vuk)=ktr(uk)0 = \operatorname{tr}(u^k v) - \operatorname{tr}(vu^k) = k \operatorname{tr}(u^k), luego tr(uk)=0\operatorname{tr}(u^k) = 0.

Paso 2: trazas de potencias nulas implican nilpotencia (sobre C\C). Sean λ1,,λr\lambda_1, \dots, \lambda_r los valores propios no nulos distintos de uu, con multiplicidades m1,,mrm_1, \dots, m_r (en el polinomio característico, que se escinde sobre C\C; véase el Capítulo 3). Las trazas de las potencias son tr(uk)=imiλik\operatorname{tr}(u^k) = \sum_i m_i \lambda_i^k (trigonalícese: la diagonal de la potencia kk-ésima de una matriz triangular está formada por las potencias kk-ésimas). El sistema imiλik=0\sum_i m_i \lambda_i^k = 0 para k=1,,rk = 1, \dots, r es invertible de tipo Vandermonde en las incógnitas miλim_i\lambda_i (matriz (λik1)(\lambda_i^{k-1}) por la diagonal λi\lambda_i, con todos los λi0\lambda_i \neq 0 distintos): todo miλi=0m_i \lambda_i = 0, imposible con mi1m_i \geq 1 salvo que r=0r = 0. Luego uu no tiene ningún valor propio no nulo: su polinomio característico es (X)n(-X)^n y Cayley–Hamilton (Capítulo 3) da un=0u^n = 0, es decir, uu es nilpotente.

Ejercicio 2.11 ★★

(Vandermonde) Para a0,,anKa_0, \dots, a_n \in K, demuestra

det(111a0a1ana0na1nann)=0i<jn(ajai).\det\begin{pmatrix} 1 & 1 & \cdots & 1\\ a_0 & a_1 & \cdots & a_n\\ \vdots & \vdots & & \vdots\\ a_0^n & a_1^n & \cdots & a_n^n \end{pmatrix} = \prod_{0 \leq i < j \leq n} (a_j - a_i).

(Considera el determinante como polinomio en ana_n: identifica su grado, sus raíces y su coeficiente director; después induce.)

Solución

Solución de Ejercicio 2.11.

Escribamos V(a0,,an)V(a_0, \dots, a_n) para el determinante e induzcamos sobre nn; el caso inicial es V(a0)=1V(a_0) = 1. Fijemos a0,,an1a_0, \dots, a_{n-1} y consideremos D(T)=V(a0,,an1,T)D(T) = V(a_0, \dots, a_{n-1}, T), el determinante con última columna (1,T,,Tn)(1, T, \dots, T^n): desarrollando por esa columna, DD es un polinomio de grado n\leq n en TT cuyo coeficiente en TnT^n es el menor V(a0,,an1)V(a_0, \dots, a_{n-1}). Supongamos primero que a0,,an1a_0, \dots, a_{n-1} son distintos. Para cada T=aiT = a_i (i<ni < n) dos columnas coinciden, luego D(ai)=0D(a_i) = 0: con nn raíces distintas y grado n\leq n,

D(T)=V(a0,,an1)i=0n1(Tai),D(T) = V(a_0, \dots, a_{n-1}) \prod_{i=0}^{n-1}(T - a_i),

y tomando T=anT = a_n junto con la hipótesis de inducción se obtiene la fórmula del producto. Si dos de los a0,,an1a_0, \dots, a_{n-1} coinciden, ambos miembros valen 00 (columnas repetidas; factor repetido) y la fórmula se cumple trivialmente.

Ejercicio 2.12 ★★★

Sean A,B,C,DMn(K)A, B, C, D \in \mathcal{M}_n(K) con KK infinito, y supongamos CD=DCCD = DC. Demuestra que

det(ABCD)=det(ADBC).\det\begin{pmatrix} A & B\\ C & D\end{pmatrix} = \det(AD - BC).

(Trata primero el caso DD invertible, multiplicando por la derecha por (I0D1CI)\left(\begin{smallmatrix} I & 0\\ -D^{-1}C & I\end{smallmatrix}\right); después sustituye DD por D+tID + tI y compara dos polinomios en tt.)

Solución

Solución de Ejercicio 2.12.

Caso DD invertible. Multipliquemos por la derecha por la matriz por bloques T=(I0D1CI)T = \left(\begin{smallmatrix} I & 0\\ -D^{-1}C & I\end{smallmatrix}\right), triangular por bloques con diagonal de unos, de detT=1\det T = 1 (su determinante, por la fórmula de las permutaciones, solo recoge los bloques diagonales: la regla por bloques del Ejercicio 2.2):

(ABCD)T=(ABD1CBCDD1CD)=(ABD1CB0D),\begin{pmatrix} A & B\\ C & D\end{pmatrix} T = \begin{pmatrix} A - BD^{-1}C & B\\ C - DD^{-1}C & D\end{pmatrix} = \begin{pmatrix} A - BD^{-1}C & B\\ 0 & D\end{pmatrix},

cuyo determinante vale det(ABD1C)detD=det((ABD1C)D)=det(ADBD1CD)\det(A - BD^{-1}C)\det D = \det\bigl((A - BD^{-1}C)D\bigr) = \det(AD - BD^{-1}CD). Como CD=DCCD = DC, se tiene BD1CD=BCBD^{-1}CD = BC: el determinante es det(ADBC)\det(AD - BC).

Caso general. Sea Dt=D+tID_t = D + tI; entonces sigue siendo CDt=DtCCD_t = D_tC. Ambas funciones

f(t)=det(ABCDt)yg(t)=det(ADtBC)f(t) = \det\begin{pmatrix} A & B\\ C & D_t\end{pmatrix} \qquad\text{y}\qquad g(t) = \det(AD_t - BC)

son polinómicas en tt. El polinomio det(D+tI)\det(D + tI) es mónico de grado nn, luego tiene a lo sumo nn raíces: para todos los tt salvo un número finito, DtD_t es invertible y f(t)=g(t)f(t) = g(t) por el primer caso. Dos polinomios sobre un cuerpo infinito que coinciden en infinitos puntos son iguales: f=gf = g, y t=0t = 0 concluye.

2.5 Problema: la alternativa de Fredholm

¿Cuándo tiene solución el sistema lineal u(x)=bu(x) = b? La respuesta completa es un enunciado de dualidad: exactamente cuando bb es anulado por toda forma lineal que anula la imagen de uu, y esas formas son calculables, pues son el núcleo de la traspuesta. Este problema de fin de semana construye el diccionario completo de la dualidad en dimensión finita (factorización de formas, bidualidad, cálculo de anuladores, traspuesta), demuestra la alternativa de Fredholm en dimensión finita y termina con la forma traza y una caracterización: la traza es el único invariante lineal de la semejanza. En todo el problema, EE y FF son KK-espacios vectoriales de dimensión finita y n=dimEn = \dim E.

Problema 2.1

Problema de fin de semana — la dualidad en dimensión finita y la alternativa de Fredholm

Notación: para SES \subseteq E^*, el preanulador es S={xE:φ(x)=0 para toda φS}S_\circ = \{x \in E : \varphi(x) = 0 \text{ para toda } \varphi \in S\}; los anuladores FF^\circ y las traspuestas uTu^{\mathsf T} son los de la Definición 2.5 y de la Definición 2.8.

Parte I — El lema de factorización. Sean φ1,,φp,φE\varphi_1, \dots, \varphi_p, \varphi \in E^*.

  1. Sea Φ ⁣:EKp\Phi \colon E \to K^p, x(φ1(x),,φp(x))x \mapsto (\varphi_1(x), \dots, \varphi_p(x)). Identifica kerΦ\ker\Phi, prueba que ΦT\Phi^{\mathsf T} envía las formas coordenadas de KpK^p a los φi\varphi_i, y deduce

    dim(kerφ1kerφp)=ndimVect(φ1,,φp).\dim \bigl(\ker\varphi_1 \cap \dots \cap \ker\varphi_p\bigr) = n - \dim \operatorname{Vect}(\varphi_1, \dots, \varphi_p).
  2. (Lema de factorización) Demuestra la equivalencia:

    φVect(φ1,,φp)    kerφ1kerφpkerφ.\varphi \in \operatorname{Vect}(\varphi_1, \dots, \varphi_p) \iff \ker\varphi_1 \cap \dots \cap \ker\varphi_p \subseteq \ker\varphi .
  3. Deduce que (φ1,,φp)(\varphi_1, \dots, \varphi_p) es libre si y solo si ikerφi\bigcap_i \ker\varphi_i tiene dimensión npn - p; y que un subespacio de codimensión pp es intersección de pp hiperplanos, nunca de menos.
  4. En R4\R^4, sean φ1=x+yz\varphi_1 = x + y - z, φ2=y+zt\varphi_2 = y + z - t, ψ=x+2yt\psi = x + 2y - t y ψ=x+y+t\psi' = x + y + t. Decide, mediante el lema de factorización, si ψ\psi y ψ\psi' pertenecen a Vect(φ1,φ2)\operatorname{Vect}(\varphi_1, \varphi_2).
  5. Sobre E=R2[X]E = \R_2[X], prueba que ψ0 ⁣:PP(0)\psi_0 \colon P \mapsto P(0), ψ1 ⁣:PP(1)\psi_1 \colon P \mapsto P(1) y ψ2 ⁣:P01P(t) ⁣dt\psi_2 \colon P \mapsto \int_0^1 P(t)\dd t forman una base de EE^*, calcula la base (P0,P1,P2)(P_0, P_1, P_2) de EE de la que es la dual, y halla el único PR2[X]P \in \R_2[X] con P(0)=1P(0) = 1, P(1)=2P(1) = 2, 01P=32\int_0^1 P = \frac32.

Parte II — Bidualidad y cálculo de anuladores.

  1. Prueba que la aplicación de evaluación J ⁣:EEJ \colon E \to E^{**}, J(x)(φ)=φ(x)J(x)(\varphi) = \varphi(x), es lineal e inyectiva, y por tanto un isomorfismo en dimensión finita.
  2. (Doble anulador) Prueba que J(F)=F:=(F)J(F) = F^{\circ\circ} := (F^\circ)^\circ para todo subespacio FEF \subseteq E: mediante la identificación JJ, el anulador del anulador es el propio subespacio.
  3. Demuestra el cálculo de anuladores: (F+G)=FG(F + G)^\circ = F^\circ \cap G^\circ y (FG)=F+G(F \cap G)^\circ = F^\circ + G^\circ.
  4. Deduce (y vuelve a demostrar directamente) que dos formas no nulas con el mismo núcleo son proporcionales.
  5. (Base antedual) Prueba que para toda base (φ1,,φn)(\varphi_1, \dots, \varphi_n) de EE^* existe una única base (u1,,un)(u_1, \dots, u_n) de EE con φi(uj)=δij\varphi_i(u_j) = \delta_{ij}.

Parte III — El cálculo con traspuestas.

  1. Prueba que uuTu \mapsto u^{\mathsf T} es una biyección lineal de L(E,F)\mathcal{L}(E, F) sobre L(F,E)\mathcal{L}(F^*, E^*), y que (u1)T=(uT)1(u^{-1})^{\mathsf T} = (u^{\mathsf T})^{-1} cuando uu es invertible.
  2. (Naturalidad) Prueba que uTTJE=JFuu^{\mathsf T\mathsf T} \circ J_E = J_F \circ u: mediante los isomorfismos de evaluación, la doble traspuesta es uu.
  3. Prueba que uu es sobreyectiva si y solo si uTu^{\mathsf T} es inyectiva, y que uu es inyectiva si y solo si uTu^{\mathsf T} es sobreyectiva.
  4. Para uL(E)u \in \mathcal{L}(E): un subespacio FF es estable por uu si y solo si FF^\circ es estable por uTu^{\mathsf T}.
  5. Prueba que ker(uTλidE)=(im(uλidE))\ker(u^{\mathsf T} - \lambda\, \mathrm{id}_{E^*}) = \bigl(\operatorname{im}(u - \lambda\, \mathrm{id}_E)\bigr)^\circ, y deduce que uu y uTu^{\mathsf T} tienen los mismos valores propios con las mismas multiplicidades geométricas.

Parte IV — La alternativa de Fredholm.

  1. Demuestra que imu=(keruT)\operatorname{im} u = (\ker u^{\mathsf T})_\circ para uL(E,F)u \in \mathcal{L}(E, F), y deduce la alternativa de Fredholm en dimensión finita: la ecuación u(x)=bu(x) = b tiene solución si y solo si toda ψF\psi \in F^* con uTψ=0u^{\mathsf T}\psi = 0 cumple ψ(b)=0\psi(b) = 0.
  2. Forma matricial: para AMm,n(K)A \in \mathcal{M}_{m,n}(K) y bKmb \in K^m se cumple exactamente una de las dos afirmaciones siguientes: (i) Ax=bAx = b tiene solución; (ii) existe yKmy \in K^m con ATy=0A^{\mathsf T}y = 0 e yTb=1y^{\mathsf T}b = 1. Demuestra tanto el “a lo sumo una” como el “al menos una”.
  3. Halla todos los bR3b \in \R^3 para los que el sistema

    x+y=b1,y+z=b2,x+2y+z=b3x + y = b_1, \qquad y + z = b_2, \qquad x + 2y + z = b_3

    tiene solución, calculando el núcleo de la matriz traspuesta.

  4. (Un problema de Neumann discreto) Sobre E=RnE = \R^n (n3n \geq 3), definamos LL por (Lx)k=xk12(xk1+xk+1)(Lx)_k = x_k - \frac12(x_{k-1} + x_{k+1}), con índices módulo nn. Prueba que LT=LL^{\mathsf T} = L (con las identificaciones canónicas), que kerL\ker L es la recta de los vectores constantes (examina una coordenada máxima), y concluye: Lx=bLx = b tiene solución si y solo si kbk=0\sum_k b_k = 0.

Parte V — La forma traza y el teorema de invariancia. Recordemos del Ejercicio 2.9 que Atr(A)A \mapsto \operatorname{tr}(A\,\cdot) identifica Mn(K)\mathcal{M}_n(K) con su dual. Supongamos charK=0\operatorname{char} K = 0 (por ejemplo, K=Q,R,CK = \Q, \R, \C).

  1. Con esa identificación, prueba que el anulador del subespacio Sn\mathcal{S}_n de las matrices simétricas es el subespacio An\mathcal{A}_n de las antisimétricas, y recíprocamente.
  2. Prueba que el anulador del hiperplano sln={M:trM=0}\mathfrak{sl}_n = \{M : \operatorname{tr} M = 0\} es la recta KInK I_n; equivalentemente, que toda forma lineal que se anula sobre las matrices de traza nula es múltiplo de la traza.
  3. Prueba que toda matriz de Mn(K)\mathcal{M}_n(K) es suma de dos matrices invertibles.
  4. (La traza es el único invariante lineal de la semejanza) Sea tt una forma lineal sobre Mn(K)\mathcal{M}_n(K) con t(PMP1)=t(M)t(PMP^{-1}) = t(M) para toda MM y toda PP invertible. Prueba primero que t(PX)=t(XP)t(PX) = t(XP) para PP invertible, después que t(BX)=t(XB)t(BX) = t(XB) para toda BB, y concluye que t=ctrt = c \operatorname{tr} para cierto cKc \in K.
  5. Prueba que rkur\operatorname{rk} u \leq r si y solo si uu es suma de rr aplicaciones de rango 1\leq 1, es decir, u=i=1rψi()fiu = \sum_{i=1}^{r} \psi_i(\cdot)\,f_i con ψiE\psi_i \in E^*, fiFf_i \in F; deduce rk(u+v)rku+rkv\operatorname{rk}(u + v) \leq \operatorname{rk} u + \operatorname{rk} v.
  6. (Síntesis) Redacta el diccionario demostrado en este problema: subespacios frente a anuladores, sumas frente a intersecciones, aplicaciones frente a traspuestas, resolubilidad frente a ortogonalidad al núcleo traspuesto, traza frente a semejanza. Cita para cada entrada la pregunta que la demostró y di en una frase qué sustituye a los recuentos de dimensiones cuando la dimensión pasa a ser infinita (el volumen del tercer año lo precisa sobre espacios de Hilbert).
Solución

Solución de Problema 2.1.

1. Φ\Phi es lineal y kerΦ=ikerφi\ker\Phi = \bigcap_i \ker\varphi_i (una pp-tupla se anula si y solo si lo hace cada entrada). Para las formas coordenadas εi\varepsilon_i de KpK^p: ΦT(εi)=εiΦ=φi\Phi^{\mathsf T}(\varepsilon_i) = \varepsilon_i \circ \Phi = \varphi_i, luego imΦTVect(φi)\operatorname{im}\Phi^{\mathsf T} \supseteq \operatorname{Vect}(\varphi_i); recíprocamente, imΦT\operatorname{im}\Phi^{\mathsf T} está generado por los ΦT(εi)\Phi^{\mathsf T}(\varepsilon_i) (los εi\varepsilon_i generan (Kp)(K^p)^*). Así pues, rkΦ=rkΦT=dimVect(φ1,,φp)=:r\operatorname{rk}\Phi = \operatorname{rk} \Phi^{\mathsf T} = \dim\operatorname{Vect}(\varphi_1, \dots, \varphi_p) =: r (Proposición 2.10), y el teorema del rango da dimikerφi=nr\dim\bigcap_i\ker\varphi_i = n - r.

2. (\Leftarrow) Conservemos una subfamilia libre maximal, por ejemplo φ1,,φr\varphi_1, \dots, \varphi_r, que genera el mismo espacio (de modo que la hipótesis sigue leyéndose irkerφikerφ\bigcap_{i \leq r}\ker\varphi_i \subseteq \ker\varphi: la intersección sobre todos los ii coincide con la que va hasta iri \leq r, pues cada forma descartada es combinación de las demás). La aplicación Ψ=(φ1,,φr) ⁣:EKr\Psi = (\varphi_1, \dots, \varphi_r) \colon E \to K^r es sobreyectiva (pregunta 1: su rango es rr). Si Ψ(x)=Ψ(y)\Psi(x) = \Psi(y), entonces xykerΨkerφx - y \in \ker\Psi \subseteq \ker\varphi, luego φ(x)=φ(y)\varphi(x) = \varphi(y): φ\varphi factoriza como φ=λΨ\varphi = \lambda \circ \Psi con λ ⁣:KrK\lambda \colon K^r \to K bien definida; λ\lambda es lineal porque Ψ\Psi lo es y es sobreyectiva (para t=Ψ(x)t = \Psi(x), t=Ψ(x)t' = \Psi(x'): λ(t+αt)=φ(x+αx)=λ(t)+αλ(t)\lambda(t + \alpha t') = \varphi(x + \alpha x') = \lambda(t) + \alpha\lambda(t')). Escribiendo λ=ciεi\lambda = \sum c_i \varepsilon_i resulta φ=irciφi\varphi = \sum_{i \leq r} c_i\varphi_i. (\Rightarrow) Si φ=ciφi\varphi = \sum c_i \varphi_i, todo xx que anule cada φi\varphi_i anula φ\varphi.

3. Por la pregunta 1, dimkerφi=nr\dim\bigcap\ker\varphi_i = n - r con r=dimVect(φi)pr = \dim\operatorname{Vect}(\varphi_i) \leq p, y r=pr = p si y solo si la familia es libre. Sea FF un subespacio de codimensión pp: su anulador tiene dimensión pp (Teorema 2.6); una base (φ1,,φp)(\varphi_1, \dots, \varphi_p) de FF^\circ da F=ikerφiF = \bigcap_i\ker\varphi_i (fórmula de recuperación). Con menos hiperplanos no se puede: una intersección de qq hiperplanos tiene dimensión nq>np\geq n - q > n - p por la pregunta 1.

4. Calculemos kerφ1kerφ2\ker\varphi_1 \cap \ker\varphi_2: de x+yz=0x + y - z = 0 e y+zt=0y + z - t = 0, parametrizando por (y,z)(y, z) se obtiene x=zyx = z - y, t=y+zt = y + z, es decir, el plano de los vectores (zy,  y,  z,  y+z)(z - y,\; y,\; z,\; y + z). Sobre él, ψ=x+2yt=(zy)+2y(y+z)=0\psi = x + 2y - t = (z - y) + 2y - (y + z) = 0: por el lema de factorización, ψVect(φ1,φ2)\psi \in \operatorname{Vect}(\varphi_1, \varphi_2) —en efecto, ψ=φ1+φ2\psi = \varphi_1 + \varphi_2—. En cambio, ψ=x+y+t=(zy)+y+(y+z)=y+2z\psi' = x + y + t = (z - y) + y + (y + z) = y + 2z no es idénticamente nula allí (para y=1,z=0y = 1, z = 0 vale 11): ψVect(φ1,φ2)\psi' \notin \operatorname{Vect}(\varphi_1, \varphi_2).

5. Tres formas sobre un espacio de dimensión 33: basta con la libertad. Si aψ0+bψ1+cψ2=0a\psi_0 + b\psi_1 + c\psi_2 = 0, comprobemos en 1,X,X21, X, X^2: a+b+c=0a + b + c = 0, b+c2=0b + \frac c2 = 0, b+c3=0b + \frac c3 = 0; restando las dos últimas se obtiene c=0c = 0, después b=0b = 0 y a=0a = 0. Base antedual: escribiendo P=α+βX+γX2P = \alpha + \beta X + \gamma X^2 y resolviendo ψi(Pj)=δij\psi_i(P_j) = \delta_{ij} (P(0)=αP(0) = \alpha, P(1)=α+β+γP(1) = \alpha + \beta + \gamma, 01P=α+β2+γ3\int_0^1 P = \alpha + \frac\beta2 + \frac\gamma3):

P0=14X+3X2,P1=2X+3X2,P2=6X6X2.P_0 = 1 - 4X + 3X^2, \qquad P_1 = -2X + 3X^2, \qquad P_2 = 6X - 6X^2 .

(Comprobación, por ejemplo: 01P2=32=1\int_0^1 P_2 = 3 - 2 = 1, P2(0)=P2(1)=0P_2(0) = P_2(1) = 0.) El problema de interpolación se resuelve mediante las coordenadas en la base antedual:

P=1P0+2P1+32P2=1+XP = 1\cdot P_0 + 2\cdot P_1 + \tfrac32\, P_2 = 1 + X

(coeficiente en XX: 44+9=1-4 - 4 + 9 = 1; coeficiente en X2X^2: 3+69=03 + 6 - 9 = 0); en efecto, P(0)=1P(0) = 1, P(1)=2P(1) = 2, 01P=32\int_0^1 P = \frac32.

6. Linealidad: para toda φ\varphi, J(x+αy)(φ)=φ(x+αy)=J(x)(φ)+αJ(y)(φ)J(x + \alpha y)(\varphi) = \varphi(x + \alpha y) = J(x)(\varphi) + \alpha J(y)(\varphi), es decir, J(x+αy)=J(x)+αJ(y)J(x + \alpha y) = J(x) + \alpha J(y). Inyectividad: si x0x \neq 0, complétese x=e1x = e_1 hasta una base; la forma coordenada e1e_1^* cumple J(x)(e1)=10J(x)(e_1^*) = 1 \neq 0. Como dimE=dimE=dimE\dim E^{**} = \dim E^* = \dim E, inyectiva implica biyectiva.

7. Inclusión: para xFx \in F y φF\varphi \in F^\circ, J(x)(φ)=φ(x)=0J(x)(\varphi) = \varphi(x) = 0, luego J(F)FJ(F) \subseteq F^{\circ\circ}. Dimensiones (aplicando dos veces el Teorema 2.6):

dimF=dimEdimF=n(ndimF)=dimF=dimJ(F),\dim F^{\circ\circ} = \dim E^* - \dim F^\circ = n - (n - \dim F) = \dim F = \dim J(F),

por ser JJ inyectiva. Por tanto J(F)=FJ(F) = F^{\circ\circ}.

8. Primera identidad: φ\varphi anula F+GF + G si y solo si anula a la vez FF y GG (anula las sumas si y solo si anula las piezas): (F+G)=FG(F+G)^\circ = F^\circ \cap G^\circ. Segunda: la inclusión F+G(FG)F^\circ + G^\circ \subseteq (F \cap G)^\circ es clara (cada sumando anula FGF \cap G). Dimensiones, usando la primera identidad y la fórmula de Grassmann:

dim(F+G)=dimF+dimGdim(FG)=(ndimF)+(ndimG)(ndim(F+G)),\dim(F^\circ + G^\circ) = \dim F^\circ + \dim G^\circ - \dim(F^\circ \cap G^\circ) = (n - \dim F) + (n - \dim G) - \bigl(n - \dim(F + G)\bigr),

que por Grassmann en EE vale ndim(FG)=dim(FG)n - \dim(F \cap G) = \dim(F \cap G)^\circ: hay igualdad.

9. Mediante el lema: kerψkerφ\ker\psi \subseteq \ker\varphi con p=1p = 1 da φVect(ψ)\varphi \in \operatorname{Vect}(\psi), y φ0\varphi \neq 0 hace que el escalar sea no nulo. Directamente: tómese x0x_0 con ψ(x0)0\psi(x_0) \neq 0; todo xx se escribe x=(xψ(x)ψ(x0)x0)+ψ(x)ψ(x0)x0x = \bigl(x - \frac{\psi(x)}{\psi(x_0)}x_0\bigr) + \frac{\psi(x)}{\psi(x_0)} x_0, con el primer término en kerψ=kerφ\ker\psi = \ker\varphi; aplicando φ\varphi: φ(x)=φ(x0)ψ(x0)ψ(x)\varphi(x) = \frac{\varphi(x_0)}{\psi(x_0)}\psi(x).

10. Tomemos la base dual (φ1,,φn)(\varphi_1^*, \dots, \varphi_n^*) de (φ1,,φn)(\varphi_1, \dots, \varphi_n) dentro de EE^{**} (Definición 2.1 aplicada a EE^*) y pongamos uj=J1(φj)u_j = J^{-1}(\varphi_j^*): se obtiene una base de EE (JJ es un isomorfismo, pregunta 6), con φi(uj)=J(uj)(φi)=φj(φi)=δij\varphi_i(u_j) = J(u_j)(\varphi_i) = \varphi_j^*(\varphi_i) = \delta_{ij}. Unicidad: las condiciones φi(uj)=δij\varphi_i(u_j) = \delta_{ij} determinan J(uj)J(u_j) sobre la base (φi)(\varphi_i) y, por tanto, determinan uju_j.

11. Linealidad: (u+αv)Tψ=ψ(u+αv)=uTψ+αvTψ(u + \alpha v)^{\mathsf T}\psi = \psi \circ (u + \alpha v) = u^{\mathsf T}\psi + \alpha\, v^{\mathsf T}\psi. Inyectividad: si u0u \neq 0, tómese xx con u(x)0u(x) \neq 0 y ψ\psi con ψ(u(x))0\psi(u(x)) \neq 0 (el truco de la forma coordenada de la pregunta 6): entonces uTψ0u^{\mathsf T}\psi \neq 0. Los espacios L(E,F)\mathcal{L}(E,F) y L(F,E)\mathcal{L}(F^*, E^*) tienen ambos dimensión dimEdimF\dim E \dim F: es biyectiva. Si uu es invertible, la regla de inversión (vu)T=uTvT(vu)^{\mathsf T} = u^{\mathsf T}v^{\mathsf T} da uT(u1)T=(u1u)T=idEu^{\mathsf T}(u^{-1})^{\mathsf T} = (u^{-1}u)^{\mathsf T} = \mathrm{id}_{E^*} y (u1)TuT=(uu1)T=idF(u^{-1})^{\mathsf T}u^{\mathsf T} = (uu^{-1})^{\mathsf T} = \mathrm{id}_{F^*}, luego (uT)1=(u1)T(u^{\mathsf T})^{-1} = (u^{-1})^{\mathsf T}.

12. Para xEx \in E y ψF\psi \in F^*:

(uTT(JEx))(ψ)=(JEx)(uTψ)=(uTψ)(x)=ψ(u(x))=(JF(u(x)))(ψ).\bigl(u^{\mathsf T\mathsf T}(J_E x)\bigr)(\psi) = (J_E x)\bigl(u^{\mathsf T}\psi\bigr) = (u^{\mathsf T}\psi)(x) = \psi\bigl(u(x)\bigr) = \bigl(J_F(u(x))\bigr)(\psi).

Como ψ\psi es arbitraria, uTTJE=JFuu^{\mathsf T\mathsf T} \circ J_E = J_F \circ u.

13. Por la Proposición 2.10: keruT=(imu)\ker u^{\mathsf T} = (\operatorname{im} u)^\circ, luego uu es sobreyectiva     imu=F    (imu)={0}\iff \operatorname{im} u = F \iff (\operatorname{im}u)^\circ = \{0\} (Teorema 2.6)     uT\iff u^{\mathsf T} es inyectiva. Y imuT=(keru)\operatorname{im} u^{\mathsf T} = (\ker u)^\circ, así que uu es inyectiva     keru={0}    (keru)=E\iff \ker u = \{0\} \iff (\ker u)^\circ = E^*     uT\iff u^{\mathsf T} es sobreyectiva.

14. Si u(F)Fu(F) \subseteq F y φF\varphi \in F^\circ, entonces (uTφ)(x)=φ(u(x))=0(u^{\mathsf T}\varphi)(x) = \varphi(u(x)) = 0 para xFx \in F, luego uTφFu^{\mathsf T}\varphi \in F^\circ. Recíprocamente, si u(F)⊈Fu(F) \not\subseteq F, tómese xFx \in F con u(x)Fu(x) \notin F; por la fórmula de recuperación del Teorema 2.6 existe φF\varphi \in F^\circ con φ(u(x))0\varphi(u(x)) \neq 0: entonces (uTφ)(x)0(u^{\mathsf T}\varphi)(x) \neq 0 pese a que xFx \in F, luego uTφFu^{\mathsf T}\varphi \notin F^\circ y FF^\circ no es estable.

15. uTλidE=(uλidE)Tu^{\mathsf T} - \lambda\,\mathrm{id}_{E^*} = (u - \lambda\,\mathrm{id}_E)^{\mathsf T} (la trasposición es lineal y idT=id\mathrm{id}^{\mathsf T} = \mathrm{id}), de modo que su núcleo es (im(uλid))(\operatorname{im}(u - \lambda\,\mathrm{id}))^\circ (Proposición 2.10), de dimensión

nrk(uλid)=dimker(uλid)n - \operatorname{rk}(u - \lambda\,\mathrm{id}) = \dim\ker(u - \lambda\,\mathrm{id})

por el teorema del rango. En particular, uno de los núcleos es no nulo si y solo si lo es el otro: mismos valores propios y mismas multiplicidades geométricas.

16. Inclusión: si b=u(x)b = u(x) y uTψ=0u^{\mathsf T}\psi = 0, entonces ψ(b)=ψ(u(x))=(uTψ)(x)=0\psi(b) = \psi(u(x)) = (u^{\mathsf T}\psi)(x) = 0; luego imu(keruT)\operatorname{im} u \subseteq (\ker u^{\mathsf T})_\circ. Dimensiones: para un subespacio SFS \subseteq F^* se tiene S=JF1(S)S_\circ = J_F^{-1}(S^\circ) (desplegando: ySy \in S_\circ si y solo si toda ψS\psi \in S anula yy, si y solo si JF(y)SJ_F(y) \in S^\circ), de modo que dimS=dimFdimS\dim S_\circ = \dim F - \dim S. Con S=keruTS = \ker u^{\mathsf T}:

dim(keruT)=dimFdimkeruT=rkuT=rku:\dim(\ker u^{\mathsf T})_\circ = \dim F - \dim\ker u^{\mathsf T} = \operatorname{rk} u^{\mathsf T} = \operatorname{rk} u :

las dimensiones coinciden, luego imu=(keruT)\operatorname{im} u = (\ker u^{\mathsf T})_\circ. Reformulado: bimub \in \operatorname{im} u si y solo si ψ(b)=0\psi(b) = 0 para toda ψ\psi con uTψ=0u^{\mathsf T}\psi = 0; esa es la alternativa de Fredholm.

17. Identifiquemos (Km)(K^m)^* con KmK^m mediante yψyy \mapsto \psi_y, ψy(v)=yTv\psi_y(v) = y^{\mathsf T}v; entonces (uTψy)(x)=yTAx=(ATy)Tx(u^{\mathsf T}\psi_y)(x) = y^{\mathsf T}Ax = (A^{\mathsf T}y)^{\mathsf T}x, luego uTψy=ψATyu^{\mathsf T}\psi_y = \psi_{A^{\mathsf T}y}: la traspuesta es la matriz traspuesta. A lo sumo una: si Ax=bAx = b y ATy=0A^{\mathsf T}y = 0, entonces yTb=yTAx=(ATy)Tx=01y^{\mathsf T}b = y^{\mathsf T}Ax = (A^{\mathsf T}y)^{\mathsf T}x = 0 \neq 1. Al menos una: si (i) falla, la pregunta 16 proporciona una ψy\psi_y con ATy=0A^{\mathsf T}y = 0 e yTb0y^{\mathsf T}b \neq 0; reescálese yy para que valga 11.

18. A=(110011121)A = \left(\begin{smallmatrix} 1 & 1 & 0\\ 0 & 1 & 1\\ 1 & 2 & 1\end{smallmatrix}\right) (tercera fila == primera ++ segunda, luego AA es singular). Resolvamos ATy=0A^{\mathsf T}y = 0: y1+y3=0y_1 + y_3 = 0, y1+y2+2y3=0y_1 + y_2 + 2y_3 = 0, y2+y3=0y_2 + y_3 = 0 dan y1=y2=y3y_1 = y_2 = -y_3: la recta generada por y=(1,1,1)y = (1, 1, -1). Fredholm: hay solución si y solo si yTb=b1+b2b3=0y^{\mathsf T}b = b_1 + b_2 - b_3 = 0, es decir, b3=b1+b2b_3 = b_1 + b_2; visiblemente la condición correcta, puesto que la tercera ecuación es la suma de las dos primeras.

19. La matriz de LL tiene 11 en la diagonal y 12-\frac12 en las posiciones (k,k±1)(k, k\pm1) (módulo nn): es simétrica, luego LT=LL^{\mathsf T} = L con la identificación de la pregunta 17. Núcleo: si Lx=0Lx = 0, entonces cada xk=12(xk1+xk+1)x_k = \frac12(x_{k-1} + x_{k+1}). Sea k0k_0 un índice que maximiza xkx_k; la media de los dos vecinos, ambos xk0\leq x_{k_0}, vale xk0x_{k_0} solo si ambos valen xk0x_{k_0}; propagando a lo largo del ciclo, xx es constante. Recíprocamente, las constantes están en el núcleo. Así pues, kerLT=kerL=R(1,,1)\ker L^{\mathsf T} = \ker L = \R(1, \dots, 1), y la alternativa de Fredholm se lee: Lx=bLx = b tiene solución si y solo si (1,,1)Tb=kbk=0(1,\dots,1)^{\mathsf T} b = \sum_k b_k = 0, la condición discreta de compatibilidad: una “distribución de calor” sobre un anillo es realizable por un potencial si y solo si su flujo total se anula.

20. Si AA es antisimétrica y SS simétrica:

tr(AS)=tr((AS)T)=tr(STAT)=tr(SA)=tr(AS),\operatorname{tr}(AS) = \operatorname{tr}\bigl((AS)^{\mathsf T}\bigr) = \operatorname{tr}(S^{\mathsf T}A^{\mathsf T}) = -\operatorname{tr}(SA) = -\operatorname{tr}(AS),

luego 2tr(AS)=02\operatorname{tr}(AS) = 0 y (si charK2\operatorname{char} K \neq 2) tr(AS)=0\operatorname{tr}(AS) = 0: AnSn\mathcal{A}_n \subseteq \mathcal{S}_n^\circ (identificando el dual con las matrices). Dimensiones: dimSn=n2n(n+1)2=n(n1)2=dimAn\dim\mathcal{S}_n^\circ = n^2 - \frac{n(n+1)}2 = \frac{n(n-1)}2 = \dim\mathcal{A}_n: hay igualdad. Intercambiando los papeles (mismo cálculo), An=Sn\mathcal{A}_n^\circ = \mathcal{S}_n.

21. tr(InM)=trM=0\operatorname{tr}(I_nM) = \operatorname{tr} M = 0 para MslnM \in \mathfrak{sl}_n: la recta KInKI_n está en el anulador, cuya dimensión es n2(n21)=1n^2 - (n^2 - 1) = 1: hay igualdad. Traducido mediante el isomorfismo Atr(A)A \mapsto \operatorname{tr}(A\,\cdot): una forma que se anula sobre sln\mathfrak{sl}_n es tr(λIn)=λtr\operatorname{tr}(\lambda I_n\,\cdot) = \lambda\operatorname{tr}.

22. Sea MMn(K)M \in \mathcal{M}_n(K). El polinomio tdet(MtI)t \mapsto \det(M - tI) es no nulo de grado nn, luego tiene a lo sumo nn raíces; KK tiene característica 00 y por tanto es infinito: tómese λ0\lambda \neq 0 que no sea raíz. Entonces M=(MλI)+λIM = (M - \lambda I) + \lambda I escribe MM como suma de dos matrices invertibles.

23. Paso 1: para PP invertible y XX arbitraria, apliquemos la invariancia a M=XPM = XP: t(P(XP)P1)=t(XP)t(P(XP)P^{-1}) = t(XP), es decir, t(PX)=t(XP)t(PX) = t(XP). Paso 2: fijemos XX; ambos miembros de t(BX)=t(XB)t(BX) = t(XB) son lineales en BB y coinciden sobre las BB invertibles; por la pregunta 22, toda BB es suma de dos invertibles, luego coinciden en todas partes. Paso 3: tt anula todo conmutador XBBXXB - BX; los conmutadores generan sln\mathfrak{sl}_n (visto en la demostración de la Proposición 2.22), así que tt se anula sobre sln\mathfrak{sl}_n y la pregunta 21 da t=ctrt = c\operatorname{tr}. (Recíprocamente, toda ctrc\operatorname{tr} es invariante por semejanza: la traza es el invariante lineal de la semejanza.)

24. Si rku=rr\operatorname{rk} u = r' \leq r: tómese una base (f1,,fr)(f_1, \dots, f_{r'}) de imu\operatorname{im} u y escríbase u(x)=i=1rψi(x)fiu(x) = \sum_{i=1}^{r'} \psi_i(x) f_i; cada coordenada ψi(x)\psi_i(x) de u(x)u(x) es lineal en xx (composición de uu con una forma coordenada), de modo que uu es suma de rrr' \leq r aplicaciones de rango 1\leq 1 (complétese con ceros). Recíprocamente, si u=i=1rψi()fiu = \sum_{i=1}^r \psi_i(\cdot)f_i, entonces imuVect(f1,,fr)\operatorname{im} u \subseteq \operatorname{Vect}(f_1, \dots, f_r): rkur\operatorname{rk} u \leq r. Subaditividad: escríbase uu con rku\operatorname{rk} u términos y vv con rkv\operatorname{rk} v términos; la suma tiene rku+rkv\operatorname{rk} u + \operatorname{rk} v términos, luego rk(u+v)rku+rkv\operatorname{rk}(u + v) \leq \operatorname{rk} u + \operatorname{rk} v.

25. El diccionario: a un subespacio FF le corresponde FF^\circ, de dimensión complementaria (Teorema 2.6), y se vuelve atrás por bidualidad (preguntas 6–7); las sumas se intercambian con las intersecciones (pregunta 8); a una aplicación uu le corresponde uTu^{\mathsf T}, con keruT=(imu)\ker u^{\mathsf T} = (\operatorname{im}u)^\circ, imuT=(keru)\operatorname{im}u^{\mathsf T} = (\ker u)^\circ, rangos iguales, inyectividad y sobreyectividad intercambiadas, y subespacios estables y valores propios que se corresponden (preguntas 11–15); la ecuación u(x)=bu(x) = b tiene solución si y solo si bb es ortogonal a keruT\ker u^{\mathsf T} (preguntas 16–19); y sobre Mn\mathcal{M}_n el emparejamiento traza realiza todo el diccionario de forma concreta, con la traza como único invariante lineal de la semejanza (preguntas 20–23) y el rango como longitud mínima de una descomposición en tensores elementales (pregunta 24). En dimensión infinita los recuentos de dimensiones fallan y se sustituyen por hipótesis de cerrado sobre las imágenes y por la completitud; sobre espacios de Hilbert esto se convierte en el teorema de representación de Riesz y en la teoría de Fredholm para operadores compactos, demostrada honestamente en el volumen del tercer año.

Términos definidos en este capítulo

Ver los 395 términos del glosario