function y = casfiltr(b0,B,A,x);
%  IIR 和 FIR濾波器的級聯(lián)型的實現(xiàn)
% -----------------------------------------------
%  y = casfiltr(b0,B,A,x);
%  y = 輸出序列
% b0 = 級聯(lián)型的增益系數(shù)
%  B = 包含各bk的K乘3維實系數(shù)矩陣
%  A = 包含各ak的K乘3維實系數(shù)矩陣
%  x = 輸入序列
%
[K,L] = size(B);
N = length(x);
w = zeros(K+1,N);
w(1,:) = x;
for i = 1:1:K
        w(i+1,:) = filter(B(i,:),A(i,:),w(i,:));
end
y = b0*w(K+1,:);