FFT for non-periodic signal
显示 更早的评论
Assume that I have a sequence of N different samples, so my signal is non-periodic. Is it ok I take the FFT for the whole sequence at a time after using a Hanning window which length is also N? Is the result different if I use part of my signal and 'make' it periodic?
采纳的回答
Wayne King
2013-12-19
1 个投票
It is fine to use fft() on non-periodic data. The DFT (discrete Fourier transform) works just fine on non-periodic data. However, it is implicit in the DFT that the signal is extended periodically.
Look at the equation for the inverse DFT, the signal is N-periodic just as the DFT is N-periodic.
2 个评论
piyush arya
2016-6-21
Shouldn't the signal be 0 at endpoints. My signal is o at the beginning at non-zero at the end.
Chris Turnes
2016-6-21
No; there's nothing about the DFT or inverse DFT that imposes a restriction on the sample values. Like Wayne said, the DFT essentially implicitly assumes that the signal is extended periodically; that doesn't mean the samples on the ends have to match in value.
更多回答（2 个）
Wayne King
2013-12-19
编辑：Wayne King
2013-12-19
Ngoc, you have to compensate for the effect of the window. Do you have the Signal Processing Toolbox? If so, let periodogram() with the 'power' option take care of that for you
Fs = 1000;
t = 0:1/Fs:1-1/Fs;
x = 2*cos(2*pi*100*t);
[Sxx,F] = periodogram(x,flattopwin(length(x)),length(x),Fs,'power');
plot(F,Sxx)
Note the power is correctly reported as 2, which is A^2/2
If you have non-periodic data,then periodogram() with the 'power' option defaults to 'psd', which is more appropriate.
In other words, the 'power' option is appropriate when you have a line spectrum consisting of sinusoidal components.
1 个评论
Ngoc
2013-12-23
Oh thank you Wayne. Now I need to plot the spectrum results. Do you know how I can get this attached color map? I already have k for x axis, f for y axis, and now I need to plot the matrix which is the result of the fft.
类别
在 帮助中心 和 File Exchange 中查找有关 Discrete Fourier and Cosine Transforms 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
