🧮 Computational Mathematics

Inicio · Parte 13 — Teoría de la información, señales y series

279 — Espectro y densidad espectral

avanzado clase 19 de 20 4 horas demostración power_spectrum

La potencia va con el cuadrado de la amplitud: doble amplitud es cuádruple potencia.

Fórmulas

PSD[k] = |X[k]|² / N
potencia relativa = PSD[k] / Σ PSD
amplitud ×2  ⟹  potencia ×4

Desarrollo

La densidad espectral de potencia describe cómo se reparte la energía de una señal entre sus frecuencias. Se obtiene del cuadrado de la magnitud de la transformada, y a diferencia del espectro de amplitud descarta la fase, quedándose solo con cuánta energía hay en cada banda.

La relación cuadrática entre amplitud y potencia tiene consecuencias que sorprenden al leer un espectro. Una componente con el doble de amplitud aporta cuatro veces más potencia, y por tanto domina el reparto mucho más de lo que sugiere el espectro de amplitud. Las componentes pequeñas quedan aplastadas, y por eso la PSD suele graficarse en decibelios, que es escala logarítmica.

La PSD es la herramienta habitual para caracterizar ruido y para diagnosticar sistemas. El ruido blanco tiene potencia uniforme en todas las frecuencias, el ruido rosa decae como 1/f, y muchos procesos naturales exhiben ese comportamiento. Comparar la PSD medida con la esperada revela componentes anómalas: un pico a 50 Hz delata acoplamiento con la red eléctrica.

En aprendizaje automático, la PSD por bandas es un descriptor clásico. En señales biomédicas las bandas alfa, beta y theta del EEG se definen así; en audio, el reparto de potencia por bandas es la base de los coeficientes MFCC que alimentaron el reconocimiento de voz durante décadas.

Ejemplo trabajado

Dos componentes con amplitudes 2 y 1, y su reparto de potencia.

128 muestras
  componente de 5 Hz  con amplitud 2,0
  componente de 20 Hz con amplitud 1,0

bins dominantes: [5, 20]                             ✓

potencia relativa 5 Hz:   80,0 %
potencia relativa 20 Hz:  20,0 %
razón de potencias: 4,0

Las amplitudes están en razón 2:1
Las potencias están en razón 4:1                     ✓

Aunque la segunda componente tiene la mitad de amplitud,
solo aporta la cuarta parte de la energía.

Qué calcula el laboratorio

Densidad espectral de potencia y reparto de la energía.

python classes/part-13-teoria-de-la-informacion-senales-y-series/279-espectro-y-densidad-espectral/lab.py
compmath run 279

Salidas del laboratorio (9)

Muestra de la ejecución real

{
  "muestras": 128,
  "componentes": {
    "5 Hz": 2.0,
    "20 Hz": 1.0
  },
  "bins_dominantes": [
    5,
    20
  ],
  "potencia_relativa_5Hz_%": 80.0,
  "potencia_relativa_20Hz_%": 20.0,
  "razon_de_potencias": 4.0
}

Errores comunes

Dónde se usa

Caracterización de ruido, análisis de EEG por bandas, MFCC en reconocimiento de voz y diagnóstico de maquinaria por vibraciones.

Idea rectora de la parte

Nyquist fija la frecuencia mínima de muestreo; por debajo hay aliasing irreversible.

Error a evitar

Calcular log(0) sin epsilon de estabilidad.

Conexión con IA

La función de pérdida de casi todo clasificador es entropía cruzada; el VAE optimiza un ELBO con un término KL; las CNN son convoluciones aprendidas.

Bibliografía de la clase

Archivos de la clase