周期图法、Yule-walker方程进行功率谱估计的性能比较

利用周期图法进行谱估计，并绘制结果，窗函数采用矩形窗。利用Levinson-Durbin递推法求解Yule-walker方程，进行AR(6)的建模。与Matlab中periodogram（周期图）和pyulear（Yule-walker方程）中相应方法的结果进行比较和分析。

clc;
clear;

%正弦信号
f1=60;
f2=300;
f3=350;
Fs=1000;%采样频率Hz
N=1000;%采样点数
t=（0:N-1）/Fs;%采样时间s
s1=sin（2*pi*f1*t）;
s2=2*cos（2*pi*f2*t）;
s3=2.5*sin（2*pi*f3*t）;
%加入白噪声
x=awgn（s1+s2+s34）;

NFFT=1024;

X = fft（xNFFT）;
Pxx = abs（X）.^2/N;
index = 0:round（NFFT/2 -1）;
k = index*Fs/NFFT;
P = 10*log（Pxx（index+1））;
figure（1）;
plot（kP）;
title（‘直接计算周期图法‘）;

window = boxcar（N）;     %矩形窗
[Pxx1f]=periodogram（xwindowNFFTFs）;  %将结果转换为dB
figure（2）;
plot（f10*log10（Pxx1））;
hold on;
plot（kP）;
legend（‘MATLAB函数周期图法‘‘无MATLAB函数周期图法‘）;
xlabel（‘Frequency（Hz）‘）;
ylabel（‘PSD （dB/Hz）‘）;
hold off;

noverlap=10;
[Pxx2f1]=pwelch（xboxcar（200）noverlapNFFTFs）;
figure（3）;
plot（f110*log10（Pxx2））;
xlabel（‘Frequency（Hz）‘）;
ylabel（‘PSD （dB/Hz）‘）;
title（‘平均周期图法（Bartlett法）‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

    .......      1655  2020-06-05 10:40  estimation2.m

    .......       584  2020-06-04 17:19  fir_rec.m

    .......       586  2020-06-04 21:06  fir_rec1.m

----------- ---------  ---------- -----  ----

                 2825                    3