请问在下面程序中功率谱对应的频率是多少?
clear;clc;
N=1024;
randn('state',0)
x=randn(1,N);
nfft=1024;
Sr=abs(fft(x,nfft)).^2/N;
index=0:nfft-1;
k=index*N/nfft;
plot_Sr=10*log10(Sr);%对数功率谱曲线
plot(k,plot_Sr);
是1024Hz吗? N是点数吧?
回复 #2 HolySaint 的帖子
是的回复 #1 dandanzhang000 的帖子
发帖前可以先看看相关讨论http://forum.vibunion.com/forum/viewthread.php?tid=23614&extra=page%3D1%26amp%3Bfilter%3Ddigest
回复 #4 zhlong 的帖子
谢谢zhlong看了相关的贴子,我理解如下:
只有给出采样频率,才能求出其对应的频率,而上面程序并没有给出采样频率,
所以也就无法求出功率谱所对应的频率了吧
请问我的理解对不对呢?
[ 本帖最后由 dandanzhang000 于 2007-8-6 14:18 编辑 ]
回复 #5 dandanzhang000 的帖子
对啊,没有采样频率对应,你画出的频谱图里的每个振动幅值对应的横坐标是第几频率点,如果有采样频率fs,那么用fs/nfft,index=(fs/nfft):(fs/nfft):(fs/nfft*(nfft-1))plot(index,plot_Sr);就可以了
回复 #5 dandanzhang000 的帖子
对的,应该要给出采样频率。
页:
[1]