longteng2.m

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

%周期图法是把随机序列x(n)的N个观测数据视为一能量有限的序列，
%直接计算x(n)的离散傅立叶变换，得X(k)，然后再取其幅值的平方，
%并除以N，作为序列x(n)真实功率谱的估计。
% author:闫朝星 2007.12.18
% version 0.1
%---------------------------------
Fs=600;%采样频率
n=0:1/Fs:1;%产生含有噪声的序列
xn=cos(2*pi*40*n)+3*cos(2*pi*90*n)+0.1*randn(size(n));
window=boxcar(length(xn));%矩形窗
nfft=512;
[Pxx,f]=periodogram(xn,window,nfft,Fs);%直接法
plot(f,10*log10(Pxx));
window=boxcar(length(xn));%矩形窗
nfft=1024;
[Pxx,f]=periodogram(xn,window,nfft,Fs);%直接法
figure(1)
plot(f,10*log10(Pxx));