資源簡介
利用周期圖法進行譜估計，并繪制結果，窗函數采用矩形窗。利用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
評論
共有 條評論