levinson.m

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

function [a] = Levinson(k)

len=length(k);
for i=1:len
	a(i,i)=k(i);
	for j=1:i-1,
		a(i,j)=a(i-1,j)-k(i)*a(i-1,i-j);
	end
end
a=[1 a(len,:)];