xinhao.m

主要是BURG法对声信号进行谱估计的程序。 还有利用BURG法求AR模型系数

%x(n)+x(n-1)+x(n-2)+x(n-3)=v(n)
clear
Fs=1000;
n=0:1/Fs:.3;
delta=0;A1=5;
omig=pi/3;
e=delta+omig*j;
xn=A1*exp(e*n);
% xn_magnitude=abs(xn)

p=floor(length(xn)/3)+1;
nfft=1024;
[xpsd,f]=pburg(xn,p,nfft,Fs,'half');
pmax=max(xpsd);
xpsd=xpsd/pmax;
xpsd=10*log10(xpsd+0.000001);
subplot(212);
plot(f,xpsd);
title('?ù??Burg???¨???????×????');
ylabel('?????×????(dB)');
xlabel('????(Hz)');
grid on;
ymin=min(xpsd)-2;
ymax=max(xpsd)+2;
axis([0 Fs/2 ymin ymax]);