par2dir.m

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

function [b,a]=par2dir(C,B,A)
%用并联型传输系数计算直接型传输系数
%C为length(b)>=length(a)时的多项式部分
%B级联型传输函数的分子；A级联型传输函数的分母；
%b直接型传输函数的分子；a直接型传输函数的分母
[K,L]=size(A);R=[];P=[];
for i=1:1:K
    [r,p,k]=residuez(B(i,:),A(i,:));
    R=[R;r];P=[P;p];
end
[b,a]=residuez(R,P,C);
b=b(:)';a=a(:)';