Laboratorio 1 · Cátedra Acha · Departamento de Física, FCEyN, UBA

Herramientas de análisis

Respuestas · Señales periódicas: cruces, ajuste senoidal y FFT

Respuesta corta y desarrollo completo de las preguntas de la simulación.

Cada pregunta viene con una respuesta corta y con el desarrollo completo. Están plegadas a propósito: intentá contestarla vos primero — con la simulación abierta al lado, moviendo los deslizadores — y recién después abrí el desplegable. Leer la respuesta sin haber intentado no deja nada.

Volver a la simulación · Todas las respuestas

(a) Duplicá la duración del registro: ¿cuánto mejora cada método? ¿Y si duplicás fs?
Respuesta corta
Duplicar T mejora muchísimo los tres (la FFT afina Δf a la mitad y el ajuste mejora como T−3/2). Duplicar fs mejora poco: sólo como √2, por tener más muestras.
El desarrollo

La asimetría es el punto central del panel. La frecuencia se determina por cuántos ciclos completos observaste, no por cuán finamente muestreaste cada ciclo:

  • FFT: Δf = 1/Ttotal, punto. Duplicar T parte la resolución al medio; duplicar fs no la toca en absoluto (sólo sube el techo de Nyquist).
  • Ajuste senoidal: σf ∝ σ/(A·√N·T). Duplicar T da un factor 2 por T y √2 por N: casi 3. Duplicar fs da sólo √2.
  • Cruces con ajuste lineal: el número de cruces crece con T, y σb ∝ 1/√(Σ(k−k̄)²) ∝ N−3/2: también domina T.

Consecuencia experimental directa: si querés un período preciso, dejá corriendo la medición. Cambiar la placa de adquisición por una más rápida es tirar plata. Es la misma estructura que «medir N períodos del péndulo» y que «ampliar el rango de L»: la información está en la extensión de la variable independiente.

(b) Subí el ruido hasta 1.0: ¿cuál de los tres se rompe primero y por qué?
Respuesta corta
Los cruces por cero. Con σ comparable a la amplitud aparecen cruces espurios cerca del nivel medio y el método detecta períodos que no existen.
El desarrollo

El método de cruces mira sólo lo que pasa en un entorno del nivel medio, que es justamente donde la señal tiene menos amplitud y el ruido pesa más. Con A/σ ≈ 1 puede haber tres o cuatro cruces por período real, y la serie tk vs k deja de ser una recta.

La FFT aguanta bastante más: el ruido blanco se reparte en todas las frecuencias mientras la señal se concentra en una, así que la relación señal/ruido en el espectro mejora como √N respecto de la del tiempo. Es el mismo principio de la detección sincrónica.

El ajuste senoidal es el más robusto de los tres, porque usa todos los puntos y toda la forma del modelo. Su límite es distinto: necesita que el modelo sea correcto. Si la señal tiene armónicos, deriva de línea de base o amplitud que decae, el ajuste de una senoidal pura sesga la frecuencia, y la FFT —que no supone forma— pasa a ser preferible.

(c) Poné f = 95 Hz y fs = 100 Hz: ¿qué vas a medir? ¿Cómo te darías cuenta midiendo?
Respuesta corta
Vas a medir 5 Hz. Y te das cuenta cambiando fs: la frecuencia verdadera no se mueve, el alias sí.
El desarrollo

El alias está en |f − k·fs| = |95 − 100| = 5 Hz. Los datos muestreados son literalmente compatibles con las dos señales: no hay procesamiento posterior que pueda distinguirlas, porque la información se perdió en el instante del muestreo.

Una comprobación práctica consiste en: repetir con otra frecuencia de muestreo. Si medís de nuevo a 110 Hz, una señal real de 5 Hz sigue apareciendo en 5 Hz, mientras que el alias de 95 Hz salta a |95 − 110| = 15 Hz. Cualquier pico que se mueve al cambiar fs es un alias.

La prevención de fábrica es un filtro antialias analógico antes del conversor, que corta todo lo que esté por encima de fs/2. Las placas de adquisición decentes lo traen; las improvisadas, no. Y con Tracker el problema existe igual: el obturador de la cámara es tu frecuencia de muestreo, y filmar a 30 fps una vibración de 50 Hz produce exactamente este fenómeno — es por lo que en el cine las ruedas de las carretas giran al revés.

(d) El nivel de continua, ¿afecta a alguno de los tres métodos?
Respuesta corta
A los tres los deja indiferentes, si los implementás bien: cruces por el nivel medio y no por cero, ajuste con offset libre, y FFT restando la media antes de transformar.
El desarrollo

Pero cada uno falla de una forma distinta si te olvidás:

  • Cruces: si buscás cruces por cero y la señal tiene continua mayor que la amplitud, puede no haber cruces y el método deja de producir una estimación. Con continua intermedia, los cruces se corren hacia una zona de menor pendiente y la incerteza empeora.
  • Ajuste: si no ponés el término constante, el ajuste intenta compensar el offset deformando amplitud y fase, y sesga la frecuencia. Es exactamente el mismo error que forzar Tamb = 0 en el enfriamiento.
  • FFT: la continua aparece como un pico enorme en f = 0 cuyas colas —por el efecto de ventana— contaminan las frecuencias bajas y pueden tapar el pico real si f es chica. Por eso se resta la media antes de transformar, que es lo que hace esta simulación.

Ninguno de los tres errores es sutil una vez que lo viste. Los tres son frecuentísimos en informes de primer año.

Si alguna respuesta no te resulta clara o te parece discutible, llevala a clase. Varias preguntas admiten más de un enfoque, y en todos los casos conviene justificar el razonamiento.

Volver a la simulación · Índice general