dir2latc.m

来自「原名：数字信号处理基础及MATLAB实现,解压缩密码：lfjk2002@163.」· M 代码 · 共 16 行

function [K]=dir2latc(b)
%FIR的直接型到全零点格型的转换
%K格型滤波器反射系数；b直接型滤波器多项式系数。
M=length(b);
K=zeros(1,M);
b1=b(1);
if b1==0
    error('b(1)等于零')
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