sdir2cas.m

原名：数字信号处理基础及MATLAB实现,解压缩密码：[e-mail]

function [C,B,A]=sdir2cas(b,a)
%s平面上直接型到级联型的转换
%C增益系数；B级联子系统分子矩阵；A级联子系统分母矩阵；
%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));K=floor(Na/2);
if K*2==Na       %当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);
    end
elseif Na==1     %当Na=1时计算
    A=[0 real(poly(p))];
else             %当Na为奇及>1时计算
    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         %当Nb=0时计算
    B=[0 0 poly(z)];
elseif K*2==Nb   %当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    %当Nb=1时计算
    B=[0 real(poly(z))];
else            %当Nb为奇且>1时计算
    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