par2dir.m

来自「西交大刘树棠写的《数字信号处理》那本书的所有源代码」· M 代码 · 共 8 行

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(:)';