casfiltr.m

原名：数字信号处理基础及MATLAB实现,解压缩密码：[e-mail]

function y=casfiltr(b0,B,A,x)
%IIR和FIR的级联型的实现
%y为输出序列
%b0为增益常数;B级联型传输函数的分子；A级联型传输函数的分母；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,:);