durbin.m

来自「这是一个用于语音信号处理的工具箱」· M 代码 · 共 10 行

function [b, sig2]= durbin(x,n,iq,l)
      
[a sig2] = autocorr(x,n,l);
l1=l+1;
ap(1)=1;
for i=2:l1
	ap(i)=a(i-1);
end	
[b pb] = autocorr(ap,l1,iq);