function [b,a] = par2tf(C,B,A);
%《數字信號處理教程——MATLAB釋義與實現》
% 電子工業出版社出版  陳懷琛編著 2004年9月
%
% 并聯型到直接型的轉換子程序
% [b,a] = par2tf(C,B,A)
% ----------------------------------
%  b = 直接型的分子多項式系數
%  a = 直接型的分母多項式系數
%  C = 并聯型的多項式部分
%  B = 包含各bk的K乘2維實系數矩陣
%  A = 包含各ak的K乘3維實系數矩陣
%
[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(:)';