dir2par.m

《数字信号处理实验(MATLAB版)》这

function [C,B,A]=dir2par(num,den)
%直接型到并联型的转换
M=length(num);N=length(den);
[r1,p1,C]=residuez(num,den);%先求系统的单根p1，对应的留数r1及直接项C
p=cplxpair(p1,10000000*eps);%用配对函数cplxpair由p1找共轭复根p
I=cplxcomp(p1,p);             %找p1变为p时的排序变化
r=r1(I);                        %让r1的排序变化为r，保持与极点对应
%变换为二阶子系统
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);%取Brow的实部,放入系数矩阵B的相应行
      A(fix((i+1)/2),:)=real(Arow);%取Arow的实部,放入系数矩阵A的相应行
   end;
   [Brow,Arow]=residuez(r(N-1),p(N-1),[]);%处理实单根
   B(K,:)=[real(Brow),0];A(K,:)=[real(Arow),0];
else                       %N为奇数，A(z)的次数为偶,所有因式都是二阶的
   for i=1:2:N-1
      Brow=r(i:1:i+1,:);     %取出一对留数
      Arow=p(i:1:i+1,:);     %取出一对对应的极点
%二个留数极点转为二阶子系统分子分母系数
      [Brow,Arow]=residuez(Brow,Arow,[]); 
      B(fix((i+1)/2),:)=real(Brow);%取Brow的实部,放入系数矩阵B的相应行
      A(fix((i+1)/2),:)=real(Arow);%取Arow的实部,放入系数矩阵A的相应行
   end
end