yametw 发表于 2009-6-3 18:38

fft画功率谱密度图有问题

以下是两个不同频率的正弦激励生成功率谱，一个是345，一个是350，为什么样子明显不一样？
clear;figure;
a=sin((1:1000)*0.0005*345*2*pi);
NFFT=1000;
c=fftshift(fft(a));
b=abs(c.^2)/1000^2;
ff=10*log10(b(500:1000));
ff1 = 2000/2*linspace(0,1,NFFT/2+1);
plot(ff1,ff)

figure;
a=sin((1:1000)*0.0005*350*2*pi);
NFFT=1000;
c=fftshift(fft(a));
b=abs(c.^2)/1000^2;
ff=10*log10(b(500:1000));
ff1 = 2000/2*linspace(0,1,NFFT/2+1);
plot(ff1,ff)

无水1324 发表于 2009-6-3 21:20

不用fftshift呢

yametw 发表于 2009-6-3 21:29

不用fftshift也是一样

ChaChing 发表于 2009-6-3 22:43

估计这是相同问题!
http://forum.vibunion.com/forum/viewthread.php?tid=82766&page=1#pid426435

yametw 发表于 2009-6-3 23:07

不是那个问题，因为我虽然只有500步的序列，计算fft时也只用了500步，并没有补零

ChaChing 发表于 2009-6-4 08:00

不是补零问题, 频域在f1/f2上刚好无点!

yametw 发表于 2009-6-4 13:44

恩，就是，原来的频率间隔太大了
页: [1]
查看完整版本: fft画功率谱密度图有问题