江汉 发表于 2006-8-5 21:50

两种单边自功率谱密度计算的幅值不一样

呵呵,我就是因为这个问题才找到这个论坛的,这也是我第一次发帖,盼望好心人的帮助,谢谢!

第一种方法我用的是matlab里的函数psd
n=1:5000;
gc=sin(40*pi/(5000*0.08)*n*0.08);
Hs=spectrum.periodogram;
psd(Hs,gc,'Fs',1/0.08,'NFFT',512)

第二种方法我先求出时域信号的频谱,再代入△t=1时的双侧谱估计函数,最后乘上我的采样间隔0.08s,再*2倍得到单侧谱
A=fft(gc,512);
for n=0:256
s(n+1)=0.08*(abs(A(n+1)))*(abs(A(n+1)))/(2*pi*5000)*2*2*pi;
end
plot((0:256)/0.16/256,s);grid

两种方法得到的结果,频率一样,就是幅值不一样.我实在检查不出错误了,请各位师兄师姐帮忙,感谢

simon21 发表于 2006-8-6 08:45

这个问题你可以先看一下帖子http://forum.vibunion.com/forum/viewthread.php?tid=17307的讨论,如果还有问题大家可以继续讨论

江汉 发表于 2006-8-6 21:36

谢谢simon21了,原来只能得到相对值,我就先这么理解吧
页: [1]
查看完整版本: 两种单边自功率谱密度计算的幅值不一样