par2dir.m

西交大刘树棠写的《数字信号处理》那本书的所有源代码

function [b,a]=par2dir(C,B,A);
[K,L]=size(A);R=[];P=[];
for i=1:1:K
[r,p,k]=residuez(B(i,:),A(i,:));
R=[R;r];P=[P;p];
end
[b,a]=residuez(R,P,C);
b=b(:)';a=a(:)';