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.
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.
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.
| Semana | Unidad | Temas | Lectura |
|---|---|---|---|
| 1 | Vectores: combinaciones lineales, longitudes y productos escalares | Vectores, combinaciones lineales, espacio generado, longitud | MIT 18.06 L1; ILA 2.1–2.2 |
| 2 | Vectores: combinaciones lineales, longitudes y productos escalares | Productos escalares, ángulos, ortogonalidad; Práctica A; inicio de la eliminación | ILA 6.1 |
| 3 | Sistemas de ecuaciones lineales y eliminación gaussiana | Eliminación, formas escalonadas, conjuntos solución; Práctica B | MIT 18.06 L2; ILA 1.1–1.3 |
| 4 | Álgebra de matrices: multiplicación, inversas y factorización LU | Multiplicación de matrices, inversas, A = LU; Práctica C | MIT 18.06 L3–4; ILA 3.4–3.5 |
| 5 | Modelos lineales: economía de insumo–producto (Leontief), redes y balanceo | Modelos de insumo–producto, redes, balanceo | Hefferon Topics; EPA USEEIO |
| 6 | Espacios vectoriales, subespacios, independencia, base, rango y los cuatro subespacios fundamentales | Subespacios, espacio columna, espacio nulo, Ax = 0 | MIT 18.06 L5–7; ILA 2.6 |
| 7 | Espacios vectoriales, subespacios, independencia, base, rango y los cuatro subespacios fundamentales | Independencia, base, dimensión, rango, cuatro subespacios; práctica de examen parcial | MIT 18.06 L8–10; ILA 2.5, 2.7, 2.9 |
| 8 | Ortogonalidad, proyecciones, mínimos cuadrados y Gram–Schmidt | Ortogonalidad y proyecciones | MIT 18.06 L14–15; ILA 6.3 |
| 9 | Ortogonalidad, proyecciones, mínimos cuadrados y Gram–Schmidt | Mínimos cuadrados, Gram–Schmidt, A = QR; Práctica D | MIT 18.06 L16–17; ILA 6.4–6.5 |
| 10 | Determinantes: propiedades, cofactores, regla de Cramer y volumen | Determinantes, cofactores, regla de Cramer, volumen | MIT 18.06 L18–20; ILA 4.1–4.3 |
| 11 | Valores propios y vectores propios | Valores propios, vectores propios, método de la potencia, modos | MIT 18.06 L21; ILA 5.1–5.2, 5.5 |
| 12 | Diagonalización, potencias de matrices, cadenas de Markov y sistemas dinámicos lineales | Diagonalización, cadenas de Markov, sistemas dinámicos; Práctica E | MIT 18.06 L22–25; ILA 5.4, 5.6 |
| 13 | La descomposición en valores singulares, el análisis de componentes principales y las matrices definidas positivas | Matrices definidas positivas y la SVD | MIT 18.06 L27–29 |
| 14 | La descomposición en valores singulares, el análisis de componentes principales y las matrices definidas positivas | Aproximación de rango bajo, PCA, número de condición; Práctica F | MIT 18.06 L29; sitio del libro de Strang |
| 15 | Repaso | Repaso y serie de práctica final | Todas las unidades |
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.
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
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
- ‖v‖ = √(3² + 4²) = √25 = 5.
- 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
- v·w = 3 + 8 = 11.
- ‖v‖ = 5, ‖w‖ = √5 = 2.236.
- 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
- 2(1, 0, 1) + 3(0, 1, 1) = (2, 3, 5).
- 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.
Halla ‖(1, 2, 2)‖.
Ver la respuesta
√(1 + 4 + 4) = 3.
Halla (1, 2, 3)·(4, −5, 6).
Ver la respuesta
4 − 10 + 18 = 12.
¿Son ortogonales (2, −1) y (1, 2)?
Ver la respuesta
Yes: 2 − 2 = 0.
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
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
- Interactive Linear Algebra, 2.1 Vectors: Texto gratuito con figuras interactivas.
- Interactive Linear Algebra, 6.1 Dot products and orthogonality: Longitudes, ángulos y ortogonalidad.
- MIT OpenCourseWare 18.06 Linear Algebra (Strang, Spring 2010): Clase 1: la geometría de las ecuaciones lineales.
- PhET: Vector Addition: Interactivo.
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
Ejemplos resueltos paso a paso
Ejemplo 2.1 Resuelve x + y + z = 6, 2x + y − z = 1, x − y + 2z = 5.
Ver los pasos
- R₂ − 2R₁: −y − 3z = −11. R₃ − R₁: −2y + z = −1.
- R₃ − 2R₂: 7z = 21, así que z = 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
- Ecuación 1 menos ecuación 2: x² + y² − [(x − 10)² + y²] = 25 − 65, así que 20x − 100 = −40.
- Ecuación 1 menos ecuación 3: 20y − 100 = −20.
- 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
- R₂ − 2R₁ da 0 = 0 en el primer sistema: un pivote, una variable libre.
- En el segundo, R₂ − 2R₁ da 0 = 1.
Respuesta: Infinitas; ninguna.
Problemas de práctica
Resuelve cada uno primero en papel y luego comprueba.
Resuelve x + y = 5, x − y = 1.
Ver la respuesta
Sumando se obtiene 2x = 6: x = 3, y = 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.
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.
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
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
- Interactive Linear Algebra, 1.1 Systems of linear equations: Sección de texto gratuito.
- Interactive Linear Algebra, 1.2 Row reduction: Formas escalonadas paso a paso.
- MIT OpenCourseWare 18.06 Linear Algebra (Strang, Spring 2010): Clase 2: eliminación con matrices.
- Jim Hefferon, Linear Algebra (free text): Capítulo uno: sistemas lineales.
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
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
- AB: fila 1 · columnas = (2, 1); fila 2 = (4, 3).
- 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
- ad − bc = 24 − 14 = 10.
- 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
- Multiplicador ℓ₂₁ = 6/2 = 3.
- 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.
¿Qué tamaño tiene una matriz 3 × 2 por una 2 × 4?
Ver la respuesta
3 × 4.
Invierte [[2, 1], [5, 3]].
Ver la respuesta
det = 1, inversa [[3, −1], [−5, 2]].
¿Aproximadamente cuántas operaciones requiere la eliminación para n = 1000?
Ver la respuesta
n³/3 ≈ 3.33 × 10⁸.
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
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
- Interactive Linear Algebra, 3.4 Matrix multiplication: Sección de texto gratuito.
- Interactive Linear Algebra, 3.5 Matrix inverses: Inversas y Gauss–Jordan.
- MIT OpenCourseWare 18.06 Linear Algebra (Strang, Spring 2010): Clases 3–4: multiplicación, inversas y A = LU.
- MIT OpenCourseWare 18.06SC Linear Algebra (self-paced, Fall 2011): Vídeos de clases prácticas sobre LU.
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
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
- I − A = [[0.8, −0.3], [−0.1, 0.6]], det = 0.48 − 0.03 = 0.45.
- (I − A)⁻¹ = (1/0.45)[[0.6, 0.3], [0.1, 0.8]].
- 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
- 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
- C: x = z. H: 4x = 2w. O: 2y = 2z + w.
- 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.
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.
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.
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.
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
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
- U.S. EPA: US Environmentally-Extended Input-Output (USEEIO) models: Los modelos insumo–producto ampliados ambientalmente de la EPA.
- Jim Hefferon, Linear Algebra (free text): Tema: análisis insumo–producto de Leontief.
- MIT OpenCourseWare 18.06 Linear Algebra (Strang, Spring 2010): Clases y listas de problemas de aplicaciones.
- Interactive Linear Algebra, 1.2 Row reduction: Resolver los sistemas.
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
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
- R₂ − 2R₁ = 0; R₃ − R₁ = [0, −2, −2]. Dos pivotes, así que r = 2.
- Ax = 0: x + z = 0 y −2y − 2z = 0, con z libre.
- 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
- 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.
- 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
- C(A) y C(Aᵀ): r = 3.
- 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.
¿Es la recta y = 2x + 1 un subespacio de ℝ²?
Ver la respuesta
No: no contiene (0, 0).
Una matriz 4 × 4 tiene un espacio nulo de dimensión 2. ¿Cuál es su rango?
Ver la respuesta
4 − 2 = 2.
¿Cuál es la dimensión del espacio de matrices 2 × 2?
Ver la respuesta
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
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
- Interactive Linear Algebra, 2.6 Subspaces: Sección de texto gratuito.
- Interactive Linear Algebra, 2.5 Linear independence: Pruebas de independencia.
- Interactive Linear Algebra, 2.9 The rank theorem: Rango y nulidad.
- MIT OpenCourseWare 18.06 Linear Algebra (Strang, Spring 2010): Clases 5–10: espacios, rango y los cuatro subespacios.
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
Ejemplos resueltos paso a paso
Ejemplo 6.1 Proyecta b = (2, 3) sobre a = (1, 1).
Ver los pasos
- a·b = 5, a·a = 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
- A tiene filas (1, 0), (1, 1), (1, 2); AᵀA = [[3, 3], [3, 5]], Aᵀb = (7, 10).
- 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
- q₁ = (1, 1, 0)/√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.
Proyecta (3, 4) sobre (1, 0).
Ver la respuesta
(3, 0).
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.
¿Es ortogonal Q = [[0, 1], [1, 0]]?
Ver la respuesta
Yes: QᵀQ = I.
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
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
- Interactive Linear Algebra, 6.3 Orthogonal projection: Sección de texto gratuito.
- Interactive Linear Algebra, 6.5 The method of least squares: Mínimos cuadrados con imágenes.
- Interactive Linear Algebra, 6.4 Orthogonal sets: Gram–Schmidt.
- MIT OpenCourseWare 18.06 Linear Algebra (Strang, Spring 2010): Clases 14–17: ortogonalidad, proyecciones, mínimos cuadrados, Gram–Schmidt.
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
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 − 3 × 0) = 1.
- − 2 × (0 × 1 − 3 × 2) = −2 × (−6) = 12.
- + 0 × (…) = 0.
Respuesta: det = 13.
Ejemplo 7.2 Resuelve 2x + y = 5, x + 3y = 10 con la regla de Cramer.
Ver los pasos
- det A = 6 − 1 = 5.
- x = det[[5, 1], [10, 3]]/5 = (15 − 10)/5.
- 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
- Á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.
¿det [[3, 1], [2, 4]]?
Ver la respuesta
12 − 2 = 10.
¿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.
Si det A = 4 para una A 3 × 3, ¿cuánto vale det(2A)?
Ver la respuesta
2³ × 4 = 32.
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
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
- Interactive Linear Algebra, 4.1 Determinants: definition: Sección de texto gratuito.
- Interactive Linear Algebra, 4.3 Determinants and volumes: Significado geométrico.
- MIT OpenCourseWare 18.06 Linear Algebra (Strang, Spring 2010): Clases 18–20: determinantes, cofactores, regla de Cramer.
- Jim Hefferon, Linear Algebra (free text): Capítulo cuatro: determinantes.
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
Ejemplos resueltos paso a paso
Ejemplo 8.1 Halla los valores y vectores propios de A = [[2, 1], [1, 2]].
Ver los pasos
- det(A − λI) = (2 − λ)² − 1 = λ² − 4λ + 3 = 0.
- λ = 3: A − 3I = [[−1, 1], [1, −1]] da v = (1, 1). λ = 1: v = (1, −1).
- 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
- λ² − 7λ + 10 = 0 (traza 7, det 10).
- (λ − 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
- Valores propios de [[2, −1], [−1, 1]]: λ = (3 ± √5)/2 = 0.382 y 2.618.
- ω₁² = 100 × 0.382 = 38.2 s⁻², ω₁ = 6.18 rad/s.
- 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.
¿Valores propios de [[5, 0], [0, −2]]?
Ver la respuesta
5 y −2 (diagonal).
Una matriz 3 × 3 tiene valores propios 1, 2, 3. Halla su traza y su determinante.
Ver la respuesta
Traza 6, det 6.
¿Valores propios de la rotación de 90° [[0, −1], [1, 0]]?
Ver la respuesta
λ² + 1 = 0: λ = ±i.
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
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
- Interactive Linear Algebra, 5.1 Eigenvalues and eigenvectors: Sección de texto gratuito.
- Interactive Linear Algebra, 5.2 The characteristic polynomial: Hallar valores propios.
- Interactive Linear Algebra, 5.5 Complex eigenvalues: Rotación y escalado.
- MIT OpenCourseWare 18.06 Linear Algebra (Strang, Spring 2010): Clase 21: valores propios y vectores propios.
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
Ejemplos resueltos paso a paso
Ejemplo 9.1 Usa la diagonalización para hallar A¹⁰ para A = [[2, 1], [1, 2]].
Ver los pasos
- De la unidad 8: λ = 3, 1 con v = (1, 1), (1, −1). P = [[1, 1], [1, −1]], P⁻¹ = ½[[1, 1], [1, −1]].
- 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
- (P − I)q = 0: −0.2q₁ + 0.3q₂ = 0, así que q₁ = 1.5q₂.
- 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
- x₁ = (0.8, 0.2); x₂ = (0.64 + 0.06, 0.16 + 0.14).
- 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.
¿Es [[0.5, 0.2], [0.5, 0.8]] una matriz de Markov?
Ver la respuesta
Sí: entradas no negativas, columnas que suman 1.
Halla su estado estacionario.
Ver la respuesta
0.5q₁ = 0.2q₂, así que q = (2/7, 5/7) = (0.286, 0.714).
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ᵏ.
¿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
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
- Interactive Linear Algebra, 5.4 Diagonalization: Sección de texto gratuito.
- Interactive Linear Algebra, 5.6 Stochastic matrices: Cadenas de Markov y estados estacionarios.
- MIT OpenCourseWare 18.06 Linear Algebra (Strang, Spring 2010): Clases 22–25: diagonalización, matrices de Markov, matrices simétricas.
- MIT OpenCourseWare 18.06SC Linear Algebra (self-paced, Fall 2011): Clases prácticas sobre matrices de Markov.
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
Ejemplos resueltos paso a paso
Ejemplo 10.1 Halla los valores singulares de A = [[3, 0], [4, 5]].
Ver los pasos
- AᵀA = [[25, 20], [20, 25]], valores propios 45 y 5.
- σ₁ = √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
- σ₁²/(σ₁² + σ₂²) = 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
- Cada término necesita un u (1000), un v (1000) y un σ: 2001 números.
- 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.
¿Es [[2, 1], [1, 2]] definida positiva?
Ver la respuesta
Sí: los valores propios 3 y 1 son positivos.
Halla el número de condición de diag(10, 0.1).
Ver la respuesta
10/0.1 = 100.
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).
¿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
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
- MIT OpenCourseWare 18.06 Linear Algebra (Strang, Spring 2010): Clases 27–29: matrices definidas positivas y la SVD.
- MIT OpenCourseWare 18.06SC Linear Algebra (self-paced, Fall 2011): Sesión sobre la descomposición en valores singulares.
- Gilbert Strang, Introduction to Linear Algebra (book site): Notas de Strang sobre la SVD y los datos.
- Jim Hefferon, Linear Algebra (free text): Temas sobre matrices ortogonales y proyecciones.
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.
Procedimiento
- 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).
- Calcula la longitud de cada píxel con SQRT(SUMSQ(…)) y el producto escalar de cada par con SUMPRODUCT.
- Halla el ángulo entre cada par con ACOS y DEGREES.
- Calcula el NDVI de cada píxel. Opcional: descarga una escena Landsat gratuita y repite con píxeles reales.
Tabla de datos
| Píxel | Verde | Red | NIR | NDVI |
|---|---|---|---|---|
| Vegetación | 0.10 | 0.08 | 0.40 | 0.667 |
| Suelo | 0.20 | 0.25 | 0.30 | 0.091 |
| Agua | 0.06 | 0.04 | 0.02 | −0.333 |
Reflectancias de ejemplo elegidas para parecerse a vegetación, suelo desnudo y agua típicos.
Cálculo de ejemplo
- NDVI de la vegetación = (0.40 − 0.08)/(0.40 + 0.08) = 0.667.
- NDVI del suelo = 0.05/0.55 = 0.091; NDVI del agua = −0.02/0.06 = −0.333.
- Á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.
Procedimiento
- 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.
- 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.
- Traza circunferencias con una cuerda y un lápiz y marca dónde se cortan.
- 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
| Baliza | Posición (cm) | Distancia (cm) | d² |
|---|---|---|---|
| 1 | (0, 0) | 5.000 | 25 |
| 2 | (10, 0) | 8.062 | 65 |
| 3 | (0, 10) | 6.708 | 45 |
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
- Eq. 1 − Eq. 2: 20x − 100 = −40, x = 3. Eq. 1 − Eq. 3: 20y − 100 = −20, y = 4.
- 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|.
Procedimiento
- Dibuja el cuadrado unidad con vértices (0, 0), (1, 0), (1, 1), (0, 1).
- 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).
- Dibuja cada imagen y halla su área contando cuadros o con la fórmula del área de Gauss.
- 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
| Matriz | det | Área de la imagen | Orientación |
|---|---|---|---|
| A = [[2, 1], [0, 1]] | 2 | 2 | conservada |
| R (rotación de 90°) | 1 | 1 | conservada |
| F (reflexión) | −1 | 1 | invertida |
| AR | 2 | 2 | conservada |
Valores exactos para el cuadrado unidad.
Cálculo de ejemplo
- det A = 2 × 1 − 1 × 0 = 2; la imagen es un paralelogramo de base 2 y altura 1.
- 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.
Procedimiento
- 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.
- 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.
- Compara con LINEST o con una línea de tendencia del gráfico. Representa los residuos b − Ax̂.
- 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) | 0 | 8 | 16 | 24 | 32 |
|---|---|---|---|---|---|
| h (mm) | 0 | 22 | 48 | 78 | 114 |
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
- t̄ = 16, h̄ = 52.4; Σ(t − t̄)(h − h̄) = 2272, Σ(t − t̄)² = 640.
- D = 2272/640 = 3.55 mm/yr; C = 52.4 − 3.55 × 16 = −4.4 mm.
- 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.
Procedimiento
- 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).
- Construye P en una hoja de cálculo y empieza con x₀ = (1, 0): todo el material en uso.
- 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.
- Repite con r = 0.5 y compara los años medios en uso, 1/(f(1 − r)).
Tabla de datos
| Tasa de recuperación r | Probabilidad de seguir en uso | En uso tras 10 años | Años medios en uso |
|---|---|---|---|
| 0.321 | 0.8642 | 0.232 | 7.36 |
| 0.50 | 0.9000 | 0.349 | 10.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
- Probabilidad de seguir = 1 − 0.2 × 0.679 = 0.8642; 0.8642¹⁰ = 0.232.
- Años medios en uso = 1/0.1358 = 7.36; con r = 0.5 es 1/0.1 = 10.0 años.
- λ = 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.
Procedimiento
- 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.
- Crea una malla 8 × 8 de niveles de gris (0–255) para un dibujo sencillo, o carga una foto en Python (gratuito) con NumPy.
- Calcula la SVD (numpy.linalg.svd) y reconstruye la imagen con k = 1, 2 y 4 términos.
- Anota el almacenamiento k(8 + 8 + 1) y la parte de Σσ² conservada, y juzga cuándo la imagen se ve bien.
Tabla de datos
| k | Números almacenados (imagen 8 × 8) | Parte de 64 |
|---|---|---|
| 1 | 17 | 26.6% |
| 2 | 34 | 53.1% |
| 4 | 68 | 106.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
- 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.
- 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
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)
Halla el ángulo entre (1, 0, 1) y (0, 1, 1).
Ver la solución completa
- Producto escalar = 1; longitudes √2 cada una; cos θ = 1/2.
Respuesta: 60°.
Resuelve x + 2y = 4, 3x − y = 5.
Ver la solución completa
- R₂ − 3R₁: −7y = −7.
Respuesta: x = 2, y = 1.
Invierte [[3, 2], [1, 1]].
Ver la solución completa
- det = 1.
Respuesta: [[1, −2], [−1, 3]].
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
- x = 120/0.6.
Respuesta: 200.
Una matriz 6 × 9 tiene rango 4. Halla dim N(A) y dim N(Aᵀ).
Ver la solución completa
- n − r y m − r.
Respuesta: 5 y 2.
Serie de práctica final (todas las unidades)
Ajusta y = C + Dt a (0, 0), (1, 2), (2, 2) por mínimos cuadrados.
Ver la solución completa
- AᵀA = [[3, 3], [3, 5]], Aᵀb = (4, 6); det 6.
Respuesta: C = 0.333, D = 1.
Halla det [[2, 1, 0], [1, 2, 1], [0, 1, 2]].
Ver la solución completa
- 2(4 − 1) − 1(2 − 0) + 0.
Respuesta: 4.
Halla los valores propios de [[6, 2], [2, 3]].
Ver la solución completa
- Traza 9, det 14: λ² − 9λ + 14 = 0.
Respuesta: 7 y 2.
Halla el estado estacionario de la matriz de Markov [[0.9, 0.2], [0.1, 0.8]].
Ver la solución completa
- 0.1q₁ = 0.2q₂.
Respuesta: (0.667, 0.333).
Halla los valores singulares de [[2, 0], [0, −3]].
Ver la solución completa
- AᵀA = diag(4, 9).
Respuesta: 3 y 2.
¿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
- 20(500 + 800 + 1).
Respuesta: 26020, aproximadamente el 6.5 %.
Hoja de fórmulas
Constantes
| Constante | Valor |
|---|---|
| Velocidad de la luz (usada para las distancias GPS) | 299,792,458 m/s |
| 1 milla | 1.609344 km |
| 1 pie | 0.3048 m |
| 1 pulgada | 25.4 mm |
| Flops para resolver un sistema n × n por eliminación | unas n³/3 multiplicaciones–sumas (Strang) |
| Matriz identidad I | Unos en la diagonal, ceros en el resto: AI = IA = A |
Fuente: NIST Guide to the SI, Appendix B.9: conversion factors
Textos y clases gratuitos
Usa estos libros de texto, clases y herramientas gratuitos junto con el curso.
- MIT OpenCourseWare 18.06 Linear Algebra (Strang, Spring 2010): Las clases en vídeo completas, listas de problemas y exámenes de Strang.
- MIT OpenCourseWare 18.06SC Linear Algebra (self-paced, Fall 2011): Versión a tu ritmo con vídeos de clases prácticas y problemas resueltos.
- Interactive Linear Algebra (Margalit and Rabinoff, Georgia Tech): Libro de texto interactivo gratuito usado en Georgia Tech.
- Jim Hefferon, Linear Algebra (free text): Libro de texto completo y gratuito con respuestas.
- Gilbert Strang, Introduction to Linear Algebra (book site): El sitio del libro de Strang con notas y figuras.
- PhET: Vector Addition: Vectores y componentes interactivos.
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.
Rúbrica de autoevaluación
No se otorga nota ni crédito. Úsala para valorar con honestidad dónde estás en cada unidad.
| Nivel | Cómo se ve |
|---|---|
| 4 · Dominio | Elige 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 · Competente | Resuelve los problemas de práctica con algún error aritmético ocasional. |
| 2 · En desarrollo | Sigue los ejemplos resueltos pero confunde filas y columnas o los cuatro subespacios. |
| 1 · Inicial | Manipula símbolos sin comprobar tamaños ni respuestas; necesita repasar la unidad. |
Ponderación sugerida si quieres una nota
| Parte | Peso |
|---|---|
| Problemas de práctica | 25% |
| Prácticas | 20% |
| Serie de práctica de mitad de curso | 25% |
| Serie de práctica final | 30% |
Apunta al nivel 3 o superior en cada unidad antes de pasar al siguiente curso.



