function [b,a] = pfe2rfz(R,P,C)
%在 z域中作部分分式擴展所得有理函數
% --------------------------------------------------------------
% [b,a] = pfe2rfz(R,P,C)
%  b = 有理函數的分子多項式系數 [b0,b1,...,bM],
%  a = 有理函數的分母多項式系數 [a0,a1,...,aN] ,
%  R =  p向量中各極點的留數行向量
%  p = 包含有理函數極點(a的根)的行向量
%  C = (M-N) 長(或零長)的包含多項式系數的向量
%
R = R'; p = p';
[b1,a1]=residue(R,P,[]);
b1=fliplr(b1); a1=fliplr(a1);
[r,p,k]=residue(b1,a1);
[b,a]=residue(r,p,fliplr(C));
b=fliplr(b); a=fliplr(a);
b=real(b); a=real(a);
b=b/a(1); a=a/a(1);