matlab量化误差分析函数

function [CBA]=tf2par（ba）
%直接型向并联型转换
M=length（b）;N=length（a）;
[r1p1C]=residuez（ba）;%先求系统的单根p1对应的留数r1及直接项C
p=cplxpair（p10.000000001）;%先配对函数cplxpair由p1找共轭复根p，le-9的误差
I=cplxcomp（p1p）;%找p1变为p的排序变化
r=r1（I）;
%变为2阶子系统
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:）;%取出一对极点
        ri=r（i:i+1:）;%取出一对对应留数
        [BiAi]=residue