casfiltr.m

来自「关于数字信号处理的一些matlab例程」· M 代码 · 共 10 行

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,:)