Error getting Fourier transform of a 1d signal
2 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
Anshuman Pal
el 2 de Mzo. de 2020
Comentada: Star Strider
el 5 de Mzo. de 2020
Hello,
I have an almost-sinusoidal signal (see attached file), which I believe is composed of a mixture of maybe 2 (or 3) Fourier modes:
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/274746/image.png)
I wanted to verify this by plotting its Fourier transform. I expected a flat signal with a bunch of peaks at the signal harmonics. However, using fft() on the signal gives me a 2d plot that makes absolutely no sense to me. Can someone please help me to understand what is going on?
4 comentarios
Star Strider
el 2 de Mzo. de 2020
Remember that you need to plot the absolute value (use the abs function) of the fresult fo the fft.
Also, since the signal has a D-C offset, subtract the mean of the signal before taking the fft in order to see the peaks more clearly:
FTsignal = fft(signal - mean(signal))/numel(signal);
assuming that ‘signal’ is a vector.
Respuesta aceptada
Star Strider
el 3 de Mzo. de 2020
The easiest way to determine the location of the peaks is to use the findpeaks function (with the appropriate name-value pairs to isolate your peaks of interest). I cannot guarantee that the frequencies will be integer multiples of each other, however findpeaks is the best way to identify them.
Improve the frequency resolution by increasing the length of the fft itself. Specify an appropriately large value for ‘n’ (the second argument) to increase the frequency resolution to whatever you want.
6 comentarios
Star Strider
el 5 de Mzo. de 2020
As always, my pleasure!
I experimented with them until I got the result I wanted.
With respect to peak prominences, it is possible to get those directly from findpeaks by requesting extra outputs, then choosing the minimum that returns what you want. Here, I just experimented.
The same for histfit. I experimented with different bin numbers until I got what appeared to me to be an acceptable result.
Más respuestas (0)
Ver también
Categorías
Más información sobre Spectral Measurements en Help Center y File Exchange.
Productos
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!