burg.m

burg算法的matlab实现

function  res=burg(sound,p)      %sound表示采集的语音信号，p是期望的阶                                 数,res是结果

size=length(sound);              %变量大小定义
e=zeros(size,p);
b=zeros(size,p);
a=zeros(size,p);
k=zeros(1,p);

e(1:size,1)=sound;               %burg算法的初始化，e,b表示格形滤波器的各                                 阶状态，算法中的重要参数
b(1:size,1)=sound;

a(1,1)=(2*e(2:size,1)'*b(1:size-1,1))/(e(1:size,1)'*e(1:size,1)+b(1:size-1,1)'*b(1:size-1,1));
k(1,1)=a(1,1);
e(2:size,2)=e(2:size,1)-k(1,1)*b(1:size-1,1);
b(2:size,2)=b(1:size-1,1)-k(1,1)*e(2:size,1);

for  i=2:p                      %计算p阶滤波器参数，实现burg算法，用递归                                 实现
   a(i,i)=(2*e(i+1:size,i)'*b(i:size-1,i))/(e(i+1:size,i)'*e(i+1:size,i)+b(i:size-1,i)'*b(i:size-1,i));
   k(1,i)=a(i,i);

   e(i+1:size,i+1)=e(i+1:size,i)-k(1,i)*b(i:size-1,i);
   b(i+1:size,i+1)=b(i:size-1,i)-k(1,i)*e(i+1:size,i);
   
   for j=i-1:-1:1
       a(j,i)=a(j,i-1)-k(1,i)*a(i-j,i-1);
   end
end
res=a;