dir2ladr.m

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

function [K,C]=dir2ladr(b,a)
%直接型到零极点格型梯形的转换
%C梯形系数；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