parfilt.m

digital signal processing-MATLAB

function y = parfilt(C,B,A,x);
%《数字信号处理教程——MATLAB释义与实现》子程序
% 并联型IIR 滤波器对信号的滤波计算
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
%  [y] = parfilt(C,B,A,x);
% ----------------------------------------
%  y = 输出序列
%  C =当 M >= N时(FIR) 的多项式部分
%  B = 包含各bk的K乘2维实系数矩阵
%  A = 包含各ak的K乘3维实系数矩阵
%  x = 输入序列
%
[K,L] = size(B); N = length(x);
w = zeros(K+1,N);			% 对各子系统输出矩阵初始化
w(1,:) = filter(C,1,x);		% 直接部分的输出
for i = 1:1:K
     w(i+1,:) = filter(B(i,:),A(i,:),x); 	% 各二阶子系统的输出
end
y = sum(w);				% 求总输出