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

江汉

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

第一种方法我用的是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＝１时的双侧谱估计函数，最后乘上我的采样间隔０.０８s，再*２倍得到单侧谱
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

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

江汉

谢谢simon21了,原来只能得到相对值，我就先这么理解吧