为什么这样求得的自功率谱密度不对呢？

xurundong

用如下的方法求解功率谱密度
clc;
clear;
Fs=1024;
n=0:1/Fs:1;
xn=10*cos(2*pi*40*n);  %+3*cos(2*pi*100*n);
nfft=1024;
window=hamming(100);
noverlap=20;
range='half';
[Pxx,f]=pwelch(xn,window,noverlap,nfft,Fs,range);
figure(1)
plot(f,Pxx);

我认为正确的结果为：10×10/(2*pi)=15.915
而得到的确实如下图