在Matlab上用FFT计算和绘制信号的频谱
提问于 2017-07-26 14:19:05
在Matlab上使用FFT计算信号的频谱时，我遗漏了一些东西。我的代码：
%% compute the spectrum of the data (data(t))
L = length(time); % length of the sample
NFFT = 2^(nextpow2(L)-1); % Next power of 2 from length of y
Y = fft(data,NFFT);%/NFFT;%L;
Fs = 1/(mean(time(2:end)-time(1:end-1))); % compute the sampling frequency
f = Fs/2*linspace(0,1,NFFT/2+1);
loglog(f,2*abs(Y(1:NFFT/2+1)))
title('Single-Sided Amplitude Spectrum of My Data')
xlabel('Frequency (Hz)')
ylabel('|Y(f)|')你能告诉我我哪里搞砸了吗？
我尝试使用相同信号的这两个采样(相同的采样频率；在0-10和0-100两个不同的时间范围内)来检查算法是否有效：
fs=1000;
time10 = [0:1/fs:10];
time100 = [0:1/fs:100];
data10 = sin(2*pi*0.23 .*time10)+cos(2*pi*12 .*time10);
data100 = sin(2*pi*0.23 .*time100)+cos(2*pi*12 .*time100);我想这两种光谱应该是相互支持的，但他们没有...如下所示：https://www.dropbox.com/s/wfols9o409pr94u/FFT_spectrum_StackOverflow.png?dl=0 https://www.dropbox.com/s/a8vmzwto6x4130w/FFT_spectrum_StackOverflow.fig?dl=0
谢谢
回答 1
Stack Overflow用户
回答已采纳
发布于 2017-07-27 12:21:02
好消息是，你对频谱的计算本身并没有什么问题。
问题是，通过查看不同长度的样本，您实际上是在查看两个完全不同的样本。在时域中，它们可以看作是无限长的正弦与不同长度的矩形窗口相乘的结果。
在频域中，无限长连续时间正弦信号的频谱与矩形窗口的频谱进行卷积。随着窗口长度的不同，这些窗口的相应频谱具有不同的宽度(较长的矩形窗口的频谱较窄)。因此，无限长的正弦信号频谱中的尖峰将分布在不同的带宽上。这就是你所看到的。
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/45329651
复制相关文章
点击加载更多
腾讯云开发者
