Hola,
Estuvimos viendo que la parte de FFT quedó más complicada de lo que esperábamos. Para darle un cierre al tema, subí un zip en la sección Laboratorios que contiene un script y un par de imágenes. En el script está bastante comentado el paso a paso, como para que puedan entender lo que se hace.
Algunas aclaraciones:
- Interpolación: una de las señales quedó muestreada con la mitad de fs que las otras, por eso hay unos factores de interpolación para las señales (en el tiempo), con el objetivo de corregir esto. Además, la ventana temporal es pequeña como para tener buena resolución en frecuencia, entonces lo que hacemos es interpolar en frecuencia, que es equivalente a agrandar la ventana temporal colocando ceros al final de la señal. Supongo que esto se explicó en teórico y lo habrán visto en cursos anteriores, pero pueden consultar si no queda claro. Si estamos en las hipótesis que nos permiten interpolar en frecuencia, podemos interpolar con tantos puntos como necesitemos. En este caso, con 50 puntos el espectro queda suave.
- Elección de la parte de la señal a transformar: no queremos ver la serie de Fourier asociada a la periodicidad de la onda cuadrada, por lo que es necesario cortar la señal. ¿Con qué parte nos quedamos? A priori habíamos dicho que tomaran sólo medio periodo. Procesando los datos, observamos que si hacemos eso perdemos información del comportamiento de la punta en frecuencias altas, que en realidad es la información de interés, así que concluimos que es mejor cortar la señal en un periodo completo. Es necesario que ese periodo termine en 0, para poder interpolar en frecuencia sin agregar nueva información (nuevos saltos).
- Valor de continua: si miran la ecuación de la punta de prueba, en continua debería comportarse como un divisor x10 ideal, por lo que sería esperable que la componente espectral en continua estuviera superpuesta. Se eliminó el factor de energía que podría estar molestando, normalizando todos los espectros por la energía del espectro del caso compensado, pero aún se aprecia una diferencia.
Saludos,
Florencia