ar.m

含有噪声的余弦序列

%含有噪声的余弦序列，分别采用周‘期图法’与‘改进协方差法’估计序的功率谱
Fs=1000; %采样频率
%产生序列
n=0:1/Fs:1;
xn=sin(2*pi*50*n)+randn(size(n))
    %周期图法
figure(1)
%周期图法的参数设置
window=hamming(length(xn));
nfft=1024;%FFT的点数
%周期图法计算功率谱
[Pxx,f]=periodogram(xn,window,nfft,Fs);
%绘制并标注图形
plot(f,10*log10(Pxx)) ,grid
xlabel('Frequency(Hz)')
ylabel('Power Spectral Density(dB/Hz)')
title('Periodogram PSD Estimate')
%burg法
figure(2);
%burg法的参数设置
order=20;%AR模型的阶数
%burg法估计并绘制功率谱
pburg(xn,order,nfft,Fs)