latc2dir.m

用DSP译文版课本的程序,相信对学DSP的同学会有所帮助

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);