casfiltr.m

西交大刘树棠写的《数字信号处理》那本书的所有源代码

function y=casfiltr(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,:);