par2dir.m

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

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