He intentado diferentes códigos de transformada de Fourier en ondas sinusoidales, y todos producen un espectro distribuido con una resonancia en la frecuencia de la señal cuando teóricamente deberían mostrar una sola barra.
La frecuencia de muestreo tiene poco efecto (10 kHz aquí), sin embargo, el número de ciclos:
Un ciclo:
100 ciclos:
100000 ciclos:
Parece que la transformada de Fourier converge solo por un número infinito de ciclos, ¿por qué es eso? ¿No debería una ventana de tiempo de exactamente un ciclo traer los mismos resultados que la de N ciclos?
Aplicación: Esto es por curiosidad y también porque quiero obtener cuánto la respuesta escalonada de un sistema de primer orden excitará la resonancia de un ensamblaje mecánico. Por lo tanto, necesito una transformación precisa de Fourier de la respuesta ... en la que ya no confío. ¿Qué podría hacer para mejorar la precisión en base al caso de "onda sinusoidal"?
PD: Estas capturas de pantalla particulares se basan en el código aquí .
fuente
Respuestas:
Este es un artefacto de ventanas.
El código vinculado rellena una señal de muestra de 10,000 con ceros para que la longitud sea una potencia de dos.
Tenga en cuenta que en el código anterior, la FFT se toma con el tamaño de FFT,
NFFT
que es la siguiente potencia de 2 más grande que la longitud de la señal (en este caso, 16,384). De la documentación de Mathworksfft()
:Esto significa que en realidad no está tomando una FFT de una 'onda sinusoidal pura', está tomando la FFT de una onda sinusoidal con una señal plana después de ella.
Esto es equivalente a tomar la FFT de una onda sinusoidal multiplicada por una función de ventana cuadrada. El espectro FFT es entonces la convolución del espectro de frecuencia de onda sinusoidal (una función de impulso) con el espectro de frecuencia de onda cuadrada (sinc (f)).
Si cambia
L = 16,384
para que no haya relleno de cero de la señal, observará unaperfect
FFT.Otras palabras clave de búsqueda: "Fuga espectral", "Función de ventana", "Ventana de Hamming".
Editar: Limpié parte del material que escribí sobre este tema en la universidad, que entra en detalles sustancialmente más. Lo publiqué en mi blog .
fuente