latc2dir.m

来自「stft(含:短时傅立叶变换 Gabor变换等)」· M 代码 · 共 15 行

function [b] = latc2dir(K)
% All-Zero Lattice form to FIR Direct form Conversion
% ---------------------------------------------------
% [b] = latc2dir(K)
%  b = FIR direct form coefficients (impulse response)
%  K = Lattice filter coefficients (reflection coefficients)
%
M = length(K);
J = 1; A = 1;
for m=2:1:M
	A = [A,0]+conv([0,K(m)],J);
	J = fliplr(A);
end
b=A*K(1);