这几个求功率谱的函数，幅值的差别是因为什么啊？

lfxiansheng

t=0:1/1000:1;
x=sin(2*pi*100*t)+sin(2*pi*200*t)+sin(2*pi*400*t)+randn(size(t));
Len=length(x);NFFT=1024;Fs=1000;
[Pxx F]=periodogram(x,hamming(Len),NFFT,Fs);
figure,plot(F,Pxx.*Fs);xlabel('periodogram方法求自功率谱');ylabel('幅值');
[Pxx F]=pwelch(x,100,80,1024,1000);
figure,plot(F,Pxx.*Fs);xlabel('pwelch方法1求自功率谱');ylabel('幅值');   %这个幅值如何解释？？？
[Pxx F]=pwelch(x,length(x),[],1024,1000);
figure,plot(F,Pxx.*Fs);xlabel('pwelch方法2求自功率谱');ylabel('幅值');   %为什么方法1与方法2有这么大的差别 那%个是正确的啊
[Pxx,F]=psd(x,1024,1000,hamming(length(x)),0);
figure,plot(F,Pxx.*2);xlabel('psd方法求自功率谱');ylabel('幅值');

[ 本帖最后由 lfxiansheng 于 2007-9-24 11:11 编辑 ]

马社

我自己编了求功率谱的函数，我不相信这些函数，太乱了。

hahaer

功率谱 是 估计量
本身 就有 误差的 存在
所以不同的 方法 有 差别是很正常的
通常情况下 平均周期 图 法认为 是 比较 准确的

马社

功率谱 不是 估计量，算出来是多少就多少

94117239

楼上看看数理统计就知道什么是估计量了，不是估计值

VibrationMaster

1.先把所有的随机干扰去掉, 2.用整周期采样(千万要这样). 比较几次就发现规律了