power_s_e_index.asv

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

clear
Fs=500;
nfft=1024;
N=0.1:0.2:1.1
for i=1:length(N)
n=0:1/Fs:N(i);
w1=2*100*pi;A1=5;
w2=2*150*pi;A2=12;
xn=A1*exp(j*w1*n)+A2*exp(j*w2*n)+1*rands(1);
% subplot(211);
% plot(n,xn);
% xlabel('n');
% ylabel('xn');
% title('xn=A1*exp(jw1n)+A2*exp(jw2n)+e(n)');
% ymax_xn=max(xn)+0.2;
% ymin_xn=min(xn)-0.2;
% axis([0 0.1 ymin_xn ymax_xn]);
p=floor(length(xn)/3);
[xpsd,f]=pburg(xn,p,nfft,Fs);
pmax=max(xpsd);
xpsd=xpsd/pmax;
xpsd=10*log10(xpsd);
subplot(3,2,i);
plot(f,xpsd);
title('Power Spectral estimate with burg');
ylabel('Power Spectral estimate(dB)');
xlabel('f(Hz)');
grid on;
ymin_psd=min(xpsd)-1;
ymax_psd=max(xpsd)+1;
axis([20 Fs/2 ymin_psd ymax_psd]);
end