puguji.m

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

clear all;
clc
close all;
clf
%产生AR信号，以及谐波信号
n=129;

%产生x1信号即AR信号
noise=rands(1,n);
x1=zeros(1,n);
x1(1)=0;x1(2)=0;
rad=0.98;
thita=pi/4;
a1=rad*exp(j*thita);
a2=rad*exp(-j*thita);
for m=3:n
    x1(m)=noise(m)-a1*x1(m-1)-a2*x1(m-2);
end
%产生x2信号即谐波信号
x2=zeros(1,n);
fs=100;f1=35;f2=36;
w1=2*pi*f1/fs;
w2=2*pi*f2/fs;
A1=3;A2=4;
for l=1:n
    x2(l)=A1*exp(j*w1*l)+A2*exp(j*w2*l);
end

%用burg算法得到a1和a2
[a1,gan1,e1]=fburg_n(x1,n,10);
[a2,gan2,e2]=fburg_n(x2,n,10);

%算两者的功率谱
fftx1=fft2(a1);
fftx2=fft2(a2);

H1=abs(fftx1);
H2=abs(fftx2);
p1=H1.*H1;
p2=H2.*H2;
plot(p1); title('p1');
figure;plot(p2);title('p2');