dubin.m

dubin算法的matlab实现

function  res=dubin(sound,p)           %sound表示语音信号，p为计算阶数，res为结果
size=length(sound);                    %变量大小定义
rela=zeros(size+p-1,1);
rela(1:size,1)=sound;                 %计算相关系数矩阵，用FFT实现
s=fft(rela);                  
s=abs(s);
s=s.*s
rela=ifft(s)

e=zeros(size,p);
a=zeros(p,p);
k=zeros(p,1);

e(1,1)=rela(1,1);                   %dubin递归算法初始化
k(1,1)=rela(2,1)/e(1,1);
a(1,1)=k(1,1);
e(2,1)=(1-k(1,1)*k(1,1))*e(1,1);


for i=2:p                         %递归算法到p阶
    jieguo=0;
    for j=1:i
        num=a(j,i-1)*rela(i-j+1,1);;
        jieguo=jieguo+num;
    end
    k(i,1)=(rela(i+1,1)-jieguo)/e(i,1);
    a(i,i)=k(i,1);
    for j=1:i-1
       a(j,i)=a(j,i-1)-k(i,1)*a(i-j,i-1);
    end
    e(i+1,1)=(1-k(i,1)*k(i,1))*e(i,1);
    res=a;
end