%《MATLAB在電子信息課程中的應用》第七章dir2par子程序
% 直接型到級聯型轉換
% 電子工業出版社出版  陳懷琛 吳大正 高西全合著 2001年10月

function [C,B,A] = dir2par(b,a);
% 直接型到并聯型的轉換
% --------------------------------------
% [C,B,A] = dir2par(b,a)
%  C = 當length(b) >= length(a)時的多項式部分
%  B = 包含各bk的K乘2 階實系數矩陣
%  A = 包含各ak的K乘3 階實系數矩陣
%  b = 直接型的分子多項式系數
%  a = 直接型的分母多項式系數
%

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) 的階次奇數d, one factor is first order
        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