dir2ladr.m

来自「<Matlab 7辅助信号处理技术与应用》的源码」· M 代码 · 共 36 行

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

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

a1=a(1);
a=a/a1;
b=b/a1;
M=length(b);
N=length(a);
if M>N
    error('***向量b的长度必须<=向量a的长度***')
end
b=[b,zeros(1,N-M)];
K=zeros(1,N-1);
A=zeros(N-1,N-1);
C=b;
for m=N-1:-1:1
    A(m,1:m)=-a(2:m+1)*C(m+1);
    K(m)=a(m+1);
    J=fliplr(a);
    a=(a-K(m)*J)/(1-K(m)*K(m));
    a=a(1:m);
    C(m)=b(m)+sum(diag(A(m:N-1,1:N-m)));
end