d_k1.m

来自「采用MELP编制的压缩与解码程序」· M 代码 · 共 11 行

%compute the No.1 reflection coefficient
function k=d_k1(p)
k=-p(10);
for n=1:9
   j=11-n;
   for i=1:j-1
      pp(i)=(p(i)-k*p(j-i))/(1-k^2);  %overflow
   end
   p=pp;
   k=-p(j-1);
end