function [C,B,A] = dir2par(b,a);
% 直接型到并聯(lián)型的轉(zhuǎn)換
% --------------------------------------
% [C,B,A] = dir2par(b,a)
%  C = 當(dāng)length(b) >= length(a)時(shí)的多項(xiàng)式部分
%  B = 包含各bk的K乘2維實(shí)系數(shù)矩陣
%  A = 包含各ak的K乘3維實(shí)系數(shù)矩陣
%  b = 直接型的分子多項(xiàng)式系數(shù)
%  a = 直接型的分母多項(xiàng)式系數(shù)
%
M = length(b); N = length(a);

[r1,p1,C] = residuez(b,a);
p = cplxpair(p1,10000000*eps);
I = cplxcomp(p1,p);
r = r1(I);

K = floor(N/2); B = zeros(K,2); A = zeros(K,3);
if K*2 == N; %N 為偶, A(z) 的次數(shù)為奇,有一個(gè)因子是一階的 
	for i=1:2:N-2
		Brow = r(i:1:i+1,:);
		Arow = p(i:1:i+1,:);
		[Brow,Arow] = residuez(Brow,Arow,[]);
		B(fix((i+1)/2),:) = real(Brow');
		A(fix((i+1)/2),:) = real(Arow');
	end
	[Brow,Arow] = residuez(r(N-1),p(N-1),[]);
	B(K,:) = [real(Brow') 0]; A(K,:) = [real(Arow') 0];
else
        for i=1:2:N-1
		Brow = r(i:1:i+1,:);
		Arow = p(i:1:i+1,:);
		[Brow,Arow] = residuez(Brow,Arow,[]);
		B(fix((i+1)/2),:) = real(Brow');
		A(fix((i+1)/2),:) = real(Arow');
	end     
end