sdir2cas.m

来自「这个是我数字信号处理这门课的课程设计,文件包括用matlab编写的代码和设计报告」· M 代码 · 共 56 行

function [C,B,A]=sdir2cas(b,a)
%直接型转换成级联型
%b/a->C*(B/A)
Na=length(a)-1;
Nb=length(b)-1;

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%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
    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