par2tf.m

matlab数字信号处理代码

function [b,a] = par2tf(C,B,A);
%《数字信号处理教程——MATLAB释义与实现》
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
% 并联型到直接型的转换子程序
% [b,a] = par2tf(C,B,A)
% ----------------------------------
%  b = 直接型的分子多项式系数
%  a = 直接型的分母多项式系数
%  C = 并联型的多项式部分
%  B = 包含各bk的K乘2维实系数矩阵
%  A = 包含各ak的K乘3维实系数矩阵
%
[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(:)';