function y=casfiltr(b0,B,A,x)
%IIR和FIR的級聯(lián)型的實(shí)現(xiàn)
%y為輸出序列
%b0為增益常數(shù);B級聯(lián)型傳輸函數(shù)的分子；A級聯(lián)型傳輸函數(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,:);