short time fourier transform
3 次查看(过去 30 天)
显示 更早的评论
I would like to perform a short time fourier transform (STFT) to a synthetic data series. I can compute the fourier transform by:
fs = 40;
t = 0:(1 / fs):4;
y1 = [sin(2 * pi * 5 * t(t <= 2)), sin(2 * pi * 10 * t(t > 2))];
subplot(311);
plot(t,y1,'k');
Fy1 = abs(ifft(y1));
N = numel(t);
idx = 1:numel(Fy1) / 2; % # Indices of half the spectrum
f = fs * (0:(N - 1)) / N; % # Actual frequencies
subplot(312);
plot(f(idx),2*Fy1(idx),'k');
As a third subplot, I would now like to compute the STFT otherwise known as the windowed fourier transform, which will show hoe the frequency of the signal varies in time. How can this be done?
0 个评论
回答(3 个)
Wayne King
2012-12-13
Why not just use spectrogram if you have the Signal Processing Toolbox?
You can output the STFT and/or the short-time periodograms (PSD estimates)
Wayne King
2012-12-13
Then just use the 'yaxis' option
t = 0:0.001:1-0.001;
x=chirp(t,0,1,150);
spectrogram(x,'yaxis')
If you actually output arguments from spectrogram, you have much more control over the plot. See the help.
0 个评论
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Time-Frequency Analysis 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!