function [b,a] = pfe2rfz(R,P,C)
%在 z域中作部分分式擴(kuò)展所得有理函數(shù)
% --------------------------------------------------------------
% [b,a] = pfe2rfz(R,P,C)
%  b = 有理函數(shù)的分子多項(xiàng)式系數(shù) [b0,b1,...,bM],
%  a = 有理函數(shù)的分母多項(xiàng)式系數(shù) [a0,a1,...,aN] ,
%  R =  p向量中各極點(diǎn)的留數(shù)行向量
%  p = 包含有理函數(shù)極點(diǎn)(a的根)的行向量
%  C = (M-N) 長(zhǎng)(或零長(zhǎng))的包含多項(xiàng)式系數(shù)的向量
%
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);