Álgebra lineal

Vectores, matrices y sistemas de ecuaciones: eliminación, LU, inversas, modelos de insumo–producto, los cuatro subespacios fundamentales, proyecciones y mínimos cuadrados, determinantes, valores propios, cadenas de Markov y la SVD, con seis prácticas. Gratis, a tu ritmo y vinculado a GPS, Landsat, satélites del nivel del mar, modelos meteorológicos, reciclaje y cartografía mineral.

NASA/Kim Shiflett (KSC-20210927-PH-KLS04_0013, Sept. 27, 2021)NASA/JPL-Caltech (PIA22431)CNES (PIA25623)NASA Blue Marble

Cursos › MATH 230

Descargar la guía de estudio en PDF (inglés)

Sobre este curso

Este es el esquema de un curso gratuito de autoestudio escrito por Saving Planets. No es un curso universitario, no otorga créditos académicos y ninguna institución participa.

Su alcance sigue un curso de álgebra lineal de un semestre de segundo año. Se basa en los temas y la profundidad de MIT 18.06 Linear Algebra impartido por Gilbert Strang (gratis en MIT OpenCourseWare) y en los cursos de álgebra lineal de segundo año de Georgia Tech y Purdue, usando el texto gratuito de Georgia Tech Interactive Linear Algebra de Dan Margalit y Joseph Rabinoff y el libro gratuito Linear Algebra de Jim Hefferon, que enlazamos donde ayudan. Ninguno de esos autores o instituciones respalda esta página.

Lo que añadimos: cada unidad se vincula a trabajo real sobre el planeta con una fuente pública citada (GPS, Landsat 9, satélites del nivel del mar de la NASA, supercomputadoras meteorológicas de la NOAA, datos de insumo–producto y reciclaje de la EPA, estudios de edificios del USGS, espectrómetro mineral EMIT de la NASA), cada número de la página se recalcula con un script de verificación antes de publicarse, y cada práctica se hace en una hoja de cálculo gratuita o en papel.

Síguelo en orden. Cada unidad tiene objetivos, ideas clave, ecuaciones, ejemplos paso a paso, práctica con respuestas completas, errores comunes, tarjetas de estudio y un cuestionario. Tu progreso se guarda solo en este navegador.

Cómo ayuda este curso a salvar el planeta

El álgebra lineal es el lenguaje de los datos sobre el planeta. Un píxel de satélite es un vector, una posición GPS es la solución de un sistema lineal, un pronóstico del tiempo son millones de ecuaciones resueltas una y otra vez, y la contaminación de una economía puede rastrearse con una sola matriz.

Las mismas pocas ideas, eliminación, proyección, valores propios y la SVD, convierten mediciones en bruto en mapas de bosques, tendencias del nivel del mar y polvo mineral. Cada unidad termina con una breve sección con fuentes. Los valores que no provienen de una fuente se marcan como valores de ejemplo.

Volver arriba

Antes de empezar

Requisitos previos

  • Álgebra de secundaria: resolver dos ecuaciones con dos incógnitas y trabajar con funciones.
  • Cálculo I ayuda con límites y derivadas en algunos ejemplos, pero no se necesita cálculo en la mayor parte del curso.
  • Soltura con una hoja de cálculo; Introducción a la programación ayuda en las versiones opcionales de las prácticas en Python.
  • Cálculo III y Ecuaciones diferenciales combinan bien con este curso pero no son obligatorios.

Cómo estudiar

  • Haz primero cada cálculo a mano con ejemplos pequeños (2 × 2 y 3 × 3) y luego compruébalo con una hoja de cálculo o calculadora.
  • Piensa en Ax de dos formas: como productos escalares de las filas con x, y como una combinación de las columnas de A.
  • Comprueba cada respuesta: sustituye x en Ax = b, multiplica A por su inversa o verifica Av = λv.
  • Controla los tamaños: una matriz m × n por una n × p da m × p, y tamaños que no coinciden indican un error.
  • Haz dibujos en 2D y 3D: rectas, planos, proyecciones y estiramientos hacen las reglas fáciles de recordar.
  • Haz cada ejercicio en papel antes de abrir la respuesta, y repite los cuestionarios unos días después.

Material de laboratorio (para los laboratorios presenciales)

  • Un ordenador con una hoja de cálculo gratuita (sus funciones MMULT, MINVERSE, MDETERM y LINEST hacen el trabajo matricial).
  • Papel cuadriculado, una regla, un transportador y un lápiz.
  • Opcional: Python gratuito con NumPy para las versiones más grandes de las prácticas.
  • Un cordel y chinchetas o cinta adhesiva para la práctica de trilateración, y un mapa impreso o una hoja de papel.
  • Una foto en escala de grises o una pequeña cuadrícula de números para la práctica de compresión de imágenes.

Volver arriba

Calendario sugerido de 15 semanas

Unas 6 a 9 horas por semana: lectura, problemas de práctica y, en las semanas de laboratorio, el laboratorio.

SemanaUnidadTemasLectura
1Vectores: combinaciones lineales, longitudes y productos escalaresVectores, combinaciones lineales, espacio generado, longitudMIT 18.06 L1; ILA 2.1–2.2
2Vectores: combinaciones lineales, longitudes y productos escalaresProductos escalares, ángulos, ortogonalidad; Práctica A; inicio de la eliminaciónILA 6.1
3Sistemas de ecuaciones lineales y eliminación gaussianaEliminación, formas escalonadas, conjuntos solución; Práctica BMIT 18.06 L2; ILA 1.1–1.3
4Álgebra de matrices: multiplicación, inversas y factorización LUMultiplicación de matrices, inversas, A = LU; Práctica CMIT 18.06 L3–4; ILA 3.4–3.5
5Modelos lineales: economía de insumo–producto (Leontief), redes y balanceoModelos de insumo–producto, redes, balanceoHefferon Topics; EPA USEEIO
6Espacios vectoriales, subespacios, independencia, base, rango y los cuatro subespacios fundamentalesSubespacios, espacio columna, espacio nulo, Ax = 0MIT 18.06 L5–7; ILA 2.6
7Espacios vectoriales, subespacios, independencia, base, rango y los cuatro subespacios fundamentalesIndependencia, base, dimensión, rango, cuatro subespacios; práctica de examen parcialMIT 18.06 L8–10; ILA 2.5, 2.7, 2.9
8Ortogonalidad, proyecciones, mínimos cuadrados y Gram–SchmidtOrtogonalidad y proyeccionesMIT 18.06 L14–15; ILA 6.3
9Ortogonalidad, proyecciones, mínimos cuadrados y Gram–SchmidtMínimos cuadrados, Gram–Schmidt, A = QR; Práctica DMIT 18.06 L16–17; ILA 6.4–6.5
10Determinantes: propiedades, cofactores, regla de Cramer y volumenDeterminantes, cofactores, regla de Cramer, volumenMIT 18.06 L18–20; ILA 4.1–4.3
11Valores propios y vectores propiosValores propios, vectores propios, método de la potencia, modosMIT 18.06 L21; ILA 5.1–5.2, 5.5
12Diagonalización, potencias de matrices, cadenas de Markov y sistemas dinámicos linealesDiagonalización, cadenas de Markov, sistemas dinámicos; Práctica EMIT 18.06 L22–25; ILA 5.4, 5.6
13La descomposición en valores singulares, el análisis de componentes principales y las matrices definidas positivasMatrices definidas positivas y la SVDMIT 18.06 L27–29
14La descomposición en valores singulares, el análisis de componentes principales y las matrices definidas positivasAproximación de rango bajo, PCA, número de condición; Práctica FMIT 18.06 L29; sitio del libro de Strang
15RepasoRepaso y serie de práctica finalTodas las unidades

Volver arriba

Tu progreso

Marca las unidades y los laboratorios al terminarlos y haz los cuestionarios. El progreso y las puntuaciones se guardan solo en este navegador y nunca se nos envían.

Activa JavaScript para seguir tu progreso.

Volver arriba

Unidad 1 · Semanas 1–2

Vectores: combinaciones lineales, longitudes y productos escalares

Vectores en dos, tres y n dimensiones, suma y escalado, combinaciones lineales y espacio generado, longitud, vectores unitarios, producto escalar, ángulo entre vectores, ortogonalidad y las desigualdades de Cauchy–Schwarz y triangular.

Objetivos de aprendizaje

Al terminar esta unidad podrás:

  • Sumar, escalar y combinar vectores a mano y con dibujos.
  • Calcular longitudes, vectores unitarios y productos escalares.
  • Hallar el ángulo entre dos vectores y comprobar la ortogonalidad.
  • Describir el espacio generado por uno, dos o tres vectores.

Conceptos clave

  • Un vector en ℝⁿ es una lista ordenada de n números. En 2D y 3D es una flecha; en dimensiones mayores es una lista, como los 11 valores de banda de un píxel de satélite.
  • Una combinación lineal c₁v₁ + c₂v₂ + … escala vectores y los suma. El espacio generado es el conjunto de todas las combinaciones lineales: una recta, un plano o todo el espacio.
  • Longitud ‖v‖ = √(v·v). Dividir por la longitud da un vector unitario u = v/‖v‖ de longitud 1.
  • Producto escalar v·w = v₁w₁ + v₂w₂ + … = ‖v‖‖w‖cos θ. Si v·w = 0 los vectores son ortogonales (perpendiculares).
  • Cauchy–Schwarz: |v·w| ≤ ‖v‖‖w‖. Desigualdad triangular: ‖v + w‖ ≤ ‖v‖ + ‖w‖.

Ecuaciones clave

∥v∥=v12+v22+⋯+vn2\|v\|=\sqrt{v_1^2+v_2^2+\cdots+v_n^2}
Longitud
v⋅w=v1w1+⋯+vnwnv\cdot w=v_1w_1+\cdots+v_nw_n
Producto escalar
cos⁡θ=v⋅w∥v∥ ∥w∥\cos\theta=\frac{v\cdot w}{\|v\|\,\|w\|}
Ángulo entre vectores
∣v⋅w∣≤∥v∥ ∥w∥|v\cdot w|\le\|v\|\,\|w\|
Desigualdad de Cauchy–Schwarz

Ejemplos resueltos paso a paso

Ejemplo 1.1 Halla la longitud de v = (3, 4) y el vector unitario en su dirección.

Ver los pasos
  1. ‖v‖ = √(3² + 4²) = √25 = 5.
  2. u = v/5.

Respuesta: ‖v‖ = 5, u = (0.6, 0.8).

Ejemplo 1.2 Halla el ángulo entre v = (3, 4) y w = (1, 2).

Ver los pasos
  1. v·w = 3 + 8 = 11.
  2. ‖v‖ = 5, ‖w‖ = √5 = 2.236.
  3. cos θ = 11/(5 × 2.236) = 0.9839.

Respuesta: θ = 10.30°.

Ejemplo 1.3 Escribe 2(1, 0, 1) + 3(0, 1, 1) como un solo vector. ¿Está (2, 3, 4) en el espacio generado por estos dos vectores?

Ver los pasos
  1. 2(1, 0, 1) + 3(0, 1, 1) = (2, 3, 5).
  2. Toda combinación c₁(1, 0, 1) + c₂(0, 1, 1) = (c₁, c₂, c₁ + c₂), así que la tercera componente debe ser la suma de las dos primeras.

Respuesta: (2, 3, 5); (2, 3, 4) no está en el espacio generado porque 2 + 3 ≠ 4.

Problemas de práctica

Resuelve cada uno primero en papel y luego comprueba.

  1. Halla ‖(1, 2, 2)‖.

    Ver la respuesta

    √(1 + 4 + 4) = 3.

  2. Halla (1, 2, 3)·(4, −5, 6).

    Ver la respuesta

    4 − 10 + 18 = 12.

  3. ¿Son ortogonales (2, −1) y (1, 2)?

    Ver la respuesta

    Yes: 2 − 2 = 0.

  4. Halla el vector unitario en la dirección de (6, 8).

    Ver la respuesta

    ‖(6, 8)‖ = 10, así que (0.6, 0.8).

Ideas erróneas comunes

  • Mito: El producto escalar es un vector.

    En realidad: Es un solo número (un escalar).

  • Mito: ‖v + w‖ = ‖v‖ + ‖w‖ siempre.

    En realidad: Solo cuando v y w apuntan en el mismo sentido; si no, es menor.

  • Mito: Los vectores solo tienen sentido en 2D o 3D.

    En realidad: Los vectores de datos pueden tener cientos de componentes; se aplican las mismas reglas.

Tarjetas de estudio

  • Combinación lineal

    c₁v₁ + c₂v₂ + …

  • Espacio generado

    Todas las combinaciones lineales de un conjunto de vectores.

  • Longitud

    ‖v‖ = √(v·v).

  • Producto escalar

    v·w = Σ vᵢwᵢ = ‖v‖‖w‖cos θ.

  • Ortogonal

    v·w = 0.

  • Vector unitario

    v/‖v‖.

Cuestionario de autoevaluación

1. ‖(3, 4)‖ es igual a…

Ver la respuesta

5. √(9 + 16) = 5.

2. v·w = 0 significa que los vectores son…

Ver la respuesta

ortogonales. cos θ = 0, así que θ = 90°.

3. El espacio generado por dos vectores no paralelos en ℝ³ es…

Ver la respuesta

un plano. Dos direcciones independientes llenan un plano.

4. Duplicar v cambia ‖v‖ en un factor de…

Ver la respuesta

2. ‖cv‖ = |c|‖v‖.

Cómo ayuda esto a salvar el planeta

Cada píxel de Landsat 9 es un vector. Según el USGS, el Operational Land Imager 2 de Landsat 9 registra 9 bandas espectrales a 30 m (15 m para la banda pancromática) y su Thermal Infrared Sensor 2 registra 2 bandas térmicas a 100 m, así que cada punto del suelo se convierte en una lista de 11 números. Los datos son de 14 bits, 2¹⁴ = 16384 tonos por banda, y Landsat 9 adquiere unas 750 escenas al día, todas gratuitas.

La salud de bosques y cultivos se lee en la banda roja (banda 4, 0.64–0.67 µm) y la banda del infrarrojo cercano (banda 5, 0.85–0.88 µm). El índice de vegetación NDVI = (NIR − red)/(NIR + red) compara dos componentes del vector píxel. Con reflectancias de ejemplo de 0.40 (NIR) y 0.08 (rojo), NDVI = 0.32/0.48 = 0.667, típico de vegetación sana; el suelo desnudo da valores cercanos a cero.

En Saving Planets: Observación de la Tierra · Deforestación

Fuentes: USGS: Landsat 9 · Interactive Linear Algebra, 2.1 Vectors

Recursos de aprendizaje gratuitos

Volver arriba

Unidad 2 · Semanas 2–3

Sistemas de ecuaciones lineales y eliminación gaussiana

Sistemas lineales vistos por filas y por columnas, matrices aumentadas, eliminación hasta la forma escalonada, sustitución hacia atrás, forma escalonada reducida, pivotes y variables libres, y cuándo un sistema tiene una, ninguna o infinitas soluciones.

Objetivos de aprendizaje

Al terminar esta unidad podrás:

  • Escribir un sistema como Ax = b y como matriz aumentada.
  • Resolver sistemas por eliminación y sustitución hacia atrás.
  • Hallar la forma escalonada reducida e identificar pivotes y variables libres.
  • Decidir si un sistema tiene una, ninguna o infinitas soluciones.

Conceptos clave

  • Vista por filas: cada ecuación es una recta (2D) o un plano (3D) y la solución es donde se cortan. Vista por columnas: b es una combinación de las columnas de A, con pesos x.
  • Tres operaciones de fila no cambian las soluciones: intercambiar filas, multiplicar una fila por un número no nulo y sumar a una fila un múltiplo de otra.
  • La eliminación produce una forma triangular superior (escalonada); la sustitución hacia atrás resuelve después de abajo arriba. Los pivotes son las primeras entradas no nulas de cada fila.
  • Una columna sin pivote da una variable libre e infinitas soluciones. Una fila 0 = c con c ≠ 0 significa que no hay solución (incompatible).
  • Las ecuaciones de distancia al cuadrado como (x − a)² + (y − b)² = d² no son lineales, pero restar una de otra cancela x² e y² y deja ecuaciones lineales. Así es como los problemas de posicionamiento se convierten en sistemas lineales.

Ecuaciones clave

Ax=bA\mathbf{x}=\mathbf{b}
Forma matricial
Ri←Ri−ℓijRj,  ℓij=aij/ajjR_i\leftarrow R_i-\ell_{ij}R_j,\ \ \ell_{ij}=a_{ij}/a_{jj}
Paso de eliminación y multiplicador
solutions: one, none, or infinitely many\text{solutions: one, none, or infinitely many}
Los tres únicos casos
d=c td=c\,t
Distancia a partir del tiempo de viaje de la señal (GPS)

Ejemplos resueltos paso a paso

Ejemplo 2.1 Resuelve x + y + z = 6, 2x + y − z = 1, x − y + 2z = 5.

Ver los pasos
  1. R₂ − 2R₁: −y − 3z = −11. R₃ − R₁: −2y + z = −1.
  2. R₃ − 2R₂: 7z = 21, así que z = 3.
  3. Sustitución hacia atrás: −y − 9 = −11 da y = 2; luego x = 6 − 2 − 3 = 1.

Respuesta: (x, y, z) = (1, 2, 3).

Ejemplo 2.2 Hay balizas en (0, 0), (10, 0) y (0, 10) km. Un receptor mide distancias al cuadrado de 25, 65 y 45 km². Halla su posición.

Ver los pasos
  1. Ecuación 1 menos ecuación 2: x² + y² − [(x − 10)² + y²] = 25 − 65, así que 20x − 100 = −40.
  2. Ecuación 1 menos ecuación 3: 20y − 100 = −20.
  3. Resuelve las dos ecuaciones lineales.

Respuesta: (x, y) = (3, 4) km.

Ejemplo 2.3 ¿Cuántas soluciones tienen x + 2y = 3, 2x + 4y = 6? ¿Y x + 2y = 3, 2x + 4y = 7?

Ver los pasos
  1. R₂ − 2R₁ da 0 = 0 en el primer sistema: un pivote, una variable libre.
  2. En el segundo, R₂ − 2R₁ da 0 = 1.

Respuesta: Infinitas; ninguna.

Problemas de práctica

Resuelve cada uno primero en papel y luego comprueba.

  1. Resuelve x + y = 5, x − y = 1.

    Ver la respuesta

    Sumando se obtiene 2x = 6: x = 3, y = 2.

  2. Una señal tipo GPS tarda 0.00505783 s en llegar. ¿Qué distancia recorrió a 299792458 m/s?

    Ver la respuesta

    d = ct = 1516299 m, unos 1516 km.

  3. Halla la forma escalonada reducida de [[1, 2, 3], [2, 5, 8]] (aumentada).

    Ver la respuesta

    R₂ − 2R₁ = [0, 1, 2]; R₁ − 2R₂ = [1, 0, −1]: x = −1, y = 2.

  4. Un sistema 3 × 5 tiene 3 pivotes. ¿Cuántas variables libres hay?

    Ver la respuesta

    5 − 3 = 2.

Ideas erróneas comunes

  • Mito: Multiplicar una fila por 0 es una operación de fila válida.

    En realidad: Destruye información; el multiplicador debe ser no nulo.

  • Mito: Un sistema con más incógnitas que ecuaciones no tiene solución.

    En realidad: Si es compatible, tiene infinitas; también puede ser incompatible.

  • Mito: Una fila de ceros significa que no hay solución.

    En realidad: Solo una fila 0 = c con c ≠ 0 lo significa; 0 = 0 solo indica una variable libre.

Tarjetas de estudio

  • Matriz aumentada

    [A | b].

  • Pivote

    Primera entrada no nula de una fila en forma escalonada.

  • Variable libre

    Variable cuya columna no tiene pivote.

  • RREF

    Pivotes 1, ceros arriba y abajo.

  • Incompatible

    Una fila 0 = c, c ≠ 0.

  • Sustitución hacia atrás

    Resolver el sistema triangular de abajo arriba.

Cuestionario de autoevaluación

1. Una fila 0 0 0 | 5 en forma escalonada significa…

Ver la respuesta

ninguna solución. 0 = 5 es imposible.

2. ¿Cuál NO es una operación de fila?

Ver la respuesta

elevar una fila al cuadrado. Elevar al cuadrado no es lineal.

3. Una columna sin pivote da…

Ver la respuesta

una variable libre. Su variable puede elegirse libremente.

4. En la vista por columnas, Ax = b pregunta si b es…

Ver la respuesta

una combinación de las columnas de A. Ax = x₁a₁ + x₂a₂ + …

Cómo ayuda esto a salvar el planeta

El GPS es una máquina de álgebra lineal. GPS.gov informa que al menos 24 satélites operativos se mantienen disponibles el 95 % del tiempo y que se han operado 31 satélites durante más de una década, a unos 20200 km de altitud, cada uno dando dos vueltas a la Tierra al día. Están en seis planos orbitales igualmente espaciados para que los usuarios vean al menos cuatro satélites desde prácticamente cualquier punto del planeta.

La actividad de trilateración de GPS.gov explica la idea: distancia = velocidad × tiempo, con la velocidad igual a la de la luz, 299792458 m/s, y las distancias a al menos 3 satélites fijan una posición 3D. El reloj de un receptor real también es desconocido, así que en el modelo estándar de los libros las incógnitas son x, y, z y el desfase del reloj, cuatro incógnitas que necesitan cuatro satélites; restar las ecuaciones de distancia, como en el ejemplo resuelto, las convierte en un sistema lineal. Un error de reloj de solo un microsegundo desviaría una distancia en 299.8 m. Las posiciones precisas guían tractores que aplican menos fertilizante, barcos en rutas eficientes y las estaciones de medición que vigilan el nivel del mar y el hielo.

En Saving Planets: Observación de la Tierra · Espacio

Fuentes: GPS.gov: Space segment · GPS.gov: Trilateration (NOAA classroom activity)

Recursos de aprendizaje gratuitos

Volver arriba

Unidad 3 · Semana 4

Álgebra de matrices: multiplicación, inversas y factorización LU

Productos matriz–vector y matriz–matriz, las reglas (y no reglas) del álgebra de matrices, traspuestas, matrices de eliminación, la inversa y Gauss–Jordan, la fórmula de la inversa 2 × 2, A = LU y el coste de la eliminación.

Objetivos de aprendizaje

Al terminar esta unidad podrás:

  • Multiplicar matrices comprobando que los tamaños coinciden.
  • Hallar inversas con la fórmula 2 × 2 y con Gauss–Jordan.
  • Factorizar A = LU y usarla para resolver sistemas.
  • Estimar el coste de la eliminación al crecer n.

Conceptos clave

  • (AB)ᵢⱼ = producto escalar de la fila i de A con la columna j de B. Una matriz m × n por una n × p es m × p. En general AB ≠ BA.
  • La inversa A⁻¹ cumple AA⁻¹ = A⁻¹A = I. Existe exactamente cuando A tiene n pivotes (det A ≠ 0). (AB)⁻¹ = B⁻¹A⁻¹ y (AB)ᵀ = BᵀAᵀ.
  • Gauss–Jordan: reducir por filas [A | I] a [I | A⁻¹].
  • La eliminación sin intercambios de filas guarda sus multiplicadores en una L triangular inferior con unos en la diagonal, lo que da A = LU. Resuelve Lc = b y luego Ux = c; con L y U guardadas, cada nuevo término independiente es barato.
  • La eliminación en una matriz n × n cuesta unas n³/3 operaciones de multiplicación–suma; resolver con L y U guardadas cuesta unas n² por término independiente.

Ecuaciones clave

(AB)ij=∑kaikbkj(AB)_{ij}=\sum_k a_{ik}b_{kj}
Producto de matrices
[abcd]−1=1ad−bc[d−b−ca]\begin{bmatrix}a&b\\c&d\end{bmatrix}^{-1}=\frac{1}{ad-bc}\begin{bmatrix}d&-b\\-c&a\end{bmatrix}
Inversa 2 × 2
A=LUA=LU
Factorización LU
cost≈n3/3\text{cost}\approx n^3/3
Número de operaciones de la eliminación

Ejemplos resueltos paso a paso

Ejemplo 3.1 Halla AB y BA para A = [[1, 2], [3, 4]] y B = [[0, 1], [1, 0]].

Ver los pasos
  1. AB: fila 1 · columnas = (2, 1); fila 2 = (4, 3).
  2. BA intercambia las filas de A: [[3, 4], [1, 2]].

Respuesta: AB = [[2, 1], [4, 3]], BA = [[3, 4], [1, 2]], así que AB ≠ BA.

Ejemplo 3.2 Invierte A = [[4, 7], [2, 6]].

Ver los pasos
  1. ad − bc = 24 − 14 = 10.
  2. A⁻¹ = (1/10)[[6, −7], [−2, 4]].

Respuesta: A⁻¹ = [[0.6, −0.7], [−0.2, 0.4]].

Ejemplo 3.3 Factoriza A = [[2, 1], [6, 8]] como LU.

Ver los pasos
  1. Multiplicador ℓ₂₁ = 6/2 = 3.
  2. R₂ − 3R₁ = [0, 5].

Respuesta: L = [[1, 0], [3, 1]], U = [[2, 1], [0, 5]].

Problemas de práctica

Resuelve cada uno primero en papel y luego comprueba.

  1. ¿Qué tamaño tiene una matriz 3 × 2 por una 2 × 4?

    Ver la respuesta

    3 × 4.

  2. Invierte [[2, 1], [5, 3]].

    Ver la respuesta

    det = 1, inversa [[3, −1], [−5, 2]].

  3. ¿Aproximadamente cuántas operaciones requiere la eliminación para n = 1000?

    Ver la respuesta

    n³/3 ≈ 3.33 × 10⁸.

  4. Si A y B son invertibles, ¿cuánto vale (AB)⁻¹?

    Ver la respuesta

    B⁻¹A⁻¹.

Ideas erróneas comunes

  • Mito: AB = BA.

    En realidad: La multiplicación de matrices normalmente no es conmutativa.

  • Mito: (AB)⁻¹ = A⁻¹B⁻¹.

    En realidad: El orden se invierte: B⁻¹A⁻¹.

  • Mito: Los ordenadores resuelven Ax = b calculando A⁻¹.

    En realidad: Factorizan A = LU (o similar); formar A⁻¹ es más lento y menos preciso.

Tarjetas de estudio

  • Identidad

    AI = IA = A.

  • Inversa

    AA⁻¹ = I.

  • Inversa 2 × 2

    Intercambia a y d; cambia el signo de b y c; divide entre ad − bc.

  • Regla de la traspuesta

    (AB)ᵀ = BᵀAᵀ.

  • LU

    A = LU a partir de la eliminación.

  • Coste

    Unas n³/3 para la eliminación.

Cuestionario de autoevaluación

1. Una matriz 2 × 3 por una 3 × 5 es…

Ver la respuesta

2 × 5. Tamaños exteriores.

2. Una matriz es invertible exactamente cuando tiene…

Ver la respuesta

n pivotes. Un conjunto completo de pivotes.

3. Duplicar n multiplica el coste de la eliminación por aproximadamente…

Ver la respuesta

8. n³ crece 2³ = 8 veces.

4. L en A = LU contiene…

Ver la respuesta

los multiplicadores. Multiplicadores bajo la diagonal, unos en la diagonal.

Cómo ayuda esto a salvar el planeta

Los pronósticos meteorológicos y climáticos son, en esencia, enormes cálculos de álgebra lineal repetidos muchas veces al día. La NOAA informó en agosto de 2023 que sus dos supercomputadoras operativas gemelas, Dogwood en Manassas, Virginia, y Cactus en Phoenix, Arizona, se mejoraron un 20 % hasta 14.5 petaflops cada una, de modo que juntas el sistema de pronóstico puede procesar 29 mil billones de cálculos por segundo, y que con sus supercomputadoras de investigación su capacidad total es de 49 petaflops.

La regla n³/3 muestra por qué importa la velocidad. Como ejemplo, un sistema denso con n = 100000 incógnitas necesita unas 3.33 × 10¹⁴ operaciones; a un ritmo sostenido de 14.5 × 10¹⁵ operaciones por segundo son unos 0.023 s, pero duplicar n multiplica el trabajo por 8. Los modelos reales usan matrices dispersas y solucionadores ingeniosos, pero ese mismo crecimiento cúbico es lo que, según la NOAA, pagará la potencia extra: modelos de mayor resolución que mejoran pronósticos y alertas.

En Saving Planets: Introducción a la programación · Soluciones

Fuentes: NOAA: NOAA completes upgrade to weather and climate supercomputer system (Aug. 10, 2023) · MIT OpenCourseWare 18.06 Linear Algebra (Strang, Spring 2010)

Recursos de aprendizaje gratuitos

Volver arriba

Unidad 4 · Semana 5

Modelos lineales: economía de insumo–producto (Leontief), redes y balanceo

Construir modelos lineales a partir de sistemas reales: el modelo insumo–producto de Leontief x = Ax + d y su solución x = (I − A)⁻¹d, modelos ampliados ambientalmente que convierten el gasto en emisiones, flujos en redes y ajuste de ecuaciones químicas.

Objetivos de aprendizaje

Al terminar esta unidad podrás:

  • Plantear y resolver un modelo insumo–producto de Leontief.
  • Interpretar las entradas de (I − A)⁻¹.
  • Ampliar un modelo insumo–producto con intensidades de emisión para hallar una huella.
  • Escribir problemas de balance de flujos y de balance químico como sistemas lineales.

Conceptos clave

  • En un modelo insumo–producto, la columna j de la matriz de consumo A indica lo que el sector j compra a cada sector para producir un dólar de producción.
  • La producción total x debe cubrir lo que usan las industrias (Ax) más la demanda final d: x = Ax + d, así que (I − A)x = d y x = (I − A)⁻¹d.
  • La entrada (i, j) de la inversa de Leontief (I − A)⁻¹ es la producción total del sector i necesaria, directamente y a través de todas las cadenas de suministro, por un dólar de demanda final del sector j.
  • Los modelos ampliados ambientalmente multiplican la producción por una fila de intensidades de emisión e (por ejemplo kg CO₂ por dólar): huella = e·x.
  • Las leyes de conservación también dan ecuaciones lineales: el flujo que entra es igual al que sale en cada nodo de una red, y los átomos de cada elemento se equilibran en una reacción química.

Ecuaciones clave

x=Ax+d\mathbf{x}=A\mathbf{x}+\mathbf{d}
Balance insumo–producto
x=(I−A)−1d\mathbf{x}=(I-A)^{-1}\mathbf{d}
Solución de Leontief
footprint=e⋅x\text{footprint}=\mathbf{e}\cdot\mathbf{x}
Modelo ampliado ambientalmente
∑in=∑out\sum \text{in}=\sum \text{out}
Balance de flujo en un nodo

Ejemplos resueltos paso a paso

Ejemplo 4.1 Una economía de dos sectores tiene A = [[0.2, 0.3], [0.1, 0.4]] y demanda final d = (100, 50) (valores de ejemplo, millones de $). Halla la producción total.

Ver los pasos
  1. I − A = [[0.8, −0.3], [−0.1, 0.6]], det = 0.48 − 0.03 = 0.45.
  2. (I − A)⁻¹ = (1/0.45)[[0.6, 0.3], [0.1, 0.8]].
  3. x = (1/0.45)(60 + 15, 10 + 40).

Respuesta: x = (166.7, 111.1) millones de $.

Ejemplo 4.2 Con intensidades de emisión de ejemplo e = (0.5, 2.0) kg CO₂ por dólar, halla la huella de esa demanda final.

Ver los pasos
  1. huella = 0.5 × 166.7 + 2.0 × 111.1 (en kg por $ × millones de $).

Respuesta: 305.6 millones de kg = 305.6 kt CO₂ (ejemplo).

Ejemplo 4.3 Ajusta x CH₄ + y O₂ → z CO₂ + w H₂O.

Ver los pasos
  1. C: x = z. H: 4x = 2w. O: 2y = 2z + w.
  2. Set x = 1: z = 1, w = 2, y = 2.

Respuesta: CH₄ + 2 O₂ → CO₂ + 2 H₂O.

Problemas de práctica

Resuelve cada uno primero en papel y luego comprueba.

  1. En el modelo resuelto, ¿cuánta producción del sector 1 requiere en total 1 $ de demanda final del sector 2?

    Ver la respuesta

    Entrada (1, 2) de (I − A)⁻¹ = 0.3/0.45 = 0.667.

  2. Una economía de un solo sector usa 0.25 de su propia producción. ¿Qué producción total satisface una demanda de 300?

    Ver la respuesta

    x = 300/(1 − 0.25) = 400.

  3. Un nodo tiene entradas de 40 y 25 y una salida f₁ más una salida de 30. Halla f₁.

    Ver la respuesta

    f₁ = 65 − 30 = 35.

  4. Ajusta H₂ + O₂ → H₂O.

    Ver la respuesta

    2 H₂ + O₂ → 2 H₂O.

Ideas erróneas comunes

  • Mito: La producción total es igual a la demanda final.

    En realidad: Las industrias también consumen la producción de las demás, así que x es mayor que d.

  • Mito: La huella proviene solo del sector al que compras.

    En realidad: La inversa de Leontief suma todos los proveedores aguas arriba.

  • Mito: Cualquier A da una economía sensata.

    En realidad: Las sumas de las columnas deben ser menores que 1 para que (I − A)⁻¹ tenga entradas no negativas.

Tarjetas de estudio

  • Matriz de consumo

    Insumos por dólar de producción.

  • Modelo de Leontief

    x = Ax + d.

  • Inversa de Leontief

    (I − A)⁻¹, requerimientos totales.

  • Modelo ampliado

    huella = e·x.

  • Balance de flujo

    Entrada = salida en cada nodo.

  • Sistema homogéneo

    Ax = 0.

Cuestionario de autoevaluación

1. La solución de Leontief es…

Ver la respuesta

x = (I − A)⁻¹d. A partir de x = Ax + d.

2. La columna j de A describe…

Ver la respuesta

lo que compra el sector j por dólar de producción. Insumos por unidad de producción.

3. En un modelo ampliado, una huella es…

Ver la respuesta

e·x. Intensidad por producción, sumada.

4. Ajustar una ecuación química es un…

Ver la respuesta

sistema lineal homogéneo. El recuento de átomos da Ax = 0.

Cómo ayuda esto a salvar el planeta

La EPA de EE. UU. usa exactamente este modelo para rastrear la contaminación a través de la economía. Según la EPA, sus modelos USEEIO combinan datos de transacciones económicas entre 389 sectores industriales con datos ambientales de esos sectores, que abarcan el uso de suelo, agua, energía y minerales, la contaminación del aire, los nutrientes y los tóxicos, para evaluar los impactos de bienes y servicios. La EPA los describe como útiles para calcular la huella de carbono, y respaldan su trabajo de gestión sostenible de materiales.

Una matriz de consumo de 389 sectores tiene 389² = 151321 entradas, demasiadas para trabajar a mano pero fáciles para los métodos de eliminación y LU de la unidad 3. La inversa de Leontief es lo que permite a una ciudad, empresa u hogar ver las emisiones ocultas en sus cadenas de suministro, no solo en la puerta de la fábrica.

En Saving Planets: Soluciones · Química General I

Fuentes: U.S. EPA: US Environmentally-Extended Input-Output (USEEIO) models · Interactive Linear Algebra, 1.1 Systems of linear equations

Recursos de aprendizaje gratuitos

Volver arriba

Unidad 5 · Semanas 6–7

Espacios vectoriales, subespacios, independencia, base, rango y los cuatro subespacios fundamentales

Espacios vectoriales y subespacios, espacio columna y espacio nulo, resolución completa de Ax = 0 y Ax = b, independencia lineal, base y dimensión, rango, el teorema del rango y los cuatro subespacios fundamentales de Strang.

Objetivos de aprendizaje

Al terminar esta unidad podrás:

  • Comprobar si un conjunto es un subespacio.
  • Hallar bases del espacio columna y del espacio nulo de una matriz.
  • Comprobar la independencia lineal de vectores.
  • Usar el teorema del rango y las dimensiones de los cuatro subespacios fundamentales.

Conceptos clave

  • Un subespacio contiene el vector cero y es cerrado bajo la suma y la multiplicación por escalares. Las rectas y planos que pasan por el origen son subespacios; los que no pasan por el origen no lo son.
  • Espacio columna C(A): todos los Ax, el conjunto de b para los que Ax = b tiene solución. Espacio nulo N(A): todos los x con Ax = 0.
  • Los vectores son linealmente independientes si la única combinación que da cero es la de todos los coeficientes nulos. Una base es un conjunto generador independiente; toda base de un espacio tiene el mismo número de vectores, su dimensión.
  • Rango r = número de pivotes = dim C(A). Teorema del rango: r + dim N(A) = n, el número de columnas.
  • Cuatro subespacios fundamentales de una matriz m × n de rango r: C(A) (dim r) y N(Aᵀ) (dim m − r) en ℝᵐ; C(Aᵀ) (dim r) y N(A) (dim n − r) en ℝⁿ. N(A) es ortogonal al espacio fila.

Ecuaciones clave

C(A)={Ax},  N(A)={x:Ax=0}C(A)=\{A\mathbf{x}\},\ \ N(A)=\{\mathbf{x}:A\mathbf{x}=0\}
Espacio columna y espacio nulo
r+dim⁡N(A)=nr+\dim N(A)=n
Teorema del rango
x=xp+xn\mathbf{x}=\mathbf{x}_p+\mathbf{x}_n
Solución completa de Ax = b
dim⁡N(AT)=m−r\dim N(A^T)=m-r
Espacio nulo izquierdo

Ejemplos resueltos paso a paso

Ejemplo 5.1 Halla el rango y una base del espacio nulo de A = [[1, 2, 3], [2, 4, 6], [1, 0, 1]].

Ver los pasos
  1. R₂ − 2R₁ = 0; R₃ − R₁ = [0, −2, −2]. Dos pivotes, así que r = 2.
  2. Ax = 0: x + z = 0 y −2y − 2z = 0, con z libre.
  3. Toma z = 1: x = −1, y = −1.

Respuesta: r = 2, N(A) generado por (−1, −1, 1); dim N(A) = 3 − 2 = 1.

Ejemplo 5.2 ¿Son independientes (1, 2, 3), (4, 5, 6), (7, 8, 9)?

Ver los pasos
  1. Reduce por filas la matriz con esas columnas: R₂ − 2R₁ y R₃ − 3R₁ dan las filas [0, −3, −6] y [0, −6, −12], que son múltiplos.
  2. Solo dos pivotes.

Respuesta: No: (7, 8, 9) = 2(4, 5, 6) − (1, 2, 3).

Ejemplo 5.3 Una matriz 5 × 8 tiene rango 3. Da las dimensiones de sus cuatro subespacios fundamentales.

Ver los pasos
  1. C(A) y C(Aᵀ): r = 3.
  2. N(A): n − r = 8 − 3; N(Aᵀ): m − r = 5 − 3.

Respuesta: 3, 3, 5 y 2.

Problemas de práctica

Resuelve cada uno primero en papel y luego comprueba.

  1. ¿Es la recta y = 2x + 1 un subespacio de ℝ²?

    Ver la respuesta

    No: no contiene (0, 0).

  2. Una matriz 4 × 4 tiene un espacio nulo de dimensión 2. ¿Cuál es su rango?

    Ver la respuesta

    4 − 2 = 2.

  3. ¿Cuál es la dimensión del espacio de matrices 2 × 2?

    Ver la respuesta

    4.

  4. ¿Pueden ser independientes 4 vectores de ℝ³?

    Ver la respuesta

    No: como mucho 3.

Ideas erróneas comunes

  • Mito: Cualquier plano de ℝ³ es un subespacio.

    En realidad: Solo los planos que pasan por el origen.

  • Mito: El rango cuenta las entradas no nulas.

    En realidad: Cuenta los pivotes (columnas independientes).

  • Mito: El espacio nulo vive en el mismo espacio que b.

    En realidad: N(A) está en ℝⁿ (el lado de x); b está en ℝᵐ.

Tarjetas de estudio

  • Subespacio

    Contiene 0; cerrado bajo + y escalado.

  • Espacio columna

    All Ax.

  • Espacio nulo

    Todos los x con Ax = 0.

  • Base

    Conjunto generador independiente.

  • Rango

    Número de pivotes.

  • Teorema del rango

    r + dim N(A) = n.

Cuestionario de autoevaluación

1. El teorema del rango dice que r + dim N(A) es igual al…

Ver la respuesta

n. Número de columnas.

2. Ax = b tiene solución exactamente cuando b está en…

Ver la respuesta

C(A). b debe ser una combinación de las columnas.

3. El espacio nulo es ortogonal al…

Ver la respuesta

el espacio fila. Cada fila por x da 0.

4. ¿Cuántos vectores tiene una base de ℝ³?

Ver la respuesta

3. Dimensión 3.

Cómo ayuda esto a salvar el planeta

La dimensión cuenta incógnitas, y mejores pronósticos necesitan más. La NOAA dijo que la mejora de sus supercomputadoras en 2023 permitiría al Global Forecast System de EE. UU. mejorar su resolución horizontal de 13 km a 9 km, lo que ayuda a ver y simular fenómenos de menor escala.

Una malla con espaciado h tiene unas (1/h)² columnas por unidad de área, así que pasar de 13 km a 9 km multiplica el número de columnas de la malla por (13/9)² = 2.09. Cada una de esas columnas lleva temperatura, viento, humedad y presión en muchos niveles, de modo que el vector de estado de un pronóstico vive en un espacio de dimensión muchísimo mayor, y las matrices que lo actualizan crecen en consecuencia.

En Saving Planets: Cálculo III · Ecuaciones diferenciales

Fuentes: NOAA: NOAA completes upgrade to weather and climate supercomputer system (Aug. 10, 2023) · Interactive Linear Algebra, 2.7 Basis and dimension

Recursos de aprendizaje gratuitos

Volver arriba

Unidad 6 · Semanas 8–9

Ortogonalidad, proyecciones, mínimos cuadrados y Gram–Schmidt

Vectores y subespacios ortogonales, proyección sobre una recta y sobre un subespacio, las ecuaciones normales AᵀAx̂ = Aᵀb, ajuste de rectas y curvas a datos por mínimos cuadrados, bases ortonormales, Gram–Schmidt y A = QR.

Objetivos de aprendizaje

Al terminar esta unidad podrás:

  • Proyectar un vector sobre una recta y sobre un espacio columna.
  • Ajustar una recta a datos con las ecuaciones normales.
  • Construir una base ortonormal con Gram–Schmidt.
  • Explicar por qué los mínimos cuadrados son una proyección.

Conceptos clave

  • La proyección de b sobre la recta que pasa por a es p = (a·b/a·a) a; el error e = b − p es ortogonal a a.
  • Cuando Ax = b no tiene solución (más ecuaciones que incógnitas, datos con ruido), elige x̂ para que ‖b − Ax̂‖ sea lo más pequeño posible. El error debe ser ortogonal a C(A), lo que da las ecuaciones normales AᵀAx̂ = Aᵀb.
  • Para una recta y = C + Dt por los puntos (tᵢ, yᵢ), las columnas de A son una de unos y la de los valores tᵢ.
  • Los vectores ortonormales qᵢ cumplen qᵢ·qⱼ = 0 (i ≠ j) y ‖qᵢ‖ = 1. Si Q tiene columnas ortonormales, QᵀQ = I y las proyecciones se vuelven p = QQᵀb.
  • Gram–Schmidt resta a cada vector nuevo sus proyecciones sobre los anteriores y luego normaliza. El resultado es A = QR con R triangular superior.

Ecuaciones clave

p=a⋅ba⋅a a\mathbf{p}=\frac{\mathbf{a}\cdot\mathbf{b}}{\mathbf{a}\cdot\mathbf{a}}\,\mathbf{a}
Proyección sobre una recta
ATAx^=ATbA^TA\hat{\mathbf{x}}=A^T\mathbf{b}
Ecuaciones normales
P=A(ATA)−1ATP=A(A^TA)^{-1}A^T
Matriz de proyección
A=QRA=QR
Factorización de Gram–Schmidt

Ejemplos resueltos paso a paso

Ejemplo 6.1 Proyecta b = (2, 3) sobre a = (1, 1).

Ver los pasos
  1. a·b = 5, a·a = 2.
  2. p = 2.5(1, 1); e = b − p = (−0.5, 0.5), y e·a = 0.

Respuesta: p = (2.5, 2.5).

Ejemplo 6.2 Ajusta y = C + Dt a (0, 1), (1, 2), (2, 4) por mínimos cuadrados.

Ver los pasos
  1. A tiene filas (1, 0), (1, 1), (1, 2); AᵀA = [[3, 3], [3, 5]], Aᵀb = (7, 10).
  2. det = 6; C = (5 × 7 − 3 × 10)/6, D = (3 × 10 − 3 × 7)/6.

Respuesta: C = 0.833, D = 1.5: y = 0.833 + 1.5t.

Ejemplo 6.3 Aplica Gram–Schmidt a a₁ = (1, 1, 0), a₂ = (1, 0, 1).

Ver los pasos
  1. q₁ = (1, 1, 0)/√2.
  2. a₂ − (a₂·q₁)q₁ = (1, 0, 1) − ½(1, 1, 0) = (0.5, −0.5, 1), de longitud √1.5 = 1.225.

Respuesta: q₁ = (0.707, 0.707, 0), q₂ = (0.408, −0.408, 0.816).

Problemas de práctica

Resuelve cada uno primero en papel y luego comprueba.

  1. Proyecta (3, 4) sobre (1, 0).

    Ver la respuesta

    (3, 0).

  2. Escribe las ecuaciones normales para ajustar y = Dt (sin ordenada en el origen) a (1, 2), (2, 3).

    Ver la respuesta

    Σt² D = Σty: 5D = 8, D = 1.6.

  3. ¿Es ortogonal Q = [[0, 1], [1, 0]]?

    Ver la respuesta

    Yes: QᵀQ = I.

  4. Si b ya está en C(A), ¿cuál es el error de mínimos cuadrados?

    Ver la respuesta

    Cero: p = b.

Ideas erróneas comunes

  • Mito: Los mínimos cuadrados resuelven Ax = b exactamente.

    En realidad: Hallan el Ax̂ más cercano cuando no existe solución exacta.

  • Mito: Cualquier base facilita las proyecciones.

    En realidad: Las bases ortonormales sí; si no, necesitas (AᵀA)⁻¹.

  • Mito: Un ajuste por una recta significa que los datos son lineales.

    En realidad: Revisa los residuos; una curva puede ajustar mejor (añade una columna t²).

Tarjetas de estudio

  • Proyección sobre una recta

    p = (a·b/a·a)a.

  • Ecuaciones normales

    AᵀAx̂ = Aᵀb.

  • Matriz de proyección

    P = A(AᵀA)⁻¹Aᵀ.

  • Ortonormal

    qᵢ·qⱼ = 0, ‖qᵢ‖ = 1.

  • Gram–Schmidt

    Restar proyecciones, normalizar.

  • QR

    A = QR.

Cuestionario de autoevaluación

1. Las ecuaciones normales son…

Ver la respuesta

AᵀAx̂ = Aᵀb. Error ortogonal a C(A).

2. El vector de error de mínimos cuadrados es ortogonal al…

Ver la respuesta

espacio columna. Eso es lo que lo hace el más corto.

3. Para Q con columnas ortonormales, QᵀQ =…

Ver la respuesta

I. Productos escalares 0 o 1.

4. Gram–Schmidt produce…

Ver la respuesta

A = QR. Q ortonormal, R triangular.

Cómo ayuda esto a salvar el planeta

La subida del nivel del mar se mide ajustando curvas a datos satelitales con ruido. La NASA informa que los altímetros satelitales, TOPEX/Poseidon, Jason-1, -2 y -3 y ahora Sentinel-6 Michael Freilich, han seguido el nivel medio global del mar desde 1993, que ha subido 114.4 mm desde entonces, y que el ritmo de subida ha pasado de 0.08 pulgadas por año (0.20 cm por año) en 1993 a 0.17 pulgadas por año (0.44 cm por año) ahora, unas 2.2 veces más rápido.

Un ajuste lineal por mínimos cuadrados da un ritmo medio; añadir una columna t² a A ajusta una curva cuya pendiente variable mide la aceleración que describe la NASA. La Práctica D ajusta datos de ejemplo con la forma del registro de la NASA usando exactamente las ecuaciones normales de esta unidad.

En Saving Planets: Océanos · Observación de la Tierra

Fuentes: NASA: Global mean sea level (Vital Signs) · Interactive Linear Algebra, 6.5 The method of least squares

Recursos de aprendizaje gratuitos

Volver arriba

Unidad 7 · Semana 10

Determinantes: propiedades, cofactores, regla de Cramer y volumen

El determinante a partir de sus propiedades, su cálculo por eliminación y por desarrollo por cofactores, det(AB) = det A det B, el criterio de invertibilidad, la regla de Cramer y el determinante como factor de escala de área o volumen de una aplicación lineal.

Objetivos de aprendizaje

Al terminar esta unidad podrás:

  • Calcular determinantes 2 × 2 y 3 × 3 por cofactores y por eliminación.
  • Usar las propiedades de los determinantes bajo operaciones de fila.
  • Resolver sistemas pequeños con la regla de Cramer.
  • Interpretar |det A| como un factor de escala de área o volumen.

Conceptos clave

  • det A ≠ 0 exactamente cuando A es invertible. El determinante de una matriz triangular es el producto de sus entradas diagonales, así que det A = ± (producto de los pivotes).
  • Operaciones de fila: un intercambio cambia el signo, multiplicar una fila por c multiplica det por c, y sumar a una fila un múltiplo de otra deja det sin cambios.
  • det(AB) = det A · det B, det(Aᵀ) = det A y det(A⁻¹) = 1/det A.
  • Desarrollo por cofactores a lo largo de una fila: det A = Σ aᵢⱼCᵢⱼ con Cᵢⱼ = (−1)ⁱ⁺ʲ por el menor. Regla de Cramer: xᵢ = det(Bᵢ)/det A, donde Bᵢ tiene b en la columna i.
  • Una aplicación lineal x ↦ Ax multiplica las áreas (2D) o los volúmenes (3D) por |det A|; un determinante negativo además invierte la orientación.

Ecuaciones clave

det⁡[abcd]=ad−bc\det\begin{bmatrix}a&b\\c&d\end{bmatrix}=ad-bc
Determinante 2 × 2
det⁡(AB)=det⁡A det⁡B\det(AB)=\det A\,\det B
Regla del producto
xi=det⁡Bidet⁡Ax_i=\frac{\det B_i}{\det A}
Regla de Cramer
area(AS)=∣det⁡A∣ area(S)\text{area}(AS)=|\det A|\,\text{area}(S)
Factor de escala de áreas

Ejemplos resueltos paso a paso

Ejemplo 7.1 Halla det [[1, 2, 0], [0, 1, 3], [2, 0, 1]] por cofactores a lo largo de la fila 1.

Ver los pasos
  1. 1 × (1 × 1 − 3 × 0) = 1.
  2. − 2 × (0 × 1 − 3 × 2) = −2 × (−6) = 12.
  3. + 0 × (…) = 0.

Respuesta: det = 13.

Ejemplo 7.2 Resuelve 2x + y = 5, x + 3y = 10 con la regla de Cramer.

Ver los pasos
  1. det A = 6 − 1 = 5.
  2. x = det[[5, 1], [10, 3]]/5 = (15 − 10)/5.
  3. y = det[[2, 5], [1, 10]]/5 = (20 − 5)/5.

Respuesta: x = 1, y = 3.

Ejemplo 7.3 Halla el área del paralelogramo de lados (3, 0) y (1, 2).

Ver los pasos
  1. Área = |det [[3, 1], [0, 2]]| = |6 − 0|.

Respuesta: 6 unidades cuadradas.

Problemas de práctica

Resuelve cada uno primero en papel y luego comprueba.

  1. ¿det [[3, 1], [2, 4]]?

    Ver la respuesta

    12 − 2 = 10.

  2. ¿det [[2, 0, 1], [1, 3, 2], [1, 1, 1]]?

    Ver la respuesta

    2(3 − 2) − 0 + 1(1 − 3) = 0, así que es singular.

  3. Si det A = 4 para una A 3 × 3, ¿cuánto vale det(2A)?

    Ver la respuesta

    2³ × 4 = 32.

  4. Si det A = 5, ¿cuánto vale det(A⁻¹)?

    Ver la respuesta

    1/5 = 0.2.

Ideas erróneas comunes

  • Mito: det(A + B) = det A + det B.

    En realidad: Falso en general; solo det(AB) se separa.

  • Mito: det(2A) = 2 det A.

    En realidad: Para n × n, det(2A) = 2ⁿ det A.

  • Mito: El desarrollo por cofactores es el método rápido.

    En realidad: Requiere unos n! pasos; la eliminación unos n³/3.

Tarjetas de estudio

  • det 2 × 2

    ad − bc.

  • Invertibilidad

    det A ≠ 0.

  • Regla del producto

    det(AB) = det A det B.

  • Intercambio de filas

    Cambia el signo.

  • Regla de Cramer

    xᵢ = det Bᵢ/det A.

  • Factor de área

    |det A|.

Cuestionario de autoevaluación

1. Intercambiar dos filas multiplica det por…

Ver la respuesta

−1. El signo cambia.

2. det A = 0 significa que A es…

Ver la respuesta

singular. No invertible.

3. El det de una matriz triangular es…

Ver la respuesta

el producto de la diagonal. Desarrolla a lo largo de la diagonal.

4. Una aplicación con det = 3 cambia las áreas en un factor de…

Ver la respuesta

3. |det A|.

Cómo ayuda esto a salvar el planeta

Convertir píxeles en superficie terrestre es un determinante. El USGS indica para Landsat 9 bandas multiespectrales de 30 m y una banda pancromática de 15 m. La aplicación de los índices de píxel (i, j) a coordenadas del terreno es, localmente, una matriz; para una malla alineada de 30 m es diag(30, 30) m, cuyo determinante es 900, así que cada píxel cubre 900 m². La malla pancromática de 15 m tiene determinante 225 m², de modo que cuatro píxeles pancromáticos caben en un píxel multiespectral.

Contar los píxeles clasificados como bosque y multiplicar por este factor de escala es como los registros satelitales miden la pérdida de bosque en hectáreas (1 ha = 10000 m², unos 11.1 píxeles Landsat). Cuando una imagen se rota o se cizalla para encajar en un mapa, el determinante de la nueva matriz mantiene correcta el área.

En Saving Planets: Deforestación · Observación de la Tierra

Fuentes: USGS: Landsat 9 · Interactive Linear Algebra, 4.3 Determinants and volumes

Recursos de aprendizaje gratuitos

Volver arriba

Unidad 8 · Semana 11

Valores propios y vectores propios

Los vectores propios como direcciones que una matriz solo estira, la ecuación característica det(A − λI) = 0, las comprobaciones con traza y determinante, valores propios de matrices especiales, valores propios complejos y rotación, el método de la potencia y los modos naturales de estructuras que vibran.

Objetivos de aprendizaje

Al terminar esta unidad podrás:

  • Hallar valores propios con det(A − λI) = 0 y vectores propios en el espacio nulo de A − λI.
  • Comprobar valores propios con la traza y el determinante.
  • Reconocer valores propios reales, repetidos y complejos.
  • Usar el método de la potencia para hallar el valor propio dominante.

Conceptos clave

  • Av = λv con v ≠ 0: el vector propio v mantiene su dirección y se escala por el valor propio λ.
  • Los valores propios son las raíces de det(A − λI) = 0, un polinomio de grado n. Para cada λ, los vectores propios llenan el espacio nulo de A − λI.
  • Comprobaciones: los valores propios suman la traza (suma de la diagonal) y su producto es det A.
  • Las matrices simétricas tienen valores propios reales y vectores propios ortogonales. Las matrices de rotación tienen valores propios complejos. Las matrices triangulares muestran sus valores propios en la diagonal.
  • Método de la potencia: multiplica repetidamente un vector inicial por A y reescálalo; se alinea con el vector propio del mayor |λ|. En vibraciones, Kv = ω²Mv da las frecuencias naturales ω y las formas modales v.

Ecuaciones clave

Av=λvA\mathbf{v}=\lambda\mathbf{v}
Ecuación de valores propios
det⁡(A−λI)=0\det(A-\lambda I)=0
Ecuación característica
λ1+⋯+λn=tr⁡A,  λ1⋯λn=det⁡A\lambda_1+\cdots+\lambda_n=\operatorname{tr}A,\ \ \lambda_1\cdots\lambda_n=\det A
Comprobaciones con traza y determinante
Kv=ω2MvK\mathbf{v}=\omega^2M\mathbf{v}
Modos naturales de una estructura

Ejemplos resueltos paso a paso

Ejemplo 8.1 Halla los valores y vectores propios de A = [[2, 1], [1, 2]].

Ver los pasos
  1. det(A − λI) = (2 − λ)² − 1 = λ² − 4λ + 3 = 0.
  2. λ = 3: A − 3I = [[−1, 1], [1, −1]] da v = (1, 1). λ = 1: v = (1, −1).
  3. Comprobación: 3 + 1 = 4 = traza; 3 × 1 = 3 = det.

Respuesta: λ = 3 con (1, 1); λ = 1 con (1, −1).

Ejemplo 8.2 Halla los valores propios de [[4, 1], [2, 3]].

Ver los pasos
  1. λ² − 7λ + 10 = 0 (traza 7, det 10).
  2. (λ − 5)(λ − 2) = 0.

Respuesta: λ = 5 y λ = 2.

Ejemplo 8.3 Un edificio de ejemplo de dos plantas tiene masas de piso iguales m y rigidez de entrepiso k, con K/m = (k/m)[[2, −1], [−1, 1]] y k/m = 100 s⁻². Halla la frecuencia fundamental.

Ver los pasos
  1. Valores propios de [[2, −1], [−1, 1]]: λ = (3 ± √5)/2 = 0.382 y 2.618.
  2. ω₁² = 100 × 0.382 = 38.2 s⁻², ω₁ = 6.18 rad/s.
  3. f₁ = ω₁/(2π), T₁ = 1/f₁.

Respuesta: f₁ = 0.984 Hz, periodo T₁ = 1.02 s (valores de ejemplo).

Problemas de práctica

Resuelve cada uno primero en papel y luego comprueba.

  1. ¿Valores propios de [[5, 0], [0, −2]]?

    Ver la respuesta

    5 y −2 (diagonal).

  2. Una matriz 3 × 3 tiene valores propios 1, 2, 3. Halla su traza y su determinante.

    Ver la respuesta

    Traza 6, det 6.

  3. ¿Valores propios de la rotación de 90° [[0, −1], [1, 0]]?

    Ver la respuesta

    λ² + 1 = 0: λ = ±i.

  4. Si Av = 4v, ¿cuánto vale A²v?

    Ver la respuesta

    16v.

Ideas erróneas comunes

  • Mito: El vector cero es un vector propio.

    En realidad: Los vectores propios deben ser no nulos.

  • Mito: La reducción por filas conserva los valores propios.

    En realidad: Los cambia; usa det(A − λI).

  • Mito: Toda matriz real tiene valores propios reales.

    En realidad: Las rotaciones dan complejos; las matrices simétricas son el caso seguro.

Tarjetas de estudio

  • Vector propio

    Av = λv, v ≠ 0.

  • Ecuación característica

    det(A − λI) = 0.

  • Comprobación con la traza

    Σλ = traza.

  • Comprobación con el determinante

    Πλ = det A.

  • Método de la potencia

    La multiplicación repetida halla el vector propio dominante.

  • Modos naturales

    Kv = ω²Mv.

Cuestionario de autoevaluación

1. Los valores propios de una matriz triangular son…

Ver la respuesta

sus entradas diagonales. det(A − λI) es un producto de términos diagonales.

2. La suma de los valores propios es igual a…

Ver la respuesta

la traza. Comprobación con la traza.

3. El método de la potencia halla…

Ver la respuesta

el valor propio dominante. Gana el mayor |λ|.

4. Las matrices simétricas tienen…

Ver la respuesta

valores propios reales y vectores propios ortogonales. Teorema espectral.

Cómo ayuda esto a salvar el planeta

Los ingenieros sísmicos leen los edificios a través de sus valores propios. Un estudio del USGS de Mehmet Çelebi (2022) usó los periodos fundamentales y los porcentajes de amortiguamiento crítico de 41 edificios altos de varios países. Encontró una clara tendencia lineal entre la altura del edificio y el periodo fundamental para edificios de acero y de hormigón armado, aunque con gran dispersión por las diferencias de diseño estructural, y que el amortiguamiento crítico era claramente inferior al 3 %, en línea con otros estudios recientes sobre edificios altos.

El periodo fundamental lo fija el menor valor propio del problema rigidez–masa, como en el ejemplo del edificio de dos plantas. Un amortiguamiento bajo significa que un edificio alto puede vibrar mucho tiempo cuando ondas sísmicas de periodo largo coinciden con ese periodo; por eso los edificios instrumentados y sus modos medidos orientan ciudades más seguras y duraderas.

En Saving Planets: Ecuaciones diferenciales · Dinámica

Fuentes: USGS: Preliminary study of correlation of natural periods and damping percentages of tall buildings in several countries (Çelebi, 2022) · Interactive Linear Algebra, 5.1 Eigenvalues and eigenvectors

Recursos de aprendizaje gratuitos

Volver arriba

Unidad 9 · Semana 12

Diagonalización, potencias de matrices, cadenas de Markov y sistemas dinámicos lineales

Diagonalizar A = PDP⁻¹, calcular rápidamente potencias Aᵏ, matrices simétricas y el teorema espectral A = QΛQᵀ, matrices de Markov (estocásticas) y estados estacionarios, y el comportamiento a largo plazo de xₖ₊₁ = Axₖ.

Objetivos de aprendizaje

Al terminar esta unidad podrás:

  • Diagonalizar una matriz con n vectores propios independientes.
  • Calcular Aᵏ y predecir el comportamiento a largo plazo de xₖ₊₁ = Axₖ.
  • Hallar el estado estacionario de una cadena de Markov.
  • Usar el teorema espectral para matrices simétricas.

Conceptos clave

  • Si A tiene n vectores propios independientes (las columnas de P), entonces A = PDP⁻¹ con los valores propios en la diagonal de D, y Aᵏ = PDᵏP⁻¹.
  • Escribe x₀ = c₁v₁ + … + cₙvₙ. Entonces xₖ = Aᵏx₀ = c₁λ₁ᵏv₁ + … + cₙλₙᵏvₙ; el mayor |λ| domina al crecer k.
  • Una matriz de Markov tiene entradas no negativas y columnas que suman 1. Siempre tiene λ = 1, y los demás valores propios cumplen |λ| ≤ 1. El estado estacionario es el vector propio de λ = 1, escalado para que sus entradas sumen 1.
  • El segundo mayor |λ| fija lo rápido que una cadena de Markov olvida su inicio: la diferencia se reduce como |λ₂|ᵏ.
  • Teorema espectral: una matriz simétrica real puede escribirse A = QΛQᵀ con vectores propios ortonormales en Q.

Ecuaciones clave

A=PDP−1,  Ak=PDkP−1A=PDP^{-1},\ \ A^k=PD^kP^{-1}
Diagonalización y potencias
xk=c1λ1kv1+⋯+cnλnkvn\mathbf{x}_k=c_1\lambda_1^k\mathbf{v}_1+\cdots+c_n\lambda_n^k\mathbf{v}_n
Solución de x_{k+1} = Ax_k
Pq=q,  ∑qi=1P\mathbf{q}=\mathbf{q},\ \ \textstyle\sum q_i=1
Estado estacionario de Markov
A=QΛQTA=Q\Lambda Q^T
Teorema espectral

Ejemplos resueltos paso a paso

Ejemplo 9.1 Usa la diagonalización para hallar A¹⁰ para A = [[2, 1], [1, 2]].

Ver los pasos
  1. De la unidad 8: λ = 3, 1 con v = (1, 1), (1, −1). P = [[1, 1], [1, −1]], P⁻¹ = ½[[1, 1], [1, −1]].
  2. A¹⁰ = P diag(3¹⁰, 1) P⁻¹; cada entrada es (3¹⁰ ± 1)/2.

Respuesta: A¹⁰ = [[29525, 29524], [29524, 29525]].

Ejemplo 9.2 Una matriz de Markov es P = [[0.8, 0.3], [0.2, 0.7]] (valores de ejemplo). Halla el estado estacionario.

Ver los pasos
  1. (P − I)q = 0: −0.2q₁ + 0.3q₂ = 0, así que q₁ = 1.5q₂.
  2. q₁ + q₂ = 1 da q₂ = 0.4.

Respuesta: q = (0.6, 0.4).

Ejemplo 9.3 Partiendo de x₀ = (1, 0), halla x₂ para la misma P, y el xₖ general.

Ver los pasos
  1. x₁ = (0.8, 0.2); x₂ = (0.64 + 0.06, 0.16 + 0.14).
  2. Los valores propios son 1 y 0.5 (traza 1.5, det 0.5), así que xₖ = (0.6, 0.4) + 0.5ᵏ(0.4, −0.4).

Respuesta: x₂ = (0.70, 0.30); xₖ → (0.6, 0.4).

Problemas de práctica

Resuelve cada uno primero en papel y luego comprueba.

  1. ¿Es [[0.5, 0.2], [0.5, 0.8]] una matriz de Markov?

    Ver la respuesta

    Sí: entradas no negativas, columnas que suman 1.

  2. Halla su estado estacionario.

    Ver la respuesta

    0.5q₁ = 0.2q₂, así que q = (2/7, 5/7) = (0.286, 0.714).

  3. Si A = PDP⁻¹ con D = diag(0.9, 0.5), ¿qué pasa con Aᵏx para k grande?

    Ver la respuesta

    Decae a cero, dominado por 0.9ᵏ.

  4. ¿Se puede diagonalizar [[1, 1], [0, 1]]?

    Ver la respuesta

    No: λ = 1 doble pero solo un vector propio independiente.

Ideas erróneas comunes

  • Mito: Toda matriz se puede diagonalizar.

    En realidad: Necesita n vectores propios independientes; [[1, 1], [0, 1]] solo tiene uno.

  • Mito: El estado estacionario de una cadena de Markov depende del inicio.

    En realidad: En una cadena regular no; todo inicio converge al mismo q.

  • Mito: Aᵏ significa elevar cada entrada a la potencia k.

    En realidad: Solo las matrices diagonales funcionan así; por eso diagonalizamos.

Tarjetas de estudio

  • Diagonalización

    A = PDP⁻¹.

  • Potencia de una matriz

    Aᵏ = PDᵏP⁻¹.

  • Matriz de Markov

    No negativa, columnas que suman 1.

  • Estado estacionario

    Pq = q, entradas que suman 1.

  • Teorema espectral

    A = QΛQᵀ.

  • Matriz defectuosa

    Demasiados pocos vectores propios para diagonalizar.

Cuestionario de autoevaluación

1. Toda matriz de Markov tiene el valor propio…

Ver la respuesta

1. Las columnas suman 1.

2. Aᵏ para A = PDP⁻¹ es…

Ver la respuesta

PDᵏP⁻¹. Los P⁻¹P interiores se cancelan.

3. La velocidad de convergencia de una cadena de Markov la fija…

Ver la respuesta

|λ₂|. La diferencia se reduce como |λ₂|ᵏ.

4. Las matrices simétricas se diagonalizan como…

Ver la respuesta

QΛQᵀ. Teorema espectral.

Cómo ayuda esto a salvar el planeta

El destino de los materiales es una cadena de Markov. La EPA informa que Estados Unidos generó 292.4 millones de toneladas de residuos sólidos urbanos en 2018, 4.9 libras por persona y día. Unos 69 millones de toneladas se reciclaron y 25 millones de toneladas se compostaron, una tasa combinada del 32.1 %; unos 35 millones de toneladas (11.8 %) se quemaron con recuperación de energía; y más de 146 millones de toneladas (50 %) fueron a vertedero.

Si las proporciones de cada año se mantuvieran fijas, un material dejaría de usarse hacia esos destinos con esas probabilidades, y solo la parte reciclada volvería a la economía. La Práctica E construye una matriz de transición de ejemplo con una proporción de recuperación del 32.1 % y muestra, con potencias de matrices y el vector propio de λ = 1, cómo aumentar el reciclaje cambia cuánto tiempo sigue en uso el material y cuánto llega al vertedero.

En Saving Planets: Soluciones

Fuentes: U.S. EPA: National overview, facts and figures on materials, wastes and recycling · Interactive Linear Algebra, 5.6 Stochastic matrices

Recursos de aprendizaje gratuitos

Volver arriba

Unidad 10 · Semanas 13–14

La descomposición en valores singulares, el análisis de componentes principales y las matrices definidas positivas

Matrices definidas positivas y AᵀA, la SVD A = UΣVᵀ, valores singulares y la geometría de una matriz como rotar–estirar–rotar, la mejor aproximación de rango bajo y la compresión de imágenes, el número de condición y el análisis de componentes principales de datos.

Objetivos de aprendizaje

Al terminar esta unidad podrás:

  • Comprobar si una matriz simétrica es definida positiva.
  • Calcular valores singulares a partir de los valores propios de AᵀA.
  • Construir una aproximación de rango k y medir lo que conserva.
  • Realizar un pequeño análisis de componentes principales.

Conceptos clave

  • Una matriz simétrica es definida positiva cuando xᵀAx > 0 para todo x ≠ 0; equivalentemente, todos sus valores propios (o todos sus pivotes) son positivos. AᵀA es siempre simétrica y semidefinida positiva.
  • Toda matriz m × n tiene una SVD A = UΣVᵀ: V y U tienen columnas ortonormales, y Σ contiene los valores singulares σ₁ ≥ σ₂ ≥ … ≥ 0, las raíces cuadradas de los valores propios de AᵀA.
  • A = σ₁u₁v₁ᵀ + σ₂u₂v₂ᵀ + …. Conservar los k primeros términos da la mejor aproximación de rango k (Eckart–Young); el error en la norma 2 es σₖ₊₁.
  • El número de condición σ₁/σₙ mide cuánto pueden amplificarse los errores en b al resolver Ax = b.
  • PCA: centra los datos (resta las medias de las columnas); entonces los vectores singulares derechos (vectores propios de la matriz de covarianza) son las direcciones principales, y σᵢ² mide la varianza que explica cada una.

Ecuaciones clave

A=UΣVTA=U\Sigma V^T
Descomposición en valores singulares
σi=λi(ATA)\sigma_i=\sqrt{\lambda_i(A^TA)}
Valores singulares
Ak=∑i=1kσiuiviTA_k=\sum_{i=1}^{k}\sigma_i\mathbf{u}_i\mathbf{v}_i^T
Mejor aproximación de rango k
κ=σ1/σn\kappa=\sigma_1/\sigma_n
Número de condición

Ejemplos resueltos paso a paso

Ejemplo 10.1 Halla los valores singulares de A = [[3, 0], [4, 5]].

Ver los pasos
  1. AᵀA = [[25, 20], [20, 25]], valores propios 45 y 5.
  2. σ₁ = √45, σ₂ = √5. Comprobación: σ₁σ₂ = 15 = |det A|.

Respuesta: σ₁ = 6.708, σ₂ = 2.236.

Ejemplo 10.2 ¿Qué parte de la «energía» Σσ² de AᵀA conserva la aproximación de rango 1 de esa A?

Ver los pasos
  1. σ₁²/(σ₁² + σ₂²) = 45/50.

Respuesta: 90%.

Ejemplo 10.3 Una imagen en escala de grises de 1000 × 1000 se guarda como una SVD de rango 50. ¿Cuántos números se guardan, comparado con la imagen completa?

Ver los pasos
  1. Cada término necesita un u (1000), un v (1000) y un σ: 2001 números.
  2. 50 × 2001 = 100050 frente a 1000000.

Respuesta: 100050 números, aproximadamente el 10.0 % del original.

Problemas de práctica

Resuelve cada uno primero en papel y luego comprueba.

  1. ¿Es [[2, 1], [1, 2]] definida positiva?

    Ver la respuesta

    Sí: los valores propios 3 y 1 son positivos.

  2. Halla el número de condición de diag(10, 0.1).

    Ver la respuesta

    10/0.1 = 100.

  3. Unos datos 2D centrados tienen covarianza [[2, 1], [1, 2]]. ¿Qué parte de la varianza explica la primera componente principal?

    Ver la respuesta

    Valores propios 3 y 1: 3/4 = 75 %, a lo largo de (1, 1).

  4. ¿Cuáles son los valores singulares de una matriz ortogonal?

    Ver la respuesta

    Todos iguales a 1.

Ideas erróneas comunes

  • Mito: Los valores singulares son los valores propios de A.

    En realidad: Son las raíces cuadradas de los valores propios de AᵀA, y nunca son negativos.

  • Mito: Solo las matrices cuadradas tienen SVD.

    En realidad: Todas las matrices la tienen.

  • Mito: El PCA funciona sobre datos sin procesar.

    En realidad: Centra primero los datos, y escala las variables medidas en unidades distintas.

Tarjetas de estudio

  • SVD

    A = UΣVᵀ.

  • Valor singular

    σᵢ = √λᵢ(AᵀA).

  • Definida positiva

    xᵀAx > 0, todos los λ > 0.

  • Eckart–Young

    La SVD truncada es la mejor aproximación de rango k.

  • Número de condición

    σ₁/σₙ.

  • PCA

    SVD de los datos centrados.

Cuestionario de autoevaluación

1. Los valores singulares son…

Ver la respuesta

no negativos. Raíces cuadradas de los valores propios de AᵀA, ≥ 0.

2. La mejor aproximación de rango k conserva…

Ver la respuesta

los k mayores términos singulares. Eckart–Young.

3. Un número de condición grande significa…

Ver la respuesta

que los errores pueden amplificarse. Casi singular.

4. En PCA, la primera dirección principal es…

Ver la respuesta

el primer vector singular derecho de los datos centrados. La dirección de mayor varianza.

Cómo ayuda esto a salvar el planeta

El espectrómetro de imágenes EMIT del JPL de la NASA, en la Estación Espacial Internacional, cartografía los minerales de las regiones de origen del polvo de la Tierra. El JPL dice que recoge más de 100000 espectros por segundo, es decir 3.6 × 10⁸ espectros por hora, y que su misión se ha prolongado al menos hasta 2026. El polvo importa para el clima: el JPL señala que los minerales oscuros absorben la luz solar y calientan el aire, mientras que los claros la reflejan y enfrían.

Cada espectro es un vector largo de reflectancias en muchas longitudes de onda, y las longitudes de onda vecinas están muy correlacionadas. El PCA y la SVD comprimen esos datos en unas pocas componentes, separan la señal del ruido y aíslan las firmas de minerales concretos: las mismas matemáticas que la práctica de compresión de imágenes.

En Saving Planets: Observación de la Tierra · Espacio

Fuentes: NASA JPL: EMIT, Earth Surface Mineral Dust Source Investigation · MIT OpenCourseWare 18.06 Linear Algebra (Strang, Spring 2010)

Recursos de aprendizaje gratuitos

Volver arriba

El laboratorio (Computation lab)

El laboratorio es la parte práctica del curso. Aquí es gratuito y solo para practicar, sin créditos.

Cada práctica se hace en una hoja de cálculo gratuita o en papel; las versiones en Python son opcionales. Haz pausas regulares de la pantalla y maneja alfileres y tijeras con cuidado.

Referencia de seguridad: OSHA: Computer workstations eTool

Práctica A: vectores de píxeles y un índice de vegetación

Semana 2

Objetivo: Tratar los píxeles de satélite como vectores: calcular longitudes, ángulos entre píxeles y el NDVI, y ver que la vegetación, el suelo y el agua apuntan en direcciones distintas.

cos⁡θ=p⋅q∥p∥∥q∥,NDVI=NIR−redNIR+red\cos\theta=\frac{\mathbf{p}\cdot\mathbf{q}}{\|\mathbf{p}\|\|\mathbf{q}\|},\quad \text{NDVI}=\frac{\text{NIR}-\text{red}}{\text{NIR}+\text{red}}
Ángulo entre píxeles y NDVI

Procedimiento

  1. En una hoja de cálculo, introduce tres píxeles de ejemplo como filas de reflectancias (verde, rojo, NIR): vegetación (0.10, 0.08, 0.40), suelo (0.20, 0.25, 0.30), agua (0.06, 0.04, 0.02).
  2. Calcula la longitud de cada píxel con SQRT(SUMSQ(…)) y el producto escalar de cada par con SUMPRODUCT.
  3. Halla el ángulo entre cada par con ACOS y DEGREES.
  4. Calcula el NDVI de cada píxel. Opcional: descarga una escena Landsat gratuita y repite con píxeles reales.

Tabla de datos

PíxelVerdeRedNIRNDVI
Vegetación0.100.080.400.667
Suelo0.200.250.300.091
Agua0.060.040.02−0.333

Reflectancias de ejemplo elegidas para parecerse a vegetación, suelo desnudo y agua típicos.

Cálculo de ejemplo

  1. NDVI de la vegetación = (0.40 − 0.08)/(0.40 + 0.08) = 0.667.
  2. NDVI del suelo = 0.05/0.55 = 0.091; NDVI del agua = −0.02/0.06 = −0.333.
  3. Ángulo vegetación–suelo: producto escalar = 0.160, longitudes 0.420 y 0.4387, cos θ = 0.8683, θ = 29.7°.

Seguridad

  • Haz una breve pausa de la pantalla cada media hora.

Alternativa en casa o virtual

  • Toda la práctica se hace en una hoja de cálculo gratuita.

PhET: Vector Addition · USGS: Landsat 9

Vínculo con el planeta: El USGS indica las bandas roja (banda 4) e infrarroja cercana (banda 5) de Landsat 9, y sus datos son gratuitos, así que cualquiera puede cartografiar la vegetación de esta forma. USGS: Landsat 9

Práctica B: trilateración, de cuerdas a un sistema lineal

Semana 3

Objetivo: Hallar una posición a partir de distancias, primero con cuerdas sobre un mapa como en la actividad de GPS.gov y luego convirtiendo las ecuaciones de circunferencias en un sistema lineal y resolviéndolo por eliminación.

2(a2−a1)x+2(b2−b1)y=d12−d22+a22−a12+b22−b122(a_2-a_1)x+2(b_2-b_1)y=d_1^2-d_2^2+a_2^2-a_1^2+b_2^2-b_1^2
Diferencia de dos ecuaciones de circunferencias

Procedimiento

  1. Imprime el mapa de trilateración de GPS.gov (sin escalar) o marca tres «satélites» en papel cuadriculado en (0, 0), (10, 0) y (0, 10) cm.
  2. Convierte los tiempos de señal en distancias con d = ct (c = 299792458 m/s); para la versión en papel usa las distancias de ejemplo de abajo.
  3. Traza circunferencias con una cuerda y un lápiz y marca dónde se cortan.
  4. Resta las ecuaciones de las circunferencias por pares para obtener dos ecuaciones lineales, resuélvelas por eliminación en una hoja de cálculo y compara con tu dibujo.

Tabla de datos

BalizaPosición (cm)Distancia (cm)d²
1(0, 0)5.00025
2(10, 0)8.06265
3(0, 10)6.70845

Distancias de ejemplo para un receptor en (3, 4) cm. La actividad de GPS.gov usa cuatro satélites en un mapa a escala 1:21283839.

Cálculo de ejemplo

  1. Eq. 1 − Eq. 2: 20x − 100 = −40, x = 3. Eq. 1 − Eq. 3: 20y − 100 = −20, y = 4.
  2. El tiempo A de GPS.gov = 0.00505783 s da d = 1516299 m; a 1:21283839 son 7.12 cm en el mapa.

Seguridad

  • Maneja chinchetas y tijeras con cuidado.

Alternativa en casa o virtual

  • La versión con cuerda y papel es la práctica para casa.

GPS.gov: Trilateration activity · Interactive Linear Algebra, 1.2 Row reduction

Vínculo con el planeta: GPS.gov diseñó esta actividad con la NOAA; los receptores reales resuelven el mismo tipo de sistema, con una incógnita de reloj, muchas veces por segundo. GPS.gov: Trilateration (NOAA classroom activity) · GPS.gov: Space segment

Práctica C: matrices como transformaciones del plano

Semana 4

Objetivo: Ver cómo actúan las matrices sobre formas: estirar, cizallar, rotar y reflejar un cuadrado en papel cuadriculado, y confirmar que las áreas se escalan por |det A|.

area(AS)=∣det⁡A∣ area(S)\text{area}(AS)=|\det A|\,\text{area}(S)
Factor de escala de áreas

Procedimiento

  1. Dibuja el cuadrado unidad con vértices (0, 0), (1, 0), (1, 1), (0, 1).
  2. Multiplica cada vértice por A = [[2, 1], [0, 1]] (una cizalla y un estiramiento), por R = [[0, −1], [1, 0]] (una rotación de 90°) y por F = [[1, 0], [0, −1]] (una reflexión).
  3. Dibuja cada imagen y halla su área contando cuadros o con la fórmula del área de Gauss.
  4. Compara cada área con |det| y anota qué aplicaciones invierten la orientación (det < 0). Después comprueba det(AR) = det A det R.

Tabla de datos

MatrizdetÁrea de la imagenOrientación
A = [[2, 1], [0, 1]]22conservada
R (rotación de 90°)11conservada
F (reflexión)−11invertida
AR22conservada

Valores exactos para el cuadrado unidad.

Cálculo de ejemplo

  1. det A = 2 × 1 − 1 × 0 = 2; la imagen es un paralelogramo de base 2 y altura 1.
  2. det(AR) = det A × det R = 2 × 1 = 2.

Seguridad

  • Ninguno más allá del trabajo de escritorio normal.

Alternativa en casa o virtual

  • Solo necesitas papel cuadriculado y un lápiz.

Interactive Linear Algebra, 3.1 Matrix transformations · Interactive Linear Algebra, 4.3 Determinants and volumes

Vínculo con el planeta: Las proyecciones cartográficas y el registro de imágenes usan estas mismas matrices; el determinante mantiene honestas las mediciones de área de bosques e inundaciones. USGS: Landsat 9

Práctica D: ajuste por mínimos cuadrados de datos del nivel del mar

Semana 9

Objetivo: Ajustar una recta, y luego una parábola, a datos del nivel del mar con las ecuaciones normales, y comparar con el ajuste integrado de la hoja de cálculo.

ATAx^=ATb,A=[ 1  t ]A^TA\hat{\mathbf{x}}=A^T\mathbf{b},\quad A=[\,\mathbf{1}\ \ \mathbf{t}\,]
Ecuaciones normales para y = C + Dt

Procedimiento

  1. Introduce los datos de ejemplo de abajo (o descarga el registro de la NASA desde su página Vital Signs): t en años desde 1993 y h en mm.
  2. Construye A con una columna de unos y una columna de t. Calcula AᵀA y Aᵀb con MMULT y TRANSPOSE, y luego x̂ = MINVERSE(AᵀA)·Aᵀb.
  3. Compara con LINEST o con una línea de tendencia del gráfico. Representa los residuos b − Ax̂.
  4. Añade una columna t² y vuelve a ajustar. Un coeficiente de t² positivo significa que la subida se acelera.

Tabla de datos

t (años desde 1993)08162432
h (mm)0224878114

Valores de ejemplo con la forma del registro de la NASA (que indica 114.4 mm de subida desde 1993); usa los datos reales para una respuesta real.

Cálculo de ejemplo

  1. t̄ = 16, h̄ = 52.4; Σ(t − t̄)(h − h̄) = 2272, Σ(t − t̄)² = 640.
  2. D = 2272/640 = 3.55 mm/yr; C = 52.4 − 3.55 × 16 = −4.4 mm.
  3. Los residuos son positivos en ambos extremos y negativos en el medio: una curva ajusta mejor, como sugiere el ritmo creciente que indica la NASA.

Seguridad

  • Haz una breve pausa de la pantalla cada media hora.

Alternativa en casa o virtual

  • Toda la práctica se hace en una hoja de cálculo gratuita.

NASA: Global mean sea level · Interactive Linear Algebra, 6.5 Least squares

Vínculo con el planeta: La NASA informa que el ritmo de subida pasó de 0.20 cm por año en 1993 a 0.44 cm por año ahora, una aceleración que puedes detectar con una columna t². NASA: Global mean sea level (Vital Signs)

Práctica E: un modelo de Markov de materiales y reciclaje

Semana 12

Objetivo: Construir una matriz de transición para el material en uso y el perdido en vertedero o incineración, usar potencias de matrices para seguirlo en el tiempo y ver cómo una tasa de reciclaje mayor mantiene el material en uso más tiempo.

xk+1=Pxk,P=[1−f(1−r)0f(1−r)1]\mathbf{x}_{k+1}=P\mathbf{x}_k,\quad P=\begin{bmatrix}1-f(1-r)&0\\f(1-r)&1\end{bmatrix}
Cadena de dos estados: en uso, perdido

Procedimiento

  1. Elige una tasa de descarte de ejemplo f = 0.2 por año y una tasa de recuperación r = 0.321 (la proporción de reciclaje más compostaje de la EPA en 2018).
  2. Construye P en una hoja de cálculo y empieza con x₀ = (1, 0): todo el material en uso.
  3. Usa MMULT para avanzar 10 años y compara con la fórmula (1 − f(1 − r))ᵏ para la proporción que sigue en uso.
  4. Repite con r = 0.5 y compara los años medios en uso, 1/(f(1 − r)).

Tabla de datos

Tasa de recuperación rProbabilidad de seguir en usoEn uso tras 10 añosAños medios en uso
0.3210.86420.2327.36
0.500.90000.34910.0

Modelo de ejemplo: la tasa de descarte f = 0.2 es un valor de ejemplo; r = 0.321 es la proporción de la EPA de 2018.

Cálculo de ejemplo

  1. Probabilidad de seguir = 1 − 0.2 × 0.679 = 0.8642; 0.8642¹⁰ = 0.232.
  2. Años medios en uso = 1/0.1358 = 7.36; con r = 0.5 es 1/0.1 = 10.0 años.
  3. λ = 1 corresponde al estado «perdido», que absorbe todo a largo plazo; λ = 0.8642 fija la rapidez.

Seguridad

  • Haz una breve pausa de la pantalla cada media hora.

Alternativa en casa o virtual

  • Toda la práctica se hace en una hoja de cálculo gratuita.

Interactive Linear Algebra, 5.6 Stochastic matrices · EPA: National overview of materials and waste

Vínculo con el planeta: La EPA informa de una tasa de reciclaje y compostaje del 32.1 % en 2018 y un 50 % a vertedero; modelos como este muestran cuánto vale subir esa tasa. U.S. EPA: National overview, facts and figures on materials, wastes and recycling

Práctica F: método de la potencia y compresión de imágenes con SVD

Semana 14

Objetivo: Hallar un valor propio dominante por multiplicación repetida y luego comprimir una pequeña imagen en escala de grises con una SVD truncada y ver con qué pocos números se conserva la imagen.

xk+1=Axk∥Axk∥,Ak=∑i=1kσiuiviT\mathbf{x}_{k+1}=\frac{A\mathbf{x}_k}{\|A\mathbf{x}_k\|},\quad A_k=\sum_{i=1}^k\sigma_i\mathbf{u}_i\mathbf{v}_i^T
Método de la potencia y SVD truncada

Procedimiento

  1. En una hoja de cálculo, empieza con x₀ = (1, 0) y multiplica por A = [[2, 1], [1, 2]] diez veces, reescalando cada vez; observa cómo el cociente de longitudes se acerca a 3.
  2. Crea una malla 8 × 8 de niveles de gris (0–255) para un dibujo sencillo, o carga una foto en Python (gratuito) con NumPy.
  3. Calcula la SVD (numpy.linalg.svd) y reconstruye la imagen con k = 1, 2 y 4 términos.
  4. Anota el almacenamiento k(8 + 8 + 1) y la parte de Σσ² conservada, y juzga cuándo la imagen se ve bien.

Tabla de datos

kNúmeros almacenados (imagen 8 × 8)Parte de 64
11726.6%
23453.1%
468106.3%

Recuentos exactos de almacenamiento. Para una imagen 8 × 8, k = 4 ya almacena más que los 64 números originales; el ahorro aparece en imágenes grandes (unidad 10: el rango 50 de una imagen 1000 × 1000 conserva el 10.0 %).

Cálculo de ejemplo

  1. Método de la potencia: (1, 0) → (2, 1) → (5, 4) → (14, 13); ‖(14, 13)‖/‖(5, 4)‖ = 19.10/6.403 = 2.98, cerca de λ = 3.
  2. Almacenamiento para k = 1: 8 + 8 + 1 = 17 números, 17/64 = 26.6 %.

Seguridad

  • Haz una breve pausa de la pantalla cada media hora.

Alternativa en casa o virtual

  • Las versiones en hoja de cálculo y en Python gratuito son ambas prácticas para casa.

MIT OpenCourseWare 18.06 Linear Algebra (Strang, Spring 2010) · Interactive Linear Algebra, 5.1 Eigenvectors

Vínculo con el planeta: EMIT, del JPL de la NASA, devuelve más de 100000 espectros por segundo; comprimir y separar esos datos usa exactamente estas herramientas. NASA JPL: EMIT, Earth Surface Mineral Dust Source Investigation

Volver arriba

Series de práctica para exámenes

Cronométrate: unos 10 minutos por problema. Luego compara con las soluciones completas.

Ejercicios de práctica para el parcial (unidades 1–5)

  1. Halla el ángulo entre (1, 0, 1) y (0, 1, 1).

    Ver la solución completa
    1. Producto escalar = 1; longitudes √2 cada una; cos θ = 1/2.

    Respuesta: 60°.

  2. Resuelve x + 2y = 4, 3x − y = 5.

    Ver la solución completa
    1. R₂ − 3R₁: −7y = −7.

    Respuesta: x = 2, y = 1.

  3. Invierte [[3, 2], [1, 1]].

    Ver la solución completa
    1. det = 1.

    Respuesta: [[1, −2], [−1, 3]].

  4. Una economía de un solo sector usa 0.4 de su propia producción. ¿Qué producción satisface una demanda de 120?

    Ver la solución completa
    1. x = 120/0.6.

    Respuesta: 200.

  5. Una matriz 6 × 9 tiene rango 4. Halla dim N(A) y dim N(Aᵀ).

    Ver la solución completa
    1. n − r y m − r.

    Respuesta: 5 y 2.

Serie de práctica final (todas las unidades)

  1. Ajusta y = C + Dt a (0, 0), (1, 2), (2, 2) por mínimos cuadrados.

    Ver la solución completa
    1. AᵀA = [[3, 3], [3, 5]], Aᵀb = (4, 6); det 6.

    Respuesta: C = 0.333, D = 1.

  2. Halla det [[2, 1, 0], [1, 2, 1], [0, 1, 2]].

    Ver la solución completa
    1. 2(4 − 1) − 1(2 − 0) + 0.

    Respuesta: 4.

  3. Halla los valores propios de [[6, 2], [2, 3]].

    Ver la solución completa
    1. Traza 9, det 14: λ² − 9λ + 14 = 0.

    Respuesta: 7 y 2.

  4. Halla el estado estacionario de la matriz de Markov [[0.9, 0.2], [0.1, 0.8]].

    Ver la solución completa
    1. 0.1q₁ = 0.2q₂.

    Respuesta: (0.667, 0.333).

  5. Halla los valores singulares de [[2, 0], [0, −3]].

    Ver la solución completa
    1. AᵀA = diag(4, 9).

    Respuesta: 3 y 2.

  6. ¿Cuántos números almacena una SVD de rango 20 de una imagen 500 × 800, y qué parte de la imagen completa es eso?

    Ver la solución completa
    1. 20(500 + 800 + 1).

    Respuesta: 26020, aproximadamente el 6.5 %.

Volver arriba

Hoja de fórmulas

∥v∥=v⋅v, cos⁡θ=v⋅w∥v∥∥w∥\|v\|=\sqrt{v\cdot v},\ \cos\theta=\frac{v\cdot w}{\|v\|\|w\|}
Ax=b,  A=LUA\mathbf{x}=\mathbf{b},\ \ A=LU
[abcd]−1=1ad−bc[d−b−ca]\begin{bmatrix}a&b\\c&d\end{bmatrix}^{-1}=\frac{1}{ad-bc}\begin{bmatrix}d&-b\\-c&a\end{bmatrix}
(AB)−1=B−1A−1, (AB)T=BTAT(AB)^{-1}=B^{-1}A^{-1},\ (AB)^T=B^TA^T
x=(I−A)−1d\mathbf{x}=(I-A)^{-1}\mathbf{d}
r+dim⁡N(A)=nr+\dim N(A)=n
p=a⋅ba⋅aa,  ATAx^=ATb\mathbf{p}=\frac{a\cdot b}{a\cdot a}a,\ \ A^TA\hat{x}=A^Tb
A=QR, QTQ=IA=QR,\ Q^TQ=I
det⁡(AB)=det⁡Adet⁡B, xi=det⁡Bi/det⁡A\det(AB)=\det A\det B,\ x_i=\det B_i/\det A
det⁡(A−λI)=0, ∑λ=tr⁡A, ∏λ=det⁡A\det(A-\lambda I)=0,\ \sum\lambda=\operatorname{tr}A,\ \prod\lambda=\det A
A=PDP−1, Ak=PDkP−1, A=QΛQTA=PDP^{-1},\ A^k=PD^kP^{-1},\ A=Q\Lambda Q^T
A=UΣVT, σi=λi(ATA), κ=σ1/σnA=U\Sigma V^T,\ \sigma_i=\sqrt{\lambda_i(A^TA)},\ \kappa=\sigma_1/\sigma_n

Constantes

ConstanteValor
Velocidad de la luz (usada para las distancias GPS)299,792,458 m/s
1 milla1.609344 km
1 pie0.3048 m
1 pulgada25.4 mm
Flops para resolver un sistema n × n por eliminaciónunas n³/3 multiplicaciones–sumas (Strang)
Matriz identidad IUnos en la diagonal, ceros en el resto: AI = IA = A

Fuente: NIST Guide to the SI, Appendix B.9: conversion factors

Volver arriba

Textos y clases gratuitos

Usa estos libros de texto, clases y herramientas gratuitos junto con el curso.

Volver arriba

Glosario

Base
Un conjunto linealmente independiente que genera un espacio.
Espacio columna
Todas las combinaciones Ax de las columnas de A.
Número de condición
σ₁/σₙ, cuánto pueden crecer los errores al resolver Ax = b.
Determinante
Factor de escala de volumen con signo de A; es cero exactamente cuando A es singular.
Diagonalización
A = PDP⁻¹ con los vectores propios en P.
Dimensión
Número de vectores de cualquier base.
Valor propio
λ con Av = λv para algún v no nulo.
Factorización LU
A = LU a partir de la eliminación.
Mínimos cuadrados
El x̂ que minimiza ‖b − Ax‖.
Inversa de Leontief
(I − A)⁻¹, producción total necesaria por unidad de demanda final.
Independencia lineal
Solo la combinación nula da el vector cero.
Matriz de Markov
Entradas no negativas, columnas que suman 1.
Espacio nulo
Todos los x con Ax = 0.
Ortonormal
Vectores unitarios ortogonales entre sí.
Pivote
Primera entrada no nula de una fila en forma escalonada.
Definida positiva
Simétrica con xᵀAx > 0 para x ≠ 0.
Proyección
Punto más cercano a b en un subespacio.
Rango
Número de pivotes; dimensión del espacio columna.
Valor singular
σᵢ = √λᵢ(AᵀA).
Espacio generado
Todas las combinaciones lineales de un conjunto de vectores.
Subespacio
Un conjunto que contiene el 0 y es cerrado bajo la suma y el escalado.

Volver arriba

Rúbrica de autoevaluación

No se otorga nota ni crédito. Úsala para valorar con honestidad dónde estás en cada unidad.

NivelCómo se ve
4 · DominioElige la herramienta adecuada (eliminación, proyección, valores propios o SVD), calcula a mano los casos pequeños, comprueba las respuestas (Ax = b, traza y determinante) y explica el significado con dibujos y en contexto.
3 · CompetenteResuelve los problemas de práctica con algún error aritmético ocasional.
2 · En desarrolloSigue los ejemplos resueltos pero confunde filas y columnas o los cuatro subespacios.
1 · InicialManipula símbolos sin comprobar tamaños ni respuestas; necesita repasar la unidad.

Ponderación sugerida si quieres una nota

PartePeso
Problemas de práctica25%
Prácticas20%
Serie de práctica de mitad de curso25%
Serie de práctica final30%

Apunta al nivel 3 o superior en cada unidad antes de pasar al siguiente curso.

Volver arriba

Ayudar / Donar