function [b,a]=par2dir(C,B,A);
%變并聯(lián)形式為直接形式
%[b,a]=par2dir(C,B,A);
%C=當(dāng)分子多項(xiàng)式階數(shù)大于分母多項(xiàng)式階數(shù)時(shí)產(chǎn)生的多項(xiàng)式
%B=K列3行bk系數(shù)矩陣
%A=K列3行ak系數(shù)矩陣
%a=直接型分子多項(xiàng)式系數(shù)
%b=直接型分母多項(xiàng)式系數(shù)
[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(:)';