资源简介

利用周期图法进行谱估计,并绘制结果,窗函数采用矩形窗。利用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


评论

共有 条评论