ladr2dir.m

matlab应用实例。如：FIR、IIR数字滤波器设计等。

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

% 极点－零点Lattice结构到直接型结构的转换
% ------------------------------------------------------------------------------
% [b,a]=ladr2dir(K,C)
% b=直接型的分子多项式系数
% a＝直接型的分母多项式系数
% K=格型系数(反射系数),[K1,K2,...,KN]
% C=梯形系数,[C1,C2,...,CN]
% ------------------------------------------------------------------------------

N=length(K);
M=length(C);
C=[C,zeros(1,N-M+1)];
J=1;
a=1;
A=zeros(N,N);
for m=1:1:N
    a=[a,0]+conv([0,K(m)],J);
    A(m,1:m)=-a(2:m+1);
    J=fliplr(a);
end
b(N+1)=C(N+1);
for m=N:-1:1
    A(m,1:m)=A(m,1:m)*C(m+1);
    b(m)=C(m)-sum(diag(A(m:N,1:N-m+1)));
end