dir2ladr.m

Digital Signal Processing Using MATLAB V.4随书代码

function [K,C] = dir2ladr(b,a)

% IIR Direct form to pole-zero Lattice/Ladder form Conversion

% -----------------------------------------------------------

% [K,C] = dir2ladr(b,a)

%  K = Lattice coefficients (reflection coefficients), [K1,...,KN]

%  C = Ladder Coefficients, [C0,...,CN]

%  b = Numerator polynomial coefficients (deg <= Num deg)

%  a = Denominator polynomial coefficients

%

a1 = a(1); a = a/a1; b = b/a1;

M = length(b); N = length(a);

if M > N

     error('   *** length of b must be <= length of a ***')

end

b = [b, zeros(1,N-M)]; K = zeros(1,N-1);

A = zeros(N-1,N-1); C = b;

for m = N-1:-1:1

     A(m,1:m) = -a(2:m+1)*C(m+1);

     K(m) = a(m+1);

     J = fliplr(a);

     a = (a-K(m)*J)/(1-K(m)*K(m));

     a = a(1:m);

     C(m) = b(m) + sum(diag(A(m:N-1,1:N-m)));

end