Mathematics · Book 3 · Bachelor Year 1

Matemáticas universitarias — Grado 1

Matemáticas universitarias — Grado 1 · Bachelor Year 1

22Determinantes y sistemas lineales

El determinante condensa en un escalar la respuesta a "¿son estos?" nn vectoriza un base?” — y, geométricamente, mide el volumen ellos engendrado. Lo caracterizamos por sus propiedades (multilineal, alternante, normalizado), calcularlo en las dimensiones 22 y 33 y por expansión de cofactor en general, y ponerlo a funcionar en sistemas lineales, junto con el algoritmo multiuso: eliminación gaussiana.

22.1 el determinante

Teorema 22.1 (Caracterización)

Hay exactamente una aplicación det ⁣:Mn(K)K\det \colon \mathcal{M}_n(K) \to K, visto como una función de las columnas nn, es decir:

  1. lineal in each column (los demás arreglados);
  2. alterno: intercambiar dos columnas cambia el signo (por lo tanto, dos columnas iguales dan 00);
  3. normalizado: detIn=1\det I_n = 1.

Para n=2n = 2 y 33:

abcd=adbc,abcdefghi=aei+bfg+cdhcegbdiafh\begin{vmatrix} a & b\\ c & d\end{vmatrix} = ad - bc, \qquad \begin{vmatrix} a & b & c\\ d & e & f\\ g & h & i\end{vmatrix} = aei + bfg + cdh - ceg - bdi - afh

(la regla 3×33 \times 3 de Sarrus: productos de diagonales descendentes menos los ascendentes).

Demostración. Admitido a este nivel.

Observación 22.2

Para n=2n = 2: expandir por bilinealidad en las columnas canónicas da la fórmula, que a la inversa satisface los axiomas — una completa prueba; n=3n = 3 es idéntico a más términos. El caso general (existencia por la suma sobre permutaciones, unicidad por el mismo expansión) exige la firma de un permutación y se difiere a el segundo año; Usamos libremente los axiomas y las consecuencias. abajo.

La expansión n=2n = 2 al completo, ya que es la plantilla: con columnas C1=ae1+ce2C_1 = a\,e_1 + c\,e_2 y C2=be1+de2C_2 = b\,e_1 + d\,e_2, la bilinealidad da

det(C1,C2)=abdet(e1,e1)+addet(e1,e2)+cbdet(e2,e1)+cddet(e2,e2),\det(C_1, C_2) = ab\det(e_1, e_1) + ad\det(e_1, e_2) + cb\det(e_2, e_1) + cd\det(e_2, e_2),

y la alternancia mata los pares repetidos mientras se voltea det(e2,e1)=det(e1,e2)\det(e_2, e_1) = -\det(e_1, e_2): todo la aplicación colapsa para (adbc)det(e1,e2)=adbc(ad - bc)\det(e_1, e_2) = ad - bc por normalización. Unicidad es visible en el cálculo mismo — los axiomas no dejaron elección en cualquier paso — y esta es exactamente la singularidad escalada hecho utilizado en la prueba de la regla del producto a continuación.

Teorema 22.3 (Propiedades)

Para A,BMn(K)A, B \in \mathcal{M}_n(K):

  1. agregar a una columna un múltiplo de otra no cambia la determinante; multiplicar una columna por λ\lambda la multiplica por λ\lambda(entonces det(λA)=λndetA\det(\lambda A) = \lambda^n \det A);
  2. det(AB)=detAdetB\det(AB) = \det A\, \det B;
  3. AA es invertible     \iff detA0\det A \neq 0     \iff las columnas formar un base de KnK^n; y luego det(A1)=(detA)1\det(A^{-1}) = (\det A)^{-1};
  4. det(AT)=detA\det(A^{\mathsf T}) = \det A — entonces cada regla de columna es también una regla de fila;
  5. el determinante de una matriz triangular es el producto de su entradas diagonales.

Demostración. (1) Por linealidad, det(,Ci+λCj,)=detA+λdet(,Cj,)\det(\dots, C_i + \lambda C_j, \dots) = \det A + \lambda\det(\dots, C_j, \dots) donde el segundo determinante tiene dos columnas iguales: cero.

(2) Repare AA y considere φ(B)=det(AB)\varphi(B) = \det(AB) como una función de las columnas de BB: como ABAB tiene columnas ABjAB_j,φ\varphi es multilineal y alterno en el BjB_j. Admitimos, con Teorema 22.1, su singularidad enunciado en forma escalada: cada multilineal alternante aplicación φ\varphi de las columnas es igual a φ(In)det\varphi(I_n) \cdot \det. Aquí φ(In)=detA\varphi(I_n) = \det A, entonces det(AB)=detAdetB\det(AB) = \det A \cdot \det B.

(3) Si AA es invertible:detAdetA1=detI=1\det A\,\det A^{-1} = \det I = 1, entonces detA0\det A \neq 0 y se cumple la fórmula inversa. SiAA no es invertible, sus columnas están vinculadas (Corolario 20.9 y Proposición 20.2); expresando una columna a través los demás y expandir por linealidad deja a los determinantes con dos columnas iguales: detA=0\det A = 0. El base enunciado es Proposición 19.8.

(4) Admitido con la construcción general (es inmediato al fórmula permutación); lo grabamos para usar operaciones de fila.

(5) Si alguna entrada diagonal desaparece, las primeras columnas kk son vinculado para algunos productos kk(consideraciones de clasificación) y det=0=\det = 0 =. De lo contrario, borre cada columna debajo y a la izquierda por operaciones de tipo (1) — posible en forma triangular — llegando a la diagonal matriz, cuyo determinante es el producto de las entradas por multilinealidad de InI_n.

Ejemplo 22.4 (Las reglas, comprobadas en números.)

Tome A=(1234)A = \begin{pmatrix} 1 & 2\\ 3 & 4\end{pmatrix}(detA=2\det A = -2) y B=(0111)B = \begin{pmatrix} 0 & 1\\ 1 & 1\end{pmatrix}(detB=1\det B = -1). entonces

AB=(2347),det(AB)=1412=2=(2)(1);det(AT)=1324=2=detA.AB = \begin{pmatrix} 2 & 3\\ 4 & 7\end{pmatrix}, \quad \det(AB) = 14 - 12 = 2 = (-2)(-1) ; \qquad \det(A^{\mathsf T}) = \begin{vmatrix} 1 & 3\\ 2 & 4 \end{vmatrix} = -2 = \det A .

Multiplicatividad e invariancia de transposición confirmadas — mientras la aditividad FALSO falla en el mismo par:

det(A+B)=1345=7detA+detB=3.\det(A + B) = \begin{vmatrix} 1 & 3\\ 4 & 5\end{vmatrix} = -7 \neq \det A + \det B = -3 .

Treinta segundos de aritmética de este tipo, después de cualquier determinante Se invoca la identidad, es el seguro contra errores más barato disponible.

Ejemplo 22.5 (Determinantes como áreas)

El paralelogramo abarcado por u=(2,0)u = (2, 0) y v=(1,3)v = (1, 3) tiene base 22 y altura 33: área 66. y

2103=6:\begin{vmatrix} 2 & 1\\ 0 & 3\end{vmatrix} = 6 :

el determinante 2×22\times2is el área firmada del paralelogramo de sus columnas. Los axiomas vuelven a contar la geometría: agregar un múltiplo de una columna a la otra es cortar, deslizar el paralelogramo paralelo a un lado sin cambiar base o altura (operación (1) de Teorema 22.3); escalar una columna escala el área; intercambiando las columnas volteadas la orientación, de donde el signo, det(v,u)=6\det(v, u) = -6. EnR3\R^3 la misma lectura da volúmenes firmados, y det\abs{\det} se convierte en el factor de escala de volumen universal de aplicaciones lineales — el hecho detrás de la fórmula de cambio de variables para múltiples integrales en el volumen Año 2.

Teorema 22.6 (Expansión de cofactor)

Sean AMn(K)A \in \mathcal{M}_n(K) y Δij\Delta_{ij} el determinante de AA con la fila ii y la columna jj eliminadas. Entonces, para cualquier columna fija jj(o fila, por transposición):

detA=i=1n(1)i+jaijΔij.\det A = \sum_{i=1}^{n} (-1)^{i+j}\, a_{ij}\, \Delta_{ij} .

Demostración. Admitido a este nivel.

Ejemplo 22.7

Expandiendo a lo largo de la primera columna:

210121012=2211211012=2×32=4.\begin{vmatrix} 2 & 1 & 0\\ 1 & 2 & 1\\ 0 & 1 & 2 \end{vmatrix} = 2\begin{vmatrix} 2 & 1\\ 1 & 2\end{vmatrix} - 1\begin{vmatrix} 1 & 0\\ 1 & 2\end{vmatrix} = 2 \times 3 - 2 = 4 .

Estrategia: crear ceros primero (operaciones de fila/columna), luego expandir a lo largo de la línea más vacía.

Ejemplo 22.8 (El cofactor inverso, una vez a mano)

Para A=(110011101)A = \begin{pmatrix} 1 & 1 & 0\\ 0 & 1 & 1\\ 1 & 0 & 1\end{pmatrix}:detA=1(1)1(1)+0=2\det A = 1(1) - 1(-1) + 0 = 2. los nueve cofactores (1)i+jΔij(-1)^{i+j}\Delta_{ij} se ensamblan en

Com(A)=(111111111),A1=1detACom(A)T=12(111111111),\operatorname{Com}(A) = \begin{pmatrix} 1 & 1 & -1\\ -1 & 1 & 1\\ 1 & -1 & 1 \end{pmatrix}, \qquad A^{-1} = \frac{1}{\det A}\operatorname{Com}(A)^{\mathsf T} = \frac12\begin{pmatrix} 1 & -1 & 1\\ 1 & 1 & -1\\ -1 & 1 & 1 \end{pmatrix},

la fórmula citada en Ejercicio 22.8. Marque una fila-columna par: (fila 11 de AA)(columna 11 de A1A^{-1})=12(1+1+0)=1= \frac12(1 + 1 + 0) = 1, y contra la columna 22:12(1+1+0)=0\frac12(-1 + 1 + 0) = 0. Nueve determinantes 2×22\times2 para una inversa 3×33\times3: ya en este tamaño, la reducción de filas (Ejercicio 22.3) es más barata — El valor de la fórmula del cofactor es teórico (integralidad en Ejercicio 22.8, diferenciabilidad de la inversa en posteriores volúmenes), no computacionales.

Ejemplo 22.9 (La regla del bloque triangular, en tamaño 44)

Reclamación: det(MN0P)=detMdetP\det\begin{pmatrix} M & N\\ 0 & P\end{pmatrix} = \det M\,\det P para bloques 2×22\times2. Borre el bloque NN por columna operaciones: agregar a las columnas 3,43, 4 combinaciones adecuadas de columnas 1,21, 2 elimina NNcuando MM es invertible (resolver MΛ=NM\Lambda = -N para los coeficientes de combinación Λ\Lambda), dejando det(M00P)\det\begin{pmatrix} M & 0\\ 0 & P\end{pmatrix}; entonces expansión de cofactor a lo largo de la primera columna, dos veces, da detMdetP\det M\det P para esta forma de bloque diagonal. SiMM no es invertible, sus columnas están vinculadas, por lo que las dos primeras columnas de la matriz grande están vinculadas (sus mitades inferiores son cero): ambas los lados desaparecen. La regla se extiende a cualquier tamaño de bloque por el mismo argumento de dos casos — y es el motor de Ejercicio 22.10.

Ejemplo 22.10 (Un determinante 4×44 \times 4, estratégicamente)

Δ=1234234134124123.\Delta = \begin{vmatrix} 1 & 2 & 3 & 4\\ 2 & 3 & 4 & 1\\ 3 & 4 & 1 & 2\\ 4 & 1 & 2 & 3 \end{vmatrix}.

Cada fila suma 1010: la operación C1C1+C2+C3+C4C_1 \leftarrow C_1 + C_2 + C_3 + C_4 hace que la primera columna sea constante y la factorización 1010 deja unos. Entonces LiLiL1L_i \leftarrow L_i - L_1(i2i \geq 2) se borra la primera columna:

Δ=101234011302220111=10113222111=10×16=160,\Delta = 10\begin{vmatrix} 1 & 2 & 3 & 4\\ 0 & 1 & 1 & -3\\ 0 & 2 & -2 & -2\\ 0 & -1 & -1 & -1 \end{vmatrix} = 10\begin{vmatrix} 1 & 1 & -3\\ 2 & -2 & -2\\ -1 & -1 & -1 \end{vmatrix} = 10 \times 16 = 160,

el último determinante 3×33\times3 se expande a lo largo de su primera fila: 1(22)1(22)+(3)(22)=0+4+12=161(2 - 2) - 1(-2 - 2) + (-3)(-2 - 2) = 0 + 4 + 12 = 16. Moraleja: una operación bien elegida (detectar la suma constante de filas) supera dieciséis cofactores.

Método 22.11 (Elegir una estrategia determinante)

Escanee la matriz antes de calcular nada.

  1. Sumas constantes de filas o columnas: agrega todo en una línea, factoriza el valor común (Ejemplo 22.10, Ejercicio 22.7).
  2. Estructura repetitiva: restar filas vecinas o columnas para crear ceros; Los patrones de escaleras colapsan hacia la forma triangular, cuyo determinante se lee en la diagonal.
  3. Ceros aislados: expandir a lo largo de la línea más vacía (Ejemplo 22.7); familias recursivas (tridiagonal, Ejercicio 22.6) produce recurrencias este manera.
  4. un parametro: el determinante es un polinomio en eso; encontrar sus raíces detectando los valores degenerados (filas iguales, columnas proporcionales), luego fije el polinomio descendido por grado y coeficiente principal. Para la matriz de Ejercicio 22.7: m=1m = 1 da tres filas iguales (rango 11, una raíz doble),m=2m = -2 hace que las filas suman cero (una raíz más); el determinante tiene grado 33 en mm con término principalm3-m^3(el producto anti-diagonal mmmm\cdot m\cdot m, cuyo Sarrus el signo es 1-1), por lo que debe ser (m+2)(m1)2-(m+2)(m-1)^2 — no Se necesita expansión y los dos métodos se controlan entre sí.

Ejemplo 22.12 (determinante de Vandermonde)

Para escalares x1,,xnx_1, \dots, x_n:

V(x1,,xn)=1x1x12x1n11x2x22x2n11xnxn2xnn1=1i<jn(xjxi).V(x_1, \dots, x_n) = \begin{vmatrix} 1 & x_1 & x_1^2 & \cdots & x_1^{n-1}\\ 1 & x_2 & x_2^2 & \cdots & x_2^{n-1}\\ \vdots & & & & \vdots\\ 1 & x_n & x_n^2 & \cdots & x_n^{n-1} \end{vmatrix} = \prod_{1 \leq i < j \leq n} (x_j - x_i) .

Bosquejo de prueba (detallado en Ejercicio 22.5): operaciones de columna CkCkx1Ck1C_k \leftarrow C_k - x_1 C_{k-1} desde la derecha despeja el primero fila, y factorizar cada fila restante se reduce a V(x2,,xn)V(x_2, \dots, x_n). Distinto de cero si el xix_i son distintos por pares — el determinante detrás de interpolación de Lagrange (Ejemplo 20.10).

22.2 Sistemas lineales

Definición 22.13

Un sistema lineal de ecuaciones nn en pp las incógnitas son AX=BAX = B con AMn,p(K)A \in \mathcal{M}_{n,p}(K),BKnB \in K^n; es homogéneo cuando B=0B = 0. Su solución conjunto, cuando no vacío, es X0+kerAX_0 + \ker A: una solución particular más la general solución homogénea — un afín subespacio de dimensión prkAp - \operatorname{rk} A(rango - nulidad).

Ejemplo 22.14 (La estructura afín, hecha visible)

resolver

{x+y+z=3xy+2z=2.\begin{cases} x + y + z = 3\\ x - y + 2z = 2 . \end{cases}

Restando las ecuaciones: 2yz=12y - z = 1, entonces z=2y1z = 2y - 1 y x=3yz=43yx = 3 - y - z = 4 - 3y. Las soluciones forman la línea.

(x,y,z)=(4, 0, 1)+y(3, 1, 2)(yR):(x, y, z) = (4,\ 0,\ -1) + y\,(-3,\ 1,\ 2) \qquad (y \in \R):

la solución particular X0=(4,0,1)X_0 = (4, 0, -1)(la eleccióny=0y = 0) más la línea del kernel kerA=Vect(3,1,2)\ker A = \operatorname{Vect}(-3, 1, 2) de el sistema homogéneo asociado — comprobar: (3)+1+2=0(-3) + 1 + 2 = 0 y (3)1+4=0(-3) - 1 + 4 = 0. Geométricamente, dos planos no paralelos de R3\R^3 se cruza a lo largo de una línea, y el recuento de dimensiones prkA=32=1p - \operatorname{rk} A = 3 - 2 = 1 lo sabía antes de que resolviéramos cualquier cosa. Cambiar la solución particular (digamos y=1y = 1:X0=(1,1,1)X_0' = (1, 1, 1)) cambia la descripción, no la línea: un afín subespacio tiene muchos orígenes y una dirección.

Teorema 22.15 (Sistemas Cramer cuadrados)

Si AGLn(K)A \in GL_n(K), el sistema AX=BAX = B tiene la solución única X=A1BX = A^{-1}B, cuyos coordenadas son

xj=detAjdetA,Aj=A with column j replaced by B.x_j = \frac{\det A_j}{\det A}, \qquad A_j = A \text{ with column } j \text{ replaced by } B .

Demostración. La unicidad y la existencia son la invertibilidad. Para la fórmula: escriba B=kxkCkB = \sum_k x_k C_k(columnas de AA); entonces, por multilinealidad y alternancia,

detAj=det(C1,,kxkCk,,Cn)=kxkdet(C1,,Ck,,Cn)=xjdetA,\det A_j = \det\Bigl(C_1, \dots, \sum_k x_k C_k, \dots, C_n\Bigr) = \sum_k x_k \det(C_1, \dots, C_k, \dots, C_n) = x_j \det A ,

todos los términos excepto k=jk = j tienen una columna repetida.

Ejemplo 22.16 (Cramer con un parámetro, completo)

Para mRm \in \R, resuelva

{x+my=1mx+y=2.\begin{cases} x + m y = 1\\ m x + y = 2 . \end{cases}

El determinante es 1m21 - m^2. Caso genéricom±1m \neq \pm1: Cramer da

x=1m211m2=12m1m2,y=11m21m2=2m1m2,x = \frac{\begin{vmatrix} 1 & m\\ 2 & 1\end{vmatrix}}{1 - m^2} = \frac{1 - 2m}{1 - m^2}, \qquad y = \frac{\begin{vmatrix} 1 & 1\\ m & 2\end{vmatrix}}{1 - m^2} = \frac{2 - m}{1 - m^2},

una solución limpia para cada mm admisible (verificar en m=0m = 0: (1,2)(1, 2), obviamente cierto). Casos degenerados: en m=1m = 1 las ecuaciones dicen x+y=1x + y = 1 y x+y=2x + y = 2: incompatibles; en m=1m = -1 leen xy=1x - y = 1 y x+y=2-x + y = 2, es decir xy=1x - y = 1 y xy=2x - y = -2: incompatibles nuevamente. La desaparición del determinante anuncia que algo degenera, pero nunca dice qué — vacío o infinito debe decidirse mirando en el lado derecho. Observe también cómo las fórmulas señalan su límites propios: como m1m \to 1^{-},x=12m1m2x = \frac{1 - 2m}{1 - m^2} \to -\infty; el punto de solución se escapa cuando las dos líneas se vuelven paralelo.

Método 22.17 (Eliminación gaussiana en sistemas)

Reducir por filas la matriz aumentada (AB)(A \mid B) a escalón formulario.

  1. Si aparece un pivote en la última columna (línea 0=10 = 1): no solución.
  2. De lo contrario, las incógnitas se dividen en incógnitas de pivote y libre unknowns (parámetros); sustitución hacia atrás expresa lo primero en lo segundo: la solución conjunto es una afín subespacio de dimensión == número de libre incógnitas.

Las fórmulas de Cramer son para teoría y sistemas pequeños; la eliminación es El algoritmo práctico.

Ejemplo 22.18 (Una discusión con el parámetro)

Para mRm \in \R, considere

{x+y+mz=1x+my+z=1mx+y+z=1.\begin{cases} x + y + mz = 1\\ x + my + z = 1\\ mx + y + z = 1 . \end{cases}

La matriz tiene el determinante (m+2)(m1)2-(m+2)(m-1)^2(calculado en Ejercicio 22.7 sumando todas las columnas a la primera). Para m1,2m \neq 1, -2: solución única x=y=z=1m+2x = y = z = \frac{1}{m+2}(por simetría). Para m=1m = 1: una ecuación repetida tres veces, un plano de soluciones. Para m=2m = -2: la suma de las tres ecuaciones da 0=30 = 3, sin solución.

Observación 22.19 (Errores comunes)

The determinant is not lineal in the matrix: det(A+B)detA+detB\det(A + B) \neq \det A + \det B(ya det(I2+I2)=42\det(I_2 + I_2) = 4 \neq 2); eso es lineal en cada columna por separado, que es un diferente cosa por completo. Escalada: det(λA)=λndetA\det(\lambda A) = \lambda^n\det A, no λdetA\lambda\det A— cada una de las columnas nn está escalada. Operaciones de fila are not all libre: LiLi+λLjL_i \leftarrow L_i + \lambda L_j preserva el determinante, pero un intercambio cambia el firmar y LiλLiL_i \leftarrow \lambda L_i lo multiplica por λ\lambda — Los errores contables aquí son la fuente clásica de errores. signos en cálculos basados en eliminación. Un determinante cero es el comienzo, no el final.: dice “rango <n< n” pero no qué rango; sólo trabajo adicional (escalón formulario, o los menores de Ejercicio 22.12) lo localiza — cf. el caso m=1m = 1 vs m=2m = -2 en Ejemplo 22.18. Cramer necesita invertibilidad: cuando detA=0\det A = 0 el Las fórmulas xj=detAj/detAx_j = \det A_j/\det A no tienen sentido y el sistema perfectamente puede tener (infinitas) soluciones. Sólo las matrices cuadradas tienen determinantes: para un rectangular sistema, la eliminación es la única herramienta.

Observación 22.20 (Adónde van los determinantes)

Tres vidas aguardan a este escalar. Geométrico: det\abs{\det} es el factor de escala de área o volumen del aplicación asociado — hecho preciso para el plano en Capítulo 23 y, como el jacobiano de un cambio de variables, en el múltiplo integrales del Año 2 volumen. Algebraico: det(AλI)\det(A - \lambda I), el característica polinomio, abre teoría de valores propios en el año 2 — la identidad A2(trA)A+(detA)I=0A^2 - (\operatorname{tr} A)A + (\det A)I = 0 de el problema del fin de semana Capítulo 21 es su primera sombra. Analítico: determinantes de matrices especiales (Vandermonde, Cauchy, Gram) deciden cuando la interpolación, la descomposición y proyección los problemas están bien planteados; el problema del fin de semana a continuación evalúa las dos primeras familias por completo.

Observación 22.21 (Perspectivas dentro del Libro 3)

Este capítulo cierra la columna vertebral de álgebra lineal del volumen y sus dos capítulos restantes cobran los dividendos. en Capítulo 23: la matriz de Gram (vi,vj)\bigl(\langle v_i, v_j\rangle\bigr) prueba la libertad mediante un determinante (Ejercicio 23.11), y las isometrías planas dividirse en rotaciones y reflexiones a lo largo del signo de su determinante — la clasificación del problema del fin de semana allí corre sobre él. En Capítulo 25: la cantidad de Monge rts2rt - s^2 es el determinante de la matriz simétrica de segunda derivados, y las ecuaciones normales de mínimos cuadrados son una Sistema Cramer cuya matriz es una matriz de Gram (por lo tanto, momento) — invertible precisamente según los criterios con sabor a Vandermonde establecido aquí. Cuando esos capítulos afirman "invertible" o “positivo”, los recibos están en este.

22.3 Ceremonias

Ejercicio 22.1

Calcular:

3152,123456789,111124139.\begin{vmatrix} 3 & 1\\ 5 & 2 \end{vmatrix}, \qquad \begin{vmatrix} 1 & 2 & 3\\ 4 & 5 & 6\\ 7 & 8 & 9\end{vmatrix}, \qquad \begin{vmatrix} 1 & 1 & 1\\ 1 & 2 & 4\\ 1 & 3 & 9\end{vmatrix}.
Solución

Solución de Ejercicio 22.1.

3×21×5=13 \times 2 - 1 \times 5 = 1.

Segundo: L2L2L1L_2 \leftarrow L_2 - L_1,L3L3L2L_3 \leftarrow L_3 - L_2(en las filas originales) dan filas (1,2,3),(3,3,3),(3,3,3)(1,2,3), (3,3,3), (3,3,3): dos iguales filas, determinante 00. (Sarrus confirma:45+84+961054872=045 + 84 + 96 - 105 - 48 - 72 = 0.)

Tercero: es Vandermonde con x=1,2,3x = 1, 2, 3 (Ejemplo 22.12): (21)(31)(32)=2(2-1)(3-1)(3-2) = 2.

Ejercicio 22.2

¿Para qué λR\lambda \in \R es la familia ((1,1,λ),(1,λ,1),(λ,1,1))\bigl((1, 1, \lambda), (1, \lambda, 1), (\lambda, 1, 1)\bigr) una base de R3\R^3?

Solución

Solución de Ejercicio 22.2.

El determinante es igual (suma todas las columnas al primer factor) (λ+2)(\lambda + 2) veces

11λ1λ1111=(λ1)2\begin{vmatrix} 1 & 1 & \lambda\\ 1 & \lambda & 1\\ 1 & 1 & 1 \end{vmatrix} = -(\lambda - 1)^2

(borrar con L1L1L3L_1 \leftarrow L_1 - L_3,L2L2L3L_2 \leftarrow L_2 - L_3 y expandir), dando det=(λ+2)(λ1)2\det = -(\lambda+2)(\lambda-1)^2. Base     det0    λ{1,2}\iff \det \neq 0 \iff \lambda \notin \{1, -2\}.

Ejercicio 22.3

Resuelva por la regla de Cramer:

{2x+y=53x2y=4,then{x+y+z=6xy+z=22x+yz=1.\begin{cases} 2x + y = 5\\ 3x - 2y = 4 , \end{cases} \qquad\text{then}\qquad \begin{cases} x + y + z = 6\\ x - y + z = 2\\ 2x + y - z = 1 . \end{cases}
Solución

Solución de Ejercicio 22.3.

Primer sistema: det=7\det = -7; x=175142=147=2x = \frac{1}{-7}\begin{vmatrix} 5 & 1\\ 4 & -2\end{vmatrix} = \frac{-14}{-7} = 2, y=172534=77=1y = \frac{1}{-7}\begin{vmatrix} 2 & 5\\ 3 & 4\end{vmatrix} = \frac{-7}{-7} = 1. Verificar:2(2)+1=52(2) + 1 = 5;3(2)2=43(2) - 2 = 4.

Segundo sistema: después de L2L1L_2 - L_1 y L32L1L_3 - 2L_1, las filas pasan a ser (1,1,1)(1,1,1),(0,2,0)(0,-2,0),(0,1,3)(0,-1,-3), entonces

detA=111111211=1×2013=6.\det A = \begin{vmatrix} 1&1&1\\ 1&-1&1\\ 2&1&-1\end{vmatrix} = 1 \times \begin{vmatrix} -2 & 0\\ -1 & -3\end{vmatrix} = 6 .

Cramer, reemplazando columnas por (6,2,1)T(6,2,1)^{\mathsf T}:

x=66=1,y=126=2,z=186=3x = \frac{6}{6} = 1, \qquad y = \frac{12}{6} = 2, \qquad z = \frac{18}{6} = 3

(los numeradores se calculan de la misma manera). Verificar: 1+2+3=61 + 2 + 3 = 6;12+3=21 - 2 + 3 = 2;2+23=12 + 2 - 3 = 1.

Ejercicio 22.4

Resuelva por eliminación gaussiana, describiendo la solución conjunto:

{x+2yz+t=12x+4y+zt=5x+2y+2z2t=4.\begin{cases} x + 2y - z + t = 1\\ 2x + 4y + z - t = 5\\ x + 2y + 2z - 2t = 4 . \end{cases}
Solución

Solución de Ejercicio 22.4.

Reducir la matriz aumentada: L2L22L1L_2 \leftarrow L_2 - 2L_1,L3L3L1L_3 \leftarrow L_3 - L_1:

(121110033300333)(121110011100000).\begin{pmatrix} 1 & 2 & -1 & 1 & 1\\ 0 & 0 & 3 & -3 & 3\\ 0 & 0 & 3 & -3 & 3 \end{pmatrix} \to \begin{pmatrix} 1 & 2 & -1 & 1 & 1\\ 0 & 0 & 1 & -1 & 1\\ 0 & 0 & 0 & 0 & 0 \end{pmatrix}.

Pivote desconocido x,zx, z; libre desconocidosy,ty, t. Sustitución hacia atrás:z=1+tz = 1 + t,x=12y+zt=22yx = 1 - 2y + z - t = 2 - 2y. Solución conjunto:

{(22y,  y,  1+t,  t):y,tR}=(2,0,1,0)+Vect((2,1,0,0),(0,0,1,1)),\{(2 - 2y,\; y,\; 1 + t,\; t) : y, t \in \R\} = (2, 0, 1, 0) + \operatorname{Vect}\bigl((-2,1,0,0),\, (0,0,1,1)\bigr),

un plano afín (dimensión 2=4rk22 = 4 - \operatorname{rk} 2) de R4\R^4.

Ejercicio 22.5 ★★

Demuestre la fórmula de Vandermonde de Ejemplo 22.12 mediante inducción en nn, con las operaciones de columna CkCkx1Ck1C_k \leftarrow C_k - x_1 C_{k-1} realizadas desde k=nk = n hasta k=2k = 2.

Solución

Solución de Ejercicio 22.5.

Inducción; n=1n = 1 es el producto vacío=1= 1. Para el paso, realice CkCkx1Ck1C_k \leftarrow C_k - x_1 C_{k-1} para k=n,n1,,2k = n, n-1, \dots, 2(en este orden, por lo que cada operación utiliza una columna aún no modificada). el la primera fila pasa a ser (1,0,,0)(1, 0, \dots, 0); en la fila i2i \geq 2, el kk-ésimo la entrada se convierte en xik1x1xik2=xik2(xix1)x_i^{k-1} - x_1 x_i^{k-2} = x_i^{k-2}(x_i - x_1). Expandiendo a lo largo de la primera fila y factorizando (xix1)(x_i - x_1) de cada fila ii:

V(x1,,xn)=i=2n(xix1)V(x2,,xn),V(x_1, \dots, x_n) = \prod_{i=2}^{n} (x_i - x_1)\cdot V(x_2, \dots, x_n),

y la hipótesis de inducción completa el producto i<j(xjxi)\prod_{i<j}(x_j - x_i).

Ejercicio 22.6 ★★

(Tridiagonal) Sea DnD_n el determinante n×nn \times n con 22 en la diagonal, 11 en las dos diagonales adyacentes,00 en el resto. Expandiendo a lo largo de la primera fila, demuestre Dn=2Dn1Dn2D_n = 2D_{n-1} - D_{n-2} y calcular DnD_n(D1=2D_1 = 2,D2=3D_2 = 3).

Solución

Solución de Ejercicio 22.6.

Expandiendo DnD_na lo largo de la primera fila:Dn=2Dn1110Dn2-blockD_n = 2 D_{n-1} - 1\cdot\begin{vmatrix} 1 & \ast\\ 0 & D_{n-2}\text{-block} \end{vmatrix}; el segundo determinante, expandido a lo largo de su primera columna, es Dn2D_{n-2}. Entonces Dn=2Dn1Dn2D_n = 2D_{n-1} - D_{n-2}, es decir DnDn1=Dn1Dn2D_n - D_{n-1} = D_{n-1} - D_{n-2}: las diferencias son constantes, iguales a D2D1=1D_2 - D_1 = 1. Por lo tanto Dn=D1+(n1)=n+1D_n = D_1 + (n - 1) = n + 1. (Compruebe: D2=3D_2 = 3, y el caso 3×33\times3 es Ejemplo 22.7: D3=4D_3 = 4.)

Ejercicio 22.7 ★★

Completar Ejemplo 22.18: calcular el determinante 11m1m1m11\begin{vmatrix} 1 & 1 & m\\ 1 & m & 1\\ m & 1 & 1\end{vmatrix} por la operación C1C1+C2+C3C_1 \leftarrow C_1 + C_2 + C_3, y realizar la discusión completa del sistema.

Solución

Solución de Ejercicio 22.7.

C1C1+C2+C3C_1 \leftarrow C_1 + C_2 + C_3 hace que la primera columna sea constante (m+2)(m+2); factorizarlo:

det=(m+2)11m1m1111=L1L3, L2L3(m+2)00m10m10111=(m+2)((m1)2)\det = (m+2)\begin{vmatrix} 1 & 1 & m\\ 1 & m & 1\\ 1 & 1 & 1 \end{vmatrix} \overset{L_1 - L_3,\ L_2 - L_3}{=} (m+2)\begin{vmatrix} 0 & 0 & m-1\\ 0 & m-1 & 0\\ 1 & 1 & 1 \end{vmatrix} = (m+2)\cdot\bigl(-(m-1)^2\bigr)

(expandir a lo largo de la primera columna: la entrada única 11 lleva el signo ++, y el determinante restante 2×22 \times 2 es 00(m1)(m1)=(m1)20 \cdot 0 - (m-1)(m-1) = -(m-1)^2).

Discusión. m{1,2}m \notin \{1, -2\}: Sistema Cramer; por la simetría de las ecuaciones, x=y=zx = y = z, y cada ecuación da (m+2)x=1(m + 2)x = 1: solución única (1m+2,1m+2,1m+2)\bigl(\frac{1}{m+2}, \frac{1}{m+2}, \frac{1}{m+2}\bigr).m=1m = 1: las tres ecuaciones se leen x+y+z=1x + y + z = 1: las soluciones forman el plano afín x+y+z=1x + y + z = 1.m=2m = -2: sumando las tres ecuaciones se obtiene 0=30 = 3: solución vacía conjunto.

Ejercicio 22.8 ★★

Sea AMn(R)A \in \mathcal{M}_n(\R) con entradas entero. probar que AA tiene una inversa con entradas enteras si y sólo si detA=±1\det A = \pm 1. (For the direct sense, take determinants; for the converse, admit — or prove for n3n \leq 3 via cofactors — that A1=1detACom(A)TA^{-1} = \frac{1}{\det A}\,\operatorname{Com}(A)^{\mathsf T} with integer cofactor matrix.)

Solución

Solución de Ejercicio 22.8.

(\Rightarrow) SiA1A^{-1} tiene entradas enteras:detAdetA1=1\det A \cdot \det A^{-1} = 1 con ambos determinantes enteros (sumas de productos de entradas): dos números enteros con producto 11 son ambos ±1\pm1.

(\Leftarrow) La fórmula del cofactor A1=1detACom(A)TA^{-1} = \frac{1}{\det A}\operatorname{Com}(A)^{\mathsf T}(comprobada para n3n \leq 3 por expansión directa, admitida en general) tiene Com(A)\operatorname{Com}(A) con entradas enteras (cada cofactor es un determinante entero); dividir por detA=±1\det A = \pm 1 mantiene los números enteros.

Ejercicio 22.9 ★★★

Calcular el determinante n×nn \times n de la matriz aI+bJaI + bJ (Ejercicio 21.9), es decir con a+ba + b en la diagonal y bb en otro lugar. (Suma todas las columnas a la primera, factoriza, luego claro.) Recuperar la condición de invertibilidad a0a \neq 0,a+nb0a + nb \neq 0.

Solución

Solución de Ejercicio 22.9.

Agregue todas las columnas a la primera: cada entrada de la nueva primera columna es a+nba + nb; Factorízalo, de modo que la primera columna sea todo unos. entonces el operaciones de fila LiLiL1L_i \leftarrow L_i - L_1(i2i \geq 2) claro cada entrada debajo de 11 en la parte superior izquierda y deje aa en la diagonal y 00 en otras partes de esas filas: la matriz es triangular superior con diagonal (1,a,,a)(1, a, \dots, a). Por lo tanto

det(aI+bJ)=(a+nb)an1,\det(aI + bJ) = (a + nb)\, a^{\,n-1} ,

distinto de cero si a0a \neq 0 y a+nb0a + nb \neq 0: la condición de Ejercicio 21.9.

Ejercicio 22.10 ★★★

Sea A,BMn(R)A, B \in \mathcal{M}_n(\R). demostrar que

det(ABBA)=det(A+B)det(AB),\det\begin{pmatrix} A & B\\ B & A \end{pmatrix} = \det(A + B)\,\det(A - B),

por columna de bloque y operaciones de fila (C1C1+C2C_1 \leftarrow C_1 + C_2, luego L2L2L1L_2 \leftarrow L_2 - L_1, en forma de bloque), suponiendo el natural regla triangular de bloque det(MN0P)=detMdetP\det\begin{pmatrix} M & N\\ 0 & P\end{pmatrix} = \det M \det P— probada para 2×22 \times 2 bloques en Ejemplo 22.9.

Solución

Solución de Ejercicio 22.10.

Operaciones de bloque (cada una de las cuales es una composición del escalar nn correspondiente operaciones, permitidas por Teorema 22.3 (1)):

ABBA=C1C1+C2A+BBA+BA=L2L2L1A+BB0AB=det(A+B)det(AB),\begin{vmatrix} A & B\\ B & A\end{vmatrix} \overset{C_1 \leftarrow C_1 + C_2}{=} \begin{vmatrix} A + B & B\\ A + B & A\end{vmatrix} \overset{L_2 \leftarrow L_2 - L_1}{=} \begin{vmatrix} A + B & B\\ 0 & A - B\end{vmatrix} = \det(A+B)\,\det(A-B),

usando la regla del bloque triangular para el último paso.

Ejercicio 22.11 ★★

(Circulante de orden 33) Sea a,b,cCa, b, c \in \C y

Δ=abccabbca.\Delta = \begin{vmatrix} a & b & c\\ c & a & b\\ b & c & a \end{vmatrix}.

Demuestre que Δ=(a+b+c)(a2+b2+c2abbcca)\Delta = (a + b + c)(a^2 + b^2 + c^2 - ab - bc - ca) y factorice completamente C\C usando j=e2iπ/3j = \eu^{2\iu\pi/3}:

Δ=(a+b+c)(a+jb+j2c)(a+j2b+jc).\Delta = (a + b + c)(a + jb + j^2c)(a + j^2b + jc) .

(Start with C1C1+C2+C3C_1 \leftarrow C_1 + C_2 + C_3; for the complex form, note that the column (1,j,j2)T(1, j, j^2)^{\mathsf T} almost behaves like an eigenvector.)

Solución

Solución de Ejercicio 22.11.

C1C1+C2+C3C_1 \leftarrow C_1 + C_2 + C_3 hace que la primera columna sea constante (a+b+c)(a + b + c); factorizarlo, entonces L2L2L1L_2 \leftarrow L_2 - L_1, L3L3L1L_3 \leftarrow L_3 - L_1:

Δ=(a+b+c)1bc0abbc0cbac=(a+b+c)[(ab)(ac)+(bc)2],\Delta = (a+b+c)\begin{vmatrix} 1 & b & c\\ 0 & a - b & b - c\\ 0 & c - b & a - c \end{vmatrix} = (a+b+c)\bigl[(a-b)(a-c) + (b-c)^2\bigr],

y ampliando, (ab)(ac)+(bc)2=a2+b2+c2abbcca(a-b)(a-c) + (b-c)^2 = a^2 + b^2 + c^2 - ab - bc - ca. Sobre C\C, con j3=1j^3 = 1 y 1+j+j2=01 + j + j^2 = 0:

(a+jb+j2c)(a+j2b+jc)=a2+b2+c2+(j+j2)(ab+bc+ca)=a2+b2+c2abbcca,\begin{align*} (a + jb + j^2c)(a + j^2b + jc) &= a^2 + b^2 + c^2 + (j + j^2)(ab + bc + ca)\\ &= a^2 + b^2 + c^2 - ab - bc - ca , \end{align*}

de ahí la factorización completa. (Estructuralmente: la columna (1,j,j2)T(1, j, j^2)^{\mathsf T} satisface M(1,j,j2)T=(a+jb+j2c)(1,j,j2)TM\,(1, j, j^2)^{\mathsf T} = (a + jb + j^2c)(1, j, j^2)^{\mathsf T}, y de manera similar para j2j^2 y 11: los tres factores son los tres "valores propios" del circulante, una historia sistematizada en el volumen Año 2.)

Ejercicio 22.12 ★★★

(Rango y menores) Let AMn,p(K)A \in \mathcal{M}_{n,p}(K). demostrar que rkA\operatorname{rk} A equivale al tamaño más grande rr de un submatriz r×rr \times r invertible de AA(una submatriz mantiene la entradas en los cruces de rr filas elegidas y rr elegidas columnas). (If rkA=r\operatorname{rk} A = r, select rrlibre columns, then rrlibre rows of the resulting n×rn \times r block; conversely, an invertible submatrix forces the corresponding columns of AA to be libre.)

Solución

Solución de Ejercicio 22.12.

Escribe r=rkAr = \operatorname{rk} A.

An invertible r×rr \times r submatrix exists. Elige rrlibre columnas de AA y dejemos que BMn,rB \in \mathcal{M}_{n,r} sea la matriz que forman: rkB=r\operatorname{rk} B = r. Desde el rango de fila es igual al rango de columna (Teorema 21.13), BB tiene rrlibre filas; mantener esas filas produce una submatriz r×rr \times r de AA de rango rr, es decir invertible.

Ninguno más grande lo hace. Sea SS un s×ss \times s invertible submatriz, tomada de las columnas j1,,jsj_1, \dots, j_s y las filas i1,,isi_1, \dots, i_s de AA. Si una combinación kλkCjk=0\sum_k \lambda_k C_{j_k} = 0 de las columnas lleno correspondientes desaparece, luego, al leer solo las filas i1,,isi_1, \dots, i_s se obtiene kλkSk=0\sum_k \lambda_k S_k = 0 en las columnas de SS, por lo que todo λk=0\lambda_k = 0 (SS invertible): las columnas Cj1,,CjsC_{j_1}, \dots, C_{j_s} de AA son libre y srkA=rs \leq \operatorname{rk} A = r.

Por lo tanto rkA\operatorname{rk} A es exactamente el tamaño más grande de un submatriz reversible.

22.4 Problema: la doble alternante de Cauchy

Problema 22.1

Dos determinantes rigen las aplicaciones de este capítulo: la Determinante de Vandermonde, evaluado en Ejercicio 22.5, y el determinante de cauchy det(1ai+bj)\det\bigl(\frac{1}{a_i + b_j}\bigr), evaluado aquí. A su alrededor se acumula este problema. la caja de herramientas alternativa: polinomio trucos de columna, interpolación por Cramer, el matriz de hilabert, el discriminante de una cúbica, y el método de alternancia polinomios. En todo momento, V(x1,,xn)=i<j(xjxi)V(x_1, \dots, x_n) = \prod_{i < j}(x_j - x_i) denota Vandermonde. valor.

Parte I — The Vandermonde toolbox.

  1. Calcule V(1,2,3,4)V(1, 2, 3, 4) y recuerde por qué la interpolación en nn nodos distintos por pares es un sistema Cramer.
  2. (Polinomio alternante) Sea P0,,Pn1P_0, \dots, P_{n-1}mónico con degPk=k\deg P_k = k. probar

    det(Pj1(xi))1i,jn=V(x1,,xn):\det\bigl(P_{j-1}(x_i)\bigr)_{1 \leq i, j \leq n} = V(x_1, \dots, x_n) :

    Las operaciones de columna reemplazan cada columna de potencia por cualquier mónico. escalera, libre.

  3. Aplicar la pregunta 2 al binomio polinomios Bk=X(X1)(Xk+1)k!B_k = \frac{X(X-1)\cdots(X-k+1)}{k!}: demostrar que para números enterosm1<m2<<mnm_1 < m_2 < \dots < m_n,

    V(m1,,mn)0!1!2!(n1)!N:\frac{V(m_1, \dots, m_n)}{0!\,1!\,2!\cdots(n-1)!} \in \N :

    el producto de todas las diferencias por pares de números enteros nn es divisible por el superfactorial 0!1!(n1)!0!\,1!\cdots(n-1)!.

  4. Probar det(xij)1i,jn=x1x2xnV(x1,,xn)\det\bigl(x_i^{\,j}\bigr)_{1 \leq i, j \leq n} = x_1 x_2 \cdots x_n\, V(x_1, \dots, x_n)(poderes ahora a partir de 11).
  5. (Matriz de momentos) Sea S=(pi+j2)1i,jnS = \bigl(p_{i+j-2}\bigr)_{1 \leq i, j \leq n} donde pk=x1k++xnkp_k = x_1^k + \dots + x_n^k. probar que S=WTWS = W^{\mathsf T} W para la matriz W=(xij1)ijW = (x_i^{\,j-1})_{ij}, deduzca

    detS=V(x1,,xn)2,\det S = V(x_1, \dots, x_n)^2 ,

    y concluya: nnreal los números están en pares distintos si y sólo si su matriz de momentos es invertible y detS0\det S \geq 0 siempre.

Parte II — Interpolation revisited. Nodos x1<<xnx_1 < \dots < x_n, valores y1,,yny_1, \dots, y_n.

  1. Escriba las condiciones “P=c0+c1X++cn1Xn1P = c_0 + c_1X + \dots + c_{n-1}X^{n-1} interpola” como sistema lineal en el ckc_k con la matriz WW, y recuperar de detW=V0\det W = V \neq 0 la existencia y unicidad del interpolante (compare las dos pruebas anteriores, Teorema 8.23 y Ejemplo 20.10).
  2. Por la regla de Cramer y expansión de cofactor de la pertinente determinante a lo largo de su última columna, demuestre que el El coeficiente principal del interpolante es

    cn1=i=1nyiji(xixj).c_{n-1} = \sum_{i=1}^{n} \frac{y_i}{\prod_{j \neq i}(x_i - x_j)} .
  3. (Confluente Vandermonde) Computación

    1x1x12012x11x2x22=(x2x1)2,\begin{vmatrix} 1 & x_1 & x_1^2\\ 0 & 1 & 2x_1\\ 1 & x_2 & x_2^2 \end{vmatrix} = (x_2 - x_1)^2 ,

    e interpretar: los datos (P(x1),P(x1),P(x2))\bigl(P(x_1), P'(x_1), P(x_2)\bigr) determinan un PR2[X]P \in \R_2[X] único cuando x1x2x_1 \neq x_2(interpolación de Hermite).

  4. Encuentre el PR2[X]P \in \R_2[X]único con P(0)=1P(0) = 1,P(0)=0P'(0) = 0,P(1)=2P(1) = 2 y compare su respuesta con pregunta 8.

Parte III — The determinante de cauchy. Sean a1,,ana_1, \dots, a_n y b1,,bnb_1, \dots, b_n escalares con ai+bj0a_i + b_j \neq 0 para todos i,ji, j, y

Cn=det(1ai+bj)1i,jn.C_n = \det\Bigl(\frac{1}{a_i + b_j}\Bigr)_{1 \leq i, j \leq n} .
  1. Calcule C2C_2 a mano y póngalo en la forma “productos de diferencias sobre productos de sumas”.
  2. Para n2n \geq 2, realice LiLiLnL_i \leftarrow L_i - L_n(i<ni < n) y factorice filas y columnas para demostrar

    Cn=i<n(anai)j(an+bj)  detM,C_n = \frac{\prod_{i<n}(a_n - a_i)}{\prod_{j}(a_n + b_j)}\;\det M,

    donde MM concuerda con la matriz de Cauchy en las filas i<ni < n y tiene la última fila (1,1,,1)(1, 1, \dots, 1).

  3. Realice CjCjCnC_j \leftarrow C_j - C_n(j<nj < n) en MM, factorizar nuevamente y concluir por inducción La doble alternante de Cauchy:

    Cn=1i<jn(ajai)(bjbi)i,j(ai+bj).C_n = \frac{\prod_{1 \leq i < j \leq n}(a_j - a_i)(b_j - b_i)}{\prod_{i, j}(a_i + b_j)} .
  4. Deducir el criterio de invertibilidad (el aia_i por pares distinto y el bjb_j distinto por pares). Para el matriz de hilabert Hn=(1i+j1)H_n = \bigl(\frac{1}{i + j - 1}\bigr): calcula detH2\det H_2 y detH3\det H_3 a partir del fórmula y verifique que H21H_2^{-1} tenga entradas enteras.
  5. Demuestre que para bjb_j distintos por pares y cualquier mano derecha lado, el sistema jcjai+bj=yi\sum_j \frac{c_j}{a_i + b_j} = y_i (i=1,,ni = 1, \dots, n) tiene una solución única y se relaciona esto a la existencia y unicidad de la fracción parcial descomposiciones con polos simples (Teorema 9.5).

Parte IV — The discriminant of a cubic. Deja λ1,λ2,λ3\lambda_1, \lambda_2, \lambda_3 sean las raíces (en C\C) de X3+pX+qX^3 + pX + q y pk=λ1k+λ2k+λ3kp_k = \lambda_1^k + \lambda_2^k + \lambda_3^k.

  1. Usando λ3=pλq\lambda^3 = -p\lambda - q en cada raíz y Vieta (p1=0p_1 = 0), calcule p2=2pp_2 = -2p,p3=3qp_3 = -3q y p4=2p2p_4 = 2p^2.
  2. Con la pregunta 5 (sobre C\C, manteniendo detS=V2\det S = V^2), calcular

    disc=V(λ1,λ2,λ3)2=302p02p3q2p3q2p2=4p327q2.\operatorname{disc} = V(\lambda_1, \lambda_2, \lambda_3)^2 = \begin{vmatrix} 3 & 0 & -2p\\ 0 & -2p & -3q\\ -2p & -3q & 2p^2 \end{vmatrix} = -4p^3 - 27q^2 .
  3. Deducir: X3+pX+qX^3 + pX + q tiene una raíz repetida si y sólo si 4p3+27q2=04p^3 + 27q^2 = 0; consulte X33X+2=(X1)2(X+2)X^3 - 3X + 2 = (X - 1)^2(X + 2).
  4. Supongamos que p,qp, q es real. Demostrar que la cúbica tiene tres. raíces reales distintas si y sólo si disc>0\operatorname{disc} > 0, y una raíz conjugada real más dos no reales si y sólo si disc<0\operatorname{disc} < 0. (If λ3=λ2λ2\lambda_3 = \conj{\lambda_2} \neq \lambda_2 and λ1R\lambda_1 \in \R, show VV is purely imaginary.)

Part V — Dividends, and the alternating method.

  1. Para 0<a1<a2<<an0 < a_1 < a_2 < \dots < a_n, mostrar det(1ai+aj)>0\det\bigl(\frac{1}{a_i + a_j}\bigr) > 0.
  2. Calcular det((mij1))1i,j3\det\bigl(\binom{m_i}{j-1}\bigr)_{1 \leq i, j \leq 3} para (m1,m2,m3)=(2,4,7)(m_1, m_2, m_3) = (2, 4, 7), primero por preguntas 2–3, luego por expansión directa.
  3. Sea λ1,,λn\lambda_1, \dots, \lambda_n distinto por pares y distinto de cero. Usando una matriz de Vandermonde invertible, demostrar nuevamente que las sucesiones geométricas ((λik)k0)1in\bigl((\lambda_i^{\,k})_{k \geq 0}\bigr)_{1 \leq i \leq n} forma un familia libre del espacio de secuencias.
  4. Calcular det(1i+j)1i,j3\det\bigl(\frac{1}{i + j}\bigr)_{1 \leq i, j \leq 3} a partir de la doble alternante.
  5. (Alternando polinomios) Llame a polinomio FF en x1,,xnx_1, \dots, x_nalterno al intercambiar dos variable cambia de signo. Demuestre que un FF alterno desaparece cada vez que xi=xjx_i = x_j(iji \neq j), y deduce — una variable a la vez, según el teorema del factor — que FF es divisible por i<j(xjxi)\prod_{i<j}(x_j - x_i).
  6. Utilice la pregunta 23 para refutar la fórmula de Vandermonde con sin inducción: el determinante det(xij1)\det(x_i^{\,j-1}) es un alternando polinomio de grado total (n2)\binom n2, por tanto, un múltiplo constante de i<j(xjxi)\prod_{i<j}(x_j - x_i); Identifica la constante comparando un monomio.
  7. Síntesis, en cuatro frases: ¿qué propiedad única de el determinante (qué axioma) genera todos los factorizaciones de este problema; ¿Por qué la matriz de momentos? La identidad de la pregunta 5 genera un enunciado sobre complejo distinción en un computable real prueba de signos; cuales dos matrices clásicas eran evaluados completamente aquí y qué problemas lineal gobernar; y cómo el método alterno de preguntas 23–24 explica, de un plumazo, por qué sigue apareciendo i<j(xjxi)\prod_{i<j}(x_j - x_i). Nombra el teorema de la Parte III.
Solución

Solución de Problema 22.1.

1. V(1,2,3,4)=(21)(31)(41)(32)(42)(43)=123121=12V(1,2,3,4) = (2-1)(3-1)(4-1)(3-2)(4-2)(4-3) = 1 \cdot 2\cdot 3\cdot 1\cdot 2\cdot 1 = 12. Interpolación en distintos nodos solicita los coeficientes de PP resolviendo Wc=yW c = y con W=(xij1)W = (x_i^{\,j-1}) y detW=V0\det W = V \neq 0: un Cramer sistema.

2. Trabaje en las columnas de izquierda a derecha. C1C_1 es la columna constante P0(xi)=1P_0(x_i) = 1(P0P_0mónico de grado 00). Supongamos que las columnas 1,,j11, \dots, j-1 ya se han reducido. a los poderes puros 1,xi,,xij21, x_i, \dots, x_i^{\,j-2}. Desde Pj1=Xj1+k<j1αkXkP_{j-1} = X^{j-1} + \sum_{k < j-1}\alpha_k X^k, restando de CjC_j la combinación kαk(column of xik)\sum_k \alpha_k\,(\text{column of } x_i^k) — una operación que no cambia el determinante — deja la columna de pura potencia xij1x_i^{\,j-1}. Después de la última columna, La matriz es la matriz de Vandermonde: det=V(x1,,xn)\det = V(x_1, \dots, x_n).

3. Los polinomios (j1)!Bj1(j-1)!\,B_{j-1} son mónico de grado j1j - 1, por lo que la pregunta 2 da

det(Bj1(mi))=V(m1,,mn)0!1!(n1)!.\det\bigl(B_{j-1}(m_i)\bigr) = \frac{V(m_1, \dots, m_n)}{0!\,1!\cdots(n-1)!} .

El lado izquierdo es el determinante de una matriz con entero entradas (BkB_k tiene un valor entero en Z\Z: preguntas 16–17 de el problema del fin de semana Problema 18.1), por lo tanto, un número entero; y es positivo desde V(m1,,mn)>0V(m_1, \dots, m_n) > 0 para m1<<mnm_1 < \dots < m_n. Entonces el superfactorial divide es el producto de todos diferencias por pares.

4. Factoriza xix_i de cada fila ii: det(xij)j=1..n=x1xndet(xij1)=x1xnV\det(x_i^{\,j})_{j = 1..n} = x_1\cdots x_n\, \det(x_i^{\,j-1}) = x_1\cdots x_n\,V.

5. (WTW)ij=kxki1xkj1=pi+j2(W^{\mathsf T}W)_{ij} = \sum_k x_k^{\,i-1} x_k^{\,j-1} = p_{i+j-2}:S=WTWS = W^{\mathsf T}W. Por lo tanto detS=det(WT)detW=V2\det S = \det(W^{\mathsf T})\det W = V^2 (Teorema 22.3 (2),(4)). De verdad xix_i:detS=V20\det S = V^2 \geq 0 y SS son invertibles si V0V \neq 0 y si xix_i son distinto por pares — una prueba de signo definido computable a partir de la sumas de potencia solas.

6. Las condiciones de interpolación kckxik=yi\sum_{k} c_k\,x_i^{\,k} = y_i forman el sistema Wc=yWc = y;detW=V0\det W = V \neq 0 da existencia y unicidad a la vez. Este es el tercero prueba en el libro: fórmula explícita en Teorema 8.23, argumento del núcleo en Ejemplo 20.10, aquí Cramer.

7. Cramer: cn1=detW/detWc_{n-1} = \det W'/\det W donde WW' es WW con su última columna reemplazada por yy. Expandiendo detW\det W' a lo largo esa columna:

detW=i=1n(1)i+nyiV(x1,,xi^,,xn).\det W' = \sum_{i=1}^n (-1)^{i+n} y_i\,V(x_1, \dots, \widehat{x_i}, \dots, x_n) .

Ahora V=V(i)j<i(xixj)j>i(xjxi)V = V(\setminus i)\cdot\prod_{j<i}(x_i - x_j)\prod_{j>i} (x_j - x_i), y convirtiendo los costos del segundo producto. (1)ni(-1)^{n-i}:

(1)i+nV(i)V=(1)i+n(1)niji(xixj)=1ji(xixj),(-1)^{i+n}\,\frac{V(\setminus i)}{V} = \frac{(-1)^{i+n}(-1)^{n-i}}{\prod_{j\neq i}(x_i - x_j)} = \frac{1}{\prod_{j\neq i}(x_i - x_j)} ,

de donde cn1=iyi/ji(xixj)c_{n-1} = \sum_i y_i/\prod_{j \neq i}(x_i - x_j) — la fórmula de diferencia dividida nuevamente.

8. L3L3L1L_3 \leftarrow L_3 - L_1 proporciona las filas (1,x1,x12)(1, x_1, x_1^2),(0,1,2x1)(0, 1, 2x_1),(0, x2x1, (x2x1)(x2+x1))(0,\ x_2 - x_1,\ (x_2-x_1)(x_2+x_1)); expandiendo a lo largo de la primera columna y factorizando (x2x1)(x_2 - x_1):

(x2x1)12x11x2+x1=(x2x1)(x2x1)=(x2x1)2.(x_2 - x_1)\begin{vmatrix} 1 & 2x_1\\ 1 & x_2 + x_1 \end{vmatrix} = (x_2 - x_1)(x_2 - x_1) = (x_2 - x_1)^2 .

No nulo para x1x2x_1 \neq x_2: el sistema lineal expresa P(x1)=uP(x_1) = u,P(x1)=vP'(x_1) = v,P(x2)=wP(x_2) = w sobre los coeficientes de PR2[X]P \in \R_2[X] es Cramer — interpolación de Hermite con un El nudo duplicado está bien planteado.

9. P=a+bX+cX2P = a + bX + cX^2 con a=P(0)=1a = P(0) = 1,b=P(0)=0b = P'(0) = 0,a+b+c=P(1)=2a + b + c = P(1) = 2:c=1c = 1, entonces P=1+X2P = 1 + X^2, único. Consistencia: aquí x1=0x_1 = 0,x2=1x_2 = 1 y el determinante de La pregunta 8 es (10)2=10(1 - 0)^2 = 1 \neq 0.

10. Cálculo directo:

C2=1(a1+b1)(a2+b2)1(a1+b2)(a2+b1)=(a1+b2)(a2+b1)(a1+b1)(a2+b2)i,j(ai+bj),C_2 = \frac{1}{(a_1+b_1)(a_2+b_2)} - \frac{1}{(a_1+b_2)(a_2+b_1)} = \frac{(a_1+b_2)(a_2+b_1) - (a_1+b_1)(a_2+b_2)} {\prod_{i,j}(a_i+b_j)} ,

y el numerador se expande a a1b1+a2b2a1b2a2b1=(a2a1)(b2b1)a_1b_1 + a_2b_2 - a_1b_2 - a_2b_1 = (a_2 - a_1)(b_2 - b_1): diferencias sobre sumas.

11. Para i<ni < n, la nueva entrada de la fila ii es

1ai+bj1an+bj=anai(ai+bj)(an+bj).\frac{1}{a_i + b_j} - \frac{1}{a_n + b_j} = \frac{a_n - a_i}{(a_i + b_j)(a_n + b_j)} .

Factorice (anai)(a_n - a_i) de cada fila i<ni < n, luego 1an+bj\frac1{a_n + b_j} de cada columna jj: lo que queda tiene entradas 1ai+bj\frac1{a_i + b_j} en filas i<ni < n y constante 11 en fila nn— la matriz MM, con el prefactor anunciado.

12. En MM, para j<nj < n la operación CjCjCnC_j \leftarrow C_j - C_nconvierte la fila nn en (0,,0,1)(0, \dots, 0, 1) y, en la fila i<ni < n,

1ai+bj1ai+bn=bnbj(ai+bj)(ai+bn).\frac{1}{a_i + b_j} - \frac{1}{a_i + b_n} = \frac{b_n - b_j}{(a_i + b_j)(a_i + b_n)} .

Factorice (bnbj)(b_n - b_j) de cada columna j<nj < n y 1ai+bn\frac1{a_i + b_n} de cada fila i<ni < n, luego expanda a lo largo de la última fila (signo (1)n+n=+1(-1)^{n+n} = +1): el determinante restante es Cn1C_{n-1}. Recopilando los factores de las preguntas 11–12:

Cn=i<n(anai)j<n(bnbj)j(an+bj)i<n(ai+bn)  Cn1,C_n = \frac{\prod_{i<n}(a_n - a_i)\,\prod_{j<n}(b_n - b_j)} {\prod_{j}(a_n + b_j)\,\prod_{i<n}(a_i + b_n)}\;C_{n-1},

y la inducción (base C1=1a1+b1C_1 = \frac1{a_1+b_1}) ensambla exactamente la doble alternante de Cauchy: los factores (ajai)(bjbi)(a_j - a_i)(b_j - b_i) para todos los pares, sobre todas las sumas (ai+bj)(a_i + b_j).

13. La fórmula desaparece si solo si aj=aia_j = a_i o bj=bib_j = b_i: la matriz de Cauchy es invertible si ambas familias son distintos por pares. Hilbert: ai=ia_i = i,bj=j1b_j = j - 1. Para n=2n = 2: numerador (21)(10)=1(2-1)(1-0) = 1, denominador 1223=121\cdot2\cdot2\cdot3 = 12:detH2=112\det H_2 = \frac1{12}. Para n=3n = 3: numerador [(1)(2)(1)]2=4\bigl[(1)(2)(1)\bigr]^2 = 4, denominador (123)(234)(345)=62460=8640(1\cdot2\cdot3) (2\cdot3\cdot4)(3\cdot4\cdot5) = 6\cdot24\cdot60 = 8640: detH3=48640=12160\det H_3 = \frac{4}{8640} = \frac1{2160}. Inverso para n=2n = 2:

H21=12(1312121)=(46612),H_2^{-1} = 12\begin{pmatrix} \frac13 & -\frac12\\[2pt] -\frac12 & 1\end{pmatrix} = \begin{pmatrix} 4 & -6\\ -6 & 12 \end{pmatrix},

todos los números enteros (un fenómeno cierto para cada nn).

14. La matriz del sistema es la matriz de Cauchy, invertible por la pregunta 13 cuando el bjb_j(y el aia_i) son distinta por pares: solución única. Interpretación: un racional La función R=jcjX+bjR = \sum_j \frac{c_j}{X + b_j} con polos simple es determinado por nn de sus valores R(a1),,R(an)R(a_1), \dots, R(a_n), y por el contrario, cualquier hoja de datos de este tipo se realiza exactamente una vez — la contraparte muestral de la existencia de fracción parcial y teorema de unicidad (Teorema 9.5).

15. Vieta para X3+pX+qX^3 + pX + q:λ1+λ2+λ3=0\lambda_1 + \lambda_2 + \lambda_3 = 0,i<jλiλj=p\sum_{i<j}\lambda_i\lambda_j = p, entonces p1=0p_1 = 0 y p2=p122p=2pp_2 = p_1^2 - 2p = -2p. Cada raíz satisface λ3=pλq\lambda^3 = -p\lambda - q; sumando:p3=pp13q=3qp_3 = -p\,p_1 - 3q = -3q. multiplicando por λ\lambda y sumando:p4=pp2qp1=2p2p_4 = -p\,p_2 - q\,p_1 = 2p^2.

16. Por pregunta 5 (la identidad S=WTWS = W^{\mathsf T}W y detS=V2\det S = V^2 son válidos sobre C\C),

V2=302p02p3q2p3q2p2=3(4p39q2)+(2p)(04p2)=4p327q2,V^2 = \begin{vmatrix} 3 & 0 & -2p\\ 0 & -2p & -3q\\ -2p & -3q & 2p^2 \end{vmatrix} = 3\bigl(-4p^3 - 9q^2\bigr) + (-2p)\bigl(0 - 4p^2\bigr) = -4p^3 - 27q^2 ,

expandiéndose a lo largo de la primera fila.

17. Una raíz repetida significa dos λi\lambda_i iguales, es decir V=0V = 0, es decir disc=4p327q2=0\operatorname{disc} = -4p^3 - 27q^2 = 0. Para X33X+2X^3 - 3X + 2:4(3)3+274=108+108=04(-3)^3 + 27\cdot4 = -108 + 108 = 0, coincidente la doble raíz 11 de (X1)2(X+2)(X-1)^2(X+2).

18. Las raíces no reales de un cúbico real vienen en conjugado pares, por lo que ocurren exactamente dos casos cuando disc0\operatorname{disc} \neq 0. Tres raíces reales distintas:VV es real y distinta de cero, por lo que disc=V2>0\operatorname{disc} = V^2 > 0. Una raíz real λ1\lambda_1 y λ3=λ2R\lambda_3 = \conj{\lambda_2} \notin \R: entonces

(λ2λ1)(λ3λ1)=λ2λ12>0,λ3λ2=2iImλ20,(\lambda_2 - \lambda_1)(\lambda_3 - \lambda_1) = \abs{\lambda_2 - \lambda_1}^2 > 0, \qquad \lambda_3 - \lambda_2 = -2\iu\,\operatorname{Im}\lambda_2 \neq 0,

entonces VV es un número puramente imaginario distinto de cero y disc=V2<0\operatorname{disc} = V^2 < 0. Los dos signos caracterizan el dos casos.

19. Tome bi=aib_i = a_i en la doble alternante: el el numerador es i<j(ajai)2>0\prod_{i<j}(a_j - a_i)^2 > 0 y el denominador i,j(ai+aj)>0\prod_{i,j}(a_i + a_j) > 0(todas las entradas positivas): el determinante es positivo. (En lenguaje posterior: el núcleo 1x+y\frac1{x+y} es positivo definido.)

20. Por las preguntas 2–3, el determinante es igual V(2,4,7)/(0!1!2!)=(42)(72)(74)2=302=15V(2,4,7)/(0!\,1!\,2!) = \frac{(4-2)(7-2)(7-4)}{2} = \frac{30}{2} = 15. Directamente, la matriz es

(1211461721),det=(8442)2(216)+(74)=4230+3=15.\begin{pmatrix} 1 & 2 & 1\\ 1 & 4 & 6\\ 1 & 7 & 21 \end{pmatrix}, \qquad \det = (84 - 42) - 2(21 - 6) + (7 - 4) = 42 - 30 + 3 = 15 .

21. Supongamos ici(λik)k=0\sum_i c_i\,(\lambda_i^{\,k})_{k} = 0 como una secuencia. Al leer k=0,1,,n1k = 0, 1, \dots, n-1 se obtiene WTc=0W^{\mathsf T}c = 0 con W=(λij1)W = (\lambda_i^{\,j-1}) invertible (det=V0\det = V \neq 0, distinto λi\lambda_i):c=0c = 0. Las secuencias geométricas son libre.

22. a=b=(1,2,3)a = b = (1, 2, 3): numerador [(21)(31)(32)]2=4\bigl[(2-1)(3-1) (3-2)\bigr]^2 = 4; denominador i,j(i+j)=(234)(345)(456)=2460120=172800\prod_{i,j}(i + j) = (2\cdot3\cdot4)(3\cdot4\cdot5)(4\cdot5\cdot6) = 24\cdot60\cdot120 = 172800. Por lo tanto det(1i+j)=4172800=143200\det\bigl(\frac1{i+j}\bigr) = \frac{4}{172800} = \frac1{43200}.

23. Si xi=xjx_i = x_j, el intercambio de las dos variables soluciona el punto pero debe cambiar el signo de FF:F=FF = -F, por lo F=0F = 0 allí. Divisibilidad: ver FF como polinomio en el single variable xnx_n con coeficientes en las demás variables; eso desaparece en los “valores” n1n - 1 x1,,xn1x_1, \dots, x_{n-1}, por lo que La factorización repetida (Teorema 8.7) da F=i<n(xnxi)GF = \prod_{i<n}(x_n - x_i)\cdot G con GGpolinomio. el el prefactor es invariante bajo swaps de dos índices i,j<ni, j < n, por lo que GG se alterna en x1,,xn1x_1, \dots, x_{n-1} y la inducción completa: i<j(xjxi)\prod_{i<j}(x_j - x_i)divideFF.

24. D=det(xij1)D = \det(x_i^{\,j-1}) es un polinomio en el xix_i; intercambiar dos variables intercambia dos filas, por lo que DD es alternando, y por la pregunta 23, D=ci<j(xjxi)D = c\,\prod_{i<j}(x_j - x_i) para algunos polinomio cc. Grados totales:DD tiene grado 0+1++(n1)=(n2)\leq 0 + 1 + \dots + (n-1) = \binom n2, el producto tiene grado exactamente (n2)\binom n2:cc es una constante. El monomio x2x32xnn1x_2\,x_3^2\cdots x_n^{\,n-1} tiene coeficiente 11 en DD(producto diagonal) y 11 en el producto (elija la variable de índice mayor en cada factor): c=1c = 1, y la fórmula de Vandermonde cae sin inducción.

25. (i) Alternancia — el axioma “dos columnas iguales matar el determinante” — es el motor: produjo cada factor (xjxi)(x_j - x_i),(ajai)(a_j - a_i),(bjbi)(b_j - b_i) en el problema. (ii) La identidad detS=V2\det S = V^2 reemplaza la raíces individualmente complejas e inalcanzables por sus sumas de potencias, que son reales polinomios en los coeficientes, por lo que la distinción se convierte en el signo de un número real computable. (iii) El Determinante de Vandermonde gobierna la interpolación polinomio, el determinante de cauchy gobierna fracciones parciales y muestreadas funciones racionales (con el matriz de hilabert como su más famoso caso especial). (iv) Cualquier polinomio alternativo es divisible por i<j(xjxi)\prod_{i<j}(x_j - x_i), y un conteo de grados luego fija dicho polinomio hasta una constante — razón por la cual este producto sigue reapareciendo dondequiera que un determinante desaparece coincidencias. El teorema de la Parte III es El doble de Cauchy alternante.