casfiltr.m

Matlab环境下带控制面板的低通、高通滤波器设计

function yn=casfiltr(b0,B,A,xn);
[K,L]=size(B);
N=length(xn);
w=zeros(K+1,N);
w(1,:)=xn;
for i=1:K
    w(i+1,:)=filter(B(i,:),A(i,:),w(i,:));
end
yn=b0*w(K+1,:);