dir21atc.m

来自「关于数字信号处理的一些matlab例程」· M 代码 · 共 17 行

function [K]=dir21atc(b)
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