dir2par.m

这是一个dsp的matlab包

function [C,B,A]=dir2par(b,a);
%变直接形式为级联形式
%[C,B,A]=dir2par(b,a)
%C=当分子多项式阶数大于分母多项式阶数时产生的多项式
%B=k列3行bk系数矩阵
%A=k列3行ak系数矩阵
%a=直接型分子多项式系数
%b=直接型分母多项式系数 

M=length(b);N=length(a);
%将多项式分式分解为多项式与一阶多项式的和的形式
    [r1,p1,C]=residuez(b,a);
%p1为系统极点
    p=cplxpair(p1,10000000*eps);
%重新排序后，就不一定能够将一阶多项式分式的分子分母
%刚好对应，必须对分子系数排序
%函数cplxcomp见后
    I=cplxcomp(p1,p);
    r=r1(I);
%初始化系数矩阵
    K=floor(N/2);B=zeros(K,2);A=zeros(K,3);
 if K*2==N ;
%N是偶数，则A(z)的阶数是奇数，一个因子是一阶的
%将因子合并成二阶多项式分式
        for i=1:2:N-2
           Brow=r(i:1:i+1,:);
           Arow=p(i:1:i+1,:);
           [Brow,Arow]=residuez(Brow,Arow,[]);
           B(fix((i+1)/2),:)=real(Brow);
           A(fix((i+1)/2),:)=real(Arow);
       end
       [Brow,Arow]=residuez(r(N-1),p(N-1),[]);
       A(K,:)=[real(Arow) 0];
       B(K,:)=[real(Brow) 0];
   else  
       for i=1:2:N-1
           Brow=r(i:1:i+1,:);
           Arow=p(i:1:i+1,:);  
           [Brow,Arow]=residuez(Brow,Arow,[]);
           A(fix(i+1)/2,:)=real(Arow);
           B(fix(i+1)/2,:)=real(Brow);
          
       end  
end