matlab谱估计

实验目的：

一、.理解参数模型谱估计的基本思路和方法优缺点。

二、掌握利用MATLAB编程进行AR模型谱估计的基本方法、计算步骤和函数使用。

1.产生一个多点频信号，利用上述参数谱估计方法进行功率谱估计；对比经典谱估计与参数谱估计几种方法的处理结果，并进行分析。

程序代码及编程思路：

clear all;

N=128;

p=40;

NFFT=2048;

Fs=2;

n=0:N-1;

randn(‘state’,0);%保持随机数不变

x=cos(0.3pin)+cos(0.4pin)+cos(0.44pin)+randn(size(n));

%产生3个点频的信号

[P,f]=periodogram(x,[],NFFT,2);%周期图法

[Py,fy]=pyulear(x,p,NFFT,2);%AR模型

[Px,F]=pwelch(x,hamming(65),overlap,NFFT);%Welch法

%使用汉明窗，叠合32点，FFT点数为2048

subplot(311);

plot(f,10*log§);

grid on;

title(‘周期图法’);

axis([0 1 -60 60]);

subplot(312);

plot(fy,10*log(Py));

grid on;

title([‘Yule-Walker方法，阶次为’,num2str§]);

axis([0 1 -60 60])

subplot(313)

plot(F,Px);