[求助]请高手指点一下

[hetao]

我现在做的是用直接法估计信号的功率谱密度：
程序如下：
fs=1000;
%generate the sequence with the AWGN noise
n=0:1/fs:1;
xn=cos(2*pi*40*n)+3*cos(2*pi*100*n)+randn(size(n));
%computing the DFT of the sequence
nfft=2048;
XK=fft(xn,nfft);
%computing the PSD
Pxx=abs(XK).^2/length(n);
%generate the figure
index=0:round(nfft/2-1);
k=index*fs/nfft;
k1=index/nfft;
plot_pxx=10*log10(Pxx(index+1));
plot(k1,plot_pxx)
ylabel('power spectrum density ')
xlabel('sequence fs=1000')


问题是：
1：当fs取不同的值的时候，功率谱密度的差距很大，原因？？？？？？？
2：k=index*fs/nfft这个语句我知道每次都要用的，但是是什么意思呢？？？

[songzy41]

1：当fs取不同的值的时候，功率谱密度的差距很大，原因？？？？？？？
2：k=index*fs/nfft这个语句我知道每次都要用的，但是是什么意思呢？？？

对于笫1个问题，楼主应看些数字信号处理的书，了解一下在DFT中存在“泄漏”和“栅栏”这两概念；
对笫2个问题，主要是用来计算频率作为横坐标的标度。

作者: songzy41

1：当fs取不同的值的时候，功率谱密度的差距很大，原因？？？？？？？
2：k=index*fs/nfft这个语句我知道每次都要用的，但是是什么意思呢？？？

对于笫1个问题，楼主应看些数字信号处理的书，了解一下在DFT中存在“泄漏”和“栅栏”这两概念；
对笫2个问题，主要是用来计算频率作为横...

多谢了 马上就看！
当时就没有注意到这两个概念！