dir2ladr.m

Digital Signal Processing with MATLAB是Proakis於DSP的經點著作 內涵此書的source code

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