sdir2cas.m

s平面中直接形式到级联形式的转换 ％适合模拟滤波器的 ％C为增益系数 ％B为包含各bk的K乘3维实系数矩阵 ％A为包含各ak的K乘3维实系数矩阵 ％b为直接形式的分子多项式系数 ％a为直接

function [C,B,A]=sdir2cas(b,a)
%s平面中直接形式到级联形式的转换
%适合模拟滤波器的
%C为增益系数
%B为包含各bk的K乘3维实系数矩阵
%A为包含各ak的K乘3维实系数矩阵
%b为直接形式的分子多项式系数
%a为直接形式的分母多项式系数
%
Na=length(a)-1;
Nb=length(b)-1;
%
%计算增益系数C
%
b0=b(1);
b=b/b0;
a0=a(1);
a=a/a0;
C=b0/a0;
%
%分母的二阶因子部分
%
p=cplxpair(roots(a));
%cplxpair 按共扼复数对对复数分组，排序按实部大小递增；
%而在一个共扼对中，含负虚部的元素排在前面，而且纯实数排在全部复数对的后面
K=floor(Na/2);
%向负无穷方向取整数
if K*2==Na
    A=zeros(K,3);
    for n=1:2:Na
        Arow=p(n:1:n+1,:);
        Arow=poly(Arow);
        A(fix((n+1)/2),:)=real(Arow);
        %a的二阶子系统多项式系数矩阵为A
    end
elseif Na==1
    A=[0 real(poly(p))];
else
    A=zeros(K+1,3);
    for n=1:2:2*K
        Arow=p(n:1:n+1,:);
        Arow=poly(Arow);
        A(fix((n+1)/2),:)=real(Arow);
    end
    A(K+1,:)=[0 real(poly(p(Na)))];
end
%
%分子的二阶因子部分
%
z=cplxpair(roots(b));
K=floor(Nb/2);
if Nb==0
    B=[0 0 poly(z)];
elseif K*2==Nb
    B=zeros(K,3);
    for n=1:2:Nb
        Brow=z(n:1:n+1,:);
        Brow=poly(Brow);
        B(fix((n+1)/2),:)=real(Brow);
    end
elseif Nb==1
    B=[0 real(poly(z))];
else
    B=zeros(K+1,3);
    for n=1:2:2*K
        Brow=z(n:1:n+1,:)
        Brow=poly(Brow);
        B(fix((n+1)/2),:)=real(Brow);
    end
    B(K+1,:)=[0 real(poly(z(Nb)))];
end