dir2par.m

来自「原名：数字信号处理基础及MATLAB实现,解压缩密码：lfjk2002@163.」· M 代码 · 共 31 行

function [C,B,A]=dir2par(b,a)
%用直接型传输系数计算并联型传输系数
%C为length(b)>=length(a)时的多项式部分
%B级联型传输函数的分子；A级联型传输函数的分母；
%b直接型传输函数的分子；a直接型传输函数的分母
M=length(b);N=length(a);
[r1,p1,C]=residuez(b,a);          %用residuez函数计算部分分式展开
p=cplxpair(p1,10000000*eps);      %用cplxpair函数实现复数数组分类为复共轭对
I=cplxcomp(p1,p);                 %cplxcomp是扩展函数，配合cplxpair使复数数组互相对应
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),[]);
    B(K,:)=[real(Brow) 0];A(K,:)=[real(Arow) 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,[]);
        B(fix((i+1)/2),:)=real(Brow);
        A(fix((i+1)/2),:)=real(Arow);
    end
end