dir2latc.m

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

function [K]=dir2latc(b)
%[K]=dir2latc(b);
%K=格型滤波器的系数矩阵
%b=FIR滤波器的直接形式系数矩阵
M=length(b);
K=zeros(1,M);
b1=b(1);
if b1==0
   error('b(1) is equal to zero')
end
K(1)=b1;A=b/b1;
for m=M:-1:2
   K(m)=A(m);
   j=fliplr(A);
   A=(A-K(m)*j)/(1-K(m)*K(m));
   A=A(1:m-1);
End