par2dir.m

来自「matlab应用实例。如：FIR、IIR数字滤波器设计等。」· M 代码 · 共 20 行

function [b,a]=par2dir(C,B,A);
%变并联形式为直接形式
%[b,a]=par2dir(C,B,A);
%C=当分子多项式阶数大于分母多项式阶数时产生的多项式
%B=K列3行bk系数矩阵
%A=K列3行ak系数矩阵
%a=直接型分子多项式系数
%b=直接型分母多项式系数
[K,L]=size(A);
R=[];
P=[];
for i=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(:)';