function?[CBA]=tf2par（ba）
%直接型向并聯(lián)型轉(zhuǎn)換
M=length（b）;N=length（a）;
[r1p1C]=residuez（ba）;%先求系統(tǒng)的單根p1對(duì)應(yīng)的留數(shù)r1及直接項(xiàng)C
p=cplxpair（p10.000000001）;%先配對(duì)函數(shù)cplxpair由p1找共軛復(fù)根p，le-9的誤差
I=cplxcomp（p1p）;%找p1變?yōu)閜的排序變化
r=r1（I）;
%變?yōu)?階子系統(tǒng)
K=floor（N/2）;B=zeros（K2）;A=zeros（K3）;
if?K*2==N
????for?i=1:2:N-2
????????pi=p（i:i+1:）;%取出一對(duì)極點(diǎn)
????????ri=r（i:i+1:）;%取出一對(duì)對(duì)應(yīng)留數(shù)
????????[BiAi]=residue