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