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;