longteng1.m

相关函数法该方法先由序列估计出自相关函数R

%相关函数法：该方法先由序列估计出自相关函数R，
%然后对其进行傅立叶变换，得到功率谱估计。
%当延时和数据长度相比很小的时候可以有良好的估计精度。
% author:闫朝星 2007.12.18
% version 0.1
%---------------------------------
Fs=500;%采样频率
n=0:1/Fs:1;
%产生含有噪声的序列
xn=cos(2*pi*40*n)+3*cos(2*pi*90*n)+randn(size(n));
nfft=512;
cxn=xcorr(xn,'unbiased');%计算序列的自相关函数
CXk=fft(cxn,nfft);
Pxx=abs(CXk);
index=0:round(nfft/2-1);
k=index*Fs/nfft;
plot_Pxx=10*log10(Pxx(index+1));
figure(1)
plot(k,plot_Pxx);