Diagonalización y modos normales

Desacoplar sistemas lineales y leer modos independientes. Teoría, fórmulas, ejemplos resueltos y gráficas interactivas.

Que es

Diagonalizar una matriz A significa encontrar una matriz P, formada por sus autovectores como columnas, tal que P^-1 A P = D, donde D es una matriz diagonal con los autovalores de A en la diagonal. Esta operación cambia la base en la que se describe el sistema: en la base de autovectores, la transformación que antes mezclaba todas las coordenadas ahora actúa de forma independiente sobre cada una.

Cuando la matriz A describe un sistema físico acoplado, sus autovectores son precisamente los modos normales: combinaciones específicas de las coordenadas originales que se comportan como si fueran variables independientes, cada una evolucionando con su propio autovalor asociado sin interferir con las demás.

Por qué aparece en física

Aparece porque muchos sistemas físicos de interés (osciladores acoplados, sólidos vibrantes, circuitos con varios elementos, sistemas cuánticos de varios niveles) están descritos por ecuaciones donde cada coordenada depende de todas las demás. Resolver ese sistema acoplado directamente es complicado, pero si se diagonaliza la matriz que lo describe, el problema se convierte en varios problemas independientes, uno por cada modo normal, cada uno tan simple como un oscilador aislado.

También aparece al simplificar el cálculo de potencias de una matriz (A^n = P D^n P^-1, mucho más simple que multiplicar A por si misma n veces) y en la resolución de sistemas de ecuaciones diferenciales lineales acoplados, donde diagonalizar desacopla las ecuaciones y permite resolver cada una por separado.

Idea visual o geométrica

La idea visual es esta: imagina un sistema de coordenadas girado y posiblemente deformado, alineado exactamente con las direcciones propias de la matriz. En esa base especial, la transformación deja de mezclar direcciones: cada eje se estira o se comprime de forma independiente, según su autovalor correspondiente. Es como encontrar los "ejes naturales" del problema, aquellos en los que la física se ve más simple.

En un sistema de masas y resortes, cada modo normal es una forma concreta de mover todas las masas a la vez, de manera que, si el sistema empieza exactamente en esa configuración, oscila con una única frecuencia sin cambiar de forma, aunque la amplitud crezca o decrezca.

Lenguaje simbólico mínimo

Si A tiene autovectores v1, ..., vn linealmente independientes con autovalores λ1, ..., λn, entonces P = [v1 | v2 | ... | vn] y D = diag(λ1, ..., λn) satisfacen AP = PD, es decir, A = P D P^-1. Cualquier vector x se puede escribir en la base de autovectores como x = c1 v1 + ... + cn vn, y la acción de A sobre x se reduce a Ax = c1 λ1 v1 + ... + cn λn vn.

Para un sistema de ecuaciones diferenciales lineales dx/dt = Ax, el cambio a la base de autovectores separa el sistema en n ecuaciones independientes dy_i/dt = λ_i y_i, cada una con solución exponencial y_i(t) = y_i(0) e^(λ_i t).

Que significa el resultado

El resultado significa que el comportamiento completo y aparentemente complicado del sistema acoplado se puede escribir como una superposición de comportamientos simples independientes: cada modo normal evoluciona a su propio ritmo, determinado por su autovalor, sin intercambiar energía o informacion con los demás modos mientras el sistema permanezca lineal.

Físicamente, cualquier estado inicial del sistema puede descomponerse en una combinación de modos normales, y la evolución temporal completa se obtiene evolucionando cada modo por separado y sumando el resultado al final, en la base original.

Cuando no debes usarlo

No debes intentar diagonalizar una matriz que no tiene suficientes autovectores independientes (matrices defectuosas): en ese caso, la diagonalización simple no existe y se necesita la forma canónica de Jordan, que incluye bloques no diagonales.

Tampoco tiene sentido buscar modos normales si el sistema es no lineal: la superposición de modos independientes solo funciona porque la ecuación es lineal, y una vez que aparecen términos cuadráticos o de orden superior, los modos dejan de evolucionar de forma independiente y empiezan a intercambiar energía entre si.

Resumen operativo

Encuentra los autovalores y autovectores de la matriz que describe el sistema acoplado, construye la matriz de cambio de base P con esos autovectores, y usa la relación A = P D P^-1 para convertir el problema original en un conjunto de problemas independientes en la base de modos normales. Resuelve cada modo por separado y, si es necesario, transforma el resultado de vuelta a las coordenadas originales.

Procedimiento

Antes de calcular

Comprueba primero que la matriz tiene suficientes autovectores linealmente independientes para diagonalizarse: si es simétrica, esto está garantizado; si no lo es, verifica la multiplicidad de cada autovalor antes de continuar.

Decide también si el objetivo final es simplificar un cálculo (potencias de la matriz), resolver un sistema de ecuaciones diferenciales, o interpretar modos normales de un sistema físico acoplado, ya que esto orienta como interpretar el resultado final.

Paso 1: identifica el objeto matemático

El objeto central es la matriz de cambio de base P, formada por los autovectores como columnas, y la matriz diagonal D formada por los autovalores en el mismo orden. Antes de construirlas, asegúrate de tener ya calculados todos los autovalores y sus autovectores asociados.

Ordena los autovalores de forma consistente: la posición i de la diagonal de D debe corresponder exactamente a la columna i de P.

Paso 2: elige la operación

Si el objetivo es diagonalizar, construye P y D, y comprueba la relación AP = PD antes de invertir P. Si el objetivo es resolver un sistema dx/dt = Ax, cambia de variable a y = P^-1 x, resuelve cada ecuación independiente dy_i/dt = λ_i y_i, y después deshaz el cambio con x = Py.

Si el objetivo es calcular una potencia A^n, usa A^n = P D^n P^-1, calculando D^n simplemente elevando cada autovalor a la potencia n, mucho más simple que multiplicar A repetidamente.

Paso 3: conserva unidades y signos

Mantén la correspondencia exacta entre cada columna de P y el autovalor correspondiente en D: intercambiar el orden en uno sin cambiarlo en el otro invalida todo el cálculo posterior.

Si los autovectores representan magnitudes físicas con unidades, verifica que la matriz P mezcla correctamente esas unidades y que el resultado final, tras deshacer el cambio de base, recupera las unidades esperadas.

Paso 4: calcula

Calcula P^-1 explícitamente si necesitas expresar un vector en la base de modos normales; si la matriz es simétrica, aprovecha que P puede elegirse ortogonal, de forma que P^-1 = P^T, lo que simplifica enormemente el cálculo.

Resuelve el sistema desacoplado en la base de autovectores, uno a uno, antes de recombinar el resultado en la base original.

Paso 5: interpreta el resultado

Interpreta cada componente en la base de autovectores como la amplitud de un modo normal independiente. Un autovalor real positivo indica crecimiento exponencial de ese modo, uno negativo indica decaimiento, y uno complejo indica oscilación combinada con crecimiento o decaimiento según su parte real.

Al volver a la base original, describe el estado del sistema como una superposición de esos modos, indicando que fracción de cada modo está presente en el estado inicial dado.

Checklist de autocontrol

Comprueba que la correspondencia entre columnas de P y elementos de D es correcta, que AP = PD se cumple, que el número de autovectores independientes coincide con el tamaño de la matriz, y que puedes describir en palabras que representa físicamente cada modo normal obtenido.

Ejemplo guiado

Situación física

Un sistema de dos masas iguales conectadas por resortes tiene matriz de acoplamiento A = [[2, −1], [−1, 2]] (ver el leaf de autovalores y autovectores para el cálculo detallado). Se conocen sus autovalores λ1 = 1, λ2 = 3 con autovectores v1 = (1, 1), v2 = (1, −1). Se quiere describir la evolución completa del sistema a partir de una condición inicial arbitraria.

Traducción matemática

Diagonalizar A permite escribir cualquier estado inicial x(0) como combinación de los modos normales: x(0) = c1 v1 + c2 v2. Cada modo evoluciona de forma independiente según su propio autovalor, y la solución completa es la superposición de ambos.

Calculo paso a paso

Si la condición inicial es x(0) = (2, 0), se resuelve c1 (1,1) + c2 (1,−1) = (2, 0), lo que da c1 = 1, c2 = 1. Cada modo oscila con frecuencia proporcional a la raíz de su autovalor: el modo v1 con frecuencia proporcional a √1 y el modo v2 con frecuencia proporcional a √3. La solución completa es x(t) = 1·v1·cos(√1 t) + 1·v2·cos(√3 t) (asumiendo velocidad inicial nula).

Lectura del resultado

El estado inicial (2, 0), que no es ni puramente simétrico ni puramente antisimétrico, se descompone en partes iguales de ambos modos normales. Con el tiempo, ambos modos oscilan a frecuencias distintas, por lo que el movimiento conjunto deja de tener una forma fija y se vuelve una combinación cambiante de ambos patrones.

Error común evitado

Se evita intentar resolver las dos ecuaciones acopladas directamente sin cambiar de base, lo que requeriría manipular ambas coordenadas simultáneamente en cada paso. Al diagonalizar, cada modo se resuelve como un oscilador simple independiente.

Ejemplo inverso

Resultado dado

Un sistema tiene un único modo normal excitado, con autovalor negativo λ = −2.

Que pregunta física permite responder

Un autovalor negativo en este contexto (para una ecuación de evolución temporal de primer orden) indica que ese modo decae exponencialmente con el tiempo, no que oscila. Esto permite anticipar que, sin importar la amplitud inicial de ese modo, su contribución al estado del sistema desaparece a largo plazo, dejando dominar a los modos con autovalores de decaimiento más lento o con autovalores positivos si los hay.

Errores frecuentes

Error 1: desordenar la correspondencia entre columnas de P y elementos de D

Parece razonable escribir los autovalores y autovectores en cualquier orden, siempre que estén todos presentes. Falla porque la columna i de P debe corresponder exactamente al autovalor en la posición i de D; invertir el orden en uno sin el otro produce una matriz que no satisface AP = PD. Se detecta comprobando esa igualdad directamente. Se corrige emparejando explícitamente cada autovector con su autovalor antes de construir las matrices. Mini-ejemplo: si λ1 = 1 con v1 y λ2 = 3 con v2, escribir P = [v2 | v1] exige también D = diag(3, 1), no diag(1, 3).

Error 2: creer que todo sistema acoplado tiene modos normales sencillos

Parece razonable esperar que cualquier sistema lineal se diagonalice sin problema. Falla si la matriz es defectuosa, es decir, si no tiene tantos autovectores independientes como su tamaño. Se detecta contando los autovectores obtenidos y comparándolos con el tamaño de la matriz. Se corrige reconociendo que esos casos necesitan un tratamiento distinto (forma de Jordan) en vez de forzar una diagonalización que no existe. Mini-ejemplo: una matriz con un autovalor doble puede tener solo un autovector independiente, no dos.

Error 3: superponer modos como si el sistema fuera no lineal

Parece razonable aplicar la misma superposición de modos normales a cualquier sistema. Falla si el sistema tiene términos no lineales, donde los modos dejan de evolucionar de forma independiente y empiezan a intercambiar energía entre si. Se detecta comprobando si las ecuaciones del sistema son realmente lineales en las coordenadas usadas. Se corrige linealizando explícitamente antes de diagonalizar, y declarando el rango de validez de esa aproximación. Mini-ejemplo: un péndulo doble con oscilaciones grandes no tiene modos normales exactos, solo aproximados para pequeñas amplitudes.

Error 4: olvidar deshacer el cambio de base al final

Parece razonable dejar la solución expresada en la base de autovectores. Falla porque el problema físico original suele estar planteado en las coordenadas originales, y el resultado en la base de modos normales no siempre tiene una interpretación física directa componente a componente. Se detecta preguntando si el resultado final responde la pregunta original del problema. Se corrige aplicando x = Py para volver a las coordenadas físicas antes de dar la respuesta. Mini-ejemplo: la amplitud de un modo normal no es directamente la posición de una masa concreta, sino una combinación de varias.

Error 5: confundir autovalores complejos con un error de cálculo

Parece razonable pensar que un autovalor complejo indica un fallo en el planteamiento. Falla porque en sistemas con rotación u oscilación, los autovalores complejos son físicamente esperados y describen comportamiento oscilatorio, no un error. Se detecta comprobando si la matriz del sistema es simétrica (que garantiza autovalores reales) o no. Se corrige interpretando la parte real como crecimiento o decaimiento y la parte imaginaria como frecuencia, en vez de descartar el resultado. Mini-ejemplo: un sistema amortiguado y oscilante tiene autovalores complejos con parte real negativa, lo que describe oscilación decreciente.

Regla rápida de detección

Si AP no coincide exactamente con PD, si se cuenta menos autovectores independientes que dimensiones tiene la matriz, o si se olvida deshacer el cambio de base antes de interpretar el resultado físicamente, el error probablemente está en el planteamiento de la diagonalización, no en la aritmética final.

Aplicaciones físicas

1. Osciladores acoplados y moléculas

Aparece al describir las vibraciones de un conjunto de masas conectadas o de los átomos de una molécula. Se diagonaliza la matriz de acoplamiento para encontrar los modos normales de vibración. El resultado da las frecuencias propias y las formas de vibración independientes. El error peligroso es intentar excitar un modo normal ignorando que en la práctica un estado inicial arbitrario contiene una mezcla de varios modos. Como control final, el alumno debe explicar en que combinación de modos se descompone un estado inicial dado.

2. Vibraciones de sólidos y redes cristalinas

Aparece al estudiar como vibra una red de átomos en un sólido. Se diagonaliza la matriz que describe las fuerzas entre átomos vecinos. El resultado da los modos normales de la red, relacionados con los fonones en la descripción cuántica. El error peligroso es tratar cada átomo de forma aislada, ignorando que su movimiento está acoplado al de sus vecinos. Como control final, el alumno debe explicar por qué la frecuencia de cada modo depende de la longitud de onda asociada.

3. Sistemas cuánticos de varios niveles

Aparece al diagonalizar el operador hamiltoniano de un sistema descrito por una matriz finita. El resultado da los niveles de energía posibles (autovalores) y los estados propios asociados (autovectores). El error peligroso es olvidar que un estado cuántico general es una superposición de estados propios, no un único estado propio. Como control final, el alumno debe explicar que significa medir la energía de un sistema que se encuentra en una superposición de estados propios.

4. Estabilidad de puntos de equilibrio

Aparece al linealizar un sistema dinámico cerca de un punto de equilibrio y diagonalizar la matriz jacobiana resultante. El resultado clasifica el equilibrio en estable, inestable o mixto según el signo de los autovalores. El error peligroso es ignorar un autovalor positivo pequeño pensando que su efecto es despreciable a corto plazo. Como control final, el alumno debe explicar que ocurre con una perturbación a lo largo de cada modo normal del sistema linealizado.

5. Circuitos eléctricos acoplados

Aparece al describir varias mallas de un circuito conectadas entre si mediante inductancias o capacidades mutuas. Se diagonaliza la matriz de acoplamiento para encontrar los modos de oscilación independientes de corriente o carga. El resultado simplifica el análisis de la respuesta transitoria del circuito. El error peligroso es analizar cada malla por separado sin considerar el acoplamiento mutuo con las demás. Como control final, el alumno debe explicar como se relaciona la frecuencia de cada modo normal con los valores de inductancia y capacidad del circuito.

Practica

Nivel 1: reconocimiento

Identifica, dada una matriz con sus autovalores y autovectores ya calculados, cuales serian las columnas de P y los elementos de D antes de construirlas explícitamente.

Ejercicio: para una matriz 2×2 con autovalores 2 y 5 y sus autovectores correspondientes, escribe P y D respetando la correspondencia correcta.

Nivel 2: aplicación directa

Construye P y D para una matriz sencilla y verifica que AP = PD sustituyendo directamente.

Ejercicio: diagonaliza la matriz de un sistema de dos masas acopladas y comprueba el resultado con la verificación AP = PD.

Nivel 3: elección de herramienta

Decide si el problema requiere diagonalizar para resolver un sistema de ecuaciones diferenciales, para calcular una potencia de la matriz, o para interpretar modos normales físicos, y adapta el procedimiento a ese objetivo.

Ejercicio: plantea como usarías la diagonalización para calcular A^10 sin multiplicar la matriz diez veces.

Nivel 4: interpretación física

Interpreta una condición inicial dada como una combinación de modos normales y describe como evoluciona cada modo por separado.

Ejercicio: dado un estado inicial de un sistema de dos masas acopladas, descompónlo en los dos modos normales conocidos y describe la evolución de cada uno.

Nivel 5: error intencionado

Analiza un cálculo donde se ha desordenado la correspondencia entre autovalores y autovectores al construir P y D, o donde se ha olvidado deshacer el cambio de base. Detecta el error, corrígelo y explica la consecuencia física.

Ejercicio: inventa un error plausible al emparejar autovalores con autovectores, muestra por qué parece razonable, corrígelo y señala que resultado incorrecto habría producido.